Related papers: A Spectral Symplectic Algorithm for Cylindrical El…
Solving quaternion kinematical differential equations is one of the most significant problems in the automation, navigation, aerospace and aeronautics literatures. Most existing approaches for this problem neither preserve the norm of…
We introduce a new electrostatic particle-in-cell algorithm capable of using large timesteps compared to particle gyro-period under a uniform external magnetic field. The algorithm extends earlier electrostatic fully implicit PIC…
Geometric numerical integration has recently been exploited to design symplectic accelerated optimization algorithms by simulating the Lagrangian and Hamiltonian systems from the variational framework introduced in Wibisono et al. In this…
In this work, a kernel-based surrogate for integrating Hamiltonian dynamics that is symplectic by construction and tailored to large prediction horizons is proposed. The method learns a scalar potential whose gradient enters a…
Metriplectic dynamical systems consist of a special combination of a Hamiltonian and a (generalized) entropy-gradient flow, such that the Hamiltonian is conserved and entropy is dissipated/produced (depending on a sign convention). It is…
The intention of this article is to illustrate the use of methods from symplectic geometry for practical purposes. Our intended audience is scientists interested in orbits of Hamiltonian systems (e.g. the three-body problem). The main…
Retraction maps have been generalized to discretization maps in (Barbero Li\~n\'an and and Mart\'{\i}n de Diego, 2022). Discretization maps are used to systematically derive numerical integrators that preserve the symplectic structure, as…
A novel electromagnetic particle-in-cell algorithm has been developed for fully kinetic plasma simulations on unstructured (irregular) meshes in complex body-of-revolution geometries. The algorithm, implemented in the BORPIC++ code,…
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…
An infinite dimensional canonical symplectic structure and structure-preserving geometric algorithms are developed for the photon-matter interactions described by the Schr\"odinger-Maxwell equations. The algorithms preserve the symplectic…
We present a practical algorithm based on symplectic splitting methods to integrate numerically in time the Schr\"odinger equation. When discretized in space, the Schr\"odinger equation can be recast as a classical Hamiltonian system…
We present a symplectic integrator, based on the canonical midpoint rule, for classical spin systems in which each spin is a unit vector in $\mathbb{R}^3$. Unlike splitting methods, it is defined for all Hamiltonians, and is…
This article considers non-relativistic charged particle dynamics in both static and non-static electromagnetic fields, which are governed by nonseparable, possibly time-dependent Hamiltonians. For the first time, explicit symplectic…
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…
Representing the electrodynamics of relativistically drifting particle ensembles in discrete, co-propagating Galilean coordinates enables the derivation of a Particle-in-Cell algorithm that is intrinsically free of the Numerical Cherenkov…
Variational integrators are derived for structure-preserving simulation of stochastic Hamiltonian systems with a certain type of multiplicative noise arising in geometric mechanics. The derivation is based on a stochastic discrete…
Predicting the behaviors of Hamiltonian systems has been drawing increasing attention in scientific machine learning. However, the vast majority of the literature was focused on predicting separable Hamiltonian systems with their kinematic…
Singular theories, characterised by the presence of degeneracies in their Lagrangian or Hamiltonian descriptions, require the systematic implementation of constraints in order to obtain well-defined dynamics. While the symplectic framework…
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…
Symplectic integration methods based on operator splitting are well established in many branches of science. For Hamiltonian systems which split in more than two parts, symplectic methods of higher order have been studied in detail only for…