Darrell Duffie, Rui Kan
5 min
This paper introduces a multi-factor model of the term structure of interest rates that is both consistent with the absence of arbitrage and empirically convenient. Unlike many previous models that rely on unobservable latent variables, this framework uses the yields of zero-coupon bonds at fixed maturities as the primary state variables. By defining the model in this way, the state variables are directly observable from the current yield curve, which simplifies the interpretation and application of the model.
The authors utilize an affine term structure model, meaning that the yield of any zero-coupon bond is an affine (linear plus constant) function of the chosen basis set of yields. The paper provides the necessary and sufficient conditions for this representation, showing that the drift and diffusion functions of the underlying stochastic differential equation must also be affine. This structure ensures that bond prices can be solved using Riccati ordinary differential equations, making the model mathematically tractable.
A significant portion of the paper is dedicated to the practical implementation of the model. Because closed-form solutions are not always available, especially when incorporating stochastic volatility, the authors develop a finite difference algorithm to price interest rate derivatives. They introduce a change of variables that maps the state space into a compact set and orthogonalizes the system, which significantly improves the stability and convergence of the numerical solution. The authors also demonstrate how the model can be extended to include jump-diffusion processes to account for sudden shifts in interest rate expectations.
This work provides a unified framework that encompasses many existing term structure models as special cases. By grounding the model in observable yields rather than abstract latent factors, it offers a more direct link between theoretical pricing models and market data. The inclusion of numerical techniques for solving the associated partial differential equations makes this model a practical tool for practitioners involved in the pricing and risk management of fixed-income derivatives.
Alex: So the geometry of the state space itself is the obstacle. [[RP_SECTION:coordinate-transformation-techniques|Coordinate Transformation Techniques]]
Sam: Right. Their solution is a coordinate transformation — combining linear shifts and trigonometric functions — that maps the constrained domain onto a rectangular grid. That stabilizes the numerical solution and ensures boundary conditions are respected throughout. It doesn't distort the underlying dynamics, but it does add a layer of complexity to calibration: because the transformation is non-linear, you have to be precise when mapping observed yields back into the transformed space.
Alex: So it's a bridge — you preserve the theoretical purity of the affine framework, but use this numerical machinery to make it robust enough for actual trading constraints.
Sam: That's the right way to read it. The coordinate transformation is what turns the theoretical construct into a functional pricing engine. The authors demonstrate it remains tractable even for two-factor models, which is the real test if you need to price exotic swaptions in real time.
Alex: Where does this leave the framework when markets don't move smoothly? Affine diffusions assume continuous dynamics, but bond markets during liquidity shocks don't behave that way. [[RP_SECTION:incorporating-jump-diffusion|Incorporating Jump Diffusion]]
Sam: That's the gap they address in the extension. They incorporate a jump-diffusion process by adding an infinitesimal generator to the SDE that allows the state vector to experience sudden discrete shifts. The key to keeping this tractable is choosing jump distributions with known Laplace transforms — that choice preserves the exponential affine form, which is the engine of the model's efficiency.
Alex: So the Laplace transform condition is what keeps the ODE structure intact when a shock hits?
Sam: Exactly. You solve for the jump parameters the same way you solve for the diffusion parameters — a system of ODEs, not a fundamentally different computational problem. That means you get a realistic representation of tail risk without abandoning the computational speed the framework was built on.
Alex: Does adding jumps make calibration substantially harder?
Sam: It adds real complexity — you lose the standard staircase algorithm for inverting finite difference steps. But the authors show that for specific dimensions, numerical algorithms can still handle the jump distributions effectively. It's not a clean extension, but it's a workable one.
Alex: So the arc of the paper is essentially: start with a theoretically clean observable-factor model, build the numerical machinery to make it stable, then extend it toward the discontinuous dynamics that actually characterize stressed markets. [[RP_SECTION:framework-trade-offs-and-summary|Framework Trade-offs and Summary]]
Sam: That's the right summary. Each step involves a trade-off — tractability for flexibility, analytical cleanliness for empirical realism. What the framework buys you throughout is that the model stays anchored to observable market data. You're not propagating estimation error from latent factors into your derivative prices. For a desk pricing under time pressure, that anchor matters more than theoretical perfection.
Alex: Thanks for listening to ResearchPod.