BLACK–SCHOLES REPRESENTATION FOR ASIAN OPTIONS

Asian options are securities with a payoff that depends on the average of the underlying stock price over a certain time interval. We identify three natural assets that appear in pricing of the Asian options, namely a stock S, a zero coupon bond B^T with maturity T, and an abstract asset A (an “average asset”) that pays off a weighted average of the stock price number of units of a dollar at time T. It turns out that each of these assets has its own martingale measure, allowing us to obtain Black–Scholes type formulas for the fixed strike and the floating strike Asian options. The model independent formulas are analogous to the Black–Scholes formula for the plain vanilla options; they are expressed in terms of probabilities under the corresponding martingale measures that the Asian option will end up in the money. Computation of these probabilities is relevant for hedging. In contrast to the plain vanilla options, the probabilities for the Asian options do not admit a simple closed form solution. However, we show that it is possible to obtain the numerical values in the geometric Brownian motion model efficiently, either by solving a partial differential equation numerically, or by computing the Laplace transform. Models with stochastic volatility or pure jump models can be also priced within the Black–Scholes framework for the Asian options.     

Analytical valuation of Asian options with counterparty risk under stochastic volatility models

In this paper, we consider Asian options with counterparty risk under stochastic volatility models. We propose a simple way to construct stochastic volatility models through the market factor channel. In the proposed framework, we obtain an explicit pricing formula of Asian options with counterparty risk and illustrate the effects of systematic risk on Asian option prices. Specially, the U-shaped and inverted U-shaped curves appear when we keep the total risk of the underlying asset and the issuer's assets unchanged, respectively.     

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.     

ASYMPTOTIC BEHAVIOR OF DISTRIBUTION DENSITIES IN MODELS WITH STOCHASTIC VOLATILITY. I

We study the asymptotic behavior of distribution densities arising in stock price models with stochastic volatility. The main objects of our interest in the present paper are the density of time averages of a geometric Brownian motion and the density of the stock price process in the Hull–White model. We find explicit formulas for leading terms in asymptotic expansions of these densities and give error estimates. As an application of our results, sharp asymptotic formulas for the price of an Asian option are obtained.     

Pricing Options With Curved Boundaries1

This paper provides a general valuation method for the European options whose payoff is restricted by curved boundaries contractually set on the underlying asset price process when it follows the geometric Brownian motion. Our result is based on the generalization of the Levy formula on the Brownian motion by T. W. Anderson in sequential analysis. We give the explicit probability formula that the geometric Brownian motion reaches in an interval at the maturity date without hitting either the lower or the upper curved boundaries. Although the general pricing formulae for options with boundaries are expressed as infinite series in the general case, our numerical study suggests that the convergence of the series is rapid. Our results include the formulae for options with a lower boundary by Merton (1973), for path-dependent options by Goldman, Sossin, and Gatto (1979), and for some corporate securities as special cases.     

ANALYTIC APPROXIMATIONS FOR MULTI‐ASSET OPTION PRICING

We derive general analytic approximations for pricing European basket and rainbow options on N assets. The key idea is to express the option’s price as a sum of prices of various compound exchange options, each with different pairs of subordinate multi‐ or single‐asset options. The underlying asset prices are assumed to follow lognormal processes, although our results can be extended to certain other price processes for the underlying. For some multi‐asset options a strong condition holds, whereby each compound exchange option is equivalent to a standard single‐asset option under a modified measure, and in such cases an almost exact analytic price exists. More generally, approximate analytic prices for multi‐asset options are derived using a weak lognormality condition, where the approximation stems from making constant volatility assumptions on the price processes that drive the prices of the subordinate basket options. The analytic formulae for multi‐asset option prices, and their Greeks, are defined in a recursive framework. For instance, the option delta is defined in terms of the delta relative to subordinate multi‐asset options, and the deltas of these subordinate options with respect to the underlying assets. Simulations test the accuracy of our approximations, given some assumed values for the asset volatilities and correlations. Finally, a calibration algorithm is proposed and illustrated.     

PRICING A CLASS OF EXOTIC OPTIONS VIA MOMENTS AND SDP RELAXATIONS

We present a new methodology for the numerical pricing of a class of exotic derivatives such as Asian or barrier options when the underlying asset price dynamics are modeled by a geometric Brownian motion or a number of mean-reverting processes of interest. This methodology identifies derivative prices with infinite-dimensional linear programming problems involving the moments of appropriate measures, and then develops suitable finite-dimensional relaxations that take the form of semidefinite programs (SDP) indexed by the number of moments involved. By maximizing or minimizing appropriate criteria, monotone sequences of both upper and lower bounds are obtained. Numerical investigation shows that very good results are obtained with only a small number of moments. Theoretical convergence results are also established.     

