Pricing VIX options with volatility clustering

We investigate the valuation of volatility index (VIX) options by developing a model with a self-exciting Hawkes process that allows for clustering in the VIX. In the proposed framework, we find semianalytical expressions for the characteristic function and forward characteristic function, and then we solve the pricing problem of standard-start and forward-start options via the fast Fourier transform. The empirical results provide evidence to support the significance of accounting for volatility clustering when pricing VIX options.     

Distributions implied by American currency futures options: A ghost's smile?

A new and easily applicable method for estimating risk‐neutral distributions (RND) implied by American futures options is proposed. It amounts to inverting the Barone‐Adesi and Whaley method (BAW method) to get the BAW implied volatility smile. Extensive empirical tests show that the BAW smile is equivalent to the volatility smile implied by corresponding European options. Therefore, the procedure leads to a legitimate RND estimation method. Further, the investigation of the currency options traded on the Chicago Mercantile Exchange and OTC markets in parallel provides us with insights on the structure and interaction of the two markets. Unequally distributed liquidity in the OTC market seems to lead to price distortions and an ensuing interesting “ghost‐like” shape of the RND density implied by CME options. Finally, using the empirical results, we propose a parsimonious generalization of the existing methods for estimating volatility smiles from OTC options. A single free parameter significantly improves the fit. © 2004 Wiley Periodicals, Inc. Jrl Fut Mark 24:147–178, 2004     

DYNAMIC SPANNING: ARE OPTIONS AN APPROPRIATE INSTRUMENT?1

Ross (1976) has shown, in a static framework, how options can complete financial markets. This paper examines the possible extensions of Ross's idea in a dynamic setup. Surprisingly enough, we find that the answer is very sensitive to the choice of the stochastic model for the underlying security returns. More specifically we obtain the following results: In a discrete-time model, classical European options typically become redundant with some probability (Proposition 2.1). Obnly path dependent (“exotic”) options may generate dynamic spanning (Proposition 4.1). In a continuous-time model with stochastic volatility of the underlying security, and under reasonable assumptions, a European option is always a good instrument for completing markets (Proposition 5.2).     

OPTION HEDGING AND IMPLIED VOLATILITIES IN A STOCHASTIC VOLATILITY MODEL1

In the stochastic volatility framework of Hull and White (1987), we characterize the so-called Black and Scholes implied volatility as a function of two arguments the ratio of the strike to the underlying asset price and the instantaneous value of the volatility By studying the variation m the first argument, we show that the usual hedging methods, through the Black and Scholes model, lead to an underhedged (resp. overhedged) position for in-the-money (resp out-of the-money) options, and a perfect partial hedged position for at the-money options These results are shown to be closely related to the smile effect, which is proved to be a natural consequence of the stochastic volatility feature the deterministic dependence of the implied volatility on the underlying volatility process suggests the use of implied volatility data for the estimation of the parameters of interest A statistical procedure of filtering (of the latent volatility process) and estimation (of its parameters) is shown to be strongly consistent and asymptotically normal.     

Applications of Eigenfunction Expansions in Continuous-Time Finance

We provide exact solutions for two closely related valuation problems in continuous-time finance. The first problem is to value generalized European-style options on stocks that pay dividends at a constant dollar rate. The second problem is to find the yield curve associated with the economy of R. C. Merton's "An Asymptotic Theory of Growth Under Uncertainty." In Merton's economic growth model, the interest rate process has a volatility linear in the rate level and a linear/quadratic drift. Both problems are solved by an eigenfunction expansion technique. The main technical difficulty is handling the problem of payoff functions that are not square-integrable with respect to the natural weight function of the models.     

Efficient trinomial trees for local-volatility models in pricing double-barrier options

A local-volatility (LV) model captures the volatility smile while retaining the preference freedom of the Black–Scholes model. Past attempts to construct a smile-consistent tree for the LV surface do not guarantee validity. This paper presents an efficient and valid smile-consistent tree for the LV model. The only assumption is that the LV surface be upper- and lower-bounded. With this tree, double-barrier options can be priced with fast convergence even in the presence of volatility smile. This is confirmed numerically. An implied tree is also presented. It recovers the LV surface reasonably well.     

