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.     

Stochastic Volatility for Lévy Processes

Three processes reflecting persistence of volatility are initially formulated by evaluating three Lévy processes at a time change given by the integral of a mean-reverting square root process. The model for the mean-reverting time change is then generalized to include non-Gaussian models that are solutions to Ornstein-Uhlenbeck equations driven by one-sided discontinuous Lévy processes permitting correlation with the stock. Positive stock price processes are obtained by exponentiating and mean correcting these processes, or alternatively by stochastically exponentiating these processes. The characteristic functions for the log price can be used to yield option prices via the fast Fourier transform. In general mean-corrected exponentiation performs better than employing the stochastic exponential. It is observed that the mean-corrected exponential model is not a martingale in the filtration in which it is originally defined. This leads us to formulate and investigate the important property of martingale marginals where we seek martingales in altered filtrations consistent with the one-dimensional marginal distributions of the level of the process at each future date.     

Time Changes for Lévy Processes

The goal of this paper is to consider pure jump Lévy processes of finite variation with an infinite arrival rate of jumps as models for the logarithm of asset prices. These processes may be written as time-changed Brownian motion. We exhibit the explicit time change for each of a wide class of Lévy processes and show that the time change is a weighted price move measure of time. Additionally, we present a number of Lévy processes that are analytically tractable, in their characteristic functions and Lévy densities, and hence are relevant for option pricing.     

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.     

VASIČEK BEYOND THE NORMAL

A general Ornstein-Uhlenbeck (OU) process is obtained upon replacing the Brownian motion appearing in the defining stochastic differential equation with a general Lévy process. Certain properties of the Brownian ancestor are distribution-free and carry over to the general OU process. Explicit expressions are obtainable for expected values of a number of functionals of interest also in the general case. Special attention is paid here to gamma- and Poisson-driven OU processes. The Brownian, Poisson, and gamma versions of the OU process are compared in various respects; in particular, their aptitude to describe stochastic interest rates is discussed in view of some standard issues in financial and actuarial mathematics: prices of zero-coupon bonds, moments of present values, and probability distributions of present values of perpetuities. The problem of possible negative interest rates finds its resolution in the general setup by taking the driving Lévy process to be nondecreasing.     

Term Structure Models Driven by General Lévy Processes

As a generalization of the Gaussian Heath–Jarrow–Morton term structure model, we present a new class of bond price models that can be driven by a wide range of Lévy processes. We deduce the forward and short rate processes implied by this model and prove that, under certain assumptions, the short rate is Markovian if and only if the volatility structure has either the Vasicek or the Ho–Lee form. Finally, we compare numerically forward rates and European call option prices in a model driven by a hyperbolic Lévy motion with those in the Gaussian model.     

The Defaultable Lévy Term Structure: Ratings and Restructuring

We introduce the intensity-based defaultable Lévy term structure model. It generalizes the default-free Lévy term structure model by Eberlein and Raible, and the intensity-based defaultable Heath-Jarrow-Morton approach of Bielecki and Rutkowski. Furthermore, we include the concept of multiple defaults, based on Schönbucher, within this generalization.     

VALUATION OF FLOATING RANGE NOTES IN LÉVY TERM-STRUCTURE MODELS

Turnbull (1995) as well as Navatte and Quittard-Pinon (1999) derived explicit pricing formulae for digital options and range notes in a one-factor Gaussian Heath–Jarrow–Morton (henceforth HJM) model. Nunes (2004) extended their results to a multifactor Gaussian HJM framework. In this paper, we generalize these results by providing explicit pricing solutions for digital options and range notes in the multivariate Lévy term-structure model of Eberlein and Raible (1999) , that is, an HJM-type model driven by a d -dimensional (possibly nonhomogeneous) Lévy process. As a byproduct, we obtain a pricing formula for floating range notes in the special case of a multifactor Gaussian HJM model that is simpler than the one provided by Nunes (2004) .     

OPTIMAL CONTINUOUS-TIME HEDGING WITH LEPTOKURTIC RETURNS

