Related papers: A New Time-Reversible Integrator for Molecular Dyn…
Molecular dynamics (MD) provides insights into atomic-scale processes by integrating over time the equations that describe the motion of atoms under the action of interatomic forces. Machine learning models have substantially accelerated MD…
This work presents a novel algorithm for impulsive optimal control of linear time-varying systems with the inclusion of input magnitude constraints. Impulsive optimal control problems, where the optimal input solution is a sum of delta…
A popular version of the finite strain Maxwell fluid is considered, which is based on the multiplicative decomposition of the deformation gradient tensor. The model combines Newtonian viscosity with hyperelasticity of Mooney-Rivlin type; it…
An algorithm for a family of self-starting high-order implicit time integration schemes with controllable numerical dissipation is proposed for both linear and nonlinear transient problems. This work builds on the previous works of the…
An inversion method for time-resolved data from ultrafast experiments is introduced, based on forward-optimisation in a trajectory basis. The method is applied to experimental data from x-ray scattering of the photochemical ring-opening…
A program for molecular calculations with B functions is reported and its performance is analyzed. All the one- and two-center integrals, and the three-center nuclear attraction integrals are computed by direct procedures, using previously…
Various many-body models are treated, which describe $N$ points confined to move on a plane circle. Their Newtonian equations of motion ("accelerations equal forces") are integrable, i. e. they allow the explicit exhibition of $N$ constants…
We propose a new route to accelerate molecular dynamics through the use of velocity jump processes allowing for an adaptive time-step specific to each atom-atom pair (2-body) interactions. We start by introducing the formalism of the new…
This paper studies multistep methods for the integration of reversible dynamical systems, with particular emphasis on the planar Kepler problem. It has previously been shown by Cano & Sanz-Serna that reversible linear multisteps for…
A general procedure for constructing conservative numerical integrators for time dependent partial differential equations is presented. In particular, linearly implicit methods preserving a time discretised version of the invariant is…
In this work we illustrate the basic development of the constrained molecular dynamics applied to the N-body problem in nuclear physics. The heavy computational taskes related to quantum effects, to the presence of the "hard core" repulsive…
Two families of symplectic methods specially designed for second-order time-dependent linear systems are presented. Both are obtained from the Magnus expansion of the corresponding first-order equation, but otherwise they differ in…
In this and subsequent paper arXiv:1011.5185 we develop a recursive approach for calculating the short-time expansion of the propagator for a general quantum system in a time-dependent potential to orders that have not yet been accessible…
Simulating the dynamics of ions near polarizable nanoparticles (NPs) using coarse-grained models is extremely challenging due to the need to solve the Poisson equation at every simulation timestep. Recently, a molecular dynamics (MD) method…
We introduce a numerical method to integrate tidal effects on collisional systems, using any definition of the external potential as a function of space and time. Rather than using a linearisation of the tidal field, this new method follows…
This paper is an attempt to solve an important class of hypersingular integral equations of the second kind. To this end, we apply a new weighted and modified perturbation method which includes some special cases of the Adomian…
This paper describes a fourth-order integration algorithm for the gravitational N-body problem based on discrete Lagrangian mechanics. When used with shared timesteps, the algorithm is momentum conserving and symplectic. We generalize the…
Time integration of ODEs or time-dependent PDEs with required resolution of the fastest time scales of the system, can be very costly if the system exhibits multiple time scales of different magnitudes. If the different time scales are…
A novel numerical approach to solving the shallow-water equations on the sphere using high-order numerical discretizations in both space and time is proposed. A space-time tensor formalism is used to express the equations of motion…
In this work, an approximate family of implicit multiderivative Runge-Kutta (MDRK) time integrators for stiff initial value problems is presented. The approximation procedure is based on the recent Approximate Implicit Taylor method (Baeza…