Abstract:In this work, we develop a multifactor approximation for $d$-dimensional Stochastic Volterra Equations (SVE) with Lipschitz coefficients and kernels of completely monotone type that may be singular. First, we prove an $L^2$-estimation between two SVEs with different kernels, which provides a quantification of the error between the SVE and any multifactor Stochastic Differential Equation (SDE) approximation. For the particular rough kernel case with Hurst parameter lying in $(0,1/2)$, we propose various approximating multifactor kernels, state their rates of convergence and illustrate their efficiency for the rough Bergomi model. Second, we study a Euler discretization of the multifactor SDE and establish a convergence result towards the SVE that is uniform with respect to the approximating multifactor kernels. These obtained results lead us to build a new multifactor Euler scheme that reduces significantly the computational cost in an asymptotic way compared to the Euler scheme for SVEs. Finally, we show that our multifactor Euler scheme outperforms the Euler scheme for SVEs for option pricing in the rough Heston model.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is about how to effectively approximate Stochastic Volterra Equations (SVE), especially those with completely monotone kernels and possibly singular kernels. Specifically, the objectives of the paper can be summarized as follows:
1. **Error Estimation**: The author first proves the $L^2$ error estimate between SVEs with two different kernels, which provides a theoretical basis for quantifying the error between SVE and any multifactor Stochastic Differential Equation (SDE) approximation.
2. **Multifactor Approximation**: For the case of rough kernels with a small Hurst parameter ($H\in(0, 1/2)$), the author proposes several multifactor approximate kernels and gives their convergence rates. These approximation methods are verified for efficiency in the rough Bergomi model.
3. **Euler Discretization**: A Euler discretization method for multifactor SDE is studied, and the result of its convergence to SVE is established, which is uniform for the approximate kernels. This lays the foundation for constructing a new multifactor Euler scheme, which significantly reduces the computational cost.
4. **Computational Efficiency Improvement**: Through the combination of approximate kernels and Euler schemes, the effectiveness in calculating option prices in the rough Heston model is demonstrated. In particular, compared with the traditional SVE Euler scheme, the new scheme reduces the computing time by more than 5 times while maintaining the same accuracy.
### Main Contributions
- **Theoretical Contribution**: It provides a strong error estimate between SVE and multifactor SDE approximation, filling the gap in the existing literature.
- **Methodological Innovation**: A new multifactor Euler scheme is proposed, which significantly reduces the computational complexity and is especially suitable for high - dimensional SVE with completely monotone kernels.
- **Practical Application**: In financial applications, especially in option pricing, it shows the superior performance of the new method, especially when dealing with rough kernels with a small Hurst parameter.
### Formula Summary
- **General Form of SVE**:
\[
X_t=x_0+\int_{0}^{t}G_1(t - s)b(X_s)\,ds+\int_{0}^{t}G_2(t - s)\sigma(X_s)\,dW_s,\quad t\geq0
\]
where $x_0\in\mathbb{R}^d$, $b:\mathbb{R}^d\rightarrow\mathbb{R}^d$, $\sigma:\mathbb{R}^d\rightarrow M_d(\mathbb{R})$ are globally Lipschitz continuous coefficients, $W$ is a standard Brownian motion, and $G_1, G_2:\mathbb{R}_+^*\rightarrow M_d(\mathbb{R})$ are kernel functions.
- **Form of Kernel Functions**:
\[
G_j(t)=\int_{\mathbb{R}_+}e^{-\rho t}M_j(\rho)\lambda(d\rho),\quad t\in(0,+\infty)
\]
where $M_1, M_2:\mathbb{R}_+\rightarrow M_d(\mathbb{R})$ are bounded measurable functions, and $\lambda$ is a measure satisfying $\int_{\mathbb{R}_+}e^{-\rho t}\lambda(d\rho)<+\infty$.
- **Truncated Kernel Functions**:
\[
G_K^j(t)=\int_{[0, K)}e^{-\rho t}M_j(\rho)\lambda(d\rho),\quad t\geq0
\]
- **Error Estimation**:
\[
\mathbb{E}[|X_t-\hat{X}_t|^2]\leq C\left(r(