Symplectic reversible integrators: Predictor-corrector methods

Glenn J. Martyna, Mark Tuckerman

Research output: Contribution to journalArticle

Abstract

A new fourth order predictor-corrector integration scheme is presented. The unique feature of the new algorithm and what distinguishes it from a Gear predictor-corrector is that the method is derived from the Trotter decomposition of a specially formulated evolution operator and as such, is both symplectic and reversible. In addition, the method retains the useful property of Gear methods that only one force evaluation per time step is required. The new integrator is tested on a harmonic plus quartic oscillator and the Henon-Heiles system. Comparisons are made to the second order velocity Verlet integrator, the true fourth order Yoshida/Suzuki schemes and fourth order Gear. In all cases, the new method works well, giving energy conservation and trajectories of much better quality than velocity Verlet and of comparable quality to the results of the true fourth order schemes for the same computational cost as velocity Verlet.

Original languageEnglish (US)
Pages (from-to)8071-8077
Number of pages7
JournalThe Journal of chemical physics
Volume102
Issue number20
StatePublished - 1995

Fingerprint

predictor-corrector methods
integrators
Gears
energy conservation
predictions
Energy conservation
Trajectories
oscillators
trajectories
Decomposition
costs
harmonics
decomposition
operators
evaluation
Costs

ASJC Scopus subject areas

  • Atomic and Molecular Physics, and Optics

Cite this

Symplectic reversible integrators : Predictor-corrector methods. / Martyna, Glenn J.; Tuckerman, Mark.

In: The Journal of chemical physics, Vol. 102, No. 20, 1995, p. 8071-8077.

Research output: Contribution to journalArticle

@article{eb4a2cbb8b1c4e489290064c8aa9dc25,
title = "Symplectic reversible integrators: Predictor-corrector methods",
abstract = "A new fourth order predictor-corrector integration scheme is presented. The unique feature of the new algorithm and what distinguishes it from a Gear predictor-corrector is that the method is derived from the Trotter decomposition of a specially formulated evolution operator and as such, is both symplectic and reversible. In addition, the method retains the useful property of Gear methods that only one force evaluation per time step is required. The new integrator is tested on a harmonic plus quartic oscillator and the Henon-Heiles system. Comparisons are made to the second order velocity Verlet integrator, the true fourth order Yoshida/Suzuki schemes and fourth order Gear. In all cases, the new method works well, giving energy conservation and trajectories of much better quality than velocity Verlet and of comparable quality to the results of the true fourth order schemes for the same computational cost as velocity Verlet.",
author = "Martyna, {Glenn J.} and Mark Tuckerman",
year = "1995",
language = "English (US)",
volume = "102",
pages = "8071--8077",
journal = "Journal of Chemical Physics",
issn = "0021-9606",
publisher = "American Institute of Physics Publising LLC",
number = "20",

}

TY - JOUR

T1 - Symplectic reversible integrators

T2 - Predictor-corrector methods

AU - Martyna, Glenn J.

AU - Tuckerman, Mark

PY - 1995

Y1 - 1995

N2 - A new fourth order predictor-corrector integration scheme is presented. The unique feature of the new algorithm and what distinguishes it from a Gear predictor-corrector is that the method is derived from the Trotter decomposition of a specially formulated evolution operator and as such, is both symplectic and reversible. In addition, the method retains the useful property of Gear methods that only one force evaluation per time step is required. The new integrator is tested on a harmonic plus quartic oscillator and the Henon-Heiles system. Comparisons are made to the second order velocity Verlet integrator, the true fourth order Yoshida/Suzuki schemes and fourth order Gear. In all cases, the new method works well, giving energy conservation and trajectories of much better quality than velocity Verlet and of comparable quality to the results of the true fourth order schemes for the same computational cost as velocity Verlet.

AB - A new fourth order predictor-corrector integration scheme is presented. The unique feature of the new algorithm and what distinguishes it from a Gear predictor-corrector is that the method is derived from the Trotter decomposition of a specially formulated evolution operator and as such, is both symplectic and reversible. In addition, the method retains the useful property of Gear methods that only one force evaluation per time step is required. The new integrator is tested on a harmonic plus quartic oscillator and the Henon-Heiles system. Comparisons are made to the second order velocity Verlet integrator, the true fourth order Yoshida/Suzuki schemes and fourth order Gear. In all cases, the new method works well, giving energy conservation and trajectories of much better quality than velocity Verlet and of comparable quality to the results of the true fourth order schemes for the same computational cost as velocity Verlet.

UR - http://www.scopus.com/inward/record.url?scp=36449009301&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=36449009301&partnerID=8YFLogxK

M3 - Article

AN - SCOPUS:36449009301

VL - 102

SP - 8071

EP - 8077

JO - Journal of Chemical Physics

JF - Journal of Chemical Physics

SN - 0021-9606

IS - 20

ER -