Related papers: A study of symplectic integrators for planetary sy…
Hamiltonian systems of ordinary and partial differential equations are fundamental mathematical models spanning virtually all physical scales. A critical property for the robustness and stability of computational methods in such systems is…
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…
Due to the chaotic nature of the Solar System, the question of its long-term stability can only be answered in a statistical sense, for instance, based on numerical ensemble integrations of nearby orbits. Destabilization of the inner…
We present a comparison of different numerical techniques for the integration of variational equations. The methods presented can be applied to any autonomous Hamiltonian system whose kinetic energy is quadratic in the generalized momenta,…
In a previous paper, second- and fourth-order explicit symplectic integrators were designed for a Hamiltonian of the Schwarzschild black hole. Following this work, we continue to trace the possibility of the construction of explicit…
In this paper, we investigate the asymptotic error distributions of symplectic methods for stochastic Hamiltonian systems and further provide Hamiltonian-specific analysis that clarifies the superiority of symplectic methods. Our…
An analytical method for investigation of the evolution of dynamical systems {\it with independent on time accuracy} is developed for perturbed Hamiltonian systems. The error-free estimation using of computer algebra enables the application…
This paper deals with the dynamics of time-reversible Hamiltonian vector fields with 2 and 3 degrees of freedom around an elliptic equilibrium point in presence of symplectic involutions. The main results discuss the existence of…
We develop a rigorous analytical Hamiltonian formalism adapted to the study of the motion of two planets in co-orbital resonance. By constructing a complex domain of holomorphy for the planetary Hamilto-nian, we estimate the size of the…
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 propose a general scheme for separation of variables in the integrable Hamiltonian systems on orbits of the loop algebra $\mathfrak{sl}(2,\Complex)\times \mathcal{P}(\lambda,\lambda^{-1})$. In particular, we illustrate the scheme by…
We give a possible splitting method to a Hamiltonian for the description of charged particles moving around the Reissner-Nordstrom-(anti)-de Sitter black hole with an external magnetic field. This Hamiltonian can be separated into six…
Modified Hamiltonian Monte Carlo (MHMC) methods combine the ideas behind two popular sampling approaches: Hamiltonian Monte Carlo (HMC) and importance sampling. As in the HMC case, the bulk of the computational cost of MHMC algorithms lies…
A new method is proposed for integrating the equations of motion of an elastic filament. In the standard finite-difference and finite-element formulations the continuum equations of motion are discretized in space and time, but it is then…
Due to the chaotic nature of the Solar System, the question of its dynamic long-term stability can only be answered in a statistical sense, e.g. based on numerical ensemble integrations of nearby orbits. Destabilization, including…
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…
Geometric aspects play an important role in the construction and analysis of structure-preserving numerical methods for a wide variety of ordinary and partial differential equations. Here we review the development and theory of symplectic…
We prove that the recently developed semiexplicit symplectic integrators for non-separable Hamiltonian systems preserve any linear and quadratic invariants possessed by the Hamiltonian systems. This is in addition to being symmetric and…
Wisdom--Holman (WH) methods are algorithms used as a basis for a wide range of codes used to solve problems in solar system and planetary dynamics. The problems range from the growth and migration of planets to the stability of the solar…
It is well known that symplectic integrators lose their near energy preservation properties when variable step sizes are used. The most common approach to combine adaptive step sizes and symplectic integrators involves the Poincar\'e…