No Arabic abstract
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 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.
We consider the stochastic volatility model $dS_t = sigma_t S_t dW_t,dsigma_t = omega sigma_t dZ_t$, with $(W_t,Z_t)$ uncorrelated standard Brownian motions. This is a special case of the Hull-White and the $beta=1$ (log-normal) SABR model, which are widely used in financial practice. We study the properties of this model, discretized in time under several applications of the Euler-Maruyama scheme, and point out that the resulting model has certain properties which are different from those of the continuous time model. We study the asymptotics of the time-discretized model in the $nto infty$ limit of a very large number of time steps of size $tau$, at fixed $beta=frac12omega^2tau n^2$ and $rho=sigma_0^2tau$, and derive three results: i) almost sure limits, ii) fluctuation results, and iii) explicit expressions for growth rates (Lyapunov exponents) of the positive integer moments of $S_t$. Under the Euler-Maruyama discretization for $(S_t,log sigma_t)$, the Lyapunov exponents have a phase transition, which appears in numerical simulations of the model as a numerical explosion of the asset price moments. We derive criteria for the appearance of these explosions.
In this article, we show how the scaling symmetry of the SABR model can be utilized to efficiently price European options. For special kinds of payoffs, the complexity of the problem is reduced by one dimension. For more generic payoffs, instead of solving the 1+2 dimensional SABR PDE, it is sufficient to solve $N_V$ uncoupled 1+1 dimensional PDEs, where $N_V$ is the number of points used to discretize one dimension. Furthermore, the symmetry argument enables us to obtain prices of multiple options, whose payoffs are related to each other by convolutions, by valuing one of them. The results of the method are compared with the Monte Carlo simulation.
We propose a hybrid method for generating arbitrage-free implied volatility (IV) surfaces consistent with historical data by combining model-free Variational Autoencoders (VAEs) with continuous time stochastic differential equation (SDE) driven models. We focus on two classes of SDE models: regime switching models and Levy additive processes. By projecting historical surfaces onto the space of SDE model parameters, we obtain a distribution on the parameter subspace faithful to the data on which we then train a VAE. Arbitrage-free IV surfaces are then generated by sampling from the posterior distribution on the latent space, decoding to obtain SDE model parameters, and finally mapping those parameters to IV surfaces.
We propose a general, very fast method to quickly approximate the solution of a parabolic Partial Differential Equation (PDEs) with explicit formulas. Our method also provides equaly fast approximations of the derivatives of the solution, which is a challenge for many other methods. Our approach is based on a computable series expansion in terms of a small parameter. As an example, we treat in detail the important case of the SABR PDE for $beta = 1$, namely $partial_{tau}u = sigma^2 big [ frac{1}{2} (partial^2_xu - partial_xu) + u rho partial_xpartial_sigma u + frac{1}{2} u^2 partial^2_sigma u , big ] + kappa (theta - sigma) partial_sigma$, by choosing $ u$ as small parameter. This yields $u = u_0 + u u_1 + u^2 u_2 + ldots$, with $u_j$ independent of $ u$. The terms $u_j$ are explicitly computable, which is also a challenge for many other, related methods. Truncating this expansion leads to computable approximations of $u$ that are in closed form, and hence can be evaluated very quickly. Most of the other related methods use the time $tau$ as a small parameter. The advantage of our method is that it leads to shorter and hence easier to determine and to generalize formulas. We obtain also an explicit expansion for the implied volatility in the SABR model in terms of $ u$, similar to Hagans formula, but including also the {em mean reverting term.} We provide several numerical tests that show the performance of our method. In particular, we compare our formula to the one due to Hagan. Our results also behave well when used for actual market data and show the mean reverting property of the volatility.
Quasi-Gaussian HJM models are a popular approach for modeling the dynamics of the yield curve. This is due to their low dimensional Markovian representation, which greatly simplifies their numerical implementation. We present a qualitative study of the solutions of the quasi-Gaussian log-normal HJM model. Using a small-noise deterministic limit we show that the short rate may explode to infinity in finite time. This implies the explosion of the Eurodollar futures prices in this model. We derive explicit explosion criteria under mild assumptions on the shape of the yield curve.