text:
# Bound State and Order Parameter Mixing Effect by Nonmagnetic Impurity Scattering in Two-band Superconductors

## I Introduction

Since the discovery of superconductivity in MgB2, great efforts have been done on this material in order to clarify the character of this 40\[K\]-superconductor. Experiments on the specific heat and penetration depth indicate a finite energy gap below the superconducting transition temperature $`T_\mathrm{c}`$. The presence of a Boron-isotope effect implies that phonon plays a crucial role in the pairing interaction. However, the origin of this superconductivity is still in controversial and various mechanism have been proposed.

Besides the mechanism, possibility of multi-gap superconductivity has been discussed in this material. Indeed, a two-gap structure has been directly observed by tunnelling experiments. The temperature dependence of the specific heat also indicates the presence of two energy gaps. Since MgB2 has two distinct Fermi surfaces originating from $`\sigma `$\- and $`\pi `$-orbitals of boron atoms, it is probable that the superconducting order parameters are different in the two bands. We also mention that multi-gap superconductivity has been investigated also in the superconducting state of transition metals.

When the order parameters are different among bands, it is known that impurities with interband scatterings strongly affect the superconducting state even if impurities are nonmagnetic. In particular, Gusman showed within the Born approximation that the order parameters in a two-band superconductor are mixed with each other by interband scatterings. In the dirty limit, only one averaged energy gap appears in the superconducting density of states. Although this impurity effect was originally considered for the superconductivity in transition metals, this effect can be expected also in MgB2 if this superconductor really has two distinct order parameters in the two bands.

So far, the mixing of order parameters by interband scatterings has been discussed within the Born approximation. With this regard, we mention that higher order scattering processes beyond the Born approximation sometimes lead to qualitatively new effects in the case when the impurity scattering affects superconductivity. (For example, formation of a bound state around a magnetic impurity in $`s`$-wave superconductivity.) Thus it is expected also in the present case that higher order scattering processes may cause new phenomena which cannot be described within the Born approximation.

In this paper we investigate nonmagnetic impurity effects in a model two-band superconductor with two different order parameters. We focus on the mixing of the order parameters by interband scatterings. In contrast to the previous works which employ the Born approximation, we take into account multi-scattering processes in terms of the inter- and intra-band scatterings using the $`t`$-matrix approximation. We clarify how the mixing effect of the order parameters which has been discussed within the Born approximation is modified by higher order scattering processes. We show that, although the interband scatterings are crucial for the mixing effect, at the same time, this effect is suppressed when the interband scattering is very strong; this result is not obtained within the Born approximation. As another new result by the $`t`$-matrix approximation, we also show that a bound state is formed around a nonmagnetic impurity depending on the detailed character of impurity scattering.

This paper is organized as follows: After explaining our model in $`\mathrm{\S }`$2, we consider a single impurity problem and discuss a bound state in $`\mathrm{\S }`$3. In $`\mathrm{\S }`$4, we investigate the case of finite impurity concentration, which is followed by summary in $`\mathrm{\S }`$5.

## II Model Two-band Superconductor

### A Hamiltonian

We consider a three-dimensional two-band $`s`$-wave superconductor. We refer to the two bands as $`a`$\- and $`b`$-band. The Hamiltonian is given by

$`H={\displaystyle \underset{\text{p},\sigma }{}}\epsilon _\text{p}^aa_{\text{p}\sigma }^{}a_{\text{p}\sigma }+{\displaystyle \underset{\text{p},\sigma }{}}\epsilon _\text{p}^bb_{\text{p}\sigma }^{}b_{\text{p}\sigma }`$ $``$ $`g_{aa}{\displaystyle \underset{\text{p},\text{p}^{}}{}}a_\text{p}^{}a_\text{p}^{}a_\text{p}^{}a_\text{p}^{}g_{bb}{\displaystyle \underset{\text{p},\text{p}^{}}{}}b_\text{p}^{}b_\text{p}^{}b_\text{p}^{}b_\text{p}^{}`$ (1)
$``$ $`g_{ab}{\displaystyle \underset{\text{p},\text{p}^{}}{}}a_\text{p}^{}a_\text{p}^{}b_\text{p}^{}b_\text{p}^{}g_{ba}{\displaystyle \underset{\text{p},\text{p}^{}}{}}b_\text{p}^{}b_\text{p}^{}a_\text{p}^{}a_\text{p}^{}`$ (2)
$`+`$ $`H_{\mathrm{imp}}.`$ (3)
Here, $`a_{\text{p}\sigma }`$ ($`b_{\text{p}\sigma }`$) represents the annihilation operator of an electron in the $`a`$-band ($`b`$-band). The kinetic energies of the $`a`$\- and $`b`$-band measured from the Fermi level are, respectively, given by $`\epsilon _\text{p}^a`$ and $`\epsilon _\text{p}^b`$. In this paper, we do not take into account the detailed band structure of MgB2 and simply assume two isotropic three-dimensional bands. The $`a`$-band ($`b`$-band) has an intraband pairing interaction described by $`g_{aa}`$ ($`g_{bb}`$). In addition, eq. (3) also includes an interband pairing interaction described by $`g_{ab}`$ and $`g_{ba}(=g_{ab}^{})`$. We do not specify the origin of the pairing interactions and, for simplicity, assume the same cut-off energy $`\omega _\mathrm{D}`$ for all the pairing interactions. We briefly mention that different cut-off energies between the two bands are crucial in considering the isotope effect as discussed by Kondo. The last term in eq. (3) represents the nonmagnetic impurity scattering:

$`H_{\mathrm{imp}}={\displaystyle \underset{\text{R}_i,\text{p},\text{p}^{},\sigma }{}}\mathrm{e}^{\mathrm{i}(\text{p}\text{p}^{})\text{R}_i}\left[v_{aa}a_{\text{p}\sigma }^{}a_{\text{p}^{}\sigma }+v_{bb}b_{\text{p}\sigma }^{}b_{\text{p}^{}\sigma }+v_{ab}a_{\text{p}\sigma }^{}b_{\text{p}^{}\sigma }+v_{ba}b_{\text{p}\sigma }^{}a_{\text{p}^{}\sigma }\right],`$ (4)
where $`v_{aa}`$ and $`v_{bb}`$, respectively, describe intraband impurity scattering in the $`a`$\- and $`b`$-band while $`v_{ab}`$ and $`v_{ba}(=v_{ab}^{})`$ represent interband scattering between the two bands. In eq. (4), we have neglected momentum dependence of scattering potentials for simplicity and retained the $`s`$-wave scattering component only. In the case of a single impurity problem, we put the position of the impurity $`R_i=0`$ and drop the summation in terms of $`R_i`$ in eq. (4).

We employ the mean field approximation and introduce the two superconducting order parameters in the $`a`$\- and $`b`$-band:

$`\{\begin{array}{c}\mathrm{\Delta }_a=g_{aa}_\text{p}a_\text{p}a_\text{p},\hfill \\ \mathrm{\Delta }_b=g_{bb}_\text{p}b_\text{p}b_\text{p}.\hfill \end{array}`$ (7)
Then the mean field Hamiltonian of eq. (3) is obtained as (we drop constant terms)

$$H_{\mathrm{MF}}=\underset{\text{p}}{}\mathrm{\Psi }_\text{p}^{}(\epsilon _\text{p}^a\rho _3\widehat{\mathrm{\Delta }}_a)\mathrm{\Psi }_\text{p}+\underset{\text{p}}{}\mathrm{\Phi }_\text{p}^{}(\epsilon _\text{p}^b\rho _3\widehat{\mathrm{\Delta }}_b)\mathrm{\Phi }_\text{p}+H_{\mathrm{imp}},$$
(8)
where $`\mathrm{\Psi }_\text{p}^{}(a_\text{p}^{},a_\text{p})`$, and $`\mathrm{\Phi }_\text{p}^{}(b_\text{p}^{},b_\text{p})`$ represent two-component Nambu field operators in the $`a`$\- and $`b`$-band, respectively; $`\rho _i`$ $`(i=1,2,3)`$ are corresponding Pauli matrices. In eq. (8), $`\widehat{\mathrm{\Delta }}_j`$ ($`j=a,b`$) are given by

$`\widehat{\mathrm{\Delta }}_j=\left(\begin{array}{cc}0& \stackrel{~}{\mathrm{\Delta }}_j\\ \stackrel{~}{\mathrm{\Delta }}_j^{}& 0\end{array}\right)=\mathrm{Re}[\stackrel{~}{\mathrm{\Delta }}_j]\rho _1\mathrm{Im}[\stackrel{~}{\mathrm{\Delta }}_j]\rho _2,`$ (11)
where $`\stackrel{~}{\mathrm{\Delta }}_j`$ are related to the order parameters as

$`\{\begin{array}{c}\stackrel{~}{\mathrm{\Delta }}_a=\mathrm{\Delta }_a+{\displaystyle \frac{g_{ab}}{g_{bb}}}\mathrm{\Delta }_b,\hfill \\ \stackrel{~}{\mathrm{\Delta }}_b=\mathrm{\Delta }_b+{\displaystyle \frac{g_{ba}}{g_{aa}}}\mathrm{\Delta }_a.\hfill \end{array}`$ (14)
The impurity term $`H_{\mathrm{imp}}`$ under the Nambu representation is given by

$$H_{\mathrm{imp}}=\underset{\text{R}_i,\text{p},\text{p}^{}}{}\mathrm{e}^{\mathrm{i}(\text{p}\text{p}^{})\text{R}_i}\left[\mathrm{\Psi }_\text{p}^{}\widehat{v}_{aa}\mathrm{\Psi }_\text{p}^{}+\mathrm{\Phi }_\text{p}^{}\widehat{v}_{bb}\mathrm{\Phi }_\text{p}^{}+\mathrm{\Psi }_\text{p}^{}\widehat{v}_{ab}\mathrm{\Phi }_\text{p}^{}+\mathrm{\Phi }_\text{p}^{}\widehat{v}_{ba}\mathrm{\Psi }_\text{p}^{}\right].$$
(15)
Here the matrix impurity potentials are defined by $`\widehat{v}_{aa}=v_{aa}\rho _3`$, $`\widehat{v}_{bb}=v_{bb}\rho _3`$ and

$`\widehat{v}_{ab}=\left(\begin{array}{cc}v_{ab}& 0\\ 0& v_{ab}^{}\end{array}\right)=\mathrm{Re}[v_{ab}]\rho _3+i\mathrm{Im}[v_{ab}]\mathrm{𝟏},\widehat{v}_{ba}=\widehat{v}_{ab}^{}.`$ (18)
We can take $`\mathrm{\Delta }_a`$ real and $`g_{ab}=g_{ba}>0`$ with use of a gauge transformation. Complex quantities in our model are then $`\mathrm{\Delta }_b=|\mathrm{\Delta }_b|\mathrm{e}^{\mathrm{i}\varphi _b}`$ and $`v_{ab}=|v_{ab}|\mathrm{e}^{\mathrm{i}\theta _{ab}}`$.

### B Clean system

In this subsection, we briefly review the theory of two-band superconductivity in the clean system in order to see how the impurity scatterings modify the superconducting state in later sections. The one-particle matrix thermal Green function in the $`j`$-band ($`j=a,b`$) is given by

$$G_j^0(p,\mathrm{i}\omega _m)=\frac{1}{\mathrm{i}\omega _m\epsilon _\text{p}^j\rho _3+\widehat{\mathrm{\Delta }}_j}.$$
(19)
Here $`\omega _m`$ is the Fermion Matsubara frequency. When we calculate the superconducting density of states (SDOS) in the $`j`$-band, we obtain

$$N_j(\omega )=\frac{1}{\pi }\underset{\text{p}}{}\mathrm{Im}G_j^0(p,i\omega _m\omega +i\delta )=N_j(0)\mathrm{Re}\left[\frac{\omega }{\sqrt{\omega ^2|\stackrel{~}{\mathrm{\Delta }}_j|^2}}\right],$$
(20)
where $`N_j(0)`$ is the density of states at the Fermi level in the normal state. Equation (20) shows that the energy gap in SDOS is governed by, not the order parameter $`\mathrm{\Delta }_j`$ itself, but $`\stackrel{~}{\mathrm{\Delta }}_j`$ defined in eq. (14). Since $`\stackrel{~}{\mathrm{\Delta }}_j`$ includes both $`\mathrm{\Delta }_a`$ and $`\mathrm{\Delta }_b`$, the superconducting states in the two bands couple with each other by the interband pairing interaction $`g_{ab}`$.

Substituting eq. (19) into the gap equation $`\mathrm{\Delta }_j=g_{jj}T_{\text{p},\omega _m}G_j^0(p,\mathrm{i}\omega _m)_{12}`$, we obtain

