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.     

RATING BASED LÉVY LIBOR MODEL

In this paper, we consider modeling of credit risk within the Libor market models. We extend the classical definition of the default‐free forward Libor rate and develop the rating based Libor market model to cover defaultable bonds with credit ratings. As driving processes for the dynamics of the default‐free and the predefault term structure of Libor rates, time‐inhomogeneous Lévy processes are used. Credit migration is modeled by a conditional Markov chain, whose properties are preserved under different forward Libor measures. Conditions for absence of arbitrage in the model are derived and valuation formulae for some common credit derivatives in this setup are presented.     

CONTINUOUSLY MONITORED BARRIER OPTIONS UNDER MARKOV PROCESSES

In this paper, we present an algorithm for pricing barrier options in one‐dimensional Markov models. The approach rests on the construction of an approximating continuous‐time Markov chain that closely follows the dynamics of the given Markov model. We illustrate the method by implementing it for a range of models, including a local Lévy process and a local volatility jump‐diffusion. We also provide a convergence proof and error estimates for this algorithm.     

MCMC ESTIMATION OF LÉVY JUMP MODELS USING STOCK AND OPTION PRICES

We examine the performances of several popular Lévy jump models and some of the most sophisticated affine jump‐diffusion models in capturing the joint dynamics of stock and option prices. We develop efficient Markov chain Monte Carlo methods for estimating parameters and latent volatility/jump variables of the Lévy jump models using stock and option prices. We show that models with infinite‐activity Lévy jumps in returns significantly outperform affine jump‐diffusion models with compound Poisson jumps in returns and volatility in capturing both the physical and risk‐neutral dynamics of the S&P 500 index. We also find that the variance gamma model of Madan, Carr, and Chang with stochastic volatility has the best performance among all the models we consider.     

HIGH‐ORDER SHORT‐TIME EXPANSIONS FOR ATM OPTION PRICES OF EXPONENTIAL LÉVY MODELS

The short‐time asymptotic behavior of option prices for a variety of models with jumps has received much attention in recent years. In this work, a novel second‐order approximation for at‐the‐money (ATM) option prices is derived for a large class of exponential Lévy models with or without Brownian component. The results hereafter shed new light on the connection between both the volatility of the continuous component and the jump parameters and the behavior of ATM option prices near expiration. In the presence of a Brownian component, the second‐order term, in time‐t, is of the form , with d₂ only depending on Y, the degree of jump activity, on σ, the volatility of the continuous component, and on an additional parameter controlling the intensity of the “small” jumps (regardless of their signs). This extends the well‐known result that the leading first‐order term is . In contrast, under a pure‐jump model, the dependence on Y and on the separate intensities of negative and positive small jumps are already reflected in the leading term, which is of the form . The second‐order term is shown to be of the form and, therefore, its order of decay turns out to be independent of Y. The asymptotic behavior of the corresponding Black–Scholes implied volatilities is also addressed. Our method of proof is based on an integral representation of the option price involving the tail probability of the log‐return process under the share measure and a suitable change of probability measure under which the pure‐jump component of the log‐return process becomes a Y‐stable process. Our approach is sufficiently general to cover a wide class of Lévy processes, which satisfy the latter property and whose Lévy density can be closely approximated by a stable density near the origin. Our numerical results show that the first‐order term typically exhibits rather poor performance and that the second‐order term can significantly improve the approximation's accuracy, particularly in the absence of a Brownian component.     

ROBUST UTILITY MAXIMIZATION WITH LÉVY PROCESSES

We study a robust portfolio optimization problem under model uncertainty for an investor with logarithmic or power utility. The uncertainty is specified by a set of possible Lévy triplets, that is, possible instantaneous drift, volatility, and jump characteristics of the price process. We show that an optimal investment strategy exists and compute it in semi‐closed form. Moreover, we provide a saddle point analysis describing a worst‐case model.     

EFFICIENT PRICING OF BARRIER OPTIONS AND CREDIT DEFAULT SWAPS IN LÉVY MODELS WITH STOCHASTIC INTEREST RATE

