SHORT-TIME IMPLIED VOLATILITY IN EXPONENTIAL LÉVY MODELS

Similar documents
M5MF6. Advanced Methods in Derivatives Pricing

Local vs Non-local Forward Equations for Option Pricing

Near-Expiry Asymptotics of the Implied Volatility in Local and Stochastic Volatility Models

Short-time asymptotics for ATM option prices under tempered stable processes

Near-expiration behavior of implied volatility for exponential Lévy models

The stochastic calculus

Pricing of some exotic options with N IG-Lévy input

UNIFORM BOUNDS FOR BLACK SCHOLES IMPLIED VOLATILITY

Large Deviations and Stochastic Volatility with Jumps: Asymptotic Implied Volatility for Affine Models

Small-time asymptotics of stopped Lévy bridges and simulation schemes with controlled bias

Financial Derivatives Section 5

Conditional Density Method in the Computation of the Delta with Application to Power Market

SADDLEPOINT APPROXIMATIONS TO OPTION PRICES 1. By L. C. G. Rogers and O. Zane University of Bath and First Chicago NBD

CRRAO Advanced Institute of Mathematics, Statistics and Computer Science (AIMSCS) Research Report. B. L. S. Prakasa Rao

Optimally Thresholded Realized Power Variations for Lévy Jump Diffusion Models

Lecture Note 8 of Bus 41202, Spring 2017: Stochastic Diffusion Equation & Option Pricing

1 Geometric Brownian motion

Estimation of Value at Risk and ruin probability for diffusion processes with jumps

The Black-Scholes Model

Risk Minimization Control for Beating the Market Strategies

MATH3075/3975 FINANCIAL MATHEMATICS TUTORIAL PROBLEMS

Tangent Lévy Models. Sergey Nadtochiy (joint work with René Carmona) Oxford-Man Institute of Quantitative Finance University of Oxford.

Lecture 4. Finite difference and finite element methods

Pricing Dynamic Solvency Insurance and Investment Fund Protection

Extended Libor Models and Their Calibration

MSc Financial Engineering CHRISTMAS ASSIGNMENT: MERTON S JUMP-DIFFUSION MODEL. To be handed in by monday January 28, 2013

Rohini Kumar. Statistics and Applied Probability, UCSB (Joint work with J. Feng and J.-P. Fouque)

1 Implied Volatility from Local Volatility

NEWCASTLE UNIVERSITY SCHOOL OF MATHEMATICS, STATISTICS & PHYSICS SEMESTER 1 SPECIMEN 2 MAS3904. Stochastic Financial Modelling. Time allowed: 2 hours

How do Variance Swaps Shape the Smile?

The Black-Scholes Model

Asset Pricing Models with Underlying Time-varying Lévy Processes

Option pricing in the stochastic volatility model of Barndorff-Nielsen and Shephard

A Continuity Correction under Jump-Diffusion Models with Applications in Finance

Universität Regensburg Mathematik

Homework Assignments

American Foreign Exchange Options and some Continuity Estimates of the Optimal Exercise Boundary with respect to Volatility

Risk Neutral Measures

Applications of short-time asymptotics to the statistical estimation and option pricing of Lévy-driven models

MASM006 UNIVERSITY OF EXETER SCHOOL OF ENGINEERING, COMPUTER SCIENCE AND MATHEMATICS MATHEMATICAL SCIENCES FINANCIAL MATHEMATICS.

The value of foresight

Using Lévy Processes to Model Return Innovations

Recovering portfolio default intensities implied by CDO quotes. Rama CONT & Andreea MINCA. March 1, Premia 14

Brownian Motion. Richard Lockhart. Simon Fraser University. STAT 870 Summer 2011

( ) since this is the benefit of buying the asset at the strike price rather

Lecture 11: Ito Calculus. Tuesday, October 23, 12

1 The continuous time limit

Change of Measure (Cameron-Martin-Girsanov Theorem)

Journal of Mathematical Analysis and Applications

Stochastic Processes and Stochastic Calculus - 9 Complete and Incomplete Market Models

The ruin probabilities of a multidimensional perturbed risk model

Introduction to Probability Theory and Stochastic Processes for Finance Lecture Notes

Optimal stopping problems for a Brownian motion with a disorder on a finite interval

1.1 Basic Financial Derivatives: Forward Contracts and Options

Chapter 3: Black-Scholes Equation and Its Numerical Evaluation

A note on the existence of unique equivalent martingale measures in a Markovian setting

Time-changed Brownian motion and option pricing

Financial Engineering. Craig Pirrong Spring, 2006

PRICING TIMER OPTIONS UNDER FAST MEAN-REVERTING STOCHASTIC VOLATILITY

