How to simulate Lévy flights in a steep potential: An explicit splitting numerical scheme
arXiv:2508.07339 · doi:10.1063/5.0301147
Abstract
We propose an effective explicit numerical scheme for simulating solutions of stochastic differential equations with confining superlinear drift terms, driven by multiplicative heavy-tailed Lévy noise. The scheme is designed to prevent explosion and accurately capture all finite moments of the solutions. In the purely Gaussian case, it correctly reproduces moments of sub-Gaussian tails of the solutions. This method is particularly well-suited for approximating statistical moments and other probabilistic characteristics of Lévy flights in steep potential landscapes.
27 pages, 10 figures