Related papers: High order symplectic integrators for perturbed Ha…
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 show here that the Hamiltonian for an electronic system may be written exactly in terms of fluctuation operators that transition constituent fragments between internally correlated states, accounting rigorously for inter-fragment…
We prove that a class of A-stable symplectic Runge--Kutta time semidiscretizations (including the Gauss--Legendre methods) applied to a class of semilinear Hamiltonian PDEs which are well-posed on spaces of analytic functions with analytic…
It is necessary to calculate the C operator for the non-Hermitian PT-symmetric Hamiltonian H=\half p^2+\half\mu^2x^2-\lambda x^4 in order to demonstrate that H defines a consistent unitary theory of quantum mechanics. However, the C…
Let M be a connected, symplectic 4-manifold. A semitoric integrable system on M essentially consists of a pair of independent, real-valued, smooth functions J and H on the manifold M, for which J generates a Hamiltonian circle action under…
Integrable Hamiltonian systems on almost-symplectic manifolds have recently drawn some attention. Under suitable properties, they have a structure analogous to those of standard symplectic-Hamiltonian completely integrable systems. Here we…
The symplectic Wisdom-Holman map revolutionized long-term integrations of planetary systems. There is freedom in such methods of how to split the Hamiltonian and which coordinate system to employ, and several options have been proposed in…
High orders of perturbation theory can be calculated by the Lipatov method, whereby they are determined by saddle-point configurations (instantons) of the corresponding functional integrals. For most field theories, the Lipatov asymptotics…
In this paper, we consider exponential integrators for semilinear Poisson systems. Two types of exponential integrators are constructed, one preserves the Poisson structure, and the other preserves energy. Numerical experiments for…
The Lie-Poisson algebra so(N+1) and some of its contractions are used to construct a family of superintegrable Hamiltonians on the ND spherical, Euclidean, hyperbolic, Minkowskian and (anti-)de Sitter spaces. We firstly present a…
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…
The idea of adaptive perturbation theory is to divide a Hamiltonian into a solvable part and a perturbation part. The solvable part contains the non-interacting sector and the diagonal elements of Fock space from the interacting terms. The…
Hill's equations are an approximation that is useful in a number of areas of astrophysics including planetary rings and planetesimal disks. We derive a symplectic method for integrating Hill's equations based on a generalized leapfrog. This…
This paper develops a structure-preserving numerical integration scheme for a class of higher-order mechanical systems. The dynamics of these systems are governed by invariant variational principles defined on higher-order tangent bundles…
The three integrable two-dimensional Henon-Heiles systems and their integrable perturbations are revisited. A family of new integrable perturbations is found, and N-dimensional completely integrable generalizations of all these systems are…
We consider symplectic manifolds with Hamiltonian torus actions which are "almost but not quite completely integrable": the dimension of the torus is one less than half the dimension of the manifold. We provide a complete set of invariants…
We show that applying any deterministic B-series method of order $p_d$ with a random step size to single integrand SDEs gives a numerical method converging in the mean-square and weak sense with order $\lfloor p_d/2\rfloor$.As an…
We present new almost time-reversible integrators for solution of planetary systems consisting of "planets" and a dominant mass ("star"). The algorithms can be considered adaptive generalizations of the Wisdom--Holman method, in which all…
In this work, we show high order splitting methods of integration without negative steps, allowing us to solve numerically irreversible problems, like reaction-diffusion equations. The methods consist in a suitable affine combinations of…
We classify the Hamiltonians $H=p_x^2+p_y^2+V(x,y)$ of all classical superintegrable systems in two dimensional complex Euclidean space with second-order constants of the motion. We similarly classify the superintegrable Hamiltonians…