Hull–White model

For the rest of this article we assume only ${\displaystyle \theta }$ has t-dependence. Neglecting the stochastic term for a moment, notice that for ${\displaystyle \alpha >0}$ the change in r is negative if r is currently "large" (greater than ${\displaystyle \theta (t)/\alpha )}$ and positive if the current value is small. That is, the stochastic process is a mean-reverting Ornstein–Uhlenbeck process.

θ is calculated from the initial yield curve describing the current term structure of interest rates. Typically α is left as a user input (for example it may be estimated from historical data). σ is determined via calibration to a set of caplets and swaptions readily tradeable in the market.

When ${\displaystyle \alpha }$, ${\displaystyle \theta }$, and ${\displaystyle \sigma }$ are constant, Itô's lemma can be used to prove that

${\displaystyle r(t)=e^{-\alpha t}r(0)+{\frac {\theta }{\alpha }}\left(1-e^{-\alpha t}\right)+\sigma e^{-\alpha t}\int _{0}^{t}e^{\alpha u}\,dW(u),}$

which has distribution

${\displaystyle r(t)\sim {\mathcal {N}}\left(e^{-\alpha t}r(0)+{\frac {\theta }{\alpha }}\left(1-e^{-\alpha t}\right),{\frac {\sigma ^{2}}{2\alpha }}\left(1-e^{-2\alpha t}\right)\right),}$

where ${\displaystyle {\mathcal {N}}(\mu ,\sigma ^{2})}$ is the normal distribution with mean ${\displaystyle \mu }$ and variance ${\displaystyle \sigma ^{2}}$.

When ${\displaystyle \theta (t)}$ is time-dependent,

${\displaystyle r(t)=e^{-\alpha t}r(0)+\int _{0}^{t}e^{\alpha (s-t)}\theta (s)ds+\sigma e^{-\alpha t}\int _{0}^{t}e^{\alpha u}\,dW(u),}$

which has distribution

${\displaystyle r(t)\sim {\mathcal {N}}\left(e^{-\alpha t}r(0)+\int _{0}^{t}e^{\alpha (s-t)}\theta (s)ds,{\frac {\sigma ^{2}}{2\alpha }}\left(1-e^{-2\alpha t}\right)\right).}$

It turns out that the time-S value of the T-maturity discount bond has distribution (note the affine term structure here!)

${\displaystyle P(S,T)=A(S,T)\exp(-B(S,T)r(S)),}$

where

${\displaystyle B(S,T)={\frac {1-\exp(-\alpha (T-S))}{\alpha }},}$

${\displaystyle A(S,T)={\frac {P(0,T)}{P(0,S)}}\exp \left(\,-B(S,T){\frac {\partial \log(P(0,S))}{\partial S}}-{\frac {\sigma ^{2}(\exp(-\alpha T)-\exp(-\alpha S))^{2}(\exp(2\alpha S)-1)}{4\alpha ^{3}}}\right).}$

Note that their terminal distribution for ${\displaystyle P(S,T)}$ is distributed log-normally.

By selecting as numeraire the time-S bond (which corresponds to switching to the S-forward measure), we have from the fundamental theorem of arbitrage-free pricing, the value at time t of a derivative which has payoff at time S.

${\displaystyle V(t)=P(t,S)\mathbb {E} _{S}[V(S)\mid {\mathcal {F}}(t)].}$

Here, ${\displaystyle \mathbb {E} _{S}}$ is the expectation taken with respect to the forward measure. Moreover, standard arbitrage arguments show that the time T forward price ${\displaystyle F_{V}(t,T)}$ for a payoff at time T given by V(T) must satisfy ${\displaystyle F_{V}(t,T)=V(t)/P(t,T)}$, thus

${\displaystyle F_{V}(t,T)=\mathbb {E} _{T}[V(T)\mid {\mathcal {F}}(t)].}$

Thus it is possible to value many derivatives V dependent solely on a single bond ${\displaystyle P(S,T)}$ analytically when working in the Hull–White model. For example, in the case of a bond put

${\displaystyle V(S)=(K-P(S,T))^{+}.}$

Because ${\displaystyle P(S,T)}$ is lognormally distributed, the general calculation used for the Black–Scholes model shows that

${\displaystyle {E}_{S}[(K-P(S,T))^{+}]=KN(-d_{2})-F(t,S,T)N(-d_{1}),}$

where

${\displaystyle d_{1}={\frac {\log(F/K)+\sigma _{P}^{2}S/2}{\sigma _{P}{\sqrt {S}}}}}$

and

${\displaystyle d_{2}=d_{1}-\sigma _{P}{\sqrt {S}}.}$

Thus today's value (with the P(0,S) multiplied back in and t set to 0) is:

${\displaystyle P(0,S)KN(-d_{2})-P(0,T)N(-d_{1}).}$

Here ${\displaystyle \sigma _{P}}$ is the standard deviation (relative volatility) of the log-normal distribution for ${\displaystyle P(S,T)}$. A fairly substantial amount of algebra shows that it is related to the original parameters via

