Graphene physics and insulator-metal transition in compressed hydrogen

Magnesium diboride differs from ordinary metallic superconductors in several important ways, including the failure of conventional models to predict accurately its unusually high transition temperature, the effects of isotope substitution on the critical transition temperature, and its anomalous specific heat. A detailed examination of the energy associated with the formation of charge-carrying pairs, referred to as the 'superconducting energy gap', should clarify why MgB(2) is different. Some early experimental studies have indicated that MgB(2) has multiple gaps, but past theoretical studies have not explained from first principles the origin of these gaps and their effects. Here we report an ab initio calculation of the superconducting gaps in MgB(2) and their effects on measurable quantities. An important feature is that the electronic states dominated by orbitals in the boron plane couple strongly to specific phonon modes, making pair formation favourable. This explains the high transition temperature, the anomalous structure in the specific heat, and the existence of multiple gaps in this material. Our analysis suggests comparable or higher transition temperatures may result in layered materials based on B, C and N with partially filled planar orbitals.

Boron in MgB_2 forms layers of honeycomb lattices with magnesium as a space filler. Band structure calculations indicate that Mg is substantially ionized, and the bands at the Fermi level derive mainly from B orbitals. Strong bonding with an ionic component and considerable metallic density of states yield a sizeable electron-phonon coupling. Using the rigid atomic sphere approximation and an analogy to Al, we estimate the coupling constant lambda to be of order 1. Together with high phonon frequencies, which we estimate via zone-center frozen phonon calculations to be between 300 and 700 cm^-1, this produces a high critical temperature, consistent with recent experiments reporting Tc=39 K (J. Akimitsu et al., to be published). Thus MgB_2 can be viewed as an analog of the long sought, but still hypothetical, superconducting metallic hydrogen.

We present a new nonempirical density functional generalized gradient approximation (GGA) that gives significant improvements for lattice constants, crystal structures, and metal surface energies over the most popular Perdew-Burke-Ernzerhof (PBE) GGA. The new functional is based on a diffuse radial cutoff for the exchange-hole in real space, and the analytic gradient expansion of the exchange energy for small gradients. There are no adjustable parameters, the constraining conditions of PBE are maintained, and the functional is easily implemented in existing codes.