Recently, advantages of conformal deformations of the contours of integration in pricing formulas for European options have been demonstrated in the context of wide classes of Lévy models, the Heston model, and other affine models. Similar deformations were used in one‐factor Lévy models to price options with barrier and lookback features and credit default swaps (CDSs). In the present paper, we generalize this approach to models, where the dynamics of the assets is modeled as , where X is a Lévy process, and the interest rate is stochastic. Assuming that X and r are independent, and , the infinitesimal generator of the pricing semigroup in the model for the short rate, satisfies weak regularity conditions, which hold for popular models of the short rate, we develop a variation of the pricing procedure for Lévy models which is almost as fast as in the case of the constant interest rate. Numerical examples show that about 0.15 second suffices to calculate prices of 8 options of same maturity in a two‐factor model with the error tolerance and less; in a three‐factor model, accuracy of order 0.001–0.005 is achieved in about 0.2 second. Similar results are obtained for quanto CDS, where an additional stochastic factor is the exchange rate. We suggest a class of Lévy models with the stochastic interest rate driven by 1–3 factors, which allows for fast calculations. This class can satisfy the current regulatory requirements for banks mandating sufficiently sophisticated credit risk models.     

Double continuation regions for American and Swing options with negative discount rate in Lévy models

In this paper, we study perpetual American call and put options in an exponential Lévy model. We consider a negative effective discount rate that arises in a number of financial applications including stock loans and real options, where the strike price can potentially grow at a higher rate than the original discount factor. We show that in this case a double continuation region arises and we identify the two critical prices. We also generalize this result to multiple stopping problems of Swing type, that is, when successive exercise opportunities are separated by i.i.d. random refraction times. We conduct an extensive numerical analysis for the Black–Scholes model and the jump‐diffusion model with exponentially distributed jumps.     

The valuation of American options in a multidimensional exponential Lévy model

We consider the problem of valuation of American options written on dividend‐paying assets whose price dynamics follow a multidimensional exponential Lévy model. We carefully examine the relation between the option prices, related partial integro‐differential variational inequalities, and reflected backward stochastic differential equations. In particular, we prove regularity results for the value function and obtain the early exercise premium formula for a broad class of payoff functions.     

Lévy betas: Static hedging with index futures

This study considers calibration to forward‐looking betas by extracting information on equity and index options from prices using Lévy models. The resulting calibrated betas are called Lévy betas. The objective of the proposed approach is to capture market expectations for future betas through option prices, as betas estimated from historical data may fail to reflect structural change in the market. By assuming a continuous‐time capital asset pricing model (CAPM) with Lévy processes, we derive an analytical solution to index and stock options, thus permitting the betas to be implied from observed option prices. One application of Lévy betas is to construct a static hedging strategy using index futures. Employing Hong Kong equity and index option data from September 16, 2008 to October 15, 2009, we show empirically that the Lévy betas during the sub‐prime mortgage crisis period were much more volatile than those during the recovery period. We also find evidence to suggest that the Lévy betas improve static hedging performance relative to historical betas and the forward‐looking betas implied by a stochastic volatility model.     

TIME‐CHANGED ORNSTEIN–UHLENBECK PROCESSES AND THEIR APPLICATIONS IN COMMODITY DERIVATIVE MODELS

This paper studies subordinate Ornstein–Uhlenbeck (OU) processes, i.e., OU diffusions time changed by Lévy subordinators. We construct their sample path decomposition, show that they possess mean‐reverting jumps, study their equivalent measure transformations, and the spectral representation of their transition semigroups in terms of Hermite expansions. As an application, we propose a new class of commodity models with mean‐reverting jumps based on subordinate OU processes. Further time changing by the integral of a Cox–Ingersoll–Ross process plus a deterministic function of time, we induce stochastic volatility and time inhomogeneity, such as seasonality, in the models. We obtain analytical solutions for commodity futures options in terms of Hermite expansions. The models are consistent with the initial futures curve, exhibit Samuelson's maturity effect, and are flexible enough to capture a variety of implied volatility smile patterns observed in commodities futures options.     

EXPECTATIONS OF FUNCTIONS OF STOCHASTIC TIME WITH APPLICATION TO CREDIT RISK MODELING

We develop two novel approaches to solving for the Laplace transform of a time‐changed stochastic process. We discard the standard assumption that the background process () is Lévy. Maintaining the assumption that the business clock () and the background process are independent, we develop two different series solutions for the Laplace transform of the time‐changed process . In fact, our methods apply not only to Laplace transforms, but more generically to expectations of smooth functions of random time. We apply the methods to introduce stochastic time change to the standard class of default intensity models of credit risk, and show that stochastic time‐change has a very large effect on the pricing of deep out‐of‐the‐money options on credit default swaps.     

Multiple curve Lévy forward price model allowing for negative interest rates

