Equity-linked annuities with multiscale hybrid stochastic and local volatility

In recent times, hybrid underlying models have become an industry standard for the pricing of derivatives and other problems in finance. This paper chooses a hybrid stochastic and local volatility model to evaluate an equity-linked annuity (ELA), which is a sort of tax-deferred annuity whose credited interest is linked to an equity index. The stochastic volatility component of the hybrid model is driven by a fast mean-reverting diffusion process while the local volatility component is given by the constant elasticity of variance (CEV) model. Since contracts of the ELA usually have long maturities over 10 years, a slowly moving factor in the stochastic volatility of stock index is expected to play a significant role in the valuation of the ELA, and thus, it is added to the aforementioned model. Based on this multiscale hybrid model, an analytic approximate formula is obtained for the price of a European option in terms of the CEV probability density function and then the result is applied to the value of the point-to-point ELA. The formula leads to the dependence structure of the ELA price on the fast and slow scale stochastic volatility and the elasticity of variance.     

An Asymptotic Expansion Approach to Currency Options with a Market Model of Interest Rates under Stochastic Volatility Processes of Spot Exchange Rates

This paper proposes an asymptotic expansion scheme of currency options with a libor market model of interest rates and stochastic volatility models of spot exchange rates. In particular, we derive closed-form approximation formulas for the density functions of the underlying assets and for pricing currency options based on a third order asymptotic expansion scheme; we do not model a foreign exchange rate’s variance such as in Heston [(1993) The Review of Financial studies, 6, 327–343], but its volatility that follows a general time-inhomogeneous Markovian process. Further, the correlations among all the factors such as domestic and foreign interest rates, a spot foreign exchange rate and its volatility, are allowed. Finally, numerical examples are provided and the pricing formula are applied to the calibration of volatility surfaces in the JPY/USD option market.     

A new approach for option pricing under stochastic volatility

We develop a new approach for pricing European-style contingent claims written on the time T spot price of an underlying asset whose volatility is stochastic. Like most of the stochastic volatility literature, we assume continuous dynamics for the price of the underlying asset. In contrast to most of the stochastic volatility literature, we do not directly model the dynamics of the instantaneous volatility. Instead, taking advantage of the recent rise of the variance swap market, we directly assume continuous dynamics for the time T variance swap rate. The initial value of this variance swap rate can either be directly observed, or inferred from option prices. We make no assumption concerning the real world drift of this process. We assume that the ratio of the volatility of the variance swap rate to the instantaneous volatility of the underlying asset just depends on the variance swap rate and on the variance swap maturity. Since this ratio is assumed to be independent of calendar time, we term this key assumption the stationary volatility ratio hypothesis (SVRH). The instantaneous volatility of the futures follows an unspecified stochastic process, so both the underlying futures price and the variance swap rate have unspecified stochastic volatility. Despite this, we show that the payoff to a path-independent contingent claim can be perfectly replicated by dynamic trading in futures contracts and variance swaps of the same maturity. As a result, the contingent claim is uniquely valued relative to its underlying’s futures price and the assumed observable variance swap rate. In contrast to standard models of stochastic volatility, our approach does not require specifying the market price of volatility risk or observing the initial level of instantaneous volatility. As a consequence of our SVRH, the partial differential equation (PDE) governing the arbitrage-free value of the contingent claim just depends on two state variables rather than the usual three. We then focus on the consistency of our SVRH with the standard assumption that the risk-neutral process for the instantaneous variance is a diffusion whose coefficients are independent of the variance swap maturity. We show that the combination of this maturity independent diffusion hypothesis (MIDH) and our SVRH implies a very special form of the risk-neutral diffusion process for the instantaneous variance. Fortunately, this process is tractable, well-behaved, and enjoys empirical support. Finally, we show that our model can also be used to robustly price and hedge volatility derivatives.     

Pricing via recursive quantization in stochastic volatility models

