No Arabic abstract
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 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.
This note shows that the cosine expansion based on the Vieta formula is equivalent to a discretization of the Parseval identity. We then evaluate the use of simple direct algorithms to compute the Shannon coefficients for the payoff. Finally, we explore the efficiency of a Filon quadrature instead of the Vieta formula for the coefficients related to the probability density function.
Following Boukai (2021) we present the Generalized Gamma (GG) distribution as a possible RND for modeling European options prices under Hestons (1993) stochastic volatility (SV) model. This distribution is seen as especially useful in situations in which the spots price follows a negatively skewed distribution and hence, Black-Scholes based (i.e. the log-normal distribution) modeling is largely inapt. We apply the GG distribution as RND to modeling current market option data on three large market-index ETFs, namely the SPY, IWM and QQQ as well as on the TLT (an ETF that tracks an index of long term US Treasury bonds). The current option chain of each of the three market-index ETFs shows of a pronounced skew of their volatility `smile which indicates a likely distortion in the Black-Scholes modeling of such option data. Reflective of entirely different market expectations, this distortion appears not to exist in the TLT option data. We provide a thorough modeling of the available option data we have on each ETF (with the October 15, 2021 expiration) based on the GG distribution and compared it to the option pricing and RND modeling obtained directly from a well-calibrated Hestons (1993) SV model (both theoretically and empirically, using Monte-Carlo simulations of the spots price). All three market-index ETFs exhibit negatively skewed distributions which are well-matched with those derived under the GG distribution as RND. The inadequacy of the Black-Scholes modeling in such instances which involve negatively skewed distribution is further illustrated by its impact on the hedging factor, delta, and the immediate implications to the retail trader. In contrast, for the TLT ETF, which exhibits no such distortion to the volatility `smile, the three pricing models (i.e. Hestons, Black-Scholes and Generalized Gamma) appear to yield similar results.
This paper introduces an intermediary between conditional expectation and conditional sublinear expectation, called R-conditioning. The R-conditioning of a random-vector in $L^2$ is defined as the best $L^2$-estimate, given a $sigma$-subalgebra and a degree of model uncertainty. When the random vector represents the payoff of derivative security in a complete financial market, its R-conditioning with respect to the risk-neutral measure is interpreted as its risk-averse value. The optimization problem defining the optimization R-conditioning is shown to be well-posed. We show that the R-conditioning operators can be used to approximate a large class of sublinear expectations to arbitrary precision. We then introduce a novel numerical algorithm for computing the R-conditioning. This algorithm is shown to be strongly convergent. Implementations are used to compare the risk-averse value of a Vanilla option to its traditional risk-neutral value, within the Black-Scholes-Merton framework. Concrete connections to robust finance, sensitivity analysis, and high-dimensional estimation are all treated in this paper.
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.