Related papers: Leapfrog methods for relativistic charged-particle…
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…
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…
A parallel implementation of coupled spin-lattice dynamics in the LAMMPS molecular dynamics package is presented. The equations of motion for both spin only and coupled spin-lattice dynamics are first reviewed, including a detailed account…
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…
This paper describes a fourth-order integration algorithm for the gravitational N-body problem based on discrete Lagrangian mechanics. When used with shared timesteps, the algorithm is momentum conserving and symplectic. We generalize the…
Hamiltonian systems of ordinary and partial differential equations are fundamental mathematical models spanning virtually all physical scales. A critical property for the robustness and stability of computational methods in such systems is…
Hamiltonian and action principle (HAP) formulations of plasma physics are reviewed for the purpose of explaining structure preserving numerical algorithms. Geometric structures associated with and emergent from HAP formulations are…
The relativistic two-body problem is considered for spinless particles subject to an external macroscopic electromagnetic field. When this field is made of the monochromatic superposition of two counter-propagating plane waves (and provided…
We develop a variational framework for addressing two-dimensional non-integrable quantum field theories through the exact structure of their integrable counterparts. Concentrating on the $\varphi^4$ Landau-Ginzburg model, we use the…
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…
Understanding the contact dynamics of nonspherical particles beyond the microscale is crucial for accurately modeling colloidal and granular systems, where shape anisotropy dictates structural organization and transport properties. In this…
The semi-implicit leapfrog time-discretization is a workhorse algorithm for initial-value MHD codes to bridge between vastly separated time scales. Inclusion of atomic interactions with neutrals breaks the functional structure of the MHD…
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…
The detection of gravitational waves has opened a new era for astronomy, allowing for the combined use of gravitational wave and electromagnetic emissions to directly probe the physics of compact objects, still poorly understood. So far,…
We study the optimal design of numerical integrators for dissipative systems, for which there exists an underlying thermodynamic structure known as GENERIC (general equation for the nonequilibrium reversible-irreversible coupling). We…
A new approach to the electroweak properties of two--particle composite systems is developed. The approach is based on the use of the instant form of relativistic Hamiltonian dynamics. The main novel feature of this approach is the new…
A complete exposition of the rest-frame instant form of dynamics for arbitrary isolated systems (particles, fields, strings, fluids)admitting a Lagrangian description is given. The starting point is the parametrized Minkowski theory…
We develop a geometric framework for the numerical integration of mechanical systems evolving on manifolds. After briefly reviewing classical numerical methods and highlighting their limitations and shortcomings in non-flat (non-Euclidean)…
Reduced magnetohydrodynamics is a simplified set of magnetohydrodynamics equations with applications to both fusion and astrophysical plasmas, possessing a noncanonical Hamiltonian structure and consequently a number of conserved…
This paper proposes a new Linear Fractional Transformation (LFT) modeling approach for uncertain Linear Parameter Varying (LPV) multibody systems with parameter-dependent equilibrium. Traditional multibody approaches, which consist in…