z-logo
open-access-imgOpen Access
On the numerical integration of motion for rigid polyatomics: The modified quaternion approach
Author(s) -
I. P. Omelyan
Publication year - 1998
Publication title -
computers in physics
Language(s) - English
Resource type - Journals
eISSN - 1558-4208
pISSN - 0894-1866
DOI - 10.1063/1.168642
Subject(s) - verlet integration , quaternion , rigidity (electromagnetism) , constraint (computer aided design) , classical mechanics , mathematics , numerical integration , molecular dynamics , physics , mathematical analysis , geometry , quantum mechanics
A revised version of the quaternion approach for numerical integration of theequations of motion for rigid polyatomic molecules is proposed. The modifiedapproach is based on a formulation of the quaternion dynamics with constraints.This allows to resolve the rigidity problem rigorously using constraint forces.It is shown that the procedure for preservation of molecular rigidity can berealized particularly simply within the Verlet algorithm in velocity form. Wedemonstrate that the presented method leads to an improved numerical stabilitywith respect to the usual quaternion rescaling scheme and it is roughly as goodas the cumbersome atomic-constraint technique.Comment: 14 pages, 2 figure

The content you want is available to Zendy users.

Already have an account? Click here to sign in.
Having issues? You can contact us here
Accelerating Research

Address

John Eccles House
Robert Robinson Avenue,
Oxford Science Park, Oxford
OX4 4GP, United Kingdom