Asymptotic methods in risk management. Advances in Financial Mathematics

NUMERICAL METHODS OF PARTIAL INTEGRO-DIFFERENTIAL EQUATIONS FOR OPTION PRICE

Logarithmic derivatives of densities for jump processes

OPTION PRICE WHEN THE STOCK IS A SEMIMARTINGALE

Deterministic Income under a Stochastic Interest Rate

Short-Time Asymptotic Methods In Financial Mathematics

Heston vs Heston. Antoine Jacquier. Department of Mathematics, Imperial College London. ICASQF, Cartagena, Colombia, June 2016

Uniform Bounds for Black Scholes Implied Volatility

Chapter 15: Jump Processes and Incomplete Markets. 1 Jumps as One Explanation of Incomplete Markets

AMH4 - ADVANCED OPTION PRICING. Contents

Math 6810 (Probability) Fall Lecture notes

Two-sided estimates for stock price distribution densities in jump-diffusion models

On the Lower Arbitrage Bound of American Contingent Claims

EFFICIENT MONTE CARLO ALGORITHM FOR PRICING BARRIER OPTIONS

How to hedge Asian options in fractional Black-Scholes model

Sensitivity of American Option Prices with Different Strikes, Maturities and Volatilities

The Use of Importance Sampling to Speed Up Stochastic Volatility Simulations

Weierstrass Institute for Applied Analysis and Stochastics Maximum likelihood estimation for jump diffusions

Applications of Lévy processes

Lecture 1: Lévy processes

Hedging Credit Derivatives in Intensity Based Models

Lecture 17. The model is parametrized by the time period, δt, and three fixed constant parameters, v, σ and the riskless rate r.

Dr. Maddah ENMG 625 Financial Eng g II 10/16/06

Advanced Topics in Derivative Pricing Models. Topic 4 - Variance products and volatility derivatives

BROWNIAN MOTION Antonella Basso, Martina Nardon

VALUATION OF FLEXIBLE INSURANCE CONTRACTS

Are stylized facts irrelevant in option-pricing?

Calculation of Volatility in a Jump-Diffusion Model

Pricing in markets modeled by general processes with independent increments

Lecture 15: Exotic Options: Barriers

Exponential utility maximization under partial information

AN ANALYTICALLY TRACTABLE UNCERTAIN VOLATILITY MODEL

arxiv: v1 [q-fin.pr] 18 Feb 2010

1 Mathematics in a Pill 1.1 PROBABILITY SPACE AND RANDOM VARIABLES. A probability triple P consists of the following components:

Asymptotic results discrete time martingales and stochastic algorithms

Economathematics. Problem Sheet 1. Zbigniew Palmowski. Ws 2 dw s = 1 t

Optimal Dividend Policy of A Large Insurance Company with Solvency Constraints. Zongxia Liang

SPDE and portfolio choice (joint work with M. Musiela) Princeton University. Thaleia Zariphopoulou The University of Texas at Austin

Non-semimartingales in finance

REFINED WING ASYMPTOTICS FOR THE MERTON AND KOU JUMP DIFFUSION MODELS

Managing Systematic Mortality Risk in Life Annuities: An Application of Longevity Derivatives

Transcription:

SHORT-TIME IMPLIED VOLATILITY IN EXPONENTIAL LÉVY MODELS ERIK EKSTRÖM1 AND BING LU Abstract. We show that a necessary and sufficient condition for the explosion of implied volatility near expiry in exponential Lévy models is the existence of jumps towards the strike price in the underlying process. When such jumps do not exist, the implied volatility converges to the volatility of the Gaussian component of the underlying Lévy process as the time to maturity tends to zero. These results are proved by comparing the short-time asymptotics of the Black-Scholes price with explicit formulas for upper and lower bounds of option prices in exponential Lévy models. 1. Introduction In financial markets, prices of European options are usually quoted in terms of their implied volatilities, which are computed through inversion of the classical Black-Scholes formula. An advantage of expressing prices in volatility is that it gives an easy comparison between options with different features. However, under a given modelling framework for the underlying asset price, there is little hope to obtain a closed-form expression of the implied volatility. Instead, a good deal of attention has been directed towards the asymptotic behavior of the volatility surface. Recently, the understanding of the short-maturity asymptotics of implied volatility has advanced, see [3] and [12] for results in local volatility models, [8], [9] and[10] for stochastic volatility models, and [1], [7], [14], [15] and [17] for exponential Lévy models. In this paper we exclusively study the short-term behavior of implied volatility in exponential Lévy models. For a nice and thorough introduction to exponential Lévy processes we refer to [5]. The small-time behavior for at-the-money options is well-understood, see [14], [15] and [17]. These references show that the implied volatility converges to the volatility of the Gaussian component of the underlying Lévy process. Moreover, in [14] it is shown how to change variables in order to have a non-trivial limit when zooming in at-the-money and close to expiry. On the other hand, the precise asymptotic behavior for out-of-the-money and in-the-money options is not clear for general exponential Lévy models. In fact, existing literature typically focuses on cases in which the implied volatility explodes close to Key words and phrases. implied volatility; exponential Lévy models; short-time asymptotic behavior. 1 Support from the Swedish Research Council (VR) is gratefully acknowledged. 1

