A bound on the superconducting transition temperature

It is notoriously difficult to make quantitative theoretical predictions of the superconducting transition temperature, Tc, either from first-principles or even from a knowledge of normal state properties. Ultimately, this reflects the fact that the energy scales involved in the superconducting state are extremely small in natural units, and that Tc depends exponentially on a subtle interplay between different interactions so that small uncertainties in microscopic processes can lead to order one effects on Tc. However, in some circumstances, it may be possible to determine (approximate) bounds on Tc. Here we propose such a bound for the conventional phonon-mediated mechanism of pairing with strongly retarded interactions, i.e. in the case in which ℏω¯≪EF\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbar \bar \omega \ll E_F$$\end{document}, where ω¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\bar \omega$$\end{document} is an appropriate characteristic phonon frequency and EF is the Fermi energy. Specifically, drawing on both empirical results (shown in Fig. 2 below) and recent results1 of determinant quantum Monte Carlo (DQMC) studies of the paradigmatic Holstein model, we propose that kBTc≤Amaxℏω¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$k_BT_c \le A_{\max }\hbar \bar \omega$$\end{document}, where Amax is a dimensionless number of order one that we estimate to be Amax ≈ 1/10. A bound for the superconducting transition temperature (Tc) is proposed for conventional superconductors with phonon-mediated pairing mechanism. I. Esterlis and S.A. Kivelson from Stanford University, and D. J. Scalapino from University of California at Santa Barbara, perform analysis on the Holstein model and propose a bound for Tc of superconductors where electron-phonon coupling (λ) plays a dominant role. They re-evaluate the Migdal-Eliashberg theory with the assumption of monotonic λ dependence of Tc and point out where the theory applies and where it doesn’t due to non-monotonic λ dependence. With this knowledge, they estimate a key dimensionless parameter, with which a bound for Tc is related. Experimentally, Tc is further related to Debye temperature, θD and a bound is proposed to be θD/10. This result adds insights into the enigmatic problem of predicting Tc.


