Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

The Journal of Chemical Physics

Quantum mechanical correlation functions, maximum entropy analytic continuation, and ring polymer molecular dynamics

Scott Habershon, Bastiaan J. Braams, David E. Manolopoulos

Published Nov 7, 2007• 69 citations
Publisher

Abstract

The maximum entropy analytic continuation (MEAC) and ring polymer molecular dynamics (RPMD) methods provide complementary approaches to the calculation of real time quantum correlation functions. RPMD becomes exact in the high temperature limit, where the thermal time betavariant Planck's over 2pi tends to zero and the ring polymer collapses to a single classical bead. MEAC becomes most reliable at low temperatures, where betavariant Planck's over 2pi exceeds the correlation time of interest and the numerical imaginary time correlation function contains essentially all of the information that is needed to recover the real time dynamics. We show here that this situation can be exploited by combining the two methods to give an improved approximation that is better than either of its parts. In particular, the MEAC method provides an ideal way to impose exact moment (or sum rule) constraints on a prior RPMD spectrum. The resulting scheme is shown to provide a practical solution to the "nonlinear operator problem" of RPMD, and to give good agreement with recent exact results for the short-time velocity autocorrelation function of liquid parahydrogen. Moreover these improvements are obtained with little extra effort, because the imaginary time correlation function that is used in the MEAC procedure can be computed at the same time as the RPMD approximation to the real time correlation function. However, there are still some problems involving long-time dynamics for which the RPMD+MEAC combination is inadequate, as we illustrate with an example application to the collective density fluctuations in liquid orthodeuterium.

DOI
10.1063/1.2786451

Related papers

Can the ring polymer molecular dynamics method be interpreted as real time quantum dynamics?

Seogjoo Jang, Anton V. Sinitskiy, Gregory A. Voth

The ring polymer molecular dynamics (RPMD) method has gained popularity in recent years as a simple approximation for calculating real time quantum correlation functions in condensed media. However, the extent to which RPMD captures real dynamical quantum effects and why it fails under certain situations have not been clearly understood. Addressing this issue has been difficult in the absence of a genuine justification for the RPMD algorithm starting from the quantum Liouville equation. To this end, a new and exact path integral formalism for the calculation of real time quantum correlation functions is presented in this work, which can serve as a rigorous foundation for the analysis of the RPMD method as well as providing an alternative derivation of the well established centroid molecular dynamics method. The new formalism utilizes the cyclic symmetry of the imaginary time path integral in the most general sense and enables the expression of Kubo-transformed quantum time correlation functions as that of physical observables pre-averaged over the imaginary time path. Upon filtering with a centroid constraint function, the formulation results in the centroid dynamics formalism. Upon filtering with the position representation of the imaginary time path integral, we obtain an exact quantum dynamics formalism involving the same variables as the RPMD method. The analysis of the RPMD approximation based on this approach clarifies that an explicit quantum dynamical justification does not exist for the use of the ring polymer harmonic potential term (imaginary time kinetic energy) as implemented in the RPMD method. It is analyzed why this can cause substantial errors in nonlinear correlation functions of harmonic oscillators. Such errors can be significant for general correlation functions of anharmonic systems. We also demonstrate that the short time accuracy of the exact path integral limit of RPMD is of lower order than those for finite discretization of path. The present quantum dynamics formulation also serves as the basis for developing new quantum dynamical methods that utilize the cyclic nature of the imaginary time path integral.

The Journal of Chemical Physics•Apr 16, 2014•64 citations

Bosonic Condensed Phase Real-time Dynamics from Ring Polymer Molecular Dynamics

Yotam M. Y. Feldman, Sourav Karmakar, Sutirtha Paul +3 more

We present a method for approximating real-time correlation functions and quantum transport coefficients of bosonic condensed phases. The direct evaluation of quantum real-time correlation functions in the path integral formulation is impossible due to a severe dynamical sign problem. Evaluating imaginary-time correlation functions and inverting them using analytic continuation is famously ill-posed. Ring polymer molecular dynamics (RPMD) provides an alternative approach resulting in accurate approximations for real-time correlation functions of condensed phases but neglects exchange effects. In this work, we develop a bosonic RPMD method, which is exact in the harmonic, short time, and high temperature limits. A critical enabling observation is that the Kubo-transformed correlation function for bosons is real only when the correlated observables are symmetric under exchange. We benchmark the method on harmonic and anharmonic model systems and then use it to directly obtain the Lieb-Liniger gas density-density real-time correlation functions at various temperatures and momentum transfers. This work enables the direct simulation of finite-temperature, real-time dynamics of large bosonic condensed phases using RPMD for the first time.

arXiv:2607.25892•Jul 28, 2026•0 citations

Approximate quantum dynamics methods for time correlation functions

Kyle K. G. Smith