$$\mathrm{\Delta }_j=g_{jj}N_j(0)_0^{\omega _\mathrm{D}}d\epsilon \frac{\stackrel{~}{\mathrm{\Delta }}_j}{\sqrt{\epsilon ^2+|\stackrel{~}{\mathrm{\Delta }}_j|^2}}\mathrm{tanh}\frac{\beta }{2}\sqrt{\epsilon ^2+|\stackrel{~}{\mathrm{\Delta }}_j|^2}(j=a,b),$$
(21)
where $`\beta =1/T`$ is the inverse of temperature. Substituting eq. (14) into eq. (21), we find $`\mathrm{\Delta }_b=C\mathrm{\Delta }_a`$, where $`C`$ is given by

$$C=\frac{1g_{aa}N_a(0)_0^{\omega _\mathrm{D}}d\epsilon \frac{1}{\sqrt{\epsilon ^2+|\stackrel{~}{\mathrm{\Delta }}_a|^2}}\mathrm{tanh}\frac{\beta }{2}\sqrt{\epsilon ^2+|\stackrel{~}{\mathrm{\Delta }}_a|^2}}{g_{aa}\frac{g_{ab}}{g_{bb}}N_a(0)_0^{\omega _\mathrm{D}}d\epsilon \frac{1}{\sqrt{\epsilon ^2+|\stackrel{~}{\mathrm{\Delta }}_a|^2}}\mathrm{tanh}\frac{\beta }{2}\sqrt{\epsilon ^2+|\stackrel{~}{\mathrm{\Delta }}_a|^2}}.$$
(22)
Thus when $`\mathrm{\Delta }_a`$ becomes finite, $`\mathrm{\Delta }_b`$ also becomes finite as far as $`g_{ab}0`$ and $`g_{bb}0`$. We also find that $`\mathrm{\Delta }_b`$ can be also taken real as well as $`\mathrm{\Delta }_a`$ in the clean system.

The superconducting transition temperature $`T_\mathrm{c}`$ is determined by

$$T_\mathrm{c}=\frac{2\gamma \omega _\mathrm{D}}{\pi }\mathrm{exp}\left[\frac{\overline{g}_{aa}+\overline{g}_{bb}\sqrt{(\overline{g}_{aa}\overline{g}_{bb})^2+4\overline{g}_{ab}\overline{g}_{ba}}}{2(\overline{g}_{aa}\overline{g}_{bb}\overline{g}_{ab}\overline{g}_{ba})}\right](\gamma =1.78),$$
(23)
where $`\overline{g}_{aa}g_{aa}N_a(0)`$, $`\overline{g}_{bb}g_{bb}N_a(0)`$, $`\overline{g}_{ab}g_{ab}N_b(0)`$ and $`\overline{g}_{ba}g_{ba}N_a(0)`$.

Temperature dependence of the order parameters in the clean system is shown in Fig. 1. When $`\overline{g}_{ab}(\overline{g}_{aa},\overline{g}_{bb}`$) (Fig. 1(a)), the overall temperature dependence is similar to the temperature dependence of the order parameter in an ordinary single-band BCS superconductor. In weak coupling BCS theory, the magnitude of the energy gap at $`T=0`$ is related to $`T_\mathrm{c}`$ as $`2\mathrm{\Delta }_{\mathrm{BCS}}(T=0)/T_\mathrm{c}=3.54`$. On the other hand, as proved by Soda and Wada, $`2\stackrel{~}{\mathrm{\Delta }}_a(T=0)/T_\mathrm{c}`$ is larger than 3.54 while $`2\stackrel{~}{\mathrm{\Delta }}_b(T=0)/T_\mathrm{c}`$ is smaller than 3.54 in Fig. 1(a). When $`\overline{g}_{ab}(\overline{g}_{aa},\overline{g}_{bb})`$ (Fig. 1(b)), the smaller order parameter ($`\mathrm{\Delta }_b`$) shows a curious temperature dependence: When $`\overline{g}_{ab}=0`$, the $`b`$-band remains in the normal state down to $`T/T_\mathrm{c}=0.513`$ (where $`T_\mathrm{c}`$ is the transition temperature in the case of $`\overline{g}_{ab}=0.002`$). Thus $`\mathrm{\Delta }_b`$ between $`0.513<T/T_\mathrm{c}<1`$ in Fig. 1(b) is due to a proximity effect caused by the interband pairing interaction $`g_{ab}`$. This induced order parameter appears even if the intraband interaction in the $`b`$-band is repulsive ($`g_{bb}<0`$). We shows an example in Fig. 1(c).

## III Bound State around a Single Impurity in Two-band Superconductor

### A Equation of bound state

Let us consider the electronic state around an impurity at $`R_i=0`$. In this case, within the neglect of the spatial variation of the order parameter, the one-particle matrix thermal Green function in the $`a`$-band $`G_a(p,p^{},i\omega _m)`$ has the form

$$G_a(p,p^{},i\omega _m)=G_a^0(p,i\omega _m)\delta _{\text{p},\text{p}^{}}+G_a^0(p,i\omega _m)\widehat{t}_{aa}(\mathrm{i}\omega )G_a^0(p^{},i\omega _m).$$
(24)
Here the $`t`$-matrix $`\widehat{t}_{aa}`$ obeys

$`\{\begin{array}{c}\widehat{t}_{aa}(\mathrm{i}\omega _m)=\widehat{v}_{aa}+\widehat{v}_{aa}_\text{p}G_a^0(p,\mathrm{i}\omega _m)\widehat{t}_{aa}(\mathrm{i}\omega _m)+\widehat{v}_{ab}_\text{p}G_b^0(p,\mathrm{i}\omega _m)\widehat{t}_{ba}(\mathrm{i}\omega _m),\hfill \\ \widehat{t}_{ba}(\mathrm{i}\omega _m)=\widehat{v}_{ba}+\widehat{v}_{ba}_\text{p}G_a^0(p,\mathrm{i}\omega _m)\widehat{t}_{aa}(\mathrm{i}\omega _m)+\widehat{v}_{bb}_\text{p}G_b^0(p,\mathrm{i}\omega _m)\widehat{t}_{ba}(\mathrm{i}\omega _m).\hfill \end{array}`$ (27)
The solution of eq.(27) is given by

$$\widehat{t}_{aa}(\mathrm{i}\omega _m)=[1\widehat{V}_{aa}(\mathrm{i}\omega _m)\underset{\text{p}}{}G_{aa}^0(p,\mathrm{i}\omega _m)]^1\widehat{V}_{aa}(\mathrm{i}\omega _m),$$
(28)
where the effective impurity potential in the $`a`$-band defined by $`\widehat{V}_{aa}(\mathrm{i}\omega _m)`$ includes interband scatterings between the two bands as

$$\widehat{V}_{aa}(\mathrm{i}\omega _m)=\widehat{v}_{aa}+\widehat{v}_{ab}\underset{\text{p}}{}G_b^0(p,\mathrm{i}\omega _m)[1\widehat{v}_{bb}\underset{\text{p}}{}G_b^0(p,\mathrm{i}\omega _m)]^1\widehat{v}_{ba}.$$
(29)
A bound state caused by the impurity scattering is obtained as a pole of the analytic-continued $`t`$-matrix of eq. (28). The equation of the bound state energy $`\omega `$ is then given by ($`\omega <\mathrm{Min}[|\stackrel{~}{\mathrm{\Delta }}_a|,|\stackrel{~}{\mathrm{\Delta }}_b|]`$)

$`{\displaystyle \frac{\omega ^2\mathrm{cos}(2\theta _{ab})\stackrel{~}{\mathrm{\Delta }}_a\stackrel{~}{\mathrm{\Delta }}_b}{\sqrt{\stackrel{~}{\mathrm{\Delta }}_a^2\omega ^2}\sqrt{\stackrel{~}{\mathrm{\Delta }}_b^2\omega ^2}}}={\displaystyle \frac{1}{2|\overline{v}_{ab}|^2}}(1+\alpha _a^2+\alpha _b^2+(|\overline{v}_{ab}|^2\alpha _a\alpha _b)^2),`$ (30)
where $`\alpha _a\pi N_a(0)v_{aa}`$ and $`\alpha _b\pi N_b(0)v_{bb}`$ represent intraband impurity potentials normalized by the density of states; $`\overline{v}_{ab}\pi \sqrt{N_a(0)N_b(0)}v_{ab}=|\overline{v}_{ab}|\mathrm{e}^{\mathrm{i}\theta _{ab}}`$ . Since we consider a single impurity problem and neglect spatial variation of the order parameters, we have taken both $`\stackrel{~}{\mathrm{\Delta }}_a`$ and $`\stackrel{~}{\mathrm{\Delta }}_b`$ real in eq. (30). The derivation of eq. (30) is summarized in Appendix A.

### B Conditions for the existence of bound state

Since the present theory is symmetric in terms of the $`a`$\- and $`b`$-band, we may consider only the case of $`|\stackrel{~}{\mathrm{\Delta }}_a||\stackrel{~}{\mathrm{\Delta }}_b|`$ without loss of generality. The left hand side (LHS) in eq. (30) equals $`\mathrm{cos}(2\theta _{ab})\mathrm{sign}[\stackrel{~}{\mathrm{\Delta }}_a\stackrel{~}{\mathrm{\Delta }}_b]`$ at $`\omega =0`$ while it becomes $`\pm \mathrm{}`$ at $`\omega =|\stackrel{~}{\mathrm{\Delta }}_b|`$ depending on the sign of $`\stackrel{~}{\mathrm{\Delta }}_b^2\mathrm{cos}(2\theta _{ab})\stackrel{~}{\mathrm{\Delta }}_a\stackrel{~}{\mathrm{\Delta }}_b`$. We then find the following three cases about the formation of the bound state:

1. $`\mathrm{cos}(2\theta _{ab})\mathrm{sign}[\stackrel{~}{\mathrm{\Delta }}_a\stackrel{~}{\mathrm{\Delta }}_b]0`$: A bound state is obtained if the right hand side (RHS) in eq. (30) is larger than $`|\mathrm{cos}(2\theta _{ab})|`$.
2. $`\mathrm{cos}(2\theta _{ab})\mathrm{sign}[\stackrel{~}{\mathrm{\Delta }}_a\stackrel{~}{\mathrm{\Delta }}_b]>0`$ and $`\stackrel{~}{\mathrm{\Delta }}_b^2\mathrm{cos}(2\theta _{ab})\stackrel{~}{\mathrm{\Delta }}_a\stackrel{~}{\mathrm{\Delta }}_b>0`$: A bound state exists because LHS in eq. (30) monotonously increases from $`\mathrm{cos}(2\theta _{ab})`$ to $`+\mathrm{}`$ as $`\omega `$ changes from $`0`$ to $`|\stackrel{~}{\mathrm{\Delta }}_b|`$.
3. $`\mathrm{cos}(2\theta _{ab})\mathrm{sign}[\stackrel{~}{\mathrm{\Delta }}_a\stackrel{~}{\mathrm{\Delta }}_b]>0`$ and $`\stackrel{~}{\mathrm{\Delta }}_b^2\mathrm{cos}(2\theta _{ab})\stackrel{~}{\mathrm{\Delta }}_a\stackrel{~}{\mathrm{\Delta }}_b0`$: A bound state is absent below $`|\stackrel{~}{\mathrm{\Delta }}_b|`$ because LHS in eq. (30) is always negative.

A bound state is not obtained when the interband potential scattering $`\overline{v}_{ab}`$ is absent because RHS in eq. (30) goes infinity. In addition, when $`\stackrel{~}{\mathrm{\Delta }}_a`$ and $`\stackrel{~}{\mathrm{\Delta }}_b`$ have the same sign and the interband potential scattering is real, one does not obtain a bound state. Thus, the phase of the interband potential scattering and the difference of the signs of $`\stackrel{~}{\mathrm{\Delta }}_a`$ and $`\stackrel{~}{\mathrm{\Delta }}_b`$ are crucial in obtaining the bound state.

Although the interband scattering is crucial for the bound state, eq. (30) also shows that the bound state disappears in the limit $`|\overline{v}_{ab}|\mathrm{}`$. (Note that RHS in eq. (30) goes infinity in this limit.) When we write the interband scattering term as $`[v_{ab}a_{0\sigma }^{}b_{0\sigma }+h.c.]`$ ($`a_{0\sigma }^{}`$ ($`b_{0\sigma }`$): creation (annihilation) operator of an electron in the $`a`$-band ($`b`$-band) at $`R=0`$) and diagonalize this term, we obtain two energy levels $`E_\pm =\pm |v_{ab}|`$. In the limit $`|v_{ab}|\mathrm{}`$, the lower level $`E_{}`$ is occupied by two electrons while the higher state $`E_+`$ becomes irrelevant owing to its extremely high energy. Then since electrons no longer come to the impurity site, the effect of the interband scattering disappears in the limit $`|\overline{v}_{ab}|\mathrm{}`$.

