Related papers: Symplectic integration for the collisional gravita…
We show that a recently discovered fourth order symplectic algorithm, which requires one evaluation of force gradient in addition to three evaluations of the force, when iterated to higher order, yielded algorithms that are far superior to…
In this work, we introduce a quadratically convergent and dynamically consistent integrator specifically designed for the replicator dynamics. The proposed scheme combines a two-stage rational approximation with a normalization step to…
Hamilton's equations of motion form a fundamental framework in various branches of physics, including astronomy, quantum mechanics, particle physics, and climate science. Classical numerical solvers are typically employed to compute the…
We design an accurate orbital integration scheme for the general N-body problem preserving all the conserved quantities but the angular momentum.This scheme is based on the chain concept (Mikkola & Aarseth 1993) and is regarded as an…
Symplectic integration of autonomous Hamiltonian systems is a well-known field of study in geometric numerical integration, but for non-autonomous systems the situation is less clear, since symplectic structure requires an even number of…
Using a Newtonian model of the Solar System with all 8 planets, we perform extensive tests on various symplectic integrators of high orders, searching for the best splitting scheme for long term studies in the Solar System. These…
We present a multiscale integrator for Hamiltonian systems with slowly varying quadratic stiff potentials that uses coarse timesteps (analogous to what the impulse method uses for constant quadratic stiff potentials). This method is based…
Two specialized algorithms for the numerical integration of the equations of motion of a Brownian walker obeying detailed balance are introduced. The algorithms become symplectic in the appropriate limits, and reproduce the equilibrium…
We present a semi-infinite program (SIP) solver for trajectory optimizations of general articulated robots. These problems are more challenging than standard Nonlinear Program (NLP) by involving an infinite number of non-convex, collision…
Accurate $N$-body simulations of multiple systems such as binaries and triples are essential for understanding the formation and evolution of interacting binaries and binary mergers, including gravitational wave sources, blue stragglers and…
In this paper we describe a new algorithm for the long-term numerical integration of the two-body problem, in which two particles interact under a Newtonian gravitational potential. Although analytical solutions exist in the unperturbed and…
We introduce a symplectic dual quaternion variational integrator(DQVI) for simulating single rigid body motion in all six degrees of freedom. Dual quaternion is used to represent rigid body kinematics and one-step Lie group variational…
We present the symplectic no-core configuration interaction (SpNCCI) framework, in which the nuclear many-body problem is solved a symmetry-adapted basis that explicitly encodes approximate symmetries associated with nuclear collectivity…
Among the single-trajectory Gaussian-based methods for solving the time-dependent Schr\"{o}dinger equation, the variational Gaussian approximation is the most accurate one. In contrast to Heller's original thawed Gaussian approximation, it…
We construct a symplectic integrator for non-separable Hamiltonian systems combining an extended phase space approach of Pihajoki and the symmetric projection method. The resulting method is semiexplicit in the sense that the main time…
We propose explicit symplectic integrators of molecular dynamics (MD) algorithms for rigid-body molecules in the canonical and isothermal-isobaric ensembles. We also present a symplectic algorithm in the constant normal pressure and lateral…
We present a symplectic integrator, based on the canonical midpoint rule, for classical spin systems in which each spin is a unit vector in $\mathbb{R}^3$. Unlike splitting methods, it is defined for all Hamiltonians, and is…
We compare the performances of symplectic and non-symplectic integrators for the computation of normal geodesics and conjugate points in sub-Riemannian geometry at the example of the Martinet case. For this case study we consider first the…
Recently, continuous-time dynamical systems have proved useful in providing conceptual and quantitative insights into gradient-based optimization, widely used in modern machine learning and statistics. An important question that arises in…
In this paper, we present a new variational integrator for problems in Lagrangian mechanics. Using techniques from Galerkin variational integrators, we construct a scheme for numerical integration that converges geometrically, and is…