Dupire Local Volatility
Dupire Local Volatility
K2
@2C
+ (rT
@K 2
qT ) C
@C
@K
rT C:
(1)
K(rT
qT )
@C
@K
1
qT C + K 2 E
2
2
T jST
@2C
:
@K 2
=K
vL =
y @w
w @y
1 @2w
2 @y 2
1
4
1
4
1
w
y2
w
@w
@y
(2)
=
(3)
1+
Ky @
@K
+2 T
@
@T
+K T
@
qT ) K @K
+ (rT
@
@K
1
4K
2
@
@K
@
+ K @K
2
(K; T ) =
@C
@T
(rT
qT ) C
1 2 @2C
2 K @K 2
@C
K @K
i:
(4)
If we set the risk-free rate rT and the dividend yield qT each equal to zero,
Equations (1) and (2) can each be solved to yield the same equation involving
local volatility, namely
2
(K; T ) =
@C
@T
1 2 @2C
K
2
@K 2
(5)
p
2 (K; T ): In this Note the derivation of
The local volatility is then vL =
these equations are all explained in detail.
=
t St dt + (St ; t)St dWt
= (rt qt ) St dt + (St ; t)St dWt :
RT
t
(6)
rs ds :
Fokker-Planck equation. Denote by f (St ; t) the probability density function of the underlying price St at time t. Then f satises the equation
@f
=
@t
@
1 @2
[ Sf (S; t)] +
@S
2 @S 2
S 2 f (S; t) :
(7)
= P (t; T )E (ST
K)+
(8)
where 1(ST >K) is the Heaviside function and where E [ ] = E [ jFt ]. In the
all the integrals in this Note, since the expectations are taken for the underlying
price at t = T it is understood that S = ST ; f (S; T ) = f (ST ; T ) and dS = dST :
We sometimes omit the subscript for notational convenience.
2
2.1
(9)
S=1
(10)
= P (t; T )f (K; T ):
We have assumed that lim f (S; T ) = 0:
S!1
2.2
@P
@T
(11)
(12)
Main Equation
2.3
2.3.1
(15)
1
f (S; T )dS:
@C
From the expression for @K
in Equation (9) we obtain
Z 1
1
@C
f (S; T )dS =
:
P
(t;
T
)
@K
K
Substitute back into Equation (15) and re-arrange terms to obtain the rst
identity
Z 1
K @C
C
:
(16)
ST f (S; T )dS =
P
(t;
T
)
P
(t;
T ) @K
K
2.3.2
Second Identity
@2C
@K 2
f (K; T ) =
2.4
1
@2C
:
P (t; T ) @K 2
(17)
First integral
=
=
[
[0
T (ST
0]
K; u0 = 1; v 0 =
S=1
@
@S
[Sf (S; T )] ; v =
Sf (S; T )dS
K @C
TC
+ T
:
P (t; T ) P (t; T ) @K
(18)
2.4.2
Second integral
K; u0 = 1; v 0 =
(ST
[0
2
=
where
S=1
@
@S
K)
0]
S 2 f (S; T )
S=K
S 2 f (S; T )
S=1
S=K
@2
@S 2
@
@S
S 2 f (S; T ) ; v =
S 2 f (S; T ) dS
K 2 f (K; T )
(K; T )2 .
@
S!1 @S
S 2 f (S; T )
= 0.
Substitute the second identity (17) for f (K; T ) to obtain the second integral I2
2
K 2 @2C
:
P (t; T ) @K 2
I2 =
2.5
(19)
1
I1 + I2 :
2
Substitute for I1 from Equation (18) and for I2 from Equation (19)
@C
+ rT C =
@T
Substitute for T = rT
Dupire equation (1)
@C
1
=
@T
2
Solve for
TC
TK
@C
1
+
@K
2
K2
@2C
@K 2
K2
@2C
+ (rT
@K 2
qT ) C
@C
@K
rT C:
@C
@T
@C
qK )K @K
+ qT C + (rT
1 2
2K
@2C
@K 2
Dupire [2] assumes zero interest rates and zero dividend yield. Hence rT =
qT = 0 so that the underlying process is dSt = (St ; t)St dWt : We obtain
(K; T )2 =
which is Equation (5).
5
@C
@T
1 2 @2C
K
2
@K 2
K)+ = 1(S>K)
@
@S 1(S>K)
= (S
@
@K (S
K)+ =
@
@K 1(S>K)
@C
@K
P (t; T )E 1(S>K)
1(S>K)
@2C
@K 2
K)
(S
K)
= P (t; T )E [ (S
K)]
In the table, ( ) denotes the Dirac delta function. Now dene the function
f (ST ; T ) as
f (ST ; T ) = P (t; T )(ST K)+ :
Recall the process for St is given by Equation (6). By Itos Lemma, f follows
the process
df =
@f
+
@T
T ST
@f
1
+
@ST
2
2
T ST
@2f
dT +
@ST2
@f
dWT :
@ST
T ST
(20)
K)+ ;
rT P (t; T )(ST
K) :
= P (t; T )
rT (ST
+P (t; T )
(21)
K)+ +
T ST 1(ST >K)
T ST 1(ST >K)
1
2
2 2
T ST
(ST
K) dT
dWT
K)+ +
T ST 1(ST >K)
=
rT (ST K)1(ST >K) + T ST 1(ST >K)
= rT K1(ST >K) qT ST 1(ST >K) :
When we take the expected value of Equation (21), the stochastic term drops
out since E [dWT ] = 0. Hence we can write the expected value of (21) as
dC
= E [df ]
= P (t; T )E rT K1(ST >K)
(22)
qT ST 1(ST >K) +
6
1
2
2 2
T ST
(ST
K) dT
so that
dC
= P (t; T )E rT K1(ST >K)
dT
qT ST 1(ST >K) +
1
2
2 2
T ST
(ST
K) .
(23)
K(rT
qT )
2 2
T ST
@C
@K
(ST
(24)
K)
1
qT C + P (t; T )E
2
2 2
T ST
(ST
K)
@C
where we have substituted @K
for P (t; T )E[1(ST >K) ]. The last term in the
last line of Equation (24) can be written
1
P (t; T )E
2
2 2
T ST
(ST
K)
=
=
=
K(rT
@2C
@K 2
qT )
1
P (t; T )E 2T ST2 jST = K E[ (ST
2
1
P (t; T )E 2T jST = K K 2 E[ (ST
2
1
@2C
E 2T jST = K K 2
2
@K 2
@C
@K
1
qT C + K 2 E
2
K)]
K)]
2
T jST
=K
@2C
:
@K 2
2
T jST
=K =
@C
@T
1 2 @2C
K
2
@K 2
which, again, is Equation (5). Hence when the dividend and interest rate are
both zero, the derivation of local volatility using Dupires approach and the
derivation using conditional expectation produce the same result.
To express local volatility in terms of implied volatility, we need the three deriv@C
@2C
atives @C
@T , @K ; and @K 2 that appear in Equation (1), but expressed in terms of
7
K
FT
RT
where FT = S0 exp 0 t dt is the forward price ( t = rt qt , risk free rate
minus dividend yield) and K is the strike price, and the "total" Black-Scholes
implied variance
w = (K; T )2 T
where (K; T ) is the implied volatility. The Black-Scholes call price can then
be written as
CBS (S0 ; K;
(K; T ) ; T )
where
d1 =
and d2 = d1
4.1
ln SK0 +
w=
yw
RT
0
1
2
= CBS (S0 ; FT ey ; w; T )
= FT fN (d1 ) ey N (d2 )g
(rt qt ) dt +
p
w
w
2
yw
1
2
(25)
1 1
+ w2
2
(26)
1 21
2w .
To express the local volatility Equation (1) in terms of y, we note that the
market call price is
C(S0 ; K; T ) = C(S0 ; FT ey ; T )
and we take derivatives. The rst derivative we need is, by the chain rule
@C
@C @K
@C
=
=
K:
@y
@K @y
@K
(27)
=
=
@
@C
@C @K
K+
@y @K
@K @y
@ 2 C 2 @C
K +
;
@K 2
@y
@A @K
@K @y ,
@C
@T
=
=
=
so that
@
@y
@C
@K
@C
@C @K
+
@T
@K @T
@C
@C
+
K T
@T
@K
@C
@C
+
@T
@y T
8
@ 2 C @K
@K 2 @y
(28)
@2C
@K 2 K.
The
(29)
RT
since K = S0 exp
that
t dt
+y
so that
@K
@T
= K
@2C
@2C 2
K =
2
@K
@y 2
T.
@C
:
@y
@C
@T
=
=
1
2
1
2
@2C
+ T C
@K 2
@C
@2C
+
2
@y
@y
K2
K
T
@C
@K
C
@C
@y
which simplies to
vL @ 2 C
@C
=
@T
2 @y 2
where vL =
[1].
4.2
@C
+
@y
TC
(30)
Before expression the local volatility Equation (1) in terms of implied volatility,
we rst derive three identities used by Gatheral [1] that help in this regard. We
use the fact that the derivatives of the standard normal cdf and pdf are, using
the chain rule, N 0 (x) = n(x)x0 and n0 (x) = xn(x)x0 . We also use the relation
n(d1 )
=
=
1
p e
2
1
p e
2
1
2
(d2 +
1
2
(d22 +2d2
= n(d2 )e d2
= n(d2 )ey :
w)
w+w)
1
2w
= FT [n(d1 )d1w
ey n(d2 )d2w ]
1
= FT n(d2 )ey d2w + w
2
h
i
1
1
=
FT ey n(d2 )w 2
2
1
2
ey n(d2 )d2w
where d1w is the rst derivative of d1 with respect to w and similarly for d2 .
The second derivative is
@ 2 CBS
@w2
1
3
1
1
FT ey
n(d2 )d2 d2w w 2
n(d2 )w 2
2
2
1
1
1
FT ey n(d2 )w 2
d2 d2w
w 1
2
2
1
3
@CBS
1 1
1
1
yw 2 + w 2
yw 2
w
@w
2
2
4
@CBS
1
1
y2
:
+
@w
8 2w 2w2
=
=
=
=
(31)
1
2
1
w
2
1 @
1
FT w 2
[ey n(d2 )]
2
@y
1
1
FT w 2 [ey n(d2 ) ey n(d2 )d2 d2y ]
2
@CBS
[1 d2 d2y ]
@w
@CBS 1
y
@w
2 w
=
=
=
=
(32)
= FT [n(d1 )d1y
= FT ey [n(d2 )d1y
=
FT ey N (d2 ):
ey N (d2 )
N (d2 )
ey n(d2 )d2y ]
n(d2 )d2y ]
FT ey N (d2 ) + FT ey n(d2 )w
@CBS
@CBS
+2
:
@y
@w
(33)
1
2
4.3
We note that when the market price C(S0 ; K; T ) is equal to the Black-Scholes
price with the implied volatility (K; T ) as the input to volatility
C(S0 ; K; T ) = CBS (S0 ; K; (K; T ); T ):
10
(34)
We can also reparameterize the Black-Scholes price in terms of the total implied
volatility w = (K; T )2 T and K = FT ey . Since w depends on K and K depends
on y, we have that w = w(y) and we can write
C(S0 ; K; T ) = CBS (S0 ; FT ey ; w(y); T ):
(35)
We need derivatives of the market call price C(S0 ; K; T ) in terms of the BlackScholes call price CBS (S0 ; FT ey ; w(y); T ). From Equation (35), the rst derivative we need is
@C
@y
@CBS
@CBS @w
+
@y
@w @y
= a(w; y) + b(w; y)c(y):
=
=
=
=
@2C
@y 2 ,
(36)
@a
@a @w
@c
@b
@b @w
+
+ b(w; y)
+
+
c(y)
(37)
@y @w @y
@y
@y @w @y
@ 2 CBS
@ 2 CBS @w @CBS @ 2 w
@ 2 CBS @w @w
@ 2 CBS
+
+
+
+
@y 2
@y@w @y
@w @y 2
@w@y
@w2 @y @y
@ 2 CBS @w @CBS @ 2 w @ 2 CBS
@ 2 CBS
+
2
+
+
@y 2
@y@w @y
@w @y 2
@w2
@w
@y
=
=
@CBS
@CBS @w
+
@T
@w @T
@CBS @w
:
TC +
@w @T
(38)
Gatheral explains that the second equality follows because the only explicit
dependence of CBS on T is through the forward price FT , even though CBS
depends implicitly on T through y and w. The reparameterized Dupire equation
(30) is reproduced here for convenience
@C
vL @ 2 C
=
@T
2 @y 2
We substitute for
tively and cancel
@CBS @w
@w @T
@C
+
@y
T C:
@C @ 2 C
@T ; @y 2 ;
and @C
@y from Equations (38), (37), and (36) respecC
from
both
sides to obtain
T
"
2
vL @ 2 CBS
@ 2 CBS @w @CBS @ 2 w @ 2 CBS @w
+
2
+
+
2
@y 2
@y@w @y
@w @y 2
@w2
@y
@CBS
@CBS @w
+
:
@y
@w @y
11
(39)
CBS @ CBS
@ CBS
Now substitute for @ @w
from the identities in Equations
2 ; @w@y ; and
@y 2
(31), (32), and (33) respectively, the idea being to end up with terms involving
@CBS
@w on the right hand side of Equation (39) that can be factored out.
"
2
@CBS @w
vL @CBS
1
y @w
1
1
y2
@w
=
2+2
+
+
@w @T
2 @w
2 w @y
8 2w 2w
@y
@2w
@y 2
@w
:
@y
BS
Remove the factor @C
@w from both sides and simplify to obtain
"
@w
1
1
y2
@w
y @w 1 @ 2 w 1
= vL 1
+
+
+
@T
w @y
2 @y 2
4
4 w
w
@y
Solve for vL to obtain the nal expression for the local volatility expressed in
2
terms of implied volatility w = (K; T ) T and the log-moneyness y = ln FKT
@w
@T
vL =
1
4.4
y @w
w @y
1 @2w
2 @y 2
1
4
1
4
1
w
y2
w
@w
@y
Alternate Derivation
2
@C @ C
In this derivation we express the derivatives @K
; @K 2 ; and @C
@T in the Dupire
equation (1) in terms of y and w = w(y), but we substitute these derivatives
directly in Equation (1) rather than in (30). This means that we take derivatives
with respect to K and T , rather than with y and T: Recall that from Equation
(35), the market call price is equal to the Black-Scholes call price with implied
volatility as input
ey N (d2 )g
p
w. The rst derivative
we need is
@C
@K
=
=
@CBS @y
@CBS @w
+
@y @K
@w @K
1 @CBS
@CBS @w
+
:
K @y
@w @K
12
(40)
1 @CBS
1 @
@CBS
+
:
K 2 @y
K @K
@y
@
@CBS @w
@CBS @ 2 w
+
+
@K
@w
@K
@w @K 2
Let A =
@C
@y
@
@K
=
=
=
@C
@y
@
@K
(41)
@A
@K
and
@A
@K
@A @y
@A @w
+
@y @K
@w @K
2
@ 2 CBS @w
@ CBS 1
+
:
2
@y K
@y@w @K
(42)
Similarly
@
@K
@CBS
@w
@ 2 CBS @w
@ 2 CBS 1
+
:
@y@w K
@w2 @K
(43)
1 @CBS
1 @ 2 CBS 1
+
2
K @y
K
@y 2 K
@ 2 CBS @w
@ 2 CBS 1
+
+
@y@w K
@w2 @K
2
1
@CBS
2
@ CBS
=
+
K2
@y 2
@y
K
=
@ 2 CBS
@w2
@w
@K
@ 2 CBS @w
@y@w @K
@CBS @ 2 w
@w
+
@K
@w @K 2
2
@ CBS @w
@y@w @K
+
(44)
@CBS @ 2 w
:
@w @K 2
=
=
@CBS
@CBS @y
@CBS @w
+
+
@T
@y @T
@w @T
@CBS
@CBS @w
+
;
T CBS +
@y T
@w @T
(45)
BS
again using the fact that @C
depends explicitly on T only through FT : Now
@T
2
@C @ C
@C
substitute for @K ; @K 2 ; and @T from Equations (40), (44), and (45) respectively
into Equation (4) for Dupire local variance, reproduced here for convenience.
@C
@T
T CBS
1 2 @2C
2 K @K 2
13
@C
K @K
BS
Applying the three useful identities in Section 4.2 allows the term @C
@w to be
factored out of the numerator and denominator. The last equation becomes
1 2
2K
2
K2
2
K
@w
@T +
@w
@K +
y
w
1
2
@w
T K @K
1
8
1
2w
@w 2
@K
y2
2w2
@2w
@K 2
i:
(46)
@
@w
@
@ w
@
@
2 T @T
+ 2 , @K
= 2 T @K
, and @K
+ @K
: Substitute into
2 = 2T
2
@K
Equation (46). The numerator in Equation (46) becomes
2
+2 T
@
+
@T
TK
@
@K
(47)
1
2
+K 2 T
@
@K
@2
+K 2 T
@K 2
Ky @
=
1
@K
=
1
y2
+
2w 2w2
1
8
T2
@
@K
1
2
@
@K
@
+ 2K 2
@K
#
@2
:
@K 2
@
@K
Replacing w with
1 + 2K T
"
y
w
y
2T
2
@
+ 2K 2
@K
#
2Ky @
@K
T2
1
8
1
2T
y2
2 4T 2
@
@K
"
K2
T2
Ky @
2
+
@K
+ 1
@
@K
Ky @
@K
K 2 y2
2
@
@K
(48)
Substituting the numerator in (47) and the denominator in (48) back to Equation (46), we obtain
@
@
+ 2 T @T
+ T K @K
h
2
@
1
@
+ K T @K
4 K T @K
2
1+
Ky @
@K
@2
@K 2
1+K T
14
@
+ K @K
2
See also the dissertation by van der Kamp [4] for additional details of this
alternate derivation.
References
[1] Gatheral, J. (2006). The Volatility Surface: A Practitioners Guide. New
York, NY: John Wiley & Sons.
[2] Dupire, B. (1994). "Pricing With a Smile." Risk 7, pp. 18-20.
[3] Derman, E., Kani, I., and M. Kamal (1996). "Trading and Hedging Local
Volatility." Goldman Sachs Quantitative Strategies Research Notes.
[4] van der Kamp, Roel (2009). "Local Volatility Modelling." M.Sc. dissertation, University of Twente, The Netherlands.
15