Hull–White one-factor model
Fitting today’s term structure while modelling mean-reverting Gaussian short-rate dynamics
01Relate Ho–Lee, Vasicek and Hull–White short-rate dynamics.
02Derive the conditional Gaussian distribution of the mean-reverting factor.
03Explain how the time-dependent drift fits the initial curve.
04Price or calibrate caps/swaptions while diagnosing negative-rate and one-factor limitations.
Identify the state variables and the behavior they add.
Hull–White adds a deterministic time-dependent drift to an Ornstein–Uhlenbeck short-rate process. That shift fits today’s discount curve exactly, while mean reversion and volatility govern future rate distributions.
Curve fit and dynamic calibration are separate tasks.
Gaussian rates permit analytic bond formulas and negative outcomes.
One factor moves the entire curve through a single source of randomness, limiting decorrelation across tenors.
Ask which instruments can identify the dynamics.
Hull–White remains a transparent benchmark for callable structures, Bermudan exercise, exposure simulation and rate-option calibration.
caps and floors
European swaptions
callable bonds
Bermudan swaptions
Rate-option calibration must state normal/lognormal/shifted volatility convention, forward, annuity, expiry, tenor, strike, settlement and collateral curves.
Write the dynamics before interpreting parameters.
Hull–White dynamics
A time-inhomogeneous Gaussian short rate mean reverts at speed a with volatility σ.
Open in AnalyticsConditional variance
Mean reversion bounds long-horizon short-rate variance when a>0.
Open in AnalyticsFull derivation
Solving the mean-reverting Gaussian factor
Apply an integrating factor to the linear short-rate SDE and separate deterministic drift from stochastic innovation.
- 01
Multiply by the integrating factor
Apply e^{at} so the mean-reversion term becomes an exact differential.
- 02
Integrate from t to T
Solve for r_T as decayed current state plus deterministic drift integral and a Gaussian stochastic integral.
- 03
Compute the variance
Itô isometry converts the stochastic-integral variance into an ordinary integral.
- 04
Fit the initial term structure
Choose θ(t), or an equivalent deterministic shift φ(t), so model P(0,T) equals the observed discount curve at every maturity.
- 05
Calibrate option dynamics
Fit a and σ—or piecewise volatility—to a governed cap/swaption basket after the curve fit.
Hull–White separates exact initial-curve fit from a parsimonious one-factor Gaussian dynamics used for optionality and scenarios.
Inputs
r_t: short ratea: mean-reversion speed\sigma: short-rate volatility\theta(t): curve-fitting drift
Assumptions and limits
- Gaussian short rates permit arbitrarily negative values.
- One factor cannot reproduce independent curve twists.
- Parameters can be weakly identified and regime dependent.
Affine bond price
Zero-coupon bond prices are exponential-affine in the short rate.
Open in AnalyticsCalibrate, compute, and challenge the dynamics.
Fit the deterministic shift to the discount curve, calibrate dynamic parameters to a stated option-volatility convention and use analytic, tree or Monte Carlo engines appropriate to exercise style.
Calibrate with bid/ask-aware weights, multiple starts and stability diagnostics. Compare caps and swaption tenors not used in the fit.
Implementation with current QuantLib
Construct HullWhite from the accepted yield-curve handle, then use model-consistent engines and CalibrationHelpers for caps/swaptions. Recheck calibration after relinking the curve and inspect helper errors individually.
API authority: upstream QuantLib reference pinned in the source registry.06PYTHON IMPLEMENTATIONOpen the implementation and checks.
- Keep market conventions and quote lineage at the boundary.
- Solve curves and dynamics in framework-free deterministic kernels.
- Return residuals, state and sensitivities with every value.
- Test analytical limits, reconstruction identities and failure domains.
Exact Hull–White factor transition
Simulate the mean-zero Ornstein–Uhlenbeck factor with the exact Gaussian variance.
from __future__ import annotations import mathimport numpy as np def hw_transition(x: np.ndarray, a: float, sigma: float, dt: float, z: np.ndarray) -> np.ndarray: if a <= 0 or sigma < 0 or dt <= 0 or x.shape != z.shape: raise ValueError("invalid Hull-White transition") decay = math.exp(-a * dt) variance = sigma*sigma * (1.0 - math.exp(-2.0*a*dt)) / (2.0*a) return decay*x + math.sqrt(variance)*z rng = np.random.default_rng(7)z = rng.standard_normal(200_000)x = hw_transition(np.zeros_like(z), 0.08, 0.01, 1.0, z)theory = 0.01**2 * (1-math.exp(-0.16)) / 0.16assert abs(np.var(x) - theory) / theory < 0.015print(f"sample variance={np.var(x):.8f} | theory={theory:.8f}")Shock one parameter and trace the full response.
Hull–White dynamics laboratory
Move mean reversion, volatility and policy shock; inspect short-rate distribution, bond response and curve-factor decay.
Slow reversion: Shocks persist across the curve.
- shock persistence
- short-rate std dev
Use Left/Right or Up/Down arrows to inspect values; Home and End jump to the bounds.
Where the model meets the book.
“Hull–White fits today’s curve by construction; calibration quality is about tomorrow’s distribution and option basket.”
discount curve
swaption/cap volatility basket
volatility convention
mean reversion
model/engine discretisation
Fit curve deterministically, then calibrate dynamics with residual, parameter-bound and out-of-sample diagnostics.
RISKmean-reversion risk
volatility buckets
one-factor basis
exercise/model risk
- freeze curve/vol snapshot
- fit deterministic shift
- calibrate a/σ
- validate helpers
- run scenarios and hedges
Production failure modes
- normal/lognormal mismatch
- curve relink without recalibration
- parameter boundary
- tree time-grid misalignment
09MACRO CONNECTIONOpen the transmission channel.
Mean reversion and policy regimes
Policy shocks move the short end immediately; the mean-reversion parameter determines how quickly model scenarios pull that disturbance back toward the fitted term structure.
transmitsmoves current short-rate state
transmitscontrols persistence
transmitsrespond through affine dynamics
outputchanges exercise and convexity
10COMMON PITFALLSOpen the failure checklist.
Calling exact initial-curve fit a successful option calibration.
Comparing normal and lognormal swaption vols directly.
Ignoring negative-rate tails of a Gaussian model.
Using one factor for curve-basis risk it cannot span.
11SOURCES / FURTHER READINGOpen sources and continue the track.
Curve construction, multi-curve, short-rate and HJM lectures
Research map for term-structure theory and numerical experiments; all platform explanations and code are original.
- Source
- Financial Engineering: Interest Rates & xVA
- Author
- L. A. Grzelak
- Ref
- main
Monte Carlo, stochastic calculus and calibration lectures
Mathematical and numerical cross-reference for model dynamics and diagnostics.
- Source
- Computational Finance Course
- Author
- L. A. Grzelak
- Ref
- main
Current bootstrapping, interpolation, curve, model, cap/floor and swaption tests
Implementation authority for production object boundaries and regression-test patterns.
- Source
- QuantLib upstream
- Author
- QuantLib contributors
- Ref
- v1.42.1