A revised version of the quaternion approach for numerical integration of the equations of motion for rigid polyatomic molecules is proposed. The modified approach is based on a formulation of the quaternion dynamics with constraints. This allows one to resolve the rigidity problem rigorously using constraint forces. It is shown that the procedure for preservation of molecular rigidity can be realized particularly simply within the Verlet algorithm in velocity form. We demonstrate that the method presented leads to an improved numerical stability with respect to the usual quaternion rescaling scheme and it is roughly as good as the cumbersome atomic-constraint technique. © 1998 American Institute of Physics.
REFERENCES
1.
2.
3.
H. Goldstein, Classical Mechanics, 2nd ed. (Addison–Wesley, Reading, MA, 1967).
4.
J.
Barojas
, D.
Levesque
, and B.
Quentrec
, Phys. Rev. A
7
, 1092
(1973
).5.
D.
Levesque
, J. J.
Weis
, and G. N.
Patey
, Mol. Phys.
51
, 333
(1984
).6.
7.
K.
Singer
, A.
Taylor
, and J. V. L.
Singer
, Mol. Phys.
33
, 1757
(1977
).8.
9.
10.
11.
12.
13.
M. P. Allen and D. J. Tildesley, Computer Simulation of Liquids (Clarendon, Oxford, 1987).
14.
D. C. Rapaport, The Art of Molecular Dynamics Simulation (Cambridge University Press, Cambridge, 1995).
15.
J. P.
Ryckaert
, G.
Ciccotti
and H. J. C.
Berendsen
, J. Comput. Phys.
23
, 327
(1977
).16.
G.
Ciccotti
, J. P.
Ryckaert
and M.
Ferrario
, Mol. Phys.
47
, 1253
(1982
).17.
C. W. Gear, ANL Report No. ANL-7126, 1966.
18.
C. W. Gear, Numerical Initial Value Problems in Ordinary Differential Equations (Prentice–Hall, Englewood Cliffs, NJ, 1971).
19.
20.
R. W. Hockney and J. W. Eastwood, Computer Simulation Using Particles (McGraw–Hill, New York, 1981).
21.
W. C.
Swope
, H. C.
Andersen
, P. H.
Berens
, and K. R.
Wilson
, J. Chem. Phys.
76
, 637
(1982
).22.
23.
J. M. Haile, Molecular Dynamics Simulations: Elementary Methods (Wiley-Interscience, New York, 1992).
24.
H. Gould and J. Tobochnik, An Introduction to Computer Simulation Methods: Applications to Physical Systems (Addison–Wesley, Reading, MA, 1988) Part 1.
25.
26.
27.
This content is only available via PDF.
© 1998 American Institute of Physics.
1998
American Institute of Physics