|Home||CHRNS||CHRNS Highlights||CHRNS Education||CHRNS-Devel||SiteMap|
GIANT ANHARMONICITY AND ELECTRON-PHONON MEDIATED SUPERCONDUCTIVITY
IN MgB2 AT 40 K
The recent surprise discovery of superconductivity in MgB2 at 39 K has stimulated a great deal of research on this intercalated graphite-like system (see Fig. 1). Sparked by this discovery, we set out to unlock the structural secrets and, in particular, to reveal the origin of the high Tc in MgB2: an electron-phonon or other exotic mechanism? To answer this fundamental question, we calculated Tc , its pressure dependence, and its isotope effect from the electronic band structure and lattice dynamics of MgB2 using density functional theory within the generalized gradient approximation [1-3].
Figure 2 shows that the features in the calculated phonon density of states (DOS) are in excellent agreement with the neutron data (GDOS), giving confidence that the calculations provide a sound description of the physical properties of the system. The DOS consists of two bands of phonons, one below 40 meV corresponding primarily to acoustic phonon modes, and one above 50 meV that mostly involves the boron motions. Inspection of the calculated phonon-dispersion curves that make up the high-energy band in the DOS reveals that the in-plane boron phonons (as depicted in the inset to Fig. 3) are very anharmonic. To demonstrate this, in Fig.3 we plot the total energy as boron atoms move according to one of these zone-center in-plane phonons with E2g symmetry. The potential indicates a very strong anharmonic term. Numerically solving the Schrödinger equation for this anharmonic potential yields a phonon energy of hw(E2g) = 74.5 meV, a 25 % enhancement over the harmonic value of 60.3 meV. This value is in good agreement with recent Raman measurements. The giant anharmonicity revealed gives the first hint that the in-plane modes are strongly coupled to the px1y bonding orbitals of boron, as shown schematically in Fig.1. This coupling is also evident from the splitting of the boron bands (red lines) with the E2g phonons (see Fig. 3). Note that the other bands are not affected by the E2g phonons.
The splitting of the boron bands, when averaged over the Fermi surface, gives an electron-phonon (EP) coupling constant 1. Using this value in the McMillan expression for Tc with w (E2g) and taking a typical value for the Coulomb repulsion µ* = 0.15, we obtain a Tc of 39.4 K, in excellent agreement with experiments. We also solved the Schrödinger equation for the potential shown in Fig. 3 for different boron masses and obtained w(E2g) = 291.8 M-0.575 and = 0.615l M 0.169, which yields the Tc -M curve shown in Fig. 4 and a boron isotope effect = 0.21, in good agreement with the experimental value of 0.26 ± 0.03.
Since the pressure dependence of Tc puts a stringent test on any theory of superconductivity, we repeated the calculations of phonons and electronic band structures discussed above for isotropic, uniaxial (along c-axis), and biaxial (in the ab-plane) pressures . We find that while w(E2g) increases with increasing pressure, the density of states at the Fermi energy decreases. The EP constant also shows significant changes with pressure. Inserting all these competitive effects into the McMillan formula yields the pressure dependence of Tc shown in Fig. 4.
For isotropic pressure Tc decreases with increasing pressure almost linearly at a rate of -1.0 K/GPa, in excellent agreement with the experimental value of -1.1 K/GPa. We also predict a cusp in the Tc-P curve around P 20 GPa. A similar cusp was recently observed experimentally. Our calculations indicate that Tc should increase first and then decrease with increasing c-compression, while it should decrease rapidly with ab-compression. Hence, when single crystal samples of MgB2 become available, measurements of the ab- and c-compression dependence of Tc should provide a critical test of our theory.
MgB2 may be the ultimate BCS s-wave superconductor, with parameters controlling Tc fully optimized to yield the highest possible Tc. However, even if Tc cannot be increased further, the low cost, light mass, easy fabrication, nearly isotropic high conductivity of MgB2, which also has a large critical current, will no doubt find many important technological applications in the near future.
 For full details of this work, see the website:
 T. Yildirim, O. Gulseren, J.F. Lynn, C.M. Brown, T.J. Udovic, Q. Huang, N. Rogado, . A. Regan, M.A. Hayward, J.S. Slusky, T. He, M.K. Haas, P. Khalifah, K. Inumaru, and R. J. Cava, Phys. Rev. Lett. 87, 37001 (2001).
 T. Yildirim and O. Gulseren, J. Phys.Chem. Solids (in press, 2001).
Last modified 19-March-2003 by website owner: NCNR (attn: John Copley)