Derivatives PricingBlack-ScholesGreeksPDE Methods

Black-Scholes: Derivation, Greeks, Limitations

Module 1 of 422 min readLevel: Medium

Setup

Market Assumptions

Let (Ω,F,P)(\Omega, \mathcal{F}, \mathbb{P}) be a filtered probability space supporting a standard Brownian motion (Wt)t0(W_t)_{t \geq 0}. The stock price (St)t0(S_t)_{t \geq 0} follows geometric Brownian motion (GBM) under the physical measure P\mathbb{P}:

dSt=μStdt+σStdWt,S0>0.dS_t = \mu S_t \, dt + \sigma S_t \, dW_t, \qquad S_0 > 0.

The following assumptions are in force throughout. Each is a modelling choice, not a physical law — and each is violated in practice:

  1. Constant volatility. σ>0\sigma > 0 is a deterministic constant.
  2. Constant risk-free rate. The risk-free rate r0r \geq 0 is constant and continuously compounded.
  3. No dividends. The stock pays no cash dividends.
  4. Continuous trading. Portfolio rebalancing can occur at every instant at zero cost.
  5. No transaction costs, no bid-ask spread. All trades execute at the mid-price.
  6. No short-selling constraints. The stock can be sold short without restriction.
  7. Markets are frictionless and complete. Every contingent claim is replicable.

Under these assumptions, the Black-Scholes framework is internally consistent. The question is not whether it is true — it is not — but whether it is useful as a baseline and how its failure modes manifest in practice.


Derivation via Delta Hedging

The original Black-Scholes (1973) derivation proceeds by constructing a locally riskless portfolio.

Let C(t,S)=C(t,St)C(t, S) = C(t, S_t) be the price at time tt of a European call option with strike KK and maturity TT. Assume CC1,2([0,T)×(0,))C \in C^{1,2}([0,T) \times (0,\infty)).

Step 1: Apply Itô's lemma. Since StS_t satisfies the GBM SDE:

dC=(Ct+μSCS+12σ2S22CS2)dt+σSCSdWt.dC = \left(\frac{\partial C}{\partial t} + \mu S \frac{\partial C}{\partial S} + \frac{1}{2}\sigma^2 S^2 \frac{\partial^2 C}{\partial S^2}\right) dt + \sigma S \frac{\partial C}{\partial S} \, dW_t.

Step 2: Form the delta-hedged portfolio. Define the portfolio

Πt=C(t,St)ΔtSt,Δt=CS(t,St).\Pi_t = C(t, S_t) - \Delta_t \cdot S_t, \qquad \Delta_t = \frac{\partial C}{\partial S}(t, S_t).

Its differential is:

dΠt=dCΔtdSt=(Ct+12σ2S22CS2)dt.d\Pi_t = dC - \Delta_t \, dS_t = \left(\frac{\partial C}{\partial t} + \frac{1}{2}\sigma^2 S^2 \frac{\partial^2 C}{\partial S^2}\right) dt.

The dWtdW_t terms cancel exactly because Δt\Delta_t is chosen as the option's partial derivative with respect to SS. The portfolio is instantaneously riskless.

Step 3: Apply no-arbitrage. A riskless portfolio must earn the risk-free rate:

dΠt=rΠtdt=r(CΔtS)dt.d\Pi_t = r \Pi_t \, dt = r(C - \Delta_t S) \, dt.

Step 4: Equate and collect. Substituting:

Ct+12σ2S22CS2=rCrSCS.\frac{\partial C}{\partial t} + \frac{1}{2}\sigma^2 S^2 \frac{\partial^2 C}{\partial S^2} = r C - r S \frac{\partial C}{\partial S}.

Rearranging gives the Black-Scholes PDE:

Ct+12σ2S22CS2+rSCSrC=0,(t,S)[0,T)×(0,),\frac{\partial C}{\partial t} + \frac{1}{2}\sigma^2 S^2 \frac{\partial^2 C}{\partial S^2} + r S \frac{\partial C}{\partial S} - r C = 0, \qquad (t, S) \in [0, T) \times (0, \infty),

