Unified understanding of superconductivity and Mott transition in alkali-doped fullerides from first principles.

Alkali-doped fullerides A 3C60 (A = K, Rb, Cs) are surprising materials where conventional phonon-mediated superconductivity and unconventional Mott physics meet, leading to a remarkable phase diagram as a function of volume per C60 molecule. We address these materials with a state-of-the-art calculation, where we construct a realistic low-energy model from first principles without using a priori information other than the crystal structure and solve it with an accurate many-body theory. Remarkably, our scheme comprehensively reproduces the experimental phase diagram including the low-spin Mott-insulating phase next to the superconducting phase. More remarkably, the critical temperatures T c's calculated from first principles quantitatively reproduce the experimental values. The driving force behind the surprising phase diagram of A 3C60 is a subtle competition between Hund's coupling and Jahn-Teller phonons, which leads to an effectively inverted Hund's coupling. Our results establish that the fullerides are the first members of a novel class of molecular superconductors in which the multiorbital electronic correlations and phonons cooperate to reach high T c s-wave superconductivity.

INTRODUCTION

Alkali-doped fullerides are unique exotic s-wave superconductors with maximum critical temperature Tc ~ 40 K, which is extremely high for the modest width of about 0.5 eV of the bands relevant to superconductivity (–). Although the evidence for a phononic pairing mechanism is robust and may suggest a conventional explanation for the superconductivity in these materials (), the experimental observation of a Mott-insulating phase next to the superconducting phase in Cs3C60 changes the perspective because the same strong correlations responsible for the Mott state are also likely to influence the superconducting state (–). This raises the fundamental question of why s-wave superconductivity, which is induced by some attractive interaction, is so resilient to (or even enhanced by) the strong repulsive interactions (). The fingerprints of phonons observed in the Mott-insulating phase, such as the low-spin state (–) and the dynamical Jahn-Teller effect (), confirm a nontrivial cooperation between the electron correlations and the electron-phonon interactions.

A variety of scenarios have been proposed to explain the superconductivity (, –). However, even the most successful theories include some adjustable parameter and/or neglect, for example, the realistic band structure, the intermolecular Coulomb interactions, and Hund’s coupling. These assumptions have prevented these theories from being conclusive on the superconducting mechanism, which arises from a subtle energy balance between the Coulomb and the electron-phonon interactions in a multiorbital system. Therefore, to pin down the pairing mechanism, a nonempirical analysis is indispensable: If fully ab initio calculations quantitatively reproduced the phase diagram and Tc, they would provide precious microscopic information to identify unambiguously the superconducting mechanism, opening the path toward a predictive theory of superconductivity in correlated molecular materials.

Here, we carry out such nonempirical analysis using a combination of density functional theory (DFT) and dynamical mean-field theory (DMFT), an approach that accurately describes the properties of many correlated materials (). DMFT is particularly accurate for compounds with short-range interactions and a large coordination number z, such as A3C60, a highly symmetric three-dimensional system with z = 12.

RESULTS

Theoretical phase diagram and comparison with experiments

In Fig. 1, we show the band structure obtained within DFT for face-centered cubic (fcc) Cs3C60. Because the t1 bands crossing the Fermi level are well isolated from the other bands, we construct a lattice Hamiltonian consisting of the t1 electron and phonon degrees of freedom. The Hamiltonian is defined in the most general form:on the fcc lattice, where each site represents a C60 molecule. , U, g, and ω are the electron one-body Hamiltonian, Coulomb interaction, electron-phonon coupling, and phonon frequency, respectively. They all have indices of the Wannier orbitals i, j, i′, j′, spin σ, σ′, momentum k,q, and the phonons also have a branch index ν. We construct maximally localized Wannier orbitals () as the basis for making the hopping and interaction parameters as short-ranged as possible. We determine all of the above parameters by state-of-the-art ab initio techniques (Materials and Methods). In particular, the recently developed constrained density-functional perturbation theory (cDFPT) () has enabled calculations of the phonon-related terms ( and ω), where the effects of the high-energy bands are incorporated into the parameter values. To solve the ab initio model, we used the extended DMFT (E-DMFT) (), which combines the ability of DMFT to treat local interactions with the inclusion of nonlocal terms. To study the superconductivity, we introduced the anomalous Green’s function, whose integral over frequency gives superconducting order parameter, into the E-DMFT equations.

Fig. 1

Band structure of fcc Cs3C60 with .

The blue dotted curves represent the Wannier-interpolated band dispersion calculated from in Eq. 1.