The dynamic structure factor of liquid para-hydrogen and ortho-deuterium in corresponding thermodynamic states, (T = 20.0 K, n = 21.24 nm⁻³) and (T = 23.0 K, n = 24.61 nm⁻³) respectively, has been computed by both the Feynman-Kleinert linearized path-integral (FK-LPI) and Ring-Polymer Molecular Dynamics (RPMD) methods and compared with Inelastic X-ray Scattering spectra. The combined use of computational and experimental methods enables a reduction in experimental uncertainties for the determination of the true sample spectrum. Furthermore, the refined experimental spectrum of para-hydrogen and ortho-deuterium is consistently reproduced by both FK-LPI and RPMD at momentum transfers lower than 12.8nm⁻¹. At larger momentum transfers the F K - LP I results agree with experiment much better for ortho-deuterium than for para-hydrogen. More specifically we found that for k ~ 20.0 nm⁻¹ para-hydrogen provides a test case for improved approximations to quantum dynamics. We meet this demand for an improved approximate quantum dynamics method by developing two classes of quasi-classical dynamics that are shown to conserve the initial quantum ensemble when used in conjunction with the Feynman-Kleinert approximation of the density operator. As shown, both classes of dynamics are able to recover the exact classical and high temperature limits of the quantum time correlation function, while a subset is able to recover the exact harmonic limit. A comparison of the approximate quantum time correlation functions obtained from both classes of dynamics are made with the exact results for the challenging model problems of the quartic and double-well potentials. It is found that this new Feynman-Kleinert Quasi-Classical Wigner (FK-QCW) method provides a great improvement over the Feynman-Kleinert implementation of the classical Wigner approximation, also known as FK-LPI, in which purely classical dynamics are used. Furthermore, it is shown that the first class of dynamics reduces to Centroid Molecular Dynamics (CMD) when used within the framework of the classical Wigner approximation for the Kubo transformed time correlation function. Finally, we apply the Feynman-Kleinert Quasi-Classical Wigner (FK- QCW) method to the same liquid para-hydrogen and ortho-deuterium system, previously studied using FK-LPI and RPMD. When applied to this challenging system, it is shown that this new FK-QCW method consistently reproduces the experimental dynamic structure factor for all momentum transfers considered. This shows that FK-QCW provides a great improvement over FK-LPI for not only model problems, but also realistic systems. Furthermore, for small momentum transfers, where RPMD is applicable, it is shown that FK-QCW provides nearly the same results as RPMD, thus suggesting that FK-QCW provides a potentially more appealing algorithm than RPMD since one is not limited to correlation functions involving linear operators. This then suggests that the FK-QCW method is a top contender in the realm of approximate quantum dynamics methods which allow for the practical evaluation of time correlation functions.

•May 1, 2014•0 citations

Real time quantum correlation functions. II. Maximum entropy numerical analytic continuation of path integral Monte Carlo and centroid molecular dynamics data

Goran Krilov, B. J. Berne

We propose a method which uses centroid molecular dynamics (CMD) [J. Cao and G. A. Voth, J. Chem. Phys. 100, 5106 (1994)] real-time data in conjunction with the imaginary-time data generated using path integral Monte Carlo simulations in a numerical analytic continuation scheme based on the maximum entropy approach. We show that significant improvement is achieved by including short-time CMD data with the imaginary-time data. In particular, for a particle bilinearly coupled to a harmonic bath, these methods lead to significant improvements over previous calculations and even allow accurate determination of transport coefficients such as the diffusion coefficient and mobility for this system. In addition we show how maximum entropy method can be used to extract accurate dynamic information from short-time CMD data, and that this approach is superior to the direct Fourier transform of long-time data for systems characterized by broad, featureless spectral distributions.

The Journal of Chemical Physics•Nov 22, 1999•54 citations

Quantum time correlation functions from complex time Monte Carlo simulations: A maximum entropy approach

Goran Krilov, Eunji Sim, B. J. Berne

We present a way of combining real-time path integral Monte Carlo simulations with a maximum entropy numerical analytic continuation scheme in a new approach for calculating time correlation functions for finite temperature many body quantum systems. The real-time dynamics is expressed in the form of the symmetrized time correlation function, which is suitable for Monte Carlo methods, and several simulation techniques are presented for evaluating this function accurately up to moderate values of time. The symmetrized time correlation function is then analytically continued in combination with imaginary time data to obtain the real-time correlation function. We test this approach on several exactly solvable problems, including two one-dimensional systems, as well two cases of vibrational relaxation of a system coupled to a dissipative environment. The computed time correlation functions are in good agreement with exact results over several multiples of the thermal time βℏ, and exhibit a significant improvement over analytic continuation of imaginary time correlation functions. Moreover, we show how the method can be systematically improved.

The Journal of Chemical Physics•Jan 15, 2001•77 citations

