---
title: Quadratic Term Structure Models
url: https://www.emergentmind.com/topics/quadratic-term-structure-models
type: topic
---

# Quadratic Term Structure Models

Searching arXiv for recent and foundational papers on quadratic term structure models.
Quadratic term structure models (QTSMs) are term-structure models in which the instantaneous short rate is specified as a quadratic function of an underlying factor process, and zero-coupon bond prices consequently take an exponential-quadratic form in the state variables. In continuous time, the factor dynamics are typically linear or Ornstein–Uhlenbeck type, while in discrete time they may be regime dependent; across these settings, the defining analytical mechanism is a system of Riccati-type equations or backward recursions for the coefficients in the bond-price exponent. The framework has been developed for zero-coupon bonds, FRAs, caplets, floors, swaptions, futures and forwards, as well as multi-curve, regime-switching, and polynomial jump-diffusion extensions [1410.1220][1512.03259][2212.04425].

## 1. Core specification and exponential-quadratic bond pricing

A standard continuous-time QTSM starts from a linear factor diffusion under a risk-neutral measure. In Hyndman–Zhou, the factor process \(X_t \in \mathbb R^n\) satisfies
\[
dX_t = (A X_t + B)\,dt + \sigma\,dW_t,
\]
and the short rate is
\[
r_t = r(X_t) = X_t^{\top}\Gamma X_t + R X_t + k,
\]
where \(\Gamma\) is symmetric and positive semidefinite, and \(k\) is chosen so that \(r(x)\ge 0\) for all \(x\). The zero-coupon bond price is then
\[
P(t,T)=E^Q\!\Bigl[\exp\!\bigl(-\int_t^T r(X_u)\,du\bigr)\Bigm|\mathcal F_t\Bigr].
\]
An equivalent generic specification in Lorig–Suaysom uses an Ornstein–Uhlenbeck-type state variable \(Y_t\in\mathbb R^d\),
\[
dY_t =(a + B\,Y_t)\,dt + \Sigma\,dW_t,
\]
with eigenvalues of \(B\) having negative real part, and a quadratic short rate
\[
R_t = q + Y_t^\top \Xi Y_t,\qquad q\ge 0,\;\Xi\succeq 0,\;\Xi_{ii}=1,
\]
so that \(R_t\ge 0\) [1410.1220][2212.04425].

The bond-pricing implication is the exponential-quadratic ansatz. In the \(Y\)-formulation,
\[
B^T_t
=\E_t\Big[e^{-\int_t^T R_s\,ds}\Big]
=\exp\Big(-F(t;T)-G(t;T)^\top Y_t-Y_t^\top H(t;T)Y_t\Big),
\]
where \((F,G,H)\) solve Riccati-type ODEs. In the \(X\)-formulation, the corresponding representation is
\[
P(t,T)
=\exp\!\Bigl[
X_t^{\top}R_2(t)X_t + R_1(t)X_t + R_0(t)
\Bigr].
\]
The matrix Riccati equation for \(R_2\), the linear equation for \(R_1\), and the scalar integral equation for \(R_0\) are obtained by matching quadratic, linear, and constant terms after applying Itô’s formula. Hyndman–Zhou further show that the solvability of these Riccati equations is linked to linear-quadratic control, and that under \(\Gamma\ge 0\) the matrix Riccati equation admits a unique continuous solution on \([0,T]\) [1410.1220].

A characteristic feature of QTSMs is that the state dependence of \(\log P(t,T)\) is quadratic rather than affine. This distinguishes them from affine term structure models even when both model classes are driven by linear Gaussian state processes. The additional quadratic term produces nonlinear Riccati structure and allows nontrivial curvature in yields while preserving closed-form or semi-closed-form bond pricing [1410.1220].

## 2. Gaussian exponentially-quadratic short-rate constructions

