Refine
Document Type
- ZIB-Report (14)
- Habilitation (1)
Language
- English (15)
Has Fulltext
- yes (15)
Is part of the Bibliography
- no (15)
Keywords
Institute
- ZIB Allgemein (15)
Symplectic methods, like the Verlet method, are a standard tool for the long term integration of Hamiltonian systems as they arise, for example, in molecular dynamics. One of the reasons for the popularity of symplectic methods is the conservation of energy over very long periods of time up to small fluctuations that scale with the order of the method. In this paper, we discuss a qualitative feature of Hamiltonian systems with separated time scales that is also preserved under symplectic discretization. Specifically, highly oscillatory degrees of freedom often lead to almost preserved quantities (adiabatic invariants). Using recent results from backward error analysis and normal form theory, we show that a symplectic method, like the Verlet method, preserves those adiabatic invariants. We also discuss step-size restrictions necessary to maintain adiabatic invariants in practical computations.
This article considers the design and implementation of variable-timestep methods for simulating holonomically constrained mechanical systems. Symplectic variable stepsizes are briefly discussed, we then consider time-reparameterization techniques employing a time-reversible (symmetric) integration method to solve the equations of motion. We give several numerical examples, including a simulation of an elastic (inextensible, unshearable) rod undergoing large deformations and collisions with the sides of a bounding box. Numerical experiments indicate that adaptive stepping can significantly smooth the numerical energy and improve the overall efficiency of the simulation.
The overall Hamiltonian structure of the Quantum-Classical Molecular Dynamics model makes - analogously to classical molecular dynamics - symplectic integration schemes the methods of choice for long-term simulations. This has already been demonstrated by the symplectic PICKABACK method. However, this method requires a relatively small step-size due to the high-frequency quantum modes. Therefore, following related ideas from classical molecular dynamics, we investigate symplectic multiple-time-stepping methods and indicate various possibilities to overcome the step-size limitation of PICKABACK.
Dynamical systems with two well-separated time-scales are investigated using normal form theory. Exponential estimates for the normal form truncation error are derived and applied to the numerical integration of differential equations (backward error analysis) and the reduction of highly oscillatory Hamiltonian systems (constrained dynamics and correcting potentials). The theoretical results are used to formulate new algorithms for the time integration of conservative Hamiltonian systems (projected multiple time stepping, soft constraints, rigid bodies, symplectic variable step-size methods).
A variable step-size, semi-explicit variant of the explicit Störmer-Verlet method has been proposed for the time-reversible integration of Newton's equations of motion by Huang & Leimkuhler. Here we propose a fully explicit version of this approach applicable to explicit and symmetric integration methods for general time-reversible differential equations. As applications, we discuss the variable step-size, time-reversible, and fully explicit integration of rigid body motion and reversible Nos\'e-Hoover dynamics.