The critical temperature () of superconductivity in compounds is generally lower smaller with alkali atoms (A). Furthermore decreases with applied pressure. In the BCS model, these trends are explained by the lower density of states at the Fermi level for a decreased lattice constant (R). There is more than one counterexample, however, suggesting that BCS does not give the whole truth. The most important one is that the compound with the largest lattice constant, , is not superconducting at all at ambient pressure. In this paper we derive a novel model where a negative lattice contribution to Hubbard U, proportional to 1/R, is taken into account. It is possible to explain why compounds with A = Li, and Na have a low or are not superconducting at all, and why is superconducting only at applied pressure and then with the highest of all alkali fullerides. It is concluded that the density of states mechanism derived in the BCS model is in doubt. Nevertheless superconductivity in depends on electron-phonon coupling. The dominating phonon is the bond stretching phonon, a breathing phonon for the whole fullerene molecular ion.

1. Introduction

The discoveries of conductivity and superconductivity (SC) in and were great events in science during the 1990s [15]. SC was later discovered in a number of other compounds, where A stand for alkali atoms: Li, Na, K, Rb, or Cs. Generally, increase of the lattice constant leads to a higher critical temperature . Some of the compounds (A = Li, Na) with the smallest lattice constants are not superconducting [68]. Only fullerides have proven to be superconducting (SC).

Remarkably, the compound with the highest lattice constant does not follow this trend [911]. Disorder-free is an antiferromagnetic insulator at ambient pressure. However, already at the quite modest pressure of 3 kbar (3000 atm), it turns into a superconductor (SC) with the highest known for any fulleride [911]. At temperatures above passes directly into a semiconducting and antiferromagnetic phase. This system thus needs an explanation model which is not depending on the existence of free electrons, as in the BCS model. Takabayashi et al. further point out [11] that transfer from SC to antiferromagnetic phase appears to be “purely electronic”, thus seemingly explicable without resorting to nuclear dynamics.

There are a number of peculiarities in experimental data in general, as summarized by Gunnarsson [12], Rosseinsky [13], Margadonna and Prassides [14], and Iwasa and Takenobu [15]. A widely accepted approach is to treat the metallic alkali fullerides as metals and apply the BCS model [1215]. In the latter model, increases if the density of states at the Fermi level increases. In the case of the Fermi level is in a narrow band, so the conclusion is that if the band width further decreases, the density of states will increase and hence also . However, contrary to this deduction, is not superconducting at all.

Expansion of the lattice by intercalation of molecules leads to a higher density of states and to a higher . Lattices of several compounds have been expanded by intercalating and this improved the SC properties [1618] for some of them, but in and SC was lost and antiferromagnetism appeared instead.

SC is a “molecular metal”. If so, one should be able to derive (Hubbard) for the solid, using molecular data. The molecular Hubbard is defined by where is the ionization energy and is the electron affinity [1215]. is the adiabatic energy [19] for charge disproportionation according to: is thus the corresponding vertical energy. Usually for a molecule or insulator, but it is still possible to have , that is, oxidation state degeneracy according to (2). The reason is that corresponds to the optimum geometry of the products, in the present case for and , in (1). means that the sites are differently occupied in a pair of electrons [charge alternant phase or charge density wave (CDW) phase].

In the CDW phase, equivalent sites may exchange electron pairs, but there may be an activation barrier. The activation barrier may be removed due to interaction with the state [19]. Interaction with is maximized when . For most estimations of end up with values larger than 1 eV. This is too large to reach . Below we will find that if the lattice contribution to is taken into account, it is possible to achieve .

As another or complimentary cause of SC, one has often referred to the Jahn-Teller effect. The Jahn-Teller effect does not include the degeneracy of (2), however. A possible rational for involving the Jahn-Teller effect at electron transport is that it is different for different site charges, but that effect must be of minor importance in comparison with the degeneracy between site oxidation states, directly related to the Mott-Hubbard problem and the structural dynamics generated by nonconserved oxidation states.

On the basis of the Jahn-Teller effect, one has concluded that a single Hg mode is important. In contradistinction, it will be argued here that the Ag breathing mode is the very cause of structural dynamics leading to SC.

