Pathwise skew-symmetric discretisation for SDEs with superlinear drift

The skew-symmetric discretisation has recently been proposed as a new robust simulation method for weakly approximating stochastic differential equations (SDEs) with non-globally Lipschitz drift. This work develops a pathwise version of the scheme by representing the noise increment as a skew-normal distribution and coupling it with the driving Brownian increments, thereby enabling its use in the multilevel Monte Carlo (MLMC) framework. Under suitable conditions, we establish strong convergence of order 1/2 in $L^2$. Subsequently, the associated MLMC estimator is shown to have computational complexity $\mathcal{O} \bigl(\varepsilon^{-2} (\log (1/\varepsilon))^2 \bigr)$ to achieve a mean-squared error $\varepsilon^2$. We then analytically compare the proposed scheme with the tamed Euler scheme, another benchmark for robust discretisation. Under a strong inward-drift regime with the current state being far from the stable region, we show that the probability of moving in the wrong direction tends to vanish in the skew-symmetric scheme, whereas the tamed Euler scheme makes such moves with a non-trivial probability. Furthermore, in the MLMC setting, employing a one-dimensional stochastic Ginzburg-Landau model, we specify the range of step sizes for which the asymptotic variance of the coupled level difference obtained via the pathwise skew-symmetric scheme is lower than that obtained via the tamed Euler scheme. Numerical experiments on several model examples support the theoretical rate of strong convergence and demonstrate the stability and effectiveness of the resulting MLMC in the superlinear drift setting.