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

TL;DR

This paper introduces a novel bead-Fourier path integral approach integrated with centroid molecular dynamics to improve the accuracy and efficiency of real-time quantum correlation function calculations, especially at low temperatures.

Contribution

The paper develops a new CMD method using bead-Fourier path integrals, reducing the number of beads needed for convergence and enhancing computational efficiency.

Findings

Achieves 4-fold reduction in beads at low temperatures

Demonstrates accuracy on 1D model systems

Method is extendable to other PI-based techniques

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 discretized 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 a 4-fold 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 (RPMD), as well as non-adiabatic PI methods.