2 ERIK EKSTRÖM AND BING LU expiry, and then derives asymptotic expansions of implied volatility. In [7] and [17], it is assumed that the Lévy measure has full support, meaning that the stock price can jump to any positive value. The condition in [15] is slightly weaker, as it merely requires a positive probability that the stock price takes a single leap across the strike level. However, the exact condition on the model to have explosion of the short-time implied volatility seems to be lacking. Moreover, very little is known about the case in which the implied volatility does not blow up for not at-the-money options. A simple example is given in [15] though, where the implied volatility converges to zero for an in-the-money option in a spectrally positive exponential Lévy model with no Gaussian component and with additional restrictions on the jump activity. The key contribution of the current paper is to give a complete picture of the short-time implied volatility behavior in a general exponential Lévy model. Our findings show that the implied volatility explodes in the shorttime limit in all exponential Lévy models except in three cases, in which it converges to the volatility of the Gaussian component of the underlying Lévy process. Denoting the strike price by K and the stock price by s, these three cases are: (i) s < K in spectrally negative exponential Lévy models; (ii) s > K in spectrally positive exponential Lévy models; (iii) s = K. In particular, if the underlying Lévy process has no Gaussian part, the implied volatility in each of these three cases converges to zero. Our results imply that neither the Gaussian component nor the absolute value of the jump sizes in the underlying Lévy process decides whether the short-time implied volatility explodes or converges. In fact, the existence of jumps towards the strike price in the underlying process is the sole cause of an implied volatility explosion. Note that even if the jump size is bounded so that it takes more than one jump to reach the strike, implied volatility explosion occurs. To prove the implied volatility explosion in cases when the underlying process have jumps towards the strike price, we use the convexity of option price in the underlying stock value to find a lower bound of the option price by removing the Brownian fluctuations as well as small jumps from the original Lévy process. For the resulting finite-activity model, the option price asymptotics can be derived explicitly, and the explosion of short-time implied volatility can thus be deduced easily by comparing the lower bound and the asymptotics of the Black-Scholes formula. To prove that the short-term implied volatility converges to the volatility of the diffusion component in the remaining cases is technically more demanding, and we deal with it in several steps. First we consider the case s < K for models of spectrally negative Lévy type with finite activity. Another application of the convexity property shows that the call price is

SHORT-TIME IMPLIED VOLATILITY IN EXPONENTIAL LÉVY MODELS 3 bounded above by a call price in a jump-to-default model, which can be computed explicitly using the Black-Scholes formula. Using this explicit upper bound, the short-term implied volatility is shown to approach the volatility of the underlying. In the next step, the class of models is extended to models with an infinite jump intensity. Here we construct a new Lévy process by replacing the infinitely many small jumps with a Brownian motion with twice the variance. Using a maximum principle for parabolic integro-differential equations, we show that the call price with the new Lévy process as underlying is greater than the call price with the original Lévy process as underlying near expiry for s < K, which, together with the result of the previous step, implies the result. Finally, in the last step we show that models of spectrally positive Lévy type can be reduced to the spectrally negative case by means of a measure change. The current paper focuses on giving a full picture of the limiting behavior of implied volatility in a general exponential Lévy model, rather than deriving asymptotic formulas in special models. However, our proofs build on various lower and upper bounds of the option price, and these bounds are of course helpful also in studies of asymptotic expansions of the implied volatility. The paper is organized as follows. In Section 2 we introduce the model and state our main results. In Section 3 we present some well-known facts on option prices that are used later. In Section 4 we deal with the cases when explosion of the short-time implied volatility occurs. Finally, Section 5 treats the cases in which the short-term implied volatility converges to the volatility of the Gaussian component of the underlying Lévy process. 2. Exponential Lévy models and the main result We consider the class of exponential Lévy processes for a stock price process S t = se Xt, where s > 0 is a constant and X = (X t ) t 0 is a Lévy process with Lévy- Khinchin representation E[exp{izX t }] = exp{tψ(z)}, z R, where ψ(z) = σ2 z 2 + iγz + (e izy 1 izy1 2 { y 1} )ν(dy). R Here the triplet (γ, σ 2, ν) is referred to as the characteristic triplet of X, σ 0 is the volatility of the Gaussian component and ν is the Lévy measure of the underlying Lévy process. We assume that the model is specified directly under the measure which is used for derivative pricing, and that the Lévy measure is such that the exponential moment E[e Xt ] is finite. More precisely, we assume that e y ν(dy) <, y >1

