An application of finite elements to option pricing

This study applied the finite element method (FEM) to pricing options. The FEM estimates the function that satisfies a governing differential equation through the assembly of piecewise continuous functions over the domain of the problem. Two common representations, a variational functional representation, and a weighted residual representation are used in the application of the method. The FEM is a versatile alternative to other popular lattice methods used in option pricing. Advantages include the abilities to directly estimate the Greeks of the option and allow nonuniform mesh construction. As an illustration of the advantages that the FEM offers, the method was used to price European put options and discrete barrier knock‐out put options. © 2001 John Wiley & Sons, Inc. Jrl Fut Mark 21:19–42, 2001     

THE BLACK-SCHOLES EQUATION REVISITED: ASYMPTOTIC EXPANSIONS AND SINGULAR PERTURBATIONS

In this paper, novel singular perturbation techniques are applied to price European, American, and barrier options. Employment of these methods leads to a significant simplification of the problem in all cases, by reducing the number of parameters. For American options, the valuation problem is reduced to a procedure that may be performed on a rudimentary handheld calculator. The method also sheds light on the evolution of option prices for all of the cases considered, the results being particularly illuminating for American and barrier options.     

Improving lattice schemes through bias reduction

Lattice schemes for option pricing, such as tree or grid/partial differential equation (p.d.e.) methods, are usually designed as a discrete version of an underlying continuous model of stock prices. The parameters of such schemes are chosen so that the discrete version “best” matches the continuous one. Only in the limit does the lattice option price model converge to the continuous one. Otherwise, a discretization bias remains. A simple modification of lattice schemes which reduces the discretization bias is proposed. The modification can, in theory, be applied to any lattice scheme. The main idea is to adjust the lattice parameters in such a way that the option price bias, not the stock price bias, is minimized. European options are used, for which the option price bias can be evaluated precisely, as a template to modify and improve American option methods. A numerical study is provided. © 2006 Wiley Periodicals, Inc. Jrl Fut Mark 26:733–757, 2006     

An efficient approximation method for American exotic options

The authors suggest a modified quadratic approximation scheme, and apply this scheme to American barrier (knock‐out) and floating‐strike lookback options. This modified scheme introduces an additional parameter into the quadratic approximation method, originally suggested by G. Barone‐Adesi and R. Whaley (1987), to reduce pricing errors. When the barrier is close to the underlying asset's current price, the approximation formula is more accurate than lattice methods because the optimal exercise boundary is independent of the underlying asset's current price. That is, the proposed method overcomes the “near‐barrier” problem that occurs in lattice methods. In addition, the pricing error decreases when the underlying asset's volatility is high. This approximation scheme is more efficient than B. Gao, J. Huang, and M. Subrahmanyam's (2000) method. As a second application of the modified approximation scheme, the authors provide an approximation formula for American floating‐strike lookback options which is the first approximation formula ever suggested in the literature. Compared to S. Babbs' (2000) binomial approach, our approximation method is more efficient after controlling for pricing errors, and is more accurate after controlling for computing time. © 2007 Wiley Periodicals, Inc. Jrl Fut Mark 27:29–59, 2007     

THE DEPENDENCE STRUCTURE OF RUNNING MAXIMA AND MINIMA: RESULTS AND OPTION PRICING APPLICATIONS

We provide general results for the dependence structure of running maxima (minima) of sets of variables in a model based on (i) Markov dynamics; (ii) no Granger causality; (iii) cross-section dependence. At the time series level, we derive recursive formulas for running minima and maxima. These formulas enable to use a "bootstrapping" technique to recursively recover the pricing kernels of barrier options from those of the corresponding European options. We also show that the dependence formulas for running maxima (minima) are completely defined from the copula function representing dependence among levels at the terminal date. The result is applied to multivariate discrete barrier digital products. Barrier Altiplanos are simply priced by (i) bootstrapping the price of univariate barrier products; (ii) evaluating a European Altiplano with these values.     

Pricing American options by canonical least‐squares Monte Carlo

Options pricing and hedging under canonical valuation have recently been demonstrated to be quite effective, but unfortunately are only applicable to European options. This study proposes an approach called canonical least‐squares Monte Carlo (CLM) to price American options. CLM proceeds in three stages. First, given a set of historical gross returns (or price ratios) of the underlying asset for a chosen time interval, a discrete risk‐neutral distribution is obtained via the canonical approach. Second, from this canonical distribution independent random samples of gross returns are taken to simulate future price paths for the underlying. Third, to those paths the least‐squares Monte Carlo algorithm is then applied to obtain early exercise strategies for American options. Numerical results from simulation‐generated gross returns under geometric Brownian motions show that the proposed method yields reasonably accurate prices for American puts. The CLM method turns out to be quite similar to the nonparametric approach of Alcock and Carmichael and simulations done with CLM provide additional support for their recent findings. CLM can therefore be viewed as an alternative for pricing American options, and perhaps could even be utilized in cases when the nature of the underlying process is not known. © 2009 Wiley Periodicals, Inc. Jrl Fut Mark 30:175–187, 2010     

