04 · Local Volatility · Volatility Surface

The Dupire
Local Volatility Model.

A study of Dupire's local volatility model, from the forward Kolmogorov equation to the closed-form local volatility formula, and its role as the unique diffusion model that reprices the entire market smile exactly.

Context Personal study — Volatility modelling
Subject Local volatility surface calibration
Model Dupire (1994)
Method Forward PDE & implied volatility surface

A diffusion that fits the whole smile.

Black-Scholes assumes a single constant volatility, which cannot reproduce the smiles and skews observed on real option markets. Stochastic volatility models such as Heston introduce an extra source of randomness to capture these effects, at the cost of an incomplete market and a harder calibration problem.

Dupire (1994) took a different route: keep the underlying as a one-factor diffusion, as in Black-Scholes, but let the volatility depend deterministically on the current spot level and on time. This single generalization is enough, in principle, to reproduce every European option price quoted on the market, for every strike and maturity, exactly.

The local volatility diffusion.

Under the risk-neutral measure \(\mathbb{Q}\), the underlying price \(S_t\) is assumed to follow:

Local volatility dynamics
\[ dS_t = (r-q)S_t\,dt + \sigma_{\mathrm{loc}}(S_t,t)\,S_t\,dW_t^{\mathbb{Q}} \]

Unlike Black-Scholes, the diffusion coefficient \(\sigma_{\mathrm{loc}}(S,t)\) is a deterministic function of the spot level and time, rather than a constant. Unlike Heston, there is still only one source of randomness driving the process: the model remains a complete market, with a unique replicating hedge, while gaining enough flexibility to match an entire volatility surface.

The forward Dupire equation.

Instead of the backward Black-Scholes PDE in \((t,S)\), Dupire works with the price of a European call \(C(K,T)\) as a function of strike \(K\) and maturity \(T\), for a fixed valuation date and spot. Applying the Fokker-Planck (forward Kolmogorov) equation to the transition density of \(S_T\) yields the Dupire forward PDE:

Dupire forward equation
\[ \frac{\partial C}{\partial T} = \frac{1}{2}\sigma_{\mathrm{loc}}^2(K,T)K^2 \frac{\partial^2 C}{\partial K^2} - (r-q)K\frac{\partial C}{\partial K} - qC \]

This equation propagates option prices forward in maturity for every strike simultaneously, in contrast with the backward equation, which propagates a single option's price backward in time. This is what makes the forward equation the natural tool for extracting local volatility from an entire quoted surface in one pass.

Dupire's local volatility formula.

Rearranging the forward PDE gives a closed-form expression for the local variance in terms of market-observable call prices:

Local volatility (price form)
\[ \sigma_{\mathrm{loc}}^2(K,T) = \frac{ \dfrac{\partial C}{\partial T} + (r-q)K\dfrac{\partial C}{\partial K} + qC }{ \dfrac{1}{2}K^2 \dfrac{\partial^2 C}{\partial K^2} } \]

The numerator is the calendar-spread sensitivity of the option price, corrected for carry, and the denominator is proportional to the butterfly sensitivity, i.e. the risk-neutral density \(\partial^2 C/\partial K^2\) implied by the market. Local volatility is therefore fully determined, at every point of the surface, by the local slope in maturity and the local convexity in strike of the quoted call prices.

Expressing local volatility from implied volatility.

Market quotes are given in Black-Scholes implied volatility \(\sigma_{\mathrm{BS}}(K,T)\) rather than in raw call prices, which is both easier to interpolate and numerically better behaved. Using log-moneyness \(k=\ln(K/F(T))\), the local variance can be written directly in terms of the implied volatility surface:

Local variance (implied volatility form)
\[ \sigma_{\mathrm{loc}}^2(K,T) = \frac{ \sigma_{\mathrm{BS}}^2 + 2\sigma_{\mathrm{BS}}T \left( \dfrac{\partial \sigma_{\mathrm{BS}}}{\partial T} + (r-q)K\dfrac{\partial \sigma_{\mathrm{BS}}}{\partial K} \right) }{ \left(1+K d_1 \sqrt{T}\,\dfrac{\partial \sigma_{\mathrm{BS}}}{\partial K}\right)^2 + K^2 T \sigma_{\mathrm{BS}} \left( \dfrac{\partial^2 \sigma_{\mathrm{BS}}}{\partial K^2} - d_1\sqrt{T} \left(\dfrac{\partial \sigma_{\mathrm{BS}}}{\partial K}\right)^2 \right) } \]

This form, due to Gatheral and others, is the one typically implemented in practice: it converts a smoothly interpolated implied volatility surface directly into a local volatility surface, without ever differentiating the noisier raw price data.

Local volatility as a conditional expectation.

A key result, also due to Dupire, is that local volatility is the risk-neutral conditional expectation of the instantaneous variance, given that the underlying is at level \(K\) at time \(T\):

Conditional expectation representation
\[ \sigma_{\mathrm{loc}}^2(K,T) = \mathbb{E}^{\mathbb{Q}} \left[ \sigma_t^2 \,\middle|\, S_T=K \right] \]

This holds for any underlying diffusion or stochastic volatility process consistent with the market smile, not only for the local volatility model itself. It is what justifies describing local volatility as the market-implied average of instantaneous variance paths, rather than as a model of the true dynamics of volatility.

Numerical pitfalls and model limitations.

01

Differentiating a noisy surface

The Dupire formula requires first and second derivatives of quoted prices or implied volatilities with respect to strike and maturity. Raw market quotes are sparse and noisy, so the surface must first be smoothed by an arbitrage-free interpolation scheme before differentiating.

Ill-conditioned derivatives
02

Static, not dynamic, calibration

Local volatility matches the smile perfectly on a given day, but the implied forward smile it generates tends to flatten unrealistically fast as time passes, unlike the persistent smiles seen historically in equity markets.

Forward smile
03

Unstable hedge ratios

Because volatility is treated as a deterministic function of spot and time rather than a genuinely stochastic factor, the model's Greeks — particularly Vega and the vanna-volga sensitivities — can behave unrealistically for exotic and barrier-type payoffs.

Exotic hedging risk
04

Denominator instability near the wings

Where the implied risk-neutral density is close to zero, in particular far in the wings of the smile or at very short maturities, the denominator of the Dupire formula becomes small and the resulting local volatility can be numerically unstable or ill-defined.

Numerical stability

From the quoted smile to a diffusion surface.

01 A one-factor diffusion with a deterministic, state-and-time-dependent volatility function.
02 Derivation of the forward PDE via the Fokker-Planck equation, propagating prices across strikes and maturities in one pass.
03 Closed-form local volatility formula from either raw call prices or, more practically, the implied volatility surface.
04 Interpretation of local variance as the risk-neutral conditional expectation of instantaneous variance given the terminal spot level.
05 Structural limitations: unrealistic forward smile dynamics, sensitivity to surface noise, and unstable hedge ratios for exotic payoffs.
Local volatility remains the industry standard for exactly reproducing a quoted vanilla surface, and is often blended with a stochastic volatility model — a stochastic-local volatility (SLV) approach — to combine exact smile fit with more realistic forward-smile dynamics.