Phonons in A3 C60 Lattice and Structural Dynamics

Abstract

The critical temperature ( _TC ) of superconductivity in _A3C60 compounds is generally lower smaller with alkali atoms (A). Furthermore _TC 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, _Cs3C60 , 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 _A3C60 compounds with A = Li, and Na have a low _TC or are not superconducting at all, and why _Cs3C60 is superconducting only at applied pressure and then with the highest _TC of all _C60 alkali fullerides. It is concluded that the density of states mechanism derived in the BCS model is in doubt. Nevertheless superconductivity in _A3C60 depends on electron-phonon coupling. The dominating phonon is the bond stretching _Ag phonon, a breathing phonon for the whole fullerene molecular ion.

Recommended by Dragan Mihailovic

1, Department of Chemistry, Chalmers University of Technology, 412 96 Göteborg, Sweden

Received 2 June 2009; Accepted 14 September 2009

1. Introduction

The discoveries of conductivity and superconductivity (SC) in _K3C60 and _Rb3C60 were great events in science during the 1990s [1-5]. SC was later discovered in a number of other _A3C60 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 _TC . Some of the _A3C60 compounds (A = Li, Na) with the smallest lattice constants are not superconducting [6-8]. Only _C60 fullerides have proven to be superconducting (SC).

Remarkably, the compound with the highest lattice constant does not follow this trend [9-11]. Disorder-free _Cs3C60 is an antiferromagnetic insulator at ambient pressure. However, already at the quite modest pressure of 3 kbar ( [approximate] 3000 atm), it turns into a superconductor (SC) with the highest _TC known for any fulleride [9-11]. At temperatures above _{TC Cs3C60} 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 _C60 fullerides as metals and apply the BCS model [12-15]. In the latter model, _TC increases if the density of states at the Fermi level increases. In the case of _A3C60 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 _TC . However, contrary to this deduction, _{Cs3 C60} is not superconducting at all.

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

SC _A3C60 is a "molecular metal". If so, one should be able to derive (Hubbard) U=0 for the solid, using molecular data. The molecular Hubbard U is defined by [figure omitted; refer to PDF] where I is the ionization energy and A is the electron affinity [12-15]. Δ^G0 is the adiabatic energy [19] for charge disproportionation according to: [figure omitted; refer to PDF] U is thus the corresponding vertical energy. Usually U>0 for a molecule or insulator, but it is still possible to have Δ^G0 =0 , that is, oxidation state degeneracy according to (2). The reason is that Δ^G0 corresponds to the optimum geometry of the products, in the present case for ^{_C60 2-} and ^{_C60 4-} , in (1). Δ^G0 <0 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 ^{_C60 3-} state [19]. Interaction with ^{_C60 3-} is maximized when Δ^G0 =0 . For _K3C60 , most estimations of U end up with values larger than 1 eV. This is too large to reach Δ^G0 =0 . Below we will find that if the lattice contribution to U is taken into account, it is possible to achieve Δ^G0 =0 .

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 H_g mode is important. In contradistinction, it will be argued here that the A_g 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) [19-22]. "Valence" here meant the charge on the _C60 molecule. The negative lattice contribution to Hubbard U is inversely proportional to the lattice constant R . Contraction of the lattice, for example by applied pressure [23-25] therefore leads to a lower U , a more negative Δ^G0 in (2), and reduced _TC . _TC in _Rb3C60 decreases almost linearly with applied pressure [25]. For _Rb3C60 Δ^G0 is negative already ambient pressure but for _Cs3C60 Δ^G0 is positive. This makes _Cs3C60 an antiferromagnet at ambient pressure, but applied pressure decreases the lattice constant enough to make Δ^G0 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 _C60 fullerides, including _Cs3C60 .

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 ANC60 N- where N is the number of electrons in LUMO. We will only be interested in the case with incompletely filled LUMO ( N<6 ). 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 2s2 2p2 configuration of carbon produces the states: 3 P, 1 D, 1 S. According to Hund's rule, the ground state has the highest possible spin, therefore 3 P. In the C60 anions the energy difference between high spin and low spin states is very small [26-28]. For C60 2- , the reason is that the exchange integrals are numerically small, of about the same size as the relevant integrals for the t1u 2 [arrow right] t2g 2 correlation effect [29-31].

