Centroid methods dynamical

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 ... 

Further improvement of the centroid method came with the introduction of centroid dynamics.Here the fundamental idea is to construct a centroid Hamiltonian in the full phase space of the system and the bath. The Boltzman factor is then the one obtained from this centroid Hamiltonian while seal time dynamics is obtained by running classical trajectories. This method has been applied to realistic systems " and recently derived from first principles.244 The main advantage of the centroid methodology is that thermodynamic quantum effects can be computed numerically exactly as it is not too difficult to converge numerically the computation of the centroid potential. 

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. 

Voth, G.A. (1996). Path-integral centroid methods in quantum statistical mechanics and dynamics. Adv. Chem. Phys. 93, 135-218... 

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... 

Takayanagi T, Shiga M. (2003) Photodissociation of CI2 in helium clusters an application of hybrid method of quantum wavepacket d5uiamics and path integral centroid molecular dynamics. Chem. Phys. Lett. 372 90-96. 

Much work has been recently done in extending path integral and centroid molecular dynamics (CMD) methods to include nuclear quantum effects in classical MD simulations [230-233]. Tachikawa etal [230] studiedp-CH O) (n= 1-3) by ab initio hybrid Monte Carlo and ab initio path integral simulations. Their simulation showed that, due to quantum effects, the average hydrogen-bonded... 

Hwang et al.131 were the first to calculate the contribution of tunneling and other nuclear quantum effects to enzyme catalysis. Since then, and in particular in the past few years, there has been a significant increase in simulations of QM-nuclear effects in enzyme reactions. The approaches used range from the quantized classical path (QCP) (e.g., Refs. 4,57,136), the centroid path integral approach,137,138 and vibrational TS theory,139 to the molecular dynamics with quantum transition (MDQT) surface hopping method.140 Most studies did not yet examine the reference water reaction, and thus could only evaluate the QM contribution to the enzyme rate constant, rather than the corresponding catalytic effect. However, studies that explored the actual catalytic contributions (e.g., Refs. 4,57,136) concluded that the QM contributions are similar for the reaction in the enzyme and in solution, and thus, do not contribute to catalysis. 

Variations on this surface hopping method that utilize Pechukas [106] formulation of mixed quantum-classical dynamics have been proposed [107,108]. Surface hopping algorithms [109-111] for non-adiabatic dynamics based on the quantum-classical Liouville equation [109,111-113] have been formulated. In these schemes the dynamics is fully prescribed by the quantum-classical Liouville operator and no additional assumptions about the nature of the classical evolution or the quantum transition probabilities are made. Quantum dynamics of condensed phase systems has also been carried out using techniques that are not based on surface hopping algorithms, in particular, centroid path integral dynamics [114] and influence functional methods [115]. 

Applications are then presented in Section IV. These examples should served as a guide as to what kinds of problems can be studied with these techniques and the limitations and possibilities for these methods. We present three examples (1) a dynamical test of the centroid quantum transition-state theory for electron transfer (ET) reactions in the crossover regime between adiabatic and nonadiabatic electron transfer, (2) the primary electron transfer reaction in bacterial photosynthesis, and (3) the diffusion kinetics of a Brownian particle in a periodic potential. Finally, Section V offers an outlook and a perspective of the current status of the field from our vantage point. 

In the early papers [4,8], the development of the CMD method was guided in part by the effective harmonic analysis and, in part, by physical reasoning. In Paper III, however, a mathematical justification of CMD was provided. In the latter analysis, it was shown that (1) CMD always yields a mathematically well-defined approximation to the quantum Kubo-transformed position or velocity correlation function, and (2) the equilibrium path centroid variable occupies an important role in the time correlation function because of the nature of the preaveraging procedure in CMD. Critical to the analysis of CMD and its justification was the phase-space centroid density formulation of Paper III, so that the momentum could be treated as an independent dynamical variable. The relationship between the centroid correlation function and the Kubo-transformed position correlation function was found to be unique if the centroid is taken as a dynamical variable. The analysis of Paper III will now be reviewed. For notational simplicity, the equations are restricted to a two-dimensional phase space, but they can readily be generalized. 

---