Real-time dynamics with bead-Fourier path integrals. II. Bead-Fourier RPMD

Feynman path integrals (PIs) have found many uses in approximate quantum dynamics methods that are able to efficiently calculate real-time quantum correlation functions. The PIs typically take the form of discrete imaginary time slices over a closed path, where the slices form the ``beads'' of a ring polymer (RP) necklace. Some methods, such as centroid molecular dynamics (CMD), use the RP to generate an effective potential for the dynamics, while others, like RP molecular dynamics (RPMD), directly utilize the RP in real-time dynamics in order to incorporate quantum effects. The standard, discretized bead forms of CMD and RPMD can require a large number of RP beads to provide accurate results for systems where quantum effects are strong, such as at low temperatures. In Paper I, we introduced the bead-Fourier (BF) CMD method, where we utilized the inclusion of a Fourier sine series to reduce the number of beads needed to converge the CMD effective potential up to eightfold. In this work, we extend RPMD to incorporate BF-PIs in the form of BF-RPMD. We study a number of different implementations of the method through the calculation of correlation functions for both linear and non-linear operators. The effectiveness of the BF-RPMD method is sensitive to both the system and form of the operators being studied, but we show that this method is able to produce results on par with standard RPMD, with at worst twofold and up to eightfold reduction in the number of beads by including two to three Fourier components.