Figure 2A shows the theoretical phase diagram as a function of temperature T and volume per C603− anion (). As a function of T and , we find a paramagnetic metal, a paramagnetic Mott insulator, and the s-wave superconducting phase where two electrons in the same orbital form the Cooper pairs. The three phases are characterized by a finite spectral weight at the Fermi level, a Mott gap, and nonzero anomalous Green’s function and order parameter, respectively. The metal-insulator transition is of the first order below 80 K, and a direct first-order superconductor-insulator transition takes place around 780 Å3 at low temperature. Between the blue solid line Vc2(T) and the black dotted line Vc1(T), we find both a metallic and an insulating solution. The first-order transition line Vc(T), where the free energies of the two solutions cross, is expected to be close to Vc2(T) at low temperature as in the pure multiorbital Hubbard model ().

Fig. 2

Theoretical and experimental phase diagrams.

(A) Phase diagram as a function of volume per C60 molecule and temperature, obtained with the DFT + E-DMFT. PM, PI, and SC denote the paramagnetic metal, the paramagnetic insulator, and the superconducting phase, respectively. Metallic and insulating E-DMFT solutions coexist between the blue solid line Vc2(T) and the black dotted line Vc1(T). The error bars for Tc originate from the statistical errors in the superconducting order parameters calculated with the quantum Monte Carlo method (see Section E in the Supplementary Materials for details). The error bars for Vc1(T) and Vc2(T) are half the interval of the volume grid in the calculation. (B) For comparison, the experimental phase diagram [adapted by Y. Kasahara from Figure 6 in Zadik et al. ()] is shown, where the region depicted in (A) corresponds to the area surrounded by the green dotted lines. AFI denotes the antiferromagnetic insulator.

In Fig. 2B, we reproduce the experimental phase diagram () for fcc A3C60 systems to highlight the impressive agreement with our calculations. In particular, the maximum Tc of ~28 K in theory is comparable to that of the experiment ~35 K (, ). The critical volumes and the slope of the metal-insulator transition line are also consistent with the experiment. It is remarkable that our fully ab initio calculation, without any empirical parameters, reproduces quantitatively the experimental phase diagram including the unconventional superconductivity and the Mott insulator.

Unusual intramolecular interactions as a key for unconventional physics

The success of our theory in reproducing the experimental phase diagram demonstrates the reliability of the whole scheme including the estimates of the interaction parameters. We are therefore in a position to disentangle the crucial ingredients leading to the proximity of the s-wave superconductivity and a Mott insulator in the same diagram. The key quantities are the effective intramolecular interactions (intraorbital and interorbital density-density interactions Ueff = U + U + Uph and U′eff = U′ + U′ + U′ph, and exchange term Jeff = J + Jph) between the t1 electrons, which are used as input interaction parameters for the E-DMFT calculation. Here, U, U′, and J are the intramolecular Coulomb interactions screened by the high-energy electrons, whereas U = U′ and J = 0 (Uph, U′ph, and Jph) represent dynamical screening contributions from the intermolecular Coulomb (electron-phonon) interactions. These dynamical screening effects make the effective interactions dependent on electron’s frequency ω.

The solid curves in Fig. 3 show Ueff(ω) and U′eff(ω), where we find U′eff(ω) > Ueff(ω) up to ω ~0.2 eV except a narrow region around ω = 0.1 eV. This remarkable inversion of the low-energy interactions is associated with the negative value of Jeff(0) discussed below because the relation U′eff(ω) ~ Ueff(ω) − 2Jeff(ω) holds. As is apparent from U > U′ and U + U > U′ + U′ (plotted by dotted and dashed curves in Fig. 3), it is the phonon contribution (Uph, U′ph, and Jph) that causes the inversion. In fact, Uph(ω) and U′ph(ω) have strongly ω-dependent structures for ω ≲ 0.2 eV because of the intramolecular Jahn-Teller phonons with the frequencies up to ~0.2 eV (), which are comparable to the t1 bandwidth ~0.5 eV.

Fig. 3

Frequency dependence of effective onsite interactions.

The effective intra- and interorbital interactions (Ueff = U + U + Uph and U′eff = U′ + U′ + U′ph, respectively) consist of the constrained random-phase approximation (cRPA) onsite Coulomb repulsion (U, U′), the dynamical screening from the off-site interactions (U = U′), and the phonon-mediated interactions (Uph, U′ph). The data are calculated for Cs3C60 with at 40 K. We assume the cRPA Coulomb interactions to be static, whose validity is substantiated in Section B in the Supplementary Materials. Inset: Frequency dependence of Ueff and U′eff along the Matsubara frequency axis.