Dependence of the bound state energy on $`|\overline{v}_{ab}|`$ is shown in Fig. 2. A bound state appears when the interband scattering becomes finite; however, the bound state energy again approaches the energy gap $`|\stackrel{~}{\mathrm{\Delta }}_b|`$ when $`|\stackrel{~}{v}_{ab}|1`$ because of the reason explained above. In Fig. 1(a), the bound state energy increases as the phase of interband impurity scattering $`\theta _{ab}`$ decreases. In the case of Fig. 1(a) ($`\stackrel{~}{\mathrm{\Delta }}_a=2`$ and $`\stackrel{~}{\mathrm{\Delta }}_b=1`$), the bound state disappears for $`\theta _{ab}\pi /6`$ because the case-(3) discussed above is satisfied. The bound state energy also increases when the intraband scattering ($`\alpha _b`$) becomes strong, as shown in Fig. 1(b). This is because electrons cannot come to the impurity site when the intraband impurity potential is strong; thus the effect of the interband scattering is weakened by the intraband impurity scattering. We also find this effect in eq. (30); when $`\alpha _a\mathrm{}`$ and/or $`\alpha _b\mathrm{}`$, RHS in eq. (30) goes infinity as in the case of $`|\overline{v}_{ab}|0`$.

## IV Order Parameter Mixing in Dirty Two-band Superconductors

### A Formulation: $`t`$-matrix approximation

In this section, we consider the case of finite impurity concentration using the $`t`$-matrix approximation. After taking spatial average in terms of the distribution of impurities, we obtain the renormalized one-particle matrix thermal Green function in the $`a`$-band as

$`G_a(p,\mathrm{i}\omega _m)`$ $`=`$ $`{\displaystyle \frac{1}{\mathrm{i}\omega _m\epsilon _\text{p}^a\rho _3+\widehat{\mathrm{\Delta }}_a\mathrm{\Sigma }_{aa}(\mathrm{i}\omega _m)}}`$ (31)
$``$ $`{\displaystyle \frac{1}{\mathrm{i}\overline{\omega }_m^a\epsilon _\text{p}^a\rho _3+\widehat{\overline{\mathrm{\Delta }}}_a}}.`$ (32)
The self-energy $`\mathrm{\Sigma }_{aa}(\mathrm{i}\omega _m)`$ within the $`t`$-matrix approximation is given by

$`\{\begin{array}{c}\mathrm{\Sigma }_{aa}(\mathrm{i}\omega _m)=n_{\mathrm{imp}}[1\overline{V}_{aa}(\mathrm{i}\omega _m)_\text{p}G_a(p,\mathrm{i}\omega _m)]^1\overline{V}_{aa},\hfill \\ \overline{V}_{aa}(\mathrm{i}\omega _m)=\widehat{v}_{aa}+\widehat{v}_{ab}_\text{p}G_b(p,\mathrm{i}\omega _m)[1\widehat{v}_{bb}G_b(p,\mathrm{i}\omega _m)]^1\widehat{v}_{ba}.\hfill \end{array}`$ (35)
Here $`n_{\mathrm{imp}}`$ is the number density of impurities. In eq. (35), $`\overline{V}_{aa}`$ has the same structure as $`\widehat{V}_{aa}`$ in eq. (29) except that the renormalized Green function is now used in place of $`G_b^0`$. The renormalized Green function in the $`b`$-band is also given by eqs. (32) and (35) where the indices ’$`a`$’ and ’$`b`$’ are interchanged. We mention that $`\mathrm{\Sigma }_{aa}`$ has the same structure as $`\widehat{t}_{aa}`$ in eq. (28).

In Eq. (32), $`\overline{\omega }_m^j`$ is a renormalized Matsubara frequency and

$`\widehat{\overline{\mathrm{\Delta }}}_j=\left(\begin{array}{cc}0& \overline{\stackrel{~}{\mathrm{\Delta }}}_j\\ \overline{\stackrel{~}{\mathrm{\Delta }}}_j^{}& 0\end{array}\right)=\mathrm{Re}[\overline{\stackrel{~}{\mathrm{\Delta }}}_j]\rho _1\mathrm{Im}[\overline{\stackrel{~}{\mathrm{\Delta }}}_j]\rho _2`$ (38)
with $`\overline{\stackrel{~}{\mathrm{\Delta }}}_j`$ being a renormalized $`\stackrel{~}{\mathrm{\Delta }}_j`$. Their expressions are derived in Appendix B and we only show the results here:

$`\{\begin{array}{c}\mathrm{i}\overline{\omega }_m^a=\mathrm{i}\omega _m+{\displaystyle \frac{u_a}{D_a}}\left[{\displaystyle \frac{|\overline{v}_{ab}|^2}{1+\alpha _b^2}}{\displaystyle \frac{\mathrm{i}\overline{\omega }_m^b}{\overline{\mathrm{\Omega }}_b}}+\eta _a{\displaystyle \frac{\mathrm{i}\overline{\omega }_m^a}{\overline{\mathrm{\Omega }}_a}}\right],\hfill \\ \mathrm{i}\overline{\omega }_m^b=\mathrm{i}\omega _m+{\displaystyle \frac{u_b}{D_b}}\left[{\displaystyle \frac{|\overline{v}_{ba}|^2}{1+\alpha _a^2}}{\displaystyle \frac{\mathrm{i}\overline{\omega }_m^a}{\overline{\mathrm{\Omega }}_a}}+\eta _b{\displaystyle \frac{\mathrm{i}\overline{\omega }_m^b}{\overline{\mathrm{\Omega }}_b}}\right],\hfill \end{array}`$ (41)
$`\{\begin{array}{c}\overline{\stackrel{~}{\mathrm{\Delta }}}_a=\stackrel{~}{\mathrm{\Delta }}_a+{\displaystyle \frac{u_a}{D_a}}\left[{\displaystyle \frac{\overline{v}_{ab}^2}{1+\alpha _b^2}}{\displaystyle \frac{\overline{\stackrel{~}{\mathrm{\Delta }}}_b}{\overline{\mathrm{\Omega }}_b}}+\eta _a{\displaystyle \frac{\overline{\stackrel{~}{\mathrm{\Delta }}}_a}{\overline{\mathrm{\Omega }}_a}}\right],\hfill \\ \overline{\stackrel{~}{\mathrm{\Delta }}}_b=\stackrel{~}{\mathrm{\Delta }}_b+{\displaystyle \frac{u_b}{D_b}}\left[{\displaystyle \frac{\overline{v}_{ba}^2}{1+\alpha _a^2}}{\displaystyle \frac{\overline{\stackrel{~}{\mathrm{\Delta }}}_a}{\overline{\mathrm{\Omega }}_a}}+\eta _b{\displaystyle \frac{\overline{\stackrel{~}{\mathrm{\Delta }}}_b}{\overline{\mathrm{\Omega }}_b}}\right],\hfill \end{array}`$ (44)
where $`\overline{\mathrm{\Omega }}_j=\sqrt{\overline{\omega }_m^j{}_{}{}^{2}+|\overline{\stackrel{~}{\mathrm{\Delta }}}_j|^2}`$, $`u_j=n_{\mathrm{imp}}/(\pi N_j(0))`$ and

$`\{\begin{array}{c}\eta _a={\displaystyle \frac{|\overline{v}_{ab}|^4}{(1+\alpha _b^2)^2}}+\left(\alpha _a{\displaystyle \frac{\alpha _b}{1+\alpha _b^2}}|\overline{v}_{ab}|^2\right)^2,\hfill \\ \eta _b={\displaystyle \frac{|\overline{v}_{ba}|^4}{(1+\alpha _a^2)^2}}+\left(\alpha _b{\displaystyle \frac{\alpha _a}{1+\alpha _a^2}}|\overline{v}_{ba}|^2\right)^2,\hfill \end{array}`$ (47)
$`\{\begin{array}{c}D_a=1+\eta _a+{\displaystyle \frac{2}{1+\alpha _b^2}}{\displaystyle \frac{1}{\overline{\mathrm{\Omega }}_a\overline{\mathrm{\Omega }}_b}}\left[\overline{\omega }_m^a\overline{\omega }_m^b|\overline{v}_{ab}|^2+\mathrm{Re}[\overline{v}_{ab}^2\overline{\stackrel{~}{\mathrm{\Delta }}}_a^{}\overline{\stackrel{~}{\mathrm{\Delta }}}_b]\right],\hfill \\ D_b=1+\eta _b+{\displaystyle \frac{2}{1+\alpha _a^2}}{\displaystyle \frac{1}{\overline{\mathrm{\Omega }}_a\overline{\mathrm{\Omega }}_b}}\left[\overline{\omega }_m^a\overline{\omega }_m^b|\overline{v}_{ba}|^2+\mathrm{Re}[\overline{v}_{ba}^2\overline{\stackrel{~}{\mathrm{\Delta }}}_b^{}\overline{\stackrel{~}{\mathrm{\Delta }}}_a]\right].\hfill \end{array}`$ (50)
### B Gap equation

The order parameter in the $`j`$-band, $`\mathrm{\Delta }_j`$, is determined by the gap equation $`\mathrm{\Delta }_j=g_{jj}T_{\text{p},\omega _m}G_j(p,\mathrm{i}\omega _m)_{12}`$:

$$\mathrm{\Delta }_j=\overline{g}_{jj}\frac{\pi }{\beta }\underset{\omega _m}{}\frac{\overline{\stackrel{~}{\mathrm{\Delta }}}_j}{\sqrt{\overline{\omega }_m^j{}_{}{}^{2}+|\overline{\stackrel{~}{\mathrm{\Delta }}}_j|^2}}(j=a,b).$$
(51)
As usual, we introduce a cutoff $`\omega _\mathrm{D}`$ in calculating the summation of the Matsubara frequency in eq. (51).

When the interband scattering $`v_{ab}`$ is absent, $`\overline{\omega }_m^j`$ and $`\overline{\stackrel{~}{\mathrm{\Delta }}}_j`$ satisfy

$$\frac{\overline{\omega }_m^j}{\overline{\stackrel{~}{\mathrm{\Delta }}}_j}=\frac{\omega _m}{\stackrel{~}{\mathrm{\Delta }}_j}.$$
(52)
In this case, the effect of impurity scatterings in eq. (51) vanishes as in the case of the nonmagnetic impurity effect in single-band $`s`$-wave superconductivity. As a result, superconductivity is unaffected by impurities. The same situation also occurs when (1) intraband scatterings or (2) interband scatterings are very strong: In the former case, since the first terms in $`[\mathrm{}]`$ in eqs. (41) and (44) are absent in the limit $`\alpha _j\mathrm{}`$, eq. (52) is satisfied. Physically, when the intraband impurity potential is strong, electrons are excluded from the impurity site. Then the effect of the interband scattering is effectively weakened, so that the situation becomes close to the case of $`v_{ab}=0`$. In the case (2), since $`\eta _j|v_{ab}|^4`$ for $`|\overline{v}_{ab}|1`$, the first terms in $`[\mathrm{}]`$ in eqs. (41) and (44) again vanish in the limit $`\overline{v}_{ab}\mathrm{}`$. As discussed in Sec. 3.2, the impurity level becomes $`\pm \mathrm{}`$ in this limit, so that the interband scattering becomes irrelevant. As a result, the situation in the limit $`v_{ab}\mathrm{}`$ is equal to the case of $`v_{ab}=0`$.

Within the Born approximation, the effect of the interband scattering is controlled by the damping rate $`\mathrm{\Gamma }_j=\pi n_{\mathrm{imp}}N_j(0)|v_{ab}|^2`$ ($`j=a,b`$) when the intraband scattering is absent. Then as shown by Gusman only one averaged energy gap appears in SDOS in the limit $`\mathrm{\Gamma }_j\mathrm{}`$ ($`j=a,b`$). (We reproduce this result based in the present model in Appendix C.) Namely, within the Born approximation, the interband scattering simply promotes the mixing of the order parameters. On the other hand, the present results show that this does not always occur; in the limit $`v_{ab}\mathrm{}`$, although $`\mathrm{\Gamma }_j\mathrm{}`$, a two-gap structure corresponding to $`|\stackrel{~}{\mathrm{\Delta }}_a|`$ and $`|\stackrel{~}{\mathrm{\Delta }}_b|`$ should be observed in SDOS as in the case of $`v_{ab}=0`$.

In solving eq. (51), we employ the following prescription: When we separate eq. (51) into the real and imaginary parts, we obtain four equations. However, since we can always take $`\mathrm{\Delta }_a`$ real, the number of variables which we have to determine is three, i.e., $`\mathrm{\Delta }_a`$, $`\mathrm{Re}[\mathrm{\Delta }_b]`$ and $`\mathrm{Im}[\mathrm{\Delta }_b]`$. Thus we have to eliminate the redundancy in eq. (51). For this purpose, we rewrite the gap equation as follows: From (41) and (44), we find that $`\overline{\omega }_m^j`$ and $`\overline{\stackrel{~}{\mathrm{\Delta }}}_j`$ can be written as