We provide the first recursive quantization-based approach for pricing options in the presence of stochastic volatility. This method can be applied to any model for which an Euler scheme is available for the underlying price process and it allows one to price vanillas, as well as exotics, thanks to the knowledge of the transition probabilities for the discretized stock process. We apply the methodology to some celebrated stochastic volatility models, including the Stein and Stein [Rev. Financ. Stud. 1991, (4), 727–752] model and the SABR model introduced in Hagan et al. [Wilmott Mag., 2002, 84–108]. A numerical exercise shows that the pricing of vanillas turns out to be accurate; in addition, when applied to some exotics like equity-volatility options, the quantization-based method overperforms by far the Monte Carlo simulation.     

General approximation schemes for option prices in stochastic volatility models

In this paper we develop a general method for deriving closed-form approximations of European option prices and equivalent implied volatilities in stochastic volatility models. Our method relies on perturbations of the model dynamics and we show how the expansion terms can be calculated using purely probabilistic methods. A flexible way of approximating the equivalent implied volatility from the basic price expansion is also introduced. As an application of our method we derive closed-form approximations for call prices and implied volatilities in the Heston [Rev. Financial Stud., 1993, 6, 327–343] model. The accuracy of these approximations is studied and compared with numerically obtained values.     

Maturity cycles in implied volatility

The skew effect in market implied volatility can be reproduced by option pricing theory based on stochastic volatility models for the price of the underlying asset. Here we study the performance of the calibration of the S&P 500 implied volatility surface using the asymptotic pricing theory under fast mean-reverting stochastic volatility described in [8]. The time-variation of the fitted skew-slope parameter shows a periodic behaviour that depends on the option maturity dates in the future, which are known in advance. By extending the mathematical analysis to incorporate model parameters which are time-varying, we show this behaviour can be explained in a manner consistent with a large model class for the underlying price dynamics with time-periodic volatility coefficients.Received: December 2003, Mathematics Subject Classification (2000): 91B70, 60F05, 60H30JEL Classification: C13, G13Jean-Pierre Fouque: Work partially supported by NSF grant DMS-0071744.Ronnie Sircar: Work supported by NSF grant DMS-0090067. We are grateful to Peter Thurston for research assistance.We thank a referee for his/her comments which improved the paper.     

Efficient estimation of drift parameters in stochastic volatility models

We study the parametric problem of estimating the drift coefficient in a stochastic volatility model , where Y is a log price process and V the volatility process. Assuming that one can recover the volatility, precisely enough, from the observation of the price process, we construct an efficient estimator for the drift parameter of the diffusion V. As an application we present the efficient estimation based on the discrete sampling with δ _n →0 and n δ _n →∞. We show that our setup is general enough to cover the case of ‘microstructure noise’ for the price process as well.      

Generic pricing of FX,inflation and stock options under stochastic interest rates and stochastic volatility

We consider the pricing of FX, inflation and stock options under stochastic interest rates and stochastic volatility, for which we use a generic multi-currency framework. We allow for a general correlation structure between the drivers of the volatility, the inflation index, the domestic (nominal) and the foreign (real) rates. Having the flexibility to correlate the underlying FX/inflation/stock index with both stochastic volatility and stochastic interest rates yields a realistic model that is of practical importance for the pricing and hedging of options with a long-term exposure. We derive explicit valuation formulas for various securities, such as vanilla call/put options, forward starting options, inflation-indexed swaps and inflation caps/floors. These vanilla derivatives can be valued in closed form under Schöbel and Zhu [Eur. Finance Rev., 1999, 4, 23–46] stochastic volatility, whereas we devise an (Monte Carlo) approximation in the form of a very effective control variate for the general Heston [Rev. Financial Stud., 1993, 6, 327–343] model. Finally, we investigate the quality of this approximation numerically and consider a calibration example to FX and inflation market data.     

A new closed-form solution as an extension of the Black–Scholes formula allowing smile curve plotting

In this paper, as a generalization of the Black–Scholes (BS) model, we elaborate a new closed-form solution for a uni-dimensional European option pricing model called the J-model. This closed-form solution is based on a new stochastic process, called the J-process, which is an extension of the Wiener process satisfying the martingale property. The J-process is based on a new statistical law called the J-law, which is an extension of the normal law. The J-law relies on four parameters in its general form. It has interesting asymmetry and tail properties, allowing it to fit the reality of financial markets with good accuracy, which is not the case for the normal law. Despite the use of one state variable, we find results similar to those of Heston dealing with the bi-dimensional stochastic volatility problem for pricing European calls. Inverting the BS formula, we plot the smile curve related to this closed-form solution. The J-model can also serve to determine the implied volatility by inverting the J-formula and can be used to price other kinds of options such as American options.     

