Explicit symplectic integrators using Hessian-vector products
- Univ. de Valladolid (Spain). Dept. de Matematica Aplicada y Computacion
- Univ. of Illinois, Urbana, IL (United States). Dept. of Computer Science and Beckman Inst.
In 1991 Rowlands proposed an effectively fourth-order, effectively two-stage, explicit symplectic integrator based on using a Hessian-vector product to modify the force evaluation in the leapfrog method, and evidence indicates that for modest accuracy this method is highly competitive. Here the authors explore the possible existence of even more efficient fourth-order explicit symplectic integrators, also based on the use of Hessian-vector products and the concept of effective order. First it is shown that the cost of a force evaluation plus a Hessian-vector product is less than twice the cost of the force alone for a sum of two-body interactions. Then a new method is found that is generally better than both the method of Rowlands and that of Calvo, according to both a theoretical measure of the error and limited numerical experiments. The basic motivation behind the new method is quite simple: do a Hessian-vector computation only every other step, significantly cutting costs while only marginally increasing the error. The idea of effective order means that they allow for both the possibility of preprocessing the initial values before application of the basic method and the possibility of postprocessing the values obtained by the basic method, but at output points only. For some applications processing is unnecessary, but in any case processing has been shown to be possible at low additional cost. The derivation of the new method illustrates how to simplify by the use of H-series the determination of parameters for methods of increased effective accuracy.
- DOE Contract Number:
- FG02-91ER25099
- OSTI ID:
- 449537
- Journal Information:
- SIAM Journal on Scientific Computing, Vol. 18, Issue 1; Other Information: PBD: Jan 1997
- Country of Publication:
- United States
- Language:
- English
Similar Records
A family of symplectic integrators: Stability, accuracy, and molecular dynamics applications
A hybrid symplectic PIC/spectral scheme for one-dimensional electrostatic simulations