Nuclear Quantum Effects in Solids Using a Colored-Noise Thermostat

We present a method, based on a non-Markovian Langevin equation, to include quantum corrections to the classical dynamics of ions in a quasiharmonic system. By properly fitting the correlation function of the noise, one can vary the fluctuations in positions and momenta as a function of the vibrational frequency, and fit them so as to reproduce the quantum-mechanical behavior, with minimal a priori knowledge of the details of the system. We discuss the application of the thermostat to diamond and to ice $\mathrm{I}h$. We find that results in agreement with path-integral methods can be obtained using only a fraction of the computational effort.