Related papers: Improving the accuracy of simulated chaotic $N$-bo…
We review a recent approach for the simulation of many-body interacting systems based on an efficient generalization of the Lanczos method for Quantum Monte Carlo simulations. This technique allows to perform systematic corrections to a…
We discuss normal forms and symplectic invariants of parabolic orbits and cuspidal tori in integrable Hamiltonian systems with two degrees of freedom. Such singularities appear in many integrable systems in geometry and mathematical physics…
We consider two strongly hyperbolic Hamiltonian formulations of general relativity and their numerical integration with a free and a partially constrained symplectic integrator. In those formulations we use hyperbolic drivers for the shift…
This paper is a continuation of our study of the dynamics of contact Hamiltonian systems in \cite{JY}, but without monotonicity assumption. Due to the complexity of general cases, we focus on the behavior of action minimizing orbits. We…
Symplectic schemes are powerful methods for numerically integrating Hamiltonian systems, and their long-term accuracy and fidelity have been proved both theoretically and numerically. However direct applications of standard symplectic…
We construct a Nekhoroshev-like result of stability with sharp constants for the planar three body problem, both in the planetary and in the restricted circular case, by using the periodic averaging technique. Our constructions can be…
We derive, and discuss the properties of, a symplectic map for the dynamics of bodies on nearly parabolic orbits. The orbits are perturbed by a planet on a circular, coplanar orbit interior to the pericenter of the parabolic orbit. The map…
We discuss an instability in the leapfrog integration algorithm, widely used in current Hybrid Monte Carlo (HMC) simulations of lattice QCD. We demonstrate the instability in the simple harmonic oscillator (SHO) system where it is manifest.…
Integrals of the Liouville $1$-form, known as the first Poincar\'e integral invariant, provide a computable figure of merit for monitoring the conservation of symplecticity in the numerical integration of Hamiltonian systems. These…
The discovery of Pluto's small moons in the last decade brought attention to the dynamics of the dwarf planet's satellites. With such systems in mind, we study a planar $N$-body system in which all the bodies are point masses, except for a…
We present a class of non-standard numerical schemes which are modifications of the discrete gradient method. They preserve the energy integral exactly (up to the round-off error). The considered class contains locally exact discrete…
In this paper we try to find examples of integrable natural Hamiltonian systems on the sphere $S^2$ with the symmetries of each Platonic polyhedra. Although some of these systems are known, their expression is extremely complicated; we try…
A goal of this work is to study the dynamics in the complement of KAM tori with focus on non-local robust transitivity. We introduce $C^r$ open sets ($r=1, 2, ..., \infty$) of symplectic diffeomorphisms and Hamiltonian systems, exhibiting…
We introduce a class of fourth order symplectic algorithms that are ideal for doing long time integration of gravitational few-body problems. These algorithms have only positive time steps, but require computing the force gradient in…
This paper summarises an investigation of the effects of low amplitude noise and periodic driving on phase space transport in 3-D Hamiltonian systems, a problem directly applicable to systems like galaxies, where such perturbations reflect…
Preintegration is a technique for high-dimensional integration over $d$-dimensional Euclidean space, which is designed to reduce an integral whose integrand contains kinks or jumps to a $(d-1)$-dimensional integral of a smooth function. The…
Generating intelligent robot behavior in contact-rich settings is a research problem where zeroth-order methods currently prevail. A major contributor to the success of such methods is their robustness in the face of non-smooth and…
We present a set of new, efficient high-order symplectic methods designed for Hamiltonian systems with cubic or quartic potentials. By demonstrating that polynomial potentials require fewer order conditions, we develop schemes that…
Long-term stability studies of nonlinear Hamiltonian systems require symplectic integration algorithms which are both fast and accurate. In this paper, we study a symplectic integration method wherein the symplectic map representing the…
Many applications in computational physics that use numerical integrators based on splitting and composition can benefit from the development of optimized algorithms and from choosing the best ordering of terms. The cost in programming and…