Numerical integration of the equations of motion for rigid polyatomics: The matrix method

✍ Scribed by Igor P. Omelyan

• Publisher: Elsevier Science
• Year: 1998
• Tongue: English
• Weight: 878 KB
• Volume: 109
• Category: Article
• ISSN: 0010-4655
• DOI: 10.1016/s0010-4655(98)00024-1

✦ Synopsis

A new scheme for numerical integration of motion for classical systems composed of rigid polyatomic molecules is proposed. The scheme is based on a matrix representation of the rotational degrees of freedom. The equations of motion are integrated within the Verlet framework in velocity form. it is shown that. contrary to previous methods, in the approach introduced the rigidity of molecules can he conserved automatically without any additional transformations. A comparison of various techniques with respect to numerical stability is made.