American option pricing under the double Heston model based on asymptotic expansion

This paper focuses on pricing American put options under the double Heston model proposed by Christoffersen et al. By introducing an explicit exercise rule, we obtain the asymptotic expansion of the solution to the partial differential equation for pricing American put options. We calculate American option price by the sum of the European option price and the early exercise premium. The early exercise premium is calculated by the difference between the American and European option prices based on asymptotic expansions. The European option price is obtained by the efficient COS method. Based on the obtained American option price, the double Heston model is calibrated by minimizing the distance between model and market prices, which yields an optimization problem that is solved by a differential evolution algorithm combined with the Matlab function fmincon.m. Numerical results show that the pricing approach is fast and accurate. Empirical results show that the double Heston model has better performance in pricing short-maturity American put options and capturing the volatility term structure of American put options than the Heston model.     

Application of Heston’s Model to the Chinese Stock Market

This article applies Heston’s (1993) stochastic volatility model to the Chinese stock market indices and subsequently assesses its pricing performance. A two-step estimation procedure is adopted to calibrate Heston’s model. First, we find that the option price is affected by both the moneyness and the maturity. Second, Heston’s model is more likely to overprice options, whereas the BS model tends to underestimate options. Finally, Heston’s model, by employing volatility as a random process, significantly improves the pricing accuracy compared to the BS model. Therefore, Heston’s model is tractable to analyze the Chinese stock market indices, and there is volatility risk that must not be overlooked in the Chinese stock market.     

Option pricing in a Garch model with tempered stable innovations

The key problem for option pricing in Garch models is that the risk-neutral distribution of the underlying at maturity is unknown. Heston and Nandi solved this problem by computing the characteristic function of the underlying by a recursive procedure. Following the same idea, Christoffersen, Heston and Jacobs proposed a Garch-like model with inverse Gaussian innovations and recently Bellini and Mercuri obtained a similar procedure in a model with Gamma innovations. We present a model with tempered stable innovations that encompasses both the CHJ and the BM models as special cases. The proposed model is calibrated on S&P500 closing option prices and its performance is compared with the CHJ, the BM and the Heston–Nandi models.     

Affine forward variance models

We introduce the class of affine forward variance (AFV) models of which both the conventional Heston model and the rough Heston model are special cases. We show that AFV models can be characterised by the affine form of their cumulant-generating function, which can be obtained as solution of a convolution Riccati equation. We further introduce the class of affine forward order flow intensity (AFI) models, which are structurally similar to AFV models, but driven by jump processes, and which include Hawkes-type models. We show that the cumulant-generating function of an AFI model satisfies a generalised convolution Riccati equation and that a high-frequency limit of AFI models converges in distribution to an AFV model.

     

A comparison principle between rough and non-rough Heston models—with applications to the volatility surface

We present a number of related comparison results, which allow one to compare moment explosion times, moment generating functions and critical moments between rough and non-rough Heston models of stochastic volatility. All results are based on a comparison principle for certain non-linear Volterra integral equations. Our upper bound for the moment explosion time is different from the bound introduced by Gerhold, Gerstenecker and Pinter [Moment explosions in the rough Heston model. Decisions in Economics and Finance, 2019, 42, 575–608] and tighter for typical parameter values. The results can be directly transferred to a comparison principle for the asymptotic slope of implied variance between rough and non-rough Heston models. This principle shows that the ratio of implied variance slopes in the rough versus non-rough Heston model increases at least with power-law behavior for small maturities.     

A numerical method to price exotic path-dependent options on an underlying described by the Heston stochastic volatility model

We consider the problem of pricing European exotic path-dependent derivatives on an underlying described by the Heston stochastic volatility model. Lipton has found a closed form integral representation of the joint transition probability density function of underlying price and variance in the Heston model. We give a convenient numerical approximation of this formula and we use the obtained approximated transition probability density function to price discrete path-dependent options as discounted expectations. The expected value of the payoff is calculated evaluating an integral with the Monte Carlo method using a variance reduction technique based on a suitable approximation of the transition probability density function of the Heston model. As a test case, we evaluate the price of a discrete arithmetic average Asian option, when the average over n = 12 prices is considered, that is when the integral to evaluate is a 2n = 24 dimensional integral. We show that the method proposed is computationally efficient and gives accurate results.     

A multivariate stochastic volatility model with applications in the foreign exchange market