We examine the behavior of optimal mean–variance hedging strategies at high rebalancing frequencies in a model where stock prices follow a discretely sampled exponential Lévy process and one hedges a European call option to maturity. Using elementary methods we show that all the attributes of a discretely rebalanced optimal hedge, i.e., the mean value, the hedge ratio, and the expected squared hedging error, converge pointwise in the state space as the rebalancing interval goes to zero. The limiting formulae represent 1-D and 2-D generalized Fourier transforms, which can be evaluated much faster than backward recursion schemes, with the same degree of accuracy. In the special case of a compound Poisson process we demonstrate that the convergence results hold true if instead of using an infinitely divisible distribution from the outset one models log returns by multinomial approximations thereof. This result represents an important extension of Cox, Ross, and Rubinstein to markets with leptokurtic returns.     

MSM Estimators of European Options on Assets with Jumps

This paper shows that, under some regularity conditions, the method of simulated moments estimator of European option pricing models developed by Bossaerts and Hillion (1993) can be extended to the case where the prices of the underlying asset follow Lévy processes, which allow for jumps, with no losses on their asymptotic properties, still allowing for the joint test of the model.     

PSEUDODIFFUSIONS AND QUADRATIC TERM STRUCTURE MODELS

The non-Gaussianity of processes observed in financial markets and the relatively good performance of Gaussian models can be reconciled by replacing the Brownian motion with Lévy processes whose Lévy densities decay as  exp(−λ| x |)  or faster, where  λ > 0  is large. This leads to asymptotic pricing models. The leading term, P ₀, is the price in the Gaussian model with the same instantaneous drift and variance. The first correction term depends on the instantaneous moments of order up to 3, that is, the skewness is taken into account, the next term depends on moments of order 4 (kurtosis) as well, etc. In empirical studies, the asymptotic formula can be applied without explicit specification of the underlying process: it suffices to assume that the instantaneous moments of order greater than 2 are small w.r.t. moments of order 1 and 2, and use empirical data on moments of order up to 3 or 4. As an application, the bond-pricing problem in the non-Gaussian quadratic term structure model is solved. For pricing of options near expiry, a different set of asymptotic formulas is developed; they require more detailed specification of the process, especially of its jump part. The leading terms of these formulas depend on the jump part of the process only, so that they can be used in empirical studies to identify the jump characteristics of the process.     

PUT-CALL SYMMETRY: EXTENSIONS AND APPLICATIONS

Classic put-call symmetry relates the prices of puts and calls at strikes on opposite sides of the forward price. We extend put-call symmetry in several directions. Relaxing the assumptions, we generalize to unified local/stochastic volatility models and time-changed Lévy processes, under a symmetry condition. Further relaxing the assumptions, we generalize to various  asymmetric  dynamics. Extending the conclusions, we take an arbitrarily given payoff of European style or single/double/sequential barrier style, and we construct a conjugate European-style claim of equal value, and thereby a semistatic hedge of the given payoff.     

POWER UTILITY MAXIMIZATION IN CONSTRAINED EXPONENTIAL LÉVY MODELS

We study power utility maximization for exponential Lévy models with portfolio constraints, where utility is obtained from consumption and/or terminal wealth. For convex constraints, an explicit solution in terms of the Lévy triplet is constructed under minimal assumptions by solving the Bellman equation. We use a novel transformation of the model to avoid technical conditions. The consequences for q‐optimal martingale measures are discussed as well as extensions to nonconvex constraints.     

Efficient static replication of European options under exponential Lévy models

This study proposes a new scheme for the static replication of European options and their portfolios. First, a general approximation formula for efficient static replication as an extension of Carr P. and Chou A. (1997, 2002) and Carr P. and Wu L. (2002) is derived. Second, a concrete procedure for implementing the scheme by applying it to plain vanilla options under exponential Lévy models is presented. Finally, numerical examples in a model developed by Carr, P., Geman, H., Madan, D., and Yor M. (2002) are used to demonstrate that the replication scheme is more efficient and more effective in practice than a standard static replication method. © 2008 Wiley Periodicals, Inc. Jrl Fut Mark 29:1–15, 2009     

