Numerical integration of the cartesian equations of motion of a system with constraints: molecular dynamics of n-alkanes

https://doi.org/10.1016/0021-9991(77)90098-5Get rights and content

Abstract

A numerical algorithm integrating the 3N Cartesian equations of motion of a system of N points subject to holonomic constraints is formulated. The relations of constraint remain perfectly fulfilled at each step of the trajectory despite the approximate character of numerical integration. The method is applied to a molecular dynamics simulation of a liquid of 64 n-butane molecules and compared to a simulation using generalized coordinates. The method should be useful for molecular dynamics calculations on large molecules with internal degrees of freedom.

References (6)

  • J. Orban et al.

    Report on CECAM Workshop

    Methods in Molecular Dynamics

    (1974)
  • J. Barojas et al.

    Phys. Rev. A

    (1973)
  • A. Rahman et al.

    J. Chem. Phys.

    (1971)
    A. Rahman et al.

    J. Chem. Phys.

    (1974)
There are more references available in the full text version of this article.

Cited by (18431)

  • Avian cryptochrome 4 binds superoxide

    2024, Computational and Structural Biotechnology Journal
View all citing articles on Scopus

Faculté des Sciences, Université Libre de Bruxelles, Brussels, Belgium.

Gruppo Nazionale di Struttura della Materia, Consiglio Nazionale delle Ricerche, and Istituto di Fisica “G. Marconi,” Università di Roma, Roma, Italy.

View full text