The main objective of this paper is to study the behavior of a daily calibration of a multivariate stochastic volatility model, namely the principal component stochastic volatility (PCSV) model, to market data of plain vanilla options on foreign exchange rates. To this end, a general setting describing a foreign exchange market is introduced. Two adequate models—PCSV and a simpler multivariate Heston model—are adjusted to suit the foreign exchange setting. For both models, characteristic functions are found which allow for an almost instantaneous calculation of option prices using Fourier techniques. After presenting the general calibration procedure, both the multivariate Heston and the PCSV models are calibrated to a time series of option data on three exchange rates—USD-SEK, EUR-SEK, and EUR-USD—spanning more than 11 years. Finally, the benefits of the PCSV model which we find to be superior to the multivariate extension of the Heston model in replicating the dynamics of these options are highlighted.     

Are classical option pricing models consistent with observed option second-order moments? Evidence from high-frequency data

As a means of validating an option pricing model, we compare the ex-post intra-day realized variance of options with the realized variance of the associated underlying asset that would be implied using assumptions as in the Black and Scholes (BS) model, the Heston, and the Bates model. Based on data for the S&P 500 index, we find that the BS model is strongly directionally biased due to the presence of stochastic volatility. The Heston model reduces the mismatch in realized variance between the two markets, but deviations are still significant. With the exception of short-dated options, we achieve best approximations after controlling for the presence of jumps in the underlying dynamics. Finally, we provide evidence that, although heavily biased, the realized variance based on the BS model contains relevant predictive information that can be exploited when option high-frequency data is not available.     

Heston stochastic vol-of-vol model for joint calibration of VIX and S&P 500 options

A parsimonious generalization of the Heston model is proposed where the volatility-of-volatility is assumed to be stochastic. We follow the perturbation technique of Fouque et al [Multiscale Stochastic Volatility for Equity, Interest Rate, and Credit Derivatives, 2011, Cambridge University Press] to derive a first-order approximation of the price of options on a stock and its volatility index. This approximation is given by Heston’s quasi-closed formula and some of its Greeks. It can be efficiently calculated since it requires to compute only Fourier integrals and the solution of simple ODE systems. We exemplify the calibration of the model with S&P 500 and VIX data.     

The microstructural foundations of leverage effect and rough volatility

We show that typical behaviors of market participants at the high frequency scale generate leverage effect and rough volatility. To do so, we build a simple microscopic model for the price of an asset based on Hawkes processes. We encode in this model some of the main features of market microstructure in the context of high frequency trading: high degree of endogeneity of market, no-arbitrage property, buying/selling asymmetry and presence of metaorders. We prove that when the first three of these stylized facts are considered within the framework of our microscopic model, it behaves in the long run as a Heston stochastic volatility model, where a leverage effect is generated. Adding the last property enables us to obtain a rough Heston model in the limit, exhibiting both leverage effect and rough volatility. Hence we show that at least part of the foundations of leverage effect and rough volatility can be found in the microstructure of the asset.     

The affine Heston model with correlated Gaussian interest rates for pricing hybrid derivatives

In this article we define a multi-factor equity–interest rate hybrid model with non-zero correlation between the stock and interest rate. The equity part is modeled by the Heston model and we use a Gaussian multi-factor short-rate process. By construction, the model fits in the framework of affine diffusion processes, allowing fast calibration to plain vanilla options. We also provide an efficient Monte Carlo simulation scheme.     

Asymptotic Expansion Formula of Option Price Under Multifactor Heston Model

The stochastic volatility model of Heston (Rev Financ Stud 6(2):327–343, 1993) has found difficulty in describing some of the important features of implied volatility dynamics, leading to a quest for multifactor extensions as well as the incorporation of time-dependent model parameters. In this paper, an asymptotic expansion approach to the multifactor Heston model with time-dependent parameters is developed. The results of Benhamou et al. (SIAM J Financ Math 1(1):289–325, 2010) are extended and it is shown that the extension to the multifactor model involves an extra expansion term that captures the interaction between variance factors. The expansion formula under constant parameters can be explicitly computed and the incorporation of time-dependent parameters is straightforward under the framework. As illustration, a two-factor model is calibrated to data of index options and variance swaps and it is found that it is possible to distinguish a short-term and long-term variance factor from the implied volatility surface and variance swap rates. Moreover, the two-factor model is able to reproduce the shapes of the implied volatility surface during various market scenarios.     

Valuing fade-in options with default risk in Heston–Nandi GARCH models

Review of Derivatives Research - In this paper, we present a pricing model to value fade-in options with default risk, where the underlying asset price is driven by the Heston–Nandi GARCH...     

Unspanned stochastic volatility and fixed income derivatives pricing