Explicit Representation of the Minimal Variance Portfolio in Markets Driven by Lévy Processes

In a market driven by a Lévy martingale, we consider a claim ξ. We study the problem of minimal variance hedging and we give an explicit formula for the minimal variance portfolio in terms of Malliavin derivatives. We discuss two types of stochastic (Malliavin) derivatives for ξ: one based on the chaos expansion in terms of iterated integrals with respect to the power jump processes and one based on the chaos expansion in terms of iterated integrals with respect to the Wiener process and the Poisson random measure components. We study the relation between these two expansions, the corresponding two derivatives, and the corresponding versions of the Clark-Haussmann-Ocone theorem.     

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.     

Recursive Formula for Arithmetic Asian Option Prices

I derive a recursive formula for arithmetic Asian option prices with finite observation times in semimartingale models. The method is based on the relationship between the risk‐neutral expectation of the quadratic variation of the return process and European option prices. The computation of arithmetic Asian option prices is straightforward whenever European option prices are available. Applications with numerical results under the Black–Scholes framework and the exponential Lévy model are proposed. © 2012 Wiley Periodicals, Inc. Jrl Fut Mark 34:220–234, 2014     

On approximating deep in‐the‐money Asian options under exponential Lévy processes

This note proposes a new approach of valuing deep in‐the‐money fixed strike and discretely monitoring arithmetic Asian options. This new approach prices Asian options whose underlying asset price evolves according to the exponential of a Lévy process as a weighted sum of European options. Numerical results from experimenting on three different types of Lévy processes—a diffusion process, a jump diffusion process, and a pure jump process—illustrate the accuracy of the approach. © 2010 Wiley Periodicals, Inc. Jrl Fut Mark     

TIME‐CHANGED MARKOV PROCESSES IN UNIFIED CREDIT‐EQUITY MODELING

This paper develops a novel class of hybrid credit‐equity models with state‐dependent jumps, local‐stochastic volatility, and default intensity based on time changes of Markov processes with killing. We model the defaultable stock price process as a time‐changed Markov diffusion process with state‐dependent local volatility and killing rate (default intensity). When the time change is a Lévy subordinator, the stock price process exhibits jumps with state‐dependent Lévy measure. When the time change is a time integral of an activity rate process, the stock price process has local‐stochastic volatility and default intensity. When the time change process is a Lévy subordinator in turn time changed with a time integral of an activity rate process, the stock price process has state‐dependent jumps, local‐stochastic volatility, and default intensity. We develop two analytical approaches to the pricing of credit and equity derivatives in this class of models. The two approaches are based on the Laplace transform inversion and the spectral expansion approach, respectively. If the resolvent (the Laplace transform of the transition semigroup) of the Markov process and the Laplace transform of the time change are both available in closed form, the expectation operator of the time‐changed process is expressed in closed form as a single integral in the complex plane. If the payoff is square integrable, the complex integral is further reduced to a spectral expansion. To illustrate our general framework, we time change the jump‐to‐default extended constant elasticity of variance model of Carr and Linetsky (2006) and obtain a rich class of analytically tractable models with jumps, local‐stochastic volatility, and default intensity. These models can be used to jointly price equity and credit derivatives.     

MARKOWITZ'S PORTFOLIO OPTIMIZATION IN AN INCOMPLETE MARKET

In this paper, for a process S , we establish a duality relation between K_p , the     - closure of the space of claims in     , which are attainable by "simple" strategies, and     , all signed martingale measures     with     , where   p ≥ 1, q ≥ 1  and     . If there exists a     with     a.s., then K_p consists precisely of the random variables     such that ϑ is predictable S -integrable and     for all     . The duality relation corresponding to the case   p = q = 2  is used to investigate the Markowitz's problem of mean–variance portfolio optimization in an incomplete market of semimartingale model via martingale/convex duality method. The duality relationship between the mean–variance efficient portfolios and the variance-optimal signed martingale measure (VSMM) is established. It turns out that the so-called market price of risk is just the standard deviation of the VSMM. An illustrative example of application to a geometric Lévy processes model is also given.