Institutional Fixed Income & Term Structure Engineering

Heath-Jarrow-Morton (1992) Forward Rate Lab

Simulate no-arbitrage instantaneous forward rate curve dynamics, multi-factor volatility structures, Monte Carlo zero-coupon bond martingale tests, and European bond option & swaption pricing under the landmark HJM (1992) framework.

Spot Short Rate r(0)

3.30%

Instantaneous f(0, 0)
10Y Forward Rate f(0, 10Y)

4.82%

Terminal yield anchor
Martingale Test MAE

1.42 bps

Zero-arbitrage error
Martingale Test RMSE

1.85 bps

Term structure dispersion
Model Parameters & Yield Curve

Initial Forward Curve $f(0, T)$ (Nelson-Siegel)

HJM Multi-Factor Volatility $\boldsymbol{\sigma}(t, T)$

Monte Carlo Simulation Grid

Derivative Pricing Calibration
Distribution of terminal forward curves $f(t_{\text{sim}}, T)$ simulated over 2,000 Monte Carlo paths with 5th, 25th, 50th, 75th, and 95th percentile uncertainty bands.
Discounted simulated bond price $\mathbb{E}^{\mathbb{Q}}[P(t, T)/B(t)]$ benchmarked against exact initial bond curve $P(0, T)$. Errors within standard error confirm the no-arbitrage condition.
Multi-factor term structure volatility decomposition: Factor 1 (Level / Mean-Reverting) and Factor 2 (Twist / Slope Hump) with the resulting endogenous no-arbitrage drift $\alpha(t, T)$.
Contingent Claims & European Option Valuation Benchmark
Zero-Coupon Call Option
HJM Monte Carlo: $0.0245
Black-76 Benchmark: $0.0248
Discrepancy: 3.1 bps
Zero-Coupon Put Option
HJM Monte Carlo: $0.0182
Black-76 Benchmark: $0.0180
Discrepancy: 2.4 bps
European Swaption ($1\text{Y} \times 4\text{Y}$)
Payer Swaption: 64.5 bps
Forward Bond $F$: 0.8412
Bond Vol $\sigma_P$: 3.82%
Sensitivity Grid: Mean-Reversion $\lambda_1$ vs Volatility $\sigma_1$ on 10Y Vol (bps)
Evaluates instantaneous 10-year forward rate volatility $\sigma_1(t, t+10) = \sigma_1 e^{-10 \lambda_1}$ in basis points.
Sensitivity Grid: Swaption Expiry vs Strike on Swaption Value (bps)
Evaluates European payer swaption premium across option expiry and strike rates for a 5-year underlying interest rate swap.
Mathematical Foundations of the Heath-Jarrow-Morton (1992) Framework
1. Instantaneous Forward Rates and the Initial Term Structure

Prior to Heath, Jarrow, and Morton (1992), interest rate modeling was dominated by short-rate models (e.g. Vasicek 1977, Cox-Ingersoll-Ross 1985, Ho-Lee 1986). A primary limitation was the difficulty of fitting the entire initial yield curve without calibrating artificial time-dependent parameters. HJM revolutionized mathematical finance by shifting the fundamental modeling variable to the instantaneous forward rate $f(t, T)$, defined as the rate contracted at time $t$ for riskless borrowing or lending over the infinitesimal interval $[T, T+dT]$:

$$f(t, T) = -\frac{\partial \ln P(t, T)}{\partial T} \quad \iff \quad P(t, T) = \exp\left(-\int_t^T f(t, u) \, du\right)$$

Because the model takes the market initial forward curve $f(0, T)$ as an exogenous initial condition, it fits the initial bond price term structure $P(0, T)$ by construction.

2. The HJM No-Arbitrage Drift Condition

Let the forward rate $f(t, T)$ evolve according to a multi-factor Itô stochastic differential equation under the empirical measure $\mathbb{P}$ or risk-neutral measure $\mathbb{Q}$:

$$df(t, T) = \alpha(t, T) \, dt + \sum_{i=1}^d \sigma_i(t, T) \, dW_i(t)$$

Under the risk-neutral measure $\mathbb{Q}$ (with the money market bank account $B(t) = \exp(\int_0^t r(s) ds)$ as numeraire), the discounted price of any zero-coupon bond $\tilde{P}(t, T) = P(t, T) / B(t)$ must be a local martingale. Applying Itô's Lemma to $P(t, T) = \exp(-\int_t^T f(t, u) du)$ and setting the drift of $\tilde{P}(t, T)$ to zero yields the celebrated HJM No-Arbitrage Drift Condition:

$$\alpha(t, T) = \sum_{i=1}^d \sigma_i(t, T) \int_t^T \sigma_i(t, u) \, du$$

Key Insight: In the risk-neutral world, there are zero free drift parameters. The drift of forward rates across all maturities is completely and uniquely dictated by the volatility structure $\boldsymbol{\sigma}(t, T)$ alone.

3. Separable Volatility and Markovian State Realizations

