Related papers: Symplectic fourth-order maps for the collisional N…
Unitary coupled cluster (UCC) theory offers a promising Hermitian alternative to conventional coupled cluster (CC) theory, but its practical implementation is hindered by the non-truncating nature of the Baker-Campbell-Hausdorff (BCH)…
For the planar $N$-body problem, we first introduce a class of moving frame suitable for orbits near central configurations, especially for total collision orbits, which is the main new ingredient of this paper. The moving frame allows us…
In this work, we utilize discrete geometric mechanics to derive a 2nd-order variational integrator so as to simulate rigid body dynamics. The developed integrator is to simulate the motion of a free rigid body and a quad-rotor. We…
The paper presents a new reduction method designed for dynamic contact problems. Recently, we have proposed an efficient reduction scheme for the node-to-node formulation, that leads to Linear Complementarity Problems (LCP). Here, we…
Conservative symmetric second-order one-step schemes are derived for dynamical systems describing various many-body systems using the Discrete Multiplier Method. This includes conservative schemes for the $n$-species Lotka-Volterra system,…
Reliable studies of the long-term dynamics of planetary systems require numerical integrators that are accurate and fast. The challenge is often formidable because the chaotic nature of many systems requires relative numerical error bounds…
By combining a standard symmetric, symplectic integrator with a new step size controller, we provide an integration scheme that is symmetric, reversible and conserves the values of the constants of motion. This new scheme is appropriate for…
We study how inexact nonlinear solvers lead to a loss of exact symplecticity in the Symplectic Euler (SE) and Stormer-Verlet (SV) schemes when applied to general nonseparable Hamiltonian systems. These schemes are implicit and require…
It has been argued that despite remarkable success, existing random matrix theories are not adequate to describe disordered conductors in the metallic regime, due to the presence of certain two-body interactions in the effective Hamiltonian…
The idea of adaptive perturbation theory is to divide a Hamiltonian into a solvable part and a perturbation part. The solvable part contains the non-interacting sector and the diagonal elements of Fock space from the interacting terms. The…
The Trotter-Suzuki approximation leads to an efficient algorithm for solving the time-dependent Schr\"odinger equation. Using existing highly optimized CPU and GPU kernels, we developed a distributed version of the algorithm that runs…
Contact-implicit motion planning-embedding contact sequencing as implicit complementarity constraints-holds the promise of leveraging continuous optimization to discover new contact patterns online. Nevertheless, the resulting optimization,…
We introduce a recent symplectic integration scheme derived for solving physically motivated systems with non-separable Hamiltonians. We show its relevance to Riemannian manifold Hamiltonian Monte Carlo (RMHMC) and provide an alternative to…
Two body tunneling problems are hard to treat analytically due to the incompatibility between tunneling and perturbation theory. The lack of classical solutions of the Euclidean Lagrangian of continuous systems further thwarts…
The simulation of multi-body systems with frictional contacts is a fundamental tool for many fields, such as robotics, computer graphics, and mechanics. Hard frictional contacts are particularly troublesome to simulate because they make the…
Symplectic integrators have made it possible to study the long-term evolution of planetary systems with direct N-body simulations. In this paper we reassess the accuracy of such simulations by running a convergence test on 20Myr…
We present a detailed comparison of several integration schemes applied to the dynamic system consisting of a charged particle on the Kerr background endowed with the axisymmetric electromagnetic test field. In particular, we compare the…
Exponential integrators based on contour integral representations lead to powerful numerical solvers for a variety of ODEs, PDEs, and other time-evolution equations. They are embarrassingly parallelizable and lead to global-in-time…
We give a thorough analysis of the convergence properties of the configuration-interaction method as applied to parabolic quantum dots among other systems, including \emph{a priori} error estimates. The method converges slowly in general,…
The St\"ackel transform is applied to the geodesic motion on Euclidean space, through the harmonic oscillator and Kepler-Coloumb potentials, in order to obtain maximally superintegrable classical systems on N-dimensional Riemannian spaces…