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.     

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.     

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.     

LÉVY PROCESSES INDUCED BY DIRICHLET (B‐)SPLINES: MODELING MULTIVARIATE ASSET PRICE DYNAMICS

We consider a new class of processes, called LG processes, defined as linear combinations of independent gamma processes. Their distributional and path‐wise properties are explored by following their relation to polynomial and Dirichlet (B‐)splines. In particular, it is shown that the density of an LG process can be expressed in terms of Dirichlet (B‐)splines, introduced independently by Ignatov and Kaishev and Karlin, Micchelli, and Rinott. We further show that the well‐known variance gamma (VG) process, introduced by Madan and Seneta, and the bilateral gamma (BG) process, recently considered by Küchler and Tappe are special cases of an LG process. Following this LG interpretation, we derive new (alternative) expressions for the VG and BG densities and consider their numerical properties. The LG process has two sets of parameters, the B‐spline knots and their multiplicities, and offers further flexibility in controlling the shape of the Levy density, compared to the VG and the BG processes. Such flexibility is often desirable in practice, which makes LG processes interesting for financial and insurance applications. Multivariate LG processes are also introduced and their relation to multivariate Dirichlet and simplex splines is established. Expressions for their joint density, the underlying LG‐copula, the characteristic, moment and cumulant generating functions are given. A method for simulating LG sample paths is also proposed, based on the Dirichlet bridge sampling of gamma processes, due to Kaishev and Dimitriva. A method of moments for estimation of the LG parameters is also developed. Multivariate LG processes are shown to provide a competitive alternative in modeling dependence, compared to the various multivariate generalizations of the VG process, proposed in the literature. Application of multivariate LG processes in modeling the joint dynamics of multiple exchange rates is also considered.     

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.     

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.     

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.