In this paper we will use a mixed-valence model for three oxidation states (MV-3) [1922]. “Valence” here meant the charge on the molecule. The negative lattice contribution to Hubbard is inversely proportional to the lattice constant . Contraction of the lattice, for example by applied pressure [2325] therefore leads to a lower , a more negative in (2), and reduced . in decreases almost linearly with applied pressure [25]. For is negative already ambient pressure but for is positive. This makes an antiferromagnet at ambient pressure, but applied pressure decreases the lattice constant enough to make negative and therefore an almost perfect SC.

Finally we will discuss electron-phonon interactions, which are of fundamental importance for SC, in MV-3 as well as BCS theory. We will see, using MV-3, that the same electron-phonon coupling is applicable to all fullerides, including .

2. Molecular Electronic and Vibrational Structure

It was early established that the valence electrons of the alkali atoms are transferred to the threefold degenerate, lowest unoccupied (LU) t1u molecular orbital (MO) of C60, thereby forming where is the number of electrons in LUMO. We will only be interested in the case with incompletely filled LUMO (). Due to the spherical symmetry, the electronic states for a single C60 may be compared to the multiplet states of a single atom in the second row with incompletely filled 2p shell. The 2s22p2 configuration of carbon produces the states: 3P, 1D, 1S. According to Hund’s rule, the ground state has the highest possible spin, therefore 3P. In the C60 anions the energy difference between high spin and low spin states is very small [2628]. For the reason is that the exchange integrals are numerically small, of about the same size as the relevant integrals for the correlation effect [2931].

The low spin ground state of , , and is a prerequisite for superconductivity according to the MV-3 model. If in addition (2) also holds with a negative , the diamagnetism of can be explained. Jahn-Teller effects and electron-phonon coupling are of minor importance in this case compared to electron correlation effects.

Another prerequisite for SC is that electron pair transfer (EPT) from a site to a site can take place without activation energy, or, in other words, that the electron pair is not trapped at the site. This has been discussed by Dušcesas and the author, and it was found that the reorganization energy is only 0.24 eV for adding two electrons to [32]. The reason for this is that the bond length change is very small.

Typical for systems is that occupation of the bonding MOs up to and including HOMO tends to shorten some of the CC bonds. The shortened bonds are said to have “double bond character”. In C60 the shorter bonds are the thirty bonds that are common to two hexagons (Figure 1). Occupation of LUMO lengthens these bonds, since LUMO tends to be antibonding in the double bonds. Full occupation of the t1u LUMO to form leads to a lengthening of the double bonds of C60 to almost equal bond lengths over the whole molecular ion [32].

We are thus looking for a mode where the t1u LUMO is stabilized on one molecule and at the same time destabilized on the other. The vibrational mode coupled to EPT between t1u orbitals on different sites is the one that shortens the double bonds and lengthens the single bonds. This is the breathing mode of the strong bonds in Figure 1.

Experimentally, breathing modes have been found to be involved in superconducting systems [33]. The theory [19, 29, 30] is equivalent to the Holstein diatomic model for a single electron [34, 35]. In neutral C60 the relevant phonon should be the Ag phonon at 1470 cm−1 [12].

Gunnarsson pointed out that it is the size, not the mass of the alkali ion that determines [12]. We therefore do not expect large contributions to the electron-phonon coupling due to the motion of the metal ions. The Ag breathing mode is very likely the phonon that couples to the electrons in SC. The Hg modes tend to increase some double bonds while decreasing the corresponding ones on the opposite side of the molecule, thus should be unimportant for EPT.

Next we will study EPT between two adjacent molecular ions A and B. The fundamental assumption in the chemical model for SC is that electrons and phonons are coupled in a motion that makes the electron pair oscillate between adjacent sites (in a time-dependent treatment). For this to happen, the activation barrier has to be sufficiently low. The exact wave function for the electron pair transferring between A and B may be complicated to write down, in particular since the relevant states are correlated and cannot be represented by single Slater determinants. If the t1u LUMO is called “a” on site A and “b” on the adjacent site B, there are three spin singlet states, where two electrons are localized on site A (aa), on site B (bb) or where the two electrons are localized on different sites (ab + ba). The symmetrized form of the former are aa + bb and aa − bb. ab + ba corresponds to the spin-coupled antiferromagnetic state [or spin density wave state (SDW)]. aa + bb and aa − bb are the necessary states for EPT, thus for SC.

It is interesting that the three A3C60 systems which are not superconducting at ambient pressure, NH3K3C60, NH3Rb3C60, and Cs3C60 with large lattice constant R, are all antiferromagnetic, thus with the ab + ba ground state [9, 36, 37]. The latter wave function cannot superconduct since the electrons are fixed one on each site.

