Markets: A Rigorous Course

8.7 Risk-neutral valuation and the Black–Scholes formula

Core about 20 min

Motivation

The PDE can be solved by transforming to the heat equation. The faster route, and the one that generalises, is to write the solution as a discounted expectation under the risk-neutral measure and compute the integral. You will do the computation once in full; it is short.

Feynman–Kac and Girsanov#

Theorem 8.7.1Feynman–Kac representation; stated

The solution of (8.18) with V(T,S)=h(S)V(T,S)=h(S) and polynomial growth is

V(t,S)=e−r(T−t) E[h(ST) ∣ St=S],(8.19)V(t,S)=e^{-r(T-t)}\,\mathbb E\big[h(S_T)\,\big|\,S_t=S\big], \tag{8.19}

where the expectation is taken under dynamics with drift r−qr-q instead of μ\mu:

dSu=(r−q)Su du+σSu dW~u,ST=Stexp⁡ ⁣((r−q−12σ2)(T−t)+σ(W~T−W~t)).(8.20)\begin{gathered} dS_u=(r-q)S_u\,du+\sigma S_u\,d\widetilde W_u,\\ S_T=S_t\exp\!\Big(\big(r-q-\tfrac12\sigma^2\big)(T-t)+\sigma(\widetilde W_T-\widetilde W_t)\Big). \end{gathered} \tag{8.20}
Why

Let Mu=e−r(u−t)V(u,Su)M_u=e^{-r(u-t)}V(u,S_u) with SS following (8.20). Itô's lemma gives dMu=e−r(u−t)[Vu+(r−q)SVS+12σ2S2VSS−rV]du+e−r(u−t)σSVS dW~udM_u=e^{-r(u-t)}\big[V_u+(r-q)SV_S+\frac12\sigma^2S^2V_{SS}-rV\big]du+e^{-r(u-t)}\sigma SV_S\,d\widetilde W_u. The bracket is zero by (8.18), so MM is a local martingale; with growth conditions it is a true martingale, and V(t,S)=Mt=E[MT]=e−r(T−t)E[h(ST)]V(t,S)=M_t=\mathbb E[M_T]=e^{-r(T-t)}\mathbb E[h(S_T)]. See Shreve II, Thm 6.4.1.∎

Girsanov's theorem (stated, not proved) The tilde process lives under a different probability measure. Let WW be a P\mathbb P-Brownian motion on [0,T][0,T] and θ\theta a constant. Define Q\mathbb Q by the density

dQdP=exp⁡ ⁣(−θWT−12θ2T).(8.21)\frac{d\mathbb Q}{d\mathbb P}=\exp\!\Big(-\theta W_T-\tfrac12\theta^2T\Big). \tag{8.21}

Then W~t=Wt+θt\widetilde W_t=W_t+\theta t is a Q\mathbb Q-Brownian motion, and Q\mathbb Q and P\mathbb P agree on which events have probability zero. Choose θ=(μ−(r−q))/σ\theta=(\mu-(r-q))/\sigma, the market price of risk: excess return per unit of volatility. Then

dS=μS dt+σS dW=μS dt+σS(dW~−θ dt)=(r−q)S dt+σS dW~,dS=\mu S\,dt+\sigma S\,dW=\mu S\,dt+\sigma S(d\widetilde W-\theta\,dt)=(r-q)S\,dt+\sigma S\,d\widetilde W,

which is (8.20). Under Q\mathbb Q the discounted gains process e−(r−q)tSte^{-(r-q)t}S_t is a martingale; this is the continuous-time version of the risk-neutral pp of 8.3.

Intuition. For one variable, the density e−θx−θ2T/2e^{-\theta x-\theta^2T/2} is exactly the likelihood ratio φ−θT,T(x)/φ0,T(x)\varphi_{-\theta T,T}(x)/\varphi_{0,T}(x) of N(−θT,T)N(-\theta T,T) to N(0,T)N(0,T). Reweighting the paths of WW by it shifts the mean of WTW_T by −θT-\theta T and leaves the variance alone. Girsanov says this works path by path for the whole process: a change of measure can change the drift of a diffusion but never its volatility. That is why option prices depend on σ\sigma (which both measures share) and not on μ\mu (which the measure change removes). The proof is in Shreve II, Thm 5.2.3, or Karatzas and Shreve (1991), 3.5.

