Tutorials › Stochastic Calculus › Girsanov Change of Measure in the Black-Scholes Model

The Black-Scholes Model · Tutorial 472 of 600

Girsanov Change of Measure in the Black-Scholes Model

Learn how exponential tilting changes the Brownian motion and converts the Black–Scholes stock drift from its physical value to the risk-free rate.

Advanced 11 min read

What You'll Learn

  • Define the constant market price of risk and the associated exponential density
  • Verify the density’s conditional-expectation martingale property and Novikov condition
  • Construct the equivalent measure on a finite time horizon
  • Show directly that the shifted process is Brownian under the new measure
  • Derive and check the stock’s risk-neutral drift
  • Distinguish the stock from its discounted value under the risk-neutral measure

From the Physical Drift to the Risk-Neutral Drift

In the Black–Scholes model, the stock’s drift under the physical probability measure describes its expected growth in that model. It need not equal the interest rate. A change of probability measure can alter the drift assigned to the Brownian driver while leaving the stock’s volatility unchanged. The key is to choose a positive density process whose weighting of Brownian paths is consistent over time.

Work on a filtered probability space \((\Omega,\mathcal{F},(\mathcal{F}_t)_{t\ge0},\mathbb{P})\) carrying a Brownian motion \(B_t\), and fix a finite horizon \(T>0\). Assume the Black–Scholes stock satisfies \(dS_t=\mu S_t\,dt+\sigma S_t\,dB_t\), with \(S_0>0\), constant \(\mu,r\in\mathbb{R}\), and constant volatility \(\sigma>0\). Here \(\mu\) is the drift under \(\mathbb{P}\), \(\sigma\) is the diffusion coefficient, and \(r\) is the constant interest rate. The strict assumption \(\sigma>0\) makes the market price of risk below well-defined.

Definition: The constant market price of risk is \(\theta=(\mu-r)/\sigma\). On \([0,T]\), define the density process $$ Z_t=\exp\left(-\theta B_t-\frac12\theta^2t\right),\qquad 0\le t\le T. $$ If \(Z_T\) has expectation one, the measure \(\mathbb{Q}\) on \(\mathcal{F}_T\) defined by \(d\mathbb{Q}=Z_T\,d\mathbb{P}\) is a probability measure. It is equivalent to \(\mathbb{P}\) when \(Z_T>0\) almost surely.

The negative sign in the exponent fixes the convention used here: the Brownian motion under \(\mathbb{Q}\) will be \(B_t^{\mathbb{Q}}=B_t+\theta t\). All stochastic integrals below are Itô integrals. Their integrands are adapted and progressively measurable, with the stated square-integrability condition whenever an \(L^2\) integral is used.

Why the Exponential Is a Valid Density

For a density process, expectation one at the terminal time is essential; the stronger conditional property \(\mathbb{E}[Z_t\mid\mathcal{F}_s]=Z_s\) for \(s\le t\) ensures that the densities at different times are consistent. Here both facts can be checked directly using independent Brownian increments. For a constant integrand, the integral \(\int_0^t\theta\,dB_s\) is defined with an adapted, progressively measurable integrand satisfying \(\mathbb{E}\int_0^t\theta^2\,ds=\theta^2t<\infty\), and equals \(\theta B_t\).

Proposition: For every finite \(T\), \(Z_t\) is a positive, true martingale under \(\mathbb{P}\) on \([0,T]\), and \(\mathbb{E}[Z_t]=1\). In particular, \(Z_T\) defines a probability measure \(\mathbb{Q}\) equivalent to \(\mathbb{P}\) on \(\mathcal{F}_T\).

Proof. Fix \(0\le s\le t\le T\), and put \(h=t-s\). Brownian increments are independent of \(\mathcal{F}_s\), with \(B_t-B_s\sim N(0,h)\). The density factors as \(Z_t=Z_s\exp(-\theta(B_t-B_s)-\tfrac12\theta^2h)\). The normal moment-generating function gives \(\mathbb{E}[\exp(-\theta(B_t-B_s))]=\exp(\tfrac12\theta^2h)\). Therefore \(\mathbb{E}[Z_t\mid\mathcal{F}_s]=Z_s\). In particular, taking \(s=0\) gives \(\mathbb{E}[Z_t]=Z_0=1\). Since \(Z_t\) is finite and strictly positive, \(Z_T\) integrates to one and is positive almost surely; hence the measure it defines is a probability measure equivalent to \(\mathbb{P}\). Integrability follows from positivity and \(\mathbb{E}[Z_t]=1\), so this is a true martingale, not merely a local martingale. \(\square\)