4 ERIK EKSTRÖM AND BING LU and that γ = γ(σ, ν) = σ2 2 (e y 1 y1 y 1 )ν(dy). R Under these assumptions, the stock price (S t ) t 0 is a martingale, and the price of a European call option is then defined by (1) C(t, s, K) = E[(S t K) + ] = E[(se Xt K) + ], where K is the strike price and t is the time to maturity. Remark 1. We emphasize that the time parameter t in the option pricing formula C(t, s, K) denotes time left to maturity, whereas the Lévy process X always starts at 0 at time 0, so se Xt denotes the position of the process after an elapsed time t. By the time-homogeneity of X, this is equivalent to starting X at 0 at time T t, and then letting it evolve until T, where T denotes the (fixed) horizon. The classical Black-Scholes formula for European call option prices in a geometric Brownian motion model is ( ln (2) C BS s K (t, s, K, σ) = sn + ) ( 1 2 σ2 t ln s σ K KN ) 1 2 σ2 t t σ t for (t, s, K, σ) (R + ) 4, where N( ) is the standard normal cumulative distribution function. With all other parameters fixed, the Black-Scholes price is a strictly increasing function of the volatility σ. Thus, for s > 0 fixed, the Black-Scholes implied volatility function ˆσ(t, K) : R + R + [0, ) of an exponential Lévy model can be defined through C(t, s, K) = C BS (t, s, K, ˆσ(t, K)). The main result of the paper is stated in the following theorem. Theorem 1. The implied volatility in an exponential Lévy model satisfies lim ˆσ(t, K) = unless (i) s < K in spectrally negative exponential Lévy models, (ii) s > K in spectrally positive exponential Lévy models, or (iii) s = K. In cases (i)-(iii), we have (3) lim ˆσ(t, K) = σ, i.e. the implied volatility converges to the volatility of the Gaussian component of the underlying Lévy process. Note that the convergence result (3) for case (iii) has been proved in [15] (see also [17], where it is proved under the additional assumption x 2 ν(dx) < ).