The Black–Scholes–Merton formula#

Theorem 8.7.2Black–Scholes–Merton formula

Let τ=T−t\tau=T-t and

d1=ln⁡(S/K)+(r−q+12σ2)τστ,d2=d1−στ=ln⁡(S/K)+(r−q−12σ2)τστ.(8.22)\begin{aligned} d_1&=\frac{\ln(S/K)+(r-q+\tfrac12\sigma^2)\tau}{\sigma\sqrt\tau},\\ d_2&=d_1-\sigma\sqrt\tau=\frac{\ln(S/K)+(r-q-\tfrac12\sigma^2)\tau}{\sigma\sqrt\tau}. \end{aligned} \tag{8.22}

Then the European call and put are worth

c=Se−qτN(d1)−Ke−rτN(d2),(8.23)c=Se^{-q\tau}N(d_1)-Ke^{-r\tau}N(d_2), \tag{8.23}
p=Ke−rτN(−d2)−Se−qτN(−d1).(8.24)p=Ke^{-r\tau}N(-d_2)-Se^{-q\tau}N(-d_1). \tag{8.24}
Proof

By (8.19)–(8.20), with Z=(W~T−W~t)/τ∼N(0,1)Z=(\widetilde W_T-\widetilde W_t)/\sqrt\tau\sim N(0,1) under Q\mathbb Q,

ST=Sexp⁡ ⁣((r−q−12σ2)τ+στZ).S_T=S\exp\!\big((r-q-\tfrac12\sigma^2)\tau+\sigma\sqrt\tau Z\big).

Step 1: the exercise region ST>KS_T>K iff στZ>ln⁡(K/S)−(r−q−12σ2)τ\sigma\sqrt\tau Z>\ln(K/S)-(r-q-\frac12\sigma^2)\tau iff Z>−d2Z>-d_2. So Q(ST>K)=Q(Z>−d2)=N(d2)\mathbb Q(S_T>K)=\mathbb Q(Z>-d_2)=N(d_2).

Step 2: split the integral

c=e−rτ∫−d2∞(Se(r−q−12σ2)τ+στz−K)φ(z) dz=I1−Ke−rτ∫−d2∞φ(z) dz=I1−Ke−rτN(d2).\begin{aligned} c&=e^{-r\tau}\int_{-d_2}^{\infty}\Big(Se^{(r-q-\frac12\sigma^2)\tau+\sigma\sqrt\tau z}-K\Big)\varphi(z)\,dz\\ &= I_1 - Ke^{-r\tau}\int_{-d_2}^\infty\varphi(z)\,dz = I_1-Ke^{-r\tau}N(d_2). \end{aligned}

Step 3: complete the square in I1I_1 Since φ(z)=e−z2/2/2π\varphi(z)=e^{-z^2/2}/\sqrt{2\pi},

e−12σ2τ+στzφ(z)=12πexp⁡ ⁣(−12(z2−2στz+σ2τ))=φ(z−στ).e^{-\frac12\sigma^2\tau+\sigma\sqrt\tau z}\varphi(z)=\frac{1}{\sqrt{2\pi}}\exp\!\Big(-\tfrac12\big(z^2-2\sigma\sqrt\tau z+\sigma^2\tau\big)\Big)=\varphi(z-\sigma\sqrt\tau).

Therefore

I1=Se−rτe(r−q)τ∫−d2∞φ(z−στ) dz=Se−qτ∫−d2−στ∞φ(y) dy=Se−qτN(d2+στ)=Se−qτN(d1),\begin{aligned} I_1&=Se^{-r\tau}e^{(r-q)\tau}\int_{-d_2}^{\infty}\varphi(z-\sigma\sqrt\tau)\,dz=Se^{-q\tau}\int_{-d_2-\sigma\sqrt\tau}^{\infty}\varphi(y)\,dy\\ &=Se^{-q\tau}N(d_2+\sigma\sqrt\tau)=Se^{-q\tau}N(d_1), \end{aligned}

