Bead-Fourier path integral molecular dynamics.

Molecular dynamics formulation of Bead-Fourier path integral method for simulation of quantum systems at finite temperatures is presented. Within this scheme, both the bead coordinates and Fourier coefficients, defining the path representing the quantum particle, are treated as generalized coordinates with corresponding generalized momenta and masses. Introduction of the Fourier harmonics together with the center-of-mass thermostating scheme is shown to remove the ergodicity problem, known to pose serious difficulties in standard path integral molecular dynamics simulations. The method is tested for quantum harmonic oscillator and hydrogen atom (Coulombic potential). The simulation results are compared with the exact analytical solutions available for both these systems. Convergence of the results with respect to the number of beads and Fourier harmonics is analyzed. It was shown that addition of a few Fourier harmonics already improves the simulation results substantially, even for a relatively small number of beads. The proposed Bead-Fourier path integral molecular dynamics is a reliable and efficient alternative to simulations of quantum systems.

[1]  S. Okazaki,et al.  Path integral molecular dynamics for Bose–Einstein and Fermi–Dirac statistics , 2000 .

[2]  Michele Parrinello,et al.  Efficient and general algorithms for path integral Car–Parrinello molecular dynamics , 1996 .

[3]  Hoover,et al.  Canonical dynamics: Equilibrium phase-space distributions. , 1985, Physical review. A, General physics.

[4]  J. A. Barker.  A quantum‐statistical Monte Carlo method; path integrals with boundary conditions , 1979 .

[5]  P. Rossky,et al.  Extracting rates of vibrational energy relaxation from centroid molecular dynamics , 2001 .

[6]  H. Metiu,et al.  The properties of electrons in sodalite saturated with alkali atoms , 1994 .

[7]  Qiang Shi,et al.  Centroid-based methods for calculating quantum reaction rate constants: Centroid sampling versus centroid dynamics , 2002 .

[8]  Bruce J. Berne,et al.  Nonergodicity in path integral molecular dynamics , 1984 .

[9]  M. Klein,et al.  Computer simulation of muonium in water , 1984 .

[10]  M. Klein,et al.  Nosé-Hoover chains : the canonical ensemble via continuous dynamics , 1992 .

[11]  D. Truhlar,et al.  A new Fourier path integral method, a more general scheme for extrapolation, and comparison of eight path integral methods for the quantum mechanical calculation of free energies , 2001 .

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

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

[14]  Jianshu Cao,et al.  A new perspective on quantum time correlation functions , 1993 .

[15]  Jianshu Cao,et al.  The formulation of quantum statistical mechanics based on the Feynman path centroid density. IV. Algorithms for centroid molecular dynamics , 1994 .

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

[17]  J. Doll,et al.  A Monte Carlo method for quantum Boltzmann statistical mechanics , 1984 .

[18]  J. D. Doll,et al.  Comment on “A comparison of the efficiency of Fourier- and discrete time-path integral Monte Carlo” [J. Chem. Phys. 109, 2123 (1998)] , 1999 .

[19]  Gregory A. Voth,et al.  The computer simulation of proton transport in water , 1999 .

[20]  M. C. Gordillo,et al.  A comparison of the efficiency of Fourier- and discrete time-path integral Monte Carlo , 1998 .

[21]  R. Coalson.  On the connection between Fourier coefficient and Discretized Cartesian path integration , 1986 .

[22]  Lyubartsev,et al.  Path-integral Monte Carlo method in quantum statistics for a system of N identical fermions. , 1993, Physical review. A, Atomic, molecular, and optical physics.

[23]  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 .

[24]  J. Doll,et al.  Cumulant methods and short time propagators , 1989 .

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

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