Abstract

The atomic motion controls important properties of materials, such as thermal transport, phase transitions, and vibrational spectra. However, simulating the ionic dynamics is exceptionally challenging when quantum fluctuations are relevant (e.g., at low temperatures or with light atoms) and the energy landscape is anharmonic. In this paper, we present the time-dependent self-consistent harmonic approximation (TDSCHA) [L. Monacelli and F. Mauri, Phys. Rev. B 103, 104305 (2021), 10.1103/PhysRevB.103.104305] in the Wigner framework, paving the way for the efficient computation of the nuclear motion in systems with sizable quantum and thermal anharmonic fluctuations. Besides the improved numerical efficiency, the Wigner formalism unveils the classical limit of TDSCHA and provides a link with the many-body perturbation theory of Feynman diagrams. We further extend the method to account for the nonlinear couplings between phonons and photons, responsible, e.g., for a nonvanishing Raman signal in high-symmetry Raman inactive crystals, firstly discussed by Rasetti and Fermi. We benchmark the method in phase III of high-pressure hydrogen ab initio. The nonlinear photon-phonon coupling reshapes the IR spectra and explains the high-frequency shoulder of the H₂ vibron observed in experiments. The Wigner TDSCHA is computationally cheap and derived from first principles. It is unbiased by assumptions on the phonon-phonon and phonon-photon scattering and does not depend on empirical parameters. Therefore, the method can be adopted in unsupervised high-throughput calculations.