Related papers: Symplectic fourth-order maps for the collisional N…
Hybrid symplectic integrators such as MERCURY are widely used to simulate complex dynamical phenomena in planetary dynamics that could otherwise not be investigated. A hybrid integrator achieves high accuracy during close encounters by…
In the Traveling Salesman Problem (TSP), a list of cities and the distances between them are given. The goal is to find the shortest possible route that visits each city exactly once and returns to the original city. The TSP has a wide…
This article considers Hamiltonian mechanical systems with potential functions admitting jump discontinuities. The focus is on accurate and efficient numerical approximations of their solutions, which will be defined via the laws of…
A system of relativistic Snyder particles with mutual two-body interaction that lives in a Non-Commutative Snyder geometry is studied. The underlying novel symplectic structure is a coupled and extended version of (single particle) Snyder…
Inverse problems are in many cases solved with optimization techniques. When the underlying model is linear, first-order gradient methods are usually sufficient. With nonlinear models, due to nonconvexity, one must often resort to…
We consider the reliable implementation of an adaptive high-order unfitted finite element method on Cartesian meshes for solving elliptic interface problems with geometrically curved singularities. We extend our previous work on the…
Quantum annealing provides a promising way to solve combinational optimization problems where the solutions correspond to the ground state of the Ising Hamiltonian. We can implement quantum annealing using the Kerr non-linear resonators,…
We investigate the performance of various methods of symplectic integration, which are based on two part splitting of the integration operator, for the numerical integration of multidimensional Hamiltonian systems. We implement these…
The optimal mass transport problem gives a geometric framework for optimal allocation, and has recently gained significant interest in application areas such as signal processing, image processing, and computer vision. Even though it can be…
This paper focuses on the parallel implementation of a direct $N$-body method~(particle-particle algorithm) and the application of multiple GPUs for galactic dynamics simulations. Application of a hybrid OpenMP-CUDA technology is considered…
We show (Theorem 3) that the symplectic reduction of the spatial $n$-body problem at non-zero angular momentum is a singular symplectic space consisting of two symplectic strata, one for spatial motions and the other for planar motions.…
Time-symmetric integration schemes share with symplectic schemes the property that their energy errors show a much better behavior than is the case for generic integration schemes. Allowing adaptive time steps typically leads to a loss of…
Two families of symplectic methods specially designed for second-order time-dependent linear systems are presented. Both are obtained from the Magnus expansion of the corresponding first-order equation, but otherwise they differ in…
Efficiently simulating many-body quantum systems with Coulomb interactions is a fundamental question in quantum physics, quantum chemistry, and quantum computing, yet it presents unique challenges: the Hamiltonian is an unbounded operator…
We consider the numerical integration of the matrix Hill's equation. Parametric resonances can appear and this property is of great interest in many different physical applications. Usually, the Hill's equations originate from a Hamiltonian…
We present a set of new, efficient high-order symplectic methods designed for Hamiltonian systems with cubic or quartic potentials. By demonstrating that polynomial potentials require fewer order conditions, we develop schemes that…
This paper introduces the Circular Restricted n-Body Problem (CRNBP), an extension of the bicircular restricted four-body problem (BCR4BP) designed to describe the dynamics of an n-body system. In the CRNBP, each massive body in the system…
Recently a new class of numerical integration methods -- ``mixed variable symplectic integrators'' -- has been introduced for studying long-term evolution in the conservative gravitational few-body problem. These integrators are an order of…
In this work, we present a new approach to the construction of variational integrators. In the general case, the estimation of the action integral in a time interval $[q_k,q_{k+1}]$ is used to construct a symplectic map $(q_k,q_{k+1})\to…
We generalize the coupled-cluster (CC) approach with singles, doubles, and the non-iterative treatment of triples termed $\Lambda$CCSD(T) to Hamiltonians containing three-body interactions. The resulting method and the underlying CC…