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

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.

[1]  K. Kinugawa,et al.  Effective potential analytic continuation approach for real time quantum correlation functions involving nonlinear operators. , 2005, The Journal of chemical physics.

[2]  Emilio Gallicchio,et al.  On the application of numerical analytic continuation methods to the study of quantum mechanical vibrational relaxation processes , 1998 .

[3]  R. Kubo.  Statistical Physics II: Nonequilibrium Statistical Mechanics , 2003 .

[4]  Thomas F. Miller,et al.  Quantum diffusion in liquid water from ring polymer molecular dynamics. , 2005, The Journal of chemical physics.

[5]  D. Chandler,et al.  Introduction To Modern Statistical Mechanics , 1987 .

[6]  B. Berne,et al.  Quantum time correlation functions from complex time Monte Carlo simulations: A maximum entropy approach , 2001 .

[7]  E. Rabani,et al.  A Short-Time Quantum Mechanical Expansion Approach to Vibrational Relaxation † , 2001 .

[8]  Gregory A Voth,et al.  A comparative study of imaginary time path integral based methods for quantum dynamics. , 2006, The Journal of chemical physics.

[9]  N. Makri,et al.  Forward-backward semiclassical dynamics for quantum fluids using pair propagators: Application to liquid para-hydrogen , 2003 .

[10]  W. T. Grandy.  Maximum entropy in action: Buck, Brian and Macaulay, Vincent A., 1991, 220 pp., Clarendon Press, Oxford, £30 pb, ISBN 0-19-8539630 , 1995 .

[11]  I. R. Mcdonald,et al.  Theory of simple liquids , 1998 .

[12]  Y. Yonetani,et al.  Transport properties of liquid para-hydrogen: The path integral centroid molecular dynamics approach , 2003 .

[13]  Mark E. Tuckerman,et al.  Molecular dynamics algorithms for path integrals at constant pressure , 1999 .

[14]  M. Parrinello,et al.  Study of an F center in molten KCl , 1984 .

[15]  John Skilling,et al.  Maximum Entropy and Bayesian Methods , 1989 .

[16]  N. Makri,et al.  Symmetrized correlation function for liquid para-hydrogen using complex-time pair-product propagators. , 2006, The Journal of chemical physics.

[17]  Thomas F. Miller,et al.  Sum rule constraints on Kubo-transformed correlation functions , 2006 .

[18]  R. Bryan,et al.  Maximum entropy analysis of oversampled data problems , 1990, European Biophysics Journal.

[19]  D. Ceperley.  Path integrals in the theory of condensed helium , 1995 .

[20]  Ian R. Craig,et al.  A refined ring polymer molecular dynamics theory of chemical reaction rates. , 2005, The Journal of chemical physics.

[21]  B. Berne,et al.  Quantum mechanical canonical rate theory: A new approach based on the reactive flux and numerical analytic continuation methods , 2000 .

[22]  David E Manolopoulos,et al.  On the short-time limit of ring polymer molecular dynamics. , 2006, The Journal of chemical physics.

[23]  William H. Press,et al.  Numerical Recipes: FORTRAN , 1988 .

[24]  B. Berne,et al.  On path integral Monte Carlo simulations , 1982 .

[25]  Tsunenobu Yamamoto,et al.  Quantum Statistical Mechanical Theory of the Rate of Exchange Chemical Reactions in the Gas Phase , 1960 .

[26]  G. Voth,et al.  A centroid molecular dynamics study of liquid para-hydrogen and ortho-deuterium. , 2004, The Journal of chemical physics.

[27]  Emilio Gallicchio,et al.  The absorption spectrum of the solvated electron in fluid helium by maximum entropy inversion of imaginary time correlation functions from path integral Monte Carlo simulations , 1994 .

[28]  Ian R. Craig,et al.  Quantum statistics and classical mechanics: real time correlation functions from ring polymer molecular dynamics. , 2004, The Journal of chemical physics.

[29]  Jarrell,et al.  Quantum Monte Carlo simulations and maximum entropy: Dynamics from imaginary-time data. , 1991, Physical review. B, Condensed matter.

[30]  Thomas F. Miller,et al.  Quantum diffusion in liquid para-hydrogen from ring-polymer molecular dynamics. , 2005, The Journal of chemical physics.

[31]  R. Kubo.  Statistical-Mechanical Theory of Irreversible Processes : I. General Theory and Simple Applications to Magnetic and Conduction Problems , 1957 .

[32]  Peter G. Wolynes,et al.  Exploiting the isomorphism between quantum theory and classical statistical mechanics of polyatomic fluids , 1981 .

[33]  F. Bermejo,et al.  Microscopic dynamics in liquid deuterium: A transition from collective to single-particle regimes , 1997 .

[34]  K. S. Singwi,et al.  Theory of Slow Neutron Scattering by Liquids. I , 1962 .

[35]  William H. Miller,et al.  Quantum mechanical rate constants for bimolecular reactions , 1983 .

[36]  A. Singh,et al.  Theory for the reorientational dynamics in glass-forming liquids , 1997 .

[37]  B. Berne,et al.  The calculation of transport properties in quantum liquids using the maximum entropy numerical analytic continuation method: Application to liquid para-hydrogen , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[38]  Kurt Kremer,et al.  Molecular dynamics simulation of a polymer chain in solution , 1993 .

[39]  B. Berne,et al.  On the Bayesian approach to calculating time correlation functions in quantum systems; reaction dynamics and spectroscopy , 2001 .

[40]  B. Berne,et al.  Quantum Rate Constants from Short-Time Dynamics: An Analytic Continuation Approach† , 2001 .

[41]  Emilio Gallicchio,et al.  On the calculation of dynamical properties of solvated electrons by maximum entropy analytic continuation of path integral Monte Carlo data , 1996 .

[42]  Victor V. Goldman,et al.  The isotropic intermolecular potential for H2 and D2 in the solid and gas phases , 1978 .

[43]  H. C. Andersen.  Molecular dynamics simulations at constant pressure and/or temperature , 1980 .

[44]  Gerhard Hummer,et al.  System-Size Dependence of Diffusion Coefficients and Viscosities from Molecular Dynamics Simulations with Periodic Boundary Conditions , 2004 .

[45]  B. Berne,et al.  Real time quantum correlation functions. II. Maximum entropy numerical analytic continuation of path integral Monte Carlo and centroid molecular dynamics data , 1999 .

[46]  Ian R. Craig,et al.  Inelastic neutron scattering from liquid para-hydrogen by ring polymer molecular dynamics , 2006 .

[47]  Ian R. Craig,et al.  Chemical reaction rates from ring polymer molecular dynamics. , 2005, The Journal of chemical physics.

[48]  The simulation of electronic absorption spectrum of a chromophore coupled to a condensed phase environment: Maximum entropy versus singular value decomposition approaches , 1997 .

[49]  P. Rossky,et al.  Quantum diffusion in liquid para-hydrogen: An application of the Feynman-Kleinert linearized path integral approximation , 2004 .

[50]  Eran Rabani,et al.  A fully self-consistent treatment of collective fluctuations in quantum liquids. , 2004, The Journal of chemical physics.

[51]  L. Hove.  Correlations in Space and Time and Born Approximation Scattering in Systems of Interacting Particles , 1954 .

[52]  E. Rabani,et al.  A self-consistent mode-coupling theory for dynamical correlations in quantum liquids: Application to liquid para-hydrogen , 2002 .