Lévy processes driven by stochastic volatility

In this paper we extend option pricing under Lévy dynamics, by assuming that the volatility of the Lévy process is stochastic. We, therefore, develop the analog of the standard stochastic volatility models, when the underlying process is not a standard (unit variance) Brownian motion, but rather a standardized Lévy process. We present a methodology that allows one to compute option prices, under virtually any set of diffusive dynamics for the parameters of the volatility process. First, we use ‘local consistency’ arguments to approximate the volatility process with a finite, but sufficiently dense Markov chain; we then use this regime switching approximation to efficiently compute option prices using Fourier inversion. A detailed example, based on a generalization of the popular stochastic volatility model of Heston (Rev Financial Stud 6 (1993) 327), is used to illustrate the implementation of the algorithms. Computer code is available at www.theponytail.net/     

Turbocharging Monte Carlo pricing for the rough Bergomi model

The rough Bergomi model, introduced by Bayer et al. [Quant. Finance, 2016, 16(6), 887–904], is one of the recent rough volatility models that are consistent with the stylised fact of implied volatility surfaces being essentially time-invariant, and are able to capture the term structure of skew observed in equity markets. In the absence of analytical European option pricing methods for the model, we focus on reducing the runtime-adjusted variance of Monte Carlo implied volatilities, thereby contributing to the model’s calibration by simulation. We employ a novel composition of variance reduction methods, immediately applicable to any conditionally log-normal stochastic volatility model. Assuming one targets implied volatility estimates with a given degree of confidence, thus calibration RMSE, the results we demonstrate equate to significant runtime reductions—roughly 20 times on average, across different correlation regimes.     

New solvable stochastic volatility models for pricing volatility derivatives

In this paper we discuss a new approach to extend a class of solvable stochastic volatility models (SVM). Usually, classical SVM adopt a CEV process for instantaneous variance where the CEV parameter γ takes just few values: 0—the Ornstein–Uhlenbeck process, 1/2—the Heston (or square root) process, 1—GARCH, and 3/2—the 3/2 model. Some other models, e.g. with γ = 2 were discovered in Henry-Labordére (Analysis, geometry, and modeling in finance: advanced methods in option pricing. Chapman & Hall/CRC Financial Mathematics Series, London, 2009) by making connection between stochastic volatility and solvable diffusion processes in quantum mechanics. In particular, he used to build a bridge between solvable superpotentials (the Natanzon superpotentials, which allow reduction of a Schrödinger equation to a Gauss confluent hypergeometric equation) and existing SVM. Here we propose some new models with ${\gamma \in \mathbb{R}}$ and demonstrate that using Lie’s symmetries they could be priced in closed form in terms of hypergeometric functions. Thus obtained new models could be useful for pricing volatility derivatives (variance and volatility swaps, moment swaps).     

Options on realized variance by transform methods: a non-affine stochastic volatility model

In this paper we study the pricing and hedging of options on realized variance in the 3/2 non-affine stochastic volatility model by developing efficient transform-based pricing methods. This non-affine model gives prices of options on realized variance that allow upward-sloping implied volatility of variance smiles. Heston's model [Rev. Financial Stud., 1993, 6, 327–343], the benchmark affine stochastic volatility model, leads to downward-sloping volatility of variance smiles—in disagreement with variance markets in practice. Using control variates, we propose a robust method to express the Laplace transform of the variance call function in terms of the Laplace transform of the realized variance. The proposed method works in any model where the Laplace transform of realized variance is available in closed form. Additionally, we apply a new numerical Laplace inversion algorithm that gives fast and accurate prices for options on realized variance, simultaneously at a sequence of variance strikes. The method is also used to derive hedge ratios for options on variance with respect to variance swaps.     

The representation of American options prices under stochastic volatility and jump-diffusion dynamics

