Cheyette (1992) Markovian HJM Model Lab

Simulate Oren Cheyette's separable volatility framework reducing Heath-Jarrow-Morton forward rates to a 2D Markovian state space \((x, y)\). Analyze exact bond pricing, swaptions, and term structure dynamics.

Model Parameters

Cheyette Volatility & Reversion
yr⁻¹
Decay speed of the forward volatility kernel.
decimal
Annualized forward rate volatility (e.g. 0.012 = 1.20% or 120 bps).

Nelson-Siegel Market Curve \(f(0, T)\)

European Swaption Details

Monte Carlo Controls
Payer Swaption (bps)

--

--
Jamshidian Critical \(x^*\)

--

Root where Swap NPV = 0
Auxiliary State \(y(T_0)\)

--

Accumulated Variance
Volatility Half-Life

--

\(t_{1/2} = \ln(2)/\kappa\)
Steady-State \(y_\infty\)

--

\(\sigma^2 / (2\kappa)\)
5Y ZCB Martingale Test

--

MC vs Exact Discrepancy

Displays instantaneous forward curves \(f(t_0, t_0 + \tau)\) under neutral state \((x=0)\) and \(\pm 2\sigma_x\) interest rate shocks. Notice the upward convexity bend caused by the deterministic auxiliary variance \(y(t_0)\).

Jamshidian (1989) decomposition breaks the payer swaption into a weighted portfolio of European put options on individual zero-coupon bonds. Bars indicate basis point contribution to total swaption value; points indicate critical strikes \(K_i\).

Monte Carlo simulation of the Markovian short rate \(r(t) = f(0, t) + x(t)\) showing the median path and 5th-95th percentile envelopes across 500 stochastic paths.

Institutional Sensitivity Matrices
Mean Reversion \(\kappa\) vs Volatility \(\sigma\) on Swaption (bps)
\(\kappa \backslash \sigma\) 0.6% 1.0% 1.4% 1.8% 2.4%
Expiry \(T_0\) vs Tenor on Convexity Bias (bps)
\(T_0 \backslash \text{Tenor}\) 1.0Y 2.0Y 5.0Y 7.0Y 10.0Y
Institutional Framework & Mathematical Derivations
1. The Heath-Jarrow-Morton (1992) Path Dependency Dilemma

In the continuous Heath-Jarrow-Morton (1992) framework, the instantaneous forward rate for maturity \(T\) follows the stochastic differential equation:

\[df(t, T) = \alpha(t, T) dt + \sigma(t, T) dW(t)\]

Under the risk-neutral martingale measure \(\mathbb{Q}\), no-arbitrage uniquely fixes the drift \(\alpha(t, T)\) to:

\[\alpha(t, T) = \sigma(t, T) \int_t^T \sigma(t, u) du\]

Integrating over time reveals that the forward rate \(f(t, T) = f(0, T) + \int_0^t \alpha(s, T) ds + \int_0^t \sigma(s, T) dW(s)\) depends on the entire historical path of Brownian increments across the continuum of maturities. Consequently, general HJM models are infinite-dimensional and non-Markovian: evaluating path-dependent options (such as Bermudan swaptions) requires simulating the entire forward curve at every time step, creating severe computational bottlenecks.

2. The Cheyette (1992) Separable Volatility Breakthrough

Oren Cheyette (1992) discovered that if the volatility structure \(\sigma(t, T)\) factorizes into an exponential separable kernel:

\[\sigma(t, T) = \sigma(t) e^{-\kappa(T - t)}\]

where \(\kappa > 0\) is the constant mean-reversion parameter and \(\sigma(t) = \sigma\), the continuous continuum of forward rates collapses into an exact 2-dimensional Markov diffusion process governed by state variables \((x(t), y(t))\):

\[\begin{aligned} dx(t) &= (y(t) - \kappa x(t)) dt + \sigma dW(t), \quad x(0) = 0 \\ dy(t) &= (\sigma^2 - 2\kappa y(t)) dt, \quad y(0) = 0 \end{aligned}\]

Here, \(x(t)\) represents the stochastic deviation of the short rate from the initial market forward rate, while \(y(t)\) represents the deterministic accumulated variance and convexity adjustment:

\[y(t) = \frac{\sigma^2}{2\kappa}\left(1 - e^{-2\kappa t}\right), \quad \lim_{t \to \infty} y(t) = \frac{\sigma^2}{2\kappa}\]
3. Exact Zero-Coupon Bond Pricing & Instantaneous Forward Rates

With the state variables \((x(t), y(t))\) established, zero-coupon bond prices \(P(t, T)\) are given in exact closed form without numerical integration:

\[P(t, T) = \frac{P^M(0, T)}{P^M(0, t)} \exp\left( -G(t, T) x(t) - \frac{1}{2} G^2(t, T) y(t) \right)\]

where \(G(t, T) = \frac{1 - e^{-\kappa(T - t)}}{\kappa}\). Notice that at time \(t = 0\), since \(x(0) = 0\) and \(y(0) = 0\), the exponential equals 1, meaning that Cheyette identically fits any arbitrary market yield curve \(P^M(0, T)\) by construction.

The corresponding instantaneous forward rate curve is:

\[f(t, T) = f(0, T) + e^{-\kappa(T - t)} x(t) + e^{-\kappa(T - t)} G(t, T) y(t)\]
4. Jamshidian (1989) Exact European Swaption Decomposition

Because zero-coupon bond prices \(P(T_0, T_i; x(T_0), y(T_0))\) are strictly monotonic decreasing functions of the single random factor \(x(T_0)\), Farshid Jamshidian's (1989) decomposition theorem applies. At swaption expiry \(T_0\), the payer swaption payoff is:

\[\left(1 - \sum_{i=1}^n c_i P(T_0, T_i)\right)^+\]

where \(c_i = K_{\text{swap}} \Delta t\) for \(i < n\) and \(c_n = 1 + K_{\text{swap}} \Delta t\). By solving for the unique root \(x^*\) such that:

\[\sum_{i=1}^n c_i P(T_0, T_i; x^*, y(T_0)) = 1\]

each bond strike is defined as \(K_i = P(T_0, T_i; x^*, y(T_0))\). The swaption decomposes into an exact portfolio of European zero-coupon bond put options:

\[\mathbf{Swaption}(T_0, T_n, K_{\text{swap}}) = \sum_{i=1}^n c_i \mathbf{ZBO}_{\text{Put}}(T_0, T_i, K_i)\]
Cheyette Model Mastery Quiz

1. What mathematical property allows Cheyette (1992) to collapse the infinite-dimensional HJM forward curve into a 2D Markovian state space?

2. What economic and mathematical quantity does the auxiliary state variable \(y(t)\) represent?

3. How does the Cheyette model value European swaptions in closed form without two-dimensional numerical quadrature?

4. Why does the Cheyette model fit the initial market term structure \(P^M(0, T)\) identically without calibration error?

Frequently Asked Questions

The Cheyette model, published by Oren Cheyette in 1992 ("Term Structure Dynamics and Derivative Pricing"), solves the infamous infinite-dimensional path dependency of the Heath-Jarrow-Morton (1992) framework. By assuming that the forward rate volatility kernel is exponentially separable (\(\sigma(t, T) = \sigma(t) e^{-\kappa(T - t)}\)), Cheyette proved that the entire continuous forward curve collapses into a finite-dimensional Markov diffusion system governed by a state variable \(x(t)\) and an accumulated variance variable \(y(t)\).

General HJM models are non-Markovian because the instantaneous forward rate drift involves path-dependent double integrals of past volatilities across all maturities. By using a separable exponential volatility function, the integral factorizes into products of time-dependent functions and maturity-dependent functions. This algebraic factorization allows the infinite continuum of forward rates to be expressed strictly in terms of two Markov state variables: \(x(t)\) and \(y(t)\).

The auxiliary state variable \(y(t)\) represents the deterministic accumulated variance and convexity adjustment of the forward rate curve. It evolves according to \(dy(t) = (\sigma^2 - 2\kappa y(t)) dt\) with \(y(0) = 0\). In zero-coupon bond pricing, \(y(t)\) drives the quadratic term \(-0.5 \cdot G^2(t, T) \cdot y(t)\) in the pricing formula, capturing the fundamental convexity bias between forward rates and zero-coupon bond yields.

Because all zero-coupon bond prices \(P(T_0, T_i; x(T_0), y(T_0))\) in the Cheyette model are strictly decreasing, monotonic functions of the single stochastic driving factor \(x(T_0)\), Jamshidian's (1989) decomposition applies. By solving for the unique critical state \(x^*\) where the underlying swap value is zero, the swaption decomposes into an exact portfolio of European put options on individual zero-coupon bonds with strikes \(K_i = P(T_0, T_i; x^*, y(T_0))\), priced via Black-type analytical formulas.

Compared to standard 1-Factor Hull-White, Cheyette provides a direct forward-rate formulation that naturally handles term-dependent volatility humps, local volatility extensions, and multi-curve discounting without recalibrating theta(t). Compared to the high-dimensional Libor Market Model (BGM), Cheyette avoids carrying hundreds of forward rates across Monte Carlo paths, enabling ultra-fast pricing of Bermudan swaptions and path-dependent exotic interest rate derivatives.