In a generic HJM specification, the spot short rate $r(t) = f(t, t)$ is non-Markovian—its evolution depends on the entire past history of Brownian increments $\int_0^t \sigma(s, t) dW(s)$. However, Carverhill (1994) and Ritchken & Sankarasubramanian (1995) established that the system admits a finite-dimensional Markovian state space if and only if the volatility structure is exponentially separable:

$$\sigma(t, T) = \sigma(t) e^{-\lambda(T - t)} \implies \alpha(t, T) = \frac{\sigma(t)^2}{\lambda} e^{-\lambda(T - t)} \left(1 - e^{-\lambda(T - t)}\right)$$

This specification corresponds exactly to the extended Vasicek / Hull-White (1990) one-factor short-rate model, proving that Hull-White is a special Markovian subset of the universal HJM framework.

4. Contingent Claims: Bond Options, Caps, and Swaptions

The HJM Monte Carlo engine prices contingent claims by simulating forward curve paths, accumulating the stochastic money market numeraire $B(t_E)$, and evaluating discounted expectation:

$$\text{Call}(t=0; t_E, T_B, K) = \mathbb{E}^{\mathbb{Q}}\left[ \frac{\max\left(P(t_E, T_B) - K, 0\right)}{B(t_E)} \right]$$

Under deterministic Gaussian volatility, this matches the analytical Black-76 bond option formula with integrated volatility $\sigma_P = \frac{\sigma}{\lambda}(1 - e^{-\lambda(T_B - t_E)})\sqrt{\frac{1 - e^{-2\lambda t_E}}{2\lambda}}$, providing an exact numerical benchmark for Monte Carlo validation.

HJM Term Structure Mastery Quiz
1. What fundamental feature distinguishes the Heath-Jarrow-Morton (HJM) model from earlier models like Vasicek or CIR?

2. Why are there zero free drift parameters in the risk-neutral HJM forward rate equation?

3. Under what condition does the HJM short rate $r(t) = f(t, t)$ become a Markov process?

4. How is the Monte Carlo martingale convergence test verified in the HJM simulator?
Frequently Asked Questions

The Heath-Jarrow-Morton (HJM) framework, published in Econometrica (1992) by David Heath, Robert Jarrow, and Andrew Morton, revolutionized term structure modeling by modeling the evolution of the entire instantaneous forward rate curve f(t, T) directly. Unlike short-rate models (e.g. Vasicek or CIR) that require calibrating artificial time-dependent drift parameters to match bond prices, HJM takes the initial forward curve f(0, T) as an exogenous input and fits market bond prices by construction.

Under the risk-neutral measure Q with the money market account as numeraire, the drift of instantaneous forward rates is uniquely determined by their volatility structure: α(t, T) = Σ σ_i(t, T) ∫_t^T σ_i(t, u) du. Because this condition is required to ensure that all discounted zero-coupon bond prices are martingales, there are zero free drift parameters in HJM; specifying the volatility structure σ(t, T) fully determines the risk-neutral drift.

In general HJM models, the short rate r(t) = f(t, t) depends on the entire past path and history of forward rate innovations. The system is Markovian in a finite number of state variables only if the volatility structure σ(t, T) satisfies special separability conditions, such as constant volatility (yielding the Ho-Lee model) or exponential decay σ(t, T) = σ exp(-λ(T - t)) (yielding the Hull-White / Vasicek model).

By applying Itô's Lemma to the zero-coupon bond pricing identity P(t, T) = exp(-∫_t^T f(t, u) du) and imposing the no-arbitrage condition, the drift of P(t, T) under the money market numeraire B(t) = exp(∫_0^t r(s) ds) equals r(t) P(t, T). Consequently, d(P(t, T)/B(t)) has zero drift, confirming that discounted bond prices are true martingales: E^Q[P(t, T)/B(t)] = P(0, T).

Hull-White (1990) can be viewed as a 1-factor Markovian special case of HJM with exponential volatility decay. However, because HJM relies on theoretical instantaneous forward rates f(t, T) that are unobservable in real trading, Brace, Gatarek, and Musiela (1997) extended HJM to discrete, observable compounded forward rates (LIBOR/SOFR), creating the LIBOR Market Model (BGM/LMM) and enabling direct Black-76 caplet calibration.
Explore Quantitative Interest Rate & Term Structure Simulations
LIBOR Market Model (BGM) Lab

Model discrete forward rates, Spot LIBOR numeraire drift stepping, and correlation decay across swaptions and ratchet caps.

Hull-White Short Rate Lab

Explore 1-factor and 2-factor Gaussian short-rate models with exact yield curve calibration and trinomial interest rate trees.

Black-Derman-Toy (BDT) Lab

Construct lognormal binomial interest rate trees fitting both term structure and yield volatilities without negative rates.

Black-Karasinski (BK) Lab

Construct calibrated lognormal short-rate trinomial trees via Arrow-Debreu state prices with negative convexity.