SHORT-TIME IMPLIED VOLATILITY IN EXPONENTIAL LÉVY MODELS 5 3. Ordering and asymptotic properties In this section we state some well-known results about the option price and implied volatility. Our first proposition gives an ordering result for models with similar characteristic triplets. Proposition 1. Assume that X i, i = 1, 2 are two Lévy processes with characteristic triplets (γ(σ 1, ν), σ 1, ν) and (γ(σ 2, ν), σ 2, ν), respectively, so that e Xi, i = 1, 2 are both martingales with the same Lévy measure ν, but with different volatilities of the Gaussian components. Let ˆσ i (t, K), i = 1, 2, be the corresponding implied volatilities. If 0 σ 1 σ 2, then (4) σ 1 ˆσ 1 (t, K) ˆσ 2 (t, K). Proof. Let X 1 be a Lévy process with a diffusion component with volatility σ 1, and let W t be a Brownian motion, independent of X 1. Then S 2 = se X2 can be represented as S 2 t = se X1 t M t, where M is the exponential martingale We then have M t = e σ 2 2 σ2 1 t+ σ 2 2 2 σ2 1 Wt. C 2 (t, s, K) = E[(se X1 t M t K) + ] = E[E[(se X1 t M t K) + X 1 t ]] E[(se X1 t K) + ] = C 1 (t, s, K), where the inequality is obtained from Jensen s inequality for conditional expectations. Consequently, ˆσ 2 (t, K) ˆσ 1 (t, K). A similar argument, in which one decomposes X 1 in a Brownian part and a jump part, gives ˆσ 1 (t, K) σ 1. Our next proposition describes the asymptotics of the option price for not at-the-money options. For two functions f and g we write f g if f(t) g(t) 1 as t 0. Proposition 2. As t 0, the Black-Scholes price defined in (2) has the asymptotics } C BS (t, s, K, σ) (s K) + C 1 σ 3 t 3/2 (ln K ln s)2 (5) exp { 2tσ 2 for s K, where C 1 = sk 2π(ln K ln s) 2. Moreover, for a given exponential Lévy model, the option price C satisfies } C(t, s, K) (s K) + C 1ˆσ 3 (t, K)t 3/2 (ln K ln s)2 (6) exp { 2tˆσ 2 (t, K) for s K.

6 ERIK EKSTRÖM AND BING LU Proof. The asymptotics (5) are well-known, see for example [12]. To prove (6), note that C(t, s, K) = E[(K S t ) + ] + s K (K s) + + s K = (s K) + as t 0 by bounded convergence. Since C BS (t, s, K, σ) only depends on t and σ through the quantity σ 2 t, this implies that ˆσ 2 (t, K)t 0 as t 0. The asymptotics (6) therefore follow from (5). 4. Implied volatility explosion This section treats the case in which the implied volatility explodes as time to maturity vanishes. To prove the following theorem, we determine an explicit lower bound of the option price by removing the Brownian motion and certain jumps using the convexity of the payoff. Theorem 2. The implied volatility tends to infinity as time to expiry tends to zero in all exponential Lévy models that do not belong to one of the cases (i), (ii) or (iii) of Theorem 1. Proof. For s > K, consider an exponential Lévy model which is not of spectrally positive type. We decompose the underlying Lévy process as the sum of two independent Lévy processes, where one of them has no Brownian fluctuations and only large negative jumps, and the other one contains the continuous fluctuations as well as positive jumps and small negative jumps. To be precise, given a Lévy process X which is not spectrally positive and with characteristic triplet (γ, σ 2, ν), take ε > 0 such that λ := ν(dy) > 0, and decompose X as where and X t = X 1 t + X 2 t, Xt 1 = t (e y 1)ν(dy) + yn(t, dy) Xt 2 = σ2 2 t + σw t t (e y 1 y1 y 1 )ν(dy) y ε + yñ(t, dy) + yn(t, dy). ε y 1 Here W is a standard Brownian motion, N is a Poisson random measure with intensity ν(dy)dt, Ñ is the compensated version of N. Note that X 1 is a compound Poisson process with drift and with jump intensity λ, and that e X1 and e X2 are independent martingales. Let β = 1 e y ν(dy) λ y>1

SHORT-TIME IMPLIED VOLATILITY IN EXPONENTIAL LÉVY MODELS 7 so that 1 β is the expected relative jump size of e X1, let and let Nt 1 = µ = λ(1 β) count the number of jumps of X 1. Then (7) N(t, dy) E[(se Xt K) + ] = E[E[(se X1 t +X2 t K) + N 1 t ]] E[(sE[e X1 t Nt 1 ] K) + ] λt (λt)n = e (se µt β n K) + n! n=0 λt (λt)n e (K se µt β n ) + s K, n! where the first inequality follows from Jensen s inequality for conditional expectations, and in the last inequality we chose the unique integer n 1 such that sβ n < K and sβ n 1 K. By (7) and Proposition 2, we have (8) lim sup lim sup e λt (λt) n (K se µt β n ) { C 1 n!ˆσ 3 (t, K)t 3/2 exp (ln K ln s)2 2tˆσ 2 (t,k) E[(se Xt K) + (s K))] { C 1ˆσ 3 (t, K)t 3/2 exp (ln K ln s)2 2ˆσ(t,K) 2 t } } = 1, which implies that ˆσ(t, K) as t 0. Finally, a change of measure technique (see e.g. [13]) can be used to reduce the case when s < K and X is a Lévy model which is not spectrally negative to the previous case. Indeed, assume that X has positive jumps and satisfies E[e Xt ] = 1. Define a probability measure Q by dq dp := ext. Note that X is a Lévy process which has negative jumps, and (e Xs ) s [0,t] is a Q-martingale. Moreover, E[(se Xt K) + ] = E Q [(s Ke Xt ) + ] = E Q [(Ke Xt s) + ] + K s, where the second equality is the put-call-parity. Now, since E Q [(Ke Xt s) + ] is a call price in an exponential Lévy model with negative jumps, with stock price K and strike price s < K, and since the Black-Scholes price satisfies C BS (t, K, s, σ) + K s = C BS (t, s, K, σ), implied volatility explosion follows from the first part of the proof.

