# zero coupon bond pricing formula using Hull White

I am having some trouble to understand the derivation of the parameters of zero coupon pricing formula using Hull White. Specifically I am trying to understand how to get

$LogA(t,T)=[tR(0,t)-TR(0,T)]+B(t,T)&space;f^M&space;(0,t)-(\sigma&space;^2&space;(1-&space;exp^(-2at)))/4a\sigma*B(t,T)^2$ --[1]

where $R(t,T)$ is the continuously compounded zero spot rate: $R(t,T)=-logA(t,T)/(T-t)+&space;B(t,T)r(t)/(T-t)$

and $f^M(0,t)$ is the initial instantaneous forward rate curve

I understand the process to get below:

$-Log(A(t,T))=\int&space;\theta&space;(u)B(u,T)du+\sigma&space;^2/2a^2*(T-t-B(t,T))-\sigma&space;^2/4a*B(t,T)^2$ --[2]

and

$\theta&space;(t)&space;=&space;af^M(0,&space;t)+\partial&space;f^M(0,T)\partial&space;t+\sigma&space;^2/2a^2*(1-e^-2at))$ ---[3]

• What is $R(0, t)$? You need to provide more background information. – Gordon Jan 20 '17 at 16:26
• You may also search this site to see whether you can find some clues. – Gordon Jan 20 '17 at 16:42
• Thanks Gordon. I added some clarifications. I did do some search but could not find the answer. – user89877 Jan 20 '17 at 16:50
• Is there a difference between $log(A(t, T))$ and $Log(A(t, T))$? It appears that the left hand side of [2] should be $-log(A(t, T))$. – Gordon Jan 20 '17 at 17:01
• See whether this question can help you. – Gordon Jan 20 '17 at 17:03

Under the Hull-White interest rate model, the short rate $r_t$ satisfies a risk-neutral SDE of the form \begin{align*} dr_t = (\theta(t)-a r_t)dt+ \sigma dW_t. \end{align*} The price at time $t$ of a zero-coupon bond with maturity $T$ and unit face value is then given by \begin{align*} P(t, T) &= A(t, T) e^{-B(t, T) r_t}, \end{align*} where \begin{align*} B(t, T) &= \frac{1}{a}\Big(1-e^{-a(T-t)} \Big),\\ A(t, T) &= \exp\left(- \int_t^T \theta(u) B(u, T) du -\frac{\sigma^2}{2a^2}\big(B(t, T) -T+t\big)-\frac{\sigma^2}{4a}B(t, T)^2\right), \tag{1}\\ \theta(t) &= af^M(0, t) +\frac{\partial f^M(0, t)}{\partial t}+\frac{\sigma^2}{2a}\left(1-e^{-2at}\right), \\ \end{align*} and \begin{align*} f^M(0, t) &= -\frac{\partial \ln P(0, t)}{\partial t}. \end{align*} Note that, \begin{align*} \ln P(0, T) = -\int_0^T f^M(0, u) du. \end{align*} Moreover, we define the yield-to-maturity $R(t, T)$ by \begin{align*} R(t, T) &=\frac{-\ln P(t, T)}{T-t}. \end{align*} We show that \begin{align*} \ln A(t, T) &= \left[tR(0, t)-TR(0, T) \right] + B(t, T)f^M(0, t) -\frac{\sigma^2}{4a}\left(1-e^{-2at} \right)B(t, T)^2.\tag{2} \end{align*}
Note that, \begin{align*} & \ \int_t^T \frac{\partial f^M(0, u)}{\partial u}B(u, T) \\ =&\ f^M(0, u)B(u, T)\, \big|_t^T - \int_t^T f^M(0, u)\frac{\partial B(u, T)}{\partial u}du\\ =&\ -f^M(0, t)B(t, T) + \int_t^T f^M(0, u) e^{-a(T-u)} du\\ =&\ -f^M(0, t)B(t, T) - a\int_t^T f^M(0, u) B(u, T) du + \int_t^T f^M(0, u)du. \end{align*} That is, \begin{align*} &\ a\int_t^T f^M(0, u) B(u, T) du + \int_t^T \frac{\partial f^M(0, u)}{\partial T}B(u, T) \\ =&\ -f^M(0, t)B(t, T)+\int_t^T f^M(0, u)du\\ =&\ -f^M(0, t)B(t, T)- \ln P(0, T)+\ln(0, t)\\ =&\ -f^M(0, t)B(t, T)+ TR(0, T)-tR(0, t). \end{align*} Then \begin{align*} &\ \int_t^T \frac{\sigma^2}{2a}\left(1-e^{-2au}\right) B(u, T) du \\ =&\ \frac{\sigma^2}{2a^2}\int_t^T \left(1-e^{-2au}\right)\left(1-e^{-a(T-u)} \right)du\\ =&\ \frac{\sigma^2}{2a^2}\int_t^T \left(1-e^{-2au} - e^{-a(T-u)} + e^{-a(T+u)}\right)du\\ =&\ \frac{\sigma^2}{2a^2}\left[T-t+\frac{1}{2a}\left(e^{-2aT}-e^{-2at} \right) + \frac{1}{a}\left(1-e^{-a(T-t)} \right) -\frac{1}{a}\left(e^{-2aT}-e^{-a(T+t)}\right) \right] \\ =&\ \frac{\sigma^2}{2a^2}\left[\big(T-t- B(t, T)\big) + \frac{1}{2a}e^{-2at}\left(-e^{-2a(T-t)}-1 +2e^{-a(T-t)} \right)\right]\\ =&\ \frac{\sigma^2}{2a^2}\big(T-t- B(t, T)\big) - \frac{\sigma^2}{4a}e^{-2at}B(t, T)^2. \end{align*} Therefore, \begin{align*} &\ \ln A(t, T) \\ =&\ - \int_t^T \theta(u) B(u, T) du -\frac{\sigma^2}{2a^2}\big(B(t, T) -T+t\big)-\frac{\sigma^2}{4a}B(t, T)^2\\ =&\ f^M(0, t)B(t, T)+ tR(0, t)- TR(0, T)- \frac{\sigma^2}{2a^2}\big(T-t- B(t, T)\big) + \frac{\sigma^2}{4a}e^{-2at}B(t, T)^2\\ &\qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad -\frac{\sigma^2}{2a^2}\big(B(t, T) -T+t\big)-\frac{\sigma^2}{4a}B(t, T)^2\\ =&\ \big(tR(0, t)- TR(0, T)\big)+f^M(0, t)B(t, T) - \frac{\sigma^2}{4a}\left(1-e^{-2at} \right)B(t, T)^2, \end{align*} which is the required Claim $(2)$ above.