Stochastic path-integral approach for predicting the superconducting temperatures of anharmonic solids

Abstract

We develop a stochastic path-integral approach for predicting the superconducting transition temperatures of anharmonic solids. By defining generalized Bloch basis, we generalize the formalism of the stochastic path-integral approach, which is originally developed for liquid systems. We implement the formalism for ab initio calculations using the projector augmented-wave method, and apply the implementation to estimate the superconducting transition temperatures of metallic deuterium and hydrogen sulfide. For metallic deuterium, which is approximately harmonic, our result coincides well with that obtained from the standard approach based on the harmonic approximation and the density functional perturbation theory. For hydrogen sulfide, we find that anharmonicity strongly suppresses the predicted superconducting transition temperature. Compared to the self-consistent harmonic approximation approach, our approach yields a transition temperature closer to the experimentally observed one.