Explicit factorization of external coordinates in constrained Statistical Mechanics models
arXiv:q-bio/0512033 · doi:10.1002/jcc.20499
Abstract
If a macromolecule is described by curvilinear coordinates or rigid constraints are imposed, the equilibrium probability density that must be sampled in Monte Carlo simulations includes the determinants of different mass-metric tensors. In this work, we explicitly write the determinant of the mass-metric tensor G and of the reduced mass-metric tensor g, for any molecule, general internal coordinates and arbitrary constraints, as a product of two functions; one depending only on the external coordinates that describe the overall translation and rotation of the system, and the other only on the internal coordinates. This work extends previous results in the literature, proving with full generality that one may integrate out the external coordinates and perform Monte Carlo simulations in the internal conformational space of macromolecules. In addition, we give a general mathematical argument showing that the factorization is a consequence of the symmetries of the metric tensors involved. Finally, the determinant of the mass-metric tensor G is computed explicitly in a set of curvilinear coordinates specially well-suited for general branched molecules.
22 pages, 2 figures, LaTeX, AMSTeX. v2: Introduccion slightly extended. Version in arXiv is slightly larger than the published one
References in corpus (1)
Cited by in corpus (6)
- A mathematical and computational review of Hartree-Fock SCF methods in Quantum Chemistry
- Introduction to protein folding for physicists
- Efficient model chemistries for peptides. I. Split-valence Gaussian basis sets and the heterolevel approximation in RHF and MP2
- The canonical equilibrium of constrained molecular models
- Exact and efficient calculation of Lagrange multipliers in constrained biological polymers: Proteins and nucleic acids as example cases
- An exact expression to calculate the derivatives of position-dependent observables in molecular simulations with flexible constraints