Breeden-Litzenberger Formula

The theoretical derivation assumes the existence of a continuous spectrum of European option prices across all possible strikes. In practice, market quotes are discrete, requiring interpolation and extrapolation techniques to construct a smooth price function before the second derivative can be taken. Gaps in strike coverage introduce truncation errors in the tails of the extracted density.

The extracted density must be strictly non-negative and integrate to unity. This imposes specific shape constraints on the option price function:

• Butterfly Spread Arbitrage: The second derivative of the call price with respect to strike must be positive. A negative value implies a negative probability, violating no-arbitrage conditions.
• Calendar Spread Arbitrage: The total variance must be monotonically increasing with time to maturity.

The formula applies strictly to European-style options that can only be exercised at expiration. American options introduce an early exercise premium that contaminates the risk-neutral density extraction. For American options on dividend-paying stocks or indices, the extracted density reflects a pseudo-distribution biased by the early exercise boundary.

The derivation assumes a constant, known risk-free interest rate for discounting. In reality, the term structure of interest rates is stochastic. Using a flat rate introduces distortion in the density's location, particularly for longer-dated expirations where rate uncertainty is material. The formula must be adjusted to incorporate deterministic discount factors from the yield curve.

For underlying assets paying dividends, the formula requires the use of de-americanized or forward-adjusted option prices. Discrete dividend payments create jumps in the underlying price process, violating the diffusion assumption. The standard approach is to convert spot prices to forward prices and apply the formula to options on the forward, extracting the density of the forward asset price at maturity.

The second derivative is highly sensitive to the extrapolation method used beyond the quoted strike range. Common approaches include:

• Flat extrapolation: Assumes constant implied volatility beyond the wings, which may underestimate tail risk.
• Parametric tail fitting: Fits extreme value distributions (e.g., Generalized Pareto) to the tails.
• SVI or SABR extrapolation: Extends the parametric volatility smile smoothly into the tails. The choice of method significantly impacts the estimated probability of extreme market moves.

The core output of the Breeden-Litzenberger formula. This is the probability distribution of future asset prices implied by option prices, assuming all assets grow at the risk-free rate. It is not the physical real-world density, but a synthetic one used for pricing derivatives.

• Extracted via the second derivative of option price w.r.t. strike: ∂²C/∂K²
• Must integrate to 1 and be non-negative everywhere
• Used to price exotic payoffs by integrating the payoff function against this density

Mathematical constraints ensuring a volatility surface is free of static arbitrage. The Breeden-Litzenberger formula provides the diagnostic tool: a valid risk-neutral density requires ∂²C/∂K² ≥ 0 for all strikes.

• Butterfly arbitrage: Violated when the density becomes negative
• Calendar arbitrage: Violated when call prices decrease with time
• A surface failing these tests cannot be used for pricing without introducing internal inconsistencies

A forward partial differential equation that derives a unique local volatility surface from a continuum of traded option prices. It is the dynamic twin of the Breeden-Litzenberger formula.

• Breeden-Litzenberger gives the marginal distribution at a fixed expiry
• Dupire links these distributions across expiries to infer the diffusion coefficient
• Together, they form the theoretical backbone of local volatility modeling

An options strategy that directly monetizes the Breeden-Litzenberger relationship. A long butterfly—buying one low-strike call, selling two middle-strike calls, and buying one high-strike call—approximates ∂²C/∂K².

• The payoff at expiry converges to a Dirac delta function centered at the middle strike
• The price of a tightly-struck butterfly divided by the strike spacing approximates the risk-neutral density
• Used to trade discrete probability views on specific price levels

A theoretical primitive that pays $1 if a specific state occurs and $0 otherwise. The Breeden-Litzenberger formula reveals that option portfolios can synthesize these securities.

• The second derivative of a call price equals the price of an Arrow-Debreu security for that strike
• This insight bridges option pricing theory and general equilibrium state-preference theory
• Forms the foundation for pricing any contingent claim via replication

The process of fitting a parametric or non-parametric model to market-quoted option prices. The Breeden-Litzenberger formula acts as a constraint and validation tool during this process.

• Calibrated models must produce non-negative risk-neutral densities across all strikes
• W-shaped or oscillating densities indicate overfitting or noisy input data
• Techniques like regularization penalize violations of the no-arbitrage conditions implied by the formula