In this paper, a new numerical method based on adaptive gradient descent optimizers is provided for computing the implied volatility from the Black-Scholes (B-S) option pricing model. It is shown that the new method is more accurate than the close form approximation. Compared with the Newton-Raphson method, the new method obtains a reliable rate of convergence and tends to be less sensitive to the beginning point.
We derive a backward and forward nonlinear PDEs that govern the implied volatility of a contingent claim whenever the latter is well-defined. This would include at least any contingent claim written on a positive stock price whose payoff at a possibly random time is convex. We also discuss suitable initial and boundary conditions for those PDEs. Finally, we demonstrate how to solve them numerically by using an iterative finite-difference approach.
We present a strikingly simple proof that two rules are sufficient to automate gradient descent: 1) dont increase the stepsize too fast and 2) dont overstep the local curvature. No need for functional values, no line search, no information about the function except for the gradients. By following these rules, you get a method adaptive to the local geometry, with convergence guarantees depending only on the smoothness in a neighborhood of a solution. Given that the problem is convex, our method converges even if the global smoothness constant is infinity. As an illustration, it can minimize arbitrary continuously twice-differentiable convex function. We examine its performance on a range of convex and nonconvex problems, including logistic regression and matrix factorization.
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.
We study the shapes of the implied volatility when the underlying distribution has an atom at zero and analyse the impact of a mass at zero on at-the-money implied volatility and the overall level of the smile. We further show that the behaviour at small strikes is uniquely determined by the mass of the atom up to high asymptotic order, under mild assumptions on the remaining distribution on the positive real line. We investigate the structural difference with the no-mass-at-zero case, showing how one can--theoretically--distinguish between mass at the origin and a heavy-left-tailed distribution. We numerically test our model-free results in stochastic models with absorption at the boundary, such as the CEV process, and in jump-to-default models. Note that while Lees moment formula tells that implied variance is at most asymptotically linear in log-strike, other celebrated results for exact smile asymptotics such as Benaim and Friz (09) or Gulisashvili (10) do not apply in this setting--essentially due to the breakdown of Put-Call duality.
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.