Table 1 shows the static (ω = 0) values of the phonon-mediated interactions for various A3C60 compounds. For all the interactions, the phonon contribution has an opposite sign with respect to U, U′, and J. Whereas |Uph(0)| and |U′ph(0)| are much smaller than U and U′ ~1 eV, |Jph(0)| ~51 meV is remarkably larger than that of the positive Hund’s coupling J ~34 meV (). As a result, an effectively negative exchange interaction Jeff(0) = J + Jph(0) ~ –17 meV is realized. The unusual sign inversion of the exchange interaction (, , , ) results from the molecular nature of the localized orbitals, which yields a small J, combined with a strong coupling between the Jahn-Teller phonons and t1 electrons (), which enhances |Jph(0)|.

Table 1

Material dependence of the static part of the phonon-mediated interactions.

Uph(0), U′ph(0), and Jph(0) are the phonon-mediated onsite intraorbital, interorbital, and exchange interaction strengths between the t1 electrons at ω = 0 calculated with the cDFPT. The energy unit is eV. The numbers just after the material names denote the volume occupied per C603− anion in Å3.

	Uph(0)	U′ph(0)	Jph(0)
K3C60 (722)	−0.15	−0.053	−0.050
Rb3C60 (750)	−0.14	−0.042	−0.051
Cs3C60 (762)	−0.11	−0.013	−0.051
Cs3C60 (784)	−0.12	−0.022	−0.051
Cs3C60 (804)	−0.13	−0.031	−0.052

Microscopic mechanism of superconductivity

We find, through a detailed analysis (Section F in the Supplementary Materials), that the unusual multiorbital interactions indeed drives the exotic s-wave superconductivity: (i) U′eff > Ueff and Jeff < 0 around ω = 0 generate a singlet pair of electrons (leading to a Cooper pair) (), which sit on the same orbital rather than on different orbitals as shown by 〈nn〉 > 〈nn〉 in Fig. 4A. (ii) Jeff further enhances the pairing through a coherent tunneling of pairs between orbitals [the Suhl-Kondo mechanism (, )].

Fig. 4

Double occupancy, size of spin, weights of intramolecular configurations, spectral functions at 40 K, and schematic pictures of representative intramolecular configurations.

(A) Volume dependence of the double-occupancy D = 〈nn〉 (red), the interorbital interspin correlation 〈nn〉 (green), and the size S of the spin per molecule (blue). (B) Spectral functions of several fcc A3C60 systems at 40 K. For comparison, we show the DFT density of states for fcc K3C60 ( = 722 Å3) as the shaded area. (C) Weights of several onsite configurations appearing in the quantum Monte Carlo simulations. (210) [(111)] generically denotes the configurations of {n1,n2,n3} = {2,1,0}, {0,2,1}, {1,0,2}, {2,0,1}, {1,2,0}, {0,1,2} [{n1,n2,n3} = {1,1,1}], with n being the occupation of orbital i. N ≠ 3 with N = n1 + n2 + n3 denotes the configurations away from half filling. (D) Illustrative pictures for the (111) and (210) configurations. The up and down arrows indicate the up- and down-spin electrons, respectively.

As we have discussed above, these unusual multiorbital interactions are caused by phonons. In this sense, phonons are necessary to explain the superconductivity. However, strong electronic correlations also play a crucial role, marking a substantial difference from conventional phonon-driven superconductivity (–): The strong correlations heavily renormalize the electronic kinetic energy (from a bare value of ~0.5 eV) so that even the small difference (~33 to 37 meV) between U′eff and Ueff is sufficient to bind the electrons into intraorbital pairs () [that is, it helps the mechanism (i)]. As a result, Tc increases with the correlation strength (see Tc versus in Fig. 2A). This analysis clearly demonstrates that the superconductivity in alkali-doped fullerides is the result of a synergy between electron-phonon interactions and electronic correlations (, ).

A peculiar low-spin Mott insulator

Also, the normal state displays remarkable consequences of the cooperative effect of the various interactions. Figure 4A shows the volume dependence of the intraorbital double occupancy D and the size S of the spin per molecule at T = 40 K. Here, we follow the metallic solution, expected to be stable in almost the entire coexistence region (). In the Mott-insulating region with the Mott gap in the spectral function (Fig. 4B), S approaches 1/2 (Fig. 4A), which is of the low-spin state, in agreement with the experimental observation (–). Another feature characteristic to the present Mott transition is the increase of D with the increase of the correlation strength, in sharp contrast with standard Mott transitions, in which D is drastically reduced.

