Anharmonic stabilization of the high-pressure simple cubic phase of calcium

The phonon spectrum of the high-pressure simple cubic phase of calcium, in the harmonic approximation, shows imaginary branches that make it mechanically unstable. In this Letter, the phonon spectrum is recalculated by using density-functional theory ab initio methods fully including anharmonic effects up to fourth order at 50 GPa. Considering that the perturbation theory cannot be employed with imaginary harmonic frequencies, a variational procedure based on the Gibbs-Bogoliubov inequality is used to estimate the renormalized phonon frequencies. The results show that strong quantum anharmonic effects make the imaginary phonons become positive even at zero temperature so that the simple cubic phase becomes mechanically stable, as experiments suggest. Moreover, our calculations find a superconducting T(c) in agreement with experiments and predict an anomalous behavior of the specific heat.