3. Hubbard and Lattice Interactions

Estimation of from molecules leads to a molecular in the range 2.7–3.1 eV [12, 13, 38, 39]. For the solid there are experimental estimates of in the range 0.8–1.27 eV [3943]. This refers to a crystal where all site charges are the same, say . The charge disproportionated CDW phase, where the charges and alternate, is favored by an additional negative contribution according to the Born model. The Born radius can be taken as the lattice constant. The decrease of the Born free energy is (atomic units) (neglecting the Born exponent). is about 20 Bohr for fullerides. (Born) is thus equal to −1/20 Hartree or −1.3 eV. This lowers the energy difference between the SDW and CDW state , thus Hubbard decreases. After this correction we obtain roughly for the fullerides. In conclusion, Hubbard is less smaller with lattice constant.

High pressure in systems leads to a smaller lattice constant. Thus the CDW phase with and is favored compared to the SDW phase with sites. If on the other hand the lattice is expanded, the SDW phase wins and becomes the ground state. SC cannot occur since the ground state is mainly ab + ba.

Figure 2 shows a total energy potential energy surfaces (PES) at the molecular level, as a function of a collective breathing mode for the carbon atoms. In the case given in Figure 2 the outer minima are lower, and this corresponds to the CDW case, the case that applies to all systems except .

The important thing is that the inner PES, corresponding to the SDW, is present slightly above the CDW-PES. In most cases there is an interaction between CDW and SDW that opens a gap. The lower state is aa + bb mixed with ab + ba. This energy level is denoted as aa + bb in Figure 2. The aa − bb is the upper state. Provided the lower state has a major component of aa + bb, the gap is a true SC gap. If, however, the SDW-PES is below the CDW-PES, the ab + ba component will be the major one in the lower state. Since ab + ba is the wave function which put an electron on A and one on B, EPT is not possible and hence there is no SC.

We arrive at the picture sketched in Figure 3. As the lattice constant () increases the CDW states, denoted aa + bb and aa − bb in Figure 3 are increasing in energy due to the numerically smaller lattice energy contribution according to (3). aa + bb interacts with ab + ba, and this interaction is maximized at the crossing point (Figure 3). The ground state is mainly of aa + bb type except for . In the ground state has a very large component of ab + ba, which means that the antiferromagnetic phase wins over SC phase. High pressure decreases the lattice constant for , changing the composition of the ground state. Interestingly, when the pressure is sufficient to change the phase into a SC phase, is the largest one, corresponding to in Figure 2.

A number of experiments have been carried out on -doped fullerides [1618]. intercalates without causing any other change than increasing the lattice constant [13]. As a function of volume/ the critical temperature increases monotonically. If the increase in lattice constant is too large, the minimum of the SDW-PES appear below the minimum of the CDW-PES and the SC is lost. Instead the antiferromagnetic phase appears. This is very clear from Figure 3. As we move to the right in the diagram the ground state wave function gets increased character of ab + ba, equivalent to antiferromagnetic coupling.

The structural change at electron transfer is associated with a total energy lowering denoted as “reorganization energy” (). Large is directly related to high effective mass in the Holstein model [34, 35]. For fullerenes has proven to be comparatively small in the cases examined, or about 0.05 eV for single electron transfer [32]. As a comparison, is about 0.22 eV in polyacetylene [44]. The Jahn-Teller effect or modifications due to chemical bonding between the fullerene molecules may modify some distances, but this is usually of little importance for .

The coupling in the case of electron pair transfer is well defined and different from the coupling in single electron transfer and, of course, different from magnetic couplings. The most important difference is that the a2 and b2 states do not couple directly, but assume the presence of the SDW ab + ba state. In principle the coupling may be calculated in a quantum chemical calculation that includes two molecules.

An important problem is to understand why A4 is an insulator. A4 with an unfilled shell-like has a different Hubbard . Equation (3) corresponds to the following equation:

We now see that the right member consists of two odd-electron systems, while A4 itself has an even number of electrons, of course. Referring to the fact that systems with an even number of electrons generally have lower energy than systems with an odd number of electrons, it is quite simple to understand why (3) has a close to zero, while for (4) is decidedly positive. In any case this may be verified in quantum chemical calculations which involve two adjacent molecules.

Knupfer and Fink have used electron energy loss spectroscopy (EELS) for A4 systems and found that Mott-Hubbard behavior with is larger than the one for the systems [45]. Iwasa applied reflection techniques to , verifying the metallic properties with [46]. Other experiments show a similar trend [47].