$`\{\begin{array}{c}\mathrm{i}\overline{\omega }_m^j=\mathrm{i}\omega _m{\displaystyle \frac{\zeta _j}{F_\omega }},\hfill \\ \overline{\stackrel{~}{\mathrm{\Delta }}}_j={\displaystyle \frac{\overline{\stackrel{~}{\mathrm{\Delta }}}_j^{}}{F_\omega }},\hfill \end{array}`$ (55)
where $`\overline{\stackrel{~}{\mathrm{\Delta }}}_j^{}`$ $`(j=a,b)`$ are related to $`\mathrm{\Delta }_a`$ and $`\mathrm{\Delta }_b`$ as

$`\{\begin{array}{c}\overline{\stackrel{~}{\mathrm{\Delta }}}_a^{}=\zeta _{aa}\mathrm{\Delta }_a+{\displaystyle \frac{\overline{g}_{ab}}{\overline{g}_{bb}}}\zeta _{ab}\mathrm{\Delta }_b,\hfill \\ \overline{\stackrel{~}{\mathrm{\Delta }}}_b^{}=\zeta _{bb}\mathrm{\Delta }_b+{\displaystyle \frac{\overline{g}_{ba}}{\overline{g}_{aa}}}\zeta _{ba}\mathrm{\Delta }_a.\hfill \end{array}`$ (58)
In eq. (55), $`F_\omega `$ is given by

$$F_\omega =\left(1\frac{u_a}{D_a}\frac{\eta _a}{\overline{\mathrm{\Omega }}_a}\right)\left(1\frac{u_b}{D_b}\frac{\eta _b}{\overline{\mathrm{\Omega }}_b}\right)\frac{|\overline{v}_{ab}|^4}{(1+\alpha _a^2)(1+\alpha _b^2)}\frac{u_au_b}{D_aD_b}\frac{1}{\overline{\mathrm{\Omega }}_a\overline{\mathrm{\Omega }}_b}.$$
(59)
In eqs. (55) and (58), $`\zeta _j`$ and $`\zeta _{jj^{}}`$ $`(j,j^{}=a,b)`$ are given by

$`\{\begin{array}{c}\zeta _a=1\left[{\displaystyle \frac{u_b}{D_b}}\eta _b{\displaystyle \frac{u_a}{D_a}}{\displaystyle \frac{|\overline{v}_{ab}|^2}{1+\alpha _b^2}}\right]{\displaystyle \frac{1}{\overline{\mathrm{\Omega }}_b}},\hfill \\ \zeta _b=1\left[{\displaystyle \frac{u_a}{D_a}}\eta _a{\displaystyle \frac{u_b}{D_b}}{\displaystyle \frac{|\overline{v}_{ba}|^2}{1+\alpha _a^2}}\right]{\displaystyle \frac{1}{\overline{\mathrm{\Omega }}_a}},\hfill \end{array}`$ (62)
$`\{\begin{array}{c}\zeta _{aa}=1\left[{\displaystyle \frac{u_b}{D_b}}\eta _b{\displaystyle \frac{\overline{g}_{ba}}{\overline{g}_{aa}}}{\displaystyle \frac{u_a}{D_a}}{\displaystyle \frac{\overline{v}_{ab}^2}{1+\alpha _b^2}}\right]{\displaystyle \frac{1}{\overline{\mathrm{\Omega }}_b}},\hfill \\ \zeta _{ab}=1\left[{\displaystyle \frac{u_b}{D_b}}\eta _b{\displaystyle \frac{\overline{g}_{bb}}{\overline{g}_{ab}}}{\displaystyle \frac{u_a}{D_a}}{\displaystyle \frac{\overline{v}_{ab}^2}{1+\alpha _b^2}}\right]{\displaystyle \frac{1}{\overline{\mathrm{\Omega }}_b}},\hfill \\ \zeta _{bb}=1\left[{\displaystyle \frac{u_a}{D_a}}\eta _a{\displaystyle \frac{\overline{g}_{ab}}{\overline{g}_{bb}}}{\displaystyle \frac{u_b}{D_b}}{\displaystyle \frac{\overline{v}_{ba}^2}{1+\alpha _a^2}}\right]{\displaystyle \frac{1}{\overline{\mathrm{\Omega }}_a}},\hfill \\ \zeta _{ba}=1\left[{\displaystyle \frac{u_a}{D_a}}\eta _a{\displaystyle \frac{\overline{g}_{aa}}{\overline{g}_{ba}}}{\displaystyle \frac{u_b}{D_b}}{\displaystyle \frac{\overline{v}_{ba}^2}{1+\alpha _a^2}}\right]{\displaystyle \frac{1}{\overline{\mathrm{\Omega }}_a}}.\hfill \end{array}`$ (67)
Substituting eq. (55) into the gap equation (51), we obtain

$`\left(\begin{array}{cc}1\overline{g}_{aa}{\displaystyle \frac{\pi }{\beta }}{\displaystyle \underset{\omega _m}{}}s(\omega _m){\displaystyle \frac{\zeta _{aa}}{\overline{\mathrm{\Omega }}_a^{}}}& \overline{g}_{aa}{\displaystyle \frac{\overline{g}_{ab}}{\overline{g}_{bb}}}{\displaystyle \frac{\pi }{\beta }}{\displaystyle \underset{\omega _m}{}}s(\omega _m){\displaystyle \frac{\zeta _{ab}}{\overline{\mathrm{\Omega }}_a^{}}}\\ \overline{g}_{bb}{\displaystyle \frac{\overline{g}_{ba}}{\overline{g}_{aa}}}{\displaystyle \frac{\pi }{\beta }}{\displaystyle \underset{\omega _m}{}}s(\omega _m){\displaystyle \frac{\zeta _{ba}}{\overline{\mathrm{\Omega }}_b^{}}}& 1\overline{g}_{bb}{\displaystyle \frac{\pi }{\beta }}{\displaystyle \underset{\omega _m}{}}s(\omega _m){\displaystyle \frac{\zeta _{bb}}{\overline{\mathrm{\Omega }}_b^{}}}\end{array}\right)\left(\begin{array}{c}\mathrm{\Delta }_a\hfill \\ \mathrm{\Delta }_b\hfill \end{array}\right)=\left(\begin{array}{c}0\hfill \\ 0\hfill \end{array}\right).`$ (74)
Here $`\overline{\mathrm{\Omega }}_j^{}=\sqrt{\zeta _j^2\omega _m^2+|\overline{\stackrel{~}{\mathrm{\Delta }}}_j^{}|^2}`$ and $`s(\omega _m)=F_\omega /|F_\omega |`$. Thus, the determinant of the $`2\times 2`$-matrix in eq. (74) must be zero in order to obtain finite $`\mathrm{\Delta }_a`$ and $`\mathrm{\Delta }_b`$:

$$\left(1\overline{g}_{aa}\frac{\pi }{\beta }\underset{\omega _m}{}s(\omega _m)\frac{\zeta _{aa}}{\overline{\mathrm{\Omega }}_a^{}}\right)\left(1\overline{g}_{bb}\frac{\pi }{\beta }\underset{\omega _m}{}s(\omega _m)\frac{\zeta _{bb}}{\overline{\mathrm{\Omega }}_b^{}}\right)\overline{g}_{ab}\overline{g}_{ba}\frac{\pi }{\beta }\underset{\omega _m}{}s(\omega _m)\frac{\zeta _{ab}}{\overline{\mathrm{\Omega }}_a^{}}\frac{\pi }{\beta }\underset{\omega _m^{}}{}s(\omega _m^{})\frac{\zeta _{ba}}{\overline{\mathrm{\Omega }}_b^{}}=0.$$
(75)
When eq. (75) is satisfied, the ”eigen-value” $`(\mathrm{\Delta }_a,\mathrm{\Delta }_b)`$ in eq. (74) obeys the equation

$$\frac{\mathrm{\Delta }_b}{\mathrm{\Delta }_a}=\frac{1\overline{g}_{aa}{\displaystyle \frac{\pi }{\beta }}{\displaystyle \underset{\omega _m}{}}s(\omega _m){\displaystyle \frac{\zeta _{aa}}{\overline{\mathrm{\Omega }}_a^{}}}}{\overline{g}_{aa}{\displaystyle \frac{\overline{g}_{ab}}{\overline{g}_{bb}}}{\displaystyle \frac{\pi }{\beta }}{\displaystyle \underset{\omega _m}{}}s(\omega _m){\displaystyle \frac{\zeta _{ab}}{\overline{\mathrm{\Omega }}_b^{}}}}.$$
(76)
The real and imaginary parts of Eq. (76) give two equations, while eq. (75) gives only one equation because it is a real equation as proved in Appendix D. Thus we can determine the order parameters $`\mathrm{\Delta }_a`$, $`\mathrm{Re}[\mathrm{\Delta }_b]`$ and $`\mathrm{Im}[\mathrm{\Delta }_b]`$ unambiguously from eqs. (75) and (76).

In the clean system, we can take $`\zeta _j=\zeta _{jj^{}}=F_\omega =1`$. Then the gap equation (74) is reduced to eq. (21) when, as usual, (1) we transform the summations in eq. (74) to complex integrations $`z`$, (2) change the variables of the integrations as $`z\epsilon \sqrt{z^2|\stackrel{~}{\mathrm{\Delta }}_j|^2}`$ and (3) introduce a cutoff $`\omega _\mathrm{D}`$.

### C Transition temperature

The transition temperature $`T_\mathrm{c}`$ is determined by eq. (75) with $`\mathrm{\Delta }_a=\mathrm{\Delta }_b=0`$:

$`0`$ $`=`$ $`\left(1\overline{g}_{aa}{\displaystyle \frac{\pi }{\beta }}{\displaystyle \underset{\omega _m}{}}s^0(\omega _m){\displaystyle \frac{\zeta _{aa}^0}{|\zeta _a^0\omega _m|}}\right)\left(1\overline{g}_{bb}{\displaystyle \frac{\pi }{\beta }}{\displaystyle \underset{\omega _m}{}}s^0(\omega _m){\displaystyle \frac{\zeta _{bb}^0}{|\zeta _b^0\omega _m|}}\right)`$ (77)
$``$ $`\overline{g}_{ab}\overline{g}_{ba}{\displaystyle \frac{\pi }{\beta }}{\displaystyle \underset{\omega _m}{}}s^0(\omega _m){\displaystyle \frac{\zeta _{ab}^0}{|\zeta _a^0\omega _m|}}{\displaystyle \frac{\pi }{\beta }}{\displaystyle \underset{\omega _m^{}}{}}s^0(\omega _m^{}){\displaystyle \frac{\zeta _{ba}^0}{|\zeta _b^0\omega _m^{}|}}.`$ (78)
Here $`\zeta _j^0`$ and $`\zeta _{jj^{}}^0`$ are given by eqs. (62) and (67) with $`\mathrm{\Delta }_a=\mathrm{\Delta }_b=0`$:

$`\{\begin{array}{c}\zeta _a^0=1\left[{\displaystyle \frac{u_b}{D_b^0}}\eta _b{\displaystyle \frac{u_a}{D_a^0}}{\displaystyle \frac{|\overline{v}_{ab}|^2}{1+\alpha _b^2}}\right]{\displaystyle \frac{1}{|\omega _m|+\mathrm{\Gamma }_b}},\hfill \\ \zeta _b^0=1\left[{\displaystyle \frac{u_a}{D_a^0}}\eta _a{\displaystyle \frac{u_b}{D_b^0}}{\displaystyle \frac{|\overline{v}_{ba}|^2}{1+\alpha _a^2}}\right]{\displaystyle \frac{1}{|\omega _m|+\mathrm{\Gamma }_a}},\hfill \end{array}`$ (81)
$`\{\begin{array}{c}\zeta _{aa}^0=1\left[{\displaystyle \frac{u_b}{D_b^0}}\eta _b{\displaystyle \frac{\overline{g}_{ba}}{\overline{g}_{aa}}}{\displaystyle \frac{u_a}{D_a^0}}{\displaystyle \frac{\overline{v}_{ab}^2}{1+\alpha _b^2}}\right]{\displaystyle \frac{1}{|\omega _m|+\mathrm{\Gamma }_b}},\hfill \\ \zeta _{ab}^0=1\left[{\displaystyle \frac{u_b}{D_b^0}}\eta _b{\displaystyle \frac{\overline{g}_{bb}}{\overline{g}_{ab}}}{\displaystyle \frac{u_a}{D_a^0}}{\displaystyle \frac{\overline{v}_{ab}^2}{1+\alpha _b^2}}\right]{\displaystyle \frac{1}{|\omega _m|+\mathrm{\Gamma }_b}},\hfill \\ \zeta _{bb}^0=1\left[{\displaystyle \frac{u_a}{D_a^0}}\eta _a{\displaystyle \frac{\overline{g}_{ab}}{\overline{g}_{bb}}}{\displaystyle \frac{u_b}{D_b^0}}{\displaystyle \frac{\overline{v}_{ba}^2}{1+\alpha _a^2}}\right]{\displaystyle \frac{1}{|\omega _m|+\mathrm{\Gamma }_a}},\hfill \\ \zeta _{ba}^0=1\left[{\displaystyle \frac{u_a}{D_a^0}}\eta _a{\displaystyle \frac{\overline{g}_{aa}}{\overline{g}_{ba}}}{\displaystyle \frac{u_b}{D_b^0}}{\displaystyle \frac{\overline{v}_{ba}^2}{1+\alpha _a^2}}\right]{\displaystyle \frac{1}{|\omega _m|+\mathrm{\Gamma }_a}}.\hfill \end{array}`$ (86)
In eq. (78), $`s^0(\omega _m)=F_\omega ^0/|F_\omega ^0|`$, where $`F_\omega ^0`$ is given by

