Black–Scholes PDE via Change of Numéraire¶
This derivation removes the physical drift \(\mu\) by unit normalization: expressing prices in units of the stock shifts the drift under Girsanov's theorem, and the martingale condition in the new units produces the PDE.
Orientation
This is the most abstract derivation in the chapter sequence. The high-level idea is: (i) take the stock \(S_t\) as the unit of account; (ii) Girsanov gives the density process \(Z_t = S_te^{-rt}/S_0\) and a new measure \(\mathbb{Q}^S\) under which \(V/S\) is a martingale; (iii) the PDE is read off by setting the drift of \(V/S\) to zero, then transforming back to \(V\). No replicating portfolio is constructed explicitly.
Recall (see § Risk-Neutral Pricing): under \(\mathbb{Q}\) the discounted stock \(e^{-rt}S_t\) is already a \(\mathbb{Q}\)-martingale, so \(Z_t = S_te^{-rt}/S_0\) defines a valid change of measure to \(\mathbb{Q}^S\). The general change-of-numéraire framework is developed in § Numéraire; this page is its concrete application to Black–Scholes.
Mechanism: Changing the Unit of Account in Two States¶
Numéraire change is fundamentally probability reweighting. Consider a one-period market with stock \(S_1 \in \{S^u, S^d\}\) and bond \(B_1 = 1+r\), both starting at \(S_0\) and \(B_0 = 1\). Pricing in dollars uses the risk-neutral measure \(\mathbb{Q}\), under which \(S_t/B_t\) is a martingale. Pricing in shares uses a new measure \(\mathbb{Q}^S\) obtained by reweighting:
The state where \(S_1\) is high gets up-weighted, the state where \(S_1\) is low gets down-weighted, by exactly the factor that makes \(1/S_t \cdot B_t = B_t/S_t\) a \(\mathbb{Q}^S\)-martingale. Asset prices in units of stock — \(V/S\) — are then the new "discounted prices." The continuous derivation below is the same operation in differential form: \(Z_t = S_te^{-rt}/S_0\) tilts \(\mathbb{Q}\) toward high-\(S_T\) paths, the new martingale condition is \(V/S\) rather than \(e^{-rt}V\), and the resulting PDE for the option price \(V\) turns out to be identical to the one obtained from the bond numéraire — the price is invariant to the chosen unit.
Setup¶
We work in the standard Black–Scholes model on a filtered probability space \((\Omega, \mathcal{F}, (\mathcal{F}_t)_{t \ge 0}, \mathbb{P})\) satisfying the usual conditions. Under the physical measure \(\mathbb{P}\):
with constant parameters \(\mu, \sigma > 0\), and a money market account \(B_t = e^{rt}\) with constant rate \(r\). Under the risk-neutral measure \(\mathbb{Q}\) (the \(B\)-martingale measure):
We take the stock \(S_t\) as numéraire. This is valid because \(S_t\) is strictly positive and pays no dividends, so its discounted price \(S_t e^{-rt}\) is a \(\mathbb{Q}\)-martingale. (When the stock pays dividends, the ex-dividend price alone is not a valid numéraire; one must use the total return process instead—see Exercise 5.) We derive the associated measure \(\mathbb{Q}^S\), the stock dynamics under \(\mathbb{Q}^S\), and ultimately the Black–Scholes PDE.
Step 1: Density Process and Girsanov Transformation¶
Radon–Nikodym Derivative¶
Recall (see § Numéraire): the density process from \(\mathbb{Q}\) to \(\mathbb{Q}^S\) is the ratio of normalized numéraires,
Dynamics of Z_t¶
Itô's product rule applied to \(Z_t = (e^{-rt}/S_0)S_t\) gives
after the \(r\,dt\) terms cancel. Thus \(Z_t = \exp(\sigma W^{\mathbb{Q}}_t - \tfrac{1}{2}\sigma^2 t)\) is a true \(\mathbb{Q}\)-martingale (constant integrand \(\sigma\) trivially satisfies Novikov).
Girsanov Transformation¶
Recall (see § Girsanov's Theorem): with density volatility \(\sigma\), the \(\mathbb{Q}^S\)-Brownian motion is
Step 2: Dynamics Under the Stock Measure¶
Stock Dynamics¶
Substitute \(dW^{\mathbb{Q}}_t = dW^S_t + \sigma\, dt\):
Under \(\mathbb{Q}^S\), the stock drift increases from \(r\) to \(r + \sigma^2\). This is natural: the density \(Z_T \propto S_T e^{-rT}\) up-weights high-\(S_T\) paths, increasing the drift under the new measure.
Verification: B_t / S_t Is a Q^S-Martingale¶
The defining property of \(\mathbb{Q}^S\) is that all asset prices normalized by \(S_t\) are martingales. We verify this for \(B_t / S_t\).
Apply Itô's formula to \(f(B_t, S_t) = B_t / S_t\), using \(dB_t = rB_t\, dt\) and the \(\mathbb{Q}^S\)-dynamics of \(S_t\):
Computing each term:
Collecting the drift:
The drift vanishes identically. Therefore
This driftless SDE shows \(B_t/S_t\) is a local \(\mathbb{Q}^S\)-martingale. To confirm it is a true martingale, solve explicitly: \(B_t/S_t = (B_0/S_0)\exp(-\sigma W_t^S - \frac{1}{2}\sigma^2 t)\), which is a Doléans-Dade exponential with constant integrand \(-\sigma\). Since \(\sigma\) is constant, \(\mathbb{E}^{\mathbb{Q}^S}[B_t/S_t] = B_0/S_0\) for all \(t\) (by the moment generating function of a Gaussian), so the process is a true \(\mathbb{Q}^S\)-martingale. \(\checkmark\)
Step 3: The Martingale Condition for the Normalized Price¶
Under \(\mathbb{Q}^S\), any traded asset price normalized by \(S_t\) is a martingale.
Assume \(V \in C^{1,2}([0,T) \times (0,\infty))\) with polynomial growth (sufficient for Itô's formula and for the local martingale below to be a true martingale, by the standard GBM moment bounds). Under \(\mathbb{Q}^S\), the normalized price
must be a \(\mathbb{Q}^S\)-martingale — this is no-arbitrage in the chosen numéraire. The martingale condition requires zero drift:
where the generator \(\mathcal{L}^S\) is:
This gives the PDE for the normalized price:
Intuition: changing units
The Radon–Nikodym derivative \(Z_T \propto S_T\) reweights paths according to the terminal stock value, tilting the probability measure toward high-\(S_T\) outcomes. In this tilted frame, normalized prices \(V/S\) are driftless (martingales), and the stock drift shifts to \(r + \sigma^2\). When we convert back to \(V\), the extra Itô terms cancel exactly with the drift shift. The precise content is the change-of-numéraire theorem, which guarantees the PDE for \(V\) is the same regardless of numéraire.
Step 4: Transform Back to V¶
The PDE \((\star)\) governs the normalized price \(u = V/S\) under the stock measure. Now we translate this back into a PDE for the option price \(V\) itself, by expressing the derivatives of \(u\) in terms of \(V\) and substituting.
Since \(V = S \cdot u\), express the derivatives of \(u\) in terms of \(V\):
Substitute into \((\star)\) and multiply through by \(S\). Collecting terms:
The \(\sigma^2\) terms from the stock-measure drift and the Itô correction cancel perfectly, leaving:
This is the Black–Scholes PDE—derived without explicitly constructing a replicating portfolio or hedging argument. (The underlying arbitrage-free pricing theory still relies on completeness and replication in the background; the point is that the derivation method avoids them.)
Pricing invariance
The normalized price \(u = V/S\) satisfies a PDE with drift coefficient \(r + \sigma^2\) (reflecting the stock-measure dynamics), but the option price \(V\) itself satisfies the standard Black–Scholes PDE with drift coefficient \(r\) (the risk-neutral drift). The cancellation of \(\sigma^2\) in the transformation is exact and guaranteed by the change-of-numéraire theorem: different numéraires give different PDEs for the normalized price, but the same PDE for \(V\).
Byproduct: Black–Scholes Formula as Two Probabilities¶
The stock-numéraire approach also yields an elegant decomposition of the Black–Scholes call price.
Derivation¶
For a European call with payoff \((S_T - K)^+\), the stock-numéraire pricing formula gives
Splitting the expectation:
For the second term, apply the abstract Bayes formula. Since \(Z_t = S_t e^{-rt}/S_0\) and \(Z_T = S_T e^{-rT}/S_0\):
Substituting \(S_0 Z_t = S_t e^{-rt}\):
Therefore:
Identification with N(d_1) and N(d_2)¶
Under \(\mathbb{Q}\), the stock drift is \(r\), so \(\ln S_T \mid \mathcal{F}_t \sim \mathcal{N}\!\left(\ln S_t + (r - \tfrac{1}{2}\sigma^2)\tau,\; \sigma^2\tau\right)\) where \(\tau = T - t\). Therefore
Under \(\mathbb{Q}^S\), the stock drift is \(r + \sigma^2\), so \(\ln S_T \mid \mathcal{F}_t \sim \mathcal{N}\!\left(\ln S_t + (r + \tfrac{1}{2}\sigma^2)\tau,\; \sigma^2\tau\right)\) (the Itô correction gives mean \(r + \sigma^2 - \frac{1}{2}\sigma^2 = r + \frac{1}{2}\sigma^2\)). Therefore
Combining:
The Black–Scholes formula emerges with a transparent probabilistic interpretation: \(N(d_1)\) is the probability of exercise under the stock measure \(\mathbb{Q}^S\), and \(N(d_2)\) is the probability of exercise under the risk-neutral measure \(\mathbb{Q}\).
Why \(d_1 \neq d_2\)
The two probabilities differ because they are computed under different measures. Under \(\mathbb{Q}^S\), paths where \(S_T\) is large are up-weighted (since \(Z_T \propto S_T\)), shifting the log-normal mean upward by \(\sigma^2 \tau\) relative to \(\mathbb{Q}\). This increases the exercise probability from \(N(d_2)\) to \(N(d_1) = N(d_2 + \sigma\sqrt{\tau})\).
Summary of the Derivation¶
The following diagram summarizes the logical flow:
flowchart LR
A["Money-market frame
(numéraire B, measure Q)
drift = r"] -->|"Radon–Nikodym
Z = Se⁻ʳᵗ/S₀"| B["Stock frame
(numéraire S, measure Q^S)
drift = r + σ²"]
B --> C["Martingale condition
u = V/S has zero drift"]
C --> D["PDE for u"]
D -->|"transform back (× S)
σ² terms cancel"| E["Black–Scholes PDE
(invariant)"]
The logical structure is:
- Choose numéraire: take \(N_t = S_t\) (the stock).
- Construct the density process: \(Z_t = S_t e^{-rt}/S_0\), verify \(dZ_t/Z_t = \sigma\, dW^{\mathbb{Q}}_t\).
- Apply Girsanov: \(W^S_t = W^{\mathbb{Q}}_t - \sigma t\) is a \(\mathbb{Q}^S\)-Brownian motion.
- Derive stock dynamics under \(\mathbb{Q}^S\): drift becomes \(r + \sigma^2\).
- Impose the martingale condition: \(u = V/S\) satisfies \(\partial_t u + \mathcal{L}^S u = 0\).
- Transform back to \(V\): the \(\sigma^2\) terms cancel, yielding the Black–Scholes PDE.
Recall (see § Self-Financing Replication): replication makes the strategy \((\alpha_t,\beta_t)\) explicit; here the same strategy is implicit, and the PDE is read off from the martingale condition for \(V/S\) rather than from \(\alpha_t = V_S\). The two methods rest on the same FTAP foundation — see § FTAP.
References¶
-
Geman, H., El Karoui, N., and Rochet, J.-C. (1995). Changes of numéraire, changes of probability measure, and option pricing. Journal of Applied Probability, 32(2), 443–458.
-
Shreve, S. E. (2004). Stochastic Calculus for Finance II: Continuous-Time Models. Springer.
-
Karatzas, I. and Shreve, S. E. (1991). Brownian Motion and Stochastic Calculus. 2nd edition, Springer.
-
Björk, T. (2009). Arbitrage Theory in Continuous Time. 3rd edition, Oxford University Press.
Exercises¶
Exercise 1. Verify that the density process \(Z_t = S_t e^{-rt}/S_0\) is a \(\mathbb{Q}\)-martingale by showing that \(dZ_t / Z_t = \sigma \, dW_t^{\mathbb{Q}}\) (i.e., the drift vanishes under \(\mathbb{Q}\)). Compute \(\mathbb{E}^{\mathbb{Q}}[Z_T]\) and confirm it equals 1.
Solution to Exercise 1
Under \(\mathbb{Q}\), the stock dynamics are \(dS_t = rS_t\, dt + \sigma S_t\, dW_t^{\mathbb{Q}}\). The density process is \(Z_t = S_t e^{-rt}/S_0\). Applying the product rule with \(e^{-rt}/S_0\) (a deterministic, bounded-variation factor):
Substituting the dynamics:
The \(r\, dt\) terms cancel exactly:
Therefore \(dZ_t / Z_t = \sigma\, dW_t^{\mathbb{Q}}\), which is driftless under \(\mathbb{Q}\), confirming \(Z_t\) is a \(\mathbb{Q}\)-local martingale.
Solving the SDE: \(Z_t = Z_0 \exp\!\bigl(\sigma W_t^{\mathbb{Q}} - \frac{1}{2}\sigma^2 t\bigr)\). Since \(Z_0 = S_0 e^0 / S_0 = 1\), we have \(Z_t = \exp\!\bigl(\sigma W_t^{\mathbb{Q}} - \frac{1}{2}\sigma^2 t\bigr)\). This is a Doléans-Dade exponential martingale. Computing the expectation:
where we used the moment generating function \(\mathbb{E}[e^{aW_T}] = e^{a^2 T/2}\). Since \(\mathbb{E}^{\mathbb{Q}}[Z_T] = 1\), the process \(Z_t\) is a true \(\mathbb{Q}\)-martingale and defines a valid probability measure \(\mathbb{Q}^S\).
Exercise 2. Under the stock measure \(\mathbb{Q}^S\), the process \(u_t = V(t, S_t)/S_t\) is a martingale. Apply Ito's lemma to \(u_t\) using the stock dynamics under \(\mathbb{Q}^S\) and show that setting the drift of \(u_t\) to zero yields the Black-Scholes PDE for \(V\).
Solution to Exercise 2
Under \(\mathbb{Q}^S\), the stock dynamics are \(dS_t = (r + \sigma^2)S_t\, dt + \sigma S_t\, dW_t^S\). The normalized price is \(u_t = V(t, S_t)/S_t\). Apply Itô's formula to \(u(t, S) = V(t, S)/S\). We need the partial derivatives:
The Itô expansion of \(u(t, S_t)\) gives:
Setting the drift to zero (the martingale condition under \(\mathbb{Q}^S\)):
Substituting the derivatives of \(u\) in terms of \(V\) and multiplying through by \(S\):
Simplifying:
The \(\sigma^2\) terms cancel in both brackets, yielding:
This is the Black–Scholes PDE.
Exercise 3. Explain why \(d_1 \neq d_2\) in the Black-Scholes formula from the change-of-numeraire perspective. Under \(\mathbb{Q}^S\), what is the distribution of \(\ln S_T\), and how does it differ from the distribution under \(\mathbb{Q}\)?
Solution to Exercise 3
The quantities \(d_1\) and \(d_2\) differ because they are computed under different probability measures, \(\mathbb{Q}^S\) and \(\mathbb{Q}\) respectively, and these measures assign different drifts to \(\ln S_T\).
Under \(\mathbb{Q}\) (risk-neutral measure): The stock dynamics are \(dS_t = rS_t\, dt + \sigma S_t\, dW_t^{\mathbb{Q}}\). By Itô's formula applied to \(\ln S_t\):
where \(\tau = T - t\). The exercise probability under \(\mathbb{Q}\) is:
Under \(\mathbb{Q}^S\) (stock measure): The stock dynamics are \(dS_t = (r + \sigma^2)S_t\, dt + \sigma S_t\, dW_t^S\). Applying Itô's formula to \(\ln S_t\):
The mean of \(\ln S_T\) under \(\mathbb{Q}^S\) is shifted upward by \(\sigma^2\tau\) compared to \(\mathbb{Q}\), because the Itô correction to the log gives \(r + \sigma^2 - \frac{1}{2}\sigma^2 = r + \frac{1}{2}\sigma^2\). The exercise probability under \(\mathbb{Q}^S\) is:
The difference \(d_1 - d_2 = \sigma\sqrt{\tau}\) arises because the Radon–Nikodym derivative \(d\mathbb{Q}^S/d\mathbb{Q} \propto S_T e^{-rT}\) up-weights paths where the stock ends high, shifting the mean of \(\ln S_T\) upward and increasing the in-the-money probability.
Exercise 4. The change-of-numeraire approach does not require constructing a self-financing portfolio. Compare the logical structure of this derivation with the delta-hedging derivation. In particular, identify which assumption (no-arbitrage, completeness, or Girsanov's theorem) plays the role that the hedging argument plays in the classical approach.
Solution to Exercise 4
Delta-hedging derivation: Constructs a portfolio \(\Pi = V - \Delta S\), applies Itô's formula under \(\mathbb{P}\), chooses \(\Delta = \partial V / \partial S\) to eliminate the \(dW\) term, then invokes the no-arbitrage condition \(d\Pi = r\Pi\, dt\). The key ingredients are: (i) a self-financing (or freeze-and-rebalance) portfolio construction, (ii) the ability to trade continuously, and (iii) the no-arbitrage principle.
Change-of-numéraire derivation: Chooses \(S_t\) as numéraire, constructs the density process \(Z_t = S_t e^{-rt}/S_0\), applies Girsanov's theorem to define \(\mathbb{Q}^S\), and requires \(V/S\) to be a \(\mathbb{Q}^S\)-martingale. No portfolio is constructed; the PDE emerges from setting the drift of the normalized price to zero.
The key structural difference is that the hedging argument in the classical approach is replaced by Girsanov's theorem in the change-of-numéraire approach. Specifically:
- In the delta-hedging approach, the choice \(\Delta = V_S\) eliminates risk and removes the physical drift \(\mu\) from the PDE. The no-arbitrage condition then determines the discount rate.
- In the change-of-numéraire approach, Girsanov's theorem absorbs the measure change from \(\mathbb{Q}\) to \(\mathbb{Q}^S\), modifying the stock drift from \(r\) to \(r + \sigma^2\). The martingale condition on \(V/S\) then yields the PDE, and the \(\sigma^2\) terms cancel when transforming back to \(V\).
Both approaches implicitly rely on completeness (to ensure uniqueness of the pricing measure) and no-arbitrage (to ensure existence of the measure). However, the change-of-numéraire approach makes the measure-theoretic structure explicit—Girsanov's theorem is the mathematical tool that plays the role of the hedging argument.
Exercise 5. Suppose the stock pays a continuous dividend yield \(q\). The ex-dividend stock price alone is not a valid numéraire (its discounted price is not a \(\mathbb{Q}\)-martingale). Instead, use the total return process \(\tilde{S}_t = e^{qt}S_t\) as numéraire, normalize by \(\tilde{S}_t\), and show that the Black–Scholes PDE becomes \(\frac{\partial V}{\partial t} + (r-q)S\frac{\partial V}{\partial S} + \frac{1}{2}\sigma^2 S^2 \frac{\partial^2 V}{\partial S^2} - rV = 0\).
Solution to Exercise 5
With a continuous dividend yield \(q\), the stock price satisfies \(dS_t = (r - q)S_t\, dt + \sigma S_t\, dW_t^{\mathbb{Q}}\) under \(\mathbb{Q}\). The total return process \(\tilde{S}_t = e^{qt} S_t\) satisfies \(d\tilde{S}_t = r\tilde{S}_t\, dt + \sigma\tilde{S}_t\, dW_t^{\mathbb{Q}}\), so its discounted price \(\tilde{S}_t e^{-rt}\) is a \(\mathbb{Q}\)-martingale, making \(\tilde{S}_t\) a valid numéraire.
Density process. The density from \(\mathbb{Q}\) to \(\mathbb{Q}^{\tilde{S}}\) is:
Since \(dZ_t/Z_t = \sigma\, dW_t^{\mathbb{Q}}\) (the \((r-q)\, dt\) terms cancel), \(Z_t\) is a \(\mathbb{Q}\)-martingale.
Girsanov transformation. Define \(W_t^{\tilde{S}} = W_t^{\mathbb{Q}} - \sigma t\), a \(\mathbb{Q}^{\tilde{S}}\)-Brownian motion.
Stock dynamics under \(\mathbb{Q}^{\tilde{S}}\). Substituting \(dW_t^{\mathbb{Q}} = dW_t^{\tilde{S}} + \sigma\, dt\):
Martingale condition. The normalized price under the numéraire \(\tilde{S}_t\) is \(u = V/\tilde{S}_t = e^{-qt}V/S\). This must be a \(\mathbb{Q}^{\tilde{S}}\)-martingale. We compute the partial derivatives of \(u(t, S) = e^{-qt}V(t,S)/S\):
Setting the drift of \(u(t, S_t)\) to zero, factoring out \(e^{-qt}\), and multiplying through by \(S\):
Collecting the \(V\) terms: \(-q - (r-q) = -r\). The \(\sigma^2\) terms cancel in both brackets, yielding:
This is the Black–Scholes PDE with continuous dividend yield \(q\). The discounting coefficient is \(-rV\) (not \(-(r-q)V\)) because the extra \(-qV\) from differentiating \(e^{-qt}\) in the normalization combines with \(-(r-q)V\) from the spatial terms to produce \(-rV\).
Exercise 6. Price a European call option directly using the stock measure \(\mathbb{Q}^S\). Show that the call price can be written as \(C = S\,\mathbb{Q}^S(S_T > K) - Ke^{-r\tau}\,\mathbb{Q}(S_T > K)\), evaluate both probabilities using the distributions of \(\ln S_T\) under each measure, and verify that the result is \(C = S\,\mathcal{N}(d_1) - Ke^{-r\tau}\,\mathcal{N}(d_2)\).
Solution to Exercise 6
The risk-neutral pricing formula for a European call with \(\tau = T - t\) is:
Second term. Under \(\mathbb{Q}\), \(\ln S_T \sim \mathcal{N}\!\left(\ln S + (r - \tfrac{1}{2}\sigma^2)\tau,\; \sigma^2\tau\right)\). Writing \(\ln S_T = \ln S + (r - \tfrac{1}{2}\sigma^2)\tau + \sigma\sqrt{\tau}\,Z\) with \(Z \sim \mathcal{N}(0,1)\):
where \(d_2 = \frac{\ln(S/K) + (r - \tfrac{1}{2}\sigma^2)\tau}{\sigma\sqrt{\tau}}\).
First term via change of numéraire. Using the density \(Z_T = S_T e^{-rT}/S_0\):
Under \(\mathbb{Q}^S\), the stock dynamics are \(dS_t = (r + \sigma^2)S_t\,dt + \sigma S_t\,dW_t^S\), so \(\ln S_T \sim \mathcal{N}\!\left(\ln S + (r + \tfrac{1}{2}\sigma^2)\tau,\; \sigma^2\tau\right)\). Therefore:
where \(d_1 = \frac{\ln(S/K) + (r + \tfrac{1}{2}\sigma^2)\tau}{\sigma\sqrt{\tau}} = d_2 + \sigma\sqrt{\tau}\).
Combining:
This is the Black–Scholes call formula. The change-of-numéraire viewpoint makes the structure transparent: \(\mathcal{N}(d_1)\) is the probability of finishing in the money under \(\mathbb{Q}^S\) (the stock measure), while \(\mathcal{N}(d_2)\) is the same probability under \(\mathbb{Q}\) (the money-market measure). The two \(d\)-values differ by \(\sigma\sqrt{\tau}\) because the stock drift is \(\sigma^2\) higher under \(\mathbb{Q}^S\).