Motivated by the challenges related to the calibration of financial models, we consider the problem of numerically solving a singular McKean–Vlasov equation \( d X_{t}= \sigma (t,X_{t}) X_{t} \frac{\sqrt{v}_{t}}{\sqrt{\mathbb{E}[v_{t}|X_{t}]}}dW_{t}, \) where $W$ is a Brownian motion and $v$ is an adapted diffusion process. This equation can be considered as a singular local stochastic volatility model. While such models are quite popular among practitioners, its well-posedness has unfortunately not yet been fully understood and in general is possibly not guaranteed at all. We develop a novel regularisation approach based on the reproducing kernel Hilbert space (RKHS) technique and show that the regularised model is well posed. Furthermore, we prove propagation of chaos. We demonstrate numerically that a thus regularised model is able to perfectly replicate option prices coming from typical local volatility models. Our results are also applicable to more general McKean–Vlasov equations.
