Related papers: N-body Integrators with Individual Time Steps from…
We present a practical algorithm based on symplectic splitting methods to integrate numerically in time the Schr\"odinger equation. When discretized in space, the Schr\"odinger equation can be recast as a classical Hamiltonian system…
A method for calculating phase-space master integrals for the decay process $1 \to n$ massless partons in QCD using integration-by-parts and differential equations techniques is discussed. The method is based on the appropriate choice of…
This paper proposes several explicit and implicit multistep frequency response optimized integrators considering first or second order derivative. A prediction-based method aiming at accelerating a novel power system transient simulation…
We propose an efficient approach for time integration of Klein-Gordon equations with highly oscillatory in time input terms. The new methods are highly accurate in the entire range, from slowly varying up to highly oscillatory regimes. Our…
Symplectic integration methods based on operator splitting are well established in many branches of science. For Hamiltonian systems which split in more than two parts, symplectic methods of higher order have been studied in detail only for…
We change a previous time-stepping algorithm for solving a multi-scale Vlasov-Poisson system within a Particle-In-Cell method, in order to do accurate long time simulations. As an exponential integrator, the new scheme allows to use large…
We present a numerical method to simulate the time evolution, according to a Hamiltonian made of local interactions, of quantum spin chains and systems alike. The efficiency of the scheme depends on the amount of the entanglement involved…
An algorithm for a family of self-starting high-order implicit time integration schemes with controllable numerical dissipation is proposed for both linear and nonlinear transient problems. This work builds on the previous works of the…
Direct N-body simulations and symplectic integrators are effective tools to study the long-term evolution of planetary systems. The Wisdom-Holman (WH) integrator in particular has been used extensively in planetary dynamics as it allows for…
We discuss the short-time perturbative expansion of the linear entropy for finite-dimensional quantum systems whose dynamics can be effectively described by a non-Hermitian Hamiltonian. We derive a timescale for the degree of mixedness for…
We address the problem of simulating pair-interaction Hamiltonians in n node quantum networks where the subsystems have arbitrary, possibly different, dimensions. We show that any pair-interaction can be used to simulate any other by…
Quantum Monte Carlo simulations are powerful and versatile tools for the quantum many-body problem. In addition to the usual calculations of energies and eigenstate observables, quantum Monte Carlo simulations can in principle be used to…
Respecting the laws of thermodynamics is crucial for ensuring that numerical simulations of dynamical systems deliver physically relevant results. In this paper, we construct a structure-preserving and thermodynamically consistent finite…
We review our results for the dynamics of isolated many-body quantum systems described by one-dimensional spin-1/2 models. We explain how the evolution of these systems depends on the initial state and the strength of the perturbation that…
Modeling self-gravity of collisionless fluids (e.g. ensembles of dark matter, stars, black holes, dust, planetary bodies) in simulations is challenging and requires some force softening. It is often desirable to allow softenings to evolve…
In this paper, we formulate and analyse exponential integrations when applied to nonlinear Schr\"{o}dinger equations in a normal or highly oscillatory regime. A kind of exponential integrators with energy preservation, optimal convergence…
Simulation of contact and friction dynamics is an important basis for control- and learning-based algorithms. However, the numerical difficulties of contact interactions pose a challenge for robust and efficient simulators. A…
We introduce a class of fourth order symplectic algorithms that are ideal for doing long time integration of gravitational few-body problems. These algorithms have only positive time steps, but require computing the force gradient in…
We present implementations of a fourth-order symplectic integrator on graphic processing units for three $N$-body models with long-range interactions of general interest: the Hamiltonian Mean Field, Ring and two-dimensional self-gravitating…
We develop a new relativistic radiation hydrodynamics code based on the Monte-Carlo algorithm. In this code, we implement a new scheme to achieve the second-order accuracy in time in the limit of a large packet number for solving the…