Fitting and testing for the implied volatility curve using parametric models

Numerous issues have arisen over the past few decades relating to the implied volatility smile in the options market; however, the extant literature reveals that relatively little effort has thus far been placed into comparing the various implied volatility models, essentially as a result of the lack of any theoretical foundation on which to base such comparative analysis. In this study, we use a comprehensive options database and employ methods of combining the various hypothesis tests to compare the different implied volatility models. To the best of our knowledge, this is the first study of its kind to address this issue using combination tests. Our empirical results reveal that the linear piecewise model is the most appropriate model for capturing the implied volatility smile, with additional robustness checks confirming the validity of this finding.     

A NEW LOOK AT SHORT‐TERM IMPLIED VOLATILITY IN ASSET PRICE MODELS WITH JUMPS

We analyze the behavior of the implied volatility smile for options close to expiry in the exponential Lévy class of asset price models with jumps. We introduce a new renormalization of the strike variable with the property that the implied volatility converges to a nonconstant limiting shape, which is a function of both the diffusion component of the process and the jump activity (Blumenthal–Getoor) index of the jump component. Our limiting implied volatility formula relates the jump activity of the underlying asset price process to the short‐end of the implied volatility surface and sheds new light on the difference between finite and infinite variation jumps from the viewpoint of option prices: in the latter, the wings of the limiting smile are determined by the jump activity indices of the positive and negative jumps, whereas in the former, the wings have a constant model‐independent slope. This result gives a theoretical justification for the preference of the infinite variation Lévy models over the finite variation ones in the calibration based on short‐maturity option prices.     

Exercise Regions And Efficient Valuation Of American Lookback Options

This paper presents an efficient method to compute the values and early exercise boundaries of American fixed strike lookback options. The method reduces option valuation to a single optimal stopping problem for standard Brownian motion and an associated path-dependent functional, indexed by one parameter in the absence of dividends and by two parameters in the presence of a dividend rate. Numerical results obtained by this method show that, after a space-time transformation, the stopping boundaries are well approximated by certain piecewise linear functions with a few pieces, leading to fast and accurate approximations for American lookback option values. An explicit decomposition formula for American lookback options is derived and applied not only to the development of these approximations but also to the asymptotic analysis of the early exercise boundary near the expiration date.     

SERIES EXPANSION OF THE SABR JOINT DENSITY

Under the SABR stochastic volatility model, pricing and hedging contracts that are sensitive to forward smile risk (e.g., forward starting options, barrier options) require the joint transition density. In this paper, we address this problem by providing closed‐form representations, asymptotically, of the joint transition density. Specifically, we construct an expansion of the joint density through a hierarchy of parabolic equations after applying total volatility‐of‐volatility scaling and a near‐Gaussian coordinate transformation. We then establish an existence result to characterize the truncation error and provide explicit joint density formulas for the first three orders. Our approach inherits the same spirit of a small total volatility‐of‐volatility assumption as in the original SABR analysis. Our results for the joint transition density serve as a basis for managing forward smile risk. Through numerical experiments, we illustrate the accuracy of our expansion in terms of joint density, marginal density, probability mass, and implied volatilities for European call options.     

BESSEL PROCESSES,STOCHASTIC VOLATILITY,AND TIMER OPTIONS

Motivated by analytical valuation of timer options (an important innovation in realized variance‐based derivatives), we explore their novel mathematical connection with stochastic volatility and Bessel processes (with constant drift). Under the Heston (1993) stochastic volatility model, we formulate the problem through a first‐passage time problem on realized variance, and generalize the standard risk‐neutral valuation theory for fixed maturity options to a case involving random maturity. By time change and the general theory of Markov diffusions, we characterize the joint distribution of the first‐passage time of the realized variance and the corresponding variance using Bessel processes with drift. Thus, explicit formulas for a useful joint density related to Bessel processes are derived via Laplace transform inversion. Based on these theoretical findings, we obtain a Black–Scholes–Merton‐type formula for pricing timer options, and thus extend the analytical tractability of the Heston model. Several issues regarding the numerical implementation are briefly discussed.     