$$F_\omega ^0=\left(1\frac{u_a}{D_a^0}\frac{\eta _a}{|\omega _m|+\mathrm{\Gamma }_a}\right)\left(1\frac{u_b}{D_b^0}\frac{\eta _b}{|\omega _m|+\mathrm{\Gamma }_b}\right)\frac{|\overline{v}_{ab}|^4}{(1+\alpha _a^2)(1+\alpha _b^2)}\frac{u_au_b}{D_a^0D_b^0}\frac{1}{(|\omega _m|+\mathrm{\Gamma }_a)(|\omega _m|+\mathrm{\Gamma }_b)}.$$
(87)
In obtaining eqs. (81)-(87), we have used the result that eq. (41) at $`T_\mathrm{c}`$ is reduced to

$`\{\begin{array}{c}\mathrm{i}\overline{\omega }_m^a=\mathrm{i}\omega _m+\mathrm{i}{\displaystyle \frac{u_a}{D_a^0}}\left[{\displaystyle \frac{|\overline{v}_{ab}|^2}{1+\alpha _b^2}}+\eta _a\right]{\displaystyle \frac{\omega _m}{|\omega _m|}}\mathrm{i}\omega _m+\mathrm{i}\mathrm{\Gamma }_a{\displaystyle \frac{\omega _m}{|\omega _m|}},\hfill \\ \mathrm{i}\overline{\omega }_m^b=\mathrm{i}\omega _m+\mathrm{i}{\displaystyle \frac{u_b}{D_b^0}}\left[{\displaystyle \frac{|\overline{v}_{ba}|^2}{1+\alpha _a^2}}+\eta _b\right]{\displaystyle \frac{\omega _m}{|\omega _m|}}\mathrm{i}\omega _m+\mathrm{i}\mathrm{\Gamma }_b{\displaystyle \frac{\omega _m}{|\omega _m|}},\hfill \end{array}`$ (90)
where

$`\{\begin{array}{c}D_a^0=\left(1+{\displaystyle \frac{|\overline{v}_{ab}|^2}{1+\alpha _b^2}}\right)^2+\left(\alpha _a{\displaystyle \frac{\alpha _b}{1+\alpha _b^2}}|\overline{v}_{ab}|^2\right)^2,\hfill \\ D_b^0=\left(1+{\displaystyle \frac{|\overline{v}_{ba}|^2}{1+\alpha _a^2}}\right)^2+\left(\alpha _b{\displaystyle \frac{\alpha _a}{1+\alpha _a^2}}|\overline{v}_{ba}|^2\right)^2.\hfill \end{array}`$ (93)
In the present case, effects of impurity scatterings are controlled by (i) the impurity potential ($`v_{aa}`$, $`v_{bb}`$, $`v_{ab}`$) and (ii) the impurity concentration described by $`u_j=n_{\mathrm{imp}}/(\pi N_j(0))`$. In Figs. 3 and 4, we show their effects on $`T_\mathrm{c}`$:

1. Figure 3 shows the effects of the interband scattering $`\overline{v}_{ab}`$ on the transition temperature $`T_\mathrm{c}`$. When $`\overline{v}_{ab}`$ becomes finite, $`T_\mathrm{c}`$ decreases with increasing $`|\overline{v}_{ab}|`$ in the region $`0<|\overline{v}_{ab}|<1`$; however, $`T_\mathrm{c}`$ again increases when $`|\overline{v}_{ab}|>1`$ reflecting that the interband scattering becomes irrelevant in the limit $`\overline{v}_{ab}\mathrm{}`$. When the interband pairing interaction $`\overline{g}_{ab}`$ is strong, the superconductong states in the two bands strongly couple with each other even in the clean system. In this case, even if transfers of Cooper-pairs occur between the two bands by interband scatterings, this effect would be weak when the phase of the interband scattering $`\theta _{ab}`$ is absent. Indeed as shown in Fig. 3(a), the suppression of $`T_\mathrm{c}`$ by $`\overline{v}_{ab}`$ is weaker for larger $`\overline{g}_{ab}`$. However, when the phase of $`\overline{v}_{ab}`$ is non-zero, the Cooper-pair acquires the phase $`2\times \theta _{ab}`$ in the transfer between the two bands. Then since the phases of the order parameters in the two bands are the same ($`\overline{g}_{bb}>0`$), the appearance of the Cooper-pair having this extra-phase would destroy superconductivity. Thus $`T_\mathrm{c}`$ is more suppressed when $`\theta _{ab}`$ is finite in Fig. 3(b). On the other hand, since the intraband scattering weakens the effect of the interband scattering, the suppression of $`T_\mathrm{c}`$ becomes weak when the intraband scattering becomes finite as shown in Fig. 3(c).
2. Effects of the impurity concentration on $`T_\mathrm{c}`$ is shown in Fig. 4. We find that the suppression of $`T_\mathrm{c}`$ continues to exist as the impurity concentration increases, which is in contrast to the case of the increase of $`|\overline{v}_{ab}|`$ in Fig. 3. In the case of $`\theta _{ab}=0`$ (Fig. 4(a)), the suppression of $`T_\mathrm{c}`$ is weaker when the interband pairing $`\overline{g}_{ab}`$ is stronger. On the other hand, the opposite tendency is obtained in the case of $`\theta _{ab}=\pi /2`$ as shown in Fig. 4(b). In this former case, when $`\overline{g}_{ab}`$ is strong, the superconducting states in the two bands strongly coupled with each other as noted in (i). Thus the mixing effect of two order parameters by interband scatterings is weaker for larger $`\overline{g}_{ab}`$ in the case of $`\theta _{ab}=0`$. However, in the case of Fig. 4(b) ($`\theta _{ab}0`$), since the Cooper-pair acquires the phase $`2\theta _{ab}`$ in the tunnelling by the interband scattering, the depairing effect becomes more remarkable when $`\overline{g}_{ab}`$ is stronger. (In the extreme case $`\overline{g}_{ab}=0`$, one can freely choose the phases of $`\mathrm{\Delta }_a`$ and $`\mathrm{\Delta }_b`$. Then the depairing effect by the tunnelling of the Cooper-pair by the interband scattering can be eliminated by choosing the difference of the phases of $`\mathrm{\Delta }_a`$ and $`\mathrm{\Delta }_b`$ as $`2\theta _{ab}`$. Indeed, one can eliminate $`\theta _{ab}`$ from the present model by an appropriate gauge transformation. This means that the depairing effect by $`\theta _{ab}`$ is absent when $`\overline{g}_{ab}=0`$.) As a result, in contrast to the case of $`\theta _{ab}=0`$, the suppression of $`T_\mathrm{c}`$ is more remarkable in Fig. 4(b) when the interband pairing $`\overline{g}_{ab}`$ is stronger.

### D Temperature dependence of $`\stackrel{~}{\mathrm{\Delta }}_a`$ and $`\stackrel{~}{\mathrm{\Delta }}_b`$

Figures 5-7 shows the temperature dependence of $`\stackrel{~}{\mathrm{\Delta }}_a`$ and $`\stackrel{~}{\mathrm{\Delta }}_b`$. In Fig. 5, when the impurity concentration ($`u_j`$) increases, $`|\stackrel{~}{\mathrm{\Delta }}_a|`$ and $`|\stackrel{~}{\mathrm{\Delta }}_b|`$ approach each other irrespective of $`\theta _{ab}`$. Although we do not show the result for $`0|\overline{v}_{ab}|<2`$ in Fig. 6, the same tendency as in Fig. 5 is obtained with increasing $`|\overline{v}_{ab}|`$ up to $`|\overline{v}_{ab}|<1`$. However, when $`|\overline{v}_{ab}|>1`$, Fig. 6 shows that difference between $`|\stackrel{~}{\mathrm{\Delta }}_a|`$ and $`|\stackrel{~}{\mathrm{\Delta }}_b|`$ becomes more remarkable for larger $`|\overline{v}_{ab}|`$ because of the same reason as that discussed in Secs. 4.2 and 4.3. This tendency is also obtained in Fig. 7 as $`\alpha _j`$ increases.

Figure 5(b) shows that $`\stackrel{~}{\mathrm{\Delta }}_b`$ becomes negative when $`u`$ becomes strong. This is because the Cooper-pair acquires the phase $`2\theta _{ab}(=\pi )`$ in the tunnelling by interband scatterings: When this tunnelling occurs frequently, it is favourable that the phases of the order parameters differ by $`\pi `$ in order to avoid the depairing effect caused by the tunnelling of the Cooper-pairs having the phase $`\pi `$. This ”$`\pi `$-phase” disappears when the interband scattering becomes very strong: In Fig. 6(b), $`\stackrel{~}{\mathrm{\Delta }}_b`$ is positive at $`\overline{v}_{ab}=10`$. In more detail, the stability of this $`\pi `$-phase is summarized in Fig. 8: The $`\pi `$-phase continues to exist when the impurity concentration is high ($`u_j>2`$) (Fig. 8(a)), while it disappears when the inter- and/or intra-band scatterings become strong ($`|\overline{v}_{ab}|>3`$, $`\alpha _j>1.5`$), as shown in panels (b) and (c).

### E Superconducting density of states

Effects of the impurity concentration and the phase/magnitude of the impurity scattering on the mixing of the order parameters have been clarified in Sec. 4.4. In this section, we study how the mixing effect is actually observed in the superconducting density of states (SDOS), because the tunnelling measurement is useful in observing this phenomenon.

The superconducting density of states in the $`j`$-band is given by

$$N_j(\omega )=N_j(0)\mathrm{Re}\left[\frac{Z_j(\omega )}{\sqrt{Z_j(\omega )^2\overline{\stackrel{~}{\mathrm{\Delta }}}_j^R(\omega )^2\overline{\stackrel{~}{\mathrm{\Delta }}}_j^I(\omega )^2}}\right]=\mathrm{Re}\left[\frac{1}{\sqrt{1\mathrm{\Xi }_j^2}}\right].$$
(94)
Here $`Z_j(\omega )=\mathrm{i}\overline{\omega }_m^j(\mathrm{i}\omega _m\omega +\mathrm{i}\delta )`$ is the analytic continued renormalized Matsubara frequency in eq. (41); $`\overline{\stackrel{~}{\mathrm{\Delta }}}_j^R(\omega )`$ and $`\overline{\stackrel{~}{\mathrm{\Delta }}}_j^I(\omega )`$ are, respectively, the analytic continued quantities of the real and imaginary part of $`\overline{\stackrel{~}{\mathrm{\Delta }}}_j`$ in eq. (44). Their expressions are given by ($`\overline{\stackrel{~}{\mathrm{\Delta }}}_j(\omega )^2\overline{\stackrel{~}{\mathrm{\Delta }}}_j^R(\omega )^2+\overline{\stackrel{~}{\mathrm{\Delta }}}_j^I(\omega )^2`$)

$`\{\begin{array}{c}Z_a(\omega )=\omega +\mathrm{i}\frac{u_a}{D_a(\omega )}\left[\frac{|\overline{v}_{ab}|^2}{1+\alpha _b^2}\frac{Z_b(\omega )}{\sqrt{Z_b(\omega )^2\overline{\stackrel{~}{\mathrm{\Delta }}}_b(\omega )^2}}+\eta _a\frac{Z_a(\omega )}{\sqrt{Z_a(\omega )^2\overline{\stackrel{~}{\mathrm{\Delta }}}_a(\omega )^2}}\right],\hfill \\ \overline{\stackrel{~}{\mathrm{\Delta }}}_a^R(\omega )=\mathrm{Re}[\stackrel{~}{\mathrm{\Delta }}_a]+\mathrm{i}\frac{u_a}{D_a(\omega )}\left[\frac{|\overline{v}_{ab}|^2}{1+\alpha _b^2}\frac{\overline{\stackrel{~}{\mathrm{\Delta }}}_b^R(\omega )\mathrm{cos}(2\theta _{ab})\overline{\stackrel{~}{\mathrm{\Delta }}}_b^I(\omega )\mathrm{sin}(2\theta _{ab})}{\sqrt{Z_b(\omega )^2\overline{\stackrel{~}{\mathrm{\Delta }}}_b(\omega )^2}}+\eta _a\frac{\overline{\stackrel{~}{\mathrm{\Delta }}}_a^R(\omega )}{\sqrt{Z_a(\omega )^2\overline{\stackrel{~}{\mathrm{\Delta }}}_a(\omega )^2}}\right],\hfill \\ \overline{\stackrel{~}{\mathrm{\Delta }}}_a^I(\omega )=\mathrm{Im}[\stackrel{~}{\mathrm{\Delta }}_a]+\mathrm{i}\frac{u_a}{D_a(\omega )}\left[\frac{|\overline{v}_{ab}|^2}{1+\alpha _b^2}\frac{\overline{\stackrel{~}{\mathrm{\Delta }}}_b^R(\omega )\mathrm{sin}(2\theta _{ab})+\overline{\stackrel{~}{\mathrm{\Delta }}}_b^I(\omega )\mathrm{cos}(2\theta _{ab})}{\sqrt{Z_b(\omega )^2\overline{\stackrel{~}{\mathrm{\Delta }}}_b(\omega )^2}}+\eta _a\frac{\overline{\stackrel{~}{\mathrm{\Delta }}}_a^I(\omega )}{\sqrt{Z_a(\omega )^2\overline{\stackrel{~}{\mathrm{\Delta }}}_a(\omega )^2}}\right],\hfill \end{array}`$ (98)
where

