Abstract

Within the generalized DMFT+Σ approach, we study disorder effects in the temperature dependence of paramagnetic critical magnetic field Hcp(T) for Hubbard model with attractive interaction. We consider the wide range of attraction potentials U—from the weak coupling limit, when superconductivity is described by BCS model, up to the limit of very strong coupling, when superconducting transition is related to Bose–Einstein condensation (BEC) of compact Cooper pairs. The growth of the coupling strength leads to the rapid growth of Hcp(T) at all temperatures. However, at low temperatures, paramagnetic critical magnetic field Hcp grows with U much more slowly, than the orbital critical field, and in BCS limit, the main contribution to the upper critical magnetic filed is of paramagnetic origin. The growth of the coupling strength also leads to the disappearance of the low temperature region of instability towards type I phase transition and Fulde–Ferrell–Larkin–Ovchinnikov (FFLO) phase, characteristic of BCS weak coupling limit. Disordering leads to the rapid drop of Hcp(T) in BCS weak coupling limit, while in BCS–BEC crossover region and BEC limit Hcp(T → 0) dependence on disorder is rather weak. Within DMFT+Σ approach, disorder influence on Hcp(T) is of universal nature at any coupling strength and is related only to disorder widening of the conduction band. In particular, this leads to the drop of the effective coupling strength with disorder, so that disordering restores the region of type I transition in the intermediate coupling region.

Notes

ACKNOWLEDGMENTS

This work was performed under the State Contract no. 0389-2014-0001 with partial support of RFBR Grant no. 17-02-00015 and the Program of Fundamental Research of the RAS Presidium no. 12 “Fundamental problems of high-temperature superconductivity.”

APPENDIX

EQUATION FOR PARAMAGNETIC CRITICAL MAGNETIC FIELD

In general case the Noziers–Schmitt-Rink approach [1] assumes, that corrections from strong pairing interaction significantly change the chemical potential of the system, but possible vertex corrections from this interaction in Cooper channel are irrelevant, so that to analyze Cooper instability we can use the weak coupling approximation (ladder approximation). In this approximation the condition of Cooper instability in disordered attractive Hubbard model is written as:

where \({{\tilde {N}}_{{0 \uparrow , \downarrow }}}\)(ε) is the “bare” (U = 0) density of states for different spin projections, “dressed” by impurity scattering. Spin splitting can be considered as an addition to chemical potential, so that introducing the “bare” density of states “dressed” by disorder in the absence of external magnetic field \({{\tilde {N}}_{0}}\)(ε), we obtain the final result for Cooper susceptibility:

In Eq. (13) energy ε is counted from the chemical potential level. If we count it from the middle of the conduction band we have to replace ε → ε – μ and the condition of Cooper instability (7) leads to the equation defining critical temperature depending on the external magnetic field, which gives the equation for paramagnetic critical magnetic filed (6). The chemical potential for different values of U and Δ should be determined from DMFT+Σ calculations, i.e. from the standard equation for electron number (band filling), which allows us to find Hcp for the wide range of model parameters, including the region of BCS–BEC crossover and the limit of strong coupling at different levels of disorder. This reflects the physical meaning of Nozieres–Scmitt-Rink approximation—in the weak coupling region the temperature of superconducting transition is controlled by the equation for Cooper instability (6), while in the strong coupling limit it is defined as the temperature of BEC, which is controlled by chemical potential. The joint solution of Eq. (6) and the equation for the chemical potential guarantees the correct interpolation for Hcp in the region of BCS–BEC crossover.