INTRODUCTION
In general, T c is a function of the phonon spectrum, the electron spectrum, the electron-phonon coupling and the electron-electron interactions. Given that the bare interactions between electrons are repulsive, it is something of a miracle that the effective interactions at low energies can be effectively attractive, i.e. produce pairing. In conventional superconductors, this miracle is a consequence of the retarded character of the phonon-induced interactions, hω ( E F . (In unconventional superconductors, it is rather the strongk dependence of the screened but still largely repulsive effective interactions that allows pairing to occur.) As we are focusing here on the case of conventional superconductors, it is reasonable to write a general expression for T c as the product of a dimensional factor hω=k B times a dimensionless function of the various dimensionless parameters that characterize the particular system in question, k B T c ¼ hω Aðw; λ; μ; Þ; where w hω=E F , λ is the (conventionally defined) dimensionless electron-phonon coupling, μ is the dimensionless electron-electron repulsion, and … represents other things, such as the form of the electron dispersion (e.g. the ratio of second to first neighbor hopping matrix elements), the phonon dispersion, etc.
Given that in typical metals, the Fermi energy is large compared to all other energies, it is reasonable to evaluate A in the w → 0 limit. In this same limit, μ is highly renormalized downward, so that μ → 0; we will therefore consider A in the limit μ = 0. Thus, the most significant parametric dependence of A concerns λ.

RESULTS AND DISCUSSION
Various approximate theoretical treatments of this problem based on Migdal-Eliashberg (ME) theory have long served as the basis for the accepted wisdom on this subject. As nominally ME theory is valid so long as wλ ( 1, it is generally believed that it can be applied even to the strong-coupling regime in which w À1 ) λ ) 1. Moreover, various approximate evaluations of the resulting self-consistency equations produce an expression for A that is a monotonically increasing function of λ; if correct, this would imply that, barring other instabilities i.e. lattice instabilities or charge-density wave (CDW) formation, the larger λ the larger T c .
We have recently shown 1 , on the basis of exact DQMC studies of the Holstein model (defined below), that this expectation is incorrect. Specifically, for a given small value of w we find that the ME approximation is extremely accurate for λ smaller than a characteristic value λ Ã $ 1, while for λ>λ Ã (even if λ ( w À1 ) ME theory is both quantitatively and qualitatively incorrect. In particular, λ * marks a crossover to a strong-coupling regime characterized by bipolaron formation, growth of commensurate CDW correlations unrelated to Fermi surface nesting, and incipient phase separation.
The Holstein model consists of a single electronic band, and a non-dispersing optical phonon coupled by the most local possible interaction to the on-site electron density. In relating model calculations with experiment, it is important to distinguish "bare parameters" (i.e., the parameters that appear in the model, such as the bare phonon frequency, ω 0 and the bare dimensionless electron-phonon coupling, λ 0 ) from "physical" quantities, such as the actual phonon dispersion, ωq, and the renormalized electronphonon coupling, λ. The tendency to phonon softening with increasing λ 0 implies that ωq<ω 0 and λ > λ 0 . Indeed, in the context of ME theory, λ diverges upon approach to a lattice instability. The breakdown of ME theory we have identified is distinct from any such lattice instability, and occurs when λ and λ 0 are both of order one. The crossover at λ * is associated with a complete rearrangement of the important low energy degrees of freedom.
To quantify this crossover, we show in Fig. 1 the T → 0 occupancy of the single particle state at the bottom of the band (nk fork ¼0), measured in DQMC and computed within ME theory. This state is far below the Fermi energy and hence, for non-interacting electrons n0 ¼ 2. On the other hand, in the strong-coupling polaronic limit, electrons are essentially localized on a lattice site, so n0 approaches the average electron density per site. For 0<λ ( 1 the electronic spectrum is perturbatively rearranged and n0 is slightly depressed. This behavior is apparent in both the ME approximation and from the DQMC results for λ<λ Ã . However, for λ>λ Ã , the DQMC results show a rapid decrease of n0, consistent with a crossover to the polaronic limit. Moreover, even though we are not able to directly compute T c (due to the difficulty in obtaining convergence of the DQMC results at low temperatures), by a series of indirect arguments, we 1 inferred that T c is maximal near the point of this crossover, λ ¼ λ Ã , and decreases dramatically when λ is either decreased or increased further. This leads us to the conclusion that there is a well-defined maximal value Max½A A max ¼ Að0; λ Ã ; 0; Þ. It turns out that a careful numerical evaluation of the full (k and ω dependent) self-consistent ME equations also leads to a nonmonotonic behavior of A, which leads to a vanishing T c for λ > λ * . However, in contrast to what is found in the DQMC, in the ME treatment the depression of T c for λ > λ * is associated with the onset of a competing incommensurate CDW order. This distinction is important, since if it were only the competition with CDW order that prevented high T c s, one could "engineer" interactions 2 that suppress CDW order so as to enhance T c .
Combining the results from the ME theory (where valid) with the DQMC results, we obtained estimates of A max for the Holstein model on the square lattice. For the range of parameters we have explored, the highest value of T c we have inferred is 0.08 times the bare phonon frequency, but because significant phonon softening occurs for λ $ λ Ã , this value of T c is 0.12 times the maximal renormalized phonon frequency. Many physically realistic generalizations of this model are possible-either by modifying the lattice structure, the electron band structure (further neighbor hopping matrix elements), the number of phonon modes and their dispersion, and the structure of the electron-phonon couplingall features represented by the … in A. There is no reason that the value of A max (…) obtained by optimizing with respect to λ, should not depend somewhat on these various features, although we have already found it to be relatively insensitive to small changes of the band structure. Still, it is an interesting exercise (which we are currently undertaking) to determine what microscopic features of the electron-phonon problem can increase A max .
For now, however, we will adopt an estimate of A max ≈ 1/10 as suggested by results for the simple Holstein model, and see how it compares with experiment in real materials.
Experimental determination of k B T c = hω In Fig. 2 we plot the superconducting T c vs. the Debye temperature, Θ D for various elemental superconductors and compounds for which data are available. In most of the data shown, Θ D is computed on the basis of the measured low temperature lattice contribution to the specific heat and the number of atoms per unit cell in the crystal structure, while in others it is inferred from e.g. low temperature resistivity. It thus represents a specific, unambiguously defined (although somewhat crude) estimate of the characteristic maximal phonon frequency hω $ k B Θ D . Also shown in the figure is the proposed bound T c ≤ Θ D /10.
Not only do we see that the bound is satisfied by all the data we have found (which is not, unfortunately, an exhaustive set), but in some cases the materials come quite close to saturating the bound, meaning that the bound may have some real significance. Specifically, for Pb (T c = 7.2 K), Nb (T c = 9.25 K), and Hg (T c = 4.15 K), three elemental superconductors known for their relatively strong electron-phonon couplings, T c /Θ D takes on the values 0.069, 0.034, and 0.058, respectively. The A-15 family of old fashion "high temperature superconductors," Nb 3 Sn (T c = 17.9 K), Nb 3 Ga (T c = 19.8 K), and Nb 3 Ge (T c = 21.8 K) have T c /Θ D equal to 0.066, 0.071, and 0.072, respectively.
At ambient pressure, the highest temperature superconductivity of a clearly conventional sort (with w ( 1) is MgB 2 , which does not appear in the figure because it has such a high Θ D = 884 K, and thus has T c /Θ D = 0.04. 3 This suggests that if a way can be found to increase the value of λ in this material, it could lead to as much as a factor of 2 enhancement of T c . The highest T c of all conventional superconductors is T c = 203 K in H 3 S at 155 GPa. 4 As far as we know, the Debye temperature has not been measured; however, if we identify ω with the largest phonon frequencies found in DFT calculations of the phonon band structure we obtain the estimate hω ¼ 0:23 eV. 5 If we accept this theoretical value, Fig. 1 Occupancy of thek ¼0 (band bottom) single-particle state as a function of the bare dimensionless electron-phonon coupling λ 0 (lower scale) and the renormalized coupling, λ (upper scale). (The definitions are as given in ref. 1 ). Solid dots are from DQMC simulations described in ref. 1 with ℏω 0 /E F = 0.1, and the data shown are at a temperature T ≈ E F /30, at which n0 has reached its asymptotic low T-value. The dashed line is the same quantity computed within ME theory. The inset shows the difference between the DQMC and ME results. The scale on the vertical axis has been chosen such that 0 is the value for non-interacting electrons and 1 is the value in the polaronic limit. (ρ = 0.8 is the average electron density per site.)  25 BKBO from ref. 6 . In most cases estimates of Θ D were obtained from lowtemperature-specific-heat measurements, in others it was estimated from e.g. low temperature resistivity. This dataset is far from exhaustive and we plan to expand on it in the future A bound on the superconducting transition temperature I Esterlis et al.
then k B T c = hω ¼ 0:08, i.e. it comes very close to saturating our bound. An especially interesting material from our perspective is Ba 1 −x K x BiO 3 (BKBO), which has an optimal T c = 32 K (T c /Θ D ≈ 0.09) for x ≈ 0.4. 6 Various features 7,8 of BKBO near optimal T c -softening of an optic phonon mode, diamagnetism above T c , proximity to a commensurate CDW phase-indicate this material may be especially relevant for studying the crossover from conventional superconductivity to strong-coupling, polaronic physics.

Further remarks
The idea of bounding T c is not new. For instance, an absolute bound for an electron-phonon mechanism T c < 30 K was proposed in ref. 9 . The loopholes in this analysis were recently summarized in ref. 10 . Conversely, a remarkable and highly influential analysis 11 of ME theory suggested that T c grows without bound with increasing λ, and can even be larger than hω=k B . It has long been recognized that this proposal was subject to the caveat that at large λ, a system may be prone to other instabilities, which could compete with superconductivity. It was shown 10 that the effect of phonon softening-still assuming the validity of ME theory-leads to a reduced prefactor in the T c expression, and thus to a bound on T c given by the bare phonon frequency, which is approached asymptotically as λ → ∞. Arguments for a bound within the context of ME theory have also been presented in refs. 12,13,14 . The breakdown of ME theory due to spontaneous breaking of translation symmetry has also been analyzed. 15,16 As the analysis in the present paper shares some features with these earlier studies, there are important ways in which it is different, both conceptually and practically. Our DQMC results show that deviations from the predictions of ME theory occur even in ranges of temperatures and λ in which no other form of order has arisen; in this range, ME theory always over-estimates the superconducting susceptibility. The non-monotonic λ dependence of T c that leads to our proposed bound is associated with a crossover at λ~λ * from a regime in which ME theory is extraordinarily accurate to a strong-coupling regime where the ME approximation breaks down entirely. This leads to a sharp drop of T c , even if the system in question has been carefully engineered to have no competing charge-ordering instabilities. The existence of an optimal λ * independent of competing instabilities, as far as we know, is inconsistent with all analysis based on ME theory, but not in conflict with any experimental observation.
Recently, it has been found that superconductivity in SrTiO 3 persists to such low electron densities that the Fermi energy is less than the typical phonon frequency, i.e. into a regime in which w ) 1. In this limit, it is far from clear that the effects of a bare repulsion, μ > 0, can still be neglected. However, as a problem in model physics, it is possible to ask whether a more general bound exists on A(w, λ, μ = 0,…). In the large w limit, it is possible to integrate out the phonons to obtain an instantaneous attractive interaction with a magnitude proportional to λ. In particular, the Holstein model in the limit w → ∞ maps onto the negative U Hubbard model with U $ ÀE F λ. It is well known that this model has an optimal T c at an intermediate value of jUj $ E F , where the maximal T c is bounded by a scale that is a fraction of E F . 17 Thus, at least in the artificial limit μ = 0, there exists a more general bound of the form Aðw; λ; μ ¼ 0; Þ Aðw; λ Ã ðwÞ; μ ¼ 0; Þ where Aðw; λ Ã ðwÞ; μ ¼ 0; Þ $ 1=w as w → ∞. There is another class of bounds on T c that can be inferred in a different manner; rather than taking normal state data, one can start with measured properties of the superconducting groundstate, from which one can attempt to bound the actual T c . Ideally, these bounds should apply to both conventional and unconventional superconductors.
For a simple BCS s-wave superconductor in the weak coupling limit, T c = BΔ 0 where Δ 0 is the zero-temperature gap and B = e γ / π ≈ 0.567. Strong-coupling effects, even in the context of BCS mean-field theory, have a tendency to decrease the value of B, and certainly fluctuational effects beyond mean-field theory will likewise decrease the value of B. Thus, taking into account the fact that the gap function can vary along the Fermi surface, one would generally expect that k B T c BjΔ max j; (2) where (with the case of unconventional superconductors in mind) Δ max is the largest value of the T = 0 gap on the Fermi surface. As far as we know, this inequality is satisfied (and often very nearly saturated) by all known crystalline superconductors. (In the presence of disorder it is, of course, possible to have gapless superconductors.) In a similar vein, a bound was proposed in ref. 18 based on the measured value of the zero-temperature superfluid stiffness, ρ s ðT ¼ 0Þ=m Ã (or equivalently from the zero-temperature value of the London penetration depth) where in a layered (quasi-2D) superconductor, a is the inter-layer separation while in a 3D superconductor, a ¼ ffiffiffi π p ξ 0 , where ξ 0 is the zero-temperature coherence length, and C ≈ 2.2. (2.2 is the ratio of T c to the zero-temperature phase stiffness of the 3D XY model on a cubic lattice.) Most superconductors satisfy this inequality but come nowhere near saturating it. However, certain high temperature superconductors, especially hole-doped cuprates, come within a factor of two or three of saturating this inequality.

DATA AVAILABILITY
The data that support the findings of this study are available from the corresponding author upon reasonable request.