Effective potential analytic continuation approach for real time quantum correlation functions involving nonlinear operators

Horikoshi, A, Kinugawa, K

We apply the effective potential analytic continuation (EPAC) method to the calculation of real time quantum correlation functions involving operators nonlinear in the position operator $\\hat{q}$. For a harmonic system the EPAC method provides the exact correlation function at all temperature ranges, while the other quantum dynamics methods, the centroid molecular dynamics and the ring polymer molecular dynamics, become worse at lower temperature. For an asymmetric anharmonic system, the EPAC correlation function is in very good agreement with the exact one at $t=0$. When the time increases from zero, the EPAC method gives good coincidence with the exact result at lower temperature. Finally, we propose a simplified version of the EPAC method to reduce the computational cost required for the calculation of the standard effective potential.

•May 6, 2005•33 citations
The maximum entropy analytic continuation (MEAC) and ring polymer molecular dynamics (RPMD) methods provide complementary approaches to the calculation of real time quantum correlation functions. RPMD becomes exact in the high temperature limit, where the thermal time betavariant Planck's over 2pi tends to zero and the ring polymer collapses to a single classical bead. MEAC becomes most reliable at low temperatures, where betavariant Planck's over 2pi exceeds the correlation time of interest and the numerical imaginary time correlation function contains essentially all of the information that is needed to recover the real time dynamics. We show here that this situation can be exploited by combining the two methods to give an improved approximation that is better than either of its parts. In particular, the MEAC method provides an ideal way to impose exact moment (or sum rule) constraints on a prior RPMD spectrum. The resulting scheme is shown to provide a practical solution to the "nonlinear operator problem" of RPMD, and to give good agreement with recent exact results for the short-time velocity autocorrelation function of liquid parahydrogen. Moreover these improvements are obtained with little extra effort, because the imaginary time correlation function that is used in the MEAC procedure can be computed at the same time as the RPMD approximation to the real time correlation function. However, there are still some problems involving long-time dynamics for which the RPMD+MEAC combination is inadequate, as we illustrate with an example application to the collective density fluctuations in liquid orthodeuterium.
The ring polymer molecular dynamics (RPMD) method has gained popularity in recent years as a simple approximation for calculating real time quantum correlation functions in condensed media. However, the extent to which RPMD captures real dynamical quantum effects and why it fails under certain situations have not been clearly understood. Addressing this issue has been difficult in the absence of a genuine justification for the RPMD algorithm starting from the quantum Liouville equation. To this end, a new and exact path integral formalism for the calculation of real time quantum correlation functions is presented in this work, which can serve as a rigorous foundation for the analysis of the RPMD method as well as providing an alternative derivation of the well established centroid molecular dynamics method. The new formalism utilizes the cyclic symmetry of the imaginary time path integral in the most general sense and enables the expression of Kubo-transformed quantum time correlation functions as that of physical observables pre-averaged over the imaginary time path. Upon filtering with a centroid constraint function, the formulation results in the centroid dynamics formalism. Upon filtering with the position representation of the imaginary time path integral, we obtain an exact quantum dynamics formalism involving the same variables as the RPMD method. The analysis of the RPMD approximation based on this approach clarifies that an explicit quantum dynamical justification does not exist for the use of the ring polymer harmonic potential term (imaginary time kinetic energy) as implemented in the RPMD method. It is analyzed why this can cause substantial errors in nonlinear correlation functions of harmonic oscillators. Such errors can be significant for general correlation functions of anharmonic systems. We also demonstrate that the short time accuracy of the exact path integral limit of RPMD is of lower order than those for finite discretization of path. The present quantum dynamics formulation also serves as the basis for developing new quantum dynamical methods that utilize the cyclic nature of the imaginary time path integral.
We present a method for approximating real-time correlation functions and quantum transport coefficients of bosonic condensed phases. The direct evaluation of quantum real-time correlation functions in the path integral formulation is impossible due to a severe dynamical sign problem. Evaluating imaginary-time correlation functions and inverting them using analytic continuation is famously ill-posed. Ring polymer molecular dynamics (RPMD) provides an alternative approach resulting in accurate approximations for real-time correlation functions of condensed phases but neglects exchange effects. In this work, we develop a bosonic RPMD method, which is exact in the harmonic, short time, and high temperature limits. A critical enabling observation is that the Kubo-transformed correlation function for bosons is real only when the correlated observables are symmetric under exchange. We benchmark the method on harmonic and anharmonic model systems and then use it to directly obtain the Lieb-Liniger gas density-density real-time correlation functions at various temperatures and momentum transfers. This work enables the direct simulation of finite-temperature, real-time dynamics of large bosonic condensed phases using RPMD for the first time.
The dynamic structure factor of liquid para-hydrogen and ortho-deuterium in corresponding thermodynamic states, (T = 20.0 K, n = 21.24 nm⁻³) and (T = 23.0 K, n = 24.61 nm⁻³) respectively, has been computed by both the Feynman-Kleinert linearized path-integral (FK-LPI) and Ring-Polymer Molecular Dynamics (RPMD) methods and compared with Inelastic X-ray Scattering spectra. The combined use of computational and experimental methods enables a reduction in experimental uncertainties for the determination of the true sample spectrum. Furthermore, the refined experimental spectrum of para-hydrogen and ortho-deuterium is consistently reproduced by both FK-LPI and RPMD at momentum transfers lower than 12.8nm⁻¹. At larger momentum transfers the F K - LP I results agree with experiment much better for ortho-deuterium than for para-hydrogen. More specifically we found that for k ~ 20.0 nm⁻¹ para-hydrogen provides a test case for improved approximations to quantum dynamics. We meet this demand for an improved approximate quantum dynamics method by developing two classes of quasi-classical dynamics that are shown to conserve the initial quantum ensemble when used in conjunction with the Feynman-Kleinert approximation of the density operator. As shown, both classes of dynamics are able to recover the exact classical and high temperature limits of the quantum time correlation function, while a subset is able to recover the exact harmonic limit. A comparison of the approximate quantum time correlation functions obtained from both classes of dynamics are made with the exact results for the challenging model problems of the quartic and double-well potentials. It is found that this new Feynman-Kleinert Quasi-Classical Wigner (FK-QCW) method provides a great improvement over the Feynman-Kleinert implementation of the classical Wigner approximation, also known as FK-LPI, in which purely classical dynamics are used. Furthermore, it is shown that the first class of dynamics reduces to Centroid Molecular Dynamics (CMD) when used within the framework of the classical Wigner approximation for the Kubo transformed time correlation function. Finally, we apply the Feynman-Kleinert Quasi-Classical Wigner (FK- QCW) method to the same liquid para-hydrogen and ortho-deuterium system, previously studied using FK-LPI and RPMD. When applied to this challenging system, it is shown that this new FK-QCW method consistently reproduces the experimental dynamic structure factor for all momentum transfers considered. This shows that FK-QCW provides a great improvement over FK-LPI for not only model problems, but also realistic systems. Furthermore, for small momentum transfers, where RPMD is applicable, it is shown that FK-QCW provides nearly the same results as RPMD, thus suggesting that FK-QCW provides a potentially more appealing algorithm than RPMD since one is not limited to correlation functions involving linear operators. This then suggests that the FK-QCW method is a top contender in the realm of approximate quantum dynamics methods which allow for the practical evaluation of time correlation functions.
We propose a method which uses centroid molecular dynamics (CMD) [J. Cao and G. A. Voth, J. Chem. Phys. 100, 5106 (1994)] real-time data in conjunction with the imaginary-time data generated using path integral Monte Carlo simulations in a numerical analytic continuation scheme based on the maximum entropy approach. We show that significant improvement is achieved by including short-time CMD data with the imaginary-time data. In particular, for a particle bilinearly coupled to a harmonic bath, these methods lead to significant improvements over previous calculations and even allow accurate determination of transport coefficients such as the diffusion coefficient and mobility for this system. In addition we show how maximum entropy method can be used to extract accurate dynamic information from short-time CMD data, and that this approach is superior to the direct Fourier transform of long-time data for systems characterized by broad, featureless spectral distributions.
We present a way of combining real-time path integral Monte Carlo simulations with a maximum entropy numerical analytic continuation scheme in a new approach for calculating time correlation functions for finite temperature many body quantum systems. The real-time dynamics is expressed in the form of the symmetrized time correlation function, which is suitable for Monte Carlo methods, and several simulation techniques are presented for evaluating this function accurately up to moderate values of time. The symmetrized time correlation function is then analytically continued in combination with imaginary time data to obtain the real-time correlation function. We test this approach on several exactly solvable problems, including two one-dimensional systems, as well two cases of vibrational relaxation of a system coupled to a dissipative environment. The computed time correlation functions are in good agreement with exact results over several multiples of the thermal time βℏ, and exhibit a significant improvement over analytic continuation of imaginary time correlation functions. Moreover, we show how the method can be systematically improved.
We apply the effective potential analytic continuation (EPAC) method to the calculation of real time quantum correlation functions involving operators nonlinear in the position operator hatq\\hat{q}hatq. For a harmonic system the EPAC method provides the exact correlation function at all temperature ranges, while the other quantum dynamics methods, the centroid molecular dynamics and the ring polymer molecular dynamics, become worse at lower temperature. For an asymmetric anharmonic system, the EPAC correlation function is in very good agreement with the exact one at t=0t=0t=0. When the time increases from zero, the EPAC method gives good coincidence with the exact result at lower temperature. Finally, we propose a simplified version of the EPAC method to reduce the computational cost required for the calculation of the standard effective potential.