$$D_a(\omega )=1+\eta _a+\frac{2|\overline{v}_{ab}|^2}{1+\alpha _b^2}\frac{Z_a(\omega )Z_b(\omega )(\overline{\stackrel{~}{\mathrm{\Delta }}}_a^R\overline{\stackrel{~}{\mathrm{\Delta }}}_b^R+\overline{\stackrel{~}{\mathrm{\Delta }}}_a^I\overline{\stackrel{~}{\mathrm{\Delta }}}_b^I)\mathrm{cos}(2\theta _{ab})(\overline{\stackrel{~}{\mathrm{\Delta }}}_b^R\overline{\stackrel{~}{\mathrm{\Delta }}}_a^I\overline{\stackrel{~}{\mathrm{\Delta }}}_b^I\overline{\stackrel{~}{\mathrm{\Delta }}}_a^R)\mathrm{sin}(2\theta _{ab})}{\sqrt{Z_a(\omega )^2\overline{\stackrel{~}{\mathrm{\Delta }}}_a(\omega )^2}\sqrt{Z_b(\omega )^2\overline{\stackrel{~}{\mathrm{\Delta }}}_b(\omega )^2}}.$$
(99)
The valuables in the $`b`$-band ($`Z_b(\omega )`$, $`\overline{\stackrel{~}{\mathrm{\Delta }}}_b^R(\omega )`$, $`\overline{\stackrel{~}{\mathrm{\Delta }}}_b^I(\omega )`$ and $`D_b(\omega )`$) are also given by eqs. (98) and (99) where the band indices ’$`a`$’ and ’$`b`$’ are interchanged. In calculating SDOS, it is convenient to use the last expression in eq. (94), where $`\mathrm{\Xi }_j=\sqrt{\overline{\stackrel{~}{\mathrm{\Delta }}}_j^R(\omega )^2+\overline{\stackrel{~}{\mathrm{\Delta }}}_j^I(\omega )^2}/Z_j(\omega )\sqrt{\overline{\mathrm{\Xi }}_j^R(\omega )^2+\overline{\mathrm{\Xi }}_j^I(\omega )^2}`$:

$`\{\begin{array}{c}\overline{\mathrm{\Xi }}_a^R={\displaystyle \frac{\mathrm{Re}[\stackrel{~}{\mathrm{\Delta }}_a]+\mathrm{i}\frac{u_a}{D_a(\omega )}\frac{|\overline{v}_{ab}|^2}{1+\alpha _b^2}\frac{\overline{\mathrm{\Xi }}_b^R\mathrm{cos}(2\theta _{ab})\overline{\mathrm{\Xi }}_b^I\mathrm{sin}(2\theta _{ab})}{\sqrt{1\mathrm{\Xi }_b^2}}}{\omega +\mathrm{i}\frac{u_a}{D_a(\omega )}\frac{|\overline{v}_{ab}|^2}{1+\alpha _b^2}\frac{1}{\sqrt{1\mathrm{\Xi }_b^2}}}},\hfill \\ \overline{\mathrm{\Xi }}_a^I={\displaystyle \frac{\mathrm{Im}[\stackrel{~}{\mathrm{\Delta }}_a]+\mathrm{i}\frac{u_a}{D_a(\omega )}\frac{|\overline{v}_{ab}|^2}{1+\alpha _b^2}\frac{\overline{\mathrm{\Xi }}_b^R\mathrm{sin}(2\theta _{ab})+\overline{\mathrm{\Xi }}_b^I\mathrm{cos}(2\theta _{ab})}{\sqrt{1\mathrm{\Xi }_b^2}}}{\omega +\mathrm{i}\frac{u_a}{D_a(\omega )}\frac{|\overline{v}_{ab}|^2}{1+\alpha _b^2}\frac{1}{\sqrt{1\mathrm{\Xi }_b^2}}}}.\hfill \end{array}`$ (102)
$`\mathrm{\Xi }_b^R`$ and $`\mathrm{\Xi }_b^I`$ are obtained from eq. (102) by interchanging the indices ’$`a`$’ and ’$`b`$’.

Figures 9 and 10 show impurity effects on SDOS in the cases of $`\theta _{ab}=0`$ and $`\theta _{ab}=\pi /2`$, respectively. When the impurity concentration increases, SDOS’s in the two bands are mixed with each other to be almost the same at $`u=20`$ in Fig. 9(a) and $`u=5`$ in Fig 10(a). In Fig. 10(a), because of the strong depairing effect by the phase of the interband scattering $`\theta _{ab}`$, the mixed SDOS at $`u=5`$ has a small energy gap compared with the case of $`u=0.1`$. This strong depairing effect also leads to a gapless behavior at $`u=1`$ in Fig. 10(a).

The same mixing effect is observed also in Figs 9(b) and 10(b): Although SDOS’s in the two bands are different at $`|\overline{v}_{ab}|=0.1`$, they become almost the same at $`|\overline{v}_{ab}|=1`$; however, since the interband impurity effect again becomes weak when $`|\overline{v}_{ab}|>1`$, SDOS’s differ from each other at $`|\overline{v}_{ab}|=100`$. Since the effect of the interband scattering is also weakened by the intraband scattering, we obtain two distinct SDOS’s at $`\alpha =10`$ in Figs. 9(c) and 10(c).

If MgB2 is really a two-band superconductor with two distinct order parameters, we expect the mixing effect in SDOS shown in Figs. 9(a) and 10(a) as the impurity concentration is increased. This effect is independent of the phase of the interband scattering as shown in Figs. 9 and 10 and, in principle, it always occurs when the impurity concentration is high enough. This effect might be useful to judge whether or not superconductivity in MgB2 has two distinct order parameters in the two bands. However, infinitely high impurity concentration is impossible experimentally. Thus when we try to observe the mixing effect, the character of impurity potential is crucial: The intraband scattering must not be strong, and the strength of the interband scattering must be moderate ($`|\overline{v}_{ab}|1`$).

## V Summary

In this paper, we have investigated nonmagnetic impurity effects on two-band $`s`$-wave superconductivity with two different order parameters. We showed that a bound state appears around a nonmagnetic impurity depending on the character of the interband scattering potential. We also considered the case of finite impurity concentration focusing on the mixing of the two order parameters by interband scatterings. In contrast to the previous works which treat the impurity scattering within the Born approximation, we took into account multi-scattering processes using the $`t`$-matrix approximation. We showed that the mixing of the two order parameters occurs when the impurity concentration becomes high, as discussed previously. However, although the interband scattering is essential for the mixing effect, we clarified that this effect disappears when the interband scattering is very strong ($`|\overline{v}_{ab}|1`$.). In addition, the intraband scattering also suppresses the mixing effect. When the magnitude of the impurity potential is very large, a clear two-gap structure is observed in the superconducting density of states as in the case of the clean system (unless the impurity concentration is very high.).

In order to observe the mixing effect of the order parameters in a two-band superconductor, we have to choose impurities with weak intraband scattering and moderate interband scattering ($`|\overline{v}_{ab}|1`$.). If one can choose this kind of impurity for MgB2, it is expected that, if the observed two-gap structure in the superconducting density of state originates from two order parameters in the two bands, the two-gap structure would gradually change into a single-gap structure as the impurity concentration increases. The observation of this impurity effect may be useful in clarifying the origin of the two-gap structure in the superconducting density of states in MgB2.

###### Acknowledgements.

The author would like to thank Professor S. Takada for fruitful discussions on the absence of the bound state in the strong interband scattering limit. He also thanks Professor A. Griffin for his hospitality during his stay in University of Toronto. He was financially supported by a Japanese Overseas Research Fellowship.

Appendices

## A Equation of Bound State Energy

From eq. (28), the energy of the bound state $`\omega `$ is determined by the equation

$`0`$ $`=`$ $`\mathrm{det}[1\widehat{V}_{aa}(\mathrm{i}\omega _m){\displaystyle \underset{\text{p}}{}}G_a^0(p,\mathrm{i}\omega _m)]`$ (A1)
$`=`$ $`\mathrm{det}[1\widehat{v}_{ab}\widehat{P}_b\widehat{v}_{ba}\widehat{P}_a]\left]\mathrm{det}\right[1\widehat{v}_{aa}{\displaystyle \underset{\text{p}}{}}G_a^0(p,\mathrm{i}\omega _m)],`$ (A2)
where $`\widehat{P}_j=_\text{p}G_j^0(p,\mathrm{i}\omega _m)[1\widehat{v}_{jj}_\text{p}G_j^0(p,\mathrm{i}\omega _m)]^1`$ ($`j=a,b`$). Since $`\widehat{v}_{aa}`$ is the intraband nonmagnetic impurity potential, no pole is obtained from the factor $`\mathrm{det}[1\widehat{v}_{aa}_\text{p}G_a^0(p,\mathrm{i}\omega _m)]`$ in eq. (A2). Hence we drop this term in what follows. The factor $`\widehat{P}_j`$ can be calculated using

$$\underset{\text{p}}{}G_j^0=\pi N_j(0)\frac{\mathrm{i}\omega _m\stackrel{~}{\mathrm{\Delta }}_j\rho _1}{\mathrm{\Omega }_j}.$$
(A3)
Here $`\mathrm{\Omega }_j=\sqrt{\omega _m^2+\stackrel{~}{\mathrm{\Delta }}_j^2}`$ and we have taken the order parameters real as $`\widehat{\mathrm{\Delta }}_j=\stackrel{~}{\mathrm{\Delta }}_j\rho _1`$. Then we obtain

$$\widehat{P}_j=\frac{\pi N_j(0)}{1+\alpha _j^2}\left[\alpha _j\rho _3+\frac{\mathrm{i}\omega _m\stackrel{~}{\mathrm{\Delta }}_j\rho _1}{\mathrm{\Omega }_j}\right].$$
(A4)
Substituting eqs.(A4) and (18) into eq. (A2), we obtain the equation of the pole as

$`0=\mathrm{det}\left(\begin{array}{cc}1A& B\\ B^{}& 1A^{}\end{array}\right)=1+|A|^2+|B|^2(A+A^{}),`$ (A7)
where the matrix elements $`A`$ and $`b`$ are given by

$`\{\begin{array}{c}A={\displaystyle \frac{1}{(1+\alpha _a^2)(1+\alpha _b^2)}}\left[|\overline{v}_{ab}|^2({\displaystyle \frac{\mathrm{i}\omega _m}{\mathrm{\Omega }_a}}+\alpha _a)({\displaystyle \frac{\mathrm{i}\omega _m}{\mathrm{\Omega }_b}}+\alpha _b)\overline{v}_{ab}^2{\displaystyle \frac{\stackrel{~}{\mathrm{\Delta }}_a}{\mathrm{\Omega }_a}}{\displaystyle \frac{\stackrel{~}{\mathrm{\Delta }}_b}{\mathrm{\Omega }_b}}\right],\hfill \\ B={\displaystyle \frac{1}{(1+\alpha _a^2)(1+\alpha _b^2)}}\left[|\overline{v}_{ab}|^2{\displaystyle \frac{\stackrel{~}{\mathrm{\Delta }}_a}{\mathrm{\Omega }_a}}({\displaystyle \frac{\mathrm{i}\omega _m}{\mathrm{\Omega }_b}}+\alpha _b)+\overline{v}_{ab}^2{\displaystyle \frac{\stackrel{~}{\mathrm{\Delta }}_b}{\mathrm{\Omega }_b}}({\displaystyle \frac{\mathrm{i}\omega _m}{\mathrm{\Omega }_a}}\alpha _a)\right].\hfill \end{array}`$ (A10)
Calculating each term in eq. (A7) and taking the analytic continuation $`\mathrm{i}\omega _m\omega +i\delta `$, we obtain eq. (30).

## B Derivation of the Renormalized Variables $`\overline{\omega }_m^j`$ and $`\overline{\stackrel{~}{\mathrm{\Delta }}}_j`$

We consider the $`a`$-band. The self-energy $`\mathrm{\Sigma }_{aa}`$ in eq. (35) can be rewritten as

