A uniformly accurate scheme for the numerical integration of penalized Langevin dynamics

10/07/2021
by   Adrien Laurent, et al.
0

In molecular dynamics, penalized overdamped Langevin dynamics are used to model the motion of a set of particles that follow constraints up to a parameter ε. The most used schemes for simulating these dynamics are the Euler integrator in ℝ^d and the constrained Euler integrator. Both have weak order one of accuracy, but work properly only in specific regimes depending on the size of the parameter ε. We propose in this paper a new consistent method with an accuracy independent of ε for solving penalized dynamics on a manifold of any dimension. Moreover, this method converges to the constrained Euler scheme when ε goes to zero. The numerical experiments confirm the theoretical findings, in the context of weak convergence and for the invariant measure, on a torus and on the orthogonal group in high dimension and high codimension.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
06/21/2021

The weak convergence order of two Euler-type discretization schemes for the log-Heston model

We study the weak convergence order of two Euler-type discretizations of...
research
08/24/2023

Linear implicit approximations of invariant measures of semi-linear SDEs with non-globally Lipschitz coefficients

This article investigates the weak approximation towards the invariant m...
research
11/09/2021

Optimal-order error estimates for Euler discretization of high-index saddle dynamics

High-index saddle dynamics provides an effective means to compute the an...
research
09/21/2020

Curved Schemes for SDEs on Manifolds

Given a stochastic differential equation (SDE) in ℝ^n whose solution is ...
research
02/22/2022

Constrained multibody kinematics and dynamics in absolute coordinates: a discussion of three approaches to representing rigid body rotation

We compare three approaches to posing the index 3 set of differential al...
research
06/02/2021

Efficient computation of statistical properties of intermittent dynamics

Intermittent maps of the interval are simple and widely-studied models f...
research
03/20/2022

Analysis of a modified Euler scheme for parabolic semilinear stochastic PDEs

We propose a modification of the standard linear implicit Euler integrat...

Please sign up or login with your details

Forgot password? Click here to reset