dorsal/arxiv
View SchemaQuasi-Hamiltonian Equations of Motion for Internal Coordinate Molecular Dynamics of Polymers
| Authors | Alexey K. Mazur |
|---|---|
| Categories | |
| ArXiv ID | physics/9703019 |
| URL | https://arxiv.org/abs/physics/9703019 |
| Journal | J.Comput.Chem.18,1354-1364,1997. |
Abstract
Conventional molecular dynamics simulations macromolecules require long computational times because the most interesting motions are very slow compared with the fast oscillations of bond lengths and bond angles that limit the integration time step. Simulation of dynamics in the space of internal coordinates, that is with bond lengths, bond angles and torsions as independent variables, gives a theoretical possibility to eliminate all uninteresting fast degrees of freedom from the system. This paper presents a new method for internal coordinate molecular dynamics simulations of macromolecules. Equations of motion are derived which are applicable to branched chain molecules with any number of internal degrees of freedom. Equations use the canonical variables and they are much simpler than existing analogs. In the numerical tests the internal coordinate dynamics are compared with the traditional Cartesian coordinate molecular dynamics in simulations a 56 residue globular protein. It is shown that the traditional and internal coordinate dynamics require the same time step size for the same accuracy and that in the standard geometry approximation of amino acids, that is with fixed bond lengths, bond angles and rigid aromatic groups, the characteristic step size is 4 fsec, that is two times higher than with fixed bond lengths only. The step size can be increased up to 11 fsec when rotation of hydrogen atoms is suppressed.
{
"annotation_id": "565c102c-1251-4e2d-9309-3b6d548c18f8",
"date_created": "2026-03-02T18:01:18.552000Z",
"date_modified": "2026-03-02T18:01:18.552000Z",
"file_hash": "3a37bf8670f8414111dd2d2af4160ab5d091d9cff3527044716ca7c12a8b5dbc",
"private": false,
"record": {
"abstract": "Conventional molecular dynamics simulations macromolecules require long\ncomputational times because the most interesting motions are very slow compared\nwith the fast oscillations of bond lengths and bond angles that limit the\nintegration time step. Simulation of dynamics in the space of internal\ncoordinates, that is with bond lengths, bond angles and torsions as independent\nvariables, gives a theoretical possibility to eliminate all uninteresting fast\ndegrees of freedom from the system. This paper presents a new method for\ninternal coordinate molecular dynamics simulations of macromolecules. Equations\nof motion are derived which are applicable to branched chain molecules with any\nnumber of internal degrees of freedom. Equations use the canonical variables\nand they are much simpler than existing analogs. In the numerical tests the\ninternal coordinate dynamics are compared with the traditional Cartesian\ncoordinate molecular dynamics in simulations a 56 residue globular protein. It\nis shown that the traditional and internal coordinate dynamics require the same\ntime step size for the same accuracy and that in the standard geometry\napproximation of amino acids, that is with fixed bond lengths, bond angles and\nrigid aromatic groups, the characteristic step size is 4 fsec, that is two\ntimes higher than with fixed bond lengths only. The step size can be increased\nup to 11 fsec when rotation of hydrogen atoms is suppressed.",
"arxiv_id": "physics/9703019",
"authors": [
"Alexey K. Mazur"
],
"categories": [
"physics.chem-ph",
"physics.bio-ph",
"physics.comp-ph"
],
"journal_ref": "J.Comput.Chem.18,1354-1364,1997.",
"title": "Quasi-Hamiltonian Equations of Motion for Internal Coordinate Molecular Dynamics of Polymers",
"url": "https://arxiv.org/abs/physics/9703019"
},
"schema_id": "dorsal/arxiv",
"source": {
"execution_id": "3ca73a17-166d-4aef-909b-1ec1bde8a77d",
"id": "arXiv Dataset IDs",
"type": "Model",
"variant": "snapshot-2026-03-01",
"version": "0.1.0"
},
"user_id": 1000002
}