This paper considers the problem of pricing American options when the dynamics of the underlying are driven by both stochastic volatility following a square-root process as used by Heston [Rev. Financial Stud., 1993, 6, 327–343], and by a Poisson jump process as introduced by Merton [J. Financial Econ., 1976, 3, 125–144]. Probability arguments are invoked to find a representation of the solution in terms of expectations over the joint distribution of the underlying process. A combination of Fourier transform in the log stock price and Laplace transform in the volatility is then applied to find the transition probability density function of the underlying process. It turns out that the price is given by an integral dependent upon the early exercise surface, for which a corresponding integral equation is obtained. The solution generalizes in an intuitive way the structure of the solution to the corresponding European option pricing problem obtained by Scott [Math. Finance, 1997, 7(4), 413–426], but here in the case of a call option and constant interest rates.     

A novel Monte Carlo approach to hybrid local volatility models

We present in a Monte Carlo simulation framework, a novel approach for the evaluation of hybrid local volatility [Risk, 1994, 7, 18–20], [Int. J. Theor. Appl. Finance, 1998, 1, 61–110] models. In particular, we consider the stochastic local volatility model—see e.g. Lipton et al. [Quant. Finance, 2014, 14, 1899–1922], Piterbarg [Risk, 2007, April, 84–89], Tataru and Fisher [Quantitative Development Group, Bloomberg Version 1, 2010], Lipton [Risk, 2002, 15, 61–66]—and the local volatility model incorporating stochastic interest rates—see e.g. Atlan [ArXiV preprint math/0604316, 2006], Piterbarg [Risk, 2006, 19, 66–71], Deelstra and Rayée [Appl. Math. Finance, 2012, 1–23], Ren et al. [Risk, 2007, 20, 138–143]. For both model classes a particular (conditional) expectation needs to be evaluated which cannot be extracted from the market and is expensive to compute. We establish accurate and ‘cheap to evaluate’ approximations for the expectations by means of the stochastic collocation method [SIAM J. Numer. Anal., 2007, 45, 1005–1034], [SIAM J. Sci. Comput., 2005, 27, 1118–1139], [Math. Models Methods Appl. Sci., 2012, 22, 1–33], [SIAM J. Numer. Anal., 2008, 46, 2309–2345], [J. Biomech. Eng., 2011, 133, 031001], which was recently applied in the financial context [Available at SSRN 2529691, 2014], [J. Comput. Finance, 2016, 20, 1–19], combined with standard regression techniques. Monte Carlo pricing experiments confirm that our method is highly accurate and fast.     

Implied liquidity: Model sensitivity

The concept of implied liquidity originates from the conic finance theory and more precisely from the law of two prices where market participants buy from the market at the ask price and sell to the market at the lower bid price. The implied liquidity λ of any financial instrument is determined such that both model prices fit as well as possible the bid and ask market quotes. It reflects the liquidity of the financial instrument: the lower the λ, the higher the liquidity. The aim of this paper is to study the evolution of the implied liquidity pre- and post-crisis under a wide range of models and to study implied liquidity time series which could give an insight for future stochastic liquidity modeling. In particular, we perform a maximum likelihood estimation of the CIR, Vasicek and CEV mean-reverting processes applied to liquidity and volatility time series. The results show that implied liquidity is far less persistent than implied volatility as the liquidity process reverts much faster to its long-run mean. Moreover, a comparison of the parameter estimates between the pre- and post-credit crisis periods indicates that liquidity tends to decrease and increase for long and short term options, respectively, during troubled periods.     

Pricing Asian options with stochastic volatility

Abstract

In this paper, we generalize the recently developed dimension reduction technique of Vecer for pricing arithmetic average Asian options. The assumption of constant volatility in Vecer's method will be relaxed to the case that volatility is randomly fluctuating and is driven by a mean-reverting (or ergodic) process. We then use the fast mean-reverting stochastic volatility asymptotic analysis introduced by Fouque, Papanicolaou and Sircar to derive an approximation to the option price which takes into account the skew of the implied volatility surface. This approximation is obtained by solving a pair of one-dimensional partial differential equations.