Related papers: A variational symplectic scheme based on Simpson's…
We have applied a collocation approach to obtain the numerical solution to the stationary Schr\"odinger equation for systems of coupled oscillators. The dependence of the discretized Hamiltonian on scale and angle parameters is exploited to…
We revisit an algorithm by Skeel et al. for computing the modified, or shadow, energy associated with the symplectic discretization of Hamiltonian systems. By rephrasing the algorithm as a Richardson extrapolation scheme arbitrary high…
A powerful method for calculating the eigenvalues of a Hamiltonian operator consists of converting the energy eigenvalue equation into a matrix equation by means of an appropriate basis set of functions. The convergence of the method can be…
A new symplectic time-reversible algorithm for numerical integration of the equations of motion in magnetic liquids is proposed. It is tested and applied to molecular dynamics simulations of a Heisenberg spin fluid. We show that the…
We construct several variational integrators--integrators based on a discrete variational principle--for systems with Lagrangians of the form L = L_A + epsilon L_B, with epsilon << 1, where L_A describes an integrable system. These…
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…
We derive the dynamics of several rigid bodies of arbitrary shape in a 2-dimensional inviscid and incompressible fluid, whose vorticity field is given by point vortices. We adopt the idea of Vankerschaver et al. (2009) to derive the…
Recently a new class of numerical integration methods -- ``mixed variable symplectic integrators'' -- has been introduced for studying long-term evolution in the conservative gravitational few-body problem. These integrators are an order of…
We present a brief tutorial on the nuts and bolts computation of a multisymplectic particle-in-cell algorithm using the discretized Lagrangian approach. This approach, originated by Marsden, Shadwick, and others, brings the benefits of…
We present a detailed comparison of several integration schemes applied to the dynamic system consisting of a charged particle on the Kerr background endowed with the axisymmetric electromagnetic test field. In particular, we compare the…
A class of Hamiltonian stochastic differential equations with multiplicative L\'{e}vy noise in the sense of Marcus, and the construction and numerical implementation methods of symplectic Euler scheme, are considered. A general symplectic…
In this paper, based on the weak form of the Hamiltonian formulation of the regularized long-wave equation and a novel approach of transforming the original Hamiltonian energy into a quadratic functional, a fully implicit and three…
There is a growing interest in the conservation of invariants when numerically solving a system of ordinary differential equations. Methods that exactly preserve these quantities in time are known as geometric integrators. In this paper we…
We describe an algorithm that computes the ground state energy and correlation functions for 2-local Hamiltonians in which interactions between qubits are weak compared to single-qubit terms. The running time of the algorithm is polynomial…
An appropriate iterative scheme for the minimization of the energy, based on the variational Monte Carlo (VMC) technique, is introduced and compared with existing stochastic schemes. We test the various methods for the 1D Heisenberg ring…
In this paper we study the performance of a symplectic numerical integrator based on the splitting method. This method is applied to a subtle problem i.e. higher order resonance of the elastic pendulum. In order to numerically study the…
For the calculation of the partition function $\mathcal{Z}$ of small, isolated and interacting many body systems an improvement with respect to previous formulations is presented. By including anharmonicities and employing a variational…
Multiscale dynamics are ubiquitous in applications of modern science. Because of time scale separation between relatively small set of slowly evolving variables and (typically) much larger set of rapidly changing variables, direct numerical…
In this work, we introduce a simple modification of the Monte Carlo algorithm, which we call step Monte Carlo (sMC). The sMC approach allows to simulate processes far from equilibrium and obtain information about the dynamic properties of…
We consider the numerical approximations of the Cahn-Hilliard equation with dynamic boundary conditions (C. Liu et. al., Arch. Rational Mech. Anal., 2019). We propose a first-order in time, linear and energy stable numerical scheme, which…