Related papers: Time-symmetric integration in astrophysics
We study analytically and experimentally certain symplectic and time-reversible N-body integrators which employ a Kepler solver for each pair-wise interaction, including the method of Hernandez & Bertschinger (2015). Owing to the Kepler…
This paper introduces an adaptive time splitting technique for the solution of stiff evolutionary PDEs that guarantees an effective error control of the simulation, independent of the fastest physical time scale for highly unsteady…
Hamilton's equations are fundamental for modeling complex physical systems, where preserving key properties such as energy and momentum is crucial for reliable long-term simulations. Geometric integrators are widely used for this purpose,…
Time-fractional parabolic equations with a Caputo time derivative are considered. For such equations, we explore and further develop the new methodology of the a-posteriori error estimation and adaptive time stepping proposed in [7]. We…
We show that adaptive time stepping in particle accelerator simulation is an enhancement for certain problems. The new algorithm has been implemented in the OPAL (Object Oriented Parallel Accelerator Library) framework, and is compared to…
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…
Time reversal symmetry occupies a distinctive role in quantum mechanics, fundamentally requiring an anti-unitary operator to ensure a physically consistent representation. As such, the time reversal operator combines a unitary…
Different possible sources are discussed for enhancement of the calculation time when solving ordinary differential equations systems to forecast space objects' motion. This paper presents an approach for building an integrator of ordinary…
The large dynamic range in some astrophysical N-body problems led to the use of adaptive multi-time-steps; however, the search for optimal strategies is still challenging. We numerically quantify the performance of the hierarchical…
We study the role of time-reversal symmetry on the dynamical response of nonlinear optical systems that behave as unidirectional ("one-way") devices. It is shown that lossless nonlinear materials, despite being nonreciprocal, are typically…
Geometric numerical integration has recently been exploited to design symplectic accelerated optimization algorithms by simulating the Lagrangian and Hamiltonian systems from the variational framework introduced in Wibisono et al. In this…
Multi-adaptive Galerkin methods are extensions of the standard continuous and discontinuous Galerkin methods for the numerical solution of initial value problems for ordinary or partial differential equations. In particular, the…
Domain decomposition based time integrators allow the usage of parallel and distributed hardware, making them well-suited for the temporal discretization of parabolic systems, in general, and degenerate parabolic problems, in particular.…
We report an accessible and robust tool for evaluating the effects of Coulomb collisions on a test particle in a plasma that obeys Maxwell-J\"uttner statistics. The implementation is based on the Beliaev-Budker collision integral which…
Following on our previous work [S. Delong and B. E. Griffith and E. Vanden-Eijnden and A. Donev, Phys. Rev. E, 87(3):033302, 2013], we develop temporal integrators for solving Langevin stochastic differential equations that arise in…
While the basic laws of physics seem time-reversal invariant, our understanding of the apparent irreversibility of the macroscopic world is well grounded in the notion of entropy. Because astrophysics deals with the largest structures in…
A Kepler solver is an analytical method used to solve a two-body problem. In this paper, we propose a new correction method by slightly modifying the Kepler solver. The only change to the analytical solutions is that the obtainment of the…
Many force-gradient explicit symplectic integration algorithms have been designed for the Hamiltonian $H=T (\mathbf{p})+V(\mathbf{q})$ with kinetic energy $T(\mathbf{p})=\mathbf{p}^2/2$ in the existing references. When the force-gradient…
One of the most accurate methods for solving the time-dependent Schr\"{o}dinger equation uses a combination of the dynamic Fourier method with the split-operator algorithm on a tensor-product grid. To reduce the number of required grid…
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 :…