Related papers: Symplectic Integrators in Corotating Coordinates
Relativistic dynamics of a charged particle in time-dependent electromagnetic fields has theoretical significance and a wide range of applications. It is often multi-scale and requires accurate long-term numerical simulations using…
Two specialized algorithms for the numerical integration of the equations of motion of a Brownian walker obeying detailed balance are introduced. The algorithms become symplectic in the appropriate limits, and reproduce the equilibrium…
Many important physical systems can be described as the evolution of a Hamiltonian system, which has the important property of being conservative, that is, energy is conserved throughout the evolution. Physics Informed Neural Networks and…
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 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…
The simulation of systems that act on multiple time scales is challenging. A stable integration of the fast dynamics requires a highly accurate approximation whereas for the simulation of the slow part, a coarser approximation is accurate…
As a generalization and extension of our previous paper [Escobar-Ruiz and Azuaje, J. Phys. A: Math. Theor. 57, 105202 (2024)], in this work, the notions of particular integral and particular integrability in classical mechanics are extended…
Efficient fourth order symplectic integrators are proposed for numerical integration of separable Hamiltonian systems H(p,q)=T(p)+V(q). Symmetric splitting coefficients with five to nine stages are obtained by higher order decomposition of…
Effective Liouville operators of the first- and the second-order symplectic integrators are obtained for the one-dimensional harmonic-oscillator system. The operators are defined only when the time step is less than two. Absolute values of…
In general relativity, the motion of an extended body moving in a given spacetime can be described by a particle on a (generally non-geodesic) worldline. In first approximation, this worldline is a geodesic of the underlying spacetime, and…
The construction of modified equations is an important step in the backward error analysis of symplectic integrator for Hamiltonian systems. In the context of partial differential equations, the standard construction leads to modified…
We present a non-canonically symplectic integration scheme tailored to numerically computing the post-Newtonian motion of a spinning black-hole binary. Using a splitting approach we combine the flows of orbital and spin contributions. In…
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…
We present a discrete total variation calculus in Hamiltonian formalism in this paper. Using this discrete variation calculus and generating functions for the flows of Hamiltonian systems, we derive two-step symplectic-energy integrators of…
Studying single-particle dynamics over many periods of oscillations is a well-understood problem solved using symplectic integration. Such integration schemes derive their update sequence from an approximate Hamiltonian, guaranteeing that…
Symplectic integrators offer many advantages for the numerical solution of Hamiltonian differential equations, including bounded energy error and the preservation of invariant sets. Two of the central Hamiltonian systems encountered in…
We present a multiscale integrator for Hamiltonian systems with slowly varying quadratic stiff potentials that uses coarse timesteps (analogous to what the impulse method uses for constant quadratic stiff potentials). This method is based…
Direct N-body simulations and symplectic integrators are effective tools to study the long-term evolution of planetary systems. The Wisdom-Holman (WH) integrator in particular has been used extensively in planetary dynamics as it allows for…
In this paper, an implicit nonsymplectic exact energy-preserving integrator is specifically designed for a ten-dimensional phase-space conservative Hamiltonian system with five degrees of freedom. It is based on a suitable…
In this paper, we compare the performance of different numerical schemes in approximating Pontryagin's Maximum Principle's necessary conditions for the optimal control of nonholonomic systems. Retraction maps are used as a seed to construct…