Related papers: N-body Integrators with Individual Time Steps from…
Hamiltonian Monte Carlo (HMC) is a powerful tool for Bayesian statistical inference due to its potential to rapidly explore high dimensional state space, avoiding the random walk behavior typical of many Markov Chain Monte Carlo samplers.…
We present sixth- and eighth-order Hermite integrators for astrophysical $N$-body simulations, which use the derivatives of accelerations up to second order ({\it snap}) and third order ({\it crackle}). These schemes do not require previous…
This work deals with the numerical approximation of plasmas which are confined by the effect of a fast oscillating magnetic field (see \cite{Bostan2012}) in the Vlasov model. The presence of this magnetic field induces oscillations (in…
We here investigate the efficient implementation of the energy-conserving methods named Hamiltonian Boundary Value Methods (HBVMs) recently introduced for the numerical solution of Hamiltonian problems. In this note, we describe an…
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…
Finite-dimensional non-canonical Hamiltonian systems arise naturally from Hamilton's principle in phase space. We present a method for deriving variational integrators that can be applied to perturbed non-canonical Hamiltonian systems on…
This paper presents a fast, economical particle-multiple-mesh N-body code optimized for large-N modelling of collisionless dynamical processes, such as black-hole wandering or bar-halo interactions, occurring within isolated galaxies. The…
Ehrenfest dynamics is a useful approximation for ab initio mixed quantum-classical molecular dynamics that can treat electronically nonadiabatic effects. Although a severe approximation to the exact solution of the molecular time-dependent…
We develop an explicit, second-order, variational time integrator for full body dynamics that preserves the momenta of the continuous dynamics, such as linear and angular momenta, and exhibits near-conservation of total energy over…
We propose an efficient method for active particle selection, working with Hermite Individual Time Steps (HITS) scheme in direct N-body simulation code $\varphi$GRAPE. For a simulation with $N$ particles, this method can reduce the…
Symplectic integration algorithms have become popular in recent years in long-term orbital integrations because these algorithms enforce certain conservation laws that are intrinsic to Hamiltonian systems. For problems with large variations…
Calculating the long term solution of ordinary differential equations, such as those of the $N$-body problem, is central to understanding a wide range of dynamics in astrophysics, from galaxy formation to planetary chaos. Because generally…
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…
In this paper, we develop a framework to construct energy-preserving methods for multi-components Hamiltonian systems, combining the exponential integrator and the partitioned averaged vector field method. This leads to numerical schemes…
The paper deals with numerical discretizations of separable nonlinear Hamiltonian systems with additive noise. For such problems, the expected value of the total energy, along the exact solution, drifts linearly with time. We present and…
We describe a new algorithm for the integration of self-gravitating fluid systems using SPH method. We split the Hamiltonian of a self-gravitating fluid system to the gravitational potential and others (kinetic and internal energies) and…
The kinetic term of the $N$-body Hamiltonian system defined on the surface of the sphere is non-separable. As a result, standard explicit symplectic integrators are inapplicable. We exploit an underlying hierarchy in the structure of the…
We generalize a recently developed method for accelerated Monte Carlo calculation of path integrals to the physically relevant case of generic many-body systems. This is done by developing an analytic procedure for constructing a hierarchy…
In this paper, a new class of energy-preserving integrators is proposed and analysed for Poisson systems by using functionally-fitted technology. The integrators exactly preserve energy and have arbitrarily high order. It is shown that the…
We discuss how dynamical fermion computations may be made yet cheaper by using symplectic integrators that conserve energy much more accurately without decreasing the integration step size. We first explain why symplectic integrators…