A non‐lattice pricing model of American options under stochastic volatility

In this article, an analytical approach to American option pricing under stochastic volatility is provided. Under stochastic volatility, the American option value can be computed as the sum of a corresponding European option price and an early exercise premium. By considering the analytical property of the optimal exercise boundary, the formula allows for recursive computation of the American option value. Simulation results show that a nonlattice method performs better than the lattice‐based interpolation methods. The stochastic volatility model is also empirically tested using S&P 500 futures options intraday transactions data. Incorporating stochastic volatility is shown to improve pricing, hedging, and profitability in actual trading. © 2006 Wiley Periodicals, Inc. Jrl Fut Mark 26:417–448, 2006     

PRICING DISCRETELY MONITORED BARRIER OPTIONS AND DEFAULTABLE BONDS IN LÉVY PROCESS MODELS: A FAST HILBERT TRANSFORM APPROACH

This paper presents a novel method to price discretely monitored single- and double-barrier options in Lévy process-based models. The method involves a sequential evaluation of Hilbert transforms of the product of the Fourier transform of the value function at the previous barrier monitoring date and the characteristic function of the (Esscher transformed) Lévy process. A discrete approximation with exponentially decaying errors is developed based on the Whittaker cardinal series (Sinc expansion) in Hardy spaces of functions analytic in a strip. An efficient computational algorithm is developed based on the fast Hilbert transform that, in turn, relies on the FFT-based Toeplitz matrix–vector multiplication. Our method also provides a natural framework for credit risk applications, where the firm value follows an exponential Lévy process and default occurs at the first time the firm value is below the default barrier on one of a discrete set of monitoring dates.     

Valuation Bounds on Barrier Options Under Model Uncertainty

This article investigates valuation bounds on barrier options under model uncertainty. This investigation enriches the literature on the model‐free valuation of these exotic options. It is found that with weak assumptions on underlying price processes, tight valuation bounds on barrier options can be sought from a set of European options. As a result, the numerical routine developed in this article can be reviewed as a new method for the evaluation of barrier options, which is independent of model assumptions. © 2012 Wiley Periodicals, Inc. Jrl Fut Mark 33:199–234, 2013     

A NEW METHOD OF PRICING LOOKBACK OPTIONS

A new method for pricing lookback options (a.k.a. hindsight options) is presented, which simplifies the derivation of analytical formulas for this class of exotics in the Black-Scholes framework. Underlying the method is the observation that a lookback option can be considered as an integrated form of a related barrier option. The integrations with respect to the barrier price are evaluated at the expiry date to derive the payoff of an equivalent portfolio of European-type binary options. The arbitrage-free price of the lookback option can then be evaluated by static replication as the present value of this portfolio. We illustrate the method by deriving expressions for generic, standard floating-, fixed-, and reverse-strike lookbacks, and then show how the method can be used to price the more complex partial-price and partial-time lookback options. The method is in principle applicable to frameworks with alternative asset-price dynamics to the Black-Scholes world.     

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.     

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.     

Step‐reset options: Design and valuation

This study proposes a new design of reset options in which the option's exercise price adjusts gradually, based on the amount of time the underlying spent beyond prespecified reset levels. Relative to standard reset options, a step‐reset design offers several desirable properties. First of all, it demands a lower option premium but preserves the same desirable reset attribute that appeals to market investors. Second, it overcomes the disturbing problem of delta jump as exhibited in standard reset option, and thus greatly reduces the difficulties in risk management for reset option sellers who hedge dynamically. Moreover, the step‐reset feature makes the option more robust against short‐term price movements of the underlying and removes the pressure of price manipulation often associated with standard reset options. To value this innovative option product, we develop a tree‐based valuation algorithm in this study. Specifically, we parameterize the trinomial tree model to correctly account for the discrete nature of reset monitoring. The use of lattice model gives us the flexibility to price step‐reset options with American exercise right. Finally, to accommodate the path‐dependent exercise price, we introduce a state‐to‐state recursive pricing procedure to properly capture the path‐dependent step‐reset effect and enhance computational efficiency. © 2002 John Wiley & Sons, Inc. Jrl Fut Mark 22:155–171, 2002     

Error analysis of finite difference and Markov chain approximations for option pricing

