The formulation of quantum statistical mechanics based on the Feynman path centroid density. IV. Algorithms for centroid molecular dynamics

Numerical algorithms are developed for the centroid molecular dynamics (centroid MD) method to calculate dynamical time correlation functions for general many‐body quantum systems. Approaches based on the normal mode path integral molecular dynamics and staging path integral Monte Carlo methods are described to carry out a direct calculation of the force on the centroid variables in the centroid MD algorithm. A more efficient, but approximate, scheme to compute the centroid force is devised which is based on the locally optimized harmonic approximation for the centroid potential. The centroid MD equations in the latter method can be solved with the help of an iterative procedure or through extended Lagrangian dynamics. A third algorithm introduces an effective centroid pseudopotential to approximate the full many‐body centroid mean force potential by effective pairwise centroid interactions. Numerical simulations for both prototype models and more realistic many‐body systems are performed to explore the f...

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

[2]  M. Toda,et al.  In: Statistical physics II , 1985 .

[3]  B. Berne,et al.  Many‐body dispersion forces of polarizable clusters and liquids , 1992 .

[4]  Tognetti,et al.  Effective potential for quantum correlation functions. , 1991, Physical review. A, Atomic, molecular, and optical physics.

[5]  Car,et al.  Unified approach for molecular dynamics and density-functional theory. , 1985, Physical review letters.

[6]  B. Berne,et al.  Low‐temperature variational approximation for the Feynman quantum propagator and its application to the simulation of quantum systems , 1990 .

[7]  B. Berne,et al.  Efficient molecular dynamics and hybrid Monte Carlo algorithms for path integrals , 1993 .

[8]  R. Pecora.  Dynamic Light Scattering , 1985 .

[9]  Jianshu Cao,et al.  The formulation of quantum statistical mechanics based on the Feynman path centroid density. III. Phase space formalism and analysis of centroid molecular dynamics , 1994 .

[10]  C. DeWitt-Morette,et al.  Techniques and Applications of Path Integration , 1981 .

[11]  D. Chandler,et al.  Quantum theory of polarization in liquids: Exact solution of the mean spherical and related approximations , 1982 .

[12]  B. Berne,et al.  A Born–Oppenheimer approximation for path integrals with an application to electron solvation in polarizable fluids , 1993 .

[13]  D. Remler,et al.  Molecular dynamics without effective potentials via the Car-Parrinello approach , 1990 .

[14]  Jianshu Cao,et al.  The formulation of quantum statistical mechanics based on the Feynman path centroid density. II. Dynamical properties , 1994 .

[15]  Bruce J. Berne,et al.  Path integral Monte Carlo studies of the behavior of excess electrons in simple fluids , 1987 .

[16]  Mark S. Swanson,et al.  Path integrals and quantum processes , 1992 .

[17]  Voth.  Calculation of equilibrium averages with Feynman-Hibbs effective classical potentials and similar variational approximations. , 1991, Physical review. A, Atomic, molecular, and optical physics.

[18]  R. Feynman,et al.  Effective classical partition functions. , 1986, Physical review. A, General physics.

[19]  D. Chandler,et al.  Excess electrons in simple fluids. I. General equilibrium theory for classical hard sphere solvents , 1984 .

[20]  Giachetti,et al.  Variational approach to quantum statistical mechanics of nonlinear systems with application to sine-Gordon chains. , 1985, Physical review letters.

[21]  Klein,et al.  Staging: A sampling technique for the Monte Carlo evaluation of path integrals. , 1985, Physical review. B, Condensed matter.

[22]  Bruce J. Berne,et al.  On the Simulation of Quantum Systems: Path Integral Methods , 1986 .

[23]  D. Ceperley,et al.  Simulation of quantum many-body systems by path-integral methods , 1984 .

[24]  Mark E. Tuckerman,et al.  Reversible multiple time scale molecular dynamics , 1992 .

[25]  R. W. Hall,et al.  A Path Integral Monte Carlo Study of Liquid Neon and the Quantum Effective Pair Potential , 1984 .

[26]  R. Feynman,et al.  Quantum Mechanics and Path Integrals , 1965 .

[27]  J. Doll,et al.  Partial averaging approach to Fourier coefficient path integration , 1986 .

[28]  J. Doll,et al.  Equilibrium and dynamical Fourier path integral methods , 2007 .