A splitting method for SDEs with locally Lipschitz drift : Illustration on the FitzHugh-Nagumo model
(2022) In Applied Numerical Mathematics 179. p.191-220- Abstract
In this article, we construct and analyse an explicit numerical splitting method for a class of semi-linear stochastic differential equations (SDEs) with additive noise, where the drift is allowed to grow polynomially and satisfies a global one-sided Lipschitz condition. The method is proved to be mean-square convergent of order 1 and to preserve important structural properties of the SDE. First, it is hypoelliptic in every iteration step. Second, it is geometrically ergodic and has an asymptotically bounded second moment. Third, it preserves oscillatory dynamics, such as amplitudes, frequencies and phases of oscillations, even for large time steps. Our results are illustrated on the stochastic FitzHugh-Nagumo model and compared with... (More)
In this article, we construct and analyse an explicit numerical splitting method for a class of semi-linear stochastic differential equations (SDEs) with additive noise, where the drift is allowed to grow polynomially and satisfies a global one-sided Lipschitz condition. The method is proved to be mean-square convergent of order 1 and to preserve important structural properties of the SDE. First, it is hypoelliptic in every iteration step. Second, it is geometrically ergodic and has an asymptotically bounded second moment. Third, it preserves oscillatory dynamics, such as amplitudes, frequencies and phases of oscillations, even for large time steps. Our results are illustrated on the stochastic FitzHugh-Nagumo model and compared with known mean-square convergent tamed/truncated variants of the Euler-Maruyama method. The capability of the proposed splitting method to preserve the aforementioned properties may make it applicable within different statistical inference procedures. In contrast, known Euler-Maruyama type methods commonly fail in preserving such properties, yielding ill-conditioned likelihood-based estimation tools or computationally infeasible simulation-based inference algorithms.
(Less)
- author
- Buckwar, Evelyn LU ; Samson, Adeline ; Tamborrino, Massimiliano and Tubikanec, Irene
- organization
- publishing date
- 2022
- type
- Contribution to journal
- publication status
- published
- subject
- keywords
- Ergodicity, FitzHugh-Nagumo model, Hypoellipticity, Locally Lipschitz drift, Mean-square convergence, Splitting methods, Stochastic differential equations
- in
- Applied Numerical Mathematics
- volume
- 179
- pages
- 30 pages
- publisher
- Elsevier
- external identifiers
-
- scopus:85129915562
- ISSN
- 0168-9274
- DOI
- 10.1016/j.apnum.2022.04.018
- language
- English
- LU publication?
- yes
- id
- d2e68e15-b8be-486c-8a2d-84f79cf77e1f
- date added to LUP
- 2022-07-14 13:17:56
- date last changed
- 2022-07-14 13:17:56
@article{d2e68e15-b8be-486c-8a2d-84f79cf77e1f, abstract = {{<p>In this article, we construct and analyse an explicit numerical splitting method for a class of semi-linear stochastic differential equations (SDEs) with additive noise, where the drift is allowed to grow polynomially and satisfies a global one-sided Lipschitz condition. The method is proved to be mean-square convergent of order 1 and to preserve important structural properties of the SDE. First, it is hypoelliptic in every iteration step. Second, it is geometrically ergodic and has an asymptotically bounded second moment. Third, it preserves oscillatory dynamics, such as amplitudes, frequencies and phases of oscillations, even for large time steps. Our results are illustrated on the stochastic FitzHugh-Nagumo model and compared with known mean-square convergent tamed/truncated variants of the Euler-Maruyama method. The capability of the proposed splitting method to preserve the aforementioned properties may make it applicable within different statistical inference procedures. In contrast, known Euler-Maruyama type methods commonly fail in preserving such properties, yielding ill-conditioned likelihood-based estimation tools or computationally infeasible simulation-based inference algorithms.</p>}}, author = {{Buckwar, Evelyn and Samson, Adeline and Tamborrino, Massimiliano and Tubikanec, Irene}}, issn = {{0168-9274}}, keywords = {{Ergodicity; FitzHugh-Nagumo model; Hypoellipticity; Locally Lipschitz drift; Mean-square convergence; Splitting methods; Stochastic differential equations}}, language = {{eng}}, pages = {{191--220}}, publisher = {{Elsevier}}, series = {{Applied Numerical Mathematics}}, title = {{A splitting method for SDEs with locally Lipschitz drift : Illustration on the FitzHugh-Nagumo model}}, url = {{http://dx.doi.org/10.1016/j.apnum.2022.04.018}}, doi = {{10.1016/j.apnum.2022.04.018}}, volume = {{179}}, year = {{2022}}, }