Valuation by Simulation of Contingent Claims with Multiple Early Exercise Opportunities

This paper introduces the application of Monte Carlo simulation technology to the valuation of securities that contain many (buying or selling) rights, but for which a limited number can be exercised per period, and penalties if a minimum quantity is not exercised before maturity. These securities combine the characteristics of American options, with the additional constraint that only a few rights can be exercised per period and therefore their price depends also on the number of living rights (i.e., American-Asian-style payoffs), and forward securities. These securities give flexibility-of-delivery options and are common in energy markets (e.g., take-or-pay or swing options) and as real options (e.g., the development of a mine). First, we derive a series of properties for the price and the optimal exercise frontier of these securities. Second, we price them by simulation, extending the Ibáñez and Zapatero (2004) method to this problem.     

VALUATION OF BARRIER OPTIONS VIA A GENERAL SELF‐DUALITY

Classical put–call symmetry relates the price of puts and calls under a suitable dual market transform. One well‐known application is the semistatic hedging of path‐dependent barrier options with European options. This, however, in its classical form requires the price process to observe rather stringent and unrealistic symmetry properties. In this paper, we develop a general self‐duality theorem to develop valuation schemes for barrier options in stochastic volatility models with correlation.     

Optimal No‐Arbitrage Bounds on S&P 500 Index Options and the Volatility Smile

This article shows that the volatility smile is not necessarily inconsistent with the Black–Scholes analysis. Specifically, when transaction costs are present, the absence of arbitrage opportunities does not dictate that there exists a unique price for an option. Rather, there exists a range of prices within which the option's price may fall and still be consistent with the Black–Scholes arbitrage pricing argument. This article uses a linear program (LP) cast in a binomial framework to determine the smallest possible range of prices for Standard & Poor's 500 Index options that are consistent with no arbitrage in the presence of transaction costs. The LP method employs dynamic trading in the underlying and risk‐free assets as well as fixed positions in other options that trade on the same underlying security. One‐way transaction‐cost levels on the index, inclusive of the bid–ask spread, would have to be below six basis points for deviations from Black–Scholes pricing to present an arbitrage opportunity. Monte Carlo simulations are employed to assess the hedging error induced with a 12‐period binomial model to approximate a continuous‐time geometric Brownian motion. Once the risk caused by the hedging error is accounted for, transaction costs have to be well below three basis points for the arbitrage opportunity to be profitable two times out of five. This analysis indicates that market prices that deviate from those given by a constant‐volatility option model, such as the Black–Scholes model, can be consistent with the absence of arbitrage in the presence of transaction costs. © 2001 John Wiley & Sons, Inc. Jrl Fut Mark 21:1151–1179, 2001     

Accounting for stochastic interest rates,stochastic volatility and a general correlation structure in the valuation of forward starting options

A quantitative analysis on the pricing of forward starting options under stochastic volatility and stochastic interest rates is performed. The main finding is that forward starting options not only depend on future smiles, but also directly on the evolution of the interest rates as well as the dependency structures among the underlying asset, the interest rates, and the stochastic volatility: compared to vanilla options, dynamic structures such as forward starting options are much more sensitive to model specifications such as volatility, interest rate, and correlation movements. We conclude that it is of crucial importance to take all these factors explicitly into account for a proper valuation and risk management of these securities. The performed analysis is facilitated by deriving closed‐form formulas for the valuation of forward starting options, hereby taking the stochastic volatility, stochastic interest rates as well the dependency structure between all these processes explicitly into account. The valuation framework is derived using a probabilistic approach, enabling a fast and efficient evaluation of the option price by Fourier inverting the forward starting characteristic functions. © 2010 Wiley Periodicals, Inc. Jrl Fut Mark 31:103–125, 2011     

The performance of traders' rules in options market