using the substitution y=z−στy=z-\sigma\sqrt\tau and 1−N(−x)=N(x)1-N(-x)=N(x). This proves (8.23). For (8.24), apply parity (8.4): p=c−Se−qτ+Ke−rτ=Ke−rτ(1−N(d2))−Se−qτ(1−N(d1))p=c-Se^{-q\tau}+Ke^{-r\tau}=Ke^{-r\tau}(1-N(d_2))-Se^{-q\tau}(1-N(d_1)).∎

Interpretation of N(d1)N(d_1) and N(d2)N(d_2) N(d2)N(d_2) is the risk-neutral probability that the call finishes in the money. N(d1)N(d_1) is the probability of the same event under the measure that uses the stock (with dividends reinvested) as numéraire; Step 3 is exactly that change of measure, with the density ST/EQ[ST]S_T/\mathbb E^{\mathbb Q}[S_T] tilting ZZ by στ\sigma\sqrt\tau. In the CRR closed form (8.12) these are the limits of the pp-tail and the p′p'-tail. Neither is the real-world probability of exercise; that would use μ\mu in place of r−qr-q.

Limits and reference values#

Remark 8.7.3 (limits) As τ↓0\tau\downarrow0, d1,2→±∞d_{1,2}\to\pm\infty according to the sign of ln⁡(S/K)\ln(S/K), and c→(S−K)+c\to(S-K)^+. As σ↓0\sigma\downarrow0 with τ>0\tau>0, d1,2→±∞d_{1,2}\to\pm\infty according to the sign of ln⁡(F/K)\ln(F/K), and c→e−rτ(F−K)+c\to e^{-r\tau}(F-K)^+: with no uncertainty the stock surely ends at the forward. As σ→∞\sigma\to\infty, d1→∞d_1\to\infty and d2→−∞d_2\to-\infty, so c→Se−qτc\to Se^{-q\tau}, the upper bound (8.6). The formula spans exactly the no-arbitrage band; this is the basis for 8.10.

Worked example 8.7reference values

S=K=100S=K=100, τ=1\tau=1, r=5%r=5\%, q=0q=0, σ=20%\sigma=20\%. Then d1=(0+0.05+0.02)/0.2=0.35d_1=(0+0.05+0.02)/0.2=0.35, d2=0.15d_2=0.15, N(0.35)=0.636831N(0.35)=0.636831, N(0.15)=0.559618N(0.15)=0.559618, and

c=100×0.636831−100e−0.05×0.559618=63.6831−53.2325=10.4506,c=100\times0.636831-100e^{-0.05}\times0.559618=63.6831-53.2325=10.4506,

precisely 10.450583572210.4505835722, and p=5.5735260223p=5.5735260223. Check parity: 10.4506−5.5735=4.8771=100−100e−0.0510.4506-5.5735=4.8771=100-100e^{-0.05}. Hull's example (S=42S=42, K=40K=40, r=10%r=10\%, σ=20%\sigma=20\%, τ=0.5\tau=0.5) gives d1=0.7693d_1=0.7693, d2=0.6278d_2=0.6278, c=4.759422c=4.759422, p=0.808599p=0.808599. A Monte Carlo with 200,000 risk-neutral draws (m08_mc.mc_price, seed 42) gives 10.4634±0.033110.4634\pm0.0331 for the first case, within half a standard error.

Practitioner notes

  • The inputs are S,K,τ,r,q,σS,K,\tau,r,q,\sigma. Only σ\sigma is not observable. In India, qq for Nifty is also awkward: index dividends are lumpy (results season) and not continuous. That is the reason to price off the future (8.8).
  • Time: use calendar days/365 for discounting and either calendar or trading days for variance, but be consistent and know which one your vendor uses. For weekly options the difference is material: a Tuesday expiry seen from Friday is 4 calendar days but 2 trading days.

Create a free account

Reading is open to everyone. An account keeps your progress, answers, notes and code drafts, and gives you the study plan, the journals and the simulator. Running code needs an account approved for it.