Saved in:
Bibliographic Details
Main Authors: Johnston, Tim, Monmarché, Pierre
Format: Preprint
Published: 2026
Subjects:
Online Access:https://arxiv.org/abs/2601.16878
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866909999112912896
author Johnston, Tim
Monmarché, Pierre
author_facet Johnston, Tim
Monmarché, Pierre
contents Many applications, such as systems of interacting particles in physics, require the simulation of diffusion processes with singular coefficients. Standard Euler schemes are then not convergent, and theoretical guarantees in this situation are scarce. In this work we introduce a Lyapunov-tamed Euler scheme, for drift coefficients for which the weak derivative is dominated by a function that obeys a certain generic Lyapunov-type condition. This allows for a range of coefficients that explode to infinity on a bounded set. We establish that, in terms of Lp-strong error, the Lyapunov-tamed scheme is consistent and moreover achieves the same order of convergence as the standard Euler scheme for Lipschitz coefficients. The general result is applied to systems of mean-field particles with singular repulsive interaction in 1D, yielding an error bound with polynomial dependency in the number of particles.
format Preprint
id arxiv_https___arxiv_org_abs_2601_16878
institution arXiv
publishDate 2026
record_format arxiv
spellingShingle A Lyapunov-tamed Euler method for singular SDEs
Johnston, Tim
Monmarché, Pierre
Probability
65C30, 65C35, 82C22, 60H35
Many applications, such as systems of interacting particles in physics, require the simulation of diffusion processes with singular coefficients. Standard Euler schemes are then not convergent, and theoretical guarantees in this situation are scarce. In this work we introduce a Lyapunov-tamed Euler scheme, for drift coefficients for which the weak derivative is dominated by a function that obeys a certain generic Lyapunov-type condition. This allows for a range of coefficients that explode to infinity on a bounded set. We establish that, in terms of Lp-strong error, the Lyapunov-tamed scheme is consistent and moreover achieves the same order of convergence as the standard Euler scheme for Lipschitz coefficients. The general result is applied to systems of mean-field particles with singular repulsive interaction in 1D, yielding an error bound with polynomial dependency in the number of particles.
title A Lyapunov-tamed Euler method for singular SDEs
topic Probability
65C30, 65C35, 82C22, 60H35
url https://arxiv.org/abs/2601.16878