Related papers: Symplectic Integrators in Corotating Coordinates
Hill's equations are an approximation that is useful in a number of areas of astrophysics including planetary rings and planetesimal disks. We derive a symplectic method for integrating Hill's equations based on a generalized leapfrog. This…
In this paper, explicit stable integrators based on symplectic and contact geometries are proposed for a non-autonomous ordinarily differential equation (ODE) found in improving convergence rate of Nesterov's accelerated gradient method.…
By combining a standard symmetric, symplectic integrator with a new step size controller, we provide an integration scheme that is symmetric, reversible and conserves the values of the constants of motion. This new scheme is appropriate for…
Time-reversible symplectic methods, which are precisely compatible with Liouville's phase-volume-conservation theorem, are often recommended for computational simulations of Hamiltonian mechanics. Lack of energy drift is an apparent…
We give a detailed study of the symplectic geometry of a family of integrable systems obtained by coupling two angular momenta in a non trivial way. These systems depend on a parameter t $\in$ [0, 1] and exhibit different behaviors…
Contact integrators are a family of geometric numerical schemes which guarantee the conservation of the contact structure. In this work we review the construction of both the variational and Hamiltonian versions of these methods. We…
A new approach is developed to integrate numerically the equations of motion for systems of interacting rigid polyatomic molecules. With the aid of a leapfrog framework, we directly involve principal angular velocities into the integration,…
Singular theories, characterised by the presence of degeneracies in their Lagrangian or Hamiltonian descriptions, require the systematic implementation of constraints in order to obtain well-defined dynamics. While the symplectic framework…
The two-dimensional n-body problem of classical mechanics is a non-integrable Hamiltonian system for n > 2. Traditional numerical integration algorithms, which are polynomials in the time step, typically lead to systematic drifts in the…
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…
This paper presents a geometric-variational approach to continuous and discrete mechanics and field theories. Using multisymplectic geometry, we show that the existence of the fundamental geometric structures as well as their preservation…
Modified Hamiltonians are used in the field of geometric numerical integration to show that symplectic schemes for Hamiltonian systems are accurate over long times. For nonlinear systems the series defining the modified Hamiltonian usually…
This article considers Hamiltonian mechanical systems with potential functions admitting jump discontinuities. The focus is on accurate and efficient numerical approximations of their solutions, which will be defined via the laws of…
Massless spinning correlators in cosmology are extremely complicated. In contrast, the scattering amplitudes of massless particles with spin are very simple. We propose that the reason for the unreasonable complexity of these correlators…
We present a collection of examples borrowed from celestial mechanics and projective dynamics. In these examples symplectic structures with singularities arise naturally from regularization transformations, Appell's transformation or…
By introducing an integration factor to the differential one-form of contact dynamics, equations of motion are derived variationally, and contact Poisson bracket and contact Lagrangian are formulated. Discrete symplectic integrator, named…
We introduce a recent symplectic integration scheme derived for solving physically motivated systems with non-separable Hamiltonians. We show its relevance to Riemannian manifold Hamiltonian Monte Carlo (RMHMC) and provide an alternative to…
The properties of a system of n = 3 coupled oscillators with linear terms in the velocities (magnetic terms) depending in two parameters are studied. We proved the existence of a bi-Hamiltonian structure arising from a non-symplectic…
We study the optimal design of numerical integrators for dissipative systems, for which there exists an underlying thermodynamic structure known as GENERIC (general equation for the nonequilibrium reversible-irreversible coupling). We…
Most numerical integration algorithms are not designed specifically for Hamiltonian systems and do not respect their characteristic properties, which include the preservation of phase space volume with time. This can lead to spurious…