Localized radial basis functions for no-arbitrage pricing of options under stochastic alpha–beta–rho dynamics

Closed-form explicit formulas for implied Black–Scholes volatilities provide a rapid evaluation method for European options under the popular stochastic alpha–beta–rho (SABR) model. However, it is well known that computed prices using the implied volatilities are only accurate for short-term maturities, but, for longer maturities, a more accurate method is required. This work addresses this accuracy problem for long-term maturities by numerically solving the no-arbitrage partial differential equation with an absorbing boundary condition at zero. Localized radial basis functions in a finite-difference mode are employed for the development of a computational method for solving the resulting two-dimensional pricing equation. The proposed method can use either multiquadrics or inverse multiquadrics, which are shown to have comparable performances. Numerical results illustrate the accuracy of the proposed method and, more importantly, that the computed risk-neutral probability densities are nonnegative. These two key properties indicate that the method of solution using localized meshless methods is a viable and efficient means for price computations under SABR dynamics. doi:10.1017/S1446181121000237

LOCALIZED RADIAL BASIS FUNCTIONS FOR NO-ARBITRAGE PRICING OF OPTIONS UNDER STOCHASTIC ALPHA–BETA–RHO DYNAMICS

Closed-form explicit formulas for implied Black–Scholes volatilities provide a rapid evaluation method for European options under the popular stochastic alpha–beta–rho (SABR) model. However, it is well known that computed prices using the implied volatilities are only accurate for short-term maturities, but, for longer maturities, a more accurate method is required. This work addresses this accuracy problem for long-term maturities by numerically solving the no-arbitrage partial differential equation with an absorbing boundary condition at zero. Localized radial basis functions in a finite-difference mode are employed for the development of a computational method for solving the resulting two-dimensional pricing equation. The proposed method can use either multiquadrics or inverse multiquadrics, which are shown to have comparable performances. Numerical results illustrate the accuracy of the proposed method and, more importantly, that the computed risk-neutral probability densities are nonnegative. These two key properties indicate that the method of solution using localized meshless methods is a viable and efficient means for price computations under SABR dynamics.

A rough SABR formula

<p style='text-indent:20px;'>Following an approach originally suggested by Balland in the context of the SABR model, we derive an ODE that is satisfied by normalized volatility smiles for short maturities under a rough volatility extension of the SABR model that extends also the rough Bergomi model. We solve this ODE numerically and further present a very accurate approximation to the numerical solution that we dub the <i>rough SABR formula</i>.</p>

ASYMPTOTICS OF THE TIME-DISCRETIZED LOG-NORMAL SABR MODEL: THE IMPLIED VOLATILITY SURFACE

We propose a novel time discretization for the log-normal SABR model which is a popular stochastic volatility model that is widely used in financial practice. Our time discretization is a variant of the Euler–Maruyama scheme. We study its asymptotic properties in the limit of a large number of time steps under a certain asymptotic regime which includes the case of finite maturity, small vol-of-vol and large initial volatility with fixed product of vol-of-vol and initial volatility. We derive an almost sure limit and a large deviations result for the log-asset price in the limit of a large number of time steps. We derive an exact representation of the implied volatility surface for arbitrary maturity and strike in this regime. Using this representation, we obtain analytical expansions of the implied volatility for small maturity and extreme strikes, which reproduce at leading order known asymptotic results for the continuous time model.