This study focuses on the usefulness of the traders' rules to predict future implied volatilities for pricing and hedging KOSPI 200 index options. There are two versions of this approach. In the “relative smile” approach, the implied volatility skew is treated as a fixed function of moneyness. In the “absolute smile” approach, the implied volatility skew is treated as a fixed function of the strike price. It is found that the “absolute smile” approach shows better performance than Black, F. and Scholes, L. ( 1973 ) model and the stochastic volatility model for both pricing and hedging options. Consistent with Jackwerth, J. C. and Rubinstein, M. (2001) and Li, M. and Pearson, N. D. (2007), the traders' rules dominate mathematically more sophisticated model, that is, the stochastic volatility model. The traders' rules can be an alternative to the sophisticated and complicated models for pricing and hedging options. © 2009 Wiley Periodicals, Inc. Jrl Fut Mark 29:999–1020, 2009     

BOUNDARY EVOLUTION EQUATIONS FOR AMERICAN OPTIONS

We consider the problem of finding optimal exercise policies for American options, both under constant and stochastic volatility settings. Rather than work with the usual equations that characterize the price exclusively, we derive and use boundary evolution equations that characterize the evolution of the optimal exercise boundary. Using these boundary evolution equations we show how one can construct very efficient computational methods for pricing American options that avoid common sources of error. First, we detail a methodology for standard static grids and then describe an improvement that defines a grid that evolves dynamically while solving the problem. When integral representations are available, as in the Black–Scholes setting, we also describe a modified integral method that leverages on the representation to solve the boundary evolution equations. Finally we compare runtime and accuracy to other popular numerical methods. The ideas and methodology presented herein can easily be extended to other optimal stopping problems.     

Currency option pricing: Mean reversion and multi‐scale stochastic volatility

This paper investigates the valuation of currency options when the underlying currency follows a mean‐reverting lognormal process with multi‐scale stochastic volatility. A closed‐form solution is derived for the characteristic function of the log‐asset price. European options are then valued by means of the Fourier inversion formula. The proposed model enables us to calibrate simultaneously to the observed currency futures and the implied volatility surface of the currency options within a unified framework. The fractional fast Fourier transform (FFT) is adopted to implement the Fourier inversion, thus ensuring that the grid spacing restriction of the standard FFT can be relaxed, which results in a more efficient computation. Using Monte Carlo simulation as a benchmark, our numerical examples show that the derived option pricing formula is accurate and efficient for practical use. © 2009 Wiley Periodicals, Inc. Jrl Fut Mark 30:938–956, 2010     

Empirical tests of canonical nonparametric American option‐pricing methods

Alcock and Carmichael (2008, The Journal of Futures Markets, 28, 717–748) introduce a nonparametric method for pricing American‐style options, that is derived from the canonical valuation developed by Stutzer (1996, The Journal of Finance, 51, 1633–1652). Although the statistical properties of this nonparametric pricing methodology have been studied in a controlled simulation environment, no study has yet examined the empirical validity of this method. We introduce an extension to this method that incorporates information contained in a small number of observed option prices. We explore the applicability of both the original method and our extension using a large sample of OEX American index options traded on the S&P100 index. Although the Alcock and Carmichael method fails to outperform a traditional implied‐volatility‐based Black–Scholes valuation or a binomial tree approach, our extension generates significantly lower pricing errors and performs comparably well to the implied‐volatility Black–Scholes pricing, in particular for out‐of‐the‐money American put options. © 2009 Wiley Periodicals, Inc. Jrl Fut Mark 30:509–532, 2010     

PRICING OF AMERICAN PATH-DEPENDENT CONTINGENT CLAIMS

We consider the problem of pricing path-dependent contingent claims. Classically, this problem can be cast into the Black-Scholes valuation framework through inclusion of the path-dependent variables into the state space. This leads to solving a degenerate advection-diffusion partial differential equation (PDE). We first estabilish necessary and sufficient conditions under which degenerate diffusions can be reduced to lower-dimensional nondegenerate diffusions. We apply these results to path-dependent options. Then, we describe a new numerical technique, called forward shooting grid (FSG) method, that efficiently copes with degenerate diffusion PDEs. Finally, we show that the FSG method is unconditionally stable and convergent. the FSG method is the first capable of dealing with the early exercise condition of American options. Several numerical examples are presented and discussed. ²