Interest Rate Models 19.1 Basic Elements in A Bond Market: T T T U T
Interest Rate Models 19.1 Basic Elements in A Bond Market: T T T U T
• Zero-coupon bond: Let B(t, τ ) be the time-t price of a bond with maturity τ (called τ -bond)
and the par value B(τ, τ ) = 1. The risk-neutral valuation principle implies
· µ Z τ ¶¯ ¸
¯
B(t, τ ) = EQ exp − ru du ¯¯ Ft .
t
−1
Y (t, τ ) = log B(t, τ ),
τ −t
which follows from
i.e. Y (t, τ ) is the constant compounding interest rate that brings B(t, τ ) to the par value.
∂
f (t, τ ) = − log B(t, τ ),
∂τ
which follows from
·Z τ ¸
B(t, τ ) exp f (t, u)du = 1,
t
1
i.e. f (t, τ ) is the instantaneous rate locked in at t for risk-free borrowing or lending at τ , with
f (t, t) = rt . Intuitively, f (t, τ ) should be interpreted as the interest rate over the infinitesimal
time interval [τ, τ + ∆τ ] as seen at time t.
To ease the notation, assume 0 ≤ τ ≤ T and consider a European option g(B(τ, T )) with
maturity τ written on a zero-coupon bond with maturity equal to the horizon T . The option g
can be hedged by a dynamic portfolio h = {(h0t , h1t ), t ∈ [0, T ]}, where h0t Bt and h1t B(t, T )
represent the time-t balance of the bank account and the time-t holding value of the zero-coupon
bond respectively. The time-t value of the portfolio is
Under certain regularity conditions, the time-t value of the option g is given by
· µ Z τ ¶ ¯ ¸
¯
Vt = EQ exp − ru du g(B(τ, T ))¯¯ Ft , (19.3)
t
which can be calculated in each problem with a specified interest rate model.
We will present three continuous-time single factor interest rate models in this section, whose
discrete-time versions appeared in Lecture 8. In each model, the SDE represents the risk neutral
dynamics under Q.
with positive constant parameters a, b and σ. Note that X = {Xt } with Xt = rt − b is an Ornstein-
Uhlenbeck (O-U) process satisfying the SDE
Hence X is a Gaussian process, which implies that Q(rt < 0) > 0 for every t ∈ [0, T ] — an
unreasonable behavior in practice. Pricing bond options would be easy due to nice properties of
the O-U process.
2
19.2.2 Cox-Ingersoll-Ross (CIR) model
√
drt = a(b − rt ) dt + σ rt dWt (19.5)
The following proposition shows that the CIR model is an important improvement of the Vasicek
model.
Proposition 19.1 There are three cases based on different ranges of the parameters:
Case 3 If 0 ≤ ab < σ 2 /2 and a < 0, then 0 < Q(Tx,0 < ∞) < 1 ∀ x > 0.
Proposition 19.1 provides us with a guideline for setting the parameters a, b, σ. {rt } in (19.5) is an
example of Feller’s square-root diffusions which enjoy great popularity in financial economics. It is
a Markov process with a closed-form transition density: the conditional density of ru given rt for
t < u is a scaled version of non-central χ2 distribution. Prices of the zero-coupon bond and some
bond options under a CIR model can be obtained with explicit expressions.
The Vasicek and CIR models, although very simple, do not model term structures directly.
When market term structure data {B(t, τ ), τ ∈ [t, T ]} are observed at time t, they are inconsistent
with the output generated by these two short rate models. In contrast, HJM proposes a general
mechanism that models term structures directly and takes observed term structures as inputs.
given the initial term structure f (0, τ ) = g(τ ), τ ∈ [0, T ], where α(·), σ(·) satisfy the HJM drift
condition
Z τ
α(t, τ ) = σ(t, τ ) σ(t, u) du (19.7)
t
3
to avoid arbitrage.
Note that in a general HJM model, the drift coefficient α(t, τ ) and volatility coefficient σ(t, τ ) are
not just functions of “current state” f (t, τ ), they also involve the history {f (u, τ ), 0 < u < t}. For
fixed τ , the process {f (t, τ ), 0 ≤ t ≤ τ } is non-Markovian. Therefore, the computation required
for the HJM model is extremely complicated, as a price to pay for greater modeling flexibility than
the Vasicek and CIR models.
Multi-factor versions of Vasicek, CIR and HJM are available that further improve the results
in empirical studies.