Related papers: An optimum time-stepping scheme for N-body simulat…
Particle-mesh simulations trade small-scale accuracy for speed compared to traditional, computationally expensive N-body codes in cosmological simulations. In this work, we show how a data-driven model could be used to learn an effective…
(Abridged) The growing role played by numerical N-body simulations in cosmological studies as a fundamental connection between theoretical modeling and direct observations has led to impressive advancements also in the development and…
Cosmological N-body simulations aim to calculate the non-linear gravitational growth of structures via particle dynamics. A crucial problem concerns the setting-up of the initial particle distribution, as standard theories of galaxy…
Beyond convergence studies and comparison of different codes, there are essentially no controls on the accuracy in the non-linear regime of cosmological N body simulations, even in the dissipationless limit. We propose and explore here a…
Mounting discoveries of debris discs orbiting newly-formed stars and white dwarfs (WDs) showcase the importance of modeling the long-term evolution of small bodies in exosystems. WD debris discs are in particular thought to form from very…
We present a new implicit asymptotic preserving time integration scheme for charged-particle orbit computation in arbitrary electromagnetic fields. The scheme is built on the Crank-Nicolson integrator and continues to recover full-orbit…
In this review we discuss Cosmological N-Body codes with a special emphasis on Particle Mesh codes. We present the mathematical model for each component of N-Body codes. We compare alternative methods for computing each quantity by…
We present the molecular hyperdynamics algorithm and its implementation to the nonorthogonal tight-binding model NTBM and the corresponding software. Due to its multiscale structure, the proposed approach provides the long time scale…
Background The development of a simulation model of full body reaching tasks that can predict endeffector trajectories and joint excursions consistent with experimental data is a non-trivial task. Because of the kinematic redundancy…
On the base of a Feynman-Kac--type formula involving Poisson stochastic processes, recently a Monte Carlo algorithm has been introduced, which describes exactly the real- or imaginary-time evolution of many-body lattice quantum systems. We…
We introduce Spiral, a third-order integration algorithm for the rotational motion of extended bodies. It requires only one force calculation per time step, does not require quaternion normalization at each time step, and can be formulated…
We present two types of meta-algorithm that can greatly improve the accuracy of existing algorithms for integrating the equations of motion of dynamical systems. The first meta-algorithm takes an integrator that is time-symmetric only for…
Binary collisions in Fermi systems obey two fundamental symmetries corresponding to the space and time inversion and to the interchange of particles and holes. We show that beyond the local and instant approximation of scattering-in and…
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…
Hard spheres are arguably one of the most fundamental model systems in soft matter physics, and hence a common topic of simulation studies. Event-driven simulation methods provide an efficient method for studying the phase behavior and…
Performing highly agile acrobatic motions with a long flight phase requires perfect timing, high accuracy, and coordination of the full-body motion. To address these challenges, we present a novel approach on timings and trajectory…
Recent detector concepts at future linear or circular $e^- e^+$ colliders emphasize the benefits of time-of-flight measurements for particle identification of long-lived charged hadrons. That method relies on a precise estimation of the…
First-order variational equations are widely used in N-body simulations to study how nearby trajectories diverge from one another. These allow for efficient and reliable determinations of chaos indicators such as the Maximal Lyapunov…
In this contribution we derive and analyze a new numerical method for kinetic equations based on a variable transformation of the moment approximation. Classical minimum-entropy moment closures are a class of reduced models for kinetic…
We present a non-iterative algorithm, FloatStepper, for coupling the motion of a rigid body and an incompressible fluid in computational fluid dynamics (CFD) simulations. The purpose of the algorithm is to remove the so-called added mass…