Electron-phonon Pairing Mechanism for Superconductivity in K_{3}C_{60} Fullerides
Daluram Yadav^{1}, Nishchhal Yadav^{1}, Sanjay K. Ghosh^{1,}
^{1}School of Studies in Physics, Vikram University, Ujjain, India
Journal of Materials Physics and Chemistry, 2013 1 (2), pp 9-18
DOI: 10.12691/jmpc-1-2-1
Received December 28, 2012; Revised February 02, 2013; Accepted February 28, 2013
Corresponding author: daluramyadav@gmail.com |
Cite This Article:
- Yadav, Daluram, Nishchhal Yadav, and Sanjay K. Ghosh. "Electron-phonon Pairing Mechanism for Superconductivity in K_{3}C_{60} Fullerides." Journal of Materials Physics and Chemistry 1.2 (2013): 9-18.
- Yadav, D. , Yadav, N. , & Ghosh, S. K. (2013). Electron-phonon Pairing Mechanism for Superconductivity in K_{3}C_{60} Fullerides. Journal of Materials Physics and Chemistry, 1(2), 9-18.
- Yadav, Daluram, Nishchhal Yadav, and Sanjay K. Ghosh. "Electron-phonon Pairing Mechanism for Superconductivity in K_{3}C_{60} Fullerides." Journal of Materials Physics and Chemistry 1, no. 2 (2013): 9-18.
Import into BibTeX | Import into EndNote | Import into RefMan | Import into RefWorks |
The electron pairing mechanism leading to superconducting state in K_{3}C_{60} fullerenes is explored within the framework of strong coupling theory. The intercage interactions between C_{60} cages and expansion of lattice due to the intercalation of K atoms are investigated using the nearest neighbor interactions. The renormalized Coulomb repulsive parameter μ* and the electron phonon coupling strength λ are obtained for the intermolecular phonon frequency ω_{er}, keeping in mind that the free electrons in lowest molecular orbital are coupled with intermolecular phonons. T_{c} is then estimated as 5 K, which is much lower as compared with the reported value. The electrons also couple with the intramolecular phonons then T_{c} enhances to 22 K. In this approach the carbon isotope effect exponent, the energy gap ratio, effect of pressure and volume on T_{c} is calculated which support the superconducting state.
fullerenes, intermolecular phonons, intramolecular phonons, transition temperature, carbon isotope effect, energy gap ratio
1. Introduction
Superconductivity in alkali metal intercalated Buckminster fullerenes with medium superconducting transition temperature (Tc ∼ 20 to 45K) [1] continue to provoke considerable interest. The attractive pairing mechanism with conventional phonon mediated electron pairing mechanisms in superconducting fullerides is a central point of several experimental studies. The normal state electronic and magnetic properties of alkali metal intercalated fullerides are still only partially understood and are reviewed at great length [2]. Superconductivity in light-element materials entered in a new era after the recent observation of a critical temperature as high as 38K in Cs_{3}C_{60} [3] and the birth of a new family of organic superconductors based on aromatic molecules, as doped picene [4], phenanthrene [5], coronene [6] and dibenzopentacene [7]. The unavoidable question about the “pairing glue” in the fullerides seems to have an answer at least for the standard members of the family, like K_{3}C_{60} or Rb_{3}C_{60}. For these compounds there is convincing evidence of an electron - phonon driven superconductivity, in which Jahn-Teller coupled intramolecular vibrations play the main role [8]. Nonetheless, the anomalous properties of “expanded” fullerides with large intermolecular distance, like ammoniated compounds [9] and most notably Cs_{3}C_{60} suggest that labeling these materials as standard superconductors described by the Bardeen-Cooper-Schrieffer (BCS) theory is at least hazardous. This compound is synthesized in two different crystal structures, the fcc which is common to most doped fullerides and a A15 structure, with a bipartite lattice [10]. A - 15 Cs_{3}C_{60} at ambient pressure is an anti ferromagnetic insulator with TN = 46 K and a spin-1/2 moment at each molecule [11], as opposed to K_{3}C_{60} and Rb_{3}C_{60}, which are superconductors in the same conditions.
The doping in C_{60} as A_{3}C_{60} (A = K, Rb, Cs) are at the centre of many novel phenomenon points to their technological and industrial importance. Neutron inelastic scattering measurements [12] and Raman spectroscopy [13] have been extensively used to probe the molecular vibration modes of C_{60} in pristine and doped C60 solids, to investigate aspects related to structural and electron-phonon interaction in superconducting fullerides. Neutron scattering [12] indicated that the vibration spectrum of K_{3}C_{60} may be divided into two regions, one of them belongs to the rotation of C_{60} molecule and the intermolecular vibrations in the range of 2.5 to 25meV and the others are the intramolecular mode with frequency 25 to 200meV. Raman scattering [13] yields the on-ball vibration modes (∼40meV - 0.2eV). In passing we note that, the phonon spectrum in AMIF consists of mode of C_{60} (26cm^{-}^{1}), the vibrations of K+ ions and the translational vibrations as a whole (∼ 45cm^{-1} - 120cm^{-}^{1}) and the on ball molecule vibrations (∼260cm^{-}^{1} -1900cm^{-1}) as revealed from lattice dynamic studies [14].
The carbon isotope effect αc = 0.30 ± 0.05 with the 100 % substitution of ^{6}C_{13} for ^{6}C_{12} [15] naturally favours the phonon mediated pairing. In the report of Fuhrer and coworkers, a reduced a, of about 0.21 ± 0.012 is documented [16]. We note that the Rb dopant isotope effect measurements verify that intermolecular phonons do not play a significant role in AMIF [17]. Later on, Burk and coworkers [18] have shown a rubidium isotope effect of αRb = -0.028 ± 0.036, a result which implies that the alkali-C60 optic phonons play at most a minor role in the pairing mechanism. However, a substantial intercalated isotope effect is expected in a situation with significant contribution from alkali- C_{60} optic phonons. We note that the value of αRB may hint at a similarity with the inverse hydrogen isotope effect in strongly anharmonic palladium hydride [19].
It is generally believed that the AMIF holds wave symmetry of the order parameter [20]. Tunneling measurements predicts a large energy gap ratio (β = 2Δ0/kBTc ) of 5.3 ± 0.2 for K_{3}C_{60} led to strong coupling regime [21]. Although tunneling measurements predicts a large magnitude of the reduced energy gap then the BCS theory (β = 3.53) support a low energy intermolecular phonons as the source of superconductivity. However, several other reports on energy gap i.e., infrared spectroscopy [22], muon spin relaxation measurement (μSR) [23] and absolute reflectivity data [24] provides energy gap value of 3.6 ∼ 4.0 close to the BCS weak coupling limit favors the intermolecular phonons. These values are also near those obtained from nuclear spin relaxation measurement [25]. Optical and tunneling studies [26] on Rb_{3}C_{60} led to β = 4.2 indicate that the superconductivity cannot be described in the weak coupling limit. Zhang and coworkers [27] first proposed that the doped alkali ion phonon mode produce a strong attraction for electrons on C_{60} using a Hubbard model to describe the superconductivity. In the theory of Varma et al. [28] the electron-phonon interaction in fullerides is induced by high frequency intramolecular vibration modes on a ball of C_{60} and it should be a key feature of superconductivity. The electron-phonon coupling parameter (λ_{ra}) is evaluated as 0.11.
The same follows from Schluter et al. [29]. The accurate results reflected from the self-consistent approach for intramolecular vibration [30]. However, the theory of Zhang and Guo led to the fact that the coupled K^{+} optical mode and the vibration of C_{60} mass centre will induce a strong electron-phonon interaction [31]. They have pointed that the value of electron-phonon coupling parameter for the intermolecular vibrations is too small to explain the superconductivity. Ivanov and Maruyama [32] have proposed a three square well model which is characterized by windows of low-frequency intermolecular phonons, high- frequency intramolecular and Coulomb energy and find that these interactions allows a coherent interpretation of superconducting fullerides properties. Furthermore, Alexandrov and Kabanov [33] have developed a non-adiabatic theory of super- conductivity.
The transition temperature Tc monotonically changes with lattice constant a. The value of a can be changed or by pressure or by substitution of different ions instead of starting. The pressure effect in fullerides is huge, dTc /dP is negative, dTc/dlna for different fullerides are close to each other. The direct measurements of the pressure dependence of the superconducting Tc for K_{3}C_{60} by Sparn and coworkers [34] noticed that Tc (= 19.3 K), decreases strongly with increasing pressure, dTc/dP = -7.8 K/GPa. Subsequently, measuring magnetic susceptibility [35] under hydrostatic pressure yields dTc /dP = -9.7 K/GPa with a Tc of 29.6 K for single-phase Rb_{3}C_{60}. In the theory of Chaban, the dependence of Tc on the pressure is argued to be connected with the chemical pressure effect in fullerides [36].
[33] Alexandrov, A. S. and Kabanov, V.V., “Theory of superconducting T_{c} of doped fullerenes,” Phys. Rev. B, 54. 3655. 1996.
It is widely believed that the high frequency (∼400 – 1400 cm^{-1}) intramolecular modes or electronic excitation would tend to favour the strong coupling theory while intermolecular modes (∼150 cm^{-1}) would require normally weak coupling theory. In the fullerenes K_{3}C_{60}, the intermolecular vibrations possess the energy in the range 2.5 to 25meV and the band structure calculations [37] have reported Fermi energy (εF) as 0.2 - 0.3eV. Looking to the various experimental reports [1,2 and 12-26, 34-37] and theoretical proposals [27,33], we have made our efforts to reveal the nature of pairing mechanism and physical properties leading to superconducting state, by considering the three dimensional K_{3}C_{60} as a diatomic lattice with C_{60} molecule and alkali metal (K) ions.
2. Methodology
The undoped C_{60} molecule is a pure carbon compounds and possesses a spherical structure. The nuclear cage of C60 have a diameter of 7.1 Å with 2 bond length (6-6 ring) 1.4 Å and (6-5 ring) 1.46 Å. The distance of the nearest approach in the molecular solid is 3.1 Å [38]. The C_{60} molecule crystallizes in the fcc phase with lattice parameter a = 14.20 Å and the distance between the centers of the nearest neighbour C60 cages is b = a/√2. Furthermore, in C_{60} molecule, out of 4 valence electrons on each carbon atom, 3 electrons participate in the 90 σ bonds along the edges of truncated isocahedron and the remaining 60 electrons of the molecule are in the orbital. In fact, the diameter of C_{60} molecule is large in such a way that the electron density is localized near the, surface of the sphere. Electronic energy band structure [39] reveals that for C_{60}, the lowest unoccupied molecular orbital is empty and can accommodate 6 electrons. The C_{60} cages they are regarded as rigid spheres and a bandwidth develop from the weak interactions between the nearest neighbor (n-n) C_{60} cages. In the following, the n-n C_{60} cages interactions are first considered and the force constants along with the phonon frequencies of the fcc lattice are then derived.
2.1. Dispersion Relation in C_{60}Consider the undoped C_{60} molecule, as a cube with N cages and volume is Na^{3}/4. When the volume is compressed from V to V+ dV the work done is:
(1) |
B being the Bulk modulus and a is the lattice parameter. For a cube V= a^{3} and ∂V/V= 3∂a/a. Keeping in mind the contribution of 6N bonds between n-n C_{60} cages, W is expressed as
(2) |
With k is the force constant for n-n cages. A comparison yields, k =3Ba /4.
A set of C60 molecule is treated as a three-dimensional monoatomic lattice with atomic mass and the position of cage as mk and u respectively. The lattice considered S unit cells numbered by an index l = 1, 2...S and for masses k (=1) [40]. The equations of motion along 110-direction are obtained from force ∂Φ / ∂u acting on the atom (k) with Φ as the potential energy of the crystal and follows
(3) |
also
(4) |
The equilibrium position of the l^{th} unit cell relative to an origin located at some atom is r(1) = l_{1}a +l_{2}b +l_{3}c, where l_{1}, l_{2}, l_{3} are integers and a, b, c are the primitive translational vectors. Using the plane wave solution as u_{j}(l) =A(q) exp i[q_{j}r(l)-ωt], with A as the amplitude, q is the wave-vector in 110 direction and ω(q) is the angular frequency, the dispersion relation reads
(5) |
At the zone boundary (q = π/b). The above is the longitudinal phonon mode, which is due to the rotational motion of the C_{60} molecule.
2.1.1. Dispersion relation in A_{3}C_{60}
The C_{60} molecule when chemically substituted with alkali metal (K) atoms, the semi conducting C_{60} becomes metallic, at the K_{3}C_{60} composition and it shows the superconducting nature. With the, doping, K atoms become fully ionized in the C_{60} crystal and give up their electrons to highly polarizable C_{60} molecule. These electrons go to de- localized lowest unoccupied molecular orbital, which can accommodate six electrons. The conduction band is empty for C_{60} crystal and for K_{3}C_{60} the conduction band is half filled up to the Fermi level. The K atoms were located in two non-equivalents tetrahedral and an octahedral position of the lattice. In doped fullerides, there exists 6N C_{60}-C_{60} bonds at a bond distance of a /, also 8N C_{60}-K bonds with bond distance as a/4 with a =14.28 Å for K_{3}C_{60}. In the true sense, a dopant entering in to parent lattice C_{60} tends to expand the lattice for creating a place for itself.
The transfer of charge will certainly result in to some contraction and the final lattice parameter of the doped material will be resultant of these effects. It is noticed that the lattice parameter for K_{3}C_{60} (a =14.28Å), Rb_{3}C_{60}^{.} (a =14.45 Å) and Cs_{3}C_{60 }(a =14.60 Å) which are more than the lattice parameter of undoped C_{60} namely_{.} a = 14.20 Å [1] and hence there is a resultant expansion. The total bond energy within the harmonic approximation is expressed as
(6) |
with, r_{CC} is the equilibrium C_{60}-C_{60} bond distance (10.04 Å) and r_{MC} denotes the metal-C_{60} bond distance and is r_{MC} = (r_{CC}/2) + r_{i}, with r_{i} = 1.63 Å is the ionic radius of K. A relationship in between the force constants K (C_{60}-C_{60}) and K' (M-C_{60}) is established when the total bond energy reach its maximum at the experimental value of lattice parameter a.
Treating K_{3}C_{60} as a three-dimensional (3D) diatomic lattice with atomic masses as m (M) and positions u (v) for C_{60}_{ }(K),_{ }the equations of motion follows
(7) |
(8) |
The plane wave solutions are used u_{i} (l) = A exp i[q_{j}r(l) - ωt] and v_{j} (l1) = B exp i [q_{j} r(l) - ωt] and k_{1} = 0.6 k'. Following harmonic approximation, all cubic and higher order terms are neglected and also the nearest neighbor interaction naturally demands the force constants and several other parameters in 3D [40]. Confining only for the longitudinal phonon modes of K_{3}C_{60} along 110-direction, the dispersion relation reads for small q value
(9) |
(10) |
Here, D_{1} = [k+(2k^{’}/3)]m^{-1}and D_{2} = [k’/3] M^{-1}. eqn. (9) is an intermolecular optical mode and acoustic characteristics are seen from eqn. (10) in the long-wavelength limit (q→0). The developed phonon modes are due to the displacement of C_{60} molecule or K^{+} ions in the K_{3}C_{60}.
2.1.2. Screening Parameter
[42] Eliashberg, G.M., “Interactions between electrons and lattice vibrations in a superconductor,” Sov. Phys. JETP, 11. 696. 1960.
The phonon mediated interaction is usually treated first and the direct Coulomb interaction is introduced in terms of pseudo potential, μ* [41,42]. In dealing with a high-T_{c}, superconductor, it is more realistic to first consider the large electron-electron interaction and later the small electron-phonon interaction strength to estimate T_{c}. The effect of screening of electrons is determined by the renormalized Coulomb repulsive parameter as [41].
(11) |
The cut-off frequency is set equal to the intermolecular optical phonon frequency. Usually, the Coulomb repulsive parameter is obtained
(12) |
Here the abbreviation Ξ^{2} = k_{s}^{2} /4k_{F}^{2} and k_{s}^{2} = 6πne^{2}/ε∞ ε_{F}.
2.1.3. Electron - Phonon Attractive Parameter
The attractive coupling constant can be calculated from the McMillan formulae [43] which is related through the density of states at fermi level N(ε_{F}) and the mean square electron-ion matrix and structure calculation of K_{3}C_{60} using full potential linear muffin-tin orbital [44] method suggest that the superconducting electrons are composed entirely Carbon 2p, the contribution of C 2s orbital is very small and the largest contributions to N (ε_{F}) are provided by the C (3) atoms which are the closest to the neighboring C_{60} molecule. The use of McMillan expression is appropriate in K_{3}C_{60} superconductors for two fold reasons: large density of states as well the more the mass difference of K and C_{60}. The electron-phonon coupling strength λ is
(13) |
<I^{2}> is a mean square electron-ion matrix element, M the molecular mass and <ω^{2}> is an averaged square molecular vibration frequency. The mean square electron-ion matrix element is
(14) |
(15) |
V (q) being the screened Coulomb potential. In terms of bare Coulomb potential V_{c} (q) one express, V (q) = V_{c} (q)/ε(q). The dielectric function ε (q,ω) is
(16) |
With V_{F} is the fermi velocity. The static dielectric function ε (q) in the long-wavelength limit is
(17) |
And the result for screened coulomb potential is
(18) |
(19) |
With is Ω is the volume of the cell and Ze denotes effective nuclear charge and is 7e. Thus, the mean square electron ion matrix element from Eq. (15) follows
(20) |
and the electron-phonon coupling strength for the intermolecular vibration frequency <ω_{er}^{2}> is
(21) |
The values of the renormalized Coulomb repulsive parameter and the attractive electron - phonon coupling constant are used while estimating the superconducting and normal state parameters.
3. Discussion and Analysis of Results
The realistic values of some physical parameters are employed following experimental data as follows. The C_{60} molecule is first analyzed using the lattice parameter a = 14.20 Å and B = 18 GPa [45]. The effective mass of C_{60} cage is 1.2 x 10^{-24} Kg and the force constant κ is obtained as 19.17 N/m. The estimated value of the longitudinal phonon frequency of the C_{60} molecule from eqn. (5) is 42.42 cm^{-1 }(5.2 meV). It is worth to mention that this phonon mode comes from 6N bonds between n-n cages of C_{60} molecule. Chemical substitution of alkali metal ion in C_{60} molecule introduces extra electrons which are accommodated in the conduction band which is half filled up to the fermi level and these conduction electrons distort the lattice of K_{3}C_{60}. The distortion leads to a net expansion as lattice parameter for K doped fullerides (14.28 Å) and parent fullerides is a = 14.20 Å. For the evaluation of force constant relationship the total bond energy is then minimized at the experimental value of a (= 14.28 Å) for K_{3}C_{60} to obtain, k^{’} as 0.17k. It appears that the magnitude of force constant for alkali metal doped fullerenes is larger than the undoped C_{60}. The dispersion relations for the intermolecular phonon modes originate due to the displacement of C_{60} molecule or K^{+} ions are then obtained.
The mass of three alkali metals (K) is used as 0.17 amu. The intermolecular alkali C_{60} optic mode (ω_{er}) is obtained as 7.8 meV at q→ 0. The wave vector q (≈ 2k_{F} sinθ) can therefore take maximum value up to 2k_{F}. For small value of q, the product qb is 10 to deduce the value of acoustic mode as 125.47 meV at 2k_{F}. The electron density n (= 4x3/a^{3}) is estimated as 4.12 x 10^{21} cm^{-3} from the lattice parameter a of K_{3}C_{60} and the band structure value of mass as 3.0 m_{e}_{ }is used [44]. The fermi energy is obtained as 0.3 eV which is consistent with the earlier data, ε_{F} = 0.2∼0.3 eV [37]. The other parameters of the electrons are the screening parameter k_{s} ≅ 0.9 Å^{-1}) and the fermi wave vector k_{F}( ≈ 0.5Å^{-1}) which are required for the determination of the coulomb repulsive parameter. We have used the value of ε∞ as 4.4 [46].
The renormalized Coulomb repulsive parameter μ* is estimated as 0.19 and is attributed to the fact that ε_{F} is higher at least an order of magnitude with the intermolecular alkali C_{60} optic phonon frequency. The electron-phonon coupling strength (λ) using eq. (21) is obtained as 1.2. We note that the dimensionless electron-phonon coupling strength as estimated by Gunnarson is about 0.5-1 [8]. Quite generally, the moderately large value of Coulomb repulsive parameter is appropriate for a narrow band width material on the metallic side of a metal-insulator transition. one can see that in dealing with K_{3}C_{60} fullerides the coupling strength (λ >1) and the renormalized coulomb repulsive parameter μ* = 0.19 for the characteristic phonon frequency (ω_{er} = 7.84 meV) is believed to distort the lattice leading to a superconducting state.
3.1. Transition TemperatureIn the regime λ >1 the strong coupling theory [47] applies and the transition temperature T_{c}, for K_{3}C_{60} is estimated following
(22) |
Where
(23) |
(24) |
The effective coupling constant λ_{eff} is deduced as 0.62 for the numerical function t(λ) value 1.07 with λ, = 1.2 and T_{c}^{er} is estimated as ≅ 5K which is much lower than the reported experimental data [1] of about 18K. The above approach with intermolecular phonons does not explain a medium T_{c}, in K_{3}C_{60} and hence it is important to look for the role of intramolecular phonons to assess the reported value.The role of collective excitation mechanism over conventional phonon mechanism in cuprates was earlier successfully explained with a generalized expression for T_{c} in a strong coupling theory [47]. The critical temperature T_{c }is obtained by both intermolecular and intramolecular phonons in fullerides as
(25) |
with
(26) |
λ_{ra} being the intramolecular coupling constant and is obtained as 0.4 from the scattering time (τ = 1.2 x 10^{-14} sec) as
(27) |
We use the reported value of ω_{ra}.= 1012 cm^{-1} from the earlier data [1], T_{c} is obtained as 22 K. The strength of coupling σ = λ_{er} / λ_{ra} is 3.0. From the analysis of superconducting T_{c}, it is natural to made two comments the high-energy intramolecular phonons are moderately coupled while to that to the low energy intermolecular phonons are strongly coupled to neighboring electrons and T_{c}, mainly arises from the significant contribution of intramolecular phonons. The dimensionless Coulomb repulsive parameter μ* is therefore believed to severely reduced and attributed to the scaling factor ln[ε_{F/}ω_{b}] appeared in denominator of the expression of μ*. As the fermi energy is larger at least an order of magnitude by intermolecular phonons, the retardation effects are expected to be very small. Generally, energy ω_{b} is presented as characterizing bosons (ω_{er}, or ω_{ra}), but opinion about its choice in screening is far from being unanimous. The only thing one must be sure of is that, if a boson is a candidate for superconductivity, one must expect μ* < μ ∼ 0.5 and of course λ - μ*> 0. To obtain some specific results, we obtain numerical results for the transition temperature with coupling parameters and are presented in the following plots for various conditions. In Figure 1, we show the result for T_{c} as a function of μ up to 0.3 for a set of parameters λ_{er} = 1.2 and λ_{ra} = 0.4. The analytic expression for T_{c} from eq. (22) clearly demonstrates that T_{c}, is strongly influenced by the Coulomb repulsive parameter and is higher for small values ofμ. For a set of coupling parameters μ* = 0.19, λ_{e}_{r} = 1.2 and λ_{ra} = 0.4 and T_{c}, of 22 K is calculated.
The present study is essentially based on the application of strong coupling theory with an additional pairing force, we have feeling that such studies describe the dependence of T_{c}, on different coupling strengths, which has not been tested carefully and is definitive to impose any constraints on theoretical approaches. An important feature of eq. (22) is displayed in Figure 2. We plot T_{c}, as a function of λ_{er}_{ }up to 2.0 for a set of parameters μ* = 0.19 and λ_{ra} = 0.4. It is noticed from the curve that T_{c}. increases initially with the enhanced value T_{c}. For the above set of parameters, if the present system is strongly coupled then one can achieve the higher values of T_{c}. The values of coupling parameters as λ_{er} (= 1.2), λ_{ra} (= 0.4) and μ* = 0.19 are quite reasonable in view of fullerides superconductors. We also investigate the effect of the intramolecular phonon coupling strength on T_{c}. In Figure 3 the variation of T_{c} for different values of λ_{ra}, varying from 0.0 to 0.7 is presented. In plotting this curve we take λ_{er}= 1.2 and μ* = 0.19.
We see from this graph that T_{c}, is highly sensitive to λ_{ra} even for moderate coupling. We now focus on the relationship in between λ_{er} and μ* within the present formalism for K_{3}C_{60} superconductors. We plot the variation of λ_{er} with μ* in Figure 4. For lower values of λ_{er}, yields unphysical values for a set of parameters (T_{c} = 22 K and λ_{ra}= 0.4).
Within the present approach the intermolecular phonons must be strongly coupled. Higher and positive μ* implies a constant repulsive interaction that is insufficient to create a superconducting state despite of the additional coupling of intramolecular phonons apart from electron-intermolecular phonon interactions.
3.1.1. Isotope Effect Coefficient
The isotope effect coefficient is
(28) |
From eq. (22) and (28) and the ω dependence of T_{c} is introduced in terms of μ* following the relation μ*dlnω = dlnμ, leads to the following result of isotope effect coefficient as
(29) |
Here ω is the average phonon frequency. The mass of carbon directly influences the intermolecular phonon frequency, as ω_{er} ∼[m M / (m + M]^{1/2} with m (M) mass of C_{60} (K) and a finite value of carbon isotope effect is expected. Equation (29) yields α ≅ 0.24 which is slightly lower than the reported value of α = 0.30 ± 0.06 with 100% substitution of ^{6}C_{13} for ^{6}C_{12} [15] and consistent with results of Fuhrer et al. [16]. The reduced inter fullerene mass [mM_{l}(m + M)]^{1/2} is definitely small in comparison with the intra fullerene mass [(m + M)/2], the carbon isotope exponent essentially causes the change in intermolecular phonon frequency and thus it is argued that although intramolecular phonons do play a significant role, however the participation of intermolecular phonons can not be ignored. Figure 5 shows the variation of estimated isotope effect from Eq. (29) with screening parameter. It is noticed that the BCS one half value is recovered for it μ* = 0.0. With the increase in μ* values the isotope effect exponent decreases for λ_{er} = 1.2 and λ_{ra} = 0.4. However, a reliable alkali metal isotope effect is not reported for K_{3}C_{60} system, but a negative α_{Rb} is documented [18] that necessarily ruled out the possibility of participation of intermolecular phonons. The justification lies in a fact that the ratio α_{c}/ακ_{ }depends on the weights of on ball carbon and carbon alkali metal dominated phonon modes that are sensitive to λ and μ*.
3.1.2. Energy Gap Parameter
The energy gap parameter β [30] is
(30) |
Where Δ (0) is the energy gap at zero temperature. The energy gap parameter provides a measurement of coupling strength for phonon mediated electron-electron interaction Substantial doubts arise on tunneling measurement data [21], which predicts significantly higher values of β in K_{3}C_{60} fullerides. The energy gap parameter (β ≅ 3.68) is estimated using ω_{er} as 63.25 cm^{-1} and T_{c} of about 5K from eq. (30). The magnitude of energy gap is larger than the BCS limit of 3.53 and thus indicates the strong coupling regime. On the other hand, the gap parameter when intramolecular phonons are considered in the pairing mechanism is estimated as 3.54 which is close to the BCS limit where the values ω_{r}_{a} = 1012 cm^{-1} and T_{c} = 22K are used. This implies that, the intramolecular phonons requires weak interactions with the conduction electrons and is consistent with the infrared spectroscopy [22], muon spin relaxation rate measurements [23], nuclear spin relaxation measurements [25] which strongly favors the participation of intramolecular phonons in the pairing mechanism. Furthermore, the absolute reflectivity measurements [24] on.K_{3}C_{60} reports β = 3.6 ∼ 4.0, which is similar to BCS weak coupling limit and the electron-phonon coupling strength cannot be as large as 2.
3.1.3. Pressure Effect on Transition Temperature
In order to analyze the pressure dependence of T_{c}, for K_{3}C_{60} system, we begin with the earlier developed Eqn. (25) according to which,
(31) |
Keeping in mind that C_{60} molecules are extremely rigid and hence pressure neglecting the pressure induced effects of intramolecular phonons,
(32) |
Henceforth, we rewrite eq. (31) as:
(33) |
Furthermore,
(34) |
Using eq. (22), we get the following results:
(35) |
Where
(36) |
We use the following
(37) |
And
(38) |
(39) |
Combining these, we write
(40) |
Here, in first and the third term are dependent on the variation of the N(ε_{F}) with pressure. In order to make progress, we have used the value B = 13.18 GPa and ε_{F} = 0.31 eV as documented in the previous section. The free electron gas result dlnN(ε_{F})/dP = (-2/3) κ_{T}, where κ_{T} (∼B^{-1}) is thermal compressibility. The calculated value of dInN (ε_{F})/dP is about - 5.058 % (Gpa)^{-1}. However, Kerkoud and coworkers reported a rate of about ∼ -10 % (Gpa)^{-1} from NMR technique [48] although of limited accuracy. From eq. (11) and (13), we write
(41) |
(42) |
Here the pressure dependence of the quantity <I^{2}> has been neglected. Using eq. (41) and (42) into eq. (40), we obtain:
(43) |
Here, we assume: dlnλ_{ra}/dP ≈ dlnN(ε_{F})/dP, dln<ω_{er}^{2}>/dP ≈ 2dln<ωer>/dP and dln<ω_{er}>/ dP ≈ dln<ω_{er}>/dP. The free electron estimate for dlnω_{p}/dP = 1/2B yields a value of about 3.79 % (GPa)^{-1}. We follow γ_{G} = B dln<ω>/dP to obtain dlnω_{er}/dP = 15.17% (GPa)^{-1}, using γ_{G} = +2 as earlier reported [35]. The parameters needed for further calculations are obtained as D_{11} =1.7 and D_{12} = -0.96 from the valus of λ = 1.2, μ =0.75 and μ* = 0.19. Using Eqn's (33-43) we find the value of dlnμ*/dP = - 0.78% (GPa)^{-1}, dlnλ_{er}/dP = -35.39 % (GPa)^{-1}, dlnλ_{ra}/dP = - 5.05 % (Gpa)^{-1}, dlnT_{c}^{er} /dP is -45.64% (Gpa)^{-1}, dlnT_{c}/dP = - 39.60% (Gpa)^{-1} and slope dT_{c}/dP= -8.7 K/Gpa which is consistent with the earlier reported value [34] of -7.8 K/GPa. Figure 6 shows the variation of T_{c} with pressure along with the experimental data.
Figure 6. Variation of T_{c} with Pressure. Closed circles represents the experimental data taken from Sparn et al. (1992) |
For fitting of the data we have employed the relation T_{c} (P) = T_{c} (0) exp [-0.3γp] using γ_{G }= +2 [35]. We follow the discussion and results of Spam et al. [34], who fit the data by using T_{c} (P) = T_{c} (0) exp [-γp] with γ = 0.44 ± 0.03 GPa^{-1}. We therefore reach the conclusion that the logarithmic derivative of the Coulomb pseudo potential with pressure, dlnμ*/dP, is smaller as compared to the other logarithmic derivatives involved in the expression of superconducting transition temperature as we dealt with. The variation of T_{c}, with pressure is thus dependent on the variation of the density of states at fermi level, the coulomb pseudo potential and the variation of the characteristic phonon frequency with pressure.
3.1.4. Volume Dependent Transition Temperature
We now explore our efforts in studying the volume derivative of T_{c}, or dT_{c}/dV. Such studies point to either weak or strong coupling of mediating bosons. We note that due to extremely rigidity of C_{60} the compression contracts the weak intermolecular bonds and high-energy intramolecular bonds are unaffected. With these considerations, the eq. (25) simplifies to the following relation,
(44) |
Leads to the following result
(45) |
With an assumption that intramolecular phonons are independent in any change of the volume [49],
(46) |
Following eq. (22), we write
(47) |
Where,
(48) |
Now, we obtain,
(49) |
Here, we define
(50) |
(51) |
(52) |
and
(53) |
Where
(54) |
Furthermore, the logarithmic volume derivative of the coupling parameter λ follows:
(55) |
Leads to the ratio
(56) |
With the values of B, dln N (ε_{F})/dP and dln <ω>/dP discussed earlier, we find Φ = 3.3 for K_{3}C_{60}. In passing we note that the value of Φ for AMIF is comparable to the value of Φ = 3.7 for usual metals as Al [49].
The physical significance of the volume derivative of T_{c} we estimate the parameters as D_{21} = 0.16, D_{22} = 70.6, D_{23} = -2.7 and D_{24} = -323 from the earlier mentioned values of various coupling strengths. Looking to eq. (49), dT_{c} /dV is being influenced by volume dependence of screening parameter, volume dependence of intermolecular vibration mode, volume dependence of intermolecular phonon coupling strength and the intramolecular phonon frequency as well various coupling parameters are independent of volume. We note that the Coulomb repulsion suppresses T_{c} as is being noticed from Figure 1. Henceforth, the contribution to dT_{c}/dV proportional to dμ*/dV is negative. We follow Crespi and Cohen [49] who argues that the logarithmic dependence of μ* on ε_{F} points to the fact that volume derivative of μ* is small in magnitude. A natural argument follows that the volume derivative of t(λ) and of μ* is small in comparison to dλ/dV, henceforth we neglect these two terms in eq. (49).
4. Conclusion
In this paper, we have devoted our efforts in revealing the attractive pairing mechanism and various physical parameters of alkali metal intercalated fullerenes superconductor within the framework of strong coupling theory. The particular system chosen for the study, because huge experimental data is available and consistent explanation is timely needed. Considering the three dimensional K_{3}C_{60} as a diatomic lattice of K and C_{60}, we derive expression for intermolecular (acoustic and optical) phonon modes. We have properly used the Bulk modulus as well as other structural parameters. The intermolecular optical phonon frequency obtained and used frequency of intramolecular phonons is consistent with the recent theoretical and experimental findings. We exercised the superconducting state parameters i.e. the transition temperature (T_{c}), the carbon isotope effect (α), the energy gap parameter (β) and the pressure as well volume effect on T_{c} describing the superconducting state in alkali metal intercalated fullerides superconductors.
The calculated value of T_{c}, α, β and dT_{c}/dP are consistent with the reported data. The deduced value of coupling strength for the intermolecular phonons (optical) favours the strong coupling mechanism. By understanding that intramolecular phonons indeed do play a significant role, we have obtained a steep increase in T_{c}. Unlikely the coupling strength for intramolecular phonons is an indicative of weak coupling theory. Various reports strongly favours the intramolecular vibrations as a key source of superconductivity and claim that intermolecular phonons will not induce the state, still we believe and argue that both the inter and intramolecular phonons participate in the superconducting state as well is essential for superconducting fullerides properties.
We argue that the reduced value of a essentially points to the large Coulomb repulsion apart from large electron-phonon coupling strength. Deduced results on α and β are consistent with the earlier reported data's. It is inferred that electron pairing with intermolecular phonons would lead to a strong coupling while to that the intramolecular phonons require weak coupling mechanism. Although several reports strongly favours the intramolecular vibrations as a key source of superconductivity and claim that intermolecular phonons will not induce the state, the present analysis reveals that both the inter- and intramolecular phonons participate in the superconducting state.
We also discuss the pressure dependence of dT_{c} /dP. We appeal from the model calculations that the pressure dependence of T_{c} demonstrates the metallic behaviour. The appreciation reveals that the developed expression yields a value of dT_{c}/dP close to the reported data. The variation of T_{c} with pressure is thus dependent on the variation of the fermi level density of states, the Coulomb pseudo potential and the variation of the phonon frequency with volume. The pressure effect in AMIF is huge and pressure derivative of T_{c}, is negative, while to that volume derivative of T_{c} is positive.
Finally, the electron pairing in fullerides superconductors and a consistent interpretation of the properties describing superconducting state as transition temperature, carbon isotope effect, energy gap parameter, pressure and volume effect on T_{c}, describing the superconducting state reveals the fact that both, the Coulomb and electron-phonon (inter- and intramolecular) are important in examining the unusual physical properties in alkali metal intercalated fullerides superconductors.
Acknowledgement
The author (Daluram Yadav) is gratefully acknowledge to Dr. N. Kaurav for valuable discussion.
References
[1] | Hebard, A.F., Rosseinsky, M.J., Haddon, R.C., Murphy, D.W., Glarum, S.H., Palstra, T.T.M., Ramirez, A.P. and Kortan, A.R., “Superconductivity at 18 K in potassium-doped C60,” Nature, 350. 600. 1991. [ CrossRef ] |
[2] | Forro, L. and Mihaly, L., “Electronic properties of doped fullerenes” Rep. Prog. Phys. 64. 649. 2001. [ CrossRef ] |
[3] | Ganin, A. Y., Takabayashi, Y., Khimyak, Y.Z., “Bulk superconductivity at 38 K in a molecular system,” Nature Mater., 7. 367-371. 2008. [ CrossRef ] [ PubMed ] |
[4] | Mitsuhashi, R., Suzuki, Y.,, Yamanari, Y., Mitamura, H., Kambe, T., Ikeda, N., Okamoto, H., Fujiwara, A., Yamaji, M., Kawasaki, N., Maniwa, Y. and Kubozono, Y., “Superconductivity in alkali metal-doped picene,” Nature, 464. 76-79. 2010. [ CrossRef ] [ PubMed ] |
[5] | Wang, X.F.,, Liu, R.H., Gui, Z., Xie, Y.L., Yan, Y.J., Ying, J.J., Luo, X.G. and Chen, X. H., “Superconductivity at 5K in alkali-metal doped phenantherene’’ Nat. Comm., 2. 507. 2011. [ CrossRef ] [ PubMed ] |
[6] | Kubozono, Y., Mitamura, H., Lee, X., He, X., Yamanari, Y., Takahashi, Y., Suzuki, Y., Kaji, Y., Eguchi, R., Akaike, K., Kambe, T., Okamoto, H., Fujiwara, A., Kato, T., Kosugi, T. and Aoki, H., “Metal-intercalated aromatic hydrocarbons: a new class of carbon-based superconductors” Phys. Chem. Chem. Phys., 13. 16476. 2011. [ CrossRef ] [ PubMed ] |
[7] | Xue, M., Cao, T., Wang, D., Wu, Y., Yang, H., Dong, X., He, J., Li, F. and Chen, G.F., “Superconductivity above 30 K in alkali metal doped hydrocarbon," Sci. Rep., 2, 389, 2012. [ CrossRef ] [ PubMed ] |
[8] | Gunnarson, O., “Superconductivity in fullerides,” Rev. Mod. Phys., 69. 575.1997. |
[9] | Ricco, M., Fumera, G., Shiroka, T., Ligabue, O., Bucci, C. and Bolzoni, F., “Metal to insulator evolution in (NH) XNaK2C60,” Phys. Rev. B, 68. 035102. 2003. [ CrossRef ] |
[10] | Ganin, A.Y., Takabayashi, Y., Jeglic, P., Arcon, D., Potocnik, A., Baker, P.J., Ohishi, Y., McDonald, M.T., Tzirakis, M.D., McLennan, A., Darling, G.R., Takata, M., Rosseinsky, M.J. and Prassides, K., “Polymorphism control of superconductivity and magnetism in Cs_{3}C_{60} close to the Mott transition,” Nature, 466. 221-225. 2010. [ CrossRef ] [ PubMed ] |
[11] | Takabayashi, Y., Ganin, A.Y., Jeglič, P., Arčon, D.,Takano, T., Iwasa, Y., Ohishi, Y., Takata, M., Takeshita, N., Prassides, K. and Rosseinsky, M. J., “The disorder-free non-BCS superconductor Cs_{3}C_{60} emerges from an antiferromagnetic insulator parent state,” Science, 323, 1585-1590, 2009. [ CrossRef ] [ PubMed ] |
[12] | Pintschovius, L., “Neutron studies of vibrations in fullerenes,” Rep. Prog. Phys., 59, 473, 1996. [ CrossRef ] |
[13] | Mitch, M.G., Chase, S. J. and Lannian, J. S., “Raman scattering and electron–phonon coupling in Rb_{x}C_{60},” Phys. Rev. Lett., 68, 883, 1992. [ CrossRef ] [ PubMed ] |
[14] | Belosludov, V.R. and Shpakov, V. P., “Lattice-dynamics study of K_{3}C_{60},” Mod. Phys. Lett., 6, 1209, 1991. [ CrossRef ] |
[15] | Chen, C.C. and Lieber, C.M., Isotope effect and superconductivity in metal doped C60,” Science, 259. 655. 1993. |
[16] | Fuhrer, S., Cherrey, K., Zettl, A. and Cohen, M. L., “Carbon isotope effect in single-crystal Rb_{3}C_{60},” Phys. Rev. Lett. 83. 404. 1999. [ CrossRef ] |
[17] | Ebbesen, T.W., Tsai, J. S., Tanigaki, K., Hiura, H., Shimakawa, Y., Kubo, Y., Hirosawa, T. and Mizuki, J., “Dopant isotope effect on superconductivity in Rb_{3}C_{60},” Physica C, 203. 163. 1992. [ CrossRef ] |
[18] | Burk, B., Crespi, V.H., Zettl, A. and Cohen, M. L., “Rubidium isotope effect in superconducting Rb_{3}C_{60},” Phys. Rev. Lett., 72. 3706. 1994. [ CrossRef ] [ PubMed ] |
[19] | Strizker, B. and Buckel, W., “Superconductivity in the palladium-hydrogen and the palladium deuterium systems,” Z. Phys. B, 257. 1. 1972. [ CrossRef ] |
[20] | Uemura, Y. J., Keren, A., Le, L.P., Luke, G.M., Sternlieb, B.J., Wu, W.D., Brewer, J.H., Whetten, R.L., Huang, S.M., Lin, S., Kaner, R.B., Diederich, F., Donovan, S., Gnmer, G. and Holczer, K., “Magnetic-field penetration depth in K_{3}C_{60} measured by muon spin relaxation,” Nature, 352. 605. 1991. [ CrossRef ] |
[21] | Zhang, Z., Chen, C.C. and Lieber, C.M., “Tunneling spectroscopy of M_{3}C_{60} superconductors: the energy gap, strong coupling, and superconductivity,” Science, 254. 1619. 1991. [ CrossRef ] [ PubMed ] |
[22] | DeGiorgi, L., Nicol, E.J., Klein, O., Gruner, G., Wachter, P., Huang, S.M., Wiley, J. and Kaner, R.B., “Optical properties of the alkali-metal-doped superconducting fullerenes: K3C60 and Rb3C60,” Phys. Rev. B., 49. 7012. 1994. [ CrossRef ] |
[23] | Kiefl, R.F., MacFarlane, W.A., Chow, K.H., Dunsiger, S., Duty, T.F., Johnston, T.M.S., Schneider, J.W., Sonier, J., Brad, L., Strongin, R.M., Fischer, J. E. and Smith III, A.B., “Coherence peak and superconducting energy gap in Rb_{3}C_{60} observed by muon spin relaxation,” Phys. Rev. Lett., 70. 3987. 1993. [ CrossRef ] [ PubMed ] |
[24] | DeGiorgi, L., Wachter, P., Gruner, G., Huang, S.M., Wiley, J. and Kaner, R. B., “Optical response of the superconducting state of K_{3}C_{60} and Rb_{3}C_{60},” Phys. Rev. Lett., 69. 2987. 1992. [ CrossRef ] [ PubMed ] |
[25] | Tycko, R., Dabbagh, G., Rosseinsky, M. J., Murphy, D.W., Ramirez, A.P. and Fleming, R.M., “Electronic properties of normal and superconducting alkali fullerides probed by 13_{C} nuclear magnetic resonance,” Phys. Rev. Lett., 68. 1912. 1992. |
[26] | Koller, D., Martin, M.C., Mihaly, L., Mihaly, G., Oszlanyi, G., Baumgartner, G. and Forro, L., “Energy gap in superconducting fullerides: optical and tunneling studies,” Phys. Rev. Lett., 77. 4082. 1996. [ CrossRef ] [ PubMed ] |
[27] | Zhang, F. C., Ogata, M. and Rice, T.M., “Attractive interaction and superconductivity for K_{3}C_{60},” Phys. Rev. Lett, 67. 3452. 1991. [ CrossRef ] [ PubMed ] |
[28] | Varma, C.M., Zaanen, J. and Raghavachari, K., “Superconductivity in the fullerenes,” Science, 254. 989. 1991. [ CrossRef ] [ PubMed ] |
[29] | Schluter, M., Lanoo, M., Needels, M., Baratt, G.A. and Tomanek, D., “Electron–phonon coupling and superconductivity in alkali-intercalated C_{60} solid,” Phys. Rev. Lett., 68. 526. 1992. [ CrossRef ] [ PubMed ] |
[30] | Kresin, V.Z., “Coupling strength in superconducting fullerenes,” Phys. Rev. B., 46, 1833, 1992. [ CrossRef ] |
[31] | Zhang, M.L. and Guo, H.Y., “Electron phonon interactions for superconductivity in alkali metal doped fullerenes,” Physica C, 227. 15-21. 1994. [ CrossRef ] |
[32] | Ivanov, V. and Maruyama, Y., “Disorder and phonon windows for superconductivity in doped fullerides,” Physica C, 247. 147. 1995. [ CrossRef ] |
[33] | Alexandrov, A. S. and Kabanov, V.V., “Theory of superconducting T_{c} of doped fullerenes,” Phys. Rev. B, 54. 3655. 1996. [ CrossRef ] |
[34] | Sparn, G., Thompson, J.D., Whetten, R.L., Huang, S.M., Kaner, R.B., Diederich, F., Gruner, G. and Holczer, K., “Pressure and field dependence of superconductivity in Rb_{3}C_{60},” Phys. Rev. Lett., 68. 1228. 1992. [ CrossRef ] [ PubMed ] |
[35] | Diederich, J., Gangopadhyay, A.K. and Schilling, J.S., “Pressure dependence of the electronic density of states and T_{c} in superconducting Rb_{3}C_{60},” Phys. Rev. B, 54. R9662. 1996. [ CrossRef ] |
[36] | Chaban, I.A., “Pressure effect in high-temperature superconductors and fullerides,” J. Supercond: Incorp. Nov. Mag., 15. 179. 2002. |
[37] | Satapthy, S., Antropov, V.O., Anderson, O. K., Jepsen, O., Gunnarson, O. and Liechtension, I., “Conduction-band structure of alkali-metal-doped C_{60},” Phys. Rev. B, 46. 1773. 1992. [ CrossRef ] |
[38] | Martins, J.L. and Troullier, N., “Structural and electronic properties of KnC60,” Phys. Rev. B, 46. 1766. 1992. [ CrossRef ] |
[39] | Erwin, S.C. and Pickett, W.E., “Theoretical normal-state transport properties of K_{3}C_{60},” Phys. Rev. B., 46. 14257. 1992. [ CrossRef ] |
[40] | Maraduddin, A.A., Montroll, E.W. and Weiss, G.H., Theory of lattice dynamics in the Harmonic approximation, Solid State Physics Suppl-3, Academic press, NewYork, 1963. |
[41] | Bogoliubov, N.N., Tolmachev, V.V. and Shirkov, D.D., “A New Method in the theory of Superconductivity,” Consultants Bureau, NewYork, 1959. |
[42] | Eliashberg, G.M., “Interactions between electrons and lattice vibrations in a superconductor,” Sov. Phys. JETP, 11. 696. 1960. |
[43] | McMillan, W.L., “Transition temperature of strong-coupled superconductors,” Phys. Rev., 167. 331. 1968. [ CrossRef ] |
[44] | Novikov, D.L., Gubanov, V.A. and Freeman, A.J., “Electronic structure, electron–phonon interaction and superconductivity in K_{3}C_{60}, Rb_{3}C_{60} and Cs_{3}C_{60},” Physica C, 191. 399. 1992. [ CrossRef ] |
[45] | Huffman, D.R., “Solid C_{60,}” Physics Today, 44. 22. 1991. [ CrossRef ] |
[46] | Hebard, A. F., Haddon, R.C., Fleming, R.M. and Kortan, A.R., Appl. Phys. Lett., 59. 2109. 1991. [ CrossRef ] |
[47] | Kresin, V.Z., “On the critical temperature for any strength of the electron–phonon coupling,” Phys. Lett. A, 122. 434. 1992. [ CrossRef ] |
[48] | Kerkoud, R., Auban Senzier, P., Jerome, D., Lambert, J.M., Zahab, A. and Bernier, P., Europhys. Lett., “13C Knight shift of the doped fullerenes K3C60,” 25. 379. 1994. |
[49] | Crespi, V.H. and Cohen, M.L., “Near constancy of the pressure dependence of Tc across families of organic and fullerene superconductors,” Phys. Rev. B, 53. 56. 1996. [ CrossRef ] |