Joint work with José Manuel Corcuera, Peter Leoni and Wim Schoutens July 15, 2009 - Eurandom
1 2 The Black-Scholes model The Lévy models 3 4 5 6 7 Delta Hedging at versus at Implied Black-Scholes Volatility The Historical Optimal P&L
introduction of the concept of implied Lévy volatility (extension of Black-Scholes implied volatility) Lévy implied time volatility Lévy implied space volatility Study of the shape of implied Lévy volatilities Model performance delta-hedging strategies (periodical rebalancing) qualitatively (Greeks) historical time-series of the S&P500 historical option prices of the Dow Jones
The Black-Scholes model The Lévy models The Black-Scholes model Diffusion part of the log-return process: modelled by geometrical Brownian motion (W t ) Definition S t = S 0 exp((r q σ 2 /2)t + σw t ), t 0. The Black-Scholes implied volatility is the volatility σ = σ(k, T ) such that the model and market option prices coincide. σ = σ(k, T ) volatility surface σ needs to be adjusted separately for each individual contract Historical stock returns: skewed and fatter tails than those of the normal distribution Development of a similar concept but now under a Lévy framework based on more empirically founded distributions
The Black-Scholes model The Lévy models The Lévy space model Lévy space stock price model: S t = S 0 exp((r q + ω)t + σx t ), t 0, where E[X 1 ] = 0, Var[X 1 ] = 1 and ω = log(φ( σi)) where φ characteristic function of X 1 : φ(u) = E [exp(iux 1 )] Note: E[X t] = 0 and Var[X t] = t Var[σX t] = σ 2 t Definition The volatility parameter σ = σ(k, T ) needed to match the model price with a given market price is called the implied Lévy space volatility.
The Black-Scholes model The Lévy models The Lévy time model Lévy time stock price model: S t = S 0 exp((r q + ωσ 2 )t + X σ2 t), t 0, where E[X 1 ] = 0, Var[X 1 ] = 1 and ω = log(φ( i)) Note: Var[X σ 2 t ] = σ 2 t Definition The volatility parameter σ = σ(k, T ) needed to match the model price with the market price is called the implied Lévy time volatility.
Pricing Vanillas under Lévy Models Carr-Madan formula in combination with FFT gives a very fast evaluation of vanillas: where C(K, T ) = exp( α log(k)) π + 0 exp( iv log(k))ϱ(v)dv, ϱ(v) = exp( rt )E[exp(i(v (α + 1)i) log(s T ))] α 2 + α v 2. + i(2α + 1)v Only dependence of the Carr-Madan formula on the model: risk neutral (i.e. under Q) characteristic function of the log-price process at maturity T : φ(u; T ) = E Q [exp(iu log(s T ))]. This characteristic function is available in closed-form for many popular Lévy processes.
Greeks under Lévy models Delta and many others Greeks can be calculated in a similar fashion: = C(K, T ) S 0 where = ϱ (v) = exp( α log(k)) π + 0 exp( rt )φ(v (α + 1)i; T S 0 (α + iv) exp( iv log(k))ϱ dv.
Option prices and Greeks computation COS method rests on Fourier-cosine series expansions and can be applied for any model if the characteristic function ψ(u; T ) is available where ψ is the characteristic function of the log-moneyness at maturity ( ( ))] ST ψ(u; T ) = E Q [exp iu log K (see Fang, F. and Oosterlee, C.W. (2008) A novel pricing method for European Options based on Fourier-cosine Series Expansions. SIAM Journal on Scientific Computing 31-2, 826-848. )
NIG Characteristic function of the normal inverse Gaussian distribution NIG(α, β, δ, µ) with parameters α > 0, β ] α, α[, δ > 0 and µ R: ( φ NIG(u; α, β, δ, µ) = exp iuµ δ ( α 2 (β + iu) 2 α 2 β 2)), u R. If the parameter β is equal to zero the distribution is symmetric around µ whereas negative and positive values of β result in negative and positive skewness NIG(α, β, δ, µ) NIG(α, 0, δ, µ) mean µ + δβ α 2 β 2 µ variance α 2 δ ( α 2 β 2) 3/2 δ α skewness 3βα 1 δ ( 1/2 α 2 β 2) 1/4 0 ) kurtosis 3 (1 + α2 +4β 2 δα 2 α 2 β 2 3 ( 1 + 1 αδ )
Meixner Characteristic function of the Meixner distribution Meixner(α, β, δ, µ) with parameters α > 0, β ] π, π[, δ > 0 and µ R: φ Meixner(u; α, β, δ, µ) = exp ( iuµ ) ( cos β 2 cosh ) ( αu iβ 2 2δ ), u R. A parameter β equal to zero indicates a symmetric distribution around µ whereas negative and positive values of β lead to negative and positive skewness mean variance skewness kurtosis Meixner(α, β, δ, µ) ) Meixner(α, 0, δ, µ) µ + αδ tan µ sin α 2 δ 2 cos 2 ( β 2 ) ( β 2 ( β 2 3 + 2 cos(β) δ α 2 δ 2 ) 2 δ 0 3 + 1 δ
Lévy waves we compute the implied Lévy space and time volatility for the NIG Lévy model for various Black-Scholes implied volatility shape and vice-versa (T = 1, r = q = 0, S 0 = 100). The symmetric cases: σ NIG Space σ NIG (κ = 0, α) for inverse square root σ 0.24 0.23 0.22 0.21 0.2 0.19 0.18 0.17 0.16 0.15 70 80 90 100 110 120 130 K σ α = 1 α = 1.25 α = 1.5 α = 1.75 α = 2 α = 2.25 α = 2.5 α = 2.75 α = 3 α = 3.25 α = 3.5 α = 3.75 α = 4 α = 4.25 α = 4.5 α = 4.75 α = 5 σ NIG Time σ NIG (κ = 0, α) for inverse square root σ 0.24 σ 0.22 α = 2.5 α = 2.75 0.2 α = 3 α = 3.25 0.18 α = 3.5 α = 3.75 α = 4 0.16 α = 4.25 α = 4.5 0.14 α = 4.75 α = 5 0.12 70 80 90 100 110 120 130 K Figure: Implied volatility for the symmetric space (left) and time (right) NIG models.
Lévy Waves Some asymmetric cases: σ NIG Space σ NIG (κ = 0.625, α) for inverse square root σ σ 0.25 0.2 0.15 70 80 90 100 110 120 130 K α = 1 α = 1.25 α = 1.5 α = 1.75 α = 2 α = 2.25 α = 2.5 α = 2.75 α = 3 α = 3.25 α = 3.5 α = 3.75 α = 4 α = 4.25 α = 4.5 α = 4.75 α = 5 Time σ NIG (κ = 0.375, α) for inverse square root σ 0.3 σ 0.28 α = 1.75 α = 2 0.26 α = 2.25 α = 2.5 0.24 α = 2.75 α = 3 α = 3.25 0.22 α = 3.5 α = 3.75 0.2 α = 4 α = 4.25 0.18 α = 4.5 α = 4.75 0.16 α = 5 70 80 90 100 110 120 130 K σ NIG Figure: Implied volatility for inverse square root volatility for asymmetric asymmetric NIG space models (κ = 0.625) (left) and asymmetric NIG time models (κ = 0.375) (right).
What is flat here is not flat there A flat Black-Scholes implied volatility curve corresponds to a reversed smiling Lévy implied volatility curve under each symmetric NIG model. σ NIG 0.215 0.21 0.205 0.2 0.195 0.19 0.185 0.18 Space σ NIG (κ = 0, α) for a flat σ 0.175 70 80 90 100 110 120 130 K σ α = 1 α = 1.25 α = 1.5 α = 1.75 α = 2 α = 2.25 α = 2.5 α = 2.75 α = 3 α = 3.25 α = 3.5 α = 3.75 α = 4 α = 4.25 α = 4.5 α = 4.75 σ NIG 0.24 0.23 0.22 0.21 0.2 0.19 0.18 0.17 0.16 Time σ NIG (κ = 0, α) for a flat σ 0.15 α = 5 70 80 90 100 110 120 130 K σ α = 2.5 α = 2.75 α = 3 α = 3.25 α = 3.5 α = 3.75 α = 4 α = 4.25 α = 4.5 α = 4.75 α = 5 Figure: Implied space (left) and time (right) volatility for a flat volatility for some symmetric NIG distributions.
What is flat here is not flat there A flat NIG implied volatility curve corresponds to a smiling implied Black-Scholes volatility curve. σ 0.225 0.22 0.215 0.21 0.205 0.2 0.195 0.19 0.185 σ for a flat space σ NIG (κ = 0, α) 0.18 70 80 90 100 110 120 130 K σ NIG α = 1 α = 1.25 α = 1.5 α = 1.75 α = 2 α = 2.25 α = 2.5 α = 2.75 α = 3 α = 3.25 α = 3.5 α = 3.75 α = 4 α = 4.25 α = 4.5 α = 4.75 σ 0.24 0.23 0.22 0.21 0.2 0.19 0.18 0.17 σ for a flat time σ NIG (κ = 0, α) 0.16 α = 5 70 80 90 100 110 120 130 K σ NIG α = 2.5 α = 2.75 α = 3 α = 3.25 α = 3.5 α = 3.75 α = 4 α = 4.25 α = 4.5 α = 4.75 α = 5 Figure: Implied volatility for a flat space (left) and time (right) volatility for some symmetric NIG distributions.
Vanna Vanna = 2 C S 0 σ = σ 0.9 0.85 0.8 0.75 Delta in function of the volatility (ATM Call) NIG space (α = 1.25) NIG space (α = 1.75) NIG time (α = 1.25) NIG time (α = 1.75) 0.5 0 0.5 Vanna in function of the volatility (ATM Call) 0.7 1 Delta 0.65 0.6 Vanna 1.5 2 0.55 0.5 0.45 2.5 3 NIG space (α = 1.25) NIG space (α = 1.75) NIG time (α = 1.25) NIG time (α = 1.75) 0.4 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 σ 3.5 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 σ Figure: Delta (left) and Vanna (right) as a function of the volatility for a maturity of 1 month (ATM option).
Charm Charm = T 0.8 0.75 Delta in function of the maturity (ATM Call) 0.4 0.35 Charm in function of the maturity (ATM Call) NIG space (α = 1.25) NIG space (α = 1.75) 0.7 0.3 0.25 NIG time (α = 1.25) NIG time (α = 1.75) Delta 0.65 Charm 0.2 0.6 NIG space (α = 1.25) 0.15 0.1 0.55 NIG space (α = 1.75) NIG time (α = 1.25) 0.05 NIG time (α = 1.75) 0.5 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 T 0 0.5 1 1.5 2 2.5 3 T Figure: Delta (left) and Charm (right) as a function of the option maturity for a volatility of 0.2 (ATM option)
Hedging error Delta Hedging at versus at Implied Black-Scholes Volat The Historical Optimal P&L Each day, we delta-hedge an ATM one-month Call option for one single day. The next day, we start with a new ATM option with again one whole month as a lifetime. Hedging indicator: HE(t 0 + t) = C t0+ t(k, T t) t0 S t0+ t+( t0 S t0 C t0 (K, T )) e r t. By considering an implied Lévy model different from the Black-Scholes model, different values of the free parameters will lead to different distributions of the hedging error (HE). Optimal free parameter set p = {p 1,... p n}: abs(µ HE ( p )) + σ HE ( p ) abs(µ HE ( p)) + σ HE ( p). Data set: VIX and S&P500 from the 2nd January 1990 to the 9th October 2008.
Evolution of the implied volatility Evolution of the implied volatility through time: Delta Hedging at versus at Implied Black-Scholes Volat The Historical Optimal P&L 0.9 0.8 Evolution of the volatility NIG space(α = 1) 0.9 0.8 Evolution of the volatility Meixner space (α = 3.45) 0.7 NIG time (α = 1.25) 0.7 Meixner time (α = 1.8) 0.6 0.6 0.5 0.5 σ 0.4 σ 0.4 0.3 0.3 0.2 0.2 0.1 0.1 0 0 500 1000 1500 2000 2500 3000 3500 4000 4500 Trading day 0 0 500 1000 1500 2000 2500 3000 3500 4000 4500 Trading day Figure: Evolution of the implied volatility for the NIG (upper) and Meixner (lower) volatility models The implied Lévy space and time volatility moves in line with the implied Black-Scholes volatility.
Evolution of the Delta Delta Hedging at versus at Implied Black-Scholes Volat The Historical Optimal P&L Evolution of the Delta through time: Evolution of the Delta Evolution of the Delta 0.6 0.55 NIG space (α = 1) NIG time (α = 1.25) 0.55 0.5 Delta 0.5 Delta 0.45 0.4 0.45 0.4 0 500 1000 1500 2000 2500 3000 3500 4000 4500 Trading day 0.35 0.3 Meixner space (α = 3.45) Meixner time (α = 1.8) 0 500 1000 1500 2000 2500 3000 3500 4000 4500 Trading day Figure: Evolution of the Delta for the NIG (upper) and Meixner (lower) volatility models The higher slope of the Vanna curve under the Lévy models explains the higher volatility of the Lévy Deltas through time in comparison with the Black-Scholes Deltas.
The Historical Optimal Implied NIG Volatility Delta Hedging at versus at Implied Black-Scholes Volat The Historical Optimal P&L HE( t) variance Variance of HE( t) NIG model (κ = 0) 3 2.8 2.6 2.4 2.2 2 1.8 1.6 1.4 NIG space NIG time 1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2 α HE( t) functional = abs(mean) + std 1.8 1.7 1.6 1.5 1.4 Functional of HE( t) NIG model (κ = 0) 1.3 1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2 α NIG space NIG time Figure: Variance (left) and functional (right) of the hedging error distribution for the NIG volatility models.
Delta Hedging at versus at Implied Black-Scholes Volat The Historical Optimal P&L The Historical Optimal Implied Meixner Volatility HE( t) variance 7 6 5 4 3 2 Variance of HE( t) Meixner model (κ = 0) Meixner space Meixner time 1 0 0.5 1 1.5 2 2.5 3 3.5 4 α HE( t) functional = abs(mean) + std Functional of HE( t) Meixner model (κ = 0) 2.8 2.6 Meixner space Meixner time 2.4 2.2 2 1.8 1.6 1.4 0 0.5 1 1.5 2 2.5 3 3.5 4 α Figure: Variance (left) and functional (right) of the hedging error distribution for the Meixner volatility models. The historical optimal model (NIG time model with α = 1.25) leads to a reduction of the HE variance and functional from around 3 to less than 1.5 and from around 1.85 to 1.3, respectively.
Dow Jones Profit and loss Delta Hedging at versus at Implied Black-Scholes Volat The Historical Optimal P&L Delta-hedging strategy: At time t 0: sell option + buy t 0 stocks At time 0 < t i < T : buy ti ti 1 stocks (rebalancing) At time T : close the option and stock positions balance at time t i amount spent until time t i to build the hedging portfolio: and Balance(t 0 ) = C t0 (K, T ) t0 S t0 Balance(t i ) = Balance(t i 1 )e r t + CF(t i ), where CF rebalance cash flow: 0 < t i T CF(t i ) = ( ti ti 1 ) Sti
Delta Hedging at versus at Implied Black-Scholes Volat The Historical Optimal P&L Dow Jones Profit and loss (cont.) Mark to Market at time t i amount spent until time t i to build the hedging portfolio after the closing of the option and stock positions: MtM(t i ) = Balance(t i ) + C ti (K, T t i ) t i S ti P&L = MtM(T ) Data : 491 liquid Put and Call option prices (different K and T ) on the Dow Jones 2 cases: sell each option once sell each option for an amount of 1$? Model which minimises the P&L variance
Dow Jones P&L results Delta Hedging at versus at Implied Black-Scholes Volat The Historical Optimal P&L 5 4.5 Variance of the P&L NIG model (κ = 0) NIG space NIG time 0.045 0.04 Weighted variance of the P&L NIG model (κ = 0) NIG space NIG time P&L variance 4 3.5 3 2.5 P&L weighted variance 0.035 0.03 0.025 0.02 2 1 1.5 2 2.5 3 3.5 4 α 0.015 1 1.5 2 2.5 3 3.5 4 α Figure: Variance (left) and weighted variance (right) of the P&L for the NIG volatility models.
Dow Jones P&L results Delta Hedging at versus at Implied Black-Scholes Volat The Historical Optimal P&L 3.5 Variance of the P&L Meixner model (κ = 0) Meixner space Meixner time Weighted variance of the P&L Meixner model (κ = 0) 0.03 Meixner space 0.028 Meixner time P&L variance 3 2.5 P&L weighted variance 0.026 0.024 0.022 0.02 2 0.5 1 1.5 2 2.5 α 0.018 0.5 1 1.5 2 2.5 α Figure: Variance (left) and weighted variance (right) of the P&L for the Meixner volatility models.
Dow Jones P&L results Delta Hedging at versus at Implied Black-Scholes Volat The Historical Optimal P&L variance of the P&L 5 4.5 4 3.5 3 2.5 2 1.5 1 0.5 Evolution of the variance of the P&L 0 03/01/00 02/01/01 02/01/02 02/01/03 02/01/04 03/01/05 03/01/06 quote date standard VG space variance of the weighted P&L 0.05 0.045 0.04 0.035 0.03 0.025 0.02 0.015 0.01 0.005 Evolution of the variance of the weighted P&L 0 3/01/2000 2/01/2001 2/01/2002 2/01/2003 2/01/2004 3/01/2005 3/01/2006 quote date standard VG space Figure: Variance (left) and weighted variance (right) of the P&L for the standard VG model volatility model.
Delta Hedging at versus at Implied Black-Scholes Volat The Historical Optimal P&L Conclusion implied Lévy space and time models obtained by replacing the normal distribution of the Black-Scholes model by a more suitable Lévy distribution Switching from the world to the Lévy world additional dof which can be used to minimize the curvature of the volatility surface any smiling or smirking volatility curve can be transformed into a flatter Lévy volatility curve under a well chosen parameter set implied Lévy models could lead to flatter volatility curves for more practical datasets minimise the absolute mean and the square root of the variance of the hedging error using the historical optimal parameters leads to a significant reduction of the hedging error minimise the variance of the P&L