Physics – Computational Physics
Scientific paper
1999-01-17
Physics
Computational Physics
27 pages, 2 figures; submitted to J. Comput. Phys
Scientific paper
A new approach is developed to integrate numerically the equations of motion for systems of interacting rigid polyatomic molecules. With the aid of a leapfrog framework, we directly involve principal angular velocities into the integration, whereas orientational positions are expressed in terms of either principal axes or quaternions. As a result, the rigidness of molecules appears to be an integral of motion, despite the atom trajectories are evaluated approximately. The algorithm derived is free of any iterative procedures and it allows to perform both energy- and temperature-conserving simulations. The corresponding integrators are time reversible but the symplectic behaviour is only achieved in mean. Symplectic versions are also described. The provide the conservation of volume in phase space precisely at each time step and, moreover, lead to exact solutions for angular velocities in the inertial-motion regime. It is shown that the algorithm exhibits excellent stability properties and conserves the energy even somewhat better than the atomic-constraint technique.
No associations
LandOfFree
A new approach to the integration of rotational motion in systems with interacting rigid bodies does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.
If you have personal experience with A new approach to the integration of rotational motion in systems with interacting rigid bodies, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and A new approach to the integration of rotational motion in systems with interacting rigid bodies will most certainly appreciate the feedback.
Profile ID: LFWR-SCP-O-575215