The reason why A4 systems are nonSC may be large Hubbard , as discussed by Rosseinsky [13]. Another important possibility may be that Jahn-Teller effect splits the t1u level as shown in STM experiments [48]. The lower level is fully occupied [49], making charge transfer impossible.

4. Discussion

In this paper we have shown that a molecular model provides a satisfactory explanation of the dependence of on lattice constant in systems. The model is based on the structural dynamics that is connected to EPT and oxidation state degeneracy. Decrease of the lattice constant in due to increased pressure leads to appearance of SC. If the pressure is further increased, increases at first, but eventually shares the fate of the other fullerides, so that decreases with increased pressure. The fact that increases with increased pressure at first is a problem not only in the BCS model but also in the present model. Since this problem has to do with the properties of the ground state wave function [(ab + ba) versus (aa + bb) character], it would be unwise to dig deeper into this problem without accurate calculations.

In the BCS model too, decrease of the lattice constant leads to a lower , but in that case the reason is decrease of the density of states. The BCS mechanism is doubtful to say the least, as has been shown during the historical development of the SC phenomenon. In the fullerides it is a quite successful model but it does not predict the nonSC behavior for . It also does not explain why and are antiferromagnets. The phase change from SC to antiferromagnetic is driven by the composition of the ground state wave function, which in its turn depends on the lattice constant.

The Jahn-Teller effect is overrun by intermolecular bonding or at least strongly modified. In the molecular SC model the importance of the Jahn-Teller effect is superceded by the oxidation state degeneracy (disproportionation). However, the Jahn-Teller effect (combined with the chemical bonding effect) may well explain why A4 systmes are not superconducting, as illustrated in the STM experiments by Wachowiak et al. [48] and further detailed by O’Shea [49]. The t1u degeneracy is split and the lower MO is fully occupied, thereby making SC impossible. An interesting discussion of the Jahn-Teller effect in fullerenes and the dependence of lattice constant on superconductivity have recently been given by Capone et al. [50]. Some conclusions are different from those of this paper, however.

Another, or contributing reason for the absence of SC in A4 compounds may be the large Hubbard for these compounds as compared to the compounds. An important task is to obtain an accurate value of , as has in fact been the goal in previous review articles [1215]. The description here introduces the nuclear coordinate into the Mott-Hubbard problem. Thereby it becomes possible to discuss the vibronic states of importance in superconductivity. The distinction between and is important.

Phonons lower the energy in the SC phase where they interact with the electrons. The motion of the electron pair is correlated with the nuclear coordinates in the breathing mode. The vibration motions of the nuclei “pump” the electron pairs between the sites without expense of energy. This happens in a quantum mechanical manner and leads to an energy lowering that stabilizes the SC phase with an energy gap to the first excited state. The antiferromagnetic phase cannot behave in the same way, since the charge distribution is inert and almost independent of the exact positions of the nuclei. The relevant breathing phonon in the case of is very likely a double-bond stretching phonon with symmetry .

In this paper the difference in lattice energy between the CDW and SDW phases has been pointed out as a general explanation of SC in . The BCS model on the other hand states that depends on the density of states, but this does not hold for . Takabayashi et al. therefore concluded that SC in is of a different type, more similar to SC in cuprates with CuO2 plane [10, 11, 51]. This is possible, not least considering the change of crystal structure from fcc in most fullerides to A15-bcc in . Of particular importance is that in the density of states at the Fermi level is higher than in the fcc fullerides and should have a higher rather than not being SC at all.

In model used here change of structure leads to different Born model contributions to the CDW state in (3) shifting the onset pressure of SC. Since it is not known which onset pressure would apply to the fcc phase, a detailed discussion of this problem would not be successful at this stage. To complicate matters further, change of structure also leads to differences in coupling between the SDW and CDW states.

The author firmly believes that the BCS model is a kind of rough limit model for the transition between ordinary metals (e.g., alkali metals and coinage metals at normal pressure) and SC metals. The presence of free electrons is there guaranteed. The BCS model should not be used in cases where insulating phases play roles [52]. In any case the BCS model appears to be too crude a model to allow adequate parametrization. In the BCS model the binding energy of a pair of electrons determines whether SC appears, whereas in reality the binding energy depends on the presence of positive nuclei and the mobility on the structural reorganization energy versus coupling [52].