Related papers: REBOUNDx: A Library for Adding Conservative and Di…
We introduce a novel numerical method to integrate partial differential equations representing the Hamiltonian dynamics of field theories. It is a multi-symplectic integrator that locally conserves the stress-energy tensor with an excellent…
The Bopp's shifts will be generalized through symplectic formalism. A special procedure, like a "diagonalization", which drives the completely deformed symplectic matrix to the standard symplectic form was found as suggested by…
In Newtonian gravity the final states of cold dissipationless collapses are characterized by several structural and dynamical properties remarkably similar to those of observed elliptical galaxies. Are these properties a peculiarity of the…
The phase space of a Hamiltonian system is symplectic. However, the post-Newtonian Hamiltonian formulation of spinning compact binaries in existing publications does not have this property, when position, momentum and spin variables $[X, P,…
Many experimental techniques aim at determining the Hamiltonian of a given system. The Hamiltonian describes the system's evolution in the absence of dissipation, and is often central to control or interpret an experiment. Here, we…
The conservation of energy, linear momentum and angular momentum are important drivers for our physical understanding of the evolution of the Universe. These quantities are also conserved in Newton's laws of motion under gravity…
The modeling and simulation of infinite-dimensional Hamiltonian systems are central problems in mathematical physics and engineering, however they pose significant computational and structural challenges for standard data-driven…
We consider the direct product of two symplectomorphisms, one of which exhibits a basic set and the other one a non-degenerate elliptic equilibrium. Under a domination condition we show that a broad class of real-analytic deformations of…
We propose a systematic method of dealing with the canonical constrained structure of reducible systems in the Dirac and symplectic approaches which involves an enlargement of phase and configuration spaces, respectively. It is not…
We present a framework for general relativistic N-body simulations in the regime of weak gravitational fields. In this approach, Einstein's equations are expanded in terms of metric perturbations about a Friedmann-Lema\^itre background,…
In a recent work of Wu, Wang, Sun and Liu, a second-order explicit symplectic integrator was proposed for the integrable Kerr spacetime geometry. It is still suited for simulating the nonintegrable dynamics of charged particles moving…
We present a fragmentation module and a composition tracking code for the $n$-body code REBOUND. Our fragmentation code utilises previous semi-analytic models and follows an implementation method similar to fragmentation for the $n$-body…
In a previous paper, second- and fourth-order explicit symplectic integrators were designed for a Hamiltonian of the Schwarzschild black hole. Following this work, we continue to trace the possibility of the construction of explicit…
Modified Hamiltonians are used in the field of geometric numerical integration to show that symplectic schemes for Hamiltonian systems are accurate over long times. For nonlinear systems the series defining the modified Hamiltonian usually…
We suggest a numerical integration procedure for solving the equations of motion of certain classical spin systems which preserves the underlying symplectic structure of the phase space. Such symplectic integrators have been successfully…
Generalized Additive Runge-Kutta schemes have shown to be a suitable tool for solving ordinary differential equations with additively partitioned right-hand sides. This work develops symplectic GARK schemes for additively partitioned…
In N-body simulations of collisionless stellar systems, the forces are softened to reduce the shot noise. Softening modifies gravity at r=|x-y| smaller than softening length epsilon and the softened forces are increasingly biased for ever…
Recently a number of analytic prescriptions for computing the non-linear matter power spectrum have appeared in the literature. These typically involve resummation or closure prescriptions which do not have a rigorous error control, thus…
We employ a port-Hamiltonian approach to model nonlinear rigid multibody systems subject to both position and velocity constraints. Our formulation accommodates Cartesian and redundant coordinates, respectively, and captures kinematic as…
The non-commutative nature of 3D rotations poses well-known challenges in generalizing planar problems to three-dimensional ones, even more so in contact-rich tasks where haptic information (i.e., forces/torques) is involved. In this sense,…