Published in

Society for Industrial and Applied Mathematics, SIAM Journal on Scientific Computing, 4(23), p. 1172-1188

DOI: 10.1137/s106482750138318x

Links

Tools

Export citation

Search in Google Scholar

Practical Construction of Modified Hamiltonians

Journal article published in 2001 by Robert D. Skeel, David J. Hardy ORCID
This paper is available in a repository.
This paper is available in a repository.

Full text: Download

Green circle
Preprint: archiving allowed
Green circle
Postprint: archiving allowed
Green circle
Published version: archiving allowed
Data provided by SHERPA/RoMEO

Abstract

One of the most fruitful ways to analyze the effects of discretization error in the numerical solution of a system of differential equations is to examine the "modified equations," which are the equations that are exactly satisfied by the (approximate) discrete solution. These do not actually exist in general, but rather are defined by an asymptotic expansion in powers of the discretization parameter. Nonetheless, if the expansion is suitably truncated, the resulting modified equations have a solution which is remarkably close to the discrete solution. In the case of a Hamiltonian system of ordinary differential equations, the modified equations are also Hamiltonian if and only if the integrator is symplectic. The existence of a modified Hamiltonian is an indicator of the validity of statistical estimates calculated from long time integration of chaotic Hamiltonian systems. Evidence for the existence of a Hamiltonian for a particular calculation is obtained by calculating modified Hami...