subject to the terminal condition C(T,S)=(SK)+C(T, S) = (S - K)^+ and boundary conditions C(t,0)=0C(t, 0) = 0, C(t,S)SC(t, S) \sim S as SS \to \infty.

Remark: The Drift μ\mu Does Not Appear

The physical drift μ\mu cancels in the portfolio construction. This is not a coincidence: it reflects the risk-neutral pricing principle derivable via Girsanov's theorem. The no-arbitrage price of any replicable claim depends only on σ\sigma and rr, not on the investor's expected return.


Derivation via Feynman-Kac

A cleaner derivation uses the Feynman-Kac representation. Under the risk-neutral measure Q\mathbb{Q}, defined by Girsanov's theorem with market price of risk λ=(μr)/σ\lambda = (\mu - r)/\sigma, the stock follows:

dSt=rStdt+σStdW~t,dS_t = r S_t \, dt + \sigma S_t \, d\widetilde{W}_t,

where W~t=Wt+λt\widetilde{W}_t = W_t + \lambda t is a Q\mathbb{Q}-Brownian motion. The no-arbitrage price is:

C(t,S)=er(Tt)EQ ⁣[(STK)+Ft].C(t, S) = e^{-r(T-t)} \mathbb{E}^{\mathbb{Q}}\!\left[(S_T - K)^+ \,\Big|\, \mathcal{F}_t\right].

Under Q\mathbb{Q}, the log-price is:

lnST=lnSt+(rσ22)τ+στZ,ZN(0,1),τ=Tt.\ln S_T = \ln S_t + \left(r - \frac{\sigma^2}{2}\right)\tau + \sigma\sqrt{\tau}\, Z, \qquad Z \sim \mathcal{N}(0,1), \quad \tau = T - t.

Evaluating the expectation by splitting the integration region {ST>K}\{S_T > K\} yields the Black-Scholes formula:

C(t,S)=SN(d1)KerτN(d2),\boxed{C(t, S) = S \, N(d_1) - K e^{-r\tau} N(d_2),}

d1=ln(S/K)+(r+σ2/2)τστ,d2=d1στ.d_1 = \frac{\ln(S/K) + (r + \sigma^2/2)\,\tau}{\sigma\sqrt{\tau}}, \qquad d_2 = d_1 - \sigma\sqrt{\tau}.

Here N()N(\cdot) denotes the standard normal CDF and n()n(\cdot) its density. By put-call parity,

P(t,S)=KerτN(d2)SN(d1).P(t, S) = K e^{-r\tau} N(-d_2) - S\, N(-d_1).

Notation

Throughout: τ=Tt\tau = T - t (time to expiry), n(x)=12πex2/2n(x) = \frac{1}{\sqrt{2\pi}}e^{-x^2/2} (standard normal pdf), N(x)=xn(u)duN(x) = \int_{-\infty}^x n(u)\,du (standard normal CDF).


The Five Greeks

The Greeks measure sensitivity of the option price to each model input. All five are derived analytically from the closed-form formula. A useful identity, easily verified from the definition of d1d_1:

Sn(d1)=Kerτn(d2).S \, n(d_1) = K e^{-r\tau} n(d_2).

This identity allows equivalent representations and simplifies many derivations.

Delta: Δ=C/S\Delta = \partial C / \partial S

Δcall=N(d1),Δput=N(d1)1=N(d1).\Delta_{\mathrm{call}} = N(d_1), \qquad \Delta_{\mathrm{put}} = N(d_1) - 1 = -N(-d_1).

Interpretation. Δ\Delta is both the hedge ratio (shares needed to replicate) and, under Q\mathbb{Q}, the probability that the call expires in the money when the underlying is a forward: Δcall=N(d1)Q(ST>K)\Delta_{\mathrm{call}} = N(d_1) \approx \mathbb{Q}(S_T > K) only if r=0r = 0; otherwise Q(ST>K)=N(d2)\mathbb{Q}(S_T > K) = N(d_2). The distinction matters at long maturities.

Gamma: Γ=2C/S2\Gamma = \partial^2 C / \partial S^2