$$\mathrm{\Sigma }_{aa}(\mathrm{i}\omega _m)=u_a\frac{1}{\frac{1}{\pi N_a(0)\overline{V}_{aa}}\frac{1}{\pi N_a(0)}_\text{p}G_a(p,\mathrm{i}\omega _m)}.$$
(B1)
The summation of the Green function in terms of $`p`$ gives

$`{\displaystyle \frac{1}{\pi N_a(0)}}{\displaystyle \underset{\text{p}}{}}G_a(p,\mathrm{i}\omega _m)={\displaystyle \frac{\mathrm{i}\overline{\omega }_m^a\widehat{\overline{\mathrm{\Delta }}}_a}{\overline{\mathrm{\Omega }}_a}}`$ $`=`$ $`{\displaystyle \frac{\mathrm{i}\overline{\omega }_m^a}{\overline{\mathrm{\Omega }}_a}}+{\displaystyle \frac{|\overline{\stackrel{~}{\mathrm{\Delta }}}_a|}{\overline{\mathrm{\Omega }}_a}}\mathrm{cos}(\overline{\stackrel{~}{\varphi }}_a)\rho _1{\displaystyle \frac{|\overline{\stackrel{~}{\mathrm{\Delta }}}_a|}{\overline{\mathrm{\Omega }}_a}}\mathrm{sin}(\overline{\stackrel{~}{\varphi }}_a)\rho _2`$ (B2)
$``$ $`\gamma _0+\gamma _1\rho _1+\gamma _2\rho _2,`$ (B3)
where $`\overline{\stackrel{~}{\varphi }}_j`$ is the phase of $`\overline{\stackrel{~}{\mathrm{\Delta }}}_j`$ ($`j=a,b`$).

Using the expression of $`\overline{V}_{aa}`$ in eq. (35) and eq. (B3), we obtain

$`\pi N_a(0)\overline{V}_{aa}`$ $`=`$ $`{\displaystyle \frac{|\overline{v}_{ab}|^2}{1+\alpha _b^2}}{\displaystyle \frac{\mathrm{i}\omega _m}{\overline{\mathrm{\Omega }}_b}}{\displaystyle \frac{|\overline{v}_{ab}|^2}{1+\alpha _b^2}}{\displaystyle \frac{|\overline{\stackrel{~}{\mathrm{\Delta }}}_b|}{\overline{\mathrm{\Omega }}_b}}\mathrm{cos}(2\theta _{ab}+\overline{\stackrel{~}{\varphi }}_b)\rho _1`$ (B4)
$`+`$ $`{\displaystyle \frac{|\overline{v}_{ab}|^2}{1+\alpha _b^2}}{\displaystyle \frac{|\overline{\stackrel{~}{\mathrm{\Delta }}}_b|}{\overline{\mathrm{\Omega }}_b}}\mathrm{sin}(2\theta _{ab}+\overline{\stackrel{~}{\varphi }}_b)\rho _2+\left(\alpha _a{\displaystyle \frac{\alpha _b}{1+\alpha _b^2}}|\overline{v}_{ab}|^2\right)\rho _3`$ (B5)
$``$ $`\lambda _0+\lambda _1\rho _1+\lambda _2\rho _2+\lambda _3\rho _3.`$ (B6)
Thus

$$\frac{1}{\pi N_a(0)\overline{V}_{aa}}=\frac{\lambda _1\rho _1+\lambda _2\rho _2+\lambda _3\rho _3\lambda _0}{\eta _a},$$
(B7)
where $`\eta _a=\lambda _1^2+\lambda _2^2+\lambda _3^2\lambda _0^2`$ is reduced to eq. (47). Substituting eqs. (B3) and (B7) into eq. (B1), we have

$$\mathrm{\Sigma }_{aa}(\mathrm{i}\omega _m)=\frac{u_a}{D_a}[(\lambda _0+\eta _a\gamma _0)+(\lambda _1\eta _a\gamma _1)\rho _1+(\lambda _2\eta _a\gamma _2)\rho _2+\lambda _3\rho _3].$$
(B8)
Here $`D_a=[(\lambda _1\eta _a\gamma _1)^2+(\lambda _2\eta _a\gamma _2)^2+\lambda _3^2(\lambda _0+\eta _a\gamma _0)^2]/\eta _a`$ equals eq. (50). Thus from eqs. (32) and (B8), we find

$`\{\begin{array}{c}\mathrm{i}\overline{\omega }_m=\mathrm{i}\omega _m\frac{u_a}{D_a}(\lambda _0+\eta _a\gamma _0),\hfill \\ |\overline{\stackrel{~}{\mathrm{\Delta }}}_a|\mathrm{cos}(\overline{\stackrel{~}{\varphi }}_a)=|\stackrel{~}{\mathrm{\Delta }}_a|\mathrm{cos}(\stackrel{~}{\varphi }_a)\frac{u_a}{D_a}(\lambda _1\eta _a\gamma _1),\hfill \\ |\overline{\stackrel{~}{\mathrm{\Delta }}}_a|\mathrm{sin}(\overline{\stackrel{~}{\varphi }}_a)=|\stackrel{~}{\mathrm{\Delta }}_a|\mathrm{sin}(\stackrel{~}{\varphi }_a)+\frac{u_a}{D_a}(\lambda _2\eta _a\gamma _2),\hfill \end{array}`$ (B12)
where $`\stackrel{~}{\varphi }_j`$ is the phase of $`\stackrel{~}{\mathrm{\Delta }}_j`$. Substituting the detailed expressions of $`\gamma _i`$ and $`\lambda _i`$ into eq. (B12), we obtain eqs. (41) and (44).

We briefly mention that the term proportional to $`\rho _3`$ in eq. (B8) gives a renormalization of the kinetic energy as

$$\stackrel{~}{\epsilon }_\text{p}^a=\epsilon _\text{p}^a+u_a\frac{\lambda _3}{D_a}.$$
(B13)
This renormalization may be crucial when we consider a correlation function; however, since it does not affect the integration in terms of the kinetic energy within the treatment in this paper, we can safely neglect it.

## C Mixing of the order parameters within the Born approximation

We consider the case of $`\alpha _a=\alpha _b=\overline{g}_{ab}=\overline{g}_{ba}=0`$. In this case, since the phase of the interband scattering $`\theta _{ab}`$ is irrelevant, we can take $`v_{ab}`$ real. In addition, $`\stackrel{~}{\mathrm{\Delta }}_j=\mathrm{\Delta }_j`$ and $`\mathrm{\Delta }_b`$ can be also taken real. When we retain scattering processes within the Born approximation, eqs. (41) and (44) are reduced to

$`\{\begin{array}{c}\mathrm{i}\overline{\omega }_m^a=\mathrm{i}\omega _m+\mathrm{\Gamma }_a{\displaystyle \frac{\mathrm{i}\overline{\omega }_m^b}{\overline{\mathrm{\Omega }}_b}},\hfill \\ \mathrm{i}\overline{\omega }_m^b=\mathrm{i}\omega _m+\mathrm{\Gamma }_b{\displaystyle \frac{\mathrm{i}\overline{\omega }_m^a}{\overline{\mathrm{\Omega }}_a}},\hfill \end{array}`$ (C3)
$`\{\begin{array}{c}\overline{\stackrel{~}{\mathrm{\Delta }}}_a=\mathrm{\Delta }_a+\mathrm{\Gamma }_a{\displaystyle \frac{\overline{\stackrel{~}{\mathrm{\Delta }}}_b}{\overline{\mathrm{\Omega }}_b}},\hfill \\ \overline{\stackrel{~}{\mathrm{\Delta }}}_b=\mathrm{\Delta }_b+\mathrm{\Gamma }_b{\displaystyle \frac{\overline{\stackrel{~}{\mathrm{\Delta }}}_a}{\overline{\mathrm{\Omega }}_a}},\hfill \end{array}`$ (C6)
where $`\mathrm{\Gamma }_a`$ and $`\mathrm{\Gamma }_b`$ are defined in eq. (90) where we take $`\alpha _a=\alpha _b=\eta _a=\eta _b=0`$ and $`D_a^0=D_b^0=1`$. Equations (C3) and (C6) can be rewritten as

$`\{\begin{array}{c}\overline{\omega }_r^a=1+\mathrm{\Gamma }_a{\displaystyle \frac{\overline{\omega }_r^b}{\overline{\mathrm{\Omega }}_b}},\hfill \\ \overline{\omega }_r^b=1+\mathrm{\Gamma }_b{\displaystyle \frac{\overline{\omega }_r^a}{\overline{\mathrm{\Omega }}_a}},\hfill \end{array}`$ (C9)
$`\{\begin{array}{c}\overline{\mathrm{\Delta }}_r^a=1+X\mathrm{\Gamma }_a{\displaystyle \frac{\overline{\mathrm{\Delta }}_r^b}{\overline{\mathrm{\Omega }}_b}},\hfill \\ \stackrel{~}{\mathrm{\Delta }}_r^b=1+{\displaystyle \frac{\mathrm{\Gamma }_b}{X}}{\displaystyle \frac{\overline{\mathrm{\Delta }}_r^a}{\overline{\mathrm{\Omega }}_a}}.\hfill \end{array}`$ (C12)
Here we have introduced $`\overline{\omega }_r^j\overline{\omega }_m^j/\omega _m`$, $`\overline{\mathrm{\Delta }}_r^j\overline{\stackrel{~}{\mathrm{\Delta }}}_j/\mathrm{\Delta }_j`$ and $`X\mathrm{\Delta }_b/\mathrm{\Delta }_a`$. When we define $`Y_j\overline{\mathrm{\Delta }}_r^j/\overline{\omega }_r^j`$, we find a relation $`Y_a=XY_b`$ in the dirty limit ($`\mathrm{\Gamma }_a1`$ and $`\mathrm{\Gamma }_b1`$). Using this relation, we find that $`Y_a`$ and $`Y_b`$ obey the following equation:

$`\left(\begin{array}{cc}1+\frac{\mathrm{\Gamma }_a}{\mathrm{\Omega }_r}& \frac{\mathrm{\Gamma }_a}{X\mathrm{\Omega }_r}\\ \frac{\mathrm{\Gamma }_b}{X\mathrm{\Omega }_r}& 1+\frac{\mathrm{\Gamma }_b}{\mathrm{\Omega }_r}\end{array}\right)\left(\begin{array}{c}Y_a\\ Y_b\end{array}\right)=\left(\begin{array}{c}1\\ 1\end{array}\right),`$ (C19)
where $`\mathrm{\Omega }_r=\sqrt{\omega _m^2+Y_b^2\stackrel{~}{\mathrm{\Delta }}_b^2}`$. The solution of eq. (C19) is given by $`Y_j=\mathrm{\Delta }_{AV}/\mathrm{\Delta }_j`$ (in the dirty limit), where

$$\mathrm{\Delta }_{AV}=\frac{\mathrm{\Delta }_aN_a(0)+\mathrm{\Delta }_bN_b(0)}{N_a(0)+N_b(0)}$$
(C20)
is an averaged order parameter. Then the gap equation (51) is reduced to

$$1=\overline{g}_{jj}Y_j\frac{\pi }{\beta }\underset{\omega _m}{}\frac{1}{\sqrt{\omega _m^2+\mathrm{\Delta }_{AV}^2}}.$$
(C21)
Substituting the explicit form of $`Y_j`$ into eq. (C21) we obtain

$`\left(\begin{array}{cc}1{\displaystyle \frac{\overline{g}_{aa}N_a(0)}{N_a(0)+N_b(0)}}{\displaystyle \frac{\pi }{\beta }}{\displaystyle \underset{\omega _m}{}}{\displaystyle \frac{1}{\sqrt{\omega _m^2+\mathrm{\Delta }_{AV}^2}}}& {\displaystyle \frac{\overline{g}_{aa}N_b(0)}{N_a(0)+N_b(0)}}{\displaystyle \frac{\pi }{\beta }}{\displaystyle \underset{\omega _m}{}}{\displaystyle \frac{1}{\sqrt{\omega _m^2+\mathrm{\Delta }_{AV}^2}}}\\ {\displaystyle \frac{\overline{g}_{bb}N_a(0)}{N_a(0)+N_b(0)}}{\displaystyle \frac{\pi }{\beta }}{\displaystyle \underset{\omega _m}{}}{\displaystyle \frac{1}{\sqrt{\omega _m^2+\mathrm{\Delta }_{AV}^2}}}& 1{\displaystyle \frac{\overline{g}_{bb}N_b(0)}{N_a(0)+N_b(0)}}{\displaystyle \frac{\pi }{\beta }}{\displaystyle \underset{\omega _m}{}}{\displaystyle \frac{1}{\sqrt{\omega _m^2+\mathrm{\Delta }_{AV}^2}}}\end{array}\right)\left(\begin{array}{c}\mathrm{\Delta }_a\\ \mathrm{\Delta }_b\end{array}\right)=0.`$ (C26)
(C27)
Equation (C27) has a solution when the following equation is satisfied:

