Related papers: Symplectic Integrator Mercury: Bug Report
Preservation of coherence is a fundamental yet subtle phenomenon in open systems. We uncover its relation to symmetries respected by the system Hamiltonian and its coupling to the environment. We discriminate between local and global…
We consider the numerical simulation of Hamiltonian systems of ordinary differential equations. Two features of Hamiltonian systems are that energy is conserved along trajectories and phase space volume is preserved by the flow. We want to…
This paper presents a general and robust method for the fluid-structure interaction of membranes and shells undergoing large displacement and large added-mass effects by coupling an immersed-boundary method with a shell finite-element…
I begin with a brief history of N-body simulation and visualization and then go on to describe various methods for creating images and animations of modern simulations in cosmology and galactic dynamics. These techniques are incorporated…
We advocate the use of on-shell constrained $M_2$ variables in order to mitigate the combinatorial problem in SUSY-like events with two invisible particles at the LHC. We show that in comparison to other approaches in the literature, the…
Symmetry is one of the most central concepts in physics, and it is no surprise that it has also been widely adopted as an inductive bias for machine-learning models applied to the physical sciences. This is especially true for models…
Optimization time integrators are effective at solving complex multi-physics problems including deformable solids with non-linear material models, contact with friction, strain limiting, etc. For challenging problems, Newton-type optimizers…
The study of controlled hybrid systems requires practical tools for approximation and comparison of system behaviors. Existing approaches to these problems impose undue restrictions on the system's continuous and discrete dynamics.…
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…
Spin models featuring infinite-range, homogeneous all-to-all interactions can be efficiently described due to the existence of a symmetry-restricted Hilbert subspace and an underlying classical phase space structure. However, when the…
We study by numerical simulations the transition from the metastable "false vacuum" to the broken symmetry phase in the preheating stage after cosmic inflation in a scalar hybrid model. We take quantum fluctuations and their back reaction…
We identify ghost trajectories of symplectic-energy-momentum (SEM) integration and show that the ghost trajectories are not time reversible. We explain how SEM integration can be regularized, in a SEM preserving manner, so that it is time…
The ongoing conjecture that the presence of horizon may induce chaos in an integrable system, is further investigated from the perspective of a uniformly accelerated frame. Particularly, we build up a model which consists of a particle…
A procedure to numerically integrate non-autonomous linear delay differential equations is presented. It is based on the use of an spectral discretization of the delayed part to transform the original problem into a matrix linear ordinary…
The problem of nonintegrability of the circular restricted three-body problem is very classical and important in the theory of dynamical systems. It was partially solved by Poincare in the nineteenth century: He showed that there exists no…
Coherent or exact equations of motion for a post-Newtonian Lagrangian formalism are the Euler-Lagrange equations without any terms truncated. They naturally conserve energy {and} angular momentum. Doubling the phase-space variables of…
The interval approach to computation of dynamics of celestial bodies in the planetary problem has been considered. It is based on the refusal from idealization of infinitely high resolving capacity of measuring tools, and forms an…
Chaotic systems such as the gravitational N-body problem are ubiquitous in astronomy. Machine learning (ML) is increasingly deployed to predict the evolution of such systems, e.g. with the goal of speeding up simulations. Strategies such as…
A general primal-dual splitting algorithm for solving systems of structured coupled monotone inclusions in Hilbert spaces is introduced and its asymptotic behavior is analyzed. Each inclusion in the primal system features compositions with…
We study numerically classical 1-dimensional Hamiltonian lattices involving inter-particle long range interactions that decay with distance like 1/r^alpha, for alpha>=0. We demonstrate that although such systems are generally characterized…