Related papers: Multiple invariants conserving Runge-Kutta type me…
Continuous Galerkin Petrov time discretization scheme is tested on some Hamiltonian systems including simple harmonic oscillator, Kepler's problem with different eccentricities and molecular dynamics problem. In particular, we implement the…
In this paper, we present a new methodology to develop arbitrary high-order structure-preserving methods for solving the quantum Zakharov system. The key ingredients of our method are: (i) the original Hamiltonian energy is reformulated…
We introduce a class of high order accurate, semi-implicit Runge-Kutta schemes in the general setting of evolution equations that arise as gradient flow for a cost function, possibly with respect to an inner product that depends on the…
In this paper, a class of high-order compact finite difference Hermite scheme is presented for the simulation of double-diffusive convection. To maintain linear stability, the convective fluxes are split into positive and negative parts,…
In this paper, we study the extended Hamilton-Jacobi Theory in the context of dynamical systems with symmetries. Given an action of a Lie group $G$ on a manifold $M$ and a $G$-invariant vector field $X$ on $M$, we construct complete…
Based on the recent development of Jacobian-free Lax-Wendroff (LW) approaches for solving hyperbolic conservation laws [Zorio, Baeza and Mulet, Journal of Scientific Computing 71:246-273, 2017], [Carrillo and Par\'es, Journal of Scientific…
In this paper we apply Kahan's nonstandard discretization to three dimensional Lotka-Volterra equations in bi-Hamiltonian form. The periodicity of the solutions and all polynomial and non-polynomial invariants are well preserved in…
Strong Stability Preserving (SSP) time integration schemes maintain stability of the forward Euler method for any initial value problem. However, only a small subset of Runge-Kutta (RK) methods are SSP, and many efficient high-order time…
We study the problem of efficient integration of variational equations in multi-dimensional Hamiltonian systems. For this purpose, we consider a Runge-Kutta-type integrator, a Taylor series expansion method and the so-called `Tangent Map'…
The variational quantum-classical algorithms are the most promising approach for achieving quantum advantage on near-term quantum simulators. Among these methods, the variational quantum eigensolver has attracted a lot of attention in…
This work concerns a study of the quantum mechanical extension of the work of Horwitz et al. [1] on the stability of classical Hamiltonian systems by geometrical methods. Simulations are carried out for several important examples, these…
In this paper, we extend the Paired-Explicit Runge-Kutta schemes by Vermeire et. al. to fourth-order of consistency. Based on the order conditions for partitioned Runge-Kutta methods we motivate a specific form of the Butcher arrays which…
We present a class of non-standard numerical schemes which are modifications of the discrete gradient method. They preserve the energy integral exactly (up to the round-off error). The considered class contains locally exact discrete…
In this letter, we revisit the IEQ method and provide a new perspective on its ability to preserve the original energy dissipation laws. The invariant energy quadratization (IEQ) method has been widely used to design energy stable numerical…
We construct complete sets of invariant quantities that are integrals of motion for two Hamiltonian systems obtained through a reduction procedure, thus proving that these systems are maximally superintegrable. We also discuss the reduction…
In this series of eight papers we present the applications of methods from wavelet analysis to polynomial approximations for a number of accelerator physics problems. In this paper we consider invariant formulation of nonlinear (Lagrangian…
It is known that HLL-type schemes are more dissipative than schemes based on characteristic decompositions. However, HLL-type methods offer greater flexibility to large systems of hyperbolic conservation laws because the eigenstructure of…
This thesis develops advanced Tensor Network (TN) methods to address Hamiltonian Lattice Gauge Theories (LGTs), overcoming limitations in real-time dynamics and finite-density regimes. A novel dressed-site formalism is introduced, enabling…
In this study we work on a novel Hamiltonian system which is Liouville integrable. In the integrable Hamiltonian model, conserved currents can be represented as Binomial polynomials in which each order corresponds to the integral of motion…
Preservation of linear and quadratic invariants by numerical integrators has been well studied. However, many systems have linear or quadratic observables that are not invariant, but which satisfy evolution equations expressing important…