Related papers: Anomaly in Symplectic Integrator
We implement several symplectic integrators, which are based on two part splitting, for studying the chaotic behavior of one- and two-dimensional disordered Klein-Gordon lattices with many degrees of freedom and investigate their numerical…
We prove some Liouville properties for sub- and supersolutions of fully nonlinear degenerate elliptic equations in the whole space. Our assumptions allow the coefficients of the first order terms to be large at infinity, provided they have…
The finite-element approach to lattice field theory is both highly accurate (relative errors $\sim1/N^2$, where $N$ is the number of lattice points) and exactly unitary (in the sense that canonical commutation relations are exactly…
We present a new class of high-order variational integrators on Lie groups. We show that these integrators are symplectic, momentum preserving, and can be constructed to be of arbitrarily high-order, or can be made to converge…
This article reviews some integrators particularly suitable for the numerical resolution of differential equations on a large time interval. Symplectic integrators are presented. Their stability on exponentially large time is shown through…
Supersymmetry transformations of first and second order are used to generate Hamiltonians with known spectra departing from the harmonic oscillator with an infinite potential barrier. It is studied also the way in which the eigenfunctions…
The finite-element approach to lattice field theory is both highly accurate (relative errors $\sim 1/N^2$, where $N$ is the number of lattice points) and exactly unitary (in the sense that canonical commutation relations are exactly…
We implement and investigate the numerical properties of a new family of integrators connecting both variants of the symplectic Euler schemes, and including an alternative to the classical symplectic mid-point scheme, with some additional…
The classical Arnold-Liouville theorem describes the geometry of an integrable Hamiltonian system near a regular level set of the moment map. Our results describe it near a nondegenerate singular level set: a tubular neighborhood of a…
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…
The out-of-time-order correlator (OTOC) of simple harmonic oscillator with extra anharmonic (quartic) interaction are calculated by the second quantization method. We obtain the analytic formulas of spectrum, Fock space states and matrix…
We suggest a numerical integration procedure for solving the equations of motion of certain classical spin systems which preserves the underlying symplectic structure of the phase space. Such symplectic integrators have been successfully…
The existence of explicit symplectic integrators for general nonseparable Hamiltonian systems is an open and important problem in both numerical analysis and computing in science and engineering, as explicit integrators are usually more…
Dependable numerical results from long-time simulations require stable numerical integration schemes. For Hamiltonian systems, this is achieved with symplectic integrators, which conserve the symplectic condition and exactly solve for the…
We consider the one-dimensional Dirac equation for the harmonic oscillator and the associated second order separated operators giving the resonances of the problem by complex dilation. The same operators have unique extensions as closed…
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…
We construct several variational integrators--integrators based on a discrete variational principle--for systems with Lagrangians of the form L = L_A + epsilon L_B, with epsilon << 1, where L_A describes an integrable system. These…
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…
In this paper, we are concerned with the construction and analysis of a new class of methods obtained as double jump compositions with complex coefficients and projection on the real axis. It is shown in particular that the new integrators…
Symplectic integrators constructed from Hamiltonian and Lie formalisms are obtained as symplectic maps whose flow follows the exact solution of a "sourrounded" Hamiltonian K = H + h^k H_1. Those modified Hamiltonians depends virtually on…