A bound on superconducting $T_c$'s

It is notoriously difficult to make quantitative theoretical predictions of the superconducting $T_c$, 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 $T_c$ depends exponentially on a subtle interplay between different interactions so that small uncertainties in microscopic processes can lead to order 1 effects on $T_c$. However, in some circumstances, it may be possible to determine (approximate) bounds on $T_c$. 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 $\hbar\bar \omega \ll E_F$ where $\bar \omega$ is an appropriate characteristic phonon frequency and $E_F$ is the Fermi energy. Specifically, drawing on both empirical results (shown in Figure 2 below) and recent results[1] of determinant quantum Monte Carlo (DQMC) studies of the paradigmatic Holstein model, we propose that \begin{equation} k_B T_c \leq A_{max} \ \hbar \bar \omega \end{equation} where $A_{max}$ is a dimensionless number of order one that we estimate to be \begin{equation} A_{max} \approx 1/10. \end{equation}

It is notoriously difficult to make quantitative theoretical predictions of the superconducting T c , 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 T c depends exponentially on a subtle interplay between different interactions so that small uncertainties in microscopic processes can lead to order 1 effects on T c . However, in some circumstances, it may be possible to determine (approximate) bounds on T c . 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 ω E F whereω is an appropriate characteristic phonon frequency and E F is the Fermi energy. Specifically, drawing on both empirical results (shown in Figure 2 below) and recent results [1] of determinant quantum Monte Carlo (DQMC) studies of the paradigmatic Holstein model, we propose that where A max is a dimensionless number of order one that we estimate to be A max ≈ 1/10.
1 Theoretical analysis 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 is 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, ω E F . (In unconventional superconductors, it is rather the strong k dependence of the screened but still largely repulsive effective interactions that allows pairing to occur.) Since 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 ω/k B times a dimensionless function of the various dimensionless parameters that characterize the particular system in question where w ≡ ω/E F , λ is the (conventionally defined) dimensionless electronphonon 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; in any case, µ > 0 will only tend to reduce A. Thus, the most significant parametric dependence of A concerns λ.
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. Since 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 electronphonon coupling, λ 0 ) from "physical" quantities, such as the actual phonon dispersion, ω q , and the renormalized electron-phonon 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 1. The crossover at λ is associated with a complete rearrangement of the important low energy degrees of freedom.
To quantify this crossover we show in Figure 1 the T → 0 occupancy of the single particle state at the bottom of the band (n k for k = 0), measured in DQMC and computed within ME theory. This state is far below the Fermi energy and hence, for non-interacting electrons n 0 = 2. On the other hand, in the strong-coupling polaronic limit, electrons are essentially localized on a lattice site, so n 0 approaches the average electron density per site. For 0 < λ 1 the electronic spectrum is perturbatively rearranged and n 0 is slightly depressed.  Figure 1: Occupancy of the k = 0 (band bottom) single-particle state as function of the bare dimensionless electron-phonon coupling λ 0 (lower scale) and the renormalized coupling, λ (upper scale). (The definitions are as given in [1]). Solid dots are from DQMC simulations described in [1] with ω 0 /E F = 0.1, and the data shown are at a temperature T ≈ E F /30, at which n 0 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.) 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 n 0 , 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 It turns out that a careful numerical evaluation of the full ( k and ω dependent) self-consistent ME equations also leads to a non-monotonic 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 coupling -all 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.
2 Experimental determination of k B T c / ω In Fig. 2 we plot the superconducting T c vs. the Debye temperature, Θ D for various elemental superconductors and compounds for which data is 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 phonon frequency ω ∼ 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 [11]. This suggests that if a way can be found to increase the value of λ in this material, it Data on the various materials were obtained from the following references: elements from [3,4]; A-15s from [5], Skutterudites from [6,7,8]; CaC 6 from [9]; BKBO from [10]. In most cases estimates of Θ D were obtained from low-temperature specific-heat measurements, in others it was estimated from e.g. low temperature resistivity. This data set is far from exhaustive and we plan to expand on it in the future.
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 [12]. As far as we know, the Debye temperature has not been measured; however, if we identifȳ ω with the largest phonon frequencies found in DFT calculations of the phonon band structure we obtain the estimate ω = 0.23 eV [13]. If we accept this theoretical value, then k B T c / ω = 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 = 32K (T c /Θ D ≈ 0.09) for x ≈ 0.4 [10]. Various features [14,15] 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 < 30K was proposed in [16]. The loopholes in this analysis were recently summarized in [17]. Conversely, a remarkable and highly influential analysis [18] of ME theory suggested that T c grows without bound with increasing λ, and can even be larger than ω/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 [17] 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 [19] and [20].
While 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 a optimal T c at an intermediate value of |U | ∼ E F [21]. Thus, at least in the artificial limit µ = 0, there exists a more general bound of the form 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 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. [22] 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 = √ πξ 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.