For comparison with a common sufficient criterion, Novikov’s condition for an integrand \(\theta_s\) on \([0,T]\) is \(\mathbb{E}[\exp(\tfrac12\int_0^T\theta_s^2\,ds)]<\infty\). It is a standard sufficient condition for the associated stochastic exponential to be a true martingale; its general proof involves technical localization and is not given here. In this constant case the condition is verified explicitly: \(\mathbb{E}[\exp(\tfrac12\int_0^T\theta^2\,ds)]=\exp(\tfrac12\theta^2T)<\infty\). The direct conditional-expectation argument above independently verifies the martingale property.

Worked Example: Constructing the Density

Take \(\mu=\tfrac15\), \(r=\tfrac1{20}\), and \(\sigma=\tfrac3{10}\), and consider the horizon \(T=2\). Find \(\theta\) and \(Z_2\), and check the density’s expectation and Novikov’s condition.

First calculate the drift difference and divide by the volatility: \(\mu-r=\tfrac15-\tfrac1{20}=\tfrac4{20}-\tfrac1{20}=\tfrac3{20}\), so \(\theta=(\tfrac3{20})/(\tfrac3{10})=\tfrac12\). Substitution into the density formula gives \(Z_2=\exp(-\tfrac12 B_2-\tfrac12(\tfrac12)^2\cdot2)=\exp(-\tfrac12B_2-\tfrac14)\).

Under \(\mathbb{P}\), \(B_2\sim N(0,2)\). Thus \(\mathbb{E}[Z_2]=e^{-1/4}\mathbb{E}[e^{-B_2/2}] =e^{-1/4}e^{(1/2)(1/4)(2)}=e^{-1/4}e^{1/4}=1\). For Novikov’s condition, \(\mathbb{E}[\exp(\tfrac12\int_0^2\theta^2\,ds)] =\exp(\tfrac12\cdot\tfrac14\cdot2)=e^{1/4}<\infty\).

Check: the directly computed expectation is exactly one, as required for a probability density. Also \(Z_2>0\) for every finite value of \(B_2\), so the resulting measure is equivalent to \(\mathbb{P}\).

The Brownian Motion Under the New Measure

Define \(B_t^{\mathbb{Q}}=B_t+\theta t\) for \(0\le t\le T\). The measure change does not make \(B_t\) itself a Brownian motion under \(\mathbb{Q}\); it is the shifted process \(B_t^{\mathbb{Q}}\) that is Brownian. The next argument verifies this directly for constant \(\theta\), rather than relying on the general Girsanov theorem.

Proposition: Under the measure \(\mathbb{Q}\) defined by \(Z_T\), the process \(B_t^{\mathbb{Q}}=B_t+\theta t\) is a Brownian motion relative to \((\mathcal{F}_t)\) on \([0,T]\).

Proof. Its paths are continuous and start at zero. To identify its increments, fix \(0\le s<t\le T\), write \(h=t-s\) and \(\Delta B=B_t-B_s\), and let \(g\) be any bounded measurable function. For an event \(A\in\mathcal{F}_s\), the density martingale property gives \(\mathbb{E}_{\mathbb{P}}[Z_T\mid\mathcal{F}_t]=Z_t\). Consequently, \(\mathbb{E}_{\mathbb{Q}}[\mathbf{1}_A g(\Delta B)] =\mathbb{E}_{\mathbb{P}}[\mathbf{1}_A Z_t g(\Delta B)]\). Factor \(Z_t=Z_s\exp(-\theta\Delta B-\tfrac12\theta^2h)\), then condition on \(\mathcal{F}_s\). Under \(\mathbb{P}\), \(\Delta B\) is independent of \(\mathcal{F}_s\) and has density \(\phi_h(x)=(2\pi h)^{-1/2}\exp(-x^2/(2h))\). Weighting this density by the increment factor gives \[ \phi_h(x)e^{-\theta x-\theta^2h/2} =(2\pi h)^{-1/2}\exp\left(-\frac{(x+\theta h)^2}{2h}\right). \] This is the normal density with mean \(-\theta h\) and variance \(h\), and integrates to one. The calculation with \(A\) shows that the conditional law of \(\Delta B\) given \(\mathcal{F}_s\) under \(\mathbb{Q}\) is this same law, independent of \(\mathcal{F}_s\). Therefore \(\Delta B+\theta h=B_t^{\mathbb{Q}}-B_s^{\mathbb{Q}}\) is independent of \(\mathcal{F}_s\) and has distribution \(N(0,h)\). Together with continuity, these increment properties establish that \(B^{\mathbb{Q}}\) is Brownian under \(\mathbb{Q}\). \(\square\)

