Related papers: Numerical integration of variational equations
A new algorithm for numerical integration of the rigid-body equations of motion is proposed. The algorithm uses the leapfrog scheme and the quantities involved are angular velocities and orientational variables which can be expressed in…
Symplectic integrators can be excellent for Hamiltonian initial value problems. Reasons for this include their preservation of invariant sets like tori, good energy behaviour, nonexistence of attractors, and good behaviour of statistical…
When numerically integrating canonical Hamiltonian systems, the long-term conservation of some of its invariants, among which the Hamiltonian function itself, assumes a central role. The classical approach to this problem has led to the…
We present a Hamiltonian Monte Carlo algorithm to sample from multivariate Gaussian distributions in which the target space is constrained by linear and quadratic inequalities or products thereof. The Hamiltonian equations of motion can be…
Generic Hamiltonian systems have a mixed phase space, where classically disjoint regions of regular and chaotic motion coexist. We present an iterative method to construct an integrable approximation, which resembles the regular dynamics of…
This letter studies symmetric and symplectic exponential integrators when applied to numerically computing nonlinear Hamiltonian systems. We first establish the symmetry and symplecticity conditions of exponential integrators and then show…
We present new splitting methods designed for the numerical integration of near-integrable Hamiltonian systems, and in particular for planetary N-body problems, when one is interested in very accurate results over a large time span. We…
A complete geometric classification of symmetries of autonomous Hamiltonian mechanical systems is established; explaining how to obtain their associated conserved quantities in all cases. In particular, first we review well-known results…
Multisymplectic variational integrators are structure preserving numerical schemes especially designed for PDEs derived from covariant spacetime Hamilton principles. The goal of this paper is to study the properties of the temporal and…
Predicting the behaviors of Hamiltonian systems has been drawing increasing attention in scientific machine learning. However, the vast majority of the literature was focused on predicting separable Hamiltonian systems with their kinematic…
We describe a new algorithm for the integration of self-gravitating fluid systems using SPH method. We split the Hamiltonian of a self-gravitating fluid system to the gravitational potential and others (kinetic and internal energies) and…
This paper is a review of results which have been recently obtained by applying mathematical concepts drawn, in particular, from differential geometry and topology, to the physics of Hamiltonian dynamical systems with many degrees of…
The existence of explicit symplectic integrators for general nonseparable Hamiltonian systems is an open and important problem in both numerical analysis and computing in science and engineering, as explicit integrators are usually more…
In this work, we present the hitherto most efficient and accurate method for the numerical integration of post-Newtonian equations of motion. We first transform the Poisson system as given by the post-Newtonian approximation to canonically…
In this paper, we develop a structure-preserving discretization of the Lagrangian framework for electromagnetism, combining techniques from variational integrators and discrete differential forms. This leads to a general family of…
In this paper, we investigate the asymptotic error distributions of symplectic methods for stochastic Hamiltonian systems and further provide Hamiltonian-specific analysis that clarifies the superiority of symplectic methods. Our…
An extended Hamiltonian approach to conduct isothermal-isobaric molecular dynamics simulations with full cell flexibility is presented. The components of the metric tensor are used as the fictitious degrees of freedom for the cell, thus…
It is well known that symplectic Runge-Kutta and Partitioned Runge-Kutta methods exactly preserve {\em quadratic} first integrals (invariants of motion) of the system being integrated. While this property is often seen as a mere curiosity…
We construct numerical integrators for Hamiltonian problems that may advantageously replace the standard Verlet time-stepper within Hybrid Monte Carlo and related simulations. Past attempts have often aimed at boosting the order of accuracy…
We describe a technique for constructing a symplectic transfer map for a charged particle moving through an accelerator component with arbitrary three-dimensional static magnetic field. The transfer map is constructed by symplectic…