1 Introduction
Stochastic interest rate models play an important role in the modeling of financial markets. The literature essentially distinguishes between short rate models, forward rate models and market models. In the sequel, we give a brief survey on the different classes of term structure models. For more detailed information, the reader is referred to the respective research articles or the textbooks [7, 21] and [26].
Widely applied short rate models are for example the Vasicek model [38], the Hull–White model [29] or the Cox–Ingersoll–Ross (CIR) model [10]. In [38] and [29] the short rate process is modeled by a stochastic differential equation (SDE) of Ornstein–Uhlenbeck (OU) type driven by a Brownian motion (BM). As a consequence, the short rate process is normally distributed in these models and may become arbitrarily negative. Both features embody severe disadvantages with view on real-world market behavior, as the distribution of interest rate data frequently deviates from the normal distribution, while interest rates do not take arbitrarily large negative values in practice. In the recent years, there indeed appeared negative interest rates from time to time, but the negative values usually were small and stayed above some lower bound. However, in [10] the short rate is modeled by a so-called square-root process. This approach leads to a mean-reverting, strictly positive and chi-square distributed short rate process. In [6] the authors propose a time-homogeneous short rate model which is extended by a deterministic shift function in order to allow for negative rates and a perfect fit to the initially observed term structure. A very detailed overview on short rate models and their properties can be found in [7, 16] and [21]. In [26] and [33] short rate models in an extended multiple-curve framework are presented. The probably most famous forward rate model is the Heath–Jarrow–Morton (HJM) model proposed in [27]. Therein, the instantaneous forward rate process is modeled directly by an arithmetic SDE driven by a BM. In [4] the HJM model is extended to a jump-diffusion setup where the forward rate process is affected by both diffusion and random jump noise. HJM type models are also treated in [7] and Chapter 7 in [28]. In [12] and [26] HJM forward rate models in an extended multi-curve framework are discussed. The class of the so-called market models was introduced in [5]. For example, the popular LIBOR model belongs to this modeling class. In most cases, market models involve geometric SDEs such that the modeled interest rates usually turn out to be strictly positive. In order to allow the modeled rates also to take small negative values, shifted market model approaches have been proposed recently. Numerous properties of affine LIBOR models are provided in [31]. Market models are also presented in [7]. In [26] LIBOR models in an extended multi-curve framework are discussed. In [17] the authors propose a Lévy forward price model in a multi-curve setup which is able to generate negative interest rates. Term structure models which are driven by Lévy processes have also been proposed in [18–20].
In the present paper, we introduce a new pure-jump multi-factor short rate model which is bounded from below by a real-valued function of time which can be chosen arbitrarily. The short rate process is modeled by a deterministic function plus a sum of pure-jump zero-reverting Ornstein–Uhlenbeck processes. It turns out that the short rate is mean-reverting and that the related bond price formula possesses an affine representation. We also provide the dynamics of the related instantaneous forward rate, the latter being of HJM type. We further derive a condition under which the forward rate model can be market-consistently calibrated. The analytical tractability of our model is illustrated by the derivation of a plain-vanilla option price formula with Fourier transform methods. With view on practical applications, we make concrete assumptions on the distribution of the jump noises and show how explicit formulas can be deduced in these cases. We conclude the paper by presenting a post-crisis extension of our short and forward rate model.
The outline of the paper is as follows: In Section 2 we introduce our new pure-jump multi-factor short rate model which is bounded from below. Section 3 is dedicated to the derivation of related bond price and forward rate representations. Section 4 is devoted to option pricing. Section 5 contains guidelines for a practical application, while putting a special focus on possible distributional choices for the modeling of the involved jump noises. In Section 6 we consider a post-crisis extension of the proposed short and forward rate model.
2 A pure-jump multi-factor short rate model
Let $\left(\Omega ,\mathbb{F},\mathcal{F}={\left({\mathcal{F}_{t}}\right)_{t\in \left[0,T\right]}},\mathbb{Q}\right)$ be a filtered probability space satisfying the usual hypotheses, i.e. ${\mathcal{F}_{t}}={\mathcal{F}_{t+}}:={\cap _{s>t}}{\mathcal{F}_{s}}$ constitutes a right-continuous filtration and $\mathbb{F}$ denotes the sigma-algebra augmented by all $\mathbb{Q}$-null sets (cf. [30, 35]). Here, $\mathbb{Q}$ is a risk-neutral probability measure and $T>0$ denotes a fixed finite time horizon. In this setup, for arbitrary $n\in \mathbb{N}$ we define the stochastic short rate process $r={\left({r_{t}}\right)_{t\in \left[0,T\right]}}$ via
where $\mu \left(t\right)$ is a differentiable real-valued deterministic ${\mathcal{L}^{1}}$-function and ${X_{t}^{k}}$ constitute pure-jump zero-reverting Ornstein–Uhlenbeck (OU) processes satisfying the SDE
with deterministic initial values ${X_{0}^{k}}:={x_{k}}\ge 0$, constant mean-reversion velocities ${\lambda _{k}}>0$ and constant volatility coefficients ${\sigma _{k}}>0$. Herein, the independent, càdlàg, increasing, pure-jump, compound Poisson Lévy processes ${L_{t}^{k}}$ are defined by
where ${D_{k}}\subseteq {\mathbb{R}^{+}}:=\left]0,\infty \right[\subset \mathbb{R}$ denote jump amplitude sets and ${N_{k}}$ constitute Poisson random measures (PRMs). Note that the processes ${X_{t}^{k}}$ and ${L_{t}^{k}}$ always jump simultaneously, while ${X_{t}^{k}}$ decays exponentially between its jumps due to the dampening linear drift term appearing in (2.2). A typical trajectory of a Lévy-driven OU process is shown in Figure 15.1 in [9]. Further note that the background-driving time-homogeneous Lévy processes ${L_{t}^{k}}$ are increasing and thus, constitute so-called subordinators. Moreover, for all $k\in \left\{1,\dots ,n\right\}$ and $\left(s,z\right)\in \left[0,T\right]\times {D_{k}}$ we define the $\mathbb{Q}$-compensated PRMs
which constitute $\left(\mathcal{F},\mathbb{Q}\right)$-martingale integrators. Herein, the positive and σ-finite Lévy measures ${\nu _{k}}$ satisfy the integrability conditions
for an arbitrary constant $\varpi \in \mathbb{R}$ (cf. [9, 17]). For all $k\in \left\{1,\dots ,n\right\}$ and $t\in \left[0,T\right]$ we obtain
(2.4)
\[ d{\tilde{N}_{k}^{\mathbb{Q}}}\left(s,z\right):=d{N_{k}}\left(s,z\right)-d{\nu _{k}}\left(z\right)ds\](2.5)
\[ {\int _{{D_{k}}}}\left(1\wedge z\right)d{\nu _{k}}\left(z\right)<\infty ,\hspace{2em}{\int _{z>1}}{e^{\varpi z}}d{\nu _{k}}\left(z\right)<\infty \]
\[ {\mathbb{E}_{\mathbb{Q}}}\big[{L_{t}^{k}}\big]=t{\int _{{D_{k}}}}zd{\nu _{k}}\left(z\right),\hspace{2em}{\mathbb{Var}_{\mathbb{Q}}}\big[{L_{t}^{k}}\big]=t{\int _{{D_{k}}}}{z^{2}}d{\nu _{k}}\left(z\right)\]
both being finite entities due to (2.5) (cf. Section 1 in [17]). We remark that the currently proposed multi-factor short rate model (2.1) has been inspired by the electricity spot price model introduced in [2]. Arithmetic multi-factor models of this type have also been investigated in [28] and Section 3.2.2 in [3].Remark 2.1.
(a) Since ${L_{t}^{k}}$ is increasing and ${X_{t}^{k}}$ is zero-reverting from above, the function $\mu \left(t\right)$ is the mean-reversion floor or lower bound of the short rate process ${r_{t}}$, i.e. it holds ${r_{t}}\ge \mu \left(t\right)\hspace{2.5pt}\mathbb{Q}$-a.s. for all $t\in \left[0,T\right]$, while ${r_{t}}$ is mean-reverting from above to $\mu \left(t\right)$. Also note that the presence of a Brownian motion (BM) as driving noise in one of the processes ${X_{t}^{1}},\dots ,{X_{t}^{n}}$ would destroy the lower boundedness of ${r_{t}}$. In contrast to the presented pure-jump approach, it appears difficult to set up (lower-) bounded processes in arithmetic BM approaches. Moreover, we recall that negative rates have been observed in real-world post-crisis interest rate markets. Such scenarios can easily be captured by our model by choosing, e.g. $\mu \left(t\right)\equiv c$, where $c<0$ is an arbitrary constant. (In practical applications, it may happen that the floor function $\mu \left(t\right)$ needs to be readjusted, if interest rates evolve lower than anticipated. This issue has been discussed in [1] in the context of the SABR model.)
(b) Our pure-jump model (2.1) is able to generate short rate trajectories which closely resemble those stemming from common Brownian motion approaches, if we allow for small jump sizes only, i.e. ${D_{k}}=\left[{\epsilon _{1}^{k}},{\epsilon _{2}^{k}}\right]$ with small constants $0<{\epsilon _{1}^{k}}<{\epsilon _{2}^{k}}$. In this context, we emphasize that the well-established pure-jump variance gamma model is likewise able to generate suitable price trajectories, although there is neither any diffusion component involved (cf. Section 2.6.3 in [3], Table 4.5 in [9], Section 5.3.7 in [37]). On top of that, our pure-jump model might even provide more flexibility concerning the modeling of distributional properties than common BM approaches, since we are able to implement tailor-made distributions via an appropriate choice of the Lévy measures ${\nu _{k}}$ which fit the empirical behavior of the rates in a best possible manner. This topic is further discussed in Section 5 below. For instance, (generalized) inverse Gaussian, tempered stable or gamma distributions might embody suitable choices (recall Appendix B.1.2 on p. 151 in [37]). We finally recall that a model of the type (2.1) has been fitted to real market data in [2] (yet in an electricity market context).
For a time partition $0\le t\le s\le T$ the solution of (2.2) under $\mathbb{Q}$ can be expressed as
where we used (2.3). The representation (2.6) implies
where $0\le t\le T$. For all $t\in \left[0,T\right]$ we next define the historical filtration
(2.6)
\[ {X_{s}^{k}}={X_{t}^{k}}{e^{-{\lambda _{k}}\left(s-t\right)}}+{\sigma _{k}}{\int _{t}^{s}}{\int _{{D_{k}}}}{e^{-{\lambda _{k}}\left(s-u\right)}}zd{N_{k}}\left(u,z\right)\](2.7)
\[ {X_{t}^{k}}={x_{k}}{e^{-{\lambda _{k}}t}}+{\sigma _{k}}{\int _{0}^{t}}{\int _{{D_{k}}}}{e^{-{\lambda _{k}}\left(t-s\right)}}zd{N_{k}}\left(s,z\right)\]Proposition 2.2.
For $0\le u\le t\le T$ we have
\[\begin{array}{l}\displaystyle {\mathbb{E}_{\mathbb{Q}}}\left({r_{t}}|{\mathcal{F}_{u}}\right)=\mu \left(t\right)+{\sum \limits_{k=1}^{n}}\left({X_{u}^{k}}{e^{-{\lambda _{k}}\left(t-u\right)}}+{\sigma _{k}}\frac{1-{e^{-{\lambda _{k}}\left(t-u\right)}}}{{\lambda _{k}}}{\int _{{D_{k}}}}zd{\nu _{k}}\left(z\right)\right),\\ {} \displaystyle {\mathbb{Var}_{\mathbb{Q}}}\left({r_{t}}|{\mathcal{F}_{u}}\right)={\sum \limits_{k=1}^{n}}{\sigma _{k}^{2}}\frac{1-{e^{-2{\lambda _{k}}\left(t-u\right)}}}{2{\lambda _{k}}}{\int _{{D_{k}}}}{z^{2}}d{\nu _{k}}\left(z\right)\end{array}\]
where the short rate process ${r_{t}}$ satisfies (2.1). Both entities are finite due to (2.5).
(Here and in what follows, we omit all proofs which are straightforward.) Taking $u=0$ in Proposition 2.2, we find for all $t\in \left[0,T\right]$
\[\begin{array}{l}\displaystyle {\mathbb{E}_{\mathbb{Q}}}\left[{r_{t}}\right]=\mu \left(t\right)+{\sum \limits_{k=1}^{n}}\left({x_{k}}{e^{-{\lambda _{k}}t}}+{\sigma _{k}}\frac{1-{e^{-{\lambda _{k}}t}}}{{\lambda _{k}}}{\int _{{D_{k}}}}zd{\nu _{k}}\left(z\right)\right),\\ {} \displaystyle {\mathbb{Var}_{\mathbb{Q}}}\left[{r_{t}}\right]={\sum \limits_{k=1}^{n}}{\sigma _{k}^{2}}\frac{1-{e^{-2{\lambda _{k}}t}}}{2{\lambda _{k}}}{\int _{{D_{k}}}}{z^{2}}d{\nu _{k}}\left(z\right).\end{array}\]
Note that it is possible to identify the entities ${\mathbb{E}_{\mathbb{Q}}}\left[{X_{t}^{k}}\right]$ and ${\mathbb{Var}_{\mathbb{Q}}}\left[{X_{t}^{k}}\right]$ inside the latter equations due to (2.1). Moreover, suppose that $\mu \left(t\right)\to \tilde{\mu }$ for $t\to \infty $ where $\tilde{\mu }\in \mathbb{R}$ is a finite constant. Then we observe
\[\begin{aligned}{}& \underset{t\to \infty }{\lim }{\mathbb{E}_{\mathbb{Q}}}\left[{r_{t}}\right]=\tilde{\mu }+{\sum \limits_{k=1}^{n}}\frac{{\sigma _{k}}}{{\lambda _{k}}}{\int _{{D_{k}}}}zd{\nu _{k}}\left(z\right),\\ {} & \underset{t\to \infty }{\lim }{\mathbb{Var}_{\mathbb{Q}}}\left[{r_{t}}\right]={\sum \limits_{k=1}^{n}}\frac{{\sigma _{k}^{2}}}{2{\lambda _{k}}}{\int _{{D_{k}}}}{z^{2}}d{\nu _{k}}\left(z\right)\end{aligned}\]
which both constitute finite constants. This limit behavior entirely stands in line with the requirements imposed on short rate models claimed on p. 46 in [7]. In the next step, we investigate the characteristic function of ${r_{t}}$ which is defined via
where $u\in \mathbb{R}$ and $t\in \left[0,T\right]$.Proposition 2.3.
For $k\in \left\{1,\dots ,n\right\}$ we define the deterministic functions
\[\begin{aligned}{}& {\Lambda _{k}}\left(s,z\right):={\sigma _{k}}{e^{-{\lambda _{k}}\left(t-s\right)}}z,\hspace{2em}{\psi _{k}}\left(t,u\right):=iu{e^{-{\lambda _{k}}t}},\\ {} & {\rho _{k}}\left(t,u\right):={\int _{0}^{t}}{\int _{{D_{k}}}}\big[{e^{iu{\Lambda _{k}}\left(s,z\right)}}-1\big]d{\nu _{k}}\left(z\right)ds.\end{aligned}\]
Then for any $u\in \mathbb{R}$ and $t\in \left[0,T\right]$ the characteristic function of ${r_{t}}$ can be decomposed as
\[ {\Phi _{{r_{t}}}}\left(u\right)={e^{iu\mu \left(t\right)}}{\prod \limits_{k=1}^{n}}{\Phi _{{X_{t}^{k}}}}\left(u\right)\]
where the characteristic function of ${X_{t}^{k}}$ is given by
with deterministic and affine characteristic exponent.
An immediate consequence of Proposition 2.3 is the subsequent affine representation
where we introduced the deterministic functions
We emphasize that ${r_{t}}$ is an affine function of the factors ${X_{t}^{1}},\dots ,{X_{t}^{n}}$ such that our model turns out to be a special case of the affine short rate models considered in Section 3.3 in [14]. To read more on affine processes we refer to [7, 14, 16, 26] and [31]. We next define the moment generating function of ${r_{t}}$ via
which implies the well-known equalities ${\Phi _{{r_{t}}}}\left(u\right)={\kappa _{{r_{t}}}}\left(iu\right)$ and ${\kappa _{{r_{t}}}}\left(v\right)={\Phi _{{r_{t}}}}\left(-iv\right)$. Note that the moment generating function ${\kappa _{{r_{t}}}}\left(v\right)$ is well-defined due to (2.5). In the sequel, we derive the time dynamics of the short rate process.
(2.8)
\[ {\Phi _{{r_{t}}}}\left(u\right)={\prod \limits_{k=1}^{n}}{e^{{\psi _{k}}\left(t,u\right){x_{k}}+{\phi _{k}}\left(t,u\right)}}\]Remark 2.5.
We recall that our model constitutes an extension of the short rate model proposed in [6], whereas we work with multiple pure-jump processes ${L_{t}^{1}},\dots ,{L_{t}^{n}}$ as driving noises instead of the single Brownian motion ${W_{t}}$ appearing in Eq. (1) in [6]. Moreover, comparing Eq. (3) in [6] with Eq. (2.1) above, we see that ${x_{t}}$ and $\varphi \left(t;\alpha \right)$ in [6] correspond in our setup to ${\textstyle\sum _{k=1}^{n}}{X_{t}^{k}}$ and $\mu \left(t\right)$, respectively.
3 Bond prices and forward rates
In this section, we derive representations for zero-coupon bond prices, forward rates and the interest rate curve related to the short rate model introduced in Section 2. To begin with, we introduce a bank account with stochastic interest rate ${r_{t}}$ satisfying
with normalized initial capital ${\beta _{0}}=1$. The solution of (3.1) reads as
where $t\in \left[0,T\right]$. In this setup, the (zero-coupon) bond price at time $t\le T$ with maturity T is given by
where $t\in \left[0,T\right]$ (cf. [6, 7, 26]). Note that $P\left(t,T\right)>0\hspace{2.5pt}\mathbb{Q}$-a.s. $\forall \hspace{2.5pt}t\in \left[0,T\right]$ by construction. Since ${r_{t}}\ge \mu \left(t\right)\hspace{2.5pt}\mathbb{Q}$-a.s. $\forall \hspace{2.5pt}t\in \left[0,T\right]$ [recall Remark 2.1 (a)], we observe
$\mathbb{Q}$-a.s. $\forall \hspace{2.5pt}t\in \left[0,T\right]$ due to (3.3) and the monotonicity of conditional expectations. The upper bound ${M_{t,T}}$ appearing in (3.4) is deterministic and strictly positive for all $0\le t\le T$. If $\mu \left(t\right)\ge 0$, then it holds $P\left(t,T\right)\le 1\hspace{2.5pt}\mathbb{Q}$-a.s. $\forall \hspace{2.5pt}t\in \left[0,T\right]$ (similar to, e.g., the CIR model [10]; also see [7, 21]). On the other hand, if $\mu \left(t\right)<0$, then we only know that ${M_{t,T}}>1$.
(3.3)
\[ P\left(t,T\right):={\beta _{t}}{\mathbb{E}_{\mathbb{Q}}}\left({\beta _{T}^{-1}}|{\mathcal{F}_{t}}\right)={\mathbb{E}_{\mathbb{Q}}}\left(\exp \left\{-{\int _{t}^{T}}{r_{s}}ds\right\}\bigg|{\mathcal{F}_{t}}\right)\]Proposition 3.1.
For $k\in \left\{1,\dots ,n\right\}$ and $t\in \left[0,T\right]$ we define the deterministic functions
Then the bond price at time $t\le T$ with maturity T possesses the affine representation
where the factors ${X_{t}^{k}}$ satisfy (2.7).
(3.5)
\[ \begin{aligned}{}& {A_{k}}\left(t,T\right):={\int _{t}^{T}}\left(-\frac{\mu \left(s\right)}{n}+{\int _{{D_{k}}}}\big[{e^{{\sigma _{k}}{B_{k}}\left(s,T\right)z}}-1\big]d{\nu _{k}}\left(z\right)\right)ds,\\ {} & {B_{k}}\left(t,T\right):=\frac{{e^{-{\lambda _{k}}\left(T-t\right)}}-1}{{\lambda _{k}}}\le 0.\end{aligned}\](3.6)
\[ P\left(t,T\right)={\prod \limits_{k=1}^{n}}{e^{{A_{k}}\left(t,T\right)+{B_{k}}\left(t,T\right){X_{t}^{k}}}}\]Proof.
First of all, we put (2.6) into (2.1) and obtain
\[ {r_{s}}=\mu \left(s\right)+{\sum \limits_{k=1}^{n}}{X_{t}^{k}}{e^{-{\lambda _{k}}\left(s-t\right)}}+{\sum \limits_{k=1}^{n}}{\int _{t}^{s}}{\int _{{D_{k}}}}{\sigma _{k}}{e^{-{\lambda _{k}}\left(s-u\right)}}zd{N_{k}}\left(u,z\right)\]
where $0\le t\le s\le T$. We next substitute the latter equation into (3.3), hereafter apply Fubini’s theorem and identify the functions ${B_{k}}\left(\boldsymbol{\cdot },T\right)$. This procedure yields
\[\begin{aligned}{}P\left(t,T\right)& =\exp \left\{-{\int _{t}^{T}}\mu \left(s\right)ds+{\sum \limits_{k=1}^{n}}{B_{k}}\left(t,T\right){X_{t}^{k}}\right\}\\ {} & \hspace{1em}\times {\mathbb{E}_{\mathbb{Q}}}\left(\exp \left\{{\sum \limits_{k=1}^{n}}{\int _{t}^{T}}{\int _{{D_{k}}}}{\sigma _{k}}{B_{k}}\left(s,T\right)zd{N_{k}}\left(s,z\right)\right\}\bigg|{\mathcal{F}_{t}}\right).\end{aligned}\]
Taking the independent increment property of the ($\mathbb{Q}$-independent) Lévy processes ${L^{1}},\dots ,{L^{n}}$ into account, we obtain
\[\begin{aligned}{}& {\mathbb{E}_{\mathbb{Q}}}\left(\exp \left\{{\sum \limits_{k=1}^{n}}{\int _{t}^{T}}{\int _{{D_{k}}}}{\sigma _{k}}{B_{k}}\left(s,T\right)zd{N_{k}}\left(s,z\right)\right\}\bigg|{\mathcal{F}_{t}}\right)\\ {} & \hspace{1em}={\prod \limits_{k=1}^{n}}{\mathbb{E}_{\mathbb{Q}}}\left[\exp \left\{{\int _{t}^{T}}{\int _{{D_{k}}}}{\sigma _{k}}{B_{k}}\left(s,T\right)zd{N_{k}}\left(s,z\right)\right\}\right]\end{aligned}\]
where $t\in \left[0,T\right]$. The usual expectations appearing here can be handled by the Lévy–Khinchin formula for additive processes (see, e.g., [9, 30, 36]) which leads us to
\[\begin{aligned}{}& {\mathbb{E}_{\mathbb{Q}}}\left[\exp \left\{{\int _{t}^{T}}{\int _{{D_{k}}}}{\sigma _{k}}{B_{k}}\left(s,T\right)zd{N_{k}}\left(s,z\right)\right\}\right]\\ {} & \hspace{1em}=\exp \left\{{\int _{t}^{T}}{\int _{{D_{k}}}}\big[{e^{{\sigma _{k}}{B_{k}}\left(s,T\right)z}}-1\big]d{\nu _{k}}\left(z\right)ds\right\}.\end{aligned}\]
Putting the latter equations together and identifying the functions ${A_{k}}\left(\boldsymbol{\cdot },T\right)$, we end up with the asserted representation (3.6). □Recall that the bond price in (3.6) is the product of exponential affine functions of the factors ${X_{t}^{1}},\dots ,{X_{t}^{n}}$ (but not of ${r_{t}}$). Also note that for all $k\in \left\{1,\dots ,n\right\}$ and $t\in \left[0,T\right]$ it holds
We remark that the functions ${B_{k}}\left(t,T\right)$ in (3.5) possess the same structure as the corresponding ones in the Vasicek model (cf. [38], or [7, 16, 21]). For all $t\in \left[0,T\right]$ Eq. (3.6) can be rewritten as
which implies $P\left(T,T\right)=1$ due to (3.7). Moreover, from (3.5) we infer the time derivatives
where ${A^{\prime }_{k}}:={\partial _{t}}{A_{k}}$ and ${B^{\prime }_{k}}:={\partial _{t}}{B_{k}}$. Hence, the functions ${B_{k}}\left(t,T\right)\le 0$ are strictly monotone increasing in t. Also note that the formulas found in (3.9) entirely stand in line with those claimed in (4.4)–(4.5) in [31]. From (3.5), (3.7) and (3.9) we deduce the following system of ordinary differential equations (ODEs)
(3.8)
\[ P\left(t,T\right)=\exp \left\{{\sum \limits_{k=1}^{n}}\left[{A_{k}}\left(t,T\right)+{B_{k}}\left(t,T\right){X_{t}^{k}}\right]\right\}\](3.9)
\[ {A^{\prime }_{k}}\left(t,T\right)=\frac{\mu \left(t\right)}{n}-{\int _{{D_{k}}}}\big[{e^{{\sigma _{k}}{B_{k}}\left(t,T\right)z}}-1\big]d{\nu _{k}}\left(z\right),\hspace{2em}{B^{\prime }_{k}}\left(t,T\right)={e^{-{\lambda _{k}}\left(T-t\right)}}>0\]
\[\begin{aligned}{}& {A_{k}}\left(t,T\right)=-{\int _{t}^{T}}{A^{\prime }_{k}}\left(s,T\right)ds,\hspace{2em}{B^{\prime }_{k}}\left(t,T\right)=1+{\lambda _{k}}{B_{k}}\left(t,T\right),\\ {} & {A_{k}}\left(T,T\right)={B_{k}}\left(T,T\right)=0\end{aligned}\]
where $t\in \left[0,T\right]$ and $k\in \left\{1,\dots ,n\right\}$. We are now prepared to derive the time dynamics of the bond price process ${\left(P\left(t,T\right)\right)_{t\in \left[0,T\right]}}$.Proposition 3.2.
For $k\in \left\{1,\dots ,n\right\}$, $t\in \left[0,T\right]$ and $z\in {D_{k}}$ we define the deterministic functions
with ${B_{k}}\left(t,T\right)$ as in (3.5). Then the bond price satisfies the t-dynamics under $\mathbb{Q}$
Recall that it holds ${\zeta _{k}}\left(t,T,z\right)\le 0$, since ${\sigma _{k}}\hspace{2.5pt}{B_{k}}\left(t,T\right)\hspace{2.5pt}z\le 0$ for all k, t and z. We stress that (3.11) possesses the same structure as the corresponding Eq. (10.9) in [37], whereas the latter stems from a Brownian motion model without jumps. In the next step, we provide the solution of the SDE (3.11).
Proposition 3.3.
For all $t\in \left[0,T\right]$ the solution of (3.11) under $\mathbb{Q}$ reads as
where the initial value $P\left(0,T\right)$ is deterministic and fulfills $P\left(0,T\right)>0$.
(3.12)
\[\begin{aligned}{}P\left(t,T\right)& =P\left(0,T\right)\exp \Bigg\{{\int _{0}^{t}}{r_{s}}ds-{\sum \limits_{k=1}^{n}}{\int _{0}^{t}}{\int _{{D_{k}}}}{\zeta _{k}}\left(s,T,z\right)d{\nu _{k}}\left(z\right)ds\\ {} & \hspace{1em}+{\sum \limits_{k=1}^{n}}{\int _{0}^{t}}{\int _{{D_{k}}}}{\sigma _{k}}{B_{k}}\left(s,T\right)zd{N_{k}}\left(s,z\right)\Bigg\}\end{aligned}\]Furthermore, for all $t\in \left[0,T\right]$ let us introduce the discontinuous Doléans-Dade exponential
where ${h_{k}}\left(s,z\right)$ is an arbitrary integrable deterministic function (which may also depend on T). We recall that ${\Xi _{0}^{k}}=1$ and that ${\left({\Xi _{t}^{k}}\right)_{t\in \left[0,T\right]}}$ constitutes a local $\mathbb{Q}$-martingale. With definition (3.13) at hand, we can express Eq. (3.12) as follows.
(3.13)
\[\begin{aligned}{}{\Xi _{t}^{k}}& :=\mathcal{E}{\big({h_{k}}\ast {\tilde{N}_{k}^{\mathbb{Q}}}\big)_{t}}:=\exp \Bigg\{{\int _{0}^{t}}{\int _{{D_{k}}}}{h_{k}}\left(s,z\right)d{\tilde{N}_{k}^{\mathbb{Q}}}\left(s,z\right)\\ {} & \hspace{1em}-{\int _{0}^{t}}{\int _{{D_{k}}}}\big[{e^{{h_{k}}\left(s,z\right)}}-1-{h_{k}}\left(s,z\right)\big]d{\nu _{k}}\left(z\right)ds\Bigg\}\end{aligned}\]Corollary 3.4.
For all $0\le t\le T$ the bond price satisfies
where ${\beta _{t}}$ is the bank account process given in (3.2), $\mathcal{E}$ denotes the Doleáns-Dade exponential defined in (3.13) and ${\xi _{k}}\left(s,z\right):={\sigma _{k}}\hspace{2.5pt}{B_{k}}\left(s,T\right)\hspace{2.5pt}z=\log \left(1+{\zeta _{k}}\left(s,T,z\right)\right)$ is a deterministic function.
(3.14)
\[ P\left(t,T\right)=P\left(0,T\right){\beta _{t}}{\prod \limits_{k=1}^{n}}\mathcal{E}{\big({\xi _{k}}\ast {\tilde{N}_{k}^{\mathbb{Q}}}\big)_{t}}\]Moreover, for all $0\le t\le T$ we define the discounted bond price
where $\hat{P}\left(0,T\right)=P\left(0,T\right)$. From (3.3) we deduce $\hat{P}\left(t,T\right)={\mathbb{E}_{\mathbb{Q}}}({\beta _{T}^{-1}}|\hspace{2.5pt}{\mathcal{F}_{t}})$ such that $\hat{P}\left(t,T\right)$ constitutes an ${\mathcal{F}_{t}}$-adapted (true) martingale under $\mathbb{Q}$, as required by the risk-neutral pricing theory. Plugging (3.14) into (3.15), for all $t\in \left[0,T\right]$ we obtain
where $P\left(0,T\right)$ is deterministic and ${\xi _{k}}$ is such as defined in Corollary 3.4. We obtain the following result.
(3.16)
\[ \hat{P}\left(t,T\right)=P\left(0,T\right){\prod \limits_{k=1}^{n}}\mathcal{E}{\big({\xi _{k}}\ast {\tilde{N}_{k}^{\mathbb{Q}}}\big)_{t}}\]Proposition 3.5.
For all $t\in \left[0,T\right]$ the discounted bond price satisfies the $\mathbb{Q}$-martingale dynamics
\[ \frac{d\hat{P}\left(t,T\right)}{\hat{P}\left(t-,T\right)}={\sum \limits_{k=1}^{n}}{\int _{{D_{k}}}}{\zeta _{k}}\left(t,T,z\right)d{\tilde{N}_{k}^{\mathbb{Q}}}\left(t,z\right)\]
where the deterministic functions ${\zeta _{k}}\left(t,T,z\right)$ are such as defined in (3.10).
With reference to [7], we define the instantaneous forward rate at time t with maturity T via
where $t\in \left[0,T\right]$ and ${\partial _{T}}$ denotes the differential operator with respect to T. Equation (3.17) is equivalent to the representation
Lemma 3.6.
For all $k\in \left\{1,\dots ,n\right\}$ and $t\in \left[0,T\right]$ it holds
\[\begin{aligned}{}{\partial _{T}}{A_{k}}\left(t,T\right)& =-\frac{\mu \left(T\right)}{n}-{\sigma _{k}}{\int _{t}^{T}}{\int _{{D_{k}}}}z{e^{{\sigma _{k}}{B_{k}}\left(s,T\right)z-{\lambda _{k}}\left(T-s\right)}}d{\nu _{k}}\left(z\right)ds,\\ {} \hspace{2.5pt}{\partial _{T}}{B_{k}}\left(t,T\right)& =-{e^{-{\lambda _{k}}\left(T-t\right)}}.\end{aligned}\]
Proof.
By the definition of ${B_{k}}\left(t,T\right)$ claimed in (3.5) we find
so that the functions ${B_{k}}\left(t,T\right)$ are strictly monotone decreasing in T. From (3.5) and (3.10) we further deduce
\[ {\partial _{T}}{A_{k}}\left(t,T\right)=-\frac{\mu \left(T\right)}{n}+{\partial _{T}}\left({\int _{t}^{T}}{\int _{{D_{k}}}}{\zeta _{k}}\left(s,T,z\right)d{\nu _{k}}\left(z\right)ds\right)\]
whereas Fubini’s theorem (see, e.g., Theorem 2.2 in [3]) leads us to
\[ {\partial _{T}}\left({\int _{t}^{T}}{\int _{{D_{k}}}}{\zeta _{k}}\left(s,T,z\right)d{\nu _{k}}\left(z\right)ds\right)={\int _{{D_{k}}}}{\partial _{T}}\left({\int _{t}^{T}}{\zeta _{k}}\left(s,T,z\right)ds\right)d{\nu _{k}}\left(z\right).\]
(We are able to apply Fubini’s theorem here, since the deterministic function ${\zeta _{k}}(s,\hspace{-0.1667em}T,\hspace{-0.1667em}z)$ is measurable and square-integrable with respect to $s\in \left[0,T\right]$ and $z\in {D_{k}}$.) The Leibniz formula for parameter integrals (see Lemma 2.4.1 on p. 13 in [28]) yields
\[\begin{aligned}{}{\partial _{T}}\left({\int _{t}^{T}}{\zeta _{k}}\left(s,T,z\right)ds\right)& ={\zeta _{k}}\left(T,T,z\right)+{\int _{t}^{T}}{\partial _{T}}{\zeta _{k}}\left(s,T,z\right)ds\\ {} & =-{\sigma _{k}}{\int _{t}^{T}}z{e^{{\sigma _{k}}{B_{k}}\left(s,T\right)z-{\lambda _{k}}\left(T-s\right)}}ds\end{aligned}\]
where we used (3.10), (3.7) and (3.19). Putting these formulas together, the proof is complete. □Proposition 3.7.
For all $t\in \left[0,T\right]$ the instantaneous forward rate can be represented as
where the factor processes ${X_{t}^{k}}$ satisfy (2.7) and ${B_{k}}\left(s,T\right)$ is like defined in (3.5).
(3.20)
\[ f\left(t,T\right)=\mu \left(T\right)+{\sum \limits_{k=1}^{n}}{\int _{t}^{T}}{\int _{{D_{k}}}}{\sigma _{k}}z{e^{{\sigma _{k}}{B_{k}}\left(s,T\right)z-{\lambda _{k}}\left(T-s\right)}}d{\nu _{k}}\left(z\right)ds+{\sum \limits_{k=1}^{n}}{X_{t}^{k}}{e^{-{\lambda _{k}}\left(T-t\right)}}\]Replacing T by t in (3.20), we immediately find $f\left(t,t\right)={r_{t}}$ due to (2.1). This equality stands in line with the usual conventions in interest rate theory (see, e.g., [7, 16, 21]).
Proposition 3.8.
For all $t\in \left[0,T\right]$ the instantaneous forward rate fulfills the pure-jump multi-factor HJM type equation
where the deterministic initial value is given by $f\left(0,T\right)=-{\partial _{T}}\log P\left(0,T\right)$.
(3.21)
\[ f\left(t,T\right)=f\left(0,T\right)+{\sum \limits_{k=1}^{n}}{\int _{0}^{t}}{\int _{{D_{k}}}}{\sigma _{k}}z{e^{-{\lambda _{k}}\left(T-s\right)}}\big\{d{N_{k}}\left(s,z\right)-{e^{{\sigma _{k}}{B_{k}}\left(s,T\right)z}}d{\nu _{k}}\left(z\right)ds\big\}\]In what follows, we illustrate how our forward rate model can be fitted to the initially observed term structure. This procedure is often called market-consistent calibration in the literature. For this purpose, we denote by ${f^{M}}\left(0,T\right)$ the deterministic initial forward rate. If $f\left(0,T\right)={f^{M}}\left(0,T\right)$ and hence, if $P\left(0,T\right)={P^{M}}\left(0,T\right)$ for all maturity times $T>0$, then the underlying model is called market-consistent.
Proposition 3.9.
The forward rate model (3.20)–(3.21) can be market-consistently calibrated to a given term structure ${f^{M}}\left(0,T\right)$ by choosing the floor function $\mu \left(\boldsymbol{\cdot }\right)\hspace{2.5pt}in$ (3.20) according to
for all maturity times $T>0$.
(3.22)
\[ \mu \left(T\right)={f^{M}}\left(0,T\right)-{\sum \limits_{k=1}^{n}}\left({x_{k}}{e^{-{\lambda _{k}}T}}+{\int _{0}^{T}}{\int _{{D_{k}}}}{\sigma _{k}}z{e^{{\sigma _{k}}{B_{k}}\left(s,T\right)z-{\lambda _{k}}\left(T-s\right)}}d{\nu _{k}}\left(z\right)ds\right)\]Note that the floor function $\mu \left(t\right)$ for all $t\in \left[0,T\right]$ can be obtained from (3.22) by replacing T by t therein. Moreover, we define the interest rate curve at time $t<T$ with maturity T via
This object is called continuously-compounded spot rate on p. 60 in [7]. It obviously holds
where $t\in \left[0,T\right[$. Comparing the exponent in (3.24) with that in (3.8), we infer
\[ R\left(t,T\right)=\frac{1}{t-T}{\sum \limits_{k=1}^{n}}\big[{A_{k}}\left(t,T\right)+{B_{k}}\left(t,T\right){X_{t}^{k}}\big]\]
where ${A_{k}}$ and ${B_{k}}$ are such as defined in (3.5). Hence, it turns out that the interest rate curve $R\left(t,T\right)$ can be represented as a sum of affine functions of the pure-jump OU factors ${X_{t}^{1}},\dots ,{X_{t}^{n}}$. In this sense, our short rate model possesses an affine term structure (cf. Section 3.2.4 in [7], or [14, 16, 21]). The latter observation entirely stands in line with (3.8).Proposition 3.10.
For all $t\in \left[0,T\right[$ the interest rate curve possesses the representation
\[\begin{aligned}{}R\left(t,T\right)& =\frac{1}{t-T}\Bigg(\log P\left(0,T\right)+{\int _{0}^{t}}{r_{s}}ds-{\sum \limits_{k=1}^{n}}{\int _{0}^{t}}{\int _{{D_{k}}}}{\zeta _{k}}\left(s,T,z\right)d{\nu _{k}}\left(z\right)ds\\ {} & \hspace{1em}+{\sum \limits_{k=1}^{n}}{\int _{0}^{t}}{\int _{{D_{k}}}}{\sigma _{k}}{B_{k}}\left(s,T\right)zd{N_{k}}\left(s,z\right)\Bigg)\end{aligned}\]
where ${\zeta _{k}}$ and ${B_{k}}$ are such as defined in (3.10) and (3.5), respectively.
Proposition 3.11.
For all $t\in \left[0,T\right]$ the integrated short rate process can be represented as
where the deterministic functions ${B_{k}}$ are such as defined in (3.5).
(3.25)
\[ {I_{t}}:={\int _{0}^{t}}{r_{s}}ds={\int _{0}^{t}}\mu \left(s\right)ds-{\sum \limits_{k=1}^{n}}{x_{k}}{B_{k}}\left(0,t\right)-{\sum \limits_{k=1}^{n}}{\int _{0}^{t}}{\int _{{D_{k}}}}{\sigma _{k}}{B_{k}}\left(s,t\right)zd{N_{k}}\left(s,z\right)\]Proof.
We substitute (2.1) and (2.7) into the definition of ${I_{t}}$ and obtain
\[ {I_{t}}={\int _{0}^{t}}\mu \left(s\right)ds-{\sum \limits_{k=1}^{n}}{x_{k}}{B_{k}}\left(0,t\right)+{\sum \limits_{k=1}^{n}}{\int _{0}^{t}}{\int _{0}^{s}}{\int _{{D_{k}}}}{\sigma _{k}}{e^{-{\lambda _{k}}\left(s-u\right)}}zd{N_{k}}\left(u,z\right)ds\]
where ${B_{k}}$ is like defined in (3.5). An application of Fubini’s theorem (see Theorem 2.2 in [3]) yields
so that the proof is complete. □4 Option pricing
In this section, we investigate the evaluation of a plain vanilla option written on the zero-coupon bond price $P\left(\boldsymbol{\cdot },T\right)$. With reference to the risk-neutral pricing theory, the price at time $t\le \tau $ of an option with payoff ${H_{\tau }}$ at the maturity time τ reads as
where β is the bank account process given in (3.2) and $\mathbb{Q}$ denotes a risk-neutral pricing measure (cf. Eq. (3.1) in [7]). We now consider a call option written on the bond price $P\left(\boldsymbol{\cdot },T\right)$ with maturity time T satisfying $T\ge \tau $. The payoff of the call option written on $P\left(\tau ,T\right)$ with deterministic strike price $K>0$ and maturity time τ is then given by
In what follows, we define the Fourier transform, respectively inverse Fourier transform, of a real-valued deterministic function $p\left(\boldsymbol{\cdot }\right)\in {\mathcal{L}^{1}}\left(\mathbb{R}\right)$ via
(4.1)
\[ {C_{t}}={\beta _{t}}{\mathbb{E}_{\mathbb{Q}}}\big({\beta _{\tau }^{-1}}{H_{\tau }}|{\mathcal{F}_{t}}\big)={\mathbb{E}_{\mathbb{Q}}}\big({e^{-{\textstyle\textstyle\int _{t}^{\tau }}{r_{s}}ds}}{H_{\tau }}|{\mathcal{F}_{t}}\big)\](4.2)
\[ {H_{\tau }}={\left[P\left(\tau ,T\right)-K\right]^{+}}:=\max \left\{0,P\left(\tau ,T\right)-K\right\}.\]
\[ \hat{p}\left(y\right):=\frac{1}{2\pi }{\int _{\mathbb{R}}}p\left(u\right){e^{-iyu}}du,p\left(u\right)={\int _{\mathbb{R}}}\hat{p}\left(y\right){e^{iyu}}dy.\]
Proposition 4.1.
[call option on bond price] For all $0\le t\le \tau \le T$ the price of a call option with payoff ${H_{\tau }}$ given in (4.2), strike price $K>0$ and maturity time τ can be expressed as
where the integrated short rate process ${I_{t}}$ is such as defined in (3.25) and
constitute deterministic functions, while $a>1$ is an arbitrary real-valued constant. Herein, the functions ${\zeta _{k}}$ and ${B_{k}}$ are such as defined in (3.10) and (3.5), respectively, while $P\left(0,T\right)$ denotes the deterministic initial bond price.
(4.3)
\[\begin{aligned}{}{C_{t}}& ={\int _{\mathbb{R}}}\hat{q}\left(y\right)\exp \Bigg\{{I_{t}}+\theta \left(\tau ,y\right)+{\sum \limits_{k=1}^{n}}{\overline{\psi }_{k}}\left(t,\tau ,y\right)\\ {} & \hspace{1em}+{\sum \limits_{k=1}^{n}}{\int _{0}^{t}}{\int _{{D_{k}}}}{\eta _{k}}\left(s,z,y\right)d{N_{k}}\left(s,z\right)\Bigg\}dy\end{aligned}\](4.4)
\[\begin{aligned}{}{\eta _{k}}\left(s,z,y\right)& :={\eta _{k}}\left(s,z,T,\tau ,y\right)\\ {} & :={\sigma _{k}}\left[\left(a+iy\right){B_{k}}\left(s,T\right)-\left(a+iy-1\right){B_{k}}\left(s,\tau \right)\right]z,\\ {} \hat{q}\left(y\right)& :=\frac{P{\left(0,T\right)^{a+iy}}}{2\pi \left(a+iy\right)\left(a+iy-1\right){K^{a+iy-1}}},\\ {} {\overline{\psi }_{k}}\left(t,\tau ,y\right)& :={\int _{t}^{\tau }}{\int _{{D_{k}}}}\big[{e^{{\eta _{k}}\left(s,z,y\right)}}-1\big]d{\nu _{k}}\left(z\right)ds,\\ {} \theta \left(\tau ,y\right)& :=\left(a+iy-1\right)\left({\int _{0}^{\tau }}\mu \left(s\right)ds-{\sum \limits_{k=1}^{n}}{x_{k}}{B_{k}}\left(0,\tau \right)\right)\\ {} & \hspace{1em}-\left(a+iy\right){\sum \limits_{k=1}^{n}}{\int _{0}^{\tau }}{\int _{{D_{k}}}}{\zeta _{k}}\left(s,T,z\right)d{\nu _{k}}\left(z\right)ds\end{aligned}\]Proof.
We substitute (4.2) and (3.12) into (4.1) and obtain
\[ {C_{t}}={\mathbb{E}_{\mathbb{Q}}}\left({e^{{I_{t}}-{I_{\tau }}}}{\big[P\left(0,T\right){e^{{G_{\tau }}}}-K\big]^{+}}|{\mathcal{F}_{t}}\right)\]
where ${I_{t}}$ denotes the integrated short rate process defined in (3.25) and
\[ {G_{\tau }}:={I_{\tau }}-{\sum \limits_{k=1}^{n}}{\int _{0}^{\tau }}{\int _{{D_{k}}}}{\zeta _{k}}\left(s,T,z\right)d{\nu _{k}}\left(z\right)ds+{\sum \limits_{k=1}^{n}}{\int _{0}^{\tau }}{\int _{{D_{k}}}}{\sigma _{k}}{B_{k}}\left(s,T\right)zd{N_{k}}\left(s,z\right)\]
is a real-valued stochastic process. For $u\in \mathbb{R}$ we introduce the deterministic function
where $a>1$ is a constant real-valued dampening parameter ensuring the integrability of the payoff function. Indeed, it holds $q\left(\boldsymbol{\cdot }\right)\in {\mathcal{L}^{1}}\left(\mathbb{R}\right)$. With the latter definition at hand, we obtain
\[ {C_{t}}={\mathbb{E}_{\mathbb{Q}}}\big({e^{{I_{t}}-{I_{\tau }}+a{G_{\tau }}}}q\left({G_{\tau }}\right)|{\mathcal{F}_{t}}\big).\]
With reference to [8] (also see [18]), we apply the inverse Fourier transform on the latter equation and hereafter, use Fubini’s theorem which leads us to
\[ {C_{t}}={\int _{\mathbb{R}}}\hat{q}(y){\mathbb{E}_{\mathbb{Q}}}\big({e^{{Z_{t,\tau }}}}|{\mathcal{F}_{t}}\big)dy\]
where we have set
for all $0\le t\le \tau $. By merging the definition of ${G_{\tau }}$ and (3.25) into the definition of ${Z_{t,\tau }}$ we deduce
\[ {Z_{t,\tau }}={I_{t}}+\theta \left(\tau ,y\right)+{\sum \limits_{k=1}^{n}}{\int _{0}^{\tau }}{\int _{{D_{k}}}}{\eta _{k}}\left(s,z,y\right)d{N_{k}}\left(s,z\right)\]
where we identified the deterministic functions $\theta \left(\tau ,y\right)$ and ${\eta _{k}}\left(s,z,y\right)$ defined in (4.4). Hence,
\[\begin{aligned}{}{\mathbb{E}_{\mathbb{Q}}}\big({e^{{Z_{t,\tau }}}}|{\mathcal{F}_{t}}\big)& =\exp \left\{{I_{t}}+\theta \left(\tau ,y\right)+{\sum \limits_{k=1}^{n}}{\int _{0}^{t}}{\int _{{D_{k}}}}{\eta _{k}}\left(s,z,y\right)d{N_{k}}\left(s,z\right)\right\}\\ {} & \hspace{1em}\times {\mathbb{E}_{\mathbb{Q}}}\left[\exp \left\{{\sum \limits_{k=1}^{n}}{\int _{t}^{\tau }}{\int _{{D_{k}}}}{\eta _{k}}\left(s,z,y\right)d{N_{k}}\left(s,z\right)\right\}\right]\end{aligned}\]
since ${I_{t}}$ is ${\mathcal{F}_{t}}$-adapted and $\theta \left(\tau ,y\right)$ is deterministic. In the derivation of the latter equation, we used the independent increment property under $\mathbb{Q}$ of the involved pure-jump integrals. We next apply the Lévy–Khinchin formula for additive processes (see, e.g., [9, 30, 36]) and derive
\[ {\mathbb{E}_{\mathbb{Q}}}\left[\exp \left\{{\sum \limits_{k=1}^{n}}{\int _{t}^{\tau }}{\int _{{D_{k}}}}{\eta _{k}}\left(s,z,y\right)d{N_{k}}\left(s,z\right)\right\}\right]=\exp \left\{{\sum \limits_{k=1}^{n}}{\overline{\psi }_{k}}\left(t,\tau ,y\right)\right\}\]
where the characteristic exponents ${\overline{\psi }_{k}}\left(t,\tau ,y\right)$ are such as defined in (4.4). Putting the latter equations together, we eventually end up with (4.3). The expression for the Fourier transform $\hat{q}\left(y\right)$ is obtained by straightforward calculations using the definition of the function $q\left(u\right)$. □Corollary 4.2.
In the special case $t=0$, the call option price formula (4.3) simplifies to
which is deterministic.
5 Practical applications
In this section, we show how the short rate model introduced in Section 2 can be implemented in practical applications. For this purpose, we now present more detailed expressions in order to prepare our model for a possible calibration of the involved parameters. First of all, let us recall that the increasing compound Poisson processes ${L_{t}^{k}}$ defined in (2.3) for every $k\in \left\{1,\dots ,n\right\}$ and $t\in \left[0,T\right]$ can be expressed as
(cf. Section 5.3.2 in [37]) where ${N_{t}^{k}}$ constitutes a standard Poisson process under $\mathbb{Q}$ with deterministic jump intensity ${\alpha _{k}}>0$. That is, ${N_{t}^{k}}\sim Poi\left({\alpha _{k}}t\right)$ such that for all $m\in {\mathbb{N}_{0}}$ it holds
\[ \mathbb{Q}\big({N_{t}^{k}}=m\big)=\frac{{\left({\alpha _{k}}t\right)^{m}}}{m!}{e^{-{\alpha _{k}}t}}.\]
The strictly positive jump amplitudes of the Lévy process ${L_{t}^{k}}$ are modeled by the i.i.d. random variables ${Y_{1}^{k}},{Y_{2}^{k}},\dots \hspace{0.1667em}$ which take values in the set ${D_{k}}\subseteq \left]0,\infty \right[$. We recall that the random variables ${Y_{1}^{k}},{Y_{2}^{k}},\dots \hspace{0.1667em}$ are independent of the Poisson processes ${N_{t}^{\overline{k}}}$ for all combinations of indices $k,\overline{k}\in \left\{1,\dots ,n\right\}$. We further put ${c_{k}}:={\mathbb{E}_{\mathbb{Q}}}[{Y_{1}^{k}}]\in {D_{k}}$ and recall that the compensated compound Poisson process ${\left({L_{t}^{k}}-{c_{k}}{\alpha _{k}}t\right)_{t\in \left[0,T\right]}}$ constitutes an $\left({\mathcal{F}_{t}},\mathbb{Q}\right)$-martingale for each k which implies
due to (2.3) and (2.4). We stress that the Poisson processes ${N_{t}^{k}}$ shall not be mixed up with the Poisson random measures $d{N_{k}}\left(t,z\right)$.In the following, we propose a number of probability distributions living on the positive half-line (recall Section B.1.2 in [37]) which constitute suitable candidates for the modeling of the jump size distribution in our new short rate model. As a first example, we propose to work with the gamma distribution and thus, assume that each random variable ${Y_{j}^{k}}$ is exponentially distributed under $\mathbb{Q}$ with parameter ${\varepsilon _{k}}>0$ for all j and k. In this case, the related Lévy measure possesses the Lebesgue density
where $z\in {D_{k}}=\left]0,\infty \right[$ and $k\in \left\{1,\dots ,n\right\}$. We find ${c_{k}}=1/{\varepsilon _{k}}$ and ${Y_{j}^{k}}\sim \Gamma \left(1,{\varepsilon _{k}}\right)$. Hence, following the notation used in Section 5.5.1 in [37], we state that we presently are in a $\Gamma \left({\alpha _{k}},{\varepsilon _{k}}\right)$-Ornstein–Uhlenbeck process setup (also see Section 8.2 in [31] and Example 15.1 in [9] in this context).
Proposition 5.1.
Suppose that the random variables ${Y_{j}^{k}}$ in (5.1) are exponentially distributed (i.e. $\Gamma \left(1,{\varepsilon _{k}}\right)$-distributed) under $\mathbb{Q}$ with parameters ${\varepsilon _{k}}>0$ for all j and k. Then, for $u\in \mathbb{R}$ and $t\in \left[0,T\right]$ the characteristic function of ${L_{t}^{k}}$ is given by
\[ {\Phi _{{L_{t}^{k}}}}\left(u\right)=\exp \left\{\frac{iu{\alpha _{k}}t}{{\varepsilon _{k}}-iu}\right\}\]
where ${\alpha _{k}}$ denotes the jump intensity of the standard Poisson process ${N_{t}^{k}}$ appearing in (5.1).
Proof.
An immediate consequence of Proposition 5.1 is the following representation for the moment generating function of ${L_{t}^{k}}$ being valid for all $v\in \mathbb{R}\setminus \left\{{\varepsilon _{k}}\right\}$
Proposition 5.2.
Assume that the random variables ${Y_{j}^{k}}$ in (5.1) are exponentially distributed (i.e. $\Gamma \left(1,{\varepsilon _{k}}\right)$-distributed) under $\mathbb{Q}$ with parameters ${\varepsilon _{k}}>0$ for all j and k. Then, for all $t\in \left[0,T\right]$, $k\in \left\{1,\dots ,n\right\}$ and $x\in \mathbb{R}$ the probability density function of ${L_{t}^{k}}$ under $\mathbb{Q}$ takes the form
Proof.
First, note that it holds
\[ {\Phi _{{L_{t}^{k}}}}\left(-u\right)={\int _{0}^{\infty }}{e^{-iux}}{f_{{L_{t}^{k}}}}\left(x\right)dx=2\pi {\hat{f}_{{L_{t}^{k}}}}\left(u\right)\]
due to the definitions of the characteristic function and the Fourier transform claimed in the sequel of (4.2). We next apply the inverse Fourier transform which yields the density function
\[ {f_{{L_{t}^{k}}}}\left(x\right)=\frac{1}{2\pi }{\int _{0}^{\infty }}{\Phi _{{L_{t}^{k}}}}\left(-u\right){e^{iux}}du.\]
We finally plug the result of Proposition 5.1 into the latter equation which completes the proof. □We stress that Eq. (5.2) can be substituted into the corresponding formulas appearing in the previous Propositions 2.2, 2.3, 3.1, 3.3, 3.7–3.10 and 4.1 yielding more explicit expressions for the involved entities, yet associated with gamma-distributed jump amplitudes in the underlying short rate model. We illustrate this statement by an application of Eq. (5.2) on Proposition 2.3. The precise result reads as follows.
Proposition 5.3.
Suppose that the random variables ${Y_{j}^{k}}\hspace{2.5pt}in$ (5.1) are exponentially distributed (i.e. $\Gamma \left(1,{\varepsilon _{k}}\right)$-distributed) under $\mathbb{Q}$ with parameters ${\varepsilon _{k}}>0$ for all j and k. Let ${\sigma _{k}}>0$ be the constant volatility coefficients introduced in (2.2). Then, for all $t\in \left[0,T\right]$ and $v\in \mathbb{R}$ with $v<{\min _{k}}\left\{{\varepsilon _{k}}/{\sigma _{k}}\right\}$, $k\in \left\{1,\dots ,n\right\}$, the moment generating function under $\mathbb{Q}$ of the short rate process ${r_{t}}$ reads as
with deterministic functions
Proof.
For each $k\in \left\{1,\dots ,n\right\}$ we define the deterministic functions ${b_{k}}\left(s,v\right):=v\hspace{2.5pt}{\sigma _{k}}\hspace{2.5pt}{e^{-{\lambda _{k}}\left(t-s\right)}}-{\varepsilon _{k}}$ which satisfy ${b_{k}}\left(s,v\right)<0$ whenever $s\in \left[0,t\right]$ and $v<{\min _{k}}\left\{{\varepsilon _{k}}/{\sigma _{k}}\right\}$. In this setting, we combine Eq. (5.2) with the definitions of ${\rho _{k}}$ and ${\Lambda _{k}}$ given in Proposition 2.3 and obtain
\[ {\rho _{k}}\left(t,-iv\right)=-{\alpha _{k}}{\int _{0}^{t}}\frac{{\varepsilon _{k}}+{b_{k}}\left(s,v\right)}{{b_{k}}\left(s,v\right)}ds.\]
We perform the integration and obtain the formula for ${\rho _{k}}$ claimed in the proposition. The representation for the moment generating function ${\kappa _{{r_{t}}}}\left(v\right)$ finally follows from Proposition 2.3. □Other distributional choices for the random variables ${Y_{j}^{k}}$ modeling the jump amplitudes would be, for example, the inverse Gaussian distribution (see Section 5.5.2 in [37]), the generalized inverse Gaussian distribution (see Section 5.3.5 in [37]) or the tempered stable distribution (see Section 5.3.6 in [37]). The related formulas for the Lebesgue density of the Lévy measure $d{\nu _{k}}\left(z\right)$ corresponding to Eq. (5.2) can be found in [37].
Remark 5.4.
We recall that the time-homogeneous compound Poisson processes ${L_{t}^{k}}$ introduced in (2.3) can be simulated according to Algorithms 6.1 and 6.2 in [9]. Further, in our model it is easily possible to calculate the moments of ${X_{t}^{k}}$ and ${r_{t}}$ (see the sequel of Proposition 2.2) so that our model can be fitted to any yield curve observed in the market by using the moment estimation method described in Section 7.2.2 in [9]. This procedure is also called moment matching, as the underlying idea is to make the empirical moments match with the theoretical moments of the model by finding a suitable parameter set.
6 A post-crisis model extension
In this section, we propose a post-crisis extension of the pure-jump lower-bounded short rate model introduced in Section 2. (To read more on post-crisis interest rate models, the reader is referred to [11–14, 17, 22–26, 32–34].) Inspired by the modeling setups presented in [33] and Chapter 2 in [26], for all $t\in \left[0,T\right]$ we define the short rate spread under $\mathbb{Q}$ by the stochastic process
showing a similar structure as (2.1). Herein, ${\mu ^{\ast }}\left(t\right)\ge 0$ constitutes an integrable real-valued deterministic function and the factors ${X_{t}^{k}}$ satisfy the SDE (2.2), but presently for indices $k\in \left\{n+1,\dots ,l\right\}$ where $l\in \mathbb{N}$ with $l>n$. Note that it holds ${s_{t}}\ge {\mu ^{\ast }}\left(t\right)\hspace{2.5pt}\mathbb{Q}$-a.s. for all $t\in \left[0,T\right]$ such that the short rate spread is bounded from below – similar to the short rate itself [recall Remark 2.1 (a)]. We interpret ${s_{t}}$ as an additive spread and therefore set for all $t\in \left[0,T\right]$
(cf. [12, 33]) where ${r_{t}}$ denotes the short rate process and ${\overline{r}_{t}}$ is called fictitious short rate, similarly to [26]. With reference to p. 46 in [26], we recall that the short rate spread ${s_{t}}$ not only incorporates credit risks, but also various other risks in the interbank sector which affect the evolution of the LIBOR rates. Let us moreover mention that the short rate ${r_{t}}$ defined in (2.1) and the short rate spread ${s_{t}}$ can be ‘correlated’ by allowing for (at least) one common factor in their respective definitions. More precisely, if the sum in the definition of ${s_{t}}$ started running from $k=n$ (instead of $k=n+1$), then the factor ${X_{t}^{n}}$ would appear both in the definition of ${r_{t}}$ and in the definition of ${s_{t}}$ such that the two latter stochastic processes would no longer be independent.
We next substitute (2.1) as well as the definition of ${s_{t}}$ into (6.1) and deduce
where we introduced the real-valued deterministic function $\overline{\mu }\left(t\right):=\mu \left(t\right)+{\mu ^{\ast }}\left(t\right)$. It obviously holds ${\overline{r}_{t}}\ge \overline{\mu }\left(t\right)\hspace{2.5pt}\mathbb{Q}$-a.s. for all $t\in \left[0,T\right]$. In accordance to Section 3.4.1 in [14], Eq. (2.35) in [26] and Section 1 in [33], we define the fictitious bond price in our post-crisis short rate model via
where $t\in \left[0,T\right]$. The object $\overline{P}\left(t,T\right)$ is sometimes called artificial bond price in the literature, as it is not physically traded in the market. Evidently, $\overline{P}\left(T,T\right)=1$. Comparing (6.2) with (2.1) and (6.3) with (3.3), we see that all our single-curve results presented in the previous sections carry over to the currently considered post-crisis case. More precisely, the entities $\mu \left(t\right)$, ${r_{t}}$, $P\left(t,T\right)$ and n emerging in the previous single-curve equations presently have to be replaced by $\overline{\mu }\left(t\right)$, ${\overline{r}_{t}}$, $\overline{P}\left(t,T\right)$ and l, respectively. Moreover, in the present case, we observe $\mathbb{Q}$-a.s. for all $t\in \left[0,T\right]$
(6.3)
\[ \overline{P}\left(t,T\right):={\mathbb{E}_{\mathbb{Q}}}\left(\exp \left\{-{\int _{t}^{T}}{\overline{r}_{u}}du\right\}\bigg|{\mathcal{F}_{t}}\right)\]
\[ 0<\overline{P}\left(t,T\right)\le \exp \left\{-{\int _{t}^{T}}\overline{\mu }\left(u\right)du\right\}\]
due to (6.3), the lower boundedness of ${\overline{r}_{t}}$ and the monotonicity of conditional expectations.Proof.
Note that taking ${\mu ^{\ast }}\left(t\right)\ge 0\hspace{2.5pt}\forall \hspace{2.5pt}t\in \left[0,T\right]$ implies ${s_{t}}\ge 0\hspace{2.5pt}\mathbb{Q}$-a.s. for all $t\in \left[0,T\right]$. In this case, we deduce ${\overline{r}_{t}}\ge {r_{t}}\hspace{2.5pt}\mathbb{Q}$-a.s. for all $t\in \left[0,T\right]$ due to (6.1). Hence, we find
\[ \exp \left\{-{\int _{t}^{T}}{\overline{r}_{u}}du\right\}\le \exp \left\{-{\int _{t}^{T}}{r_{u}}du\right\}\]
$\mathbb{Q}$-a.s. for all $t\in \left[0,T\right]$. We next take conditional expectations with respect to ${\mathcal{F}_{t}}$ and $\mathbb{Q}$ in the latter inequality, hereafter apply the monotonicity of conditional expectations and finally identify (6.3) and (3.3) in the resulting inequality. □The result obtained in Proposition 6.1 possesses the following economical interpretation: The inequality (6.4) shows that nontraded bonds are cheaper than their nonfictitious counterparts which are physically traded in the market. This feature appears economically reasonable and stands in accordance with the modeling assumptions and statements in [12], Section 2.1.3 in [26] and Section 2.1 in [33]. Moreover, combining (6.3) and (3.18), we obtain (parallel to [12])
\[ \overline{P}\left(t,T\right)=\exp \left\{-{\int _{t}^{T}}\overline{f}\left(t,u\right)du\right\}\]
where $\overline{f}$ is sometimes called fictitious/artificial forward rate in the literature. It holds $\overline{f}\left(t,t\right)={\overline{r}_{t}}$ for all $t\in \left[0,T\right]$. With reference to [11] and [12], for all $t\in \left[0,T\right]$ we introduce the forward rate spread via
so that we have not only set up a new pure-jump post-crisis short rate model, but simultaneously a new pure-jump post-crisis forward rate model of HJM-type in the current section. Recall that (6.4) is equivalent to $f\left(t,u\right)\le \overline{f}\left(t,u\right)\hspace{2.5pt}\mathbb{Q}$-a.s. for all $0\le t\le u\le T$. From this, we conclude that $g\left(t,T\right)\ge 0\hspace{2.5pt}\mathbb{Q}$-a.s. for all $t\in \left[0,T\right]$. It further holds $g\left(t,t\right)={s_{t}}$ for all $t\in \left[0,T\right]$ due to (6.1).Furthermore, in the present post-crisis setting, for a time partition $0\le t\le {T_{1}}<{T_{2}}$ we define the (forward) overnight indexed swap (OIS) rate via
\[ F\left(t,{T_{1}},{T_{2}}\right):=\frac{1}{\delta }\left(\frac{P\left(t,{T_{1}}\right)}{P\left(t,{T_{2}}\right)}-1\right)\]
(cf. Eq. (4.1) in [26]) where P denotes the zero-coupon bond price defined in (3.3) and $\delta :=\delta \left({T_{1}},{T_{2}}\right)$ is the year fraction with expiry date ${T_{1}}$ and maturity date ${T_{2}}$. With reference to [12], for $0\le t\le {T_{1}}<{T_{2}}$ we define the forward LIBOR rate via
\[ L\left(t,{T_{1}},{T_{2}}\right):=\frac{1}{\delta }\left(\frac{\overline{P}\left(t,{T_{1}}\right)}{\overline{P}\left(t,{T_{2}}\right)}-1\right)\]
where δ is the year fraction and $\overline{P}$ denotes the bond price introduced in (6.3). Note that the LIBOR rate $L\left(t,{T_{1}},{T_{2}}\right)$ shall not be mixed up with the Lévy processes ${L_{t}^{k}}$ defined in (2.3). In a pre-crisis single-curve approach, it holds $\overline{P}\left(t,T\right)=P\left(t,T\right)$ which implies $F\left(t,{T_{1}},{T_{2}}\right)=L\left(t,{T_{1}},{T_{2}}\right)\hspace{2.5pt}\mathbb{Q}$-a.s. for all t.We are now prepared to derive the dynamics of the short rate spread ${s_{t}}$, the fictitious short rate ${\overline{r}_{t}}$, the bond price $\overline{P}\left(t,T\right)$, the forward rate $\overline{f}\left(t,T\right)$, the forward rate spread $g\left(t,T\right)$ and the LIBOR rate $L\left(t,{T_{1}},{T_{2}}\right)$. The associated computations can be accomplished by similar techniques as presented in Sections 2 and 3 and thus, are not worked out explicitly. We provide as an example two results without proofs in the sequel. For all $t\in \left[0,T\right]$ it holds
\[ \frac{d\overline{P}\left(t,T\right)}{\overline{P}\left(t-,T\right)}={\overline{r}_{t}}dt+{\sum \limits_{k=1}^{l}}{\int _{{D_{k}}}}{\zeta _{k}}\left(t,T,z\right)d{\tilde{N}_{k}^{\mathbb{Q}}}\left(t,z\right)\]
where the functions ${\zeta _{k}}$ are such as claimed in (3.10). We further obtain in the post-crisis case
\[\begin{aligned}{}L\left(t,{T_{1}},{T_{2}}\right)& =\frac{1}{\delta }\Bigg(\frac{\overline{P}\left(0,{T_{1}}\right)}{\overline{P}\left(0,{T_{2}}\right)}\times {\prod \limits_{k=1}^{l}}\exp \Bigg\{{\int _{0}^{t}}{\int _{{D_{k}}}}{\Psi _{k}}\left(s,z\right)d{\nu _{k}}\left(z\right)ds\\ {} & \hspace{1em}+{\int _{0}^{t}}{\int _{{D_{k}}}}{\Sigma _{k}}\left(s,z\right)d{N_{k}}\left(s,z\right)\Bigg\}-1\Bigg)\end{aligned}\]
with deterministic functions