Related papers: Step-size effect in the time-transformed leapfrog …
Symplectic integrators separate a problem into parts that can be solved in isolation, alternately advancing these sub-problems to approximate the evolution of the complete system. Problems with a single, dominant mass can use mixed-variable…
The Serre-Green-Naghdi (SGN) equations provide a valuable framework for modelling fully nonlinear and weakly dispersive shallow-water flows. However, their elliptic formulation can considerably increase the computational cost compared to…
We study the time-asymptotic behavior of linear hyperbolic systems under partial dissipation which is localized in suitable subsets of the domain. More precisely, we recover the classical decay rates of partially dissipative systems…
We study a variation of the Trotter-Suzuki decomposition, in which a Hamiltonian exponential is approximated by an ordered product of two-qubit operator exponentials such that the Trotter step size is enhanced for a small number of terms.…
Self-force methods can be applied in calculations of the scatter angle in two-body hyperbolic encounters, working order by order in the mass ratio (assumed small) but with no recourse to a weak-field approximation. This, in turn, can inform…
It is well known that de Sitter(dS) black holes generally have a black hole horizon and a cosmological horizon, both of which have Hawking radiation. But the radiation temperature of the two horizons is generally different, so dS black…
In this paper, we develop a framework to construct energy-preserving methods for multi-components Hamiltonian systems, combining the exponential integrator and the partitioned averaged vector field method. This leads to numerical schemes…
We propose and analyze a second-order partitioned time-stepping method for a two-phase flow problem in porous media. The algorithm is based on a refactorization of Cauchy's one-leg $\theta$-method. The main part consists of the implicit…
In reinforcement learning, the TD($\lambda$) algorithm is a fundamental policy evaluation method with an efficient online implementation that is suitable for large-scale problems. One practical drawback of TD($\lambda$) is its sensitivity…
We introduce energy-preserving integrators for nonholonomic mechanical systems. We will see that the nonholonomic dynamics is completely determined by a triple $({\mathcal D}^*, \Pi, \mathcal{H})$, where ${\mathcal D}^*$ is the dual of the…
We study a simple model of interacting bosons on a d-dimensional cubic lattice whose dynamics conserves both total boson number and total boson dipole moment. This model provides a simple framework in which several remarkable consequences…
Symplectic integrators offer many advantages for the numerical solution of Hamiltonian differential equations, including bounded energy error and the preservation of invariant sets. Two of the central Hamiltonian systems encountered in…
A 3d electron topological insulator (ETI) is a phase of matter protected by particle-number conservation and time-reversal symmetry. It was previously believed that the surface of an ETI must be gapless unless one of these symmetries is…
The primary objective of this paper is to present a long-term numerical energy-preserving analysis of one-stage explicit symmetric and/or symplectic extended Runge--Kutta--Nystr\"{o}m (ERKN) integrators for highly oscillatory Hamiltonian…
In this paper, we analyze the preservation of asymptotic properties of partially dissipative hyperbolic systems when switching to a discrete setting. We prove that one of the simplest consistent and unconditionally stable numerical methods…
Current theoretical studies of electronic correlations in transition metal oxides typically only account for the local repulsion between d-electrons even if oxygen ligand p-states are an explicit part of the effective Hamiltonian.…
Dynamic simulation of elastic bodies is a longstanding task in engineering and computer graphics. In graphics, numerical integrators like implicit Euler and BDF2 are preferred due to their stability at large time steps, but they tend to…
In this work, the benefits of the phase fitting technique are embedded in high order discrete Lagrangian integrators. The proposed methodology creates integrators with zero phase lag in a test Lagrangian in a similar way used in phase…
In this work we propose a new numerical approach to distinguish between regular and chaotic orbits in Hamiltonian systems, based on the simultaneous integration of both the orbit and the deviation vectors using a symplectic scheme, hereby…
We present a new mixed variable symplectic (MVS) integrator for planetary systems, that fully resolve close encounters. The method is based on a time regularisation that allows keeping the stability properties of the symplectic integrators,…