A MULTINOMIAL APPROXIMATION FOR AMERICAN OPTION PRICES IN LÉVY PROCESS MODELS

This paper gives a tree-based method for pricing American options in models where the stock price follows a general exponential Lévy process. A multinomial model for approximating the stock price process, which can be viewed as generalizing the binomial model of Cox, Ross, and Rubinstein (1979) for geometric Brownian motion, is developed. Under mild conditions, it is proved that the stock price process and the prices of American-type options on the stock, calculated from the multinomial model, converge to the corresponding prices under the continuous time Lévy process model. Explicit illustrations are given for the variance gamma model and the normal inverse Gaussian process when the option is an American put, but the procedure is applicable to a much wider class of derivatives including some path-dependent options. Our approach overcomes some practical difficulties that have previously been encountered when the Lévy process has infinite activity.     

The valuation of European options when asset returns are autocorrelated

This article derives the closed‐form formula for a European option on an asset with returns following a continuous‐time type of first‐order moving average process, which is called an MA(1)‐type option. The pricing formula of these options is similar to that of Black and Scholes, except for the total volatility input. Specifically, the total volatility input of MA(1)‐type options is the conditional standard deviation of continuous‐compounded returns over the option's remaining life, whereas the total volatility input of Black and Scholes is indeed the diffusion coefficient of a geometric Brownian motion times the square root of an option's time to maturity. Based on the result of numerical analyses, the impact of autocorrelation induced by the MA(1)‐type process is significant to option values even when the autocorrelation between asset returns is weak. © 2006 Wiley Periodicals, Inc. Jrl Fut Mark 26:85–102, 2006     

Robustness of the Black and Scholes Formula

Consider an option on a stock whose volatility is unknown and stochastic. An agent assumes this volatility to be a specific function of time and the stock price, knowing that this assumption may result in a misspecification of the volatility. However, if the misspecified volatility dominates the true volatility, then the misspecified price of the option dominates its true price. Moreover, the option hedging strategy computed under the assumption of the misspecified volatility provides an almost sure one-sided hedge for the option under the true volatility. Analogous results hold if the true volatility dominates the misspecified volatility. These comparisons can fail, however, if the misspecified volatility is not assumed to be a function of time and the stock price. The positive results, which apply to both European and American options, are used to obtain a bound and hedge for Asian options.     

PRICING AND HEDGING DOUBLE-BARRIER OPTIONS: A PROBABILISTIC APPROACH

Barrier options have become increasingly popular over the last few years. Less expensive than standard options, they may provide the appropriate hedge in a number of risk management strategies. In the case of a single-barrier option, the valuation problem is not very difficult (see Merton 1973 and Goldman, Sosin, and Gatto 1979). the situation where the option gets knocked out when the underlying instrument hits either of two well-defined boundaries is less straightforward. Kunitomo and Ikeda (1992) provide a pricing formula expressed as the sum of an infinite series whose convergence is studied through numerical procedures and suggested to be rapid. We follow a methodology which proved quite successful in the case of Asian options (see Geman and Yor 1992,1993) and which has its roots in some fundamental properties of Brownian motion. This methodology permits the derivation of a simple expression of the Laplace transform of the double-barrir price with respect to its maturity date. the inversion of the Laplace transform using techniques developed by Geman and Eydeland (1995), is then fairly easy to perform.     

Laguerre Series for Asian and Other Options

This paper has four goals: (a) relate ladder height distributions to option values; (b) show how Laguerre expansions may be used in the computation of densities, distribution functions, and option prices; (c) derive some new results on the integral of geometric Brownian motion over a finite interval; and (d) apply the preceding results to the determination of the distribution of the integral of geometric Brownian motion and the computation of Asian option values. The usual fixed‐strike options on the average are treated, as well as options with payoffs expressed in terms of one over the average of the underlying security, which this author calls “reciprocal Asian options.” In all cases the underlying asset is represented by geometric Brownian motion, the averages are performed continuously, and the options are of European type.     

A Continuity Correction for Discrete Barrier Options

