Centroid density method

Warshel and Chu [42] and Hwang et al. [60] were the first to calculate the contribution of tunneling and other nuclear quantum effects to PT in solution and enzyme catalysis, respectively. Since then, and in particular in the past few years, there has been a significant increase in simulations of quantum mechanical-nuclear effects in enzyme and in solution reactions [16]. The approaches used range from the quantized classical path (QCP) (for example. Refs. [4, 58, 95]), the centroid path integral approach [54, 55], and variational transition state theory [96], to the molecular dynamics with quantum transition (MDQT) surface hopping method [31] and density matrix evolution [97-99]. Most studies of enzymatic reactions did not yet examine the reference water reaction, and thus could only evaluate the quantum mechanical contribution to the enzyme rate constant, rather than the corresponding catalytic effect. However, studies that explored the actual catalytic contributions (for example. Refs. [4, 58, 95]) concluded that the quantum mechanical contributions are similar for the reaction in the enzyme and in solution, and thus, do not contribute to catalysis. 

It should be noted that the imaginary time correlation function in Eq. (2.4) provides a measure of the localization of quantum particles in condensed media [17-19,52]. From this point on, the notation denotes an averaging by integrating some centroid-dependent function over the centroid position q weighted by the normalized centroid density p q ) Z. An alternative method for defining the correlation function... 

It should be noted that these equations are to be solved for each position of the centroid q. The frequency in Eq. (2.27) is the same as the effective frequency obtained for the optimized LHO reference system using the path-integral centroid density version of the Gibbs-Bogoliubov variational method [1, pp. 303-307 2, pp. 86-96], Correspondingly, Eqs. (2.27) and (2.28) are exactly the same as those in the quadratic effective potential theory [1,21-23], The derivation above does not make use of the variational principle but, instead, is the result of the vertex renormalization procedure. The diagrammatic analysis thus provides a method of systematic identification and evaluation of the corrections to the variational theory [3],... 

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. 

Linear scaling can be achieved when localized molecular orbitals (LMOs) are used to create sparsity within the Slater matrix. The density of a LMO is confined to limited region of space around its centroid, so only a few LMOs need to be evaluated for each electron. The second consideration on the way to linear scaling is to accelerate the transformation from basis function to LMO s. Linear scaling QMC methods have been a popular field of research and several algorithms have already been published. 

---