Worked Example: Tilting One Brownian Increment

Use the parameters of the first example, so \(\theta=\tfrac12\). For an interval of length \(h=2\), find the distribution under \(\mathbb{Q}\) of \(\Delta B\) and of the shifted increment \(\Delta B^{\mathbb{Q}}\).

The density tilt for \(\Delta B=x\) is \(\exp(-\tfrac12x-\tfrac12(\tfrac12)^2\cdot2)=\exp(-\tfrac12x-\tfrac14)\). Multiplying the \(N(0,2)\) density by this factor gives \((4\pi)^{-1/2}\exp(-x^2/4-x/2-1/4) =(4\pi)^{-1/2}\exp(-(x+1)^2/4)\). Thus \(\Delta B\) has distribution \(N(-1,2)\) under \(\mathbb{Q}\). Since \(\Delta B^{\mathbb{Q}}=\Delta B+\theta h=\Delta B+1\), it has distribution \(N(0,2)\).

Check: the tilted mean is \(-\theta h=-(\tfrac12)(2)=-1\), and the variance remains \(h=2\). Adding \(\theta h=1\) changes the mean to zero without changing the variance.

Transforming the Stock Dynamics

By definition, \(dB_t^{\mathbb{Q}}=dB_t+\theta\,dt\), so \(dB_t=dB_t^{\mathbb{Q}}-\theta\,dt\). Substituting into the stock equation gives \(dS_t=(\mu-\sigma\theta)S_t\,dt+\sigma S_t\,dB_t^{\mathbb{Q}}\). The choice \(\theta=(\mu-r)/\sigma\) makes \(\mu-\sigma\theta=r\). Thus the stock drift under \(\mathbb{Q}\) is \(r\), while its diffusion coefficient remains \(\sigma S_t\).

Theorem: Under the assumptions above, the Black–Scholes stock satisfies \(dS_t=rS_t\,dt+\sigma S_t\,dB_t^{\mathbb{Q}}\) under \(\mathbb{Q}\). Its discounted value \(e^{-rt}S_t\) is a true martingale on \([0,T]\).

Proof. The drift identity follows from substitution and \(\mu-\sigma\theta=\mu-(\mu-r)=r\). To verify the discounted-martingale statement, the explicit stock solution under \(\mathbb{Q}\) is \(S_t=S_0\exp((r-\tfrac12\sigma^2)t+\sigma B_t^{\mathbb{Q}})\). This formula is also obtained by applying Itô’s formula to the exponential function: for \(f(t,x)=S_0\exp((r-\tfrac12\sigma^2)t+\sigma x)\), \(f_t=(r-\tfrac12\sigma^2)f\), \(f_x=\sigma f\), and \(f_{xx}=\sigma^2f\). The Itô correction is \(\tfrac12 f_{xx}\,dt=\tfrac12\sigma^2f\,dt\); hence the total \(dt\) coefficient is \(rf\), giving \(dS_t=rS_t\,dt+\sigma S_t\,dB_t^{\mathbb{Q}}\).

The stochastic-integral integrand \(\sigma S_t\) is adapted and progressively measurable, and it is square-integrable under \(\mathbb{Q}\) on every finite horizon because \(\mathbb{E}^{\mathbb{Q}}[S_t^2]=S_0^2e^{(2r+\sigma^2)t}\), whose integral over \([0,T]\) is finite. For \(0\le s\le t\le T\), the explicit solution and the independent \(N(0,t-s)\) increment of \(B^{\mathbb{Q}}\) give \(\mathbb{E}^{\mathbb{Q}}[S_t\mid\mathcal{F}_s] =S_s e^{(r-\sigma^2/2)(t-s)}e^{\sigma^2(t-s)/2} =S_s e^{r(t-s)}\). The conditional expectation is finite, so \(\mathbb{E}^{\mathbb{Q}}[e^{-rt}S_t\mid\mathcal{F}_s]=e^{-rs}S_s\). This proves the true-martingale claim, including its integrability; no step relies only on a local-martingale property. \(\square\)