8 ERIK EKSTRÖM AND BING LU Remark 2. (Explosion rates.) To conclude that the short-term implied volatility explodes and to derive its explosion rate, existing literature uses the small-time asymptotics C(t, s, K) (s K) + (9) lim inf η > 0 t 0 t for the option price under certain conditions (e.g. full support of Lévy measure). The proof of Theorem 2 utilizes the fact that if s > K and if the underlying Lévy measure supports jumps towards the strike price, then the option price satisfies the weaker asymptotics C(t, s, K) (s K) + (10) lim inf t 0 t n η for some integer n 1 and some constant η > 0. However, as is shown in the proof above, (10) also leads to implied volatility explosion. Moreover, it is straightforward to find a lower bound for the explosion rate. Indeed, following the proof of [17, Proposition 4] gives lim inf 2nˆσ 2 (t, K)t ln 1 t (ln s K )2 1. Remark 3. (Spectrally negative models do not necessarily lead to skews.) Consider the spectrally negative exponential Lévy model in which S t = se λβt β Nt, where N is a Poisson process with intensity λ and the relative jump size 1 β (0, 1) is fixed. In this model, as t 0, the implied volatility satisfies 2nˆσ 2 (t, K)t ln 1 t (ln s K )2 for s > K, where n 1 is the unique integer such that sβ n < K and sβ n 1 K. Thus the rate of the implied volatility explosion for K just below ˆK := sβ (for which n = 2) is strictly smaller than the explosion rate for K just above ˆK (for which n = 1). An interesting consequence of this is that spectrally negative models not necessarily lead to volatility skews (compare the statement that the presence of a skew is attributed to the fear of large negative jumps, [5, page 10]). 5. Convergence of implied volatility In this section we show that the implied volatility converges to the volatility of the Gaussian component of the underlying Lévy process as the time to maturity goes to zero in cases (i) and (ii) of Theorem 1 (case (iii) is proved in [15]). The convergence results are proved in three steps. The first step treats models of spectrally negative exponential Lévy type with finite activity for s < K. In the second step we examine the case of spectrally

SHORT-TIME IMPLIED VOLATILITY IN EXPONENTIAL LÉVY MODELS 9 negative exponential Lévy processes with infinite activity for s < K. Finally, in the third step we study the case of spectrally positive exponential Lévy processes for s > K. Step 1. This step treats case (i) of Theorem 1 under the additional assumption that the underlying Lévy process has finite jump intensity. Theorem 3. In spectrally negative exponential Lévy models with finite activity, i.e. when the Lévy measure satisfies ν(0, ) = 0 and ν(, 0) <, the implied volatility satisfies for s < K. lim ˆσ(t, K) = σ Proof. First note that it follows from Proposition 1 that it suffices to consider the case σ > 0. Since the underlying spectrally negative Lévy process has finite activity, its discontinuous part is just a compound Poisson process with negative jumps. Let λ := ν(, 0) be the intensity of the compound Poisson process, denote by F (dx) = 1 λν(dx) the distribution of the jump size, let β = x<0 ex F (dx) so that 1 β is the expected relative jump size, and let µ = λ(1 β). Write C(t, s) for C(t, s, K), and note that since σ > 0, C is a classical solution of C t = 1 2 σ2 s 2 C ss + µsc s + λ C(t, e x s)f (dx) λc, see [6]. By convexity of s C(t, s) and since C(t, 0) = 0, we have C(t, e x s)f (dx) e x CF (dx) = βc, so x<0 x<0 x<0 (11) C t 1 2 σ2 s 2 C ss + µsc s + λβc λc. Let C(t, s) be the Black-Scholes option price with interest rate µ, i.e. so that C(t, s) = C BS (t, s, Ke µt, σ), (12) C t = 1 2 σ2 s 2 C ss + µsc s µc. Set V (t, s) = C(t, s) C(t, s), and note that V t 1 2 σ2 s 2 V ss + µsv s µv for t > 0 and s > 0 by (11) and (12). Since V (0, s) = 0, V (t, 0) = 0 and K V (t, s) K for all t 0 and s 0, it follows from the maximum

