A Tail-Respecting Explicit Numerical Scheme for L\'evy-Driven SDEs With Superlinear Drifts

We present an explicit numerical approximation scheme, denoted by $\{X^n\}$, for the effective simulation of solutions $X$ to a multivariate stochastic differential equation (SDE) with a superlinearly growing $\kappa$-dissipative drift, where $\kappa>1$, driven by a multiplicative heavy-tailed L\'evy process that has a finite $p$-th moment, with $p>0$. We show that the strong $L^q$-convergence $\sup_{t\in[0,T]}\mathbf E \|X^n_t-X_t\|^q=\mathcal O (h_n^{\gamma})$ holds true for any $q\in (0,p+\kappa-1)$, which is exactly the range where the $q$-moment of the solution is known to be finite. Additionally, for any $q\in (0,p)$ we establish strong uniform convergence: $\mathbf E\sup_{t\in[0,T]} \|X^n_t-X_t\|^q=\mathcal{O} ( h_n^{\delta} )$. In both cases we determine the convergence rates $\gamma$ and $\delta$. In the special case of SDEs driven solely by a Brownian motion, our numerical scheme preserves super-exponential moments of the solution. The scheme $\{X^n\}$ is realized as a combination of a well-known Euler method with a Lie-Trotter type splitting technique.