Related papers: Numerical integration of variational equations for…
An interesting family of geometric integrators for Lagrangian systems can be defined using discretizations of the Hamilton's principle of critical action. This family of geometric integrators is called variational integrators. In this…
The process of relaxation of a system of particles interacting with long-range forces is relevant to many areas of Physics. For obvious reasons, in Stellar Dynamics much attention has been paid to the case of 1/r^2 force law. However,…
We propose a variational symplectic numerical method for the time integration of dynamical systems issued from the least action principle. We assume a quadratic internal interpolation of the state and we approximate the action in a small…
We devote our studies to the subject of weakly nonintegrable dynamics of systems with a macroscopic number of degrees of freedom. Our main points of interest are the relations between the timescales of thermalization and the timescales of…
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…
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 introduce and analyze a model for the transport of particles or energy in extended lattice systems. The dynamics of the model acts on a discrete phase space at discrete times but has nonetheless some of the characteristic properties of…
A Hamiltonian formulation of generic many-particle systems with space-dependent balanced loss and gain coefficients is presented. It is shown that the balancing of loss and gain necessarily occurs in a pair-wise fashion. Further, using a…
The Hamiltonian dynamics of the classical $\phi^4$ model on a two-dimensional square lattice is investigated by means of numerical simulations. The macroscopic observables are computed as time averages. The results clearly reveal the…
Recent theoretical studies of statistical mechanical properties of systems with long range interactions are briefly reviewed. In these systems the interaction potential decays with a rate slower than 1/r^d at large distances r in d…
Variational symplectic algorithms have recently been developed for carrying out long-time simulation of charged particles in magnetic fields. As a direct consequence of their derivation from a discrete variational principle, these…
Explicit and semi-explicit geometric integration schemes for dissipative perturbations of Hamiltonian systems are analyzed. The dissipation is characterized by a small parameter $\epsilon$, and the schemes under study preserve the…
We study the problem of efficient integration of variational equations in multi-dimensional Hamiltonian systems. For this purpose, we consider a Runge-Kutta-type integrator, a Taylor series expansion method and the so-called `Tangent Map'…
We present a variational integrator based on the Lobatto quadrature for the time integration of dynamical systems issued from the least action principle. This numerical method uses a cubic interpolation of the states and the action is…
We study periodically driven closed systems with a long-ranged Hamiltonian by considering a generalized Kitaev chain with pairing terms which decay with distance as a power law characterized by exponent $\alpha$. Starting from an initial…
We consider the geometric numerical integration of Hamiltonian systems subject to both equality and "hard" inequality constraints. As in the standard geometric integration setting, we target long-term structure preservation. We…
Numerical methods that preserve geometric invariants of the system, such as energy, momentum or the symplectic form, are called geometric integrators. Variational integrators are an important class of geometric integrators. The general idea…
Using 1200 CPUs of the National Supercomputer TH-A1 and a parallel integral algorithm based on the 3500th-order Taylor expansion and the 4180-digit multiple precision data, we have done a reliable simulation of chaotic solution of Lorenz…
In this paper we consider the general setting for constructing Action Principles for three-dimensional first order autonomous equations. We present the results for some integrable and non-integrable cases of the Lotka-Volterra equation, and…
Dynamics of many-body Hamiltonian systems with long range interactions is studied, in the context of the so called $\alpha-$HMF model. Building on the analogy with the related mean field model, we construct stationary states of the…