${\displaystyle {\sqrt {S}}\sigma _{P}={\frac {\sigma }{\alpha }}(1-\exp(-\alpha (T-S))){\sqrt {\frac {1-\exp(-2\alpha S)}{2\alpha }}}.}$

Note that this expectation was done in the S-bond measure, whereas we did not specify a measure at all for the original Hull–White process. This does not matter — the volatility is all that matters and is measure-independent.

Because interest rate caps/floors are equivalent to bond puts and calls respectively, the above analysis shows that caps and floors can be priced analytically in the Hull–White model. Jamshidian's trick applies to Hull–White (as today's value of a swaption in the Hull–White model is a monotonic function of today's short rate). Thus knowing how to price caps is also sufficient for pricing swaptions. In the even that the underlying is a compounded backward-looking rate rather than a (forward-looking) LIBOR term rate, Turfus (2020) shows how this formula can be straightforwardly modified to take into account the additional convexity.

Swaptions can also be priced directly as described in Henrard (2003). Direct implementations are usually more efficient.

However, valuing vanilla instruments such as caps and swaptions is useful primarily for calibration. The real use of the model is to value somewhat more exotic derivatives such as bermudan swaptions on a lattice, or other derivatives in a multi-currency context such as Quanto Constant Maturity Swaps, as explained for example in Brigo and Mercurio (2001). The efficient and exact Monte-Carlo simulation of the Hull–White model with time dependent parameters can be easily performed, see Ostrovski (2013) and (2016). An open-source implementation of the exact Monte-Carlo simulation following Fries (2016)[1] can be found in finmath lib.[2]

Even though single factor models such as Vasicek, CIR and Hull–White model has been devised for pricing, recent research has shown their potential with regard to forecasting. In Orlando et al. (2018,[3] 2019,[4][5]) was provided a new methodology to forecast future interest rates called CIR#. The ideas, apart from turning a short-rate model used for pricing into a forecasting tool, lies in an appropriate partitioning of the dataset into subgroups according to a given distribution [6]). In there it was shown how the said partitioning enables capturing statistically significant time changes in volatility of interest rates. following the said approach, Orlando et al. (2021) [7]) compares the Hull–White model with the CIR model in terms of forecasting and prediction of interest rate directionality.

• Vasicek model
• Cox–Ingersoll–Ross model
• Black–Karasinski model

Primary references

• John Hull and Alan White, "Using Hull–White interest rate trees," Journal of Derivatives, Vol. 3, No. 3 (Spring 1996), pp. 26–36
• John Hull and Alan White, "Numerical procedures for implementing term structure models I," Journal of Derivatives, Fall 1994, pp. 7–16.
• John Hull and Alan White, "Numerical procedures for implementing term structure models II," Journal of Derivatives, Winter 1994, pp. 37–48.
• John Hull and Alan White, "The pricing of options on interest rate caps and floors using the Hull–White model" in Advanced Strategies in Financial Risk Management, Chapter 4, pp. 59–67.
• John Hull and Alan White, "One factor interest rate models and the valuation of interest rate derivative securities," Journal of Financial and Quantitative Analysis, Vol 28, No 2, (June 1993) pp. 235–254.
• John Hull and Alan White, "Pricing interest-rate derivative securities", The Review of Financial Studies, Vol 3, No. 4 (1990) pp. 573–592.

Other references

• Hull, John C. (2006). "Interest Rate Derivatives: Models of the Short Rate". Options, Futures, and Other Derivatives (6th ed.). Upper Saddle River, N.J: Prentice Hall. pp. 657–658. ISBN 0-13-149908-4. LCCN 2005047692. OCLC 60321487.
• Damiano Brigo, Fabio Mercurio (2001). Interest Rate Models — Theory and Practice with Smile, Inflation and Credit (2nd ed. 2006 ed.). Springer Verlag. ISBN 978-3-540-22149-4.
• Henrard, Marc (2003). "Explicit Bond Option and Swaption Formula in Heath–Jarrow–Morton One Factor Model," International Journal of Theoretical and Applied Finance, 6(1), 57–72. Preprint SSRN.
• Henrard, Marc (2009). Efficient swaptions price in Hull–White one factor model, arXiv, 0901.1776v1. Preprint arXiv.
• Ostrovski, Vladimir (2013). Efficient and Exact Simulation of the Hull–White Model, Preprint SSRN.
• Ostrovski, Vladimir (2016). Efficient and Exact Simulation of the Gaussian Affine Interest Rate Models., International Journal of Financial Engineering, Vol. 3, No. 02.,Preprint SSRN.
• Puschkarski, Eugen. Implementation of Hull–White's No-Arbitrage Term Structure Model, Diploma Thesis, Center for Central European Financial Markets
• Turfus, Colin (2020). Caplet Pricing with Backward-Looking Rates., Preprint SSRN.
• Letian Wang, Hull–White Model, Fixed Income Quant Group, DTCC (detailed numeric example and derivation)