10 ERIK EKSTRÖM AND BING LU principle (see [11, Theorem 9, Chapter 2]) that C C. By the definition of implied volatility we get that C BS (t, s, K, ˆσ(t, K)) (13) lim sup C(t, s) lim sup C BS (t, s, K, ˆσ(t, K)) C(t, s) = 1. On the other hand, it follows from (6) and the Black-Scholes asymptotics with positive interest rate (see [12, Lemma 2.5]) that (14) C BS (t, s, K, ˆσ(t, K)) C(t, s) ˆσ3 (t, K) σ 3 exp { µ ln K s σ 2 + (ln K s )2 2t ( ) } 1 σ 2 1 ˆσ 2 (t, K) as t 0. Combining (13), (14) and ˆσ(t, K) σ from Proposition 1 we see that ˆσ(t, K) σ as t 0, which finishes the proof. Step 2. In this step we show convergence of the short-time implied volatility when s < K in spectrally negative exponential Lévy models with infinite jump activity. Given such a model, we construct a new model by replacing the small jumps by a Brownian motion with twice the variance. This replacement technique is often used for simulation of Lévy processes, see [2], [16] and [17]. To be precise, let X be a spectrally negative Lévy process X with Lévy triplet (γ, σ 2, ν). As in Step 1, Proposition 1 yields that it suffices to consider the case σ > 0. Also, we assume that ν(( ε, 0)) > 0 for any ε > 0 (otherwise we are done by Step 1). Given ε > 0, let X ε be a Lévy process characterized by the triplet (γ ε, σε, 2 ν ε ), where σε 2 = σ 2 + 2 y 2 ν(dy), and Then e Xε ε y<0 ν ε (dy) = 1 ν(dy) γ ε = σ2 ε 2 (e y 1 y1 y 1 )ν ε (dy). R is a martingale, and σ ε > σ. For fixed K > 0, the price C ε (t, s) = E[(se Xε t K) + ] of a call option with X ε as underlying is then a classical solution of Ct ε = 1 (15) 2 σ2 εs 2 Css ε + (C ε (t, se y ) C ε s(e y 1)Cs) ε ν(dy), see [6]. Set Ω := {(t, s) : 0 < t ln K/s, 0 < s < K}, d ε

SHORT-TIME IMPLIED VOLATILITY IN EXPONENTIAL LÉVY MODELS 11 where d ε = 3 2 σ2 ε (e y 1)ν(dy). Our first result states that C ε is S-superharmonic in Ω. Lemma 1. The function C ε satisfies on Ω. Ct ε 1 2 σ2 s 2 Css ε + (C ε (t, se y ) C ε s(e y 1)Cs)ν(dy) ε y<0 Proof. Since X ε has finite jump intensity, it can be decomposed as X ε t = σ ε W t Y t + γt, where Y 0 is a compound Poisson process with only positive jumps and Denoting the density of e Xε t γ = σ2 ε 2 (e y 1)ν(dy). by f(x), we have f(x) = [ ( x E P W t ln x + Y )] t γt σ ε Y t [ = ln x+yt γt ] x E σε t ϕ(y) dy [ ( )] 1 ln x + Yt γt = E ϕ, xσ ε t σ ε t where ϕ denotes the density of the standard normal distribution. Thus differentiating C ε (t, s) = R (sx K)+ f(x) dx three times with respect to s gives C ε sss(t, s) = [ ( ) ( )] K ln(k/s) γt + Yt ln(k/s) γt + Yt E ϕ s 3 σ ε t σ ε t σεt 2 2 [ ( ) ( )] K ln(k/s) γt + Yt ln(k/s) γt E ϕ s 3 σ ε t σ ε t σεt 2 2. Consequently, C ε sss 0 in Ω.

12 ERIK EKSTRÖM AND BING LU We then have (16) (C ε (t, se y ) C ε s(e y 1)Cs)ν(dy) ε s(1 e y )(Cs ε Cs(t, ε se y ))ν(dy) ε y<0 s 2 (1 e y ) 2 Cssν(dy) ε ε y<0 s 2 Css ε y 2 ν(dy) ε y<0 ε y<0 = σ2 ε σ 2 s 2 Css ε 2 for (t, s) Ω, where the first inequality follows from C ss 0 and the second one from Csss ε 0 in Ω. Combining (15) and (16) we obtain the lemma. Next we show that C ε (t, s) C(t, s) in Ω τ := Ω {t τ} for some τ > 0. Lemma 2. For a fixed K, there exists τ > 0 such that the option prices satisfy C ε (t, s) C(t, s), for (t, s) Ω τ. Proof. We have C ε (t, Ke dεt ) C(t, Ke dεt ) = E[(exp(Xε t d ε t)k K) + ] E[(exp(X t d ε t)k K) + ] E[(exp(X ε t d ε t) 1) + ] E[(exp(σ ε W t 1 2 σ2 εt) 1) + ] E[(exp(σ ε W t 1 2 σ2 εt) 1) + ] E[(e Xt 1) +. ] The first factor on the right hand side converges to 1 as t 0 by [15, Lemma 2.7]. Using the asymptotics for at-the-money options from [15, Theorem 3.1], the second factor converges to σ ε /σ > 1 as t 0. Hence there exists τ > 0 such that C ε (t, Ke dεt ) C(t, Ke dεt ) for t τ. Let V (t, s) = C ε (t, s) C(t, s). Lemma 1 above and [6, Proposition 2] yield that (17) V t 1 2 σ2 s 2 V ss + (V (t, se y ) V + (1 e y )sv s )ν(dy) y<0 in Ω τ. On the boundaries V satisfies V (t, 0) = V (0, s) = 0 and V (t, Ke dεt ) 0 for t τ. Since X is spectrally negative, the non-local term in (17) only depends on the value of V within Ω τ. Consequently, a maximum principle similar to the one in [4] can be applied, which yields V (t, s) 0 for (t, s) Ω τ.

