Related papers: Improving the accuracy of simulated chaotic $N$-bo…
Dependable numerical results from long-time simulations require stable numerical integration schemes. For Hamiltonian systems, this is achieved with symplectic integrators, which conserve the symplectic condition and exactly solve for the…
We study the Hamiltonian dynamics and spectral theory of spin-oscillators. Because of their rich structure, spin-oscillators display fairly general properties of integrable systems with two degrees of freedom. Spin-oscillators have…
Spin-orbit coupling of planetary systems plays an important role in the dynamics and habitability of planets. However, symplectic integrators that can accurately simulate not only how orbit affects spin but also how spin affects orbit have…
We examine the non-linear stability of the Wisdom-Holman (WH) mapping applied to the integration of perturbed, highly eccentric two-body orbits. We find that the method is unstable and introduces artificial chaos into the computed…
In recent publications, the construction of explicit symplectic integrators for Schwarzschild and Kerr type spacetimes is based on splitting and composition methods for numerical integrations of Hamiltonians or time-transformed Hamiltonians…
Two new families of completely integrable perturbations of the N-dimensional isotropic harmonic oscillator Hamiltonian are presented. Such perturbations depend on arbitrary functions and N free parameters and their integrals of motion are…
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…
There is strong evidence that cosmological N-body simulations dominated by Warm Dark Matter (WDM) contain spurious or unphysical haloes, most readily apparent as regularly spaced low-mass haloes strung along filaments. We show that spurious…
We present a new mixed variable symplectic (MVS) integrator for planetary systems, that fully resolve close encounters. The method is based on a time regularisation that allows keeping the stability properties of the symplectic integrators,…
We review the implementation of individual particle time-stepping for N-body dynamics. We present a class of integrators derived from second order Hamiltonian splitting. In contrast to the usual implementation of individual time-stepping,…
Quantum mechanical systems with some degree of complexity due to multiple scattering behave as if their Hamiltonians were random matrices. Such behavior, while originally surmised for the interacting many-body system of highly excited…
Hamiltonian systems are known to conserve the Hamiltonian function, which describes the energy evolution over time. Obtaining a numerical spatio-temporal scheme that accurately preserves the discretized Hamiltonian function is often a…
Multibody dynamics simulators are an important tool in many fields, including learning and control for robotics. However, many existing dynamics simulators suffer from inaccuracies when dealing with constrained mechanical systems due to…
Dynamical systems with translational or rotational symmetry arise frequently in studies of spatially extended physical systems, such as Navier-Stokes flows on periodic domains. In these cases, it is natural to express the state of the fluid…
We present a structure preserving discretization of the fundamental spacetime geometric structures of fluid mechanics in the Lagrangian description in 2D and 3D. Based on this, multisymplectic variational integrators are developed for…
Classical Hamiltonian spin systems are continuous dynamical systems on the symplectic phase space $(S^2)^n$. In this paper we investigate the underlying geometry of a time discretization scheme for classical Hamiltonian spin systems called…
Discrete control systems, as considered here, refer to the control theory of discrete-time Lagrangian or Hamiltonian systems. These discrete-time models are based on a discrete variational principle, and are part of the broader field of…
Numerical algorithms based on variational and symplectic integrators exhibit special features that make them promising candidates for application to general relativity and other constrained Hamiltonian systems. This paper lays part of the…
A trajectory segment in an energy shell, which combines to form a closed curve with a segment in another canonically driven energy shell, adds an oscillatory semiclassical contribution to the smooth classical background of the quantum…
We design a novel, exactly energy-conserving implicit non-symplectic integration method for an eight-dimensional Hamiltonian system with four degrees of freedom. In our algorithm, each partial derivative of the Hamiltonian with respect to…