In this paper, we develop a framework for discretely compounding interest rates that is based on the forward price process approach. This approach has a number of advantages, in particular in the current market environment. Compared to the classical as well as the Lévy Libor market model, it allows in a natural way for negative interest rates and has superb calibration properties even in the presence of extremely low rates. Moreover, the measure changes along the tenor structure are significantly simplified. These properties make it an excellent base for a postcrisis multiple curve setup. Two variants for multiple curve constructions based on the multiplicative spreads are discussed. Time‐inhomogeneous Lévy processes are used as driving processes. An explicit formula for the valuation of caps is derived using Fourier transform techniques. Relying on the valuation formula, we calibrate the two model variants to market data.     

Analytical solvability and exact simulation in models with affine stochastic volatility and Lévy jumps

We investigate analytical solvability of models with affine stochastic volatility (SV) and Lévy jumps by deriving a unified formula for the conditional moment generating function of the log-asset price and providing the condition under which this new formula is explicit. The results lay a foundation for a range of valuation, calibration, and econometric problems. We then combine our theoretical results, the Hilbert transform method, various interpolation techniques, with the dimension reduction technique to propose unified simulation schemes for solvable models with affine SV and Lévy jumps. In contrast to traditional exact simulation methods, our approach is applicable to a broad class of models, maintains good accuracy, and enables efficient pricing of discretely monitored path-dependent derivatives. We analyze various sources of errors arising from the simulation approach and present error bounds. Finally, extensive numerical results demonstrate that our method is highly accurate, efficient, simple to implement, and widely applicable.     

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.     

MULTIVARIATE SUBORDINATION OF MARKOV PROCESSES WITH FINANCIAL APPLICATIONS

This paper develops the procedure of multivariate subordination for a collection of independent Markov processes with killing. Starting from d independent Markov processes with killing and an independent d‐dimensional time change , we construct a new process by time, changing each of the Markov processes with a coordinate . When is a d‐dimensional Lévy subordinator, the time changed process is a time‐homogeneous Markov process with state‐dependent jumps and killing in the product of the state spaces of . The dependence among jumps of its components is governed by the d‐dimensional Lévy measure of the subordinator. When is a d‐dimensional additive subordinator, Y is a time‐inhomogeneous Markov process. When with forming a multivariate Markov process, is a Markov process, where each plays a role of stochastic volatility of . This construction provides a rich modeling architecture for building multivariate models in finance with time‐ and state‐dependent jumps, stochastic volatility, and killing (default). The semigroup theory provides powerful analytical and computational tools for securities pricing in this framework. To illustrate, the paper considers applications to multiname unified credit‐equity models and correlated commodity models.     

DYNAMIC DEFAULTABLE TERM STRUCTURE MODELING BEYOND THE INTENSITY PARADIGM

The two main approaches in credit risk are the structural approach pioneered by Merton and the reduced‐form framework proposed by Jarrow and Turnbull and by Artzner and Delbaen. The goal of this paper is to provide a unified view on both approaches. This is achieved by studying reduced‐form approaches under weak assumptions. In particular, we do not assume the global existence of a default intensity and allow default at fixed or predictable times, such as coupon payment dates, with positive probability. In this generalized framework, we study dynamic term structures prone to default risk following the forward‐rate approach proposed by Heath, Jarrow, and Morton. It turns out that previously considered models lead to arbitrage possibilities when default can happen at a predictable time. A suitable generalization of the forward‐rate approach contains an additional stochastic integral with atoms at predictable times and necessary and sufficient conditions for an appropriate no‐arbitrage condition are given. For efficient implementations, we develop a new class of affine models that do not satisfy the standard assumption of stochastic continuity. The chosen approach is intimately related to the theory of enlargement of filtrations, for which we provide an example by means of filtering theory where the Azéma supermartingale contains upward and downward jumps, both at predictable and totally inaccessible stopping times.     

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.     

LOCAL WELL‐POSEDNESS OF MUSIELA’S SPDE WITH LÉVY NOISE

We determine sufficient conditions on the volatility coefficient of Musiela’s stochastic partial differential equation driven by an infinite dimensional Lévy process so that it admits a unique local mild solution in spaces of functions whose first derivative is square integrable with respect to a weight.     

AN EXACT FORMULA FOR DEFAULT SWAPTIONS’ PRICING IN THE SSRJD STOCHASTIC INTENSITY MODEL

We develop and test a fast and accurate semi‐analytical formula for single‐name default swaptions in the context of a shifted square root jump diffusion (SSRJD) default intensity model. The model can be calibrated to the CDS term structure and a few default swaptions, to price and hedge other credit derivatives consistently. We show with numerical experiments that the model implies plausible volatility smiles.