Abstract
A new class of predictor–corrector integration methods intended for molecular dynamics simulation is presented. The methods are derived from the original Gear methods by analysing the numerical solution of the harmonic oscillator. The corrector coefficients are chosen to improve the time-reversibility while keeping maximum stability. Tests performed on systems with forces not dependent on velocities (classical two-body problem and three simple atomistic systems) show good energy conservation and reduced deviations of measured quantities in comparison with the Verlet method. As regards the systems with velocity-dependent right-hand sides, a moderate improvement is found for the Nosé–Hoover thermostat but no improvement for fixed bonds constrained by the general constraint dynamics based on Lagrange multipliers.
Sign-in/Register to access full text options 
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a callDisclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.
