Skew Modeling: Bruno Dupire Bloomberg LP Columbia University, New York September 12, 2005
Skew Modeling: Bruno Dupire Bloomberg LP Columbia University, New York September 12, 2005
Skew Modeling: Bruno Dupire Bloomberg LP Columbia University, New York September 12, 2005
Bruno Dupire
Bloomberg LP
bdupire@bloomberg.net
Columbia University, New York
September 12, 2005
I. Generalities
Market Skews
K
Not a general phenomenon
K K
Bruno Dupire 3
Skews
• Volatility Skew: slope of implied volatility as a
function of Strike
• Link with Skewness (asymmetry) of the Risk
Neutral density function ϕ ?
Bruno Dupire 4
Why Volatility Skews?
• Market prices governed by
– a) Anticipated dynamics (future behavior of volatility or jumps)
– b) Supply and Demand
σ impl Sup
ply
and
Market Skew
Dem
and
Th. Skew
K
Bruno Dupire 5
Modeling Uncertainty
Main ingredients for spot modeling
• Many small shocks: Brownian Motion
(continuous prices) S
Bruno Dupire 6
2 mechanisms to produce Skews (1)
– b) Downward jumps
Bruno Dupire 7
2 mechanisms to produce Skews (2)
– a) Negative link between prices and volatility
S1 S2
– b) Downward jumps
S1 S2
Bruno Dupire 8
Leverage and Jumps
Dissociating Jump & Leverage effects
t0 t1 t2
x = St1-St0 y = St2-St1
• Variance : (x + y) = x + 2xy+ y
2 2 2
Bruno Dupire 10
Dissociating Jump & Leverage effects
∑ (δS )
3
• Jump: ti
i
∑ (S )( )
− S t1 δS ti
2
• Leverage: ti
i
Bruno Dupire 11
Dissociating Jump & Leverage effects
Bruno Dupire 12
Dissociating Jump & Leverage effects
Bruno Dupire 13
Break Even Volatilities
Theoretical Skew from Prices
?
=>
Problem : How to compute option prices on an underlying without
options?
For instance : compute 3 month 5% OTM Call from price history only.
1) Discounted average of the historical Intrinsic Values.
Bad : depends on bull/bear, no call/put parity.
2) Generate paths by sampling 1 day return recentered histogram.
Problem : CLT => converges quickly to same volatility for all
strike/maturity; breaks autocorrelation and vol/spot dependency.
Bruno Dupire 15
Theoretical Skew from Prices (2)
Bruno Dupire 16
Theoretical Skew
from historical prices (3)
How to get a theoretical Skew just from spot price
history? S
Example: K
ST 1
Bruno Dupire 17
Theoretical Skew
from historical prices (4)
Bruno Dupire 18
Theoretical Skew
from historical prices (4)
Bruno Dupire 19
Theoretical Skew
from historical prices (4)
Bruno Dupire 20
Theoretical Skew
from historical prices (4)
Bruno Dupire 21
Barriers as FWD Skew trades
Beyond initial vol surface fitting
Bruno Dupire 23
Barrier Static Hedging
K
L
If not touched, IV’s are equal K L K
L
Bruno Dupire 24
Static Hedging: Model Dominance
• Back to DOCK , L
L
K
Bruno Dupire 25
Skew Adjusted Barrier Hedges
dS = (aS + b )dW
aK + b
DOC K , L ↔ CK − PaL2 +b (2 L − K )
aL + b aK + b
⎛ ⎞ aK + b
↔ C K − (L − K )⎜ 2 Dig L +
a
UOC K , L CL ⎟ − C aL2 +b (2 L − K )
⎝ aL + b ⎠ aL + b aK + b
Bruno Dupire 26
Local Volatility Model
One Single Model
• We know that a model with dS = σ(S,t)dW
would generate smiles.
– Can we find σ(S,t) which fits market smiles?
– Are there several solutions?
Bruno Dupire 28
The Risk-Neutral Solution
But if drift imposed (by risk-neutrality), uniqueness of the solution
Risk
Diffusions Neutral
Processes
sought diffusion
Compatible
(obtained by integrating twice
with Smile
Fokker-Planck equation)
Bruno Dupire 29
Forward Equations (1)
• BWD Equation:
price of one option C (K 0 ,T0 ) for different (S, t )
• FWD Equation:
price of all options C (K , T ) for current (S 0 ,t0 )
• Advantage of FWD equation:
– If local volatilities known, fast computation of implied
volatility surface,
– If current implied volatility surface known, extraction of
local volatilities,
– Understanding of forward volatilities and how to lock
them.
Bruno Dupire 30
Forward Equations (2)
• Several ways to obtain them:
– Fokker-Planck equation:
• Integrate twice Kolmogorov Forward Equation
– Tanaka formula:
• Expectation of local time
– Replication
• Replication portfolio gives a much more financial
insight
Bruno Dupire 31
Fokker-Planck
• If dx = b( x, t )dW
∂ϕ 1 ∂ 2 (b 2ϕ )
• Fokker-Planck Equation: =
∂t 2 ∂x 2
∂ 2C
• Where ϕ is the Risk Neutral density. As ϕ =
∂K 2
2 ⎛ ∂C ⎞
⎛ ∂ 2C ⎞ 2⎛ 2 ∂ C ⎞
2
∂ ⎜ ⎟ ∂⎜⎜ 2 ⎟⎟ ∂ ⎜⎜ b ⎟
2 ⎟
⎝ ∂t ⎠ = ⎝ ∂K ⎠ = 1 ⎝ ∂K ⎠
∂x 2 ∂t 2 ∂x 2
∂ ∂
T 2 T
C 1 C0 2
•Ito ( ST − K ) = C0 ( S 0 ,0) + ∫ dS + ∫ σ σ
+ 0
( − 2
0 ( S t , t )) dt
∂S 2 0 ∂S 2 t
0
• Taking expectation:
1
([ ] )
C ( S 0 ,0) = C0 ( S 0 ,0) + ∫∫ Γ0 ( S , t ) Ε σ t2 St = S − σ 02 ( S , t ) ϕ ( S , t )dSdt
[ ]
2
•Equality for all (K,T) Ù Ε σ t S t = S = σ 0 ( S , t )
2 2
Bruno Dupire 33
Summary of LVM Properties
Bruno Dupire 34
Stochastic Volatility Models
Heston Model
⎡ dS
⎢ S = µ dt + v dW
⎢
( )
⎢⎣dv = κ v 2 ∞ − v dt + η v dZ dW , dZ = ρ dt
FWD
x ≡ ln τ =T −t
K
C K ,T ( x, v,τ ) = e x P1 ( x, v,τ ) − P0 ( x, v,τ )
Bruno Dupire 36
Role of parameters
• Correlation gives the short term skew
• Mean reversion level determines the long term
value of volatility
• Mean reversion strength
– Determine the term structure of volatility
– Dampens the skew for longer maturities
• Volvol gives convexity to implied vol
• Functional dependency on S has a similar effect
to correlation
Bruno Dupire 37
Spot dependency
1) σ t = xt f (S , t ), ρ (W , Z ) = 0
2) σρ (W , Z ) ≠ 0 σ
S0 ST ST
S0
• With correlation ρ
Bruno Dupire 44
Smile Dynamics
Smile dynamics: Local Vol Model (1)
• Consider, for one maturity, the smiles associated
to 3 initial spot values
Smile S −
Smile S 0
Smile S +
S − S0 S + K
Bruno Dupire 51
Smile dynamics: Local Vol Model (2)
Smile S +
Smile S −
Smile S 0
S− S0 S+ K
Bruno Dupire 52
Smile dynamics: Stoch Vol Model (1)
Skew case (r<0) Local vols
σ Smile S −
Smile S 0
Smile S +
S − S0 S+ K
σ
Local vols
Smile S − Smile S +
Smile S 0
S− S0 S+ K
Bruno Dupire 54
Smile dynamics: Jump Model
Skew case
Local vols
Smile S −
Smile S 0
Smile S +
S− S0 S + K
Bruno Dupire 55
Smile dynamics: Jump Model
Pure smile case
Local vols
Smile S 0
Smile S − Smile S +
S− S0 S+ K
Bruno Dupire 56
Smile dynamics
Weighting scheme imposes
some dynamics of the smile for S1
a move of the spot: S0 K
For a given strike K,
S↑ ⇒ σ∃K ↓
(we average lower volatilities) t
Smile today (Spot St) 26
&
Smile tomorrow (Spot St+dt) 25.5
24.5
Smile tomorrow (Spot St+dt)
if σATM=constant 24
23.5
Smile tomorrow (Spot St+dt)
in the smile model
St+dt St
Bruno Dupire 57
Volatility Dynamics of different models
Bruno Dupire 58
ATM volatility behavior
Forward Skews
a) tells that the sensitivity and the hedge ratio of vanillas depend on the
calibration to the vanilla, not on local volatility/ stochastic volatility.
To change them, jumps are needed.
But b) does not say anything on the conditional forward skews.
Bruno Dupire 60
Sensitivity of ATM volatility / S
In average, σ ATM
2
follows σ loc
2
.
Optimal hedge of vanilla under calibrated stochastic volatility corresponds to
perfect hedge ratio under LVM.
Bruno Dupire 61
Market Model of Implied Volatility
• Implied volatilities are directly observable
• Can we model directly their dynamics? (r = 0)
⎧dS
⎪⎪ S = σ dW1
⎨
⎪dσˆ = α dt + u dW + u dW
⎪⎩ σˆ 1 1 2 2
Bruno Dupire 62
Drift Condition
• Apply Ito’s lemma to C(S,σˆ , t )
• Cancel the drift term
• Rewrite derivatives of C(S,σˆ , t )
gives the condition that the drift α of dσˆ must satisfy.
σˆ
For short T, we get the Short Skew Condition (SSC):
2 2
⎛ K ⎞ ⎛ K ⎞
C ( S ,σ ,̂ t)
t →T
σˆ = ⎜ σ + u1 ln( ) ⎟ + ⎜ u2 ln( ) ⎟
2
⎝ S ⎠ ⎝ S ⎠
K
close to the money: σˆ ~ σ + u1 ln( )
2
S
Î Skew determines u1
Bruno Dupire 63
Optimal hedge ratio ∆H
Bruno Dupire 64
Optimal hedge ratio ∆H II
dσˆ .dS
∆ = CS + Cσˆ
H
(dS ) 2
⎧ dS
⎪⎪ S = σdW1
With ⎨
⎪ dσˆ = αdt + u dW + u dW
⎪⎩ σˆ 1 1 2 2
K
S0
t0 T1 T2
Bruno Dupire 67
Det. Fut. smiles & no jumps
=> = FWD smile
If ∃(S , t , K , T ) / VK ,T (S , t ) ≠ σ (K , T ) ≡ lim σ imp (K , T , K + δK , T + δT )
2 2
δK → 0
δT → 0
stripped from Smile S.t
K
Then, there exists a 2 step arbitrage:
Define ∂ 2C S0
2
( (K , T ) − V (S , t )) ∂K (S , t , K , T )
PL t ≡ σ K ,T 2 S
(
At t0 : Sell PL t ⋅ Dig S − ε ,t − Dig S + ε ,t ) t0 t T
At t: if S ∈ [S − ε , S + ε ] buy
2
CS K, T , sell σ 2
(K , T )δ K ,T
t 2
K
gives a premium = PLt at t, no loss at T
Bruno Dupire 68
Consequence of det. future smiles
Bruno Dupire 69
Example of arbitrage with Sticky Strike
Each CK,T lives in its Black-Scholes (σ impl ( K , T ) )world
C1 ≡ C K 1 ,T1 C 2 ≡ C K 2 ,T2 assume σ 1 > σ 2
P&L of Delta hedge position over dt:
δ PL (C 1 ) = 1
2
((δ S ) − σ S δ t ) Γ
2
1
2
1
δ PL (C 2 ) = 1
2
((δ S ) − σ S δ t ) Γ
2
2
2
2
Γ1C 2
Γ2 C 1
Γ1Γ2 2 2
δ PL (Γ1C 2 − Γ2 C 1 ) =
2
(
S σ 1 − σ 22 δ t > 0 ) S t1 S t + δt
(no Γ , free Θ )
! If no jump
Bruno Dupire 70
Arbitrage with Sticky Delta
• In the absence of jumps, Sticky-K is arbitrageable and Sticky-∆ even more so.
• However, it seems that quiet trending market (no jumps!) are Sticky-∆.
In trending markets, buy Calls, sell Puts and ∆-hedge.
Example:
K1 St
PF ≡ C K 2 − PK1
K2
σ 1 ,σ 2
S PF
∆-hedged PF gains
VegaK > Vega K
2 1 from S induced
σ 1 ,σ 2 volatility moves.
S PF
VegaK < Vega K
2 1
Bruno Dupire 71
Conclusion
• Both leverage and asymmetric jumps may generate skew
but they generate different dynamics