Real-time dynamics with bead-Fourier path integrals. I. Bead-Fourier CMD

Abstract

Developing new methods for the accurate and efficient calculations of real-time quantum correlation functions is deemed one of the most challenging problems of modern condensed matter theory. Many popular methods, such as centroid molecular dynamics (CMD), make use of Feynman path integrals (PIs) to efficiently introduce nuclear quantum effects into classical dynamical simulations. Conventional CMD methods use the discretized form of the PI formalism to represent a quantum particle using a series of replicas, or “beads,” connected with harmonic springs to create an imaginary time ring polymer. The alternative Fourier PI methodology, instead, represents the imaginary time path using a Fourier sine series. Presented as an intermediary between the two formalisms, bead-Fourier PIs (BF-PIs) have been shown to reduce the number of beads needed to converge equilibrium properties by including a few terms of the Fourier series. Here, a new CMD method is presented where the effective potential is calculated using BF-PIs as opposed to the typical bead PIs. We demonstrate the accuracy and efficiency of this new BF-CMD method for a series of 1D model systems and show that at low temperatures, one can achieve between a fourfold and eightfold reduction in the number of beads with the addition of a single Fourier component. The developed methodology is general and can be extended to other closely related methods, such as ring polymer molecular dynamics, as well as non-adiabatic PI methods.