SHORT-TIME IMPLIED VOLATILITY IN EXPONENTIAL LÉVY MODELS 13 We are now ready to conclude that the implied volatility converges to σ in case (i) in Theorem 1. Theorem 4. In spectrally negative exponential Lévy models, the implied volatility satisfies (18) for s < K. lim ˆσ(t, K) = σ Proof. Denote by ˆσ ε (t, K) and ˆσ(t, K) the implied volatilities in the Lévy models X ε and X, respectively. Proposition 1, Lemma 2 and Theorem 3 imply that σ lim inf ˆσ(t, K) lim sup ˆσ(t, K) lim sup ˆσ ε (t, K) = σ ε. Since ε > 0 was arbitrary, and since lim ε 0 σ 2 ε = σ 2, the theorem follows. Step 3. In this step we show that the convergence of the short-term implied volatility in case (ii) of Theorem 1 follows from Theorem 4. Theorem 5. In spectrally positive exponential Lévy models, the short-term implied volatility satisfies (19) lim ˆσ(t, K) = σ for s > K. t 0 Proof. Assume that s > K, and let X be a spectrally positive Lévy process such that S t = se Xt is a martingale. Defining a measure Q by dq dp := ext, similar arguments as in the last part of the proof of Theorem 2 show that equation (19) follows from Theorem 4. Acknowledgement. We thank two anonymous referees for their careful reading. Not only did they point out several misprints, but their suggestions also substantially simplified the proof of Theorem 5. References [1] Andersen, L. and Lipton, A. Asymptotics for exponential Lévy processes and their volatility smile: survey and new results. Int. J. Theor. Appl. Finance 16 (2013), no. 1, 1350001, 98 pp. [2] Asmussen, S. and Rosiński, J. Approximations of small jumps of Lévy processes with a view towards simulation. J. Appl. Probab. 38 (2001), no. 2, 482 493. [3] Berestycki, H., Busca, J. and Florent, I. Asymptotics and calibration of local volatility models. Quant. Finance 2 (2002), no. 1, 61-69. [4] Ciomaga, A. On the strong maximum principle for second order nonlinear parabolic integro-differential equations. Adv. Differential Equations 17 (2012), no. 7-8, 635 671. [5] Cont, R. and Tankov, P. Financial modelling with Jump Processes. [6] Cont, R. and Voltchkova, E. Integro-differential equations for option prices in exponential Lévy models. Finance Stoch. 9 (2005), no. 3, 299 325.

14 ERIK EKSTRÖM AND BING LU [7] Figueroa-López, José E. and Forde, Martin. The small-maturity smile for exponential Lévy models. SIAM J. Financial Math. 3 (2012), 33 65. [8] Forde, M. and Jacquier, A. Small-time asymptotics for Implied volatility under the Heston model. Int. J. Theor. Appl. Finance 12 (2009), no. 6, 861 876. [9] Forde, M., Jacquier, A. and Lee, R The small-time smile and term structure of implied volatility under the Heston model. SIAM J. Financial Math. 3 (2012), 690 708. [10] Forde, M., Jacquier, A. Small-time asymptotics for an uncorrelated local-stochastic volatility model. Appl. Math. Finance 18 (2011), no. 6, 517 535. [11] Friedman, A. Partial differential equations of parabolic type. Prentice-Hall, Inc., Englewood Cliffs, N.J. 1964. [12] Gatheral, J., Hsu, E.P., Laurence, P., Ouyang, C. and Wang, T.H. Asymptotics of implied volatility in local volatility models. Math. Finance 22 (2012), no. 4, 591 620. [13] Lee, R. The moment formula for implied volatility at extreme strikes. Math. Finance 14 (2004), no. 3, 469 480. [14] Mijatovic, A. and Tankov, P. A new look at short-term implied volatility in asset price models with jumps. To appear in Mathematical Finance (2014). [15] Roper, M. Implied volatility: small time-to-expiry asymptotics in exponential Lévy models. Preprint (2010). [16] Signahl, M. On error rates in normal approximations and simulation schemes for Lévy processes. Stoch. Models 19 (2003), no. 3, 287 298. [17] Tankov, P. Pricing and Hedging in Exponential Lévy Models: Review of Recent Results. Lecture Notes in Math., 2003, Springer, Berlin, 2011. E-mail address: ekstrom@math.uu.se, bing@math.uu.se Department of Mathematics, Uppsala University, Box 480, SE-751 06 Uppsala, SWEDEN