The payoff of a barrier option depends on whether or not a specified asset price, index, or rate reaches a specified level during the life of the option. Most models for pricing barrier options assume continuous monitoring of the barrier; under this assumption, the option can often be priced in closed form. Many (if not most) real contracts with barrier provisions specify discrete monitoring instants; there are essentially no formulas for pricing these options, and even numerical pricing is difficult. We show, however, that discrete barrier options can be priced with remarkable accuracy using continuous barrier formulas by applying a simple continuity correction to the barrier. The correction shifts the barrier away from the underlying by a factor of exp(bet sig sqrt dt), where bet approx 0.5826, sig is the underlying volatility, and dt is the time between monitoring instants. The correction is justified both theoretically and experimentally.     

Analytic approximation formulae for pricing forward‐starting Asian options

In this article we first identify a missing term in the Bouaziz, Briys, and Crouhy ( 1994 ) pricing formula for forward‐starting Asian options and derive the correct one. First, illustrate in certain cases that the missing term in their pricing formula could induce large pricing errors or unreasonable option prices. Second, we derive new analytic approximation formulae for valuing forward‐starting Asian options by adding the second‐order term in the Taylor series. We show that our formulae can accurately value forward‐starting Asian options with a large underlying asset's volatility or a longer time window for the average of the underlying asset prices, whereas the pricing errors for these options with the previously mentioned formula could be large. Third, we derive the hedge ratios for these options and compare their properties with those of plain vanilla options. © 2003 Wiley Periodicals, Inc. Jrl Fut Mark 23:487–516, 2003     

On peacocks and lyrebirds: Australian options,Brownian bridges,and the average of submartingales

We introduce a class of stochastic processes, which we refer to as lyrebirds. These extend a class of stochastic processes, which have recently been coined peacocks, but are more commonly known as processes that are increasing in the convex order. We show how these processes arise naturally in the context of Asian and Australian options and consider further applications, such as the arithmetic average of a Brownian bridge and the average of submartingales, including the case of Asian and Australian options where the underlying features constant elasticity of variance or is of Merton jump diffusion type.     

Pricing Discrete European Barrier Options Using Lattice Random Walks

This paper designs a numerical procedure to price discrete European barrier options in Black-Scholes model. The pricing problem is divided into a series of initial value problems, one for each monitoring time. Each initial value problem is solved by replacing the driving Brownian motion by a lattice random walk. Some results from the theory of Besov spaces show that the convergence rate of lattice methods for initial value problems depends on two factors, namely the smoothness of the initial value (or the value function) and the moments for the increments of the lattice random walk. This fact is used to obtain an efficient method to price discrete European barrier options. Numerical examples and comparisons with other methods are carried out to show that the proposed method yields fast and accurate results.     

On commodity price limits

This paper examines the behavior of futures prices and trader positions around the occurrence of price limits in commodity futures markets. We ask whether limit events are the result of shocks to fundamental volatility or the result of temporary volatility induced by the trading of noncommercial market participants (speculators). We find little evidence that limits events are the result of speculative activity, but instead associated with shocks to fundamentals that lead to persistent price changes. When futures trading halts price discovery migrates to options markets, but option prices provide a biased estimate of subsequent future prices when trading resumes.     

AN ANALYTICAL SOLUTION FOR THE TWO‐SIDED PARISIAN STOPPING TIME,ITS ASYMPTOTICS,AND THE PRICING OF PARISIAN OPTIONS

In this paper, we obtain a recursive formula for the density of the two‐sided Parisian stopping time. This formula does not require any numerical inversion of Laplace transforms, and is similar to the formula obtained for the one‐sided Parisian stopping time derived in Dassios and Lim. However, when we study the tails of the two distributions, we find that the two‐sided stopping time has an exponential tail, while the one‐sided stopping time has a heavier tail. We derive an asymptotic result for the tail of the two‐sided stopping time distribution and propose an alternative method of approximating the price of the two‐sided Parisian option.     

Enhancing managerial equity incentives with moving average payoffs

Prior research suggests that Asian stock options provide stronger managerial equity incentives than traditional stock options do, holding the cost of the option grant constant. Although this is true on the grant date, it is not over the life of the option grant. Very little of the initial advantage remains after 2 years because Asian stock options have diminishing incentive effects over time. A simple solution is to replace averaging over the option's life with averaging over a moving window. We show that moving average options do not have the diminishing incentive problem and are effective in preventing managerial gaming.     

A Simplified Quadrature Approach for Computing Bermudan Option Prices

We examine a simple quadrature approach to compute the prices of Bermudan options when the value of the corresponding European claim can be computed in closed form, one period before maturity. Using a constant grid of stock prices at early exercise time points, the known value of the European option is used as a smoothing device to enable efficient numerical integ ration with quadrature approaches. Examples with the geometric Brownian motion context and the lognormal jump‐diffusion context are provided.