Related papers: Explicit symplectic methods in black hole spacetim…
We propose to use the properties of the Lie algebra of the angular momentum to build symplectic integrators dedicated to the Hamiltonian of the free rigid body. By introducing a dependence of the coefficients of integrators on the moments…
In this manuscript, we propose efficient stochastic semi-explicit symplectic schemes tailored for nonseparable stochastic Hamiltonian systems (SHSs). These semi-explicit symplectic schemes are constructed by introducing augmented…
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…
Using a Newtonian model of the Solar System with all 8 planets, we perform extensive tests on various symplectic integrators of high orders, searching for the best splitting scheme for long term studies in the Solar System. These…
The paper intends to lay out the first steps towards constructing a unified framework to understand the symplectic and spectral theory of finite dimensional integrable Hamiltonian systems. While it is difficult to know what the best…
Our universe contains a great number of extremely compact and massive objects which are generally accepted to be black holes. Precise observations of orbital motion near candidate black holes have the potential to determine if they have the…
An approach is presented to resolve key paradoxes in black hole physics through the application of complex Riemannian spacetime. We extend the Schwarzschild metric into the complex domain, employing contour integration techniques to remove…
We construct a symplectic, globally defined, minimal-coordinate, equivariant integrator on products of 2-spheres. Examples of corresponding Hamiltonian systems, called spin systems, include the reduced free rigid body, the motion of point…
Symplectic integrators evolve dynamical systems according to modified Hamiltonians whose error terms are also well-defined Hamiltonians. The error of the algorithm is the sum of each error Hamiltonian's perturbation on the exact solution.…
We present several methods, which utilize symplectic integration techniques based on two and three part operator splitting, for numerically solving the equations of motion of the disordered, discrete nonlinear Schr\"odinger (DDNLS)…
In [8], some exact splittings are proposed for inhomogeneous quadratic differential equations including, for example, transport equations, kinetic equations, and Schr{\"o}dinger type equations with a rotation term. In this work, these exact…
For a large class of time-dependent non-Hermitain Hamiltonians expressed in terms linear and bilinear combinations of the generators for an Euclidean Lie-algebra respecting different types of PT-symmetries, we find explicit solutions to the…
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…
Generalized Additive Runge-Kutta schemes have shown to be a suitable tool for solving ordinary differential equations with additively partitioned right-hand sides. This work develops symplectic GARK schemes for additively partitioned…
Hamilton's equations of motion form a fundamental framework in various branches of physics, including astronomy, quantum mechanics, particle physics, and climate science. Classical numerical solvers are typically employed to compute the…
Most numerical integration algorithms are not designed specifically for Hamiltonian systems and do not respect their characteristic properties, which include the preservation of phase space volume with time. This can lead to spurious…
Calculating the long term solution of ordinary differential equations, such as those of the $N$-body problem, is central to understanding a wide range of dynamics in astrophysics, from galaxy formation to planetary chaos. Because generally…
This work outlines a fast, high-precision time-domain solver for scalar, electromagnetic and gravitational perturbations on hyperboloidal foliations of Kerr space-times. Time-domain Teukolsky equation solvers have typically used explicit…
We consider the numerical integration of the matrix Hill's equation. Parametric resonances can appear and this property is of great interest in many different physical applications. Usually, the Hill's equations originate from a Hamiltonian…
We propose a family of reliable symplectic integrators adapted to the Discrete Non-Linear Schr\"odinger equation; based on an idea of Yoshida (H. Yoshida, Construction of higher order symplectic integrators, Physics Letters A, 150, 5,6,7,…