Mijatovi? and Pistorius proposed an efficient Markov chain approximation method for pricing European and barrier options in general one‐dimensional Markovian models. However, sharp convergence rates of this method for realistic financial payoffs, which are nonsmooth, are rarely available. In this paper, we solve this problem for general one‐dimensional diffusion models, which play a fundamental role in financial applications. For such models, the Markov chain approximation method is equivalent to the method of lines using the central difference. Our analysis is based on the spectral representation of the exact solution and the approximate solution. By establishing the convergence rate for the eigenvalues and the eigenfunctions, we obtain sharp convergence rates for the transition density and the price of options with nonsmooth payoffs. In particular, we show that for call‐/put‐type payoffs, convergence is second order, while for digital‐type payoffs, convergence is generally only first order. Furthermore, we provide theoretical justification for two well‐known smoothing techniques that can restore second‐order convergence for digital‐type payoffs and explain oscillations observed in the convergence for options with nonsmooth payoffs. As an extension, we also establish sharp convergence rates for European options for a rich class of Markovian jump models constructed from diffusions via subordination. The theoretical estimates are confirmed using numerical examples.     

PRICING AND SEMIMARTINGALE REPRESENTATIONS OF VULNERABLE CONTINGENT CLAIMS IN REGIME‐SWITCHING MARKETS

Using a suitable change of probability measure, we obtain a Poisson series representation for the arbitrage‐free price process of vulnerable contingent claims in a regime‐switching market driven by an underlying continuous‐time Markov process. As a result of this representation, along with a short‐time asymptotic expansion of the claim's price process, we develop an efficient novel method for pricing claims whose payoffs may depend on the full path of the underlying Markov chain. The proposed approach is applied to price not only simple European claims such as defaultable bonds, but also a new type of path‐dependent claims that we term self‐decomposable, as well as the important class of vulnerable call and put options on a stock. We provide a detailed error analysis and illustrate the accuracy and computational complexity of our method on several market traded instruments, such as defaultable bond prices, barrier options, and vulnerable call options. Using again our Poisson series representation, we show differentiability in time of the predefault price function of European vulnerable claims, which enables us to rigorously deduce Feynman‐Ka? representations for the predefault pricing function and new semimartingale representations for the price process of the vulnerable claim under both risk‐neutral and objective probability measures.     

PRICING CHAINED OPTIONS WITH CURVED BARRIERS

This paper studies barrier options which are chained together, each with payoff contingent on curved barriers. When the underlying asset price hits a primary curved barrier, a secondary barrier option is given to a primary barrier option holder. Then if the asset price hits another curved barrier, a third barrier option is given, and so on. We provide explicit price formulas for these options when two or more barrier options with exponential barriers are chained together. We then extend the results to the options with general curved barriers.     

VALUATION OF CONTINUOUSLY MONITORED DOUBLE BARRIER OPTIONS AND RELATED SECURITIES

In this paper, we apply Carr's randomization approximation and the operator form of the Wiener‐Hopf method to double barrier options in continuous time. Each step in the resulting backward induction algorithm is solved using a simple iterative procedure that reduces the problem of pricing options with two barriers to pricing a sequence of certain perpetual contingent claims with first‐touch single barrier features. This procedure admits a clear financial interpretation that can be formulated in the language of embedded options. Our approach results in a fast and accurate pricing method that can be used in a rather wide class of Lévy‐driven models including Variance Gamma processes, Normal Inverse Gaussian processes, KoBoL processes, CGMY model, and Kuznetsov's β ‐class. Our method can be applied to double barrier options with arbitrary bounded terminal payoff functions, which, in particular, allows us to price knock‐out double barrier put/call options as well as double‐no‐touch options.     

Static hedging and pricing of exotic options with payoff frames

We develop a general framework for statically hedging and pricing European‐style options with nonstandard terminal payoffs, which can be applied to mixed static–dynamic and semistatic hedges for many path‐dependent exotic options including variance swaps and barrier options. The goal is achieved by separating the hedging and pricing problems to obtain replicating strategies. Once prices have been obtained for a set of basis payoffs, the pricing and hedging of financial securities with arbitrary payoff functions is accomplished by computing a set of “hedge coefficients” for that security. This method is particularly well suited for pricing baskets of options simultaneously, and is robust to discontinuities of payoffs. In addition, the method enables a systematic comparison of the value of a payoff (or portfolio) across a set of competing model specifications with implications for security design.     

THE WIENER–HOPF TECHNIQUE AND DISCRETELY MONITORED PATH‐DEPENDENT OPTION PRICING

Fusai, Abrahams, and Sgarra (2006) employed the Wiener–Hopf technique to obtain an exact analytic expression for discretely monitored barrier option prices as the solution to the Black–Scholes partial differential equation. The present work reformulates this in the language of random walks and extends it to price a variety of other discretely monitored path‐dependent options. Analytic arguments familiar in the applied mathematics literature are used to obtain fluctuation identities. This includes casting the famous identities of Baxter and Spitzer in a form convenient to price barrier, first‐touch, and hindsight options. Analyzing random walks killed by two absorbing barriers with a modified Wiener–Hopf technique yields a novel formula for double‐barrier option prices. Continuum limits and continuity correction approximations are considered. Numerically, efficient results are obtained by implementing Padé approximation. A Gaussian Black–Scholes framework is used as a simple model to exemplify the techniques, but the analysis applies to Lévy processes generally.