Accession Number : ADA191663
Title : A Constraint Algorithm for Maintaining Rigid Bonds in Molecular Dynamics Simulations of Large Molecules.
Descriptive Note : Memorandum rept.,
Corporate Author : NAVAL RESEARCH LAB WASHINGTON DC
Personal Author(s) : Lambrakos, S G ; Boris, J P ; Oran, E S ; Chandrasekhar, I ; Nagumo, M
PDF Url : ADA191663
Report Date : 04 Mar 1988
Pagination or Media Count : 39
Abstract : We present an efficient new algorithm, the Multiple Constraint Force (MCF) algorithm, for maintaining fixed distances in molecular dynamics simulations of polyatomic molecules. A reversible leapfrog integration scheme is modified to calculate the constraint forces required to maintain fixed interparticle distances (usually corresponding to chemical bonds). The procedure is exact for diatomic systems. It converges in a small number of iterations for polyatomic molecules and scales linearly with the number of bonds to be kept rigid. The MCF algorithm, like the SHAKE algorithm, iteratively corrects the constrained equations of motion. Rather than successively correcting particle positions or inverting a matrix, however, the MCF algorithm iteratively correct the constraint forces required to maintain the fixed distances. Keywords: Constraint algorithm, Molecules dynamics, Polyatomic molecules.
Descriptors : *ALGORITHMS, *DYNAMICS, *EQUATIONS OF MOTION, *MOLECULES, BONDED JOINTS, CHEMICAL BONDS, DIATOMIC MOLECULES, EFFICIENCY, INTEGRATION, ITERATIONS, MOLECULAR PROPERTIES, PARTICLES, POLYATOMIC MOLECULES, RIGIDITY, SHAKING, SIMULATION
Subject Categories : Physical Chemistry
Distribution Statement : APPROVED FOR PUBLIC RELEASE