Direct Dynamics Simulations using Hessian-based Predictor-corrector Integration Algorithms
Journal Article
·
· Journal of Chemical Physics
The research described in this product was performed in part in the Environmental Molecular Sciences Laboratory, a national scientific user facility sponsored by the Department of Energy's Office of Biological and Environmental Research and located at Pacific Northwest National Laboratory. In previous research (J. Chem. Phys. 111, 3800 (1999)) a Hessian-based integration algorithm was derived for performing direct dynamics simulations. In the work presented here, improvements to this algorithm are described. The algorithm has a predictor step based on a local second-order Taylor expansion of the potential in Cartesian coordinates, within a trust radius, and a fifth-order correction to this predicted trajectory. The current algorithm determines the predicted trajectory in Cartesian coordinates, instead of the instantaneous normal mode coordinates used previously, to ensure angular momentum conservation. For the previous algorithm the corrected step was evaluated in rotated Cartesian coordinates. Since the local potential expanded in Cartesian coordinates is not invariant to rotation, the constants of motion are not necessarily conserved during the corrector step. An approximate correction to this shortcoming was made by projecting translation and rotation out of the rotated coordinates. For the current algorithm unrotated Cartesian coordinates are used for the corrected step to assure the constants of motion are conserved. An algorithm is proposed for updating the trust radius to enhance the accuracy and efficiency of the numerical integration. This modified Hessian-based integration algorithm, with its new components, has been implemented into the VENUS/NWChem software package and compared with the velocity-Verlet algorithm for the H₂CO→H₂+CO, O₃+C₃H₆, and F-+CH₃OOH chemical reactions.
- Research Organization:
- Pacific Northwest National Laboratory (PNNL), Richland, WA (US), Environmental Molecular Sciences Laboratory (EMSL)
- Sponsoring Organization:
- USDOE
- DOE Contract Number:
- AC05-76RL01830
- OSTI ID:
- 921827
- Journal Information:
- Journal of Chemical Physics, Journal Name: Journal of Chemical Physics Journal Issue: 4 Vol. 126; ISSN 0021-9606
- Publisher:
- American Institute of Physics (AIP)
- Country of Publication:
- United States
- Language:
- English
Similar Records
An evaluation of the Gauss-Radau algorithm for the simulation of chemical dynamics
Solving Newton’s equations of motion with large timesteps using recurrent neural networks based operators
Journal Article
·
Mon Aug 01 00:00:00 EDT 1994
· Journal of Computational Physics; (United States)
·
OSTI ID:7032698
Solving Newton’s equations of motion with large timesteps using recurrent neural networks based operators
Journal Article
·
Mon Apr 04 20:00:00 EDT 2022
· Machine Learning: Science and Technology
·
OSTI ID:1861193