Related papers: JANUS: A bit-wise reversible integrator for N-body…
A new approach is developed to integrate numerically the equations of motion for systems of interacting rigid polyatomic molecules. With the aid of a leapfrog framework, we directly involve principal angular velocities into the integration,…
We present a new time-stepping criterion for N-body simulations that is based on the true dynamical time of a particle. This allows us to follow the orbits of particles correctly in all environments since it has better adaptivity than…
If one has to attain high accuracy over long timescales during the numerical computation of the N-body problem, the method called Lie-integration is one of the most effective algorithms. In this paper we present a set of recurrence…
We present new splitting methods designed for the numerical integration of near-integrable Hamiltonian systems, and in particular for planetary N-body problems, when one is interested in very accurate results over a large time span. We…
Phase fitting has been extensively used during the last years to improve the behaviour of numerical integrators on oscillatory problems. In this work, the benefits of the phase fitting technique are embedded in discrete Lagrangian…
In this paper we propose a systematic method to solve the inverse dynamical problem for a quantum system governed by the von Neumann equation: to find a class of Hamiltonians reproducing a prescribed time evolution of a pure or mixed state…
Hamiltonian trajectories are strictly time-reversible. Any time series of Hamiltonian coordinates {q} satisfying Hamilton's motion equations will likewise satisfy them when played "backwards", with the corresponding momenta changing signs :…
We construct holographic Janus solutions, which describe a conformal interface in the theory of M2-branes, in four-dimensional gauged supergravities using a perturbative method. In particular, we study three Einstein-scalar systems and…
Various many-body models are treated, which describe $N$ points confined to move on a plane circle. Their Newtonian equations of motion ("accelerations equal forces") are integrable, i. e. they allow the explicit exhibition of $N$ constants…
We consider the problem of $n$ points with positive masses interacting pairwise with forces inversely proportional to the distance between them. In particular, it is the classical gravitational, Coulomb or photo-gravitational $n$-body…
We study the axisymmetric interaction of two chemically active Janus particles. By relying on the linearity of the field equations and symmetry arguments, we derive a generic solution for the relative velocity of the particles. We show that…
We present a massive equilibrium simulation of the three-dimensional Ising spin glass at low temperatures. The Janus special-purpose computer has allowed us to equilibrate, using parallel tempering, L=32 lattices down to T=0.64 Tc. We…
Direct gravitational simulations of n-body systems have a time complexity O(n^2), which gets computationally expensive as the number of bodies increases. Distributing this workload to multiple cores significantly speeds up the computation…
We present a simple choice of integration variables that can be used to exploit the near-integrable character of problems in celestial mechanics. The approach is based on the well-known principle of variation of parameters: instead of…
Hamiltonian Monte Carlo provides efficient Markov transitions at the expense of introducing two free parameters: a step size and total integration time. Because the step size controls discretization error it can be readily tuned to achieve…
This article considers Hamiltonian mechanical systems with potential functions admitting jump discontinuities. The focus is on accurate and efficient numerical approximations of their solutions, which will be defined via the laws of…
Symplectic integrators are widely used for long-term integration of conservative astrophysical problems due to their ability to preserve the constants of motion; however, they cannot in general be applied in the presence of nonconservative…
Recently, Hamiltonian neural networks (HNN) have been introduced to incorporate prior physical knowledge when learning the dynamical equations of Hamiltonian systems. Hereby, the symplectic system structure is preserved despite the…
We present IAS15, a 15th-order integrator to simulate gravitational dynamics. The integrator is based on a Gau\ss-Radau quadrature and can handle conservative as well as non-conservative forces. We develop a step-size control that can…
The N-body problem is a classic problem involving a system of N discrete bodies mutually interacting in a dynamical system. At any moment in time there are N*(N - 1)/2 such interactions occurring. This scaling as N^2 leads to computational…