$$1=\frac{\overline{g}_{aa}N_a(0)+\overline{g}_{bb}N_b(0)}{N_a(0)+N_b(0)}\frac{\pi }{\beta }\underset{\omega _m}{}\frac{1}{\sqrt{\omega _m^2+\mathrm{\Delta }_{AV}^2}}.$$
(C28)
Equation (C28) can be interpreted as the gap equation of $`\mathrm{\Delta }_{AV}`$ with the averaged pairing interaction $`(\overline{g}_{aa}N_a(0)+\overline{g}_{bb}N_b(0))/(N_a(0)+N_b(0))`$.

In the dirty limit, we obtain the following relation:

$$\frac{\overline{\omega }_m^j}{\sqrt{\overline{\omega }_m^j{}_{}{}^{2}+\overline{\stackrel{~}{\mathrm{\Delta }}}_j^2}}=\frac{\omega _m}{\sqrt{\omega _m^2+\mathrm{\Delta }_{AV}^2}}(\omega _m>0).$$
(C29)
Since the superconducting density of states in the $`j`$-band is obtained from the analytic continuation of LHS in eq. (C29), we find that the excitation gap in both bands is equal to $`\mathrm{\Delta }_{AV}`$. Namely, in the dirty limit, only one excitation gap $`\mathrm{\Delta }_{AV}`$ appears in the superconducting density of states within the Born approximation.

## D Evaluation of Eq. (75)

When we expand eq. (75), we obtain

$`1`$ $`=`$ $`\overline{g}_{aa}{\displaystyle \frac{\pi }{\beta }}{\displaystyle \underset{\omega _m}{}}s(\omega _m){\displaystyle \frac{\zeta _{aa}}{\overline{\mathrm{\Omega }}_a^{}}}+\overline{g}_{bb}{\displaystyle \frac{\pi }{\beta }}{\displaystyle \underset{\omega _m}{}}s(\omega _m){\displaystyle \frac{\zeta _{bb}}{\overline{\mathrm{\Omega }}_b^{}}}`$ (D1)
$``$ $`\overline{g}_{aa}\overline{g}_{bb}{\displaystyle \frac{\pi }{\beta }}{\displaystyle \underset{\omega _m}{}}s(\omega _m){\displaystyle \frac{\zeta _{aa}}{\overline{\mathrm{\Omega }}_a^{}}}{\displaystyle \frac{\pi }{\beta }}{\displaystyle \underset{\omega _m^{}}{}}s(\omega _m^{}){\displaystyle \frac{\zeta _{bb}}{\overline{\mathrm{\Omega }}_b^{}}}\overline{g}_{ab}\overline{g}_{ba}{\displaystyle \frac{\pi }{\beta }}{\displaystyle \underset{\omega _m}{}}s(\omega _m){\displaystyle \frac{\zeta _{ab}}{\overline{\mathrm{\Omega }}_a^{}}}{\displaystyle \frac{\pi }{\beta }}{\displaystyle \underset{\omega _m^{}}{}}s(\omega _m^{}){\displaystyle \frac{\zeta _{ba}}{\overline{\mathrm{\Omega }}_b^{}}}.`$ (D2)
The first and the second terms in RHS in eq. (D2) include the following complex terms:

$$\frac{\pi }{\beta }\underset{\omega _m}{}s(\omega _m)\left[\overline{g}_{ba}\frac{u_a}{D_a}\frac{\overline{v}_{ab}^2}{1+\alpha _b^2}\frac{1}{\overline{\mathrm{\Omega }}_b\overline{\mathrm{\Omega }}_a^{}}+\overline{g}_{ab}\frac{u_b}{D_b}\frac{\overline{v}_{ba}^2}{1+\alpha _a^2}\frac{1}{\overline{\mathrm{\Omega }}_a\overline{\mathrm{\Omega }}_b^{}}\right].$$
(D3)
However, from eqs. (50) and (55) and the definition of $`u_j`$, we obtain $`(1+\alpha _b^2)D_a=(1+\alpha _a^2)D_b`$ and $`N_a(0)u_a=N_b(0)u_b`$. Substituting these relations into eq. (D3) and noting that $`g_{ab}=g_{ba}`$ and $`\overline{\mathrm{\Omega }}_j=\overline{\mathrm{\Omega }}_j^{}/|F_\omega |`$, we find that eq. (D3) is real. Thus the first and the second terms in eq. (D2) is also real.

The third and the fourth terms in eq. (D2) can be rewritten as

$`(\overline{g}_{aa}\overline{g}_{bb}`$ $``$ $`\overline{g}_{ab}\overline{g}_{ba}){\displaystyle \frac{\pi ^2}{\beta ^2}}{\displaystyle }_{\omega _m\omega _m^{}}{\displaystyle \frac{s(\omega _m)s(\omega _m^{})}{\overline{\mathrm{\Omega }}_a^{}(\omega _m)\overline{\mathrm{\Omega }}_b^{}(\omega _m^{})}}[(1{\displaystyle \frac{u_b}{D_b(\omega _m)}}{\displaystyle \frac{\eta _b}{\overline{\mathrm{\Omega }}_b(\omega _m)}})(1{\displaystyle \frac{u_a}{D_a(\omega _m^{})}}{\displaystyle \frac{\eta _a}{\overline{\mathrm{\Omega }}_a(\omega _m^{})}})`$ (D4)
$``$ $`{\displaystyle \frac{u_au_b|\overline{v}_{ab}|^4}{(1+\alpha _a^2)(1+\alpha _b^2)D_a(\omega _m)D_b(\omega _m^{})\overline{\mathrm{\Omega }}_a(\omega _m)\overline{\mathrm{\Omega }}_b(\omega _m^{})}}].`$ (D5)
Equation (D5) is real, so that eq. (75) is found real.

Figure Captions

1. Temperature dependence of the order parameters in the clean system. We put $`N_a(0)=N_b(0)`$. This condition is used throughout this paper. In this case, $`\overline{g}_{ab}=\overline{g}_{ba}`$ is satisfied. In panels (a) and (b), the intraband interactions $`\overline{g}_{aa}`$ and $`\overline{g}_{bb}`$ are attractive; panel (a) shows the case of $`\overline{g}_{ab}(\overline{g}_{aa},\overline{g}_{bb})`$, while panel (b) is the case of $`\overline{g}_{ab}(\overline{g}_{aa},\overline{g}_{bb})`$. Panel (c) shows the case when $`\overline{g}_{bb}`$ is repulsive.
2. Dependence of the bound state energy on the interband impurity scattering $`\overline{v}_{ab}`$. We take $`\stackrel{~}{\mathrm{\Delta }}_a=2`$, $`\stackrel{~}{\mathrm{\Delta }}_b=1`$ and $`\alpha _a=0`$. Panel (a) shows the effect of the phase of $`\overline{v}_{ab}`$ in the case of $`\alpha _b=0`$. Since $`\stackrel{~}{\mathrm{\Delta }}_b^2\mathrm{cos}(2\theta _{ab})\stackrel{~}{\mathrm{\Delta }}_a\stackrel{~}{\mathrm{\Delta }}_b0`$ for $`\theta _{ab}\pi /6`$, the bound state is absent for $`\theta _{ab}\pi /6`$. Panel (b) shows the effect of the intraband potential scattering $`\alpha _b`$ in the $`b`$-band in the case of $`\theta _{ab}=\pi /2`$.
3. Dependence of the transition temperature $`T_\mathrm{c}`$ on the interband potential scattering $`\overline{v}_{ab}`$. The transition temperature is normalized by the value at $`\overline{v}_{ab}=0`$. We take $`\overline{g}_{aa}=0.3`$, $`\overline{g}_{bb}=0.25`$, $`u_a=u_b=1`$ and $`\alpha _a=0`$. (a) Effect of the interband pairing interaction $`\overline{g}_{ab}`$ in the case of $`\theta _{ab}=0`$ and $`\alpha _b=0`$. (b) Effect of the phase of the interband potential scattering $`\theta _{ab}`$ in the case of $`\overline{g}_{ab}=0.002`$ and $`\alpha _b=0`$. (b) Effect of the intraband potential scattering in the $`b`$-band ($`\alpha _b`$). We put $`\overline{g}_{ab}=0.002`$ and $`\theta _{ab}=0`$.
4. Dependence of the transition temperature $`T_\mathrm{c}`$ on the impurity concentration $`u_j=n_{\mathrm{imp}}/(\pi N_j(0))`$. In this figure, we take $`u_a=u_b`$ and $`T_\mathrm{c}`$ is normalized by the value in clean system ($`T_\mathrm{c}^0`$). The other parameters are taken as $`\overline{g}_{aa}=0.3`$, $`\overline{g}_{bb}=0.25`$, $`\alpha _a=\alpha _b=0`$ and $`|\overline{v}_{ab}|=1`$. In panel (a), the phase of the interband scattering is taken $`\theta _{ab}=0`$, while in panel (b) $`\theta _{ab}=\pi /2`$.
5. Temperature dependence of $`\stackrel{~}{\mathrm{\Delta }}_a`$ and $`\stackrel{~}{\mathrm{\Delta }}_b`$ in the cases of (a) $`\theta _{ab}=0`$ and (b) $`\theta _{ab}=\pi /2`$. We take $`\overline{g}_{aa}=0.3`$, $`\overline{g}_{bb}=0.25`$, $`\overline{g}_{ab}=0.002`$, $`\alpha _a=\alpha _b=0`$ and $`|\overline{v}_{ab}|=1`$. The $`x`$\- and $`y`$-axis are normalized by $`T_\mathrm{c}`$ in the clean system.
6. Effect of the interband scattering on temperature dependence of $`\stackrel{~}{\mathrm{\Delta }}_a`$ and $`\stackrel{~}{\mathrm{\Delta }}_b`$ in the cases of (a) $`\theta _{ab}=0`$ and (b) $`\theta _{ab}=\pi /2`$. We take $`\overline{g}_{aa}=0.3`$, $`\overline{g}_{bb}=0.25`$, $`\overline{g}_{ab}=0.002`$, $`\alpha _a=\alpha _b=0`$ and $`u_a=u_b=1`$.
7. Effect of intraband scattering ($`\alpha _a`$ and $`\alpha _b`$) on $`T_\mathrm{c}`$. The figure shows the case of $`\theta _{ab}=0`$ ($`\overline{v}_{ab}=1`$). The other parameters are $`\overline{g}_{aa}=0.3`$, $`\overline{g}_{bb}=0.25`$, $`\overline{g}_{ab}=0.002`$, $`u_a=u_b=10`$ and $`|\overline{v}_{ab}|=1`$.
8. Stability of ”$`\pi `$-phase” caused by the interband scattering with $`\theta _{ab}=\pi /2`$. We take $`T=0.5T_\mathrm{c}`$, $`\overline{g}_{aa}=0.3`$, $`\overline{g}_{bb}=0.25`$ and $`\overline{g}_{ab}=0.1`$.
9. Superconducting density of states at $`T=0.1T_\mathrm{c}`$ in the case of $`\theta _{ab}=0`$. We take $`\overline{g}_{aa}=0.3`$, $`\overline{g}_{bb}=0.25`$ and $`\overline{g}_{ab}=0.002`$. $`N_j(\omega )`$ is normalized by the density of states at the Fermi level in the normal state ($`N_a(0)=N_b(0)N(0)`$). In this figure, ’a’ and ’b’ mean SDOS in the $`a`$\- and $`b`$-band, respectively. (a) Effect of impurity concentration ($`u=u_a=u_b`$) at $`\alpha _a=\alpha _b=0`$ and $`|\overline{v}_{ab}|=1`$. Since SDOS’s are almost the same in the two bands at $`u=5`$, we show one of them only in the figure. (b) Effect of the magnitude of the interband scattering $`|\overline{v}_{ab}|`$ at $`\alpha _a=\alpha _b=0`$ and $`u_a=u_b=10`$. (b) Effect of the intraband scattering $`\alpha (=\alpha _a=\alpha _b)`$ at $`|\overline{v}_{ab}|=1`$ and $`u_a=u_b=10`$.
10. Superconducting density of states at $`T=0.1T_\mathrm{c}`$ in the case of $`\theta _{ab}=\pi /2`$. We take $`\overline{g}_{aa}=0.3`$, $`\overline{g}_{bb}=0.25`$ and $`\overline{g}_{ab}=0.1`$. (a) Effect of $`u(=u_a=u_b)`$ at $`\alpha _a=\alpha _b=0`$ and $`|\overline{v}_{ab}|=1`$. (b) Effect of $`|\overline{v}_{ab}|`$ at $`\alpha _a=\alpha _b=0`$ and $`u_a=u_b=5`$. (c) Effect of the intraband scattering $`\alpha `$ at $`|\overline{v}_{ab}|=1`$ and $`u_a=u_b=5`$. Since $`N_a(\omega )`$ and $`N_b(\omega )`$ are almost the same at $`u=5`$ in panel (a), $`|\overline{v}_{ab}|=100`$ in panel (b) and $`\alpha =0`$ in panel (c), only $`N_a(\omega )`$ is shown in these cases.