相关论文: Anomaly in Symplectic Integrator
Symplectic integration algorithms have become popular in recent years in long-term orbital integrations because these algorithms enforce certain conservation laws that are intrinsic to Hamiltonian systems. For problems with large variations…
We show that when time-reversible symplectic algorithms are used to solve periodic motions, the energy error after one period is generally two orders higher than that of the algorithm. By use of correctable algorithms, we show that the…
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…
Explicit symplectic integrators have been important tools for accurate and efficient approximations of mechanical systems with separable Hamiltonians. For the first time, the article proposes for arbitrary Hamiltonians similar integrators,…
We present a class of symplectic integrators adapted for the integration of perturbed Hamiltonian systems of the form $H=A+\epsilon B$. We give a constructive proof that for all integer $p$, there exists an integrator with positive steps…
We study the time evolution of an ideal system composed of two harmonic oscillators coupled through a quadratic Hamiltonian with arbitrary interaction strength. We solve its dynamics analytically by employing tools from symplectic geometry.…
The evolution of any factorized time-reversible symplectic integrators, when applied to the harmonic oscillator, can be exactly solved in a closed form. The resulting modified Hamiltonians demonstrate the convergence of the Lie series…
A new family of methods involving complex coefficients for the numerical integration of differential equations is presented and analyzed. They are constructed as linear combinations of symmetric-conjugate compositions obtained from a basic…
Symplectic integration algorithms are well-suited for long-term integrations of Hamiltonian systems because they preserve the geometric structure of the Hamiltonian flow. However, this desirable property is generally lost when adaptive…
It has previously been shown that varying the numerical timestep during a symplectic orbital integration leads to a random walk in energy and angular momentum, destroying the phase space-conserving property of symplectic integrators. Here…
The isotropic harmonic oscillator in dimension 3 separates in several different coordinate systems. Separating in a particular coordinate system defines a system of three commuting operators, one of which is the Hamiltonian. We show that…
Two families of symplectic methods specially designed for second-order time-dependent linear systems are presented. Both are obtained from the Magnus expansion of the corresponding first-order equation, but otherwise they differ in…
Using the exact path integral solution for the damped harmonic oscillator it is shown that in general there does not exist an exact dissipative Liouville operator describing the dynamics of the oscillator for arbitrary initial bath…
Time-reversible symplectic methods, which are precisely compatible with Liouville's phase-volume-conservation theorem, are often recommended for computational simulations of Hamiltonian mechanics. Lack of energy drift is an apparent…
We construct symplectic integrators for Lie-Poisson systems. The integrators are standard symplectic (partitioned) Runge--Kutta methods. Their phase space is a symplectic vector space with a Hamiltonian action with momentum map $J$ whose…
In this paper, we study one-loop contributions in the double-scaling limit of the SYK model from the chord diagrams and Liouville type effective action. We compute and clarify the meaning of each component consisting of the one-loop…
An exact invariant operator of time-dependent coupled oscillators is derived using the Liouville-von Neumann equation. The unitary relation between this invariant and the invariant of two uncoupled simple harmonic oscillators is…
We present a method to construct symplecticity-preserving renormalization group maps by using the Liouville operator, and obtain correctly reduced symplectic maps describing their long-time behavior even when a resonant island chain…
Symplectic N-body integrators are widely used to study problems in celestial mechanics. The most popular algorithms are of 2nd and 4th order, requiring 2 and 6 substeps per timestep, respectively. The number of substeps increases rapidly…
In a recent work of Wu, Wang, Sun and Liu, a second-order explicit symplectic integrator was proposed for the integrable Kerr spacetime geometry. It is still suited for simulating the nonintegrable dynamics of charged particles moving…