Characterizing phonon dynamics using stochastic sampling

Predicting phonon relaxation time from molecular dynamics (MD) requires a long simulation time to compute the mode energy auto-correlation function. Here, we present an alternative approach to infer the phonon life-time from an approximate form of the energy auto-correlation function. The method requires as an input a set of sampled equilibrium configurations. A stochastic sampling method is used to generate the equilibrium configurations. We consider a truncated Taylor series expansion of the phonon energy auto-correlation function. The different terms in the truncated correlation function are obtained using the stochastic sampling approach. The expansion terms, thus, obtained are in good agreement with the corresponding values obtained using MD. We then use the approximate function to compute the phonon relaxation time. The relaxation time computed using this method is compared with that obtained from the exact correlation function. The two values are in agreement with each other.