electron scattering off magnetic impurities , unlike the potential scattering , substantially modifies the properties of _ s_-wave superconductors .
it was realized in the seminal work by abrikosov and gorkov@xcite ( ag ) that the presence of magnetic impurities reduces the superconducting gap below its bcs value and eventually may lead to the phenomenon of gapless superconductivity .
the gap suppression was investigated by means of tunneling between normal and superconducting electrodes.@xcite beyond the ag theory , it was realized that a single magnetic impurity creates a localized state within the bcs gap.@xcite a small concentration of magnetic impurities results in the formation of an impurity band@xcite which merges with the continuum as the concentration exceeds a certain value ( see ref .   for a review ) .
( the original ag theory is valid above that value . )
the impurity band was investigated by tunneling experiments with alloys ( such as pbmn)@xcite and normal metal
superconductor bilayers,@xcite and also by means of thermal transport in superconducting films.@xcite the observation of discrete levels associated with a separate magnetic impurity came later in the scanning tunneling microscopy experiments.@xcite    the effect of magnetic impurities came under scrutiny recently in connection with limits on performance of superconducting elements in qubits@xcite and photon detectors.@xcite in these applications , superconductors should perform in the ac regime .
there is also some experimental effort directly aimed at observation of the effect of magnetic impurities on the surface impedance of superconducting multilayers@xcite and unconventional superconductors.@xcite theoretically , the results on impedance@xcite were analyzed in the framework of a simple two - fluid model .
the two - fluid model assumes a simplified frequency - independent drude - type dissipative conductivity @xmath0 , while the influence of magnetic impurities on the electron spectrum should lead to a nontrivial @xmath1 dependence on the frequency @xmath2 .
a theoretical approach to studying the surface impedance in superconductors was developed in the celebrated work of mattis and bardeen,@xcite which treated the anomalous skin effect ( the skin depth is smaller than the mean free path ) taking nonmagnetic disorder into account .
independently , the clean case was considered by abrikosov , gorkov , and khalatnikov@xcite in the green - function language .
later , nam@xcite employed the same technique and demonstrated that the extreme anomalous limit is equivalent ( up to the expressions for the normal - metallic conductivities ) to the local limit , realized in the dirty regime .
the complex conductivity of a superconductor containing both potential and magnetic impurities was addressed by skalski _
_ et al.__@xcite their paper does contain general results applicable in the case of weak magnetic scattering ; however , similarly to the ag theory , scattering off magnetic impurities was considered in the born approximation , which did not allow to treat the effect of localized states on the impedance . in ref .
, larkin and ovchinnikov considered superconductors with inhomogeneous pairing , which under certain assumptions is formally equivalent to the presence of magnetic impurities in the ag regime . some of their results can therefore be applied to this problem .
however , a convenient formulation of the theory , suitable for the experimentally relevant dirty limit and making it possible to treat both the ag regime and the impurity band , as well as detailed quantitative investigation of the effect of magnetic impurities on the surface impedance in different regimes over temperature , frequency , and electron scattering rate off magnetic impurities is still lacking , to the best of our knowledge .    in this paper , we study the influence of magnetic impurities on the complex conductivity and surface impedance of a superconductor .
we are especially interested in the limit of low frequency of the impinging radiation and low temperature , and in the effect of the impurity band on the dissipation .
quantitative results are presented for the experimentally relevant dirty limit ( with respect to the nonmagnetic disorder ) .
the paper is organized as follows .
first , in sec .  [ sec : clean ] , we consider the limit of weak scattering off magnetic impurities and disregard nonmagnetic disorder ( `` clean '' limit ) . in that limit
, we find the localized states bound to magnetic impurities and discuss the formation of the impurity band at finite but small concentration of magnetic impurities and the formation of the ag density of states at higher concentrations . we also present there a qualitative consideration of the dissipative conductivity .
the main features of the developed simple picture carry over to the realistic case of a `` dirty '' superconductor doped with magnetic impurities . in sec .
[ sec : dirty ] , we consider the dirty limit for the nonmagnetic disorder in the framework of the quasiclassical formalism of the usadel equation .
we first obtain the spectrum of quasiparticle excitations in the presence of magnetic impurities ( treated beyond the born approximation ) .
then we address the effect of magnetic impurities on the low - temperature dissipative conductivity .
in addition to the low - frequency limit , we also consider in detail the frequency dependence of the dissipation near the thresholds determined by the gaps in the excitation spectrum .
our results are summarized in sec .
[ sec : concl ] .
some technical details are presented in the appendixes .
to highlight the main features of the problem related to the impurity states , we consider here a simplified case of a superconductor with magnetic impurities only .
in addition , considering the hybridization of single - impurity states into an impurity band we first assume magnetic impurities to be polarized ( magnetic moments pointing up the @xmath3 axis ) , while the impurity positions are random .
furthermore , we assume the exchange constant @xmath4 of interaction of localized classical spins @xmath5 with the spins of itinerant electrons to be small . introducing the dimensionless parameter @xmath6 describing the `` strength '' of a single magnetic impurity ,
we can write the latter requirement as @xmath7 here @xmath8 is the electron density of states ( dos ) per spin projection at the fermi level in the normal state , with @xmath9 and @xmath10 being the fermi wave vector and the fermi velocity , respectively .
weak exchange interaction will modify substantially the excitation spectrum only close to the gap . concentrating on that energy region ,
we simplify the hamiltonian of the system in the following way .
first , we express the operator of exchange interaction of electrons with an impurity ( placed at origin , @xmath11 ) in terms of the bogoliubov quasiparticles ( @xmath12 and @xmath13 are the operators of electrons and bogoliubov quasiparticles , respectively ) : @xmath14,\end{gathered}\ ] ] where @xmath15 are momenta and @xmath0 is the spin index ( @xmath16 for @xmath17 , @xmath18 for @xmath19 ) . at the same time , the standard bcs part of the hamiltonian becomes diagonal after the bogoliubov transformation , with the excitation energies @xmath20 in the superconducting state . here
@xmath21 is the bcs gap and @xmath22 are the excitation energies in the normal state .
close to the gap ( @xmath23 ) , we can approximate @xmath24 , while the bogoliubov coherence factors are simply @xmath25 , so finally the total hamiltonian takes the form @xmath26 where @xmath27 and we measure the quasiparticle energies from the edge of the continuous spectrum @xmath21 : @xmath28      it is easy to see that hamiltonian ( [ ham ] ) has a bound state for quasiparticles with spin @xmath0 satisfying the condition @xmath29 . for definiteness
, hereafter we set @xmath30 ; then the bound state exists for spin down .
we denote the wave function of the bound state in the momentum representation as @xmath31 ; the schrdinger equation for it , which follows from the form of hamiltonian ( [ ham ] ) , reads @xmath32 clearly , the wave function must have the form @xmath33 . substituting it into eq .
( [ eqphi ] ) , we find the equation for the energy of the bound state ( @xmath34 ) , @xmath35 yielding @xmath36 at @xmath37 . the power - law dependence , @xmath38 , is a direct consequence of the excitations spectrum @xmath39 . the convergence of the integral in eq .
( [ sch_e ] ) at @xmath40 validates the use of the simplified form of @xmath41 .
we emphasize that in order to obtain the bound state , it was essential to go beyond the born approximation.@xcite the wave function of the bound state in the coordinate representation has the form @xmath42 where @xmath43 is the coherence length of a clean superconductor .
the wave function of the localized state decays exponentially on the spatial scale which exceeds the coherence length by the factor @xmath44 .
clearly , the effect of multiple magnetic impurities on the quasiparticle spectrum depends crucially on the characteristic distance between them compared to @xmath45 . at lowest concentrations @xmath46 ,
separate magnetic impurities create a @xmath47-peak in the quasiparticle dos . with the increase of the concentration the peak broadens .
as long as the average inter - impurity distance is large , the broadening is due to rare occurrences of closely located impurity pairs producing split levels due to their hybridization.@xcite however , at higher concentrations , @xmath48 , the bound state on a given impurity overlaps with a large number of other localized states and forms a well - defined impurity band .
this large parameter allows one to treat the dos in that band self - consistently .
we start with an estimate of the impurity band width @xmath49 .
a quasiparticle initially localized on a given impurity may hop on @xmath50 `` nearest neighbors '' . because of their large number , the return of the quasiparticle is improbable .
therefore , we may introduce the escape rate @xmath51 for a given impurity .
it will determine the level width , and we will associate it with the width of the band , @xmath52 . on the other hand ,
the escape rate from a given site is proportional to the dos on the `` receiving '' ones .
the latter is inversely proportional to @xmath49 .
this way , the formula for the escape rate turns into a self - consistent equation for @xmath49 ( this is the essence of the self - consistent born approximation for degenerate states , analyzed in detail in the context of the broadening of high landau levels@xcite ) . to implement the scheme , we denote the hopping matrix element between two impurity states @xmath53 and @xmath54 , at positions @xmath55 and @xmath56 , as @xmath57
. the said self - consistency equation reads @xmath58 ( upon summing over positions @xmath54 , the dependence on @xmath53 disappears ) .
the estimate of @xmath57 should be obtained from a consideration of the two - impurity problem .
tunneling between two impurity sites splits the degenerate level in two .
the resulting energy splitting depends on the inter - impurity distance @xmath59 . similarly to the conventional tight - binding problem
, it is this energy splitting that should be identified with @xmath60 . this way
, we find @xmath61 .
to estimate the proportionality coefficient here , we note that the energy splitting reaches a value of the order of @xmath62 if @xmath63 .
therefore , @xmath64 .
this estimate is confirmed by a solution of the energy splitting problem , which is easily obtained by the generalization of eqs .
( [ ham])([sch_e ] ) to the case of two impurities.@xcite    summation over the random positions of the impurities @xmath54 results in @xmath65 . remarkably , we used here only the value of the wave function @xmath66 close to the impurity and the fact that @xmath66 is normalized .
this is why one may expect that the obtained estimate of @xmath67 remains valid in the presence of nonmagnetic disorder , as long as the corresponding mean free path is large compared to the electron fermi wave length .
( this expectation is indeed confirmed by the rigorous consideration of sec .
[ sec : dirty ] . ) substituting the estimate of @xmath67 into eq .
( [ w ] ) and replacing @xmath68 , we find expression for @xmath49 up to a numerical factor . together with this factor ( which will be found later ) ,
the bandwidth in terms of measurable quantities @xmath46 and @xmath69 , reads @xmath70 this expression is valid when the band is narrow , @xmath71 .
the next question is the dos @xmath72 for the bogoliubov quasiparticles ( hence the @xmath73 subscript ) inside the impurity band which can be found from the momentum - integrated green function @xmath74 as @xmath75 the green function @xmath76 of the bogoliubov quasiparticles is determined by the hamiltonian , which now , in contrast to hamiltonian ( [ ham ] ) , contains magnetic impurities of finite concentration : @xmath77 the green function can be written in terms of the self - energy @xmath78 and we treat the self - energy within the self - consistent @xmath79-matrix approximation,@xcite which diagrammatically amounts to the geometric series of fig .
[ fig : scba ] : @xmath80 ( the sign in front of @xmath4 in the self - energy corresponds to quasiparticle spin down , which is of interest for us ) . integrating eq .
( [ g ] ) over momentum , we obtain the self - consistency condition for @xmath81 : @xmath82 the impurity band corresponds to nonzero dos , that is , to the nonzero imaginary part of @xmath83 , in some domain of negative @xmath84 . analyzing eq .  ( [ g ] ) in the limit @xmath85 ( the limit of narrow band ) , we find a semicircle dos inside the impurity band : @xmath86 which is centered at @xmath69 , found earlier in eq .
( [ bs ] ) , and has the width given by eq .
( [ wfin ] ) [ note that it is the calculation leading to eq .
( [ nuom ] ) that establishes the numerical factor in eq .
( [ wfin ] ) ] .
the total dos inside the impurity band is equal to @xmath46 , demonstrating that each impurity brings one state into the system.@xcite the same eq .
( [ g ] ) can be employed to demonstrate that the lower edge of the continuous spectrum shifts from zero to @xmath87 at finite impurity concentrations , and the behavior of the dos near the edge is @xmath88 ( the positive value of @xmath89 is a manifestation of the level repulsion between the states of the continuum and the impurity states . )     corresponding to the self - consistent @xmath79-matrix approximation .
the scattering on a single impurity ( dashed lines ) is treated in all orders of the perturbation series .
the self - consistency is achieved by inserting the bold lines for the green functions , which means that they are calculated with the help of the presented self - energy . ]    as we mentioned before , a well - defined impurity band exists [ and can be described by the self - consistent @xmath79-matrix method leading to dos ( [ nuom ] ) ] if @xmath48 , which translates into @xmath90 , where @xmath91 is the fermi energy . with the increase of the impurity concentration , the impurity band broadens , @xmath92 . at @xmath93 ,
the width of the band becomes of the same order as the distance from the center of the band to the continuum edge , resulting in the merger of the impurity band with the continuum of quasiparticle states .
[ quantitatively , the analysis of eq .  ( [ g ] ) shows that the energy gap between the impurity band and continuum exists if @xmath94 . ]
the two conditions for the existence of the impurity band , @xmath48 and @xmath95 , are consistent since @xmath96 .    .
]    in the above discussion , we did not treat the self - consistency equation explicitly , assuming @xmath21 that we used to be the self - consistent order parameter corresponding to the given impurity parameters . actually , @xmath21 is suppressed by magnetic impurities
. a gapless state is realized at a critical concentration @xmath97 . at that concentration the average zeeman energy associated with the magnetic impurities ,
@xmath98 is comparable to @xmath21.@xcite figure  [ fig : axis_ns ] illustrates the results of the above discussion .
the above consideration of spin - polarized magnetic impurities can be easily generalized to the case of randomly oriented spins ( see appendix  [ app : relation ] ) .
it yields eq .
( [ app : nuom3 ] ) for the dos  now summed over the quasiparticle spin directions  inside an impurity band of width @xmath99 [ see eq .
( [ app : wfin3 ] ) ] . comparing eqs .
( [ app : nuom3 ] ) and ( [ app : wfin3 ] ) with eqs .
( [ nuom ] ) and ( [ wfin ] ) , we notice the relation @xmath100 between the widths of the band . indeed ,
randomness of impurity spin orientations leads to an additional @xmath101 factor in the hopping matrix element @xmath57 between two impurity states [ see eq .
( [ w ] ) ] , and hence to the @xmath102 factor in the average @xmath103 , finally changing the width of the band by the @xmath104 factor .
( here @xmath105 is the angle between the spins of the two impurities . )
another difference of the random - spin case from the spin - polarized case is the absence of the average zeeman energy . upon increasing @xmath46 , the gapless state
is then reached due to the electron spin - flip scattering off magnetic impurities.@xcite the latter is characterized by a scattering rate @xmath106 which , in the born approximation , is @xmath107 the transition occurs when @xmath108 , that is , at a concentration of magnetic impurities @xmath109 .      at finite but low temperatures @xmath79 ,
quasiparticles first populate the impurity band .
this results in a finite dissipative part @xmath110 of conductivity at arbitrary low frequencies .
we estimate @xmath110 at @xmath111 using a suitable version@xcite of the mattis - bardeen theory.@xcite we disregard possible spin selection rules@xcite by assuming that the spins of magnetic impurities are randomly oriented
.    we may relate @xmath112 to the absorption power per unit volume , @xmath113 , where @xmath114 is the amplitude of the applied ac electric field .
the absorption power can be evaluated as the product of the photon energy quantum @xmath2 and the rate of creation of quasiparticles .
the latter is @xmath115 and is evaluated with the help of the fermi golden rule .
the perturbation leading to the excitation of the system is @xmath116 , where @xmath117 is the operator of velocity of an electron and @xmath118 is the vector potential of the ac field ; its amplitude is @xmath119 .
the absorption power associated with the photon - induced transitions of quasiparticles within the impurity band equals@xcite @xmath120 \\
\times \nu_b ( e_1 ) \nu_b ( e_2 ) ( u_1 u_2+v_1 v_2)^2 \delta ( e_1-e_2-\omega).\end{gathered}\ ] ] writing eq .
( [ abspow2 ] ) , we accounted for the fact that the states in the impurity band do not correspond to a defined momentum , and consequently the latter is not conserved in the transitions . to estimate the average over the states @xmath121 , we employ the arguments of the work of mattis and
bardeen,@xcite which allow one to express it in terms of an analog of the normal conductivity for the impurity band .
that leads us to the estimate @xmath122 , where @xmath123 is the relaxation time and @xmath124 is a characteristic value of the dos in the impurity band . substituting @xmath125 and @xmath126
[ see eq .
( [ app : nuom3 ] ) ] , we find @xmath127 in the case of `` shallow '' levels ( @xmath37 ) we may also replace the combination of the coherence factors , @xmath128 . at lowest frequencies , @xmath129 , we may expand the difference of the quasiparticle distribution functions in eq .
( [ abspow2 ] ) as @xmath130 . as a result
, we find frequency - independent asymptotic forms for the dissipative conductivity at @xmath131 and @xmath132 , which we express in terms of the density @xmath133 of quasiparticles in the impurity band [ hence the @xmath134 superscript ] : @xmath135 in equilibrium , the density of quasiparticles is @xmath136 in the opposite limit of photon energies high compared to temperature , we may neglect @xmath137 in eq .
( [ abspow2 ] ) , finally finding @xmath138    in an ideal superconductor at zero temperature , photon absorption is possible only for photons with energy @xmath2 exceeding the threshold @xmath139 ( each photon creates a pair of quasiparticles ) .
the presence of states within the gap @xmath21 reduces the threshold frequency and modifies the dependence of the dissipative part of conductivity @xmath112 on frequency for @xmath2 exceeding the threshold . similarly to the above consideration , one may relate @xmath110 to the absorption power ( cf .  ref .  ) , @xmath140 as follows from eq .
( [ app : nuom3 ] ) , absorption is possible only if the frequency exceeds the threshold , @xmath141 to evaluate the absorption , however , the accepted approximation of shallow levels is insufficient , because it yields @xmath142 . as follows from the detailed calculation presented in sec .
[ sec : dissip_zerot ] this factor @xmath143 ; the dissipative conductivity scales as @xmath144 near the threshold .
we turn now to the practically important case of short elastic mean free time @xmath145 caused by potential scattering of electrons and consider the dirty limit , @xmath146 .
we will see that potential disorder does not affect the electron spectrum.@xcite however , its role is essential for transport properties ( complex conductivity and impedance ) .    in the dirty limit , magnetic impurities of small concentration can be described in the framework of the usadel equation.@xcite the form of the usadel equation , widely used in literature,@xcite corresponds to the born limit for magnetic scattering ( with random positions and orientations of magnetic impurities ) . at the same time , the generalization to arbitrary strength of magnetic scattering can also be obtained , as shown by marchetti and simons.@xcite employing the standard @xmath147-parametrization of the normal and anomalous quasiclassical green functions , @xmath148 we can write the generalized form as@xcite @xmath149 where @xmath150 is the diffusion constant and @xmath151 is the electron energy measured from the fermi level .
we are interested in the homogeneous solution .
following ref .  , we define a `` renormalized '' energy and order parameter , @xmath152 with these quantities , the homogeneous usadel equation becomes simple , @xmath153 and can be formally solved as @xmath154 where @xmath155 . at the same time , this is not a solution yet , but an alternative parametrization of the green functions in terms of @xmath156 .
the usadel equation actually determines the green functions dependence on energy @xmath151 ; hence , in eq .
( [ tetd ] ) @xmath156 itself depends on @xmath151 .
the latter dependence is determined by the equation , directly following from eqs .
( [ tildee ] ) and ( [ tilded ] ) : @xmath157 where @xmath158 here @xmath159 is the exact scattering rate due to magnetic impurities , normalized by @xmath21 , and @xmath160 can be shown to determine the center of the impurity band , also in units of @xmath21 . since we focus on the case when each impurity is weak , @xmath37 , we simplify the above expressions as @xmath161    once found , @xmath156 can be immediately applied for analyzing the spectrum .
the dos per spin is given by @xmath162 it turns out that the potential scattering does not influence the energy spectrum , since eqs .
( [ u_eq ] ) , ( [ gammaepsilon ] ) , and ( [ dos_u ] ) coincide@xcite with those obtained by shiba@xcite and rusinov@xcite in the clean limit ( see also appendix  [ app : relation ] ) .
finite dos corresponds to complex solutions of eq .
( [ u_eq ] ) for @xmath156 . analyzing the behavior of the function in the left - hand side ( l.h.s . ) of eq .
( [ u_eq ] ) at @xmath163 , shiba@xcite found the domain of energies [ the right - hand side ( r.h.s . ) of eq .
( [ u_eq ] ) ] that can not be matched by _
real _ values of @xmath156 . at positive energy , this domain is the impurity band centered around @xmath164 .
calculating @xmath156 inside the narrow band , we find @xmath165 for the dos and @xmath166 for the width . the dos at @xmath167 is the mirror image of eq .
( [ nuom3 ] ) .
the total number of states in the two impurity bands equals @xmath46 . the dos ( [ nuom3 ] ) coincides with the result derived in the clean limit
[ see eq .
( [ app : nuom3 ] ) ] up to a factor of 4 .
this factor accounts for the difference in the definitions of dos .
[ the dos ( [ app : nuom3 ] ) , defined for bogoliubov quasiparticles existing only at @xmath168 , corresponds to eq .
( [ nuom3 ] ) ( electron component ) folded together with its mirror image ( hole component ) and , additionally , summed over the spin projections . ]    as the concentration of magnetic impurities ( and hence @xmath159 ) grows , the upper edge of the band merges with the continuum at @xmath169 [ this result can be obtained from eq .
( [ u_eq ] ) and differs only numerically from the condition obtained for the case of polarized magnetic impurities in sec .
[ subsec : imp - band ] ] . at higher concentrations the ag regime
is realized , with the gap width reduced due to the scattering off magnetic impurities . as follows from eq .
( [ u_eq ] ) , the quasiparticle continuum starts at energy lower than the single - impurity bound - state energy @xmath170 , and the gap @xmath171 is reduced progressively with the increase of @xmath159 . at @xmath172 , one may replace @xmath173 with @xmath174 in eq .
( [ u_eq ] ) .
then it simplifies to the form @xmath175 exactly as considered by abrikosov and gorkov.@xcite note that @xmath159 itself can still be small , and we mainly focus on that case : @xmath176    the ag regime , corresponding to a single gap in the spectrum , has two characteristic features , essential for the low - frequency dissipation : first , the gap is suppressed ( the relative scale is determined by @xmath177 ) , and second , the square - root bcs singularity of the dos at @xmath178 is smeared .
the gap @xmath171 in units of @xmath21 is@xcite @xmath179 solving the usadel equation at energies close to the actual gap , @xmath180 , we obtain@xcite @xmath181 the square - root behavior of the dos [ see eq .
( [ dos_u ] ) ] in the vicinity of the gap was found in refs .  .
at low concentration , magnetic impurities substantially modify the spectrum in a relatively narrow domain of energies , @xmath182 .
we take @xmath178 as a reference point within the domain , and find @xmath183 and its derivative with respect to energy : @xmath184 the maximum of the dos is achieved slightly above @xmath178 . at the same time , considering the derivative @xmath185 and the bcs solution that makes it possible to approach the maximum from the side of higher energies , we conclude that the maximal dos is of the same order as the value determined by eq .
( [ at1a ] ) .
this result was found in ref .  .
the results for the dos are summarized in fig .
[ fig : dos ] .
note that the square - root behavior of the dos near the gap edge [ following from eq .
( [ neargap ] ) ] , being applied at @xmath178 ( which is already beyond its applicability range ) , yields @xmath186 , which only slightly differs from the accurate result @xmath187 following from eq .
( [ at1a ] ) .    in the evaluation of dissipation with the help of eq .
( [ dissip ] ) , a good approximation for the dos is given by @xmath188 at the matching points , the above expressions match by the order of magnitude , being different only by numerical factors close to 1 .
therefore , this is a rather accurate approximation for calculating the dissipation up to order - of - one numerical factors ( at the same time , in the limiting cases that we consider below , the numerical factors will be asymptotically exact ) .    at small @xmath189 and @xmath159 , magnetic impurities significantly affect the electron spectrum only at energies close to @xmath21 . in that energy range [ corresponding to eq .
( [ nuom3 ] ) and to the first two lines in eq .
( [ dos ] ) ] , we have @xmath190 and hence @xmath191 [ as follows from eq .
( [ tetd ] ) ] , so that finally @xmath192 . outside that domain , the functions @xmath193 and @xmath183 in the leading order have the conventional bcs form [ with @xmath194 instead of the last line in eq .
( [ dos ] ) ] .
note that we have discussed @xmath195 and @xmath196 at @xmath168 . at the same time , @xmath197 is even while @xmath198 is odd with respect to @xmath151 .
having found the green functions , we are ready to calculate the dissipation .
the expression for the response kernel of the current to the external electromagnetic field is derived in appendix  [ sec : app : kernel ] . from the result for the kernel
, we extract the expression for the dissipative part of conductivity : @xmath199,\end{gathered}\ ] ] where @xmath200 and @xmath201 is the drude conductivity .
this formula is the dirty - limit counterpart of eqs .
( [ abspow2 ] ) and ( [ abspow3 ] ) for the absorption power in the clean limit .
the combination of the green functions @xmath183 and @xmath193 here replaces the product of the quasiparticle densities of states and the combinations of the coherence factors in eqs .
( [ abspow2 ] ) and ( [ abspow3 ] ) .
the analysis of the present section is be based on eq .
( [ dissip ] ) . at the same time
, the results can be easily converted into the language of the surface impedance , as discussed in sec .
[ sec : app : impedance ] .
small but finite dissipative conductivity implies a small but finite real part of the surface impedance .      at zero temperature
the dissipation is absent below a certain threshold frequency . in the absence of magnetic impurities
, @xmath2 must exceed @xmath139 so that a cooper pair can be broken and two quasiparticles are created above the gap . at frequencies slightly above the threshold , @xmath203 , the dissipative part of conductivity is@xcite
@xmath204    in the presence of magnetic impurities , the bcs singularity of the green functions is smeared out and hence the threshold behavior of the dissipation changes .
the threshold frequency in the ag regime is @xmath205 , and the nonzero contribution to the integral in eq .
( [ dissip ] ) comes only from the region of energies @xmath206 , where the arguments of the green functions are close to @xmath207 so that , at first sight , we can employ the square - root form for @xmath197 and @xmath198 [ the first line in eq .
( [ dos ] ) ] and take into account that @xmath195 is even while @xmath196 is odd .
however , approximating @xmath208 at @xmath168 we find that the integrand in eq .
( [ dissip ] ) vanishes in the main order , so we should take into account a small difference between @xmath197 and @xmath198 .
this is easily done on the basis of the normalization condition @xmath209 .
note that it is the large absolute value of @xmath183 near @xmath171 [ see eq .
( [ neargap ] ) ] that leads to @xmath191 and hence to @xmath192 . at the same time
, the normalization condition immediately gives also a correction to this result : @xmath210 then we find nonzero dissipation which is quadratic in deviation from the threshold at @xmath211 : @xmath212 this result was obtain in ref .  .
the quadratic threshold behavior of @xmath112 is a consequence of the square - root threshold behavior of the dos , and experimental results of ref .
seem to confirm this @xmath110 dependence .    at lowest concentrations of magnetic impurities
, the impurity band exists inside the gap , so that the threshold frequency for absorption is determined by transitions in which a cooper pair is broken and two quasiparticles are created inside the impurity band . the threshold frequency for this process is @xmath213 , which is eq .
( [ omth ] ) written in different notation . just above the threshold , at @xmath214 , only the edge behavior ( near an edge of the impurity band ) of @xmath197 and
@xmath198 enters eq .
( [ dissip ] ) . according to eq .
( [ nuom3 ] ) , both functions near an edge grow as a square root .
we again find a cancellation in the main order and have to employ eq .
( [ imfreg ] ) .
the functional form is thus the same as in the ag regime , and the result differs from eq .
( [ sgma ] ) only by coefficients : @xmath215    next we consider the case of finite ( but low ) temperatures , when the threshold for absorption is absent since equilibrium thermally excited quasiparticles can absorb photons of arbitrarily small energy .
we start with the consideration of the ag regime in the limit of small @xmath159 , so that quasiparticle states exist only above the gap @xmath171 , which is weakly suppressed compared to @xmath21 , as described by eqs .
( [ gammas ] ) and ( [ epsilong ] ) .
most of the results of this subsection were obtained in previous publications , although sometimes in the framework of physically different problems and with the help of different techniques .
our consideration makes it possible to incorporate them in a unified manner .    in the limit that we consider , @xmath195 and @xmath196 entering eq .
( [ dissip ] ) are nonzero only above @xmath171 which is close to @xmath21 , so that at low frequencies and temperatures , @xmath216 two simplifications are possible : ( a )  we can substitute the @xmath217 s with their asymptotic exponential forms , and ( b )  since the exponentials limit the integration to the region of @xmath218 , where @xmath196 is approximately equal to @xmath195 , the contribution of the @xmath219 term is the same as the contribution from the @xmath220 term .
therefore , we can simplify eq .
( [ dissip ] ) as @xmath221 at the same time , at @xmath222 the external electromagnetic field can not excite quasiparticles across the gap ; therefore , the dissipation will be due to quasiparticles that are already excited due to low but finite temperature .
it is thus natural to express the results for the dissipation in terms of the equilibrium quasiparticle density , @xmath223 for which under the same assumptions ( @xmath224 , @xmath225 ) we obtain two limiting cases : @xmath226    below we consider various relations between the three energy scales @xmath2 , @xmath79 , and @xmath227 ( while all of them are much smaller than @xmath21 ) . starting with the case of relatively strong magnetic scattering , @xmath228 , we note that , due to the exponentials , the main contribution to the integral in eq .
( [ sigmagg ] ) comes from the narrow region above @xmath171 , where the square - root asymptote for the dos is valid [ see the first line in eq .
( [ dos ] ) ] . as a result ,
@xmath229 where @xmath230 is the macdonald function ( the modified bessel function of the second kind ) .
asymptotic forms of the macdonald function at large and small values of the argument describe the two limiting cases for the @xmath231 ratio , so that finally at @xmath232 we reproduce the result of maki,@xcite @xmath233 while at @xmath234 we find @xmath235    the result ( [ k1 ] ) was actually obtained by larkin and ovchinnikov in a physically different but formally equivalent problem.@xcite they considered a superconductor with inhomogeneous effective interaction between electrons . for small - scale inhomogeneities the problem turns out to be mathematically equivalent to the case of homogeneous effective interaction in the presence of magnetic impurities .
the results of ref .   thus apply to the problem considered in the present paper .    in the limit of weak magnetic scattering , @xmath236 , the main contribution to the integral in eq .
( [ sigmagg ] ) comes from the region of energies above the smeared bcs singularity ; in that region we can apply the unperturbed result for the dos [ see the last line in eq .
( [ dos ] ) ] .
we thus reproduce the results obtained in the absence of magnetic impurities@xcite ( we have to take into account the correspondence noted by nam@xcite between the extreme anomalous limit in the clean case@xcite and the dirty case to which our calculations refer ) : @xmath237 asymptotic forms of the macdonald function @xmath238 at large and small values of the argument describe the two limiting cases for the @xmath231 ratio , so that finally at @xmath239 we find @xmath240 ( where @xmath241 is euler s constant ) , while at @xmath242 the result is @xmath243 the logarithmic divergence in eq .
( [ dis4 ] ) at @xmath244 is due to the overlap of two square - root bcs singularities of the dos .
similarly to the case of eq .
( [ dis1 ] ) , at @xmath245 the second exponential in eq .
( [ sigmagg ] ) can be neglected while the @xmath2-shifted dos is nearly constant , @xmath246 , in the essential region of integration . as a result ,
the remaining integral in both cases is exactly the same as the one determining the quasiparticle density ( [ nqp ] ) at @xmath247 .
therefore , we obtain exactly the same result ( [ dis1 ] ) in terms of @xmath248 , although the expressions for @xmath248 itself are different at @xmath242 and @xmath245 [ eqs .
( [ nqp1 ] ) and ( [ nqp2 ] ) , respectively ] .
this result actually follows from ref .  .
finally , at @xmath249 , relatively large temperature allows rather wide region of energy integration in eq .
( [ sigmagg ] ) , while the shift of the two dos can be neglected .
similarly to the case of eq .
( [ dis4 ] ) , a logarithmic singularity appears in the integral ; however , now it is cut off not due to the finite shift @xmath2 [ as it was the case in eq .
( [ dis4 ] ) ] but due to smearing of the square - root bcs singularity of the dos by magnetic impurities . as a result ,
with the logarithmic accuracy we find @xmath250    the results for @xmath110 at @xmath222 are valid as long as the quasiparticle population may be described by the boltzmann distribution with an effective temperature .
the quasiparticle chemical potential should not necessarily be @xmath251 and the temperature @xmath79 may deviate from equilibrium as long as @xmath252 .      at equilibrium , the results for the dissipation @xmath112 calculated in sec .
[ sec : ag ] for the ag regime in the limit of low temperatures always contain an exponentially small factor @xmath253 , since the absorption of the electromagnetic field is due to the quasiparticles thermally excited above the gap . at the same time
, as we discussed above , at low concentrations of magnetic impurities a band of quasiparticles states appears below the continuum .
this means that at temperatures much smaller than the distance between the band and the continuum , the dominant contribution to the dissipation will be due to quasiparticles residing inside the impurity band .
below we consider this case , which is the main focus of our interest .
our aim now is to calculate the dissipation ( [ dissip ] ) for the case when it is due to transitions inside the impurity band ; hence , @xmath254 . at the same time
, we assume the band to be narrow , so that its width , which can be written as @xmath255 , is much smaller than the distance from the center of the band @xmath170 to the bottom of the continuum @xmath21 : since @xmath256 , the condition can be written as @xmath257 the contribution of the impurity band is dominant at @xmath258 , when the continuum contribution is exponentially smaller due to the @xmath259 factor .    since we consider energies of the order of @xmath170 , that is , of the order of @xmath21 , and @xmath260 inside the narrow band ( as discussed in sec .
[ sec : dos ] ) , we can make the same simplifications that lead us to eq .
( [ sigmagg ] ) .
the only difference is that now the region of the energy integration is confined to within the impurity band . taking into account the explicit form of the dos inside the band [ eq .
( [ nuom3 ] ) ] and making the variable of integration dimensionless , we obtain @xmath261 the integral contains the product of the two semicircle dos ( [ nuom3 ] ) shifted by the frequency @xmath2 .
the overlap exists and produces nonzero dissipation only at @xmath262 .    similarly to sec .
[ sec : ag ] , we express the final results in terms of the equilibrium quasiparticle density ( within the band ) . calculating eq .
( [ nqp ] ) for the impurity band , we find @xmath263 where @xmath264 is the modified bessel function of the first kind . in the two limiting cases for the @xmath265 ratio ,
we obtain @xmath266    below we consider various relations between the three energy scales , @xmath2 , @xmath79 , and @xmath99 ( while all of them are much smaller than @xmath21 ) . starting with the case of relatively large band width , @xmath267 , we note that , due to the exponentials , the main contribution to the integral in eq .
( [ sigmaggb ] ) comes from the narrow region above the lower limit of integration where the simplified square - root result for the dos is valid : @xmath268 .
this situation is very similar to the case of eq .
( [ k1 ] ) ( cf .  the discussion above that formula ) .
as a result , we obtain the same functional form [ with a different prefactor and with @xmath269 instead of @xmath171 in the exponent ] : @xmath270 asymptotic forms of the macdonald function at large and small values of the argument describe the two limiting cases for the @xmath231 ratio , so that finally at @xmath271 we find @xmath272 while at @xmath273 the result is @xmath274    at @xmath275 ( while @xmath99 and @xmath2 are not necessarily close to each other ) , the exponential in eq .
( [ sigmaggb ] ) can be replaced by unity in the region of integration , and the remaining integral can be written as @xmath276.\end{gathered}\ ] ] here @xmath277 and @xmath151 are the complete elliptic integrals of the first and second kind , respectively , and their argument is @xmath278 in the limit of small frequency and high temperature , @xmath279 , with the help of eq .
( [ nqp1b ] ) , expression ( [ ek ] ) can be simplified as @xmath280    there is an upper threshold for impurity band absorption at @xmath281 .
close to that threshold , eq
.  ( [ ek ] ) gives @xmath282 , where @xmath147 is the heaviside step function .
comparing eqs .
( [ eq : sigmab - dirty1 ] ) , ( [ eq : sigmab - dirty3 ] ) , and ( [ eq : sigmab - dirty2 ] ) in the dirty case with eqs .
( [ eq : sigmab - clean1 ] ) and ( [ eq : sigmab - clean2 ] ) in the clean case , we notice that the dependences of the dissipative part of conductivity on @xmath2 , @xmath79 , and @xmath99 are the same . at the same time , the overall coefficients in the two cases are different due to different limits with respect to potential scattering .      at low temperatures ,
the density of equilibrium quasiparticles becomes negligible , and in reality @xmath248 is dominated by extraneously produced quasiparticles.@xcite in quantum devices based on josephson junctions , the adverse effect of quasiparticles on the coherence may be mitigated by the inclusion of traps  sections of the superconducting leads with reduced value of the gap .
we may pose a similar question with respect to the impedance : can magnetic impurities reduce the dissipation caused by nonequilibrium quasiparticles ? such a possibility is most intriguing in the range of frequencies @xmath283 [ see eq .
( [ bs ] ) ] and quasiparticle densities @xmath284 ( see fig .  [
fig : axis_ns ] ) . in the following ,
we use the results of the two previous subsections to elucidate the dependence of the normalized dissipative conductivity @xmath285 on @xmath46 .    for estimates ,
we assume the density of quasiparticles @xmath248 fixed , and their effective temperature low , @xmath286 , so that thermal ionization to the continuum does not prevent magnetic impurities from trapping the quasiparticles .
in the absence of magnetic impurities , the dissipative response is defined by eq .
( [ dis1 ] ) . at small density , impurities reduce by @xmath46 the density of quasiparticles at the bottom of the continuous spectrum , yielding @xmath287 . in the interval of densities @xmath288 ,
the impurity band may accommodate all the quasiparticles , while being too narrow to allow absorption of the field of frequency @xmath2 .
this low - absorption regime stretches up to the densities @xmath289 . at higher @xmath46 ,
the impurity band is broad enough to allow for intraband absorption ; compared to the benchmark @xmath290 dissipation level , eq .
( [ eq : sigmab - dirty3 ] ) yields @xmath285 smaller by a factor @xmath291 . upon further increase of @xmath46 , the impurity band merges with the continuum . deep in the ag regime , @xmath292 ,
one may use eq .
( [ dis3 ] ) for the estimate of dissipation ; compared to the @xmath290 dissipation level , eq .
( [ dis3 ] ) yields @xmath285 smaller by a factor @xmath293 .    to summarize the above estimates , in the selected region of frequencies , @xmath283 , and nonequilibrium quasiparticle densities , @xmath294 , magnetic impurities indeed reduce the dissipative part of conductivity due to the reduction of the dos available for the quasiparticle transitions .
finally , we reformulate the above results for the complex conductivity @xmath1 in terms of another important experimentally measurable quantity , the surface impedance:@xcite @xmath295 the impedance can be written as @xmath296 , where @xmath297 and @xmath298 are the surface resistance and reactance , respectively .
the complex conductivity @xmath299 is given by eqs .
( [ sigma ] ) and ( [ qkeld ] ) .
since we always consider the case of weak dissipation , @xmath300 , we also obtain @xmath301 .
the response of the superconductor is then almost purely reactive with @xmath302 under assumptions of sec .
[ sec : dissip ] , we can calculate @xmath303 in the simplest zero - temperature bcs model without magnetic scattering ( since corrections due to low temperature and weak magnetic scattering , @xmath304 , are negligible ) . at the same time , the new qualitative feature arising due to finite @xmath112 ( calculated in sec .
[ sec : dissip ] ) , is the dissipation , i.e. , small but nonzero @xmath297 : @xmath305    in the near - threshold situation of sec .
[ sec : dissip_zerot ] , we are interested in frequencies @xmath306 , then@xcite @xmath307 , and we find @xmath308    in the small - frequency limit considered in secs .  [ sec : ag ] and  [ sec : impband ] , at @xmath309 , we find@xcite @xmath310 and finally obtain @xmath311    thus , with the help of eqs .
( [ rth ] ) and ( [ rlow ] ) , the results for @xmath110 ( the dissipative part of conductivity ) calculated in sec .
[ sec : dissip ] are directly translated into the surface resistance .
we have employed the quasiclassical approach to study dissipation in superconductors due to small concentration of magnetic impurities .
the superconductor was assumed to be in the dirty limit with respect to the potential ( nonmagnetic ) scattering .
the concentration of magnetic impurities was assumed to be small enough so that the gap suppression is weak . employing the extension ( proposed by marchetti and simons@xcite ) of the usadel equation beyond the born approximation for magnetic scattering , we have considered in a unified manner both the ag regime ( with a continuum of states above a gap ) and the limit of lowest magnetic impurities concentrations where the impurity band due to overlap of localized impurity states is formed below the edge of the continuum .
the deviation from the born limit was assumed to be finite but small , so that the subgap impurity states lie close to the continuum edge .
our main results refer to the limit where the temperature @xmath79 , the frequency @xmath2 , and the gap suppression @xmath227 are all much smaller than the bcs gap @xmath21 in the absence of magnetic impurities . at the same time
, the relation between @xmath79 , @xmath2 , and @xmath227 was assumed to be arbitrary , and we have obtained explicit analytical expressions in various limiting cases .
our results can be expressed in terms of quantities describing the response of superconductor to external electromagnetic field : the dissipative conductivity @xmath112 ( the real part of the complex conductivity @xmath0 ) and the surface resistance @xmath297 ( the real part of the surface impedance @xmath312 ) .    in the limit of small temperatures and small frequencies that we considered , at equilibrium ,
the dissipation is always proportional to the density of thermally excited quasiparticles and thus exponentially suppressed due to gapped character of the spectrum . on the other hand , fluctuations in positions of magnetic impurities
can lead to finite dos below the mean - field gap.@xcite investigating the effect of the `` tail '' states on the dissipation is beyond the scope of our work .
we only note that since the number of tail states is small ( compared to the number of the `` mean - field '' states considered in this paper ) , there is a wide region of applicability for our results at subgap temperatures @xmath313 , excluding only ultra - low temperatures .    _ note added . _ during preparation of this paper , we became aware of a preprint by kharitonov _
_ et al.__@xcite addressing the same topic with a complementary approach . in
that work the surface impedance of superconductors with magnetic impurities is evaluated numerically , with emphasis on the superconducting gapless regime having low energy states in the electron spectrum .
we are grateful to m.  h.  devoret , r.  j.  schoelkopf , m.  v.  feigelman , and j.  s.  meyer for helpful discussions . ya.v.f .  was supported by the rfbr ( grant no .
11 - 02 - 00077-a ) , the russian federal agency of education and science ( contract no .
p799 ) , and the program `` quantum physics of condensed matter '' of the ras .
this work was also supported by doe contract no .
de - fg02 - 08er46482 ( yale ) and the nanosciences foundation at grenoble , france .
thanks yale university for the hospitality .
his visit was made possible by the gift of victor and marina vekselberg to yale university .
in this appendix we demonstrate that the electron spectrum in the presence of magnetic impurities , regardless of their polarization , is insensitive to the nonmagnetic disorder . to see this we consider separately two different models . in sec .
[ a : clean ] we derive explicitly the equation determining the dos in clean superconductors with randomly oriented spins and we establish the direct relation with the result obtained in sec .  [
sec : dirty ] ( dirty superconductors , randomly oriented spins ) . in sec .
[ a : dirty ] we formulate the usadel equation ( dirty limit ) for superconductors with polarized magnetic impurities and find that the dos corresponds to the one obtained in sec .
[ sec : clean ] ( clean limit , polarized magnetic impurities ) .      in the clean limit , the hamiltonian describing quasiparticles near the edge of the spectrum and interacting with randomly oriented weak magnetic impurities is @xmath314 where @xmath315 is the vector of the pauli matrices acting in the spin space and @xmath316 is the ( random ) orientation of the spin of impurity @xmath54 .
[ the hamiltonian ( [ ham1 ] ) for polarized magnetic impurities corresponds to all @xmath317 . ]
the @xmath79-matrix approximation for this model yields the self - consistency equation @xmath318 for the momentum - integrated green function @xmath74 , which is now diagonal in spin space . at low concentration of weak magnetic impurities ( @xmath37 and @xmath319 ) and for energies close to the edge of the continuous spectrum ( i.e. , @xmath320 ) , eq .
( [ grandom ] ) agrees with eq .
( [ u_eq ] ) that was derived in the dirty limit provided one identifies @xmath321 [ see also eq .
( [ tetd ] ) ] .
in particular , the analysis of eq .
( [ grandom ] ) in the limit @xmath322 yields a narrow impurity band with @xmath323 for the total dos ( summed over quasiparticle spins ) and @xmath324 for the width [ which coincides with eq .
( [ wfin3 ] ) ] .      following the derivation of eq .
( [ usadel_mod ] ) by marchetti and simons,@xcite we find that the usadel equation for the case of polarized magnetic impurities is @xmath325 here @xmath326 parametrizes the quasiclassical green functions that are different for electrons with spins up and spins down . at low concentration of weak magnetic impurities ( @xmath37 and @xmath319 ) and for energies close to the edge of the continuous spectrum , eq .
( [ usadel_modup ] ) in a uniform superconducting state can be rewritten as @xmath327 in terms of the spin - resolved momentum - integrated green functions @xmath328 .
the equation for down spin is the same as eq .
( [ g ] ) derived in the clean limit , while the one for spin up does not produce an impurity band .
the usadel equation in a superconductor reads @xmath329 -i[\check\sigma_s , \check g]=0 , \label{eq : usadel}\\ \hat\delta = \begin{pmatrix } 0 & \delta \\ \delta^ * & 0 \end{pmatrix},\end{gathered}\ ] ] where the quasiclassical green s function @xmath330 is a matrix in the nambu ( pauli matrices @xmath331 ) and keldysh spaces ( pauli matrices @xmath332 ) , @xmath333 and the self - energy @xmath334 describes magnetic scattering.@xcite the covariant derivative @xmath335 $ ] depends on the vector potential @xmath336 ( we choose the gauge without the scalar potential ) .
time convolution is implicit in the matrix product .
the superconducting gap @xmath337 enters the matrix ( in the nambu space ) @xmath338 and solves the self - consistency equation @xmath339 where @xmath340 is the bcs pairing constant .      in the next order , @xmath344,\ ] ] so that @xmath345\right)\ ] ] in the same order
. if we choose the london gauge with @xmath346 ( accompanied by the requirement of @xmath347 at the surface ) , then the usadel equation in the first order reads @xmath348.\ ] ] the l.h.s .  of this equation
is a linear operator @xmath349 acting on @xmath350 , while the source term in the r.h.s .
is proportional to @xmath351 .
taking the first order also in the self - consistency equation ( [ eq : self ] ) , we see that our choice of the gauge leads to @xmath352 and @xmath353 .
( it is actually well known that @xmath354 in the london gauge.@xcite )    the current flowing in the superconductor reads @xmath355 ( where @xmath201 is the drude conductivity ) , which in the lowest order in @xmath118 reduces to @xmath356 making the fourier transform , we obtain @xmath357 with the response kernel [ we omit the @xmath358 superscript of the green function for brevity ]        the keldysh component of the green function can be expressed in terms of the retarded and advanced ones , while the advanced component can be written in terms of the retarded one , @xmath363 so that finally only the retarded component enters . in the case without the superconducting phase , it has the form @xmath364 [ @xmath183 and @xmath193 are the functions introduced in sec .  [ sec : dirty ] , see eq .
( [ gftheta ] ) ] and we obtain @xmath365 - \notag \\ - & \tanh\frac{e_+}{2 t } \bigl [ g^*(e_- ) { \mathop\mathrm{re}\nolimits}g(e_+ ) + i f^*(e_- ) { \mathop\mathrm{im}\nolimits}f(e_+ ) \bigr ] \biggr\ } , \label{qkeld}\end{aligned}\ ] ]    where @xmath366 [ note that the integrand in eq .
( [ qkeld ] ) is even with respect to @xmath151 ] .
the difference of this expression for the kernel from the results of refs .
is that our expression is written in terms of the _ quasiclassical _ green functions ( in the dirty limit ) that can be directly found from the usadel equation .
a similar expression was employed in ref .  , although it was formulated in terms of contour integrals in the plane of complex @xmath151 and the cuts in the plane were made from @xmath207 , the values of the ag gap ; thus , the ag regime with a single gap was explicitly assumed .
our expression ( [ qkeld ] ) is more general , applicable also in the case when the impurity band is present and in the case of gapless superconductivity .
defining the complex conductivity@xcite @xmath368 that determines the ( local ) response of current to electric field , @xmath369 we see that @xmath370 dissipation is determined by the real part of conductivity , @xmath371 , for which , according to eq .
( [ qkeld ] ) , we obtain eq .
( [ dissip ] ) given in the beginning of sec .
[ sec : dissip ] .
the exponential in @xmath59 factor of @xmath372 may be obtained by the standard means relating @xmath57 to the product of @xmath66 and its derivative at the midpoint between the impurities .
however , obtaining this way the preexponential factor is problematic because of the oscillatory nature of wave function ( [ wf ] ) .
the short - range nature of the potential forming @xmath66 allows us a full solution of the two - impurity problem .
the main difference of the problem considered here from the corresponding one for a single quantum particle in the field of two short - range potentials@xcite stems from the structure of wave function ( [ wf ] ) which stays finite at @xmath373 .
the usadel equation must be complemented by the self - consistency equation for the order parameter @xmath21 . as a result
, @xmath21 is suppressed by magnetic impurities .
we do not treat the self - consistency equation explicitly , assuming @xmath21 that we use to be the self - consistent order parameter corresponding to the given impurity parameters .
alternatively , eq .  ( [ u_ag ] ) can be rewritten as a fourth - order polynomial equation , and then an analytical solution is available , according to ref .
( see also the discussion of this solution in ref .  ) .
however , this formal solution is bulky and nontransparent , therefore rather than simplifying it in the corresponding limiting cases , we choose to consider those cases directly .            in sec .
va2a of ref .
, maki considers the `` low - temperature and low - frequency '' limit , which probably means that @xmath79 and @xmath2 are much smaller than both @xmath21 and @xmath227 .
the condition @xmath374 written above eq .
( 102 ) in ref .   is probably a misprint for @xmath375 .
this would mean that eq .
( 102 ) from ref .   is obtained at @xmath376 , and in the limit when the last inequality becomes strong , we obtain agreement with our eq .
( [ dis6 ] ) .
the self energy for magnetic impurities with randomly oriented spins reads@xcite @xmath377 for polarized magnetic impurities , the spin structure of the quasiclassical green function @xmath378 needs to be taken into account ( pauli matrices @xmath379 ) and the self energy reads @xmath380