Worked Example: Checking the Risk-Neutral Stock Equation

Let \(\mu=\tfrac3{10}\), \(r=\tfrac1{10}\), \(\sigma=\tfrac12\), and \(S_0=2\). Find \(\theta\) and write the stock solution under \(\mathbb{Q}\), then verify its drift by Itô’s formula.

Here \(\mu-r=\tfrac3{10}-\tfrac1{10}=\tfrac2{10}=\tfrac15\), so \(\theta=(\tfrac15)/(\tfrac12)=\tfrac25\). The transformed drift is \(\mu-\sigma\theta=\tfrac3{10}-(\tfrac12)(\tfrac25)=\tfrac3{10}-\tfrac2{10}=\tfrac1{10}=r\). The candidate solution is \(S_t=2\exp((\tfrac1{10}-\tfrac12(\tfrac12)^2)t+\tfrac12 B_t^{\mathbb{Q}}) =2\exp(-\tfrac1{40}t+\tfrac12 B_t^{\mathbb{Q}})\).

For this function \(f(t,x)=2e^{-t/40+x/2}\), the derivatives are \(f_t=-f/40\), \(f_x=f/2\), and \(f_{xx}=f/4\). The correction term is \(\tfrac12 f_{xx}\,dt=f\,dt/8\). Itô’s formula therefore gives a total drift \((-1/40+1/8)f\,dt=(4/40)f\,dt=f\,dt/10=r f\,dt\), and diffusion term \(f\,dB_t^{\mathbb{Q}}/2=\sigma f\,dB_t^{\mathbb{Q}}\).

Check: the computed drift coefficient is \(-\tfrac1{40}+\tfrac5{40}=\tfrac4{40}=\tfrac1{10}\), exactly \(r\). Also, for \(t\ge s\), the conditional expectation formula gives \(\mathbb{E}^{\mathbb{Q}}[S_t\mid\mathcal{F}_s]=S_s e^{(t-s)/10}\), confirming the growth rate independently of the differential calculation.

Common Pitfalls and Scope

Omitting the quadratic correction destroys the density. A tempting but incorrect choice is \(\widetilde Z_t=e^{-\theta B_t}\). For \(\theta\ne0\), its expectation is \(\mathbb{E}[\widetilde Z_t]=e^{\theta^2t/2}>1\) when \(t>0\), so it cannot be a probability density. The term \(-\tfrac12\theta^2t\) makes the expectation equal to one.

The sign determines the drift shift. With the convention \(B_t^{\mathbb{Q}}=B_t+\theta t\), substitution gives \(\mu-\sigma\theta\), not \(\mu+\sigma\theta\). Using the opposite sign in the density requires changing the definition of the shifted Brownian motion as well. Under the convention used here, the chosen \(\theta\) makes the drift exactly \(r\).

The measure \(\mathbb{P}\) describes the model before the change of measure; \(\mathbb{Q}\) is the equivalent risk-neutral measure on the specified finite horizon. The construction assumes constant coefficients and \(\sigma>0\). If \(\sigma=0\), the ratio defining \(\theta\) is unavailable: the stock has no Brownian risk to tilt, and a deterministic drift different from \(r\) cannot be changed by weighting Brownian paths. This tutorial establishes the measure change and the resulting stock dynamics; it does not derive prices of particular claims.

Check Your Understanding

In each question, assume the stated parameters are constant and the time horizon is finite.

  1. For \(\mu=\tfrac14\), \(r=\tfrac1{12}\), and \(\sigma=\tfrac12\), compute \(\theta\) and write \(Z_t\).
  2. For constant \(\theta\), use the normal moment-generating function to verify \(\mathbb{E}[Z_t\mid\mathcal{F}_s]=Z_s\).
  3. Under the convention \(B_t^{\mathbb{Q}}=B_t+\theta t\), what are the mean and variance under \(\mathbb{Q}\) of \(B_t-B_s\)?
  4. For \(\mu=\tfrac15\), \(r=\tfrac1{10}\), and \(\sigma=\tfrac14\), verify that the transformed drift is \(r\).
  5. Explain why proving that \(e^{-rt}S_t\) is a local martingale would not, by itself, establish the true-martingale conditional-expectation property.