We propose a parsimonious ‘unspanned stochastic volatility’ model of the term structure and study its implications for fixed-income option prices. The drift and quadratic variation of the short rate are affine in three state variables (the short rate, its long-term mean and variance) which follow a joint Markov (vector) process. Yet, bond prices are exponential affine functions of only two state variables, independent of the current interest rate volatility level. Because this result holds for an arbitrary volatility process, such a process can be calibrated to match fixed income derivative prices. Furthermore, this model can be ‘extended’ (by relaxing the time-homogeneity) to fit any arbitrary term structure. In its ‘HJM’ form, this model nests the analogous stochastic equity volatility model of Heston (1993) [Heston, S.L., 1993. A closed form solution for options with stochastic volatility. Review of Financial Studies 6, 327–343]. In particular, if the volatility process is specified to be affine, closed-form solutions for interest rate options obtain. We propose an efficient algorithm to compute these prices. An application using data on caps and floors shows that the model can capture very well the implied Black spot volatility surface, while simultaneously fitting the observed term structure.     

Valuation of FX barrier options under stochastic volatility

This paper describes European-style valuation and hedging procedures for a class of knockout barrier options under stochastic volatility. A pricing framework is established by applying mean self-financing arguments and the minimal equivalent martingale measure. Using appropriate combinations of stochastic numerical and variance reduction procedures we demonstrate that fast and accurate valuations can be obtained for down-and-out call options for the Heston model.     

An empirical comparison of GARCH option pricing models

Recent empirical studies have shown that GARCH models can be successfully used to describe option prices. Pricing such contracts requires knowledge of the risk neutral cumulative return distribution. Since the analytical forms of these distributions are generally unknown, computationally intensive numerical schemes are required for pricing to proceed. Heston and Nandi (2000) consider a particular GARCH structure that permits analytical solutions for pricing European options and they provide empirical support for their model. The analytical tractability comes at a potential cost of realism in the underlying GARCH dynamics. In particular, their model falls in the affine family, whereas most GARCH models that have been examined fall in the non-affine family. This article takes a closer look at this model with the objective of establishing whether there is a cost to restricting focus to models in the affine family. We confirm Heston and Nandi's findings, namely that their model can explain a significant portion of the volatility smile. However, we show that a simple non affine NGARCH option model is superior in removing biases from pricing residuals for all moneyness and maturity categories especially for out-the-money contracts. The implications of this finding are examined. JEL Classification G13     

Pricing Exotic Options and American Options: A Multidimensional Asymptotic Expansion Approach

This paper introduces a new method for pricing exotic options whose payoff functions depend on several stochastic indices and American options in multidimensional models. This method is based on two ideas. One is an application of the asymptotic expansion method for the law of a multidimensional diffusion process. The other is the combination of the asymptotic expansion method and the method called backward induction. The author applies the method to the problems of pricing call options on the maximum of two assets in the CEV model, average strike options in the Black–Scholes model and American options in the Heston model. Numerical examples show practical effectiveness of the proposed method.     

Gamma expansion of the Heston stochastic volatility model

We derive an explicit representation of the transitions of the Heston stochastic volatility model and use it for fast and accurate simulation of the model. Of particular interest is the integral of the variance process over an interval, conditional on the level of the variance at the endpoints. We give an explicit representation of this quantity in terms of infinite sums and mixtures of gamma random variables. The increments of the variance process are themselves mixtures of gamma random variables. The representation of the integrated conditional variance applies the Pitman–Yor decomposition of Bessel bridges. We combine this representation with the Broadie–Kaya exact simulation method and use it to circumvent the most time-consuming step in that method.     

Algorithmic market making for options

In this article, we tackle the problem of a market maker in charge of a book of options on a single liquid underlying asset. By using an approximation of the portfolio in terms of its vega, we show that the seemingly high-dimensional stochastic optimal control problem of an option market maker is in fact tractable. More precisely, when volatility is modeled using a classical stochastic volatility model—e.g. the Heston model—the problem faced by an option market maker is characterized by a low-dimensional functional equation that can be solved numerically using a Euler scheme along with interpolation techniques, even for large portfolios. In order to illustrate our findings, numerical examples are provided.     

Maximum likelihood estimation of stochastic volatility models

We develop and implement a method for maximum likelihood estimation in closed-form of stochastic volatility models. Using Monte Carlo simulations, we compare a full likelihood procedure, where an option price is inverted into the unobservable volatility state, to an approximate likelihood procedure where the volatility state is replaced by proxies based on the implied volatility of a short-dated at-the-money option. The approximation results in a small loss of accuracy relative to the standard errors due to sampling noise. We apply this method to market prices of index options for several stochastic volatility models, and compare the characteristics of the estimated models. The evidence for a general CEV model, which nests both the affine Heston model and a GARCH model, suggests that the elasticity of variance of volatility lies between that assumed by the two nested models.