Least squares Monte Carlo methods in stochastic Volterra rough volatility models

In stochastic Volterra rough volatility models, the volatility follows a truncated Brownian semi-stationary process with stochastic vol-of-vol.Recently, ecient VIX pricing Monte Carlo methods have been proposed for the case where the vol-of-vol is Markovian and independent of the volatility. Following recent empirical data, we discuss the VIX option pricing problem for a generalized framework of these models, where the vol-of-vol may depend on the volatility and/or not be Markovian. In such a setting, the aforementioned Monte Carlo methods are not valid.Moreover, the classical least squares Monte Carlo faces exponentially increasing complexity with the number of grid time steps, whilst the nested Monte Carlo method requires a prohibitive number of simulations. By exploring the innite dimensional Markovian representation of these models, we device a scalable least squares Monte Carlo for VIX option pricing. We apply our method rstly under the independence assumption for benchmarks, and then to the generalized framework. We also discuss the rough vol-of-vol setting, where Markovianity of the vol-of-vol is not present. We present simulations and benchmarks to establish the eciency of our method.