Modern Stochastics: Theory and Applications logo


  • Help
Login Register

  1. Home
  2. Issues
  3. Volume 7, Issue 2 (2020)
  4. A pure-jump mean-reverting short rate mo ...

A pure-jump mean-reverting short rate model
Volume 7, Issue 2 (2020), pp. 113–134
Markus Hess  

Authors

 
Placeholder
https://doi.org/10.15559/20-VMSTA152
Pub. online: 20 April 2020      Type: Research Article      Open accessOpen Access

Received
23 April 2019
Revised
27 March 2020
Accepted
27 March 2020
Published
20 April 2020

Abstract

A new multi-factor short rate model is presented which is bounded from below by a real-valued function of time. The mean-reverting short rate process is modeled by a sum of pure-jump Ornstein–Uhlenbeck processes such that the related bond prices possess affine representations. Also the dynamics of the associated instantaneous forward rate is provided and a condition is derived under which the model can be market-consistently calibrated. The analytical tractability of this model is illustrated by the derivation of an explicit plain vanilla option price formula. With view on practical applications, suitable probability distributions are proposed for the driving jump processes. The paper is concluded by presenting a post-crisis extension of the proposed short and forward rate model.

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
(2.1)
\[ {r_{t}}:=\mu \left(t\right)+{\sum \limits_{k=1}^{n}}{X_{t}^{k}}\]
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
(2.2)
\[ d{X_{t}^{k}}=-{\lambda _{k}}{X_{t}^{k}}dt+{\sigma _{k}}d{L_{t}^{k}}\]
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
(2.3)
\[ {L_{t}^{k}}:={\int _{0}^{t}}{\int _{{D_{k}}}}zd{N_{k}}\left(s,z\right)\]
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
(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\]
which constitute $\left(\mathcal{F},\mathbb{Q}\right)$-martingale integrators. Herein, the positive and σ-finite Lévy measures ${\nu _{k}}$ satisfy the integrability conditions
(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 \]
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
\[ {\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
(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)\]
where we used (2.3). The representation (2.6) implies
(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)\]
where $0\le t\le T$. For all $t\in \left[0,T\right]$ we next define the historical filtration
\[ {\mathcal{F}_{t}}:=\sigma \{{L_{s}^{1}},\dots ,{L_{s}^{n}}:0\le s\le t\}.\]
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
\[ {\Phi _{{r_{t}}}}\left(u\right):={\mathbb{E}_{\mathbb{Q}}}\big[{e^{iu{r_{t}}}}\big]\]
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
\[ {\Phi _{{X_{t}^{k}}}}\left(u\right)={e^{{\psi _{k}}\left(t,u\right){x_{k}}+{\rho _{k}}\left(t,u\right)}}\]
with deterministic and affine characteristic exponent.
An immediate consequence of Proposition 2.3 is the subsequent affine representation
(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)}}\]
where we introduced the deterministic functions
\[ {\phi _{k}}\left(t,u\right):={\rho _{k}}\left(t,u\right)+iu\mu \left(t\right)/n.\]
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
(2.9)
\[ {\kappa _{{r_{t}}}}\left(v\right):={\mathbb{E}_{\mathbb{Q}}}\left[{e^{v{r_{t}}}}\right]\]
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.
Proposition 2.4.
For all $t\in \left[0,T\right]$ the short rate process follows the dynamics
(2.10)
\[ d{r_{t}}=\left({\mu ^{\prime }}\left(t\right)-{\sum \limits_{k=1}^{n}}{\lambda _{k}}{X_{t}^{k}}\right)dt+{\sum \limits_{k=1}^{n}}{\sigma _{k}}{\int _{{D_{k}}}}zd{N_{k}}\left(t,z\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
(3.1)
\[ d{\beta _{t}}={r_{t}}{\beta _{t}}dt\]
with normalized initial capital ${\beta _{0}}=1$. The solution of (3.1) reads as
(3.2)
\[ {\beta _{t}}=\exp \left\{{\int _{0}^{t}}{r_{s}}ds\right\}\]
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
(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)\]
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
(3.4)
\[ P\left(t,T\right)\le {M_{t,T}}:=\exp \left\{-{\int _{t}^{T}}\mu \left(s\right)ds\right\}\]
$\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$.
Proposition 3.1.
For $k\in \left\{1,\dots ,n\right\}$ and $t\in \left[0,T\right]$ we define the deterministic functions
(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}\]
Then the bond price at time $t\le T$ with maturity T possesses the affine representation
(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}}}}\]
where the factors ${X_{t}^{k}}$ satisfy (2.7).
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
(3.7)
\[ {A_{k}}\left(t,t\right)={B_{k}}\left(t,t\right)=0.\]
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
(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\}\]
which implies $P\left(T,T\right)=1$ due to (3.7). Moreover, from (3.5) we infer the time derivatives
(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\]
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)
\[\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
(3.10)
\[ {\zeta _{k}}\left(t,T,z\right):={e^{{\sigma _{k}}{B_{k}}\left(t,T\right)z}}-1\]
with ${B_{k}}\left(t,T\right)$ as in (3.5). Then the bond price satisfies the t-dynamics under $\mathbb{Q}$
(3.11)
\[ \frac{dP\left(t,T\right)}{P\left(t-,T\right)}={r_{t}}dt+{\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).\]
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
(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}\]
where the initial value $P\left(0,T\right)$ is deterministic and fulfills $P\left(0,T\right)>0$.
Furthermore, for all $t\in \left[0,T\right]$ let us introduce the discontinuous Doléans-Dade exponential
(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}\]
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.
Corollary 3.4.
For all $0\le t\le T$ the bond price satisfies
(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}}\]
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.
Moreover, for all $0\le t\le T$ we define the discounted bond price
(3.15)
\[ \hat{P}\left(t,T\right):=\frac{P\left(t,T\right)}{{\beta _{t}}}\]
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
(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}}\]
where $P\left(0,T\right)$ is deterministic and ${\xi _{k}}$ is such as defined in Corollary 3.4. We obtain the following result.
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
(3.17)
\[ f\left(t,T\right):=-{\partial _{T}}\log P\left(t,T\right)\]
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
(3.18)
\[ P\left(t,T\right)=\exp \left\{-{\int _{t}^{T}}f\left(t,u\right)du\right\}.\]
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
(3.19)
\[ {\partial _{T}}{B_{k}}\left(t,T\right)=-{e^{-{\lambda _{k}}\left(T-t\right)}}\]
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
(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)}}\]
where the factor processes ${X_{t}^{k}}$ satisfy (2.7) and ${B_{k}}\left(s,T\right)$ is like defined in (3.5).
Proof.
We substitute (3.8) into (3.17) and obtain
\[ f\left(t,T\right)=-{\sum \limits_{k=1}^{n}}\big[{\partial _{T}}{A_{k}}\left(t,T\right)+{X_{t}^{k}}{\partial _{T}}{B_{k}}\left(t,T\right)\big].\]
Combining this equality with Lemma 3.6, we derive the claimed representation (3.20).  □
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
(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\}\]
where the deterministic initial value is given by $f\left(0,T\right)=-{\partial _{T}}\log P\left(0,T\right)$.
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
(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)\]
for all maturity times $T>0$.
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
(3.23)
\[ R\left(t,T\right):=\frac{\log P\left(t,T\right)}{t-T}.\]
This object is called continuously-compounded spot rate on p. 60 in [7]. It obviously holds
(3.24)
\[ P\left(t,T\right)={e^{-\left(T-t\right)R\left(t,T\right)}}\]
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
(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)\]
where the deterministic functions ${B_{k}}$ are such as defined in (3.5).
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
\[ {\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=-{\int _{0}^{t}}{\int _{{D_{k}}}}{\sigma _{k}}{B_{k}}\left(u,t\right)zd{N_{k}}\left(u,z\right),\]
so that the proof is complete.  □
Recall that the last jump integral in (3.25) constitutes a so-called Volterra integral, as the time parameter t appears both inside the integrand and inside the upper integration bound. Also note that it holds ${I_{t}}=\log {\beta _{t}}$ with ${I_{0}}=0$ due to (3.2).

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
(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)\]
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
(4.2)
\[ {H_{\tau }}={\left[P\left(\tau ,T\right)-K\right]^{+}}:=\max \left\{0,P\left(\tau ,T\right)-K\right\}.\]
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
\[ \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
(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}\]
where the integrated short rate process ${I_{t}}$ is such as defined in (3.25) and
(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}\]
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.
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
\[ q\left(u\right):={e^{-au}}{\left[P\left(0,T\right){e^{u}}-K\right]^{+}}\]
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
\[ {Z_{t,\tau }}:={I_{t}}-{I_{\tau }}+\left(a+iy\right){G_{\tau }}\]
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
\[ {C_{0}}={\int _{\mathbb{R}}}\hat{q}\left(y\right)\exp \left\{\theta \left(\tau ,y\right)+{\sum \limits_{k=1}^{n}}{\overline{\psi }_{k}}\left(0,\tau ,y\right)\right\}dy\]
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
(5.1)
\[ {L_{t}^{k}}={\sum \limits_{j=1}^{{N_{t}^{k}}}}{Y_{j}^{k}}\]
(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
\[ {c_{k}}{\alpha _{k}}={\int _{{D_{k}}}}zd{\nu _{k}}\left(z\right)\]
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
(5.2)
\[ d{\nu _{k}}\left(z\right)={\alpha _{k}}{\varepsilon _{k}}{e^{-{\varepsilon _{k}}z}}dz\]
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.
Successively applying the definition of the characteristic function, (2.3), the Lévy–Khinchin formula and (5.2), for $u\in \mathbb{R}$ and $t\in \left[0,T\right]$ we obtain
\[ {\Phi _{{L_{t}^{k}}}}\left(u\right)={\mathbb{E}_{\mathbb{Q}}}\big[{e^{iu{L_{t}^{k}}}}\big]=\exp \left\{{\alpha _{k}}{\varepsilon _{k}}t{\int _{0}^{\infty }}\big[{e^{iuz}}-1\big]{e^{-{\varepsilon _{k}}z}}dz\right\}.\]
We eventually perform the integration and end up with the asserted equality.  □
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\}$
\[ {\kappa _{{L_{t}^{k}}}}\left(v\right)={\Phi _{{L_{t}^{k}}}}\left(-iv\right)=\exp \left\{\frac{v{\alpha _{k}}t}{{\varepsilon _{k}}-v}\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
\[ {f_{{L_{t}^{k}}}}\left(x\right)=\frac{1}{2\pi }{\int _{0}^{\infty }}\exp \left\{iu\left(x-\frac{{\alpha _{k}}t}{{\varepsilon _{k}}+iu}\right)\right\}du.\]
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
\[ {\kappa _{{r_{t}}}}\left(v\right)={\Phi _{{r_{t}}}}\left(-iv\right)=\exp \left\{v\mu \left(t\right)+{\sum \limits_{k=1}^{n}}{\rho _{k}}\left(t,-iv\right)+{\sum \limits_{k=1}^{n}}{\psi _{k}}\left(t,-iv\right){x_{k}}\right\}\]
with deterministic functions
\[ {\psi _{k}}\left(t,-iv\right)=v{e^{-{\lambda _{k}}t}},\hspace{2em}{\rho _{k}}\left(t,-iv\right)=\frac{{\alpha _{k}}}{{\lambda _{k}}}\log \left|\frac{{\varepsilon _{k}}-v{\sigma _{k}}{e^{-{\lambda _{k}}t}}}{{\varepsilon _{k}}-v{\sigma _{k}}}\right|.\]
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
\[ {s_{t}}:={\mu ^{\ast }}\left(t\right)+{\sum \limits_{k=n+1}^{l}}{X_{t}^{k}}\]
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]$
(6.1)
\[ {\overline{r}_{t}}:={r_{t}}+{s_{t}}\]
(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
(6.2)
\[ {\overline{r}_{t}}=\overline{\mu }\left(t\right)+{\sum \limits_{k=1}^{l}}{X_{t}^{k}}\]
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
(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)\]
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]$
\[ 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.
Proposition 6.1.
It holds $\mathbb{Q}$-a.s. for all $t\in \left[0,T\right]$
(6.4)
\[ \overline{P}\left(t,T\right)\le P\left(t,T\right)\]
where $\overline{P}\left(t,T\right)$ constitutes the bond price defined in (6.3) and $P\left(t,T\right)$ is given in (3.3).
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
\[ g\left(t,T\right):=\overline{f}\left(t,T\right)-f\left(t,T\right)\]
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
\[\begin{aligned}{}{\Psi _{k}}\left(s,z\right)& :={e^{{\sigma _{k}}{B_{k}}\left(s,{T_{2}}\right)z}}-{e^{{\sigma _{k}}{B_{k}}\left(s,{T_{1}}\right)z}}<0,\\ {} {\Sigma _{k}}\left(s,z\right)& :={\sigma _{k}}\left[{B_{k}}\left(s,{T_{1}}\right)-{B_{k}}\left(s,{T_{2}}\right)\right]z>0.\end{aligned}\]

References

[1] 
Antonov, A., Konikov, M., Spector, M.: The free boundary SABR: natural extension to negative rates. Risk Magazine (September 2015).
[2] 
Benth, F., Kallsen, J., Meyer-Brandis, T.: A non-Gaussian Ornstein-Uhlenbeck Process for Electricity Spot Price Modeling and Derivatives Pricing. Appl. Math. Finance 14(2), 153–169 (2006). MR2323278. https://doi.org/10.1080/13504860600725031
[3] 
Benth, F., Saltyte-Benth, J., Koekebakker, S.: Stochastic Modeling of Electricity and Related Markets, 1st edn. World Scientific, (2008). MR2416686. https://doi.org/10.1142/9789812812315
[4] 
Björk, T., Kabanov, Y., Runggaldier, W.: Bond market structure in the presence of marked point processes. Math. Finance 7, 211–239 (1997). MR1446647. https://doi.org/10.1111/1467-9965.00031
[5] 
Brace, A., Gatarek, D., Musiela, M.: The market model of interest rate dynamics. Math. Finance 7, 127–155 (1997). MR1446645. https://doi.org/10.1111/1467-9965.00028
[6] 
Brigo, D., Mercurio, F.: A deterministic-shift extension of analytically-tractable and time-homogeneous short-rate models. Finance Stoch. 5(3), 369–387 (2001). MR1849427. https://doi.org/10.1007/PL00013541
[7] 
Brigo, D., Mercurio, F.: Interest Rate Models – Theory and Practice, 1st edn. Springer, (2001). MR1846525. https://doi.org/10.1007/978-3-662-04553-4
[8] 
Carr, P., Madan, D.: Option valuation using the fast Fourier transform. J. Comput. Finance 2(4), 61–73 (1999). https://doi.org/10.21314/JCF.1999.043
[9] 
Cont, R., Tankov, P.: Financial Modeling with Jump Processes, 1st edn. Chapman & Hall/CRC, (2004). MR2042661
[10] 
Cox, J., Ingersoll, J., Ross, S.: A theory of the term structure of interest rates. Econometrica 53, 385–407 (1985). MR0785475. https://doi.org/10.2307/1911242
[11] 
Crépey, S., Grbac, Z., Ngor, N., Skovmand, D.: A Lévy HJM multiple-curve model with application to CVA computation. Quant. Finance 15(3), 401–419 (2015). MR3313198. https://doi.org/10.1080/14697688.2014.942232
[12] 
Crépey, S., Grbac, Z., Nguyen, H.: A multiple-curve HJM model of interbank risk. Math. Financ. Econ. 6, 155–190 (2013). MR2966729. https://doi.org/10.1007/s11579-012-0083-4
[13] 
Cuchiero, C., Fontana, C., Gnoatto, A.: A General HJM Framework for Multiple Yield Curve Modeling. Finance Stoch. 20(2), 267–320 (2016). MR3479323. https://doi.org/10.1007/s00780-016-0291-5
[14] 
Cuchiero, C., Fontana, C., Gnoatto, A.: Affine Multiple Yield Curve Models. Math. Finance 29(2), 568–611 (2019). preprint: arXiv:1603.00527v2. MR3925431. https://doi.org/10.1111/mafi.12183
[15] 
Di Nunno, G., Øksendal, B., Proske, F.: Malliavin Calculus for Lévy Processes with Applications to Finance, 1st edn. Springer, (2009). MR2460554. https://doi.org/10.1007/978-3-540-78572-9
[16] 
Duffie, D., Filipovic, D., Schachermayer, W.: Affine Processes and Applications in Finance. Ann. Appl. Probab. 13, 984–1053 (2003). MR1994043. https://doi.org/10.1214/aoap/1060202833
[17] 
Eberlein, E., Gerhart, C., Grbac, Z.: Multiple curve Lévy forward price model allowing for negative interest rates. Mathematical Finance (2019). forthcoming, preprint: arXiv:1805.02605v1. https://doi.org/10.1007/978-3-030-26106-1
[18] 
Eberlein, E., Kluge, W.: Exact pricing formulae for caps and swaptions in a Lévy term structure model. J. Comput. Finance 9(2), 99–125 (2006). MR2359367. https://doi.org/10.21314/JCF.2005.158
[19] 
Eberlein, E., Kluge, W.: Valuation of floating range notes in Lévy term structure models. Math. Finance 16, 237–254 (2006). MR2212265. https://doi.org/10.1111/j.1467-9965.2006.00270.x
[20] 
Eberlein, E., Kluge, W., Papapantoleon, A.: Symmetries in Lévy Term Structure Models. Int. J. Theor. Appl. Finance 9(6), 967–986 (2006). MR2260054. https://doi.org/10.1142/S0219024906003809
[21] 
Filipovic, D.: Term-Structure Models. Springer Finance, Springer, (2009). MR2553163. https://doi.org/10.1007/978-3-540-68015-4
[22] 
Filipovic, D., Trolle, A.: The Term Structure of Interbank Risk. J. Financ. Econ. 109(3), 707–733 (2013). https://doi.org/10.1016/j.jfineco.2013.03.014
[23] 
Gallitschke, J., Seifried, S., Seifried, F.: Post-crisis Interest Rates: Xibor Mechanics and Basis Spreads. J. Bank. Finance 78, 142–152 (2017). https://doi.org/10.1016/j.jbankfin.2017.01.002
[24] 
Grasselli, M., Miglietta, G.: A Flexible Spot Multiple-curve Model. Quant. Finance 6(10), 1465–1477 (2016). MR3564921. https://doi.org/10.1080/14697688.2015.1108521
[25] 
Grbac, Z., Papapantoleon, A., Schoenmakers, J., Skovmand, D.: Affine LIBOR Models with Multiple Curves: Theory, Examples and Calibration. SIAM J. Financ. Math. 6(1), 984–1025 (2015). MR3413586. https://doi.org/10.1137/15M1011731
[26] 
Grbac, Z., Runggaldier, W.: Interest Rate Modeling: Post-Crisis Challenges and Approaches, Springer Briefs in Quantitative Finance, 1st edn. Springer, (2015). MR3467166. https://doi.org/10.1007/978-3-319-25385-5
[27] 
Heath, D., Jarrow, R., Morton, A.: Bond pricing and the term structure of interest rates: a new methodology for contingent claims valuation. Econometrica 60, 77–105 (1992). https://doi.org/10.2307/2951677
[28] 
Hess, M.: Pricing Energy, Weather and Emission Derivatives under Future Information. PhD Dissertation, University Duisburg-Essen, Germany (2013). http://duepublico.uni-duisburg-essen.de/servlets/DocumentServlet?id=31060
[29] 
Hull, J., White, A.: Pricing interest-rate-derivative securities. Rev. Financ. Stud. 3(4), 573–592 (1990). https://doi.org/10.1093/rfs/3.4.573
[30] 
Jacod, J., Shiryaev, A.: Limit Theorems for Stochastic Processes, 2nd edn. Springer, (2003). MR1943877. https://doi.org/10.1007/978-3-662-05265-5
[31] 
Keller-Ressel, M., Papapantoleon, A., Teichmann, J.: The Affine LIBOR Models. Math. Finance 23(4), 627–658 (2013). MR3094715. https://doi.org/10.1111/j.1467-9965.2012.00531.x
[32] 
Mercurio, F.: A LIBOR Market Model with a Stochastic Basis, Risk, pp. 96–101 (December 2013)
[33] 
Morino, L., Runggaldier, W.: On multicurve models for the term structure. In: Nonlinear Economic Dynamics and Financial Modeling, pp. 275–290. Springer, (2014). MR3329961
[34] 
Nguyen, T., Seifried, F.: The Multi-curve Potential Model. Int. J. Theor. Appl. Finance 18(7), 1550049 (2015). MR3423186. https://doi.org/10.1142/S0219024915500491
[35] 
Protter, P.: Stochastic Integration and Differential Equations, 2nd edn. Springer, (2005). MR2020294
[36] 
Sato, K.: Lévy Processes and Infinitely Divisible Distributions. Cambridge studies in advanced mathematics, vol. 68 (1999). MR1739520
[37] 
Schoutens, W.: Lévy Processes in Finance: Pricing Financial Derivatives. John Wiley & Sons, Ltd., (2003)
[38] 
Vasicek, O.: An equilibrium characterization of the term structure. J. Financ. Econ. 5(2), 177–188 (1977). https://doi.org/10.1016/0304-405X(77)90016-2
Exit Reading PDF XML


Table of contents
  • 1 Introduction
  • 2 A pure-jump multi-factor short rate model
  • 3 Bond prices and forward rates
  • 4 Option pricing
  • 5 Practical applications
  • 6 A post-crisis model extension
  • References

Export citation

Copy and paste formatted citation
Placeholder

Download citation in file


Share


RSS

MSTA

MSTA

  • Online ISSN: 2351-6054
  • Print ISSN: 2351-6046
  • Copyright © 2018 VTeX

About

  • About journal
  • Indexed in
  • Editors-in-Chief

For contributors

  • Submit
  • OA Policy
  • Become a Peer-reviewer

Contact us

  • ejournals-vmsta@vtex.lt
  • Mokslininkų 2A
  • LT-08412 Vilnius
  • Lithuania
Powered by PubliMill  •  Privacy policy