Related papers: On the accuracy of symplectic integrators for secu…
Long-term stability studies of nonlinear Hamiltonian systems require symplectic integration algorithms which are both fast and accurate. In this paper, we study a symplectic integration method wherein the symplectic map representing the…
This letter proposes a predictor-corrector method to strike a balance between simulation accuracy and efficiency by appropriately tuning the numerical integration step length of a power system time-domain simulation. Numerical tests…
We examine a new 2nd order integrator recently found by Omelyan et al. The integration error of the new integrator measured in the root mean square of the energy difference, $\bra\Delta H^2\ket^{1/2}$, is about 10 times smaller than that of…
Phase fitting has been extensively used during the last years to improve the behaviour of numerical integrators on oscillatory problems. In this work, the benefits of the phase fitting technique are embedded in discrete Lagrangian…
In this note, we propose a symplectic algorithm for the stable manifolds of the Hamilton-Jacobi equations combined with an iterative procedure in [Sakamoto-van~der Schaft, IEEE Transactions on Automatic Control, 2008]. Our algorithm…
The non-resonant secular dynamics of compact planetary systems are modeled by a perturbing function which is usually expanded in eccentricity and absolute inclination with respect to the invariant plane. Here, the expressions are given in a…
Structural symmetries of linear dynamical systems can be exploited for decoupling the dynamics and reducing the computational complexity of the controller implementation. However, in practical applications, inexact structural symmetries…
Symplectic quantization is a functional approach to quantum field theory that allows sampling of quantum fluctuations directly in Minkowski space time by means of a generalized Hamiltonian dynamics in an extra time variable $\tau$ which, at…
Self-interacting dark matter (SIDM) theory predicts that dark matter halos experience core-collapse, a process where the halo's inner region rapidly increases in density and decreases in size. The N-body simulations used to study this…
We propose an effective and lightweight learning algorithm, Symplectic Taylor Neural Networks (Taylor-nets), to conduct continuous, long-term predictions of a complex Hamiltonian dynamic system based on sparse, short-term observations. At…
Among the single-trajectory Gaussian-based methods for solving the time-dependent Schr\"{o}dinger equation, the variational Gaussian approximation is the most accurate one. In contrast to Heller's original thawed Gaussian approximation, it…
Computations of quantum corrections to the CMB spectrum and to scalar field dynamics during inflation very often take advantage of the "semi-classical" approach, where the metric fluctuations are simply omitted. On the other hand, a…
Large-scale cosmological simulations are an indispensable tool for modern cosmology. To enable model-space exploration, fast and accurate predictions are critical. In this paper, we show that the performance of such simulations can be…
In this paper, we consider the harmonic extension problem, which is widely used in many applications of machine learning. We find that the transitional method of graph Laplacian fails to produce a good approximation of the classical…
We design a novel, exactly energy-conserving implicit non-symplectic integration method for an eight-dimensional Hamiltonian system with four degrees of freedom. In our algorithm, each partial derivative of the Hamiltonian with respect to…
A minimal requirement for simulating multi-scale systems is to reproduce the statistical behavior of the slow variables. In particular, a good numerical method should accurately aproximate the probability density function of the…
Numerical methods that preserve geometric invariants of the system, such as energy, momentum or the symplectic form, are called geometric integrators. In this paper we present a method to construct symplectic-momentum integrators for…
We present an approach to construct appropriate and efficient emulators for Hamiltonian flow maps. Intended future applications are long-term tracing of fast charged particles in accelerators and magnetic plasma confinement configurations.…
The possibility of measuring the post-Newtonian gravitoelectric correction to the orbital period of a test particle freely orbiting a spherically symmetric mass in the Solar System is analyzed. It should be possible, in principle, to detect…
In this paper, explicit stable integrators based on symplectic and contact geometries are proposed for a non-autonomous ordinarily differential equation (ODE) found in improving convergence rate of Nesterov's accelerated gradient method.…