Peter K. Friz,William Salkeld,Thomas Wagenhofer
Abstract:We consider a class of stochastic processes with rough stochastic volatility, examples of which include the rough Bergomi and rough Stein-Stein model, that have gained considerable importance in quantitative finance.
A basic question for such (non-Markovian) models concerns efficient numerical schemes. While strong rates are well understood (order $H$), we tackle here the intricate question of weak rates. Our main result asserts that the weak rate, for a reasonably large class of test function, is essentially of order $\min \{ 3H+\tfrac12, 1 \}$ where $H \in (0,1/2]$ is the Hurst parameter of the fractional Brownian motion that underlies the rough volatility process.
Interestingly, the phase transition at $H=1/6$ is related to the correlation between the two driving factors, and thus gives additional meaning to a quantity already of central importance in stochastic volatility modelling.Our results are complemented by a lower bound which show that the obtained weak rate is indeed optimal.
What problem does this paper attempt to address?
This paper aims to study the problem of weak convergence rates in stochastic processes with rough stochastic volatility. Specifically, the paper focuses on the weak error estimation of stochastic integral functions driven by fractional Brownian motion (fBm). Such models are very important in financial quantitative analysis, such as the rough Bergomi model and the rough Stein - Stein model.
### Research Background and Motivation
Traditional numerical methods for stochastic differential equations usually focus on strong convergence rates, that is, the rates of path convergence. However, in practical applications, such as option pricing, the weak convergence rate (that is, the convergence rate of the expected value) is more important. For systems driven by fractional Brownian motion, the known strong convergence rate is \(H\), but the study of the weak convergence rate is more complex, especially in the "rough" region where \(H\in(0, 1/2]\).
### Main Contributions
The main result of the paper is to prove that, under a reasonably wide class of test functions, the weak convergence rate is essentially \(\min\{3H+\frac{1}{2}, 1\}\), where \(H\) is the Hurst parameter of fractional Brownian motion. The importance of this result lies in:
1. **Phase - transition point**: When \(H = \frac{1}{6}\), the weak convergence rate undergoes a phase transition, which is related to the correlation between the two driving factors, thus giving the Hurst parameter additional significance in stochastic volatility modeling.
2. **Optimality**: The paper also provides a lower - bound estimate, indicating that the obtained weak convergence rate is optimal.
### Technical Details
To obtain these results, the authors used an exact expansion method, based on an iterative scheme, and combined with Ito and Malliavin calculus. Key steps include:
- **Exact expansion**: For polynomial test functions, an exact expansion is carried out through an iterative scheme.
- **Kernel function modification**: The kernel function of the left - point approximation method is appropriately modified.
- **Term - by - term comparison**: Through term - by - term comparison of complex iterative integrals, the weak convergence rate is finally obtained.
### Special Cases
- When the correlation coefficient \(\rho = 0\), the weak convergence rate is increased to 1.
- For the linear wave function \(f(x)=x\) and polynomial test functions, the weak convergence rate can be further optimized.
### Conclusion
Through strict mathematical derivations, the paper provides an explicit weak convergence rate estimate for models with rough stochastic volatility. This result is not only of great significance for financial quantitative analysis, but also provides a theoretical basis for the weak convergence analysis of multi - scale stochastic systems in other fields.