The canonical equilibrium of constrained molecular models
arXiv:1105.0374 · doi:10.1140/epjst/e2011-01517-9
Abstract
In order to increase the efficiency of the computer simulation of biological molecules, it is very common to impose holonomic constraints on the fastest degrees of freedom; normally bond lengths, but also possibly bond angles. However, as any other element that affects the physical model, the imposition of constraints must be assessed from the point of view of accuracy: both the dynamics and the equilibrium statistical mechanics are model-dependent, and they will be changed if constraints are used. In this review, we investigate the accuracy of constrained models at the level of the equilibrium statistical mechanics distributions produced by the different dynamics. We carefully derive the canonical equilibrium distributions of both the constrained and unconstrained dynamics, comparing the two of them by means of a "stiff" approximation to the latter. We do so both in the case of flexible and hard constraints, i.e., when the value of the constrained coordinates depends on the conformation and when it is a constant number. We obtain the different correcting terms associated with the kinetic energy mass-metric tensor determinants, but also with the details of the potential energy in the vicinity of the constrained subspace (encoded in its first and second derivatives). This allows us to directly compare, at the conformational level, how the imposition of constraints changes the thermal equilibrium of molecular systems with respect to the unconstrained case. We also provide an extensive review of the relevant literature, and we show that all models previously reported can be considered special cases of the most general treatments presented in this work. Finally, we numerically analyze a simple methanol molecule in order to illustrate the theoretical concepts in a practical case.
arXiv comment51 pages, 7 figures, Springer LaTeX style
References in corpus (13)
- Time-reversible Born-Oppenheimer molecular dynamics
- A modified Ehrenfest formalism for efficient large-scale ab initio molecular dynamics
- Efficient formalism for large scale ab initio molecular dynamics based on time-dependent density functional theory
- A mathematical and computational review of Hartree-Fock SCF methods in Quantum Chemistry
- Introduction to protein folding for physicists
- Quantum mechanical calculation of the effects of stiff and rigid constraints in the conformational equilibrium of the Alanine dipeptide
- Efficient model chemistries for peptides. I. Split-valence Gaussian basis sets and the heterolevel approximation in RHF and MP2
- Definition of Systematic, Approximately Separable and Modular Internal Coordinates (SASMIC) for macromolecular simulation
- Constraint methods for determining pathways and free energy of activated processes
- Exact and efficient calculation of Lagrange multipliers in constrained biological polymers: Proteins and nucleic acids as example cases
- Explicit factorization of external coordinates in constrained Statistical Mechanics models
- An exact expression to calculate the derivatives of position-dependent observables in molecular simulations with flexible constraints
- A Quantum-mechanical Approach for Constrained Macromolecular Chains