Hierarchical adaptive sparse grids and quasi Monte Carlo for option pricing under the rough Bergomi model

—

Hierarchical adaptive sparse grids and quasi Monte Carlo for option pricing under the rough Bergomi model

Christian Bayer
Chiheb Ben Hammouda
Raul Tempone
Published on 12/20/2018
Options
Derivatives
Machine learning
Risk management
Volatility premium

This paper introduces novel deterministic quadrature methods for pricing European options under the rough Bergomi (rBergomi) model, a challenging non-Markovian stochastic volatility model. The authors propose two hierarchical approaches: adaptive sparse grids quadrature (ASGQ) and quasi Monte Carlo (QMC), both combined with a Brownian bridge construction and Richardson extrapolation. The Brownian bridge reduces the effective dimension of the integration problem, while Richardson extrapolation reduces the bias from time discretization, allowing coarser time grids. The methods exploit the smoothness of the integrand obtained via conditioning on the fractional Brownian motion path, which yields a Black-Scholes-type formula inside the expectation. Numerical experiments across various parameter sets, including very low Hurst parameters (H=0.05), demonstrate that ASGQ and QMC achieve substantial computational gains over standard Monte Carlo, often by one to two orders of magnitude, for relative error tolerances around 0.1%. The paper also provides numerical evidence that the weak error for both the hybrid and exact simulation schemes converges at order one, justifying the use of Richardson extrapolation. This work is the first to apply deterministic quadrature to rough volatility models, opening a new research direction beyond Monte Carlo for pricing and calibration.

Highlights

  • 1First application of deterministic quadrature methods (ASGQ and QMC) to option pricing under the rough Bergomi model.
  • 2Combines Brownian bridge construction and Richardson extrapolation to reduce effective dimension and bias.
  • 3Demonstrates substantial computational gains over standard Monte Carlo for small relative error tolerances, even for very low Hurst parameters.
  • 4Provides numerical evidence that the weak error for both hybrid and exact schemes is of order one.
  • 5Opens a new research direction for pricing under rough volatility models beyond Monte Carlo.

Methods

  • M
    Adaptive sparse grids quadrature (ASGQ) for high-dimensional integration.
  • M
    Quasi Monte Carlo (QMC) with Sobol sequences.
  • M
    Brownian bridge construction for hierarchical path generation.
  • M
    Richardson extrapolation to reduce bias from time discretization.

Results

  • R
    ASGQ and QMC achieve computational gains of up to two orders of magnitude over standard MC for relative error tolerances around 0.1%.
  • R
    The hybrid scheme exhibits consistent pre-asymptotic weak error behavior, making it suitable for Richardson extrapolation.
  • R
    The proposed methods are robust to the Hurst parameter, performing well even for H=0.05.
  • R
    Richardson extrapolation effectively reduces the number of time steps needed, lowering the effective dimension.
  • R
    The combination of Brownian bridge and deterministic quadrature significantly improves efficiency compared to standard MC.
0/5

Analyze Paper

Generate insights from "Hierarchical adaptive sparse grids and quasi Monte Carlo for...".

Suggested Actions