Centroid molecular dynamics

Cao, J., Voth, G.A. The formulation of quantum statistical mechanics based on the Feynman path centroid density. I. Equilibrium properties. J. Chem. Phys. 100 (1994) 5093-5105 II Dynamical properties. J. Chem. Phys. 100 (1994) 5106-5117 III. Phase space formalism and nalysis of centroid molecular dynamics. J. Chem. Phys. 101 (1994) 6157-6167 IV. Algorithms for centroid molecular dynamics. J. Chem. Phys. 101 (1994) 6168-6183 V. Quantum instantaneous normal mode theory of liquids. J. Chem. Phys. 101 (1994) 6184 6192. 

The Feynman-Hibbs and QFH potentials have been used extensively in simulations examining quantum effects in atomic and molecular fluids [12,15,25]. We note here that the centroid molecular dynamics method [54, 55] is related and is intermediate between a full path integral simulation and the Feynman-Hibbs approximation ... 

Diffusion constants are enhanced with the approximate inclusion of quantum effects. Changes in the ratio of diffusion constants for water and D2O with decreasing temperature are accurately reproduced with the QFF1 model. This ratio computed with the QFF1 model agrees well with the centroid molecular dynamics result at room temperature. Fully quantum path integral dynamical simulations of diffusion in liquid water are not presently possible. 

Abstract The theoretical basis for the quantum time evolution of path integral centroid variables is described, as weU as the motivation for using these variables to study condensed phase quantum dynamics. The equihbrium centroid distribution is shown to be a well-defined distribution function in the canonical ensemble. A quantum mechanical quasi-density operator (QDO) can then be associated with each value of the distribution so that, upon the application of rigorous quantum mechanics, it can be used to provide an exact definition of both static and dynamical centroid variables. Various properties of the dynamical centroid variables can thus be defined and explored. Importantly, this perspective shows that the centroid constraint on the imaginary time paths introduces a non-stationarity in the equihbrium ensemble. This, in turn, can be proven to yield information on the correlations of spontaneous dynamical fluctuations. This exact formalism also leads to a derivation of Centroid Molecular Dynamics, as well as the basis for systematic improvements of that theory. 

Following Fey nman s original work, several authors pmsued extensions of the effective potential idea to construct variational approximations for the quantum partition function (see, e g., Refs. 7,8). The importance of the path centroid variable in quantum activated rate processes was also explored and revealed, which gave rise to path integral quantum transition state theory and even more general approaches. The Centroid Molecular Dynamics (CMD) method for quantum dynamics simulation was also formulated. In the CMD method, the position centroid evolves classically on the efiective centroid potential. Various analysis and numerical tests for realistic systems have shown that CMD captures the main quantum effects for several processes in condensed matter such as transport phenomena. 

D. Marx, M. E. Tuckerman, and G. J. Martyna (1999) Quantum dynamics via adiabatic ab initio centroid molecular dynamics. Comput. Phvs. Commun. 

A disadvantage of using Car-Parrinello path integral methods is that the molecular dynamics is used only to compute averaged properties, the simulation dynamics having no direct physical meaning. A recently developed, albeit approximate method for generating fully quantum mechanical dynamics is the ab initio centroid molecular dynamics method (Marx et al., 1999 Pavese et al., 1999). Tlie application of Car-Parrinello methods to... 

General Time Correlation Functions Centroid Molecular Dynamics Method... 

---