Γ=n(d1)Sστ.\Gamma = \frac{n(d_1)}{S \sigma \sqrt{\tau}}.

Γ\Gamma is identical for calls and puts (put-call parity). It measures convexity: the rate at which the hedge ratio must be adjusted. High Γ\Gamma → frequent rebalancing → high transaction costs.

Theta: Θ=C/t\Theta = \partial C / \partial t

Θcall=Sn(d1)σ2τrKerτN(d2).\Theta_{\mathrm{call}} = -\frac{S\, n(d_1)\, \sigma}{2\sqrt{\tau}} - r K e^{-r\tau} N(d_2).

Θ\Theta is almost always negative for long calls (the option loses time value). For deep in-the-money calls, Θ\Theta can be positive (early exercise premium in American options). The Black-Scholes PDE is often written as:

Θ+12σ2S2Γ+rSΔrC=0,\Theta + \frac{1}{2}\sigma^2 S^2 \Gamma + r S \Delta - r C = 0,

which expresses the P&L decomposition: time decay offset by gamma income and financing of the delta position.

Vega: ν=C/σ\nu = \partial C / \partial \sigma

ν=Sτn(d1).\nu = S \sqrt{\tau}\, n(d_1).

Vega is positive for both calls and puts and is the same for both (put-call parity). It measures the option's sensitivity to volatility — the input that Black-Scholes treats as known but that must in practice be inferred from market prices (implied volatility).

Rho: ϱ=C/r\varrho = \partial C / \partial r

ϱcall=KτerτN(d2),ϱput=KτerτN(d2).\varrho_{\mathrm{call}} = K \tau e^{-r\tau} N(d_2), \qquad \varrho_{\mathrm{put}} = -K \tau e^{-r\tau} N(-d_2).

Rho is small for short-dated equity options but material for long-dated ones and for interest rate derivatives.


Validation

Several internal consistency checks confirm the formula is correct:

Put-call parity. CP=SKerτC - P = S - Ke^{-r\tau}. Verified directly from the formula.

Boundary behaviour. As τ0\tau \to 0: C(SK)+C \to (S - K)^+. As S0S \to 0: d1,d2d_1, d_2 \to -\infty, N(d1),N(d2)0N(d_1), N(d_2) \to 0, so C0C \to 0. As SS \to \infty: N(d1),N(d2)1N(d_1), N(d_2) \to 1, so CSKerτC \to S - Ke^{-r\tau} (intrinsic value).

Delta bounds. 0N(d1)10 \leq N(d_1) \leq 1 always; N(d1)1N(d_1) \to 1 as S/KS/K \to \infty (call becomes a forward), N(d1)0N(d_1) \to 0 as S/K0S/K \to 0 (call becomes worthless).


Limitations

Volatility Smile

The single most important empirical failure of Black-Scholes is the volatility smile: in the market, the implied volatility σ^(K,T)\hat{\sigma}(K, T) — defined as the unique σ\sigma such that BS(S,K,r,T,σ)=Cmarket\text{BS}(S, K, r, T, \sigma) = C_{\text{market}} — is not constant. It varies with strike and maturity.

In equity markets, the smile is typically a skew: implied vol is higher for low strikes (out-of-the-money puts) than high strikes, reflecting fear of large downside moves. In FX markets, the smile is more symmetric. This violates the core assumption of constant σ\sigma.

Discrete Hedging Error

Black-Scholes assumes continuous rebalancing. In practice, hedging occurs at discrete intervals {t0,t1,,tn}\{t_0, t_1, \ldots, t_n\}. The replication error per step is:

εk=12Γ(tk,Stk)[(ΔS)2σ2Stk2Δt],\varepsilon_k = \frac{1}{2}\Gamma(t_k, S_{t_k})\left[(\Delta S)^2 - \sigma^2 S_{t_k}^2 \Delta t\right],

which is non-zero when realised variance deviates from implied variance. Discrete hedging converts a continuous P&L of zero into a discrete P&L that fluctuates. The root-mean-square hedging error scales as σSΓΔt\sigma S \sqrt{\Gamma \cdot \Delta t}.