The low spin ground state of ^{_C60 2-} , ^{_C60 3-} , and ^{_C60 4-} is a prerequisite for superconductivity according to the MV-3 model. If in addition (2) also holds with a negative Δ^G0 , the diamagnetism of ^{_C60 3-} 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 ^{_C60 4-} site to a ^{_C60 2-} site can take place without activation energy, or, in other words, that the electron pair is not trapped at the ^{_C60 4-} site. This has been discussed by Duscesas and the author, and it was found that the reorganization energy is only 0.24 eV for adding two electrons to ^{_C60 2-} [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 C₆₀ the shorter bonds are the thirty (3×20/2) 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 t_1u LUMO to form ^{_C60 -6} leads to a lengthening of the double bonds of C₆₀ to almost equal bond lengths over the whole molecular ion [32].

Figure 1: C₆₀ molecule. The most relevant electron-phonon interaction mode for EPT is the in-phase bond stretch for the short and strong bonds, drawn with thick lines.

[figure omitted; refer to PDF]

We are thus looking for a mode where the t_1u LUMO is stabilized on one molecule and at the same time destabilized on the other. The vibrational mode coupled to EPT between t_1u 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 C₆₀ the relevant phonon should be the A_g phonon at 1470 cm^-1 [12].

Gunnarsson pointed out that it is the size, not the mass of the alkali ion that determines _TC [12]. We therefore do not expect large contributions to the electron-phonon coupling due to the motion of the metal ions. The A_g breathing mode is very likely the phonon that couples to the electrons in SC. The H_g 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 t_1u 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 A₃ C₆₀ systems which are not superconducting at ambient pressure, NH₃ K₃ C₆₀ , NH₃ Rb₃ C₆₀ , and Cs₃ C₆₀ 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 U and Lattice Interactions

Estimation of U from molecules leads to a molecular U in the range 2.7-3.1 eV [12, 13, 38, 39]. For the solid there are experimental estimates of U in the range 0.8-1.27 eV [39-43]. This refers to a crystal where all site charges are the same, say Z . The charge disproportionated CDW phase, where the charges (Z-1) and (Z+1) alternate, is favored by an additional negative contribution according to the Born model. The Born radius R can be taken as the lattice constant. The decrease of the Born free energy is (atomic units) [figure omitted; refer to PDF] (neglecting the Born exponent). R is about 20 Bohr for fullerides. ΔG (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 U decreases. After this correction we obtain roughly U[approximate]0 for the _A3C60 fullerides. In conclusion, Hubbard U is less smaller with lattice constant.

High pressure in _A3C60 systems leads to a smaller lattice constant. Thus the CDW phase with ^{_C60 2-} and ^{_C60 4-} is favored compared to the SDW phase with ^{_C60 3-} 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 _A3C60 systems except _Cs3C60 .

Figure 2: PESs of the MV-3 model. Curvatures are exaggerated. λ is reorganization energy. The fat horizontal bars are vibronic energy levels. aa + bb and aa - bb denote approximate wave functions for the given energy levels. For _NH3K3C60 , _NH3Rb3C60 , and _Cs3C60 , the SDW-PES is lower than CDW-PES. Fat arrows show the behavior of CDW-PES at applied pressure. |Δ^G0 | is defined by (2).

[figure omitted; refer to PDF]

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 ( R ) 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 _Cs3C60 . In _Cs3C60 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 _Cs3C60 , changing the composition of the ground state. Interestingly, when the pressure is sufficient to change the phase into a SC phase, _TC is the largest one, corresponding to Δ^G0 =0 in Figure 2.

Figure 3: Free energies for the SDW and CDW states as a function of lattices constants in (3). Arrows at the symbol for the element show the superconducting gap. _Cs3C60 is in the antiferromagnetic region (AF) where the gap is a pseudogap.

[figure omitted; refer to PDF]

A number of experiments have been carried out on _NH3 -doped fullerides [16-18]. _NH3 intercalates without causing any other change than increasing the lattice constant [13]. As a function of volume/ ^{_C60 3-} the critical temperature _TC 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 a² and b² 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 A_{4 C60} is an insulator. A_{4 C60} with an unfilled shell-like _A3C60 has a different Hubbard U . Equation (3) corresponds to the following equation: [figure omitted; refer to PDF]

We now see that the right member consists of two odd-electron systems, while A_{4 C60} 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 Δ^G0 close to zero, while Δ^G0 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 A_{4 C60} systems and found that Mott-Hubbard behavior with Δ^G0 is larger than the one for the _A3C60 systems [45]. Iwasa applied reflection techniques to _A3C60 , verifying the metallic properties with Δ^G0 [approximate]0 [46]. Other experiments show a similar trend [47].

The reason why A_{4 C60} systems are nonSC may be large Hubbard U , as discussed by Rosseinsky [13]. Another important possibility may be that Jahn-Teller effect splits the t_1u 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 _TC on lattice constant in _A3C60 systems. The model is based on the structural dynamics that is connected to EPT and oxidation state degeneracy. Decrease of the lattice constant in _Cs3C60 due to increased pressure leads to appearance of SC. If the pressure is further increased, _TC increases at first, but eventually _Cs3C60 shares the fate of the other _C60 fullerides, so that _TC decreases with increased pressure. The fact that _TC 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 _TC , 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 _A3C60 fullerides it is a quite successful model but it does not predict the nonSC behavior for _Cs3C60 . It also does not explain why _NH3K3C60 and _NH3Rb3C60 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 A_{4 C60} systmes are not superconducting, as illustrated in the STM experiments by Wachowiak et al. [48] and further detailed by O'Shea [49]. The t_1u degeneracy is split and the lower MO is fully occupied, thereby making SC impossible. An interesting discussion of the Jahn-Teller effect in _C60 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 A_{4 C60} compounds may be the large Hubbard U for these compounds as compared to the _A3C60 compounds. An important task is to obtain an accurate value of U , as has in fact been the goal in previous review articles [12-15]. 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 Δ^G0 and U 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 _A3C60 is very likely a double-bond stretching phonon with symmetry _Ag .

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 _A3C60 . The BCS model on the other hand states that _TC depends on the density of states, but this does not hold for _Cs3C60 . Takabayashi et al. therefore concluded that SC in _Cs3C60 is of a different type, more similar to SC in cuprates with CuO₂ plane [10, 11, 51]. This is possible, not least considering the change of crystal structure from fcc in most _A3C60 fullerides to A15-bcc in _Cs3C60 . Of particular importance is that in _Cs3C60 the density of states at the Fermi level is higher than in the _A3C60 fcc fullerides and should have a higher _TC 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].

References

[1] A. F. Hebard, M. J. Rosseinsky, R. C. Haddon, "Superconductivity at 18 K in potassium-doped _C60 ," Nature , vol. 350, no. 6319, pp. 600-601, 1991.

[2] M. J. Rosseinsky, A. P. Ramirez, S. H. Glarum, "Superconductivity at 28 K in _RbxC60 ," Physical Review Letters , vol. 66, no. 21, pp. 2830-2832, 1991.

[3] R. M. Fleming, A. P. Ramirez, M. J. Rosseinsky, "Relation of structure and superconducting transition temperatures in _A3C60 ," Nature , vol. 352, no. 6338, pp. 787-788, 1991.

[4] L. S. Cox, R. A. Laskey, "DNA replication occurs at discrete sites in pseudonuclei assembled from purified DNA in vitro," Cell , vol. 66, no. 2, pp. 271-275, 1991.

[5] R. C. Haddon, A. F. Hebard, M. J. Rosseinsky, "Conducting films of _C60 and _C70 by alkali-metal doping," Nature , vol. 350, no. 6316, pp. 320-322, 1991.

[6] T. Yildirim, O. Zhou, J. E. Fischer, "Intercalation of sodium heteroclusters into the _C60 lattice," Nature , vol. 360, no. 6404, pp. 568-571, 1992.

[7] K. Prassides, C. Christides, I. M. Thomas, "Crystal structure, bonding, and phase transition of the superconducting _Na2CsC60 fulleride," Science , vol. 263, no. 5149, pp. 950-954, 1994.

[8] Y. Maniwa, T. Saito, A. Ohi, "Electronic states and superconductivity in alkali-intercalated fullerides: C 13 -NMR study in _Na2RbC60 , _Na2CsC60 , _K3C60 , _K2RbC60 , _K2CsC60 , _KRbCsC60 , _Rb2CsC60 and _RbCs2C60 ," Journal of the Physical Society of Japan , vol. 52, pp. 1139-1148, 1994.

[9] T. T. M. Palstra, O. Zhou, Y. Iwasa, P. E. Sulewski, R. M. Fleming, B. R. Zegarski, "Superconductivity at 40 K in cesium doped _C60 ," Solid State Communications , vol. 93, no. 4, pp. 327-330, 1995.

[10] A. Y. Ganin, Y. Takabayashi, Y. Z. Khimyak, "Bulk superconductivity at 38 K in a molecular system," Nature Materials , vol. 7, no. 5, pp. 367-371, 2008.

[11] Y. Takabayashi, A. Y. Ganin, P. Jeglic, "The disorder-free non-BCS superconductor _Cs3C60 emerges from an antiferromagnetic insulator parent state," Science , vol. 323, no. 5921, pp. 1585-1590, 2009.

[12] O. Gunnarsson, "Superconductivity in fullerides," Reviews of Modern Physics , vol. 69, no. 2, pp. 575-606, 1997.

[13] M. J. Rosseinsky, "Recent developments in the chemistry and physics of metal fullerides," Chemistry of Materials , vol. 10, no. 10, pp. 2665-2685, 1998.

[14] S. Margadonna, K. Prassides, "Recent advances in fullerene superconductivity," Journal of Solid State Chemistry , vol. 168, no. 2, pp. 639-652, 2002.

[15] Y. Iwasa, T. Takenobu, "Superconductivity, Mott-Hubbard states, and molecular orbital order in intercalated fullerides," Journal of Physics: Condensed Matter , vol. 15, no. 13, pp. R495-R519, 2003.

[16] M. J. Rosseinsky, D. W. Murphy, R. M. Fleming, O. Zhou, "Intercalation of ammonia into _K3C60 ," Nature , vol. 364, no. 6436, pp. 425-427, 1993.

[17] O. Zhou, R. M. Fleming, D. W. Murphy, "Increased transition temperature in superconducting _Na2CsC60 by intercalation of ammonia," Nature , vol. 362, no. 6419, pp. 433-435, 1993.

[18] O. Zhou, T. T. M. Palstra, Y. Iwasa, "Structural and electronic properties of _{(NH3 )xK3C60} ," Physical Review B , vol. 52, no. 1, pp. 483-489, 1995.

[19] S. Larsson, "Localization of electrons and excitations," Chemical Physics , vol. 326, no. 1, pp. 115-122, 2006.

[20] K. Prassides, P. N. Schatz, K. Y. Wong, P. Day, "Vibronic coupling model for mixed-valence compounds. Extension to two-site two-electron systems," Journal of Physical Chemistry , vol. 90, no. 22, pp. 5588-5597, 1986.

[21] K. Prassides, P. N. Schatz, "Vibronic coupling model for mixed-valence compounds. Extension to the multimode case," Journal of Physical Chemistry , vol. 93, no. 1, pp. 83-89, 1989.

[22] S. Larsson, "Electronic structure of planar superconducting systems: from finite to extended model," Chemical Physics , vol. 236, no. 1-3, pp. 133-150, 1998.

[23] G. Sparn, J. D. Thompson, S.-M. Huang, "Pressure dependence of superconductivity in single-phase _K3C60 ," Science , vol. 252, no. 5014, pp. 1829-1831, 1991.

[24] J. E. Schirber, L. Hansen, B. Morosin, J. E. Fischer, J. D. Jorgensen, G. H. Kwei, "Effect of pressure on the superconducting primitive cubic and "polymeric" phases of _Na2CsC60 ," Physica C , vol. 260, no. 3-4, pp. 173-176, 1996.

[25] J. Diederichs, J. S. Schilling, K. W. Herwig, W. B. Yelon, "Dependence of the superconducting transition temperature and lattice parameter on hydrostatic pressure for _Rb3C60 ," Journal of Physics and Chemistry of Solids , vol. 58, no. 1, pp. 123-132, 1997.

[26] P. Bhyrappa, P. Paul, J. Stinchcombe, P. D. W. Boyd, C. A. Reed, "Synthesis and electronic characterization of discrete buckminsterfulleride salts: _C60^2- and _C60^3- ," Journal of the American Chemical Society , vol. 115, no. 23, pp. 11004-11005, 1993.

[27] P. C. Trulove, R. T. Carlin, G. R. Eaton, S. S. Eaton, "Determination of the singlet-triplet energy separation for _C60^2- in DMSO by electron paramagnetic resonance," Journal of the American Chemical Society , vol. 117, no. 23, pp. 6265-6272, 1995.

[28] S. Nakazawa, K. Sato, D. Shiomi, "Electronic and molecular structures of _C60 -based polyanionic high-spin molecular clusters: direct spin identification and electron spin transient nutation spectroscopy for high-spin chemistry," Inorganica Chimica Acta , vol. 361, no. 14-15, pp. 4031-4037, 2008.

[29] S. Larsson, P. Bernier, J. E. Fischer, S. Roth, S. A. Solin, "Electron pair transfer in _A3C60 ," Chemical Physics of Intercalation II , vol. 305, of NATO Science Series, pp. 227-231, Plenum Press, New York, NY, USA, 1993.

[30] S. Larsson, L. Rodríguez Monge, "Correlation and pairing in C[_stack60^2n- ] ions. Superconductivity of alkali and alkaline earth compounds of _C60 ," International Journal of Quantum Chemistry , vol. 27, pp. 655-665, 1993.

[31] F. Negri, G. Orlandi, F. Zerbetto, "Low-lying electronic excited states of buckminsterfullerene anions," Journal of the American Chemical Society , vol. 114, no. 8, pp. 2909-2913, 1992.

[32] G. Duskesas, S. Larsson, "Bond lengths and reorganization energies in fullerenes and their ions," Theoretical Chemistry Accounts , vol. 97, no. 1-4, pp. 110-118, 1997.

[33] D. Reznik, L. Pintschovius, M. Ito, "Electron-phonon coupling reflecting dynamic charge inhomogeneity in copper oxide superconductors," Nature , vol. 440, no. 7088, pp. 1170-1173, 2006.

[34] T. Holstein, "Studies of polaron motion--part I: the molecular-crystal model," reprinted: Annals of Physics , vol. 281, pp. 706-724, 2000 Annals of Physics , vol. 8, no. 3, pp. 325-342, 1959.

[35] T. Holstein, "Studies of polaron motion--part II: the "small" polaron," Annals of Physics , vol. 8, no. 3, pp. 343-389, 1959.

[36] K. Prassides, S. Margadonna, D. Arcon, A. Lappas, H. Shimoda, Y. Iwasa, "Magnetic ordering in the ammoniated fulleride (_ND3 )_K3C60 ," Journal of the American Chemical Society , vol. 121, no. 48, pp. 11227-11228, 1999.

[37] J. Arvanitidis, K. Papagelis, T. Takenobu, "Antiferromagnetic ordering in the expanded (_NH3 )_Rb3C60 fulleride," Physica B , vol. 326, no. 1-4, pp. 572-576, 2003., [email protected]

[38] R. L. Martin, J. P. Ritchie, "Coulomb and exchange interactions in _C60^n- ," Physical Review B , vol. 48, no. 7, pp. 4845-4849, 1993.

[39] R. W. Lof, M. A. van Veenendaal, B. Koopmans, H. T. Jonkman, G. A. Sawatzky, "Band gap, excitons, and Coulomb interaction in solid _C60 ," Physical Review Letters , vol. 68, no. 26, pp. 3924-3927, 1992.

[40] P. A. Brühweiler, A. J. Maxwell, A. Nilsson, N. Mårtensson, O. Gunnarsson, "Auger and photoelectron study of the Hubbard U in _C60 , _K3C60 , and _K6C60 ," Physical Review B , vol. 48, no. 24, pp. 18296-18299, 1993.

[41] O. Gunnarsson, E. Koch, R. M. Martin, "Mott transition in degenerate Hubbard models: application to doped fullerenes," Physical Review B , vol. 54, no. 16, pp. R11026-R11029, 1996.

[42] M. R. Pederson, A. A. Quong, "Polarizabilities, charge states, and vibrational modes of isolated fullerene molecules," Physical Review B , vol. 46, no. 20, pp. 13584-13591, 1992.

[43] V. P. Antropov, O. Gunnarsson, O. Jepsen, "Coulomb integrals and model Hamiltonians for _C60 ," Physical Review B , vol. 46, no. 20, pp. 13647-13650, 1992.

[44] S. Larsson, A. Klimkans, L. Rodríguez-Monge, G. Duskesas, "Reorganization energies in organic π systems," Journal of Molecular Structure: Theochem , vol. 425, no. 1-2, pp. 155-159, 1998.

[45] M. Knupfer, J. Fink, "Mott-Hubbard-like behavior of the energy gap of _A4C60 (A=Na,K,Rb,Cs) and _Na10C60 ," Physical Review Letters , vol. 79, no. 14, pp. 2714-2717, 1997.

[46] Y. Iwasa, K. Tanaka, T. Yasuda, T. Koda, S. Koda, "Metallic reflection spectra of _K3C60 ," Physical Review Letters , vol. 69, no. 15, pp. 2284-2287, 1992.

[47] E. Sohmen, J. Fink, W. Krätchmer, "Electronic structure studies of undoped and n -type doped fullerene _C60 ," Europhysics Letters , vol. 17, pp. 51-55, 1992.

[48] A. Wachowiak, R. Yamachika, K. H. Khoo, "Applied physics: visualization of the molecular Jahn-Teller effect in an insulating _K4C60 monolayer," Science , vol. 310, no. 5747, pp. 468-470, 2005.

[49] J. N. O'Shea, "Molecular orbitals tell the story," Science , vol. 310, no. 5747, pp. 453-454, 2005.

[50] M. Capone, M. Fabrizio, C. Castellani, E. Tosatti, "Colloquiumml: modeling the unconventional superconducting properties of expanded _A3C60 fullerides," Reviews of Modern Physics , vol. 81, no. 2, pp. 943-958, 2009.

[51] G. R. Darling, A. Y. Ganin, M. J. Rosseinsky, Y. Takabayashi, K. Prassides, "Intermolecular overlap geometry gives two classes of fulleride superconductor: electronic structure of 38 K _TcCS3C60 ," Physical Review Letters , vol. 101, no. 13, 2008., [email protected]

[52] S. Larsson accepted for publication, International Journal of Quantum Chemistry