Related papers: Extended phase-space symplectic-like integrators f…
Extended Lagrangian Born-Oppenheimer molecular dynamics [Niklasson, Phys. Rev. Lett. 100 123004 (2008)] has been generalized to the propagation of the electronic wavefunctions. The technique allows highly efficient first principles…
We prove that the recently developed semiexplicit symplectic integrators for non-separable Hamiltonian systems preserve any linear and quadratic invariants possessed by the Hamiltonian systems. This is in addition to being symmetric and…
In this paper, we consider the task of efficiently computing the numerical solution of evolutionary complex Ginzburg--Landau equations on Cartesian product domains with homogeneous Dirichlet/Neumann or periodic boundary conditions. To this…
The contour integrals, occurring in the arbitrary-order phase-integral quantization conditions given in a previous paper, are in the first- and third-order approximations expressed in terms of complete elliptic integrals in the case that…
We discuss some Lagrangian and presymplectic models concerning test particles in electromagnetic and gravitational fields, with the aim of describing an upper bound to the acceleration. Some models are based on the relativistic phase space…
We applied the method of finite-part integration [Galapon E.A Proc.R.Soc A 473, 20160567(2017)] to evaluate in closed-form the exact one-loop integral representations of the Heisenberg-Euler Lagrangian from QED for a constant magnetic field…
Today, the motion of spacecrafts is still described according to the classical Newtonian equations plus the so-called "relativistic corrections", computed with the required precision using the Post-(Post-)Newtonian formalism. The current…
Direct N-body simulations and symplectic integrators are effective tools to study the long-term evolution of planetary systems. The Wisdom-Holman (WH) integrator in particular has been used extensively in planetary dynamics as it allows for…
We introduce a new family of symplectic integrators depending on a real parameter. When the paramer is zero, the corresponding method in the family becomes the classical Gauss collocation formula of order 2s, where s denotes the number of…
The evolution of any factorized time-reversible symplectic integrators, when applied to the harmonic oscillator, can be exactly solved in a closed form. The resulting modified Hamiltonians demonstrate the convergence of the Lie series…
We present the Euler--Langrage equations for a many-body system of coupled planar pendulums. Hence, imposing initial condition data, the equations of motion are linearized and later developed in an idealized model for the pseudo-periodicity…
In this paper, exponential Runge-Kutta methods of collocation type (ERKC) which were originally proposed in (Appl Numer Math 53:323-339, 2005) are extended to semilinear parabolic problems with time-dependent delay. Two classes of the ERKC…
In the present work, we obtain the constants of motion for isoperimetric variational problems with time delay. We consider a constrained optimization problem where the Lagrangian function defining the functional depends on time delayed…
We present a new symplectic integrator designed for collisional gravitational $N$-body problems which makes use of Kepler solvers. The integrator is also reversible and conserves 9 integrals of motion of the $N$-body problem to machine…
A novel semi-Lagrangian method is introduced to solve numerically the Euler equation for ideal incompressible flow in arbitrary space dimension. It exploits the time-analyticity of fluid particle trajectories and requires, in principle,…
We investigate the computational performance of various numerical methods for the integration of the equations of motion and the variational equations for some typical classical many-body models of condensed matter physics: the…
We characterize non-degenerate Lagrangians of the form $ \int f(u_x, u_y, u_t) dx dy dt $ such that the corresponding Euler-Lagrange equations $ (f_{u_x})_x+ (f_{u_y})_y+ (f_{u_t})_t=0 $ are integrable by the method of hydrodynamic…
We consider the numerical solution of the inertial version of Landau-Lifshitz-Gilbert equation (iLLG), which describes high-frequency nutation on top of magnetization precession due to angular momentum relaxation. The iLLG equation defines…
In general, the system of $2$nd-order partial differential equations made of the Euler-Lagrange equations of classical field theories are not compatible for singular Lagrangians. This is the so-called second-order problem. The first aim of…
We construct a symplectic integrator for non-separable Hamiltonian systems combining an extended phase space approach of Pihajoki and the symmetric projection method. The resulting method is semiexplicit in the sense that the main time…