Related papers: A novel efficient structure-preserving exponential…
We develop a general framework for numerically solving differential equations while preserving invariants. As in standard projection methods, we project an arbitrary base integrator onto an invariant-preserving manifold, however, our method…
Casimir preserving integrators for stochastic Lie-Poisson equations with Stratonovich noise are developed extending Runge-Kutta Munthe-Kaas methods. The underlying Lie-Poisson structure is preserved along stochastic trajectories. A related…
Quantum Hamiltonians containing nonseparable products of non-commuting operators, such as $\hat{\bf x}^m \hat{\bf p}^n$, are problematic for numerical studies using split-operator techniques since such products cannot be represented as a…
The most advanced techniques using fault-tolerant quantum computers to estimate the ground-state energy of a chemical Hamiltonian involve compression of the Coulomb operator through tensor factorizations, enabling efficient block-encodings…
Numerical evolution of time-dependent differential equations via explicit Runge-Kutta or Taylor methods typically fails to preserve symmetries of a system. It is known that there exists no numerical integration method that in general…
This paper contributes with a new formal method of spatial discretization of a class of nonlinear distributed parameter systems that allow a port-Hamiltonian representation over a one dimensional manifold. A specific finite dimensional…
Efficient fourth order symplectic integrators are proposed for numerical integration of separable Hamiltonian systems H(p,q)=T(p)+V(q). Symmetric splitting coefficients with five to nine stages are obtained by higher order decomposition of…
Kahan introduced an explicit method of discretization for systems of first order differential equations with nonlinearities of degree at most two (quadratic vector fields). Kahan's method has attracted much interest due to the fact that it…
The goal of this project is to compare the performance of exponential time integrators with traditional methods such as diagonally implicit Runge-Kutta methods in the context of solving the system of reduced magnetohydrodynamics (RMHD). In…
In this paper, we investigate the long-time near-conservations of energy and kinetic energy by the widely used exponential integrators to highly oscillatory conservative systems. The modulated Fourier expansions of two kinds of exponential…
In this paper, we present a novel class of high-order energy-preserving schemes for solving the Zakharov-Rubenchik equations. The main idea of the scheme is first to introduce an quadratic auxiliary variable to transform the Hamiltonian…
In this paper, we develop Hamiltonian particle-in-cell methods for Vlasov-Maxwell equations by applying conforming finite element methods in space and splitting methods in time. For the spatial discretisation, the criteria for choosing…
A new family of methods involving complex coefficients for the numerical integration of differential equations is presented and analyzed. They are constructed as linear combinations of symmetric-conjugate compositions obtained from a basic…
The main objective of this paper is to present an efficient structure-preserving scheme, which is based on the idea of the scalar auxiliary variable approach, for solving the space fractional nonlinear Schr\"{o}dinger equation. First, we…
Symplectic integrators are widely implemented numerical integrators for Hamiltonian mechanics, which preserve the Hamiltonian structure (symplecticity) of the system. Although the symplectic integrator does not conserve the energy of the…
Numerical algorithms based on variational and symplectic integrators exhibit special features that make them promising candidates for application to general relativity and other constrained Hamiltonian systems. This paper lays part of the…
We present a method for the calculation of electronic structure of systems that contain tens of thousands of atoms. The method is based on the division of the system into mutually overlapping fragments and the representation of the…
In this paper, we will prove a very general result of stability for perturbations of linear integrable Hamiltonian systems, and we will construct an example of instability showing that both our result and our example are optimal. Moreover,…
Hamiltonian structure is pursued and uncovered in collisional and collisionless gyrokinetic theory. A new Hamiltonian formulation of collisionless electromagnetic theory is presented that is ideally suited to implementation on modern…
We propose a Lax equation for the non-linear sigma model which leads directly to the conserved local charges of the system. We show that the system has two infinite sets of such conserved charges following from the Lax equation, much like…