The HJM Framework
Interest Rate Modeling
Tutorial Contents
- Fixed Income Notation
- Fixed Income Probability Measures
- Multi-Currency Markets
- The HJM Framework
- Basic Fixed Income Instruments
- Caps, Floors and Swaptions
- CMS Swaps, Caps and Floors
- Bermudan Swaptions
- Structured Notes and Exotic Swaps
- Callable Libor Exotics and TARNs
- Volatility Derivatives
- Yield Curves and Curve Construction
- Managing Yield Curve Risk
- Vanilla Models With Local Volatility
- Vanilla Models With Stochastic Volatility
- Introduction to Term Structure Models
- One-Factor Short Rate Models I
- One-Factor Short Rate Models II
- Multi-Factor Short Rate Models
- The Libor Market Model I
- The Libor Market Model II
- Risk Management and Sensitivity Computations
- P&L Attribution
- Payoff Smoothing
- Vegas in the Libor Market Model
In this section, we assign dynamics to the major quantities we have defined thus far. Specifically, we look at the ZCB
prices P(t,T) and the instantaneous forward rate f(t,T). By ensuring that the dynamics of these quantities do not
admit arbitrage, we find that, when these values are driven by a finite number of Brownian Motions, the dynamics of our
various quantities must satisfy certain conditions. This is the crux of the Heath-Jarrow-Morton (HJM) framework.
The HJM framework prescribes the way in which the drift terms of our dynamics must relate to the diffusion terms in
order to satisfy no-arbitrage.
Bond Price Dynamics
Within the HJM framework, we focus on modeling an entire continuum of T-indexed bond prices P(t,T) for
t∈[0,T], starting from the known initial conditions P(0,T). We consider ZCB prices over a finite time-horizon,
i.e. T∈[0,T], T<∞, and assume that all models are generated by a d-dimensional
Brownian motion. Let W(t) represent a d-dimensional Brownian motion under the risk-neutral measure Q,
with numeraire β(t), the continuously compounded money market account. Define the discounted ZCB prices
P~(t,T)=β(t)P(t,T). In the absence of arbitrage P~(t,T) is a martingale under
Q, and the martingale representation theorem then implies that
dP~(t,T)=−P~(t,T)σP(t,T)⊺dW(t), where σP(t,T) is a d-dimensional stochastic process adapted to the filtration generated by W. We assume
that σP(t,T) is regular enough for P~(t,T) to be a square-integrable martingale, and since we have
the boundary conditions P(T,T)=1 for T∈[0,T], we impose the extra pull to par condition
σP(T,T)=0. Considering the SDE for β(t), using Ito's lemma shows that
P(t,T)dP(t,T)=P(t,T)d(β(t)P~(t,T))=P(t,T)P~(t,T)dβ(t)+β(t)dP~(t,T)=β(t)dβ(t)+P~(t,T)dP~(t,T)=r(t)dt−σP(t,T)⊺dW(t), where r(t) is the short rate. Equation
P(t,T)dP(t,T)=r(t)dt−σP(t,T)⊺dW(t) defines the class of d-dimensional HJM models.
Another application of Ito's lemma shows that forward bond prices P(t,T,T+τ)=P(t,T)P(t,T+τ) must
satisfy
dP(t,T,T+τ)=d(P(t,T)P(t,T+τ))=P(t,T)dP(t,T+τ)−P(t,T)2P(t,T+τ)dP(t,T)−P(t,T)2d⟨P(t,T),P(t,T+τ)⟩+P(t,T)3P(t,T+τ)d⟨P(t,T)⟩=P(t,T)P(t,T+τ)⋅(P(t,T+τ)dP(t,T+τ)−P(t,T)dP(t,T)−P(t,T)P(t,T+τ)d⟨P(t,T),P(t,T+τ)⟩+P(t,T)2d⟨P(t,T)⟩)=P(t,T,T+τ)⋅(−(σP(t,T+τ)−σP(t,T))⊺dW(t)−(σP(t,T+τ)−σP(t,T))⊺σP(t,T)dt). Hence, under the risk-neutral measure, forward bond prices satisfy the SDE
P(t,T,T+τ)dP(t,T,T+τ)=−(σP(t,T+τ)−σP(t,T))⊺(σP(t,T)dt+dW(t)). Since, under the T-forward measure QT P(t,T,T+τ) is a martingale, it must be driftless under the
QT measure, which immediately shows that the QT-Brownian motion is related to W(t) via
dWT(t)=dW(t)+σP(t,T)dt,(4.1) so that by Girsanov's theorem, the density process
η(t)=EtQ[dQdQT] for the measure shift
between QT and Q satisfies
η(t)dη(t)=−σP(t,T)⊺dW(t). Instantaneous Forward Rate Dynamics
While the above fully specifies the HJM dynamics of the ZCB and forward bond prices, the HJM framework is normally
expressed in terms of the instantaneous forward rates. Doing so eliminates the need to consider the short rate r
and also reveals a number of the fundamental properties of HJM models. By Ito's lemma, under the risk-neutral measure,
the quantity Y(t,T):=logP(t,T) satisfies
dY(t,T)=−σP(t,T)⊺dW(t)+g(t,T)dt, for some function g that is irrelevant for our purposes. Differentiating both sides of the above SDE by T and
taking the partial differential inside the differential term, we see that
f(t,T)=∂T∂Y(t,T) satisfies SDE
df(t,T)=μf(t,T)dt+σf(t,T)⊺dW(t), where
σf(t,T)=∂T∂σP(t,T).(4.2) Using the relationship (4.1), we find that
df(t,T)=μf(t,T)dt+σf(t,T)⊺(dWT(t)−σP(t,T)dt)=(μf(t,T)−σf(t,T)⊺σP(t,T))dt+σf(t,T)⊺dWT(t). Since f(t,T) is a QT martingale, the drift term must vanish. Furthermore, integrating (4.2) with respect
to T gives σP(t,T)=∫tTσf(t,u)du. These two results together show that
μf(t,T)=σf(t,T)⊺∫tTσf(t,u)du. This is the main result of Heath, Jarrow and Morton. Once the forward rate diffusion coefficients σf(t,T)
have been specified, the HJM model is fully specified. The HJM models take initial forward rates f(0,T)
as inputs, so that the models are automatically consistent with discount bond prices at time 0. This is true
irrespective of the choice of σf(t,T), which can be set freely.
This automatic calibration to bond prices is the main selling point of HJM models. However, there are a number of less
desirable features of HJM models. Particularly problematic is the dimensionality of the model. To describe bond
prices at some future time t, an entire continuum of forward rates {f(t,u):u∈[t,T]}. In practice, to make
this tractable requires either making assumptions on σf that permits a finite-dimensional Markovian
representation of the forward rate curve, or moving to a finite state space model that tracks forwards rates
spanning time-buckets of finite length. This second approach is explored in much more depth in the
Libor market model (LMM).
It's important to note that any arbitrage-free IR model set in a finite dimensional Brownian motion world must
be a special case of an HJM model. In particular, any such model corresponds to a particular choice of σf(t,T).
So, in some sense, the HJM framework underlies all of the models that we will cover.
Short Rate Processes
In principle specification of the short-rate process is sufficient to completely specify a full yield curve model.
Under the HJM framework, it follows that the short-rate r(t) under Q satisfies
r(t)=f(t,t)=f(0,t)+∫0tσf(u,t)⊺∫utσf(u,s)dsdu+∫0tσf(u,t)⊺dW(u). Clearly the process r(t) is not Markovian in general, as can be seen by focusing on the term
D(t)=∫0tσf(u,t)⊺dW(u), for which we have
D(T)=D(t)+∫tTσf(u,T)⊺dW(u)+[∫0tσf(u,T)⊺dW(u)−∫0tσf(u,t)⊺dW(u)]. Unless the square-bracketed terms are either non-random, or a deterministic function of D(t), we find that
EQ[D(T)∣D(t)]=EtQ[D(T)]. Examples of HJM Models
The Gaussian Model
In the Gaussian model, we assume that σP(t,T) is abounded deterministic function of t and T. It follows
that forward bond prices are then log-normal under both Q and QT. The forward rate process under
Q is
df(t,T)σf(t,T)=σf(t,T)⊺σP(t,T)dt+σf(t,T)⊺dW(t),=∂T∂σP(t,T). This implies that r(T) is Gaussian with Q moments
EtQ[f(T,T)]VartQ(f(T,T))=∫tTσf(u,T)⊺σP(u,T)du,=∫tTσf(u,T)⊺σf(u,T)du. The Gaussian short-rate model benefits from being tractable and permitting analytical price formulas for a
number of European options and futures. It is often a good starting point for the development of closed-form
approximations in other models. An example of a derivative admitting a closed-form solution under the Gaussian model is
an option on a ZCB:
Example (Option on a ZCB)
Consider a European call option with maturity S and strike price K on a T-bond. The payoff at time S is
given by
V(S)=(P(S,T)−K)+,S<T. Under the S-forward measure QS we know that
V(t)=P(t,S)EtS[(P(S,S,T)−K)+]. Applying the dynamics of P(t,S,T) we find that
V(t)=P(t,T)ϕ(d+)−KP(t,S)ϕ(d−), where
d±σ2=σlog(KP(t,S)P(t,T))±21σ2=∫tT∣σP(u,T)−σP(u,S)∣2du. Gaussian Models with Markovian Short Rate
Although tractable, the Gaussian model does not generally allow for a finite-dimensional Markovian representation,
and does not typically imply Markov behaviour of the short rate. By considering the special choice of
forward rate volatility
σf(t,T)=g(t)h(T), where h is a positive-valued real function and g(t)∈Rd, we can ensure that r(t) is Markovian. In
this case, we find that
σP(t,T)=∫tTσf(t,u)du=g(t)∫tTh(u)du(4.3) and
r(t)=f(0,t)+h(t)∫0tmf(t,u)du+h(t)∫0tg(u)⊺dW(u).(4.4) We find that the term
D(t)=∫0tσf(u,t)⊺dW(u)=h(t)∫0tg(u)⊺dW(u) is Markov, since
D(T)=h(T)∫0Tg(u)⊺dW(u)=h(t)h(T)D(t)+h(T)∫tTg(u)⊺dW(u). It quickly follows that r(t) itself is Markovian by differentiating (4.4). In fact, we have the following result.
Proposition (Markovian short rate)
In the d-dimensional Gaussian HJM model, when (4.3) holds, the short rate satisfies an SDE of the form
dr(t)=(a(t)−ξ(t)r(t))dt+σr(t)⊺dW(t), where ξ(t)∈R and σr(t)∈Rd are deterministic
functions of time, and
a(t)=∂t∂f(0,t)+ξ(t)f(0,t)+∫0tσf(u,t)⊺σf(u,t)du. Log-Normal HJM Models
To avoid the issue of negative forward rates that arises in Gaussian HJM models, it is tempting to consider forward
rates of the form
σf(t,T)=f(t,T)σ(t,T), where σ(t,T) is deterministic and bounded. In the T-forward measure,
df(t,T)=f(t,T)σ(t,T)⊺dWT(t), so that f(t,T) is log-normal. This leads to severe issues. Under Q we find that forward rates
explode to infinity with non-zero probability. Hence, all valuation formulas under Q lead to ZCB prices
being zero, implying obvious arbitrage opportunities.