Jump Risk

GBM has continuous paths. Real asset prices exhibit jumps (earnings announcements, central bank decisions, geopolitical events). A jump of size JJ at time τ\tau creates an instantaneous P&L of 12ΓJ2\frac{1}{2}\Gamma J^2 — unhedgeable in the Black-Scholes framework. Jump-diffusion models (Merton 1976, Kou 2002) incorporate this at the cost of model complexity and incomplete markets.

Lognormal Tails

GBM implies lognormal STS_T, which has thin tails relative to empirical equity return distributions. Extreme losses are systematically underpriced. This is particularly dangerous for short volatility strategies and barrier option pricing near the barrier.

Constant Rate Assumption

For long-dated options (swaptions, caps, equity LTEPs), the constant-rr assumption introduces material pricing errors. Interest rate options require models where rates themselves are stochastic (see Hull-White, LMM).


Interview Angle

L1. Derive the Black-Scholes PDE via delta hedging. State the terminal and boundary conditions for a European call. Give the formula for Δ\Delta and explain its hedge ratio interpretation.

Derivation summary. Form Π=CCSS\Pi = C - \frac{\partial C}{\partial S} \cdot S. Apply Itô's lemma to dCdC, observe that choosing Δ=C/S\Delta = \partial C/\partial S eliminates all dWdW terms. No-arbitrage forces dΠ=rΠdtd\Pi = r\Pi\, dt. Equating gives the BS PDE. Terminal condition: C(T,S)=(SK)+C(T,S) = (S-K)^+. Boundary: C(t,0)=0C(t,0) = 0. Δcall=N(d1)\Delta_{\mathrm{call}} = N(d_1): holds N(d1)N(d_1) shares to replicate the call, dynamically rebalanced.

L2. Derive all five Greeks analytically. Reconcile the Black-Scholes PDE with the P&L decomposition Θ+12σ2S2Γ+rSΔrC=0\Theta + \frac{1}{2}\sigma^2 S^2 \Gamma + rS\Delta - rC = 0. Explain why the P&L of a delta-hedged call is 12S2Γ(σ^2σR2)dt\frac{1}{2}S^2\Gamma(\hat{\sigma}^2 - \sigma_R^2)dt where σ^\hat{\sigma} is implied vol and σR\sigma_R is realised vol.

P&L of a delta-hedged call. The BS PDE uses implied vol σ^\hat{\sigma}. If realised vol is σRσ^\sigma_R \neq \hat{\sigma}, the actual instantaneous P&L of the hedged portfolio is:

dVhedged=12S2Γ(σR2σ^2)dt.dV_{\text{hedged}} = \frac{1}{2}S^2 \Gamma \left(\sigma_R^2 - \hat{\sigma}^2\right) dt.

Long gamma earns when σR>σ^\sigma_R > \hat{\sigma} (bought cheap vol), loses when σR<σ^\sigma_R < \hat{\sigma}. This is the basis of volatility trading.

L3. Critique the no-arbitrage argument in Black-Scholes: under which conditions does the replicating portfolio exist, and what breaks down if markets are incomplete? Derive put-call parity from a cash-flow argument, not from the formula. Discuss the Breeden-Litzenberger result: what does the second derivative of the call price surface with respect to strike represent, and how does it reveal the risk-neutral density?

Breeden-Litzenberger. For call prices C(K)C(K) on a fixed maturity:

2CK2=erTpSTQ(K),\frac{\partial^2 C}{\partial K^2} = e^{-rT} p_{S_T}^{\mathbb{Q}}(K),

where pQp^{\mathbb{Q}} is the risk-neutral density of STS_T. This is derived by differentiating the pricing formula C(K)=erTEQ[(STK)+]C(K) = e^{-rT}\mathbb{E}^{\mathbb{Q}}[(S_T - K)^+] twice with respect to KK. The result implies that a full call price surface in strike and maturity uniquely determines the risk-neutral transition density — this is the conceptual foundation of Dupire's local vol.

Verify your understanding before moving on.