A particularly explicit specification is the Gaussian exponentially-quadratic multi-curve model of Crépey, Grbac, Ngor, and Skovmand. Under the risk-neutral OIS-bank-account measure \(Q\), the three-dimensional factor process \(\Psi_t=(\Psi_t^1,\Psi_t^2,\Psi_t^3)'\) satisfies
\[
d\Psi_t = F\Psi_t\,dt + D\,dW_t,
\]
with
\[
F=\operatorname{diag}(-b^1,-b^2,-b^3),\qquad
D=\operatorname{diag}(\sigma^1,\sigma^2,\sigma^3),
\]
and independent Wiener components, so that
\[
d\Psi_t^i = -b^i\Psi_t^i\,dt + \sigma^i\,dW_t^i,\qquad i=1,2,3.
\]
The short rate and spread are specified as second-order polynomials of Gaussian factors:
\[
r_t = A_r + B_r'\Psi_t + \Psi_t' C_r \Psi_t,\qquad
s_t = A_s + B_s'\Psi_t + \Psi_t' C_s \Psi_t,
\]
with
\[
A_r=0,\quad B_r=e_1,\quad C_r=e_2e_2',
\]
and
\[
A_s=0,\quad B_s=\kappa e_1,\quad C_s=e_2e_2'+e_3e_3'.
\]
The parameter \(\kappa\in\mathbb R\) introduces instantaneous correlation between \(r\) and \(s\) via the shared factor \(\Psi^1\) [1512.03259].

In this Gaussian EQ class, the OIS bond price is
\[
p(t,T)
=\exp\bigl[-A(t,T)-B(t,T)'\Psi_t-\Psi_t' C(t,T)\Psi_t\bigr].
\]
Because \(F\) and \(D\) are diagonal and the coefficient structure is sparse, only selected entries are nonzero. The explicit solution includes
\[
B^1(t,T)=\frac{1}{b^1}\bigl[1-e^{-b^1(T-t)}\bigr],
\]
and
\[
C^{22}(t,T)
=
\frac{2\bigl[e^{h^2(T-t)}-1\bigr]}
{2h^2+(2b^2+h^2)\bigl(e^{h^2(T-t)}-1\bigr)},
\qquad
h^2=\sqrt{4(b^2)^2+8(\sigma^2)^2},
\]
with
\[
A(t,T)=\int_t^T\Big[(\sigma^2)^2 C^{22}(u,T)-\frac12(\sigma^1)^2(B^1(u,T))^2\Big]\,du.
\]
Hence
\[
p(t,T)
=
\exp\bigl[-A(t,T)-B^1(t,T)\Psi_t^1-C^{22}(t,T)(\Psi_t^2)^2\bigr].
\]
The synthetic Libor-curve bond
\[
\bar p(t,T)
=
E^Q\!\Big[\exp\!\Big(-\int_t^T[r_u+s_u]\,du\Big)\Bigm|\mathcal F_t\Big]
\]
has the same exponential-quadratic structure, with
\[
\bar B^1(t,T)=(1+\kappa)B^1(t,T),\qquad
\bar C^{22}=C^{22},
\]
and \(\bar C^{33}\) determined by an analogous ODE [1512.03259].

This construction is notable because the factors remain Gaussian while the short rate and spread are quadratic. The resulting model is explicitly non-affine at the short-rate level, yet analytically tractable at the bond-pricing level.

## 3. Pricing FRAs, caplets, floors, and swaptions

In the multi-curve formulation, the single-curve FRA fair rate is
\[
R^{SC}(t;T,T+\Delta)=\frac{1}{\Delta}\bigl[\nu(t,T)-1\bigr],
\qquad
\nu(t,T)=\frac{p(t,T)}{p(t,T+\Delta)}.
\]
The post-crisis multi-curve FRA rate is
\[
R^{MC}(t;T,T+\Delta)=\frac{1}{\Delta}\bigl[\bar\nu(t,T)-1\bigr],
\qquad
\bar\nu(t,T)=E^{T+\Delta}[1/\bar p(T,T+\Delta)\mid\mathcal F_t].
\]
A central identity is the factorization
\[
\bar\nu(t,T)=\nu(t,T)\cdot Ad_t^{T,\Delta}\cdot Res_t^{T,\Delta},
\]
where
\[
Ad_t^{T,\Delta}
=
E^Q\!\Big[\frac{p(T,T+\Delta)}{\bar p(T,T+\Delta)}\Bigm|\mathcal F_t\Big]
=
E^Q\!\bigl[\exp(\tilde A+\kappa B^1\Psi_T^1+\bar C^{33}(\Psi_T^3)^2)\mid\mathcal F_t\bigr],
\]
and the deterministic residual is
\[
Res_t^{T,\Delta}
=
\exp\!\Big[-\kappa \frac{(\sigma^1)^2}{2(b^1)^3}
(1-e^{-b^1\Delta})(1-e^{-b^1(T-t)})^2\Big].
\]
The adjustment factor provides an explicit route from pre-crisis single-curve values to post-crisis multi-curve values [1512.03259].

For optional claims, the same Gaussian structure yields semi-analytic formulas. A caplet on \([T,T+\Delta]\) is priced as
\[
P^{cap}(t)
=
\Delta\,p(t,T+\Delta)\,
E^{T+\Delta}\!\big[(L(T;T,T+\Delta)-R)^+\mid\mathcal F_t\big]
=
p(t,T+\Delta)\,
E^{T+\Delta}\!\big[(1/\bar p(T,T+\Delta)-(1+\Delta R))^+\mid\mathcal F_t\big].
\]
Under \(Q^{T+\Delta}\), \(\Psi_T\) is Gaussian; the three-dimensional integral is reduced to explicit one-dimensional normal-cdf expressions by introducing the set
\[
M=\{(x,y):e^{\bar A +(1+\kappa)B^1x+C^{22}y^2}\le 1+\Delta R\},
\]
identifying the two roots \(z^\pm(x,y)\) of
\[
e^{\bar A +(1+\kappa)B^1x+C^{22}y^2+\bar C^{33}z^2}=1+\Delta R,
\]
and evaluating the resulting Gaussian integrals in closed form. Swaption pricing proceeds analogously under \(Q^{T_0}\), using a region
\[
\bar M=\{(x,y):g(x,y,0)\le h(x,y)\}
\]
and two critical values \(z^\pm(x,y)\) solving \(g(x,y,z)=h(x,y)\) [1512.03259].

Lorig and Suaysom develop a complementary asymptotic methodology for caplet implied volatilities in generic QTSMs. Under the \((T+\tau)\)-forward measure, the caplet forward price \(v(t,x,y)\) solves
\[
(\partial_t+\widetilde{\mathcal A}(t))v(t,x,y)=0,
\qquad
v(T,x,y)=\tau(e^x-e^k)^+,
\]
where \((X_t,Y_t)\) is Markov with generator
\[
\widetilde{\mathcal A}(t)=\sum_{|\alpha|\le 2} a_\alpha(t,x,y)\,\partial_{(x,y)}^\alpha.
\]
They expand both the generator and the solution in powers of a book-keeping parameter \(\epsilon\),
\[
\widetilde{\mathcal A}^\epsilon(t)=\sum_{n=0}^\infty \epsilon^n\mathcal A_n(t),
\qquad
v^\epsilon(t,z)=\sum_{n=0}^\infty \epsilon^n v_n(t,z),
\]
derive explicit formulas for \(\sigma_0,\sigma_1,\sigma_2\), and prove that in the short-maturity/parabolic regime
\[
\biggl|
\sigma_{\rm true}(t,x,y;T,\tau,k)
-
\sum_{n=0}^N \sigma_n(t,x,y;T,\tau,k)
\biggr|
=
\mathcal O\bigl((T-t)^{(N+1)/2}\bigr)
\quad\text{as }(T,k)\to (t,x).
\]
In numerical experiments for the one-factor Quadratic Ornstein–Uhlenbeck model
\[
dY_t=\kappa(\theta-Y_t)\,dt+\delta\,dW_t,\qquad R_t=q+Y_t^2,
\]
the second-order approximation captures the level, slope and curvature of the exact smile extremely well for all four maturities considered, especially near-ATM, and the absolute relative errors remain below a few tenths of one percent in a neighborhood of \(k-x=0\) and small \(T\) [2212.04425].

## 4. Regime-switching and polynomial jump-diffusion generalizations

A discrete-time extension is the regime-switching QTSM of Goutte. Time is indexed by \(k=0,1,\dots,T\), and the model includes a finite-state Markov chain \(\xi_k\) with time-dependent transition matrices \(Q_k=[q_{ij}^{(k)}]\), together with a regime-dependent factor process
\[
X_{k+1}=A_k(i)X_k+b_k(i)+\Sigma_k(i)\varepsilon_{k+1},
\qquad \varepsilon_{k+1}\sim N(0,I_d),
\]
when \(\xi_k=i\). The short rate in regime \(i\) is
\[
r_k=r(k,\xi_k,X_k)
=
X_k^T Q_k(\xi_k)X_k + L_k(\xi_k)^T X_k + C_k(\xi_k).
\]
The zero-coupon bond price remains exponential-quadratic:
\[
P(k,T\mid \xi_k=i,X_k)
=
\exp\!\Bigl\{
a_{k,T}(i)+b_{k,T}(i)^T X_k + X_k^T C_{k,T}(i)X_k
\Bigr\},
\]
with coefficient families initialized at maturity by
\[
a_{T,T}(i)=0,\qquad b_{T,T}(i)=0,\qquad C_{T,T}(i)=0,
\]
and propagated backward through explicit recursions involving
\[
M_j=\bigl[I_d-2\widetilde\Sigma^T C' \widetilde\Sigma\bigr]^{-1},
\qquad
\Delta_j=\widetilde\Sigma^T b' + 2\widetilde\Sigma^T C'\widetilde b.
\]
Regularity requires \(I_d-2\Sigma_k(i)^T C_{k+1,T}(j)\Sigma_k(i)\) to be invertible and positive-definite, ensuring finiteness of the Gaussian mgf [1305.2693].

A broader continuous-time extension is the polynomial jump-diffusion framework of Filipović–Willems. The state vector \(X=(X^1,\dots,X^d)'\) evolves on a convex domain \(E\subseteq \mathbb R^d\) with generator \(\mathcal L\) satisfying \(\mathcal L\,\mathrm{Pol}_2\subseteq \mathrm{Pol}_2\). A maximally general specification is
\[
dX_t
=
(b+BX_t)\,dt
+\Sigma_0\,dW_t
+\sum_{i=1}^d X_{t,i}\Sigma_i\,dW_t
+\int_Z(\delta_0(z)+\Delta(z)X_{t-})\,\tilde N(dt,dz),
\]
so the drift is affine and diffusion-jump coefficients are at most quadratic in \(X\). With a quadratic short rate
\[
r_t=r_0+r_1'X_t+X_t'R_2X_t,
\]
bond prices again take an exponential-quadratic form
\[
P(t,T)=\exp\{A(\tau)+B(\tau)'X_t+X_t' C(\tau)X_t\},
\qquad \tau=T-t,
\]
and \((A,B,C)\) satisfy generalized Riccati ODEs that incorporate diffusion and jump contributions. The same machinery supports a quadratic dividend specification,
\[
C_t=d_0(t)+d_1(t)'X_t+X_t' D_2(t)X_t,
\]
closed-form prices for dividend futures, and a present-value formula for the stock price. For derivatives with non-exponential-quadratic payoffs, the paper proposes a moment-based approximation: compute conditional moments, fit a maximum-entropy density matching those moments, and evaluate the price by numerical quadrature. The paper also records positivity and mgf-existence restrictions, including \(R_2\succeq 0\) for nonnegative short rates and light-tail conditions on the jump measure \(\nu\) [1803.02249].

## 5. Affine comparison and the statistical-consistency obstruction

QTSMs are closely related to affine term structure models (ATSMs), but the relationship is not identity. In an ATSM, one posits an affine short rate \(r_t=a^\top X_t+\alpha\) and obtains
\[
P(t,T)=\exp\bigl(B(t,T)^\top X_t + C(t,T)\bigr),
\]
with linear Riccati equations. In a QTSM, the exponent of the bond price contains a quadratic state term and the Riccati system acquires the nonlinear contribution
\[
\frac12(R_2+R_2^\top)\sigma\sigma^\top(R_2+R_2^\top).
\]
Hyndman–Zhou also note that QTSMs can be embedded into ATSMs on an augmented state space of dimension \(n+\tfrac12 n(n+1)\) by adjoining pairwise products \(\{X_iX_j\}\), and that if \(\Gamma=0\) then the model collapses to the affine case with \(R_2\equiv 0\) [1410.1220].

A more restrictive result is provided by Benth and Detering in the HJM-diffusion setting. They consider forward curves in Musiela parametrisation,
\[
f_t(x)=g(x,Y_t),
\]
with finite-dimensional realisation
\[
dY_t=\beta_t\,dt+\sigma_t\,dW_t,
\]
and absence of arbitrage encoded by the pointwise PDE
\[
\partial_x g(x,Y_t)
=
\nabla_y g(x,Y_t)\,\beta_t
+
\frac12\sum_{i,j=1}^d
[\sigma_t\sigma_t^\top]_{ij}\,
\partial_{y_i}\partial_{y_j}g(x,Y_t).
\]
They introduce the statistical consistency condition (SCC): for every constant deterministic \(\sigma\in\mathbb R^{d\times d}\), there exists a corresponding drift process \(\beta^\sigma\) such that the resulting model remains risk-neutral with the same parametrisation \(g\). Under finite-dimensional realisation, risk-neutrality, and SCC, they prove that \(g\) must be affine in the factor variable:
\[
g(x,y)=c(x)+u(x)^\top A(y),
\]
and that \(c\) and each component of \(u\) are quasi-exponential in \(x\) [2308.02246].

This theorem directly constrains genuine quadratic term-structure parametrisations. For
\[
g(x,y)=a(x)+b(x)^\top y + y^\top C(x)y,
\]
the gradient and Hessian are
\[
\nabla_y g(x,y)=b(x)+2C(x)y,\qquad
\partial_{y_i}\partial_{y_j}g(x,y)=2C_{ij}(x),
\]
so the Hessian is independent of \(y\) while the gradient is affine in \(y\). The paper shows that the compatibility demanded by SCC can hold for all \(y\) only if \(C(x)=0\). Hence any genuine quadratic term structure with \(C(x)\not\equiv 0\) fails SCC. The paper’s stated conclusion is that a non-degenerate quadratic term structure model cannot be both arbitrage-free and statistically consistent in the diffusion setting; possible remedies are to abandon SCC and fix a single diffusion law, or to allow highly constrained state-dependent volatility, which generically forces the model back toward affine form or degeneracy [2308.02246].

## 6. Tractability, calibration, and modeling trade-offs

The principal tractability advantage emphasized in the Gaussian EQ literature is that Gaussian factors remain Gaussian under any forward-measure change, since the change only induces a linear drift shift. This has several direct consequences: log-bond and option prices reduce to normal-cdf expressions, and calibration and simulation are simpler than in exponentially affine models driven by square-root factors, where more involved distributions appear [1512.03259].

At the same time, QTSMs involve model-specific admissibility conditions. In the Gaussian EQ multi-curve model, \(r_t\) and \(s_t\) can go negative with small probability, though they remain well behaved for realistic \(\sigma\)’s and \(b\)’s. The same framework allows a deterministic shift extension \(\phi_t\) to match the initial term structure exactly while preserving Gaussianity. For optional claims, conditions such as
\[
b^3 \ge \frac{\sigma^3}{\sqrt 2}
\]
are imposed so that certain quadratic forms inside Gaussian integrals stay positive [1512.03259].

Implementation strategies depend on the extension under consideration. In polynomial jump-diffusion specifications, the paper recommends choosing \(d\) small, for example \(2\)–\(4\) factors, numerically integrating the Riccati ODEs for \((A,B,C)\), inverting the bond and dividend-future formulas to recover \(X_0\), and matching option prices by moment fits. The same paper states that matching \(4\)–\(6\) moments often yields errors well below typical bid–ask spreads. In the regime-switching discrete-time case, the one-factor special case collapses to scalar recursions, and calibration may be performed with maximum-likelihood or the extended Kalman-filter / EM algorithm under hidden regimes [1803.02249][1305.2693].

Taken together, these results place QTSMs in a distinctive position within term-structure theory. They retain closed-form or semi-closed-form pricing through exponential-quadratic transforms and Riccati systems, accommodate multi-curve spreads, caplet and swaption pricing, regime changes, and polynomial jump-diffusion effects, and permit explicit implied-volatility asymptotics. At the same time, the affine-geometry theorem shows that, once statistical consistency is imposed in a diffusion-based HJM setting, genuine quadratic yield-curve parametrisations are no longer compatible with the required no-arbitrage geometry [1410.1220][2308.02246].

Source: https://www.emergentmind.com/topics/quadratic-term-structure-models