Related papers: Improving the accuracy of simulated chaotic $N$-bo…
In recent decades, there have been many attempts to construct symplectic integrators with variable time steps, with rather disappointing results. In this paper we identify the causes for this lack of performance, and find that they fall…
This paper summarises a numerical investigation of phase mixing in time-independent Hamiltonian systems that admit a coexistence of regular and chaotic phase space regions, allowing also for low amplitude perturbations idealised as periodic…
We derive the dynamics of several rigid bodies of arbitrary shape in a 2-dimensional inviscid and incompressible fluid, whose vorticity field is given by point vortices. We adopt the idea of Vankerschaver et al. (2009) to derive the…
Variational integrators are derived for structure-preserving simulation of stochastic Hamiltonian systems with a certain type of multiplicative noise arising in geometric mechanics. The derivation is based on a stochastic discrete…
Much of standard galaxy dynamics rests on the implicit assumption that the corresponding N-body problem is (near) integrable. This notion although leading to great simplification is by no means a fact. It is therefore important to develop…
An outstanding property of any Hamiltonian system is the symplecticity of its flow, namely, the continuous trajectory preserves volume in phase space. Given a symplectic but discrete trajectory generated by a transition matrix applied at a…
Multi-symplectic integrators are typically regarded as a discretization of the Hamiltonian partial differential equations. This is due to the fact that, for generic finite-dimensional Hamiltonian systems, there exists only one independent…
The formation and evolution of protoplanetary systems, the breeding grounds of planet formation, is a complex dynamical problem that involves many orders of magnitudes. To serve this purpose, we present a new hybrid algorithm that combines…
We describe an algorithm for long-term planetary orbit integrations, including the dominant post-Newtonian effects, that employs individual timesteps for each planet. The algorithm is symplectic and exhibits short-term errors that are…
We have proposed new algorithms for the numerical integration of the equations of motion for classical spin systems. In close analogy to symplectic integrators for Hamiltonian equations of motion used in Molecular Dynamics these algorithms…
We show that every parabolic orbit of a two-degree of freedom integrable system admits a $C^\infty$-smooth Hamiltonian circle action, which is persistent under small integrable $C^\infty$ perturbations. We deduce from this result the…
We present WHFast, a fast and accurate implementation of a Wisdom-Holman symplectic integrator for long-term orbit integrations of planetary systems. WHFast is significantly faster and conserves energy better than all other Wisdom-Holman…
In this article, a brief description of Discrete Mechanics and Variational Integrators which preserve the symplectic structure of the flow will be provided and a Newton-Raphson algorithm that can be used to solve implicit equations on the…
We present a new time-stepping criterion for N-body simulations that is based on the true dynamical time of a particle. This allows us to follow the orbits of particles correctly in all environments since it has better adaptivity than…
This paper reports a theory of Koopman operators for a class of hybrid dynamical systems with globally asymptotically stable periodic orbits, called hybrid limit-cycling systems. We leverage smooth structures intrinsic to the hybrid…
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…
We describe the symplectic structure and Hamiltonian dynamics for a class of Grassmannian manifolds. Using the two dimensional sphere ($S^2$) and disc ($D^2$) as illustrative cases, we write their path integral representations using…
The solution of problems in physics is often facilitated by a change of variables. In this work we present neural transformations to learn symmetries of Hamiltonian mechanical systems. Maintaining the Hamiltonian structure requires novel…
Symplectic integrators constructed from Hamiltonian and Lie formalisms are obtained as symplectic maps whose flow follows the exact solution of a "sourrounded" Hamiltonian K = H + h^k H_1. Those modified Hamiltonians depends virtually on…
Nonholonomic systems are variational models commonly used for mechanical systems with ideal no-slip constraints. This note provides a differential-geometric derivation of the nonholonomic equations of motion for an arbitrary rigid body…