This phenomenology also descends from Ueff < U′eff and Jeff < 0, which prefer onsite low-spin configurations with two electrons on one orbital (Fig. 4D). In the metallic phase, the electron hopping tends instead to make all the different local configurations equally likely (Fig. 4C) (). As the correlation increases, the electrons gradually lose their kinetic energy, and hence, the six low-spin configurations ({n1,n2,n3} = {2,1,0}, {0,2,1}, {1,0,2}, {2,0,1}, {1,2,0}, {0,1,2}) become the majority (Fig. 4C), which leads to the decrease (increase) of S (D). Eventually, in the Mott-insulating phase, these configurations become predominant. Remarkably, the six (210) configurations, in each of which the orbital degeneracy is lifted (Fig. 4D), remain degenerate and no orbital ordering occurs. This unusual orbitally degenerate Mott state is consistent with the experimental observation of a dynamical Jahn-Teller effect ().

DISCUSSION

We have carried out a fully ab initio study of the phase diagram of alkali-doped fullerides, in which all the interactions of a multiorbital Hubbard model coupled with phonons are computed from first principles with only the information on the atomic positions. The model is then solved by means of E-DMFT. We identify that a negative Hund’s coupling arising from the electron-phonon interaction and the related condition Ueff < U′eff that favors intraorbital pairs (, , , ) underlie the existence of the superconducting state, its competition with the low-spin Mott state, and the whole phase diagram, which in turn remarkably reproduces the experimental observations.

This identifies a general condition under which a material can show a high-temperature electron-phonon superconducting phase that benefits from strong correlations (), and suggests that a wider family of molecular conductors can exhibit a similar physics as long as an inverted Hund’s exchange is stabilized by electron-phonon coupling. In view of the accuracy demonstrated for A3C60, we can apply the methodology of this article to design other similar materials, possibly tailoring their properties to maximize the transition temperature.

MATERIALS AND METHODS

We apply the DFT + DMFT method () to fcc A3C60. We start from the global band structure (Fig. 1) obtained by the DFT (see Section A in the Supplementary Materials for the calculation conditions). Considering that the low-energy properties are governed by the isolated t1 bands, we derive from first principles the lattice Hamiltonian (Eq. ) for the t1 electrons and the phonons.

As the basis for the lattice Hamiltonian, we construct the maximally localized Wannier orbitals () {ϕ} from the t1 bands. We calculate the transfer integrals by with the Kohn-Sham Hamiltonian . in Eq. is the Fourier transform of t(R), which reproduces the original DFT band dispersion well (Fig. 1). The Coulomb interaction in the low-energy model incorporates the screening effects from the high-energy bands, whereas the screening processes within the t1 bands are excluded to avoid duplication because the latter processes are taken into account when we solve the model. Such partially screened Coulomb parameters are calculated () with the cRPA (). Similarly, we calculate renormalized electron-phonon coupling and phonon frequency ω by applying the recently developed cDFPT ().

In the partition function for the lattice Hamiltonian (Eq. ) written in the coherent state path integral formalism, we can integrate out the phonon degrees of freedom. This results in an electronic model with the additional electron-electron interaction (Uph, U′ph, and Jph) mediated by phonons on top of the cRPA Coulomb interaction. The onsite phonon-mediated interaction (with Uph = V, U′ph = V, Jph = V = V) is given, on the Matsubara axis, by ():where Ω = 2nπT is the bosonic Matsubara frequency and with the number N of k-points. Here, the sum over the phonon branches runs from 1 to 189. In practice, we omit the lowest 9 branches [ν = 1 to 9, which correspond to the acoustic modes, the librations, and the alkali-ion vibrations at the octahedral sites ()]. This is justified because the coupling between these modes and the t1 electrons is tiny ().

We solve the derived lattice Hamiltonian by means of the E-DMFT (), which maps the Hamiltonian (Eq. ) onto a three-orbital single-impurity model embedded in noninteracting bath subject to self-consistent conditions. In the mapping, the E-DMFT incorporates the dynamical screening through the off-site interactions into the onsite interactions (U and U′) in a self-consistent way. To solve the impurity model, we use the continuous-time quantum Monte Carlo method based on the strong coupling expansion (). Further detail can be found in Sections C and D in the Supplementary Materials.

The frequency dependence of U and U′ in Fig. 3 is calculated with the Padé analytic continuation. We perform the analytic continuation onto ω + iη with η = 0.01 eV from the data along the Matsubara axis at 40 K for Cs3C60 with . The spectral functions in Fig. 4B are calculated with the maximum entropy method.