Rational Savings Account Models for Backward-Looking Interest Rate Benchmarks

Interest rate benchmarks are currently undergoing a major transition. The LIBOR benchmark is planned to be discontinued by the end of 2021 and 'replaced' by what ISDA calls an adjusted risk-free rate (RFR). ISDA has recently announced that the LIBOR 'replacement' will most likely be constructed from a compounded running average of RFR overnight rates over a period matching the LIBOR tenor. This new backward-looking benchmark is markedly different when compared with LIBOR. It is measurable only at the end of the term in contrast to the forward-looking LIBOR, which is measurable at the start of the term. On the other hand though, RFRs provide a simplification because the cash flows and the discount factors may be derived from the same discounting curve, thus avoiding--on a superficial level--any multi-curve complications. We develop a new class of savings account models and derive a novel interest rate system specifically designed to facilitate a high degree of tractability for the pricing of RFR-based fixed-income instruments. The rational form of the savings account models under the risk-neutral measure enables the pricing in closed form of caplets, swaptions and futures written on the backward-looking interest rate benchmark. An interesting twist is that the proposed rational savings account models are different from so-called short rate models in that they cannot necessarily be expressed as an exponentiated integral of a short rate of interest.<br>


Introduction
Since the speech of Andrew Bailey, the Chief Executive of the Financial Conduct Authority, in 2017 on the future of the London interbank offered rate (LIBOR) it has steadily become clearer that the LIBOR interest rate benchmarks are no longer going to be supported in several major currency denominations. The change is primarily justified by the observation that unsecured (term-) lending, the cost of which LIBOR is meant to measure, is no longer a significant source of financing for financial institutions. Furthermore, the current LIBOR benchmark has fallen out favour, i.a., due to its survey-based nature prone to manipulation, see e.g., Gyntelberg and Wooldridge (2008).
The move, much advocated by regulators suchlike the Financial Conduct Authority (FCA), has been towards a benchmark detached from the subjective opinions of a panel, such as the LIBOR panel of global banks, and closer to actual financial transactions. It has been decided, see ISDA (2019), that reliance on LIBOR and most of the other IBOR benchmarks should end and a transaction-based overnight rate-a so-called risk-free rate (RFR)-introduced as the alternative. In the USA, the Secured Overnight Financing Rate (SOFR)-an overnight rate derived from repurchase agreement transactions-is the chosen RFR. A similar benchmark, the reformed Sterling Overnight Index Average (SONIA), has been chosen in the UK, while in the Eurozone the Euro Short-Term Rate (eSTR) is We begin with a general description of the model, and then we present the pricing formulas in a model driven by affine processes. We wrap up with concluding remarks and an outlook on future investigations.

Derivatives on Backward-Looking Rates
Let (Ω, F , Q, (F t ) t∈[0,T] ) be a filtered probability space withT < ∞. The bank account numeraire process (B t ) t≥0 is associated with the risk-neutral measure denoted by Q. It is customary to model the bank account numeraire by where (r t ) t≥0 is a specified instantaneous short rate. Instead, we shall work with a generalised notion of the bank account referred to as an implied savings account in Döberlein and Schweizer (2001). Thus we assume that (B t ) t≥0 is a strictly positive process of finite variation, which generalises (1).
With LIBOR discontinued and superseded by an RFR benchmark, the rate-fixing and discounting are applied by using the same numeraire process (B t ) t≥0 . Hence no features typical of multi-curve interest rate systems emerge, and one is effectively back in a single-curve setting. The standard zero-coupon bond (ZCB) has a price process (P tT ) t≤T given by A European-style derivative with random payoff X T at time T has at time t ≤ T a price given by Letting δ 1d denote the day-count fraction of an overnight loan, it follows from standard interest rate theory that the discrete overnight rate is R 1d (t) := 1 δ 1d (1/P t,t+δ 1d − 1). This is what ISDA refers to as the Risk-Free Rate (RFR), a benchmark rate that may correspond to SOFR, SONIA, eSTR etc., in their respective currencies. The proposed LIBOR replacement is the following rate referred to as the backward-looking benchmark: where T ≤ T 1 , . . . , T n δ ≤ U. A standard approximation applied in the literature, and in particular in Lyashenko and Mercurio (2019), is to replace discrete by continuous compounding and hence usē We will use this approximation throughout the paper. 2 We can compare this rate to a classical discrete forward-looking term rate F(T, U) := (1/P TU − 1)/δ, 2 After consecutive applications of the tower property one sees that, for , and it thus follows that pricing linear derivatives, i.e., any derivative that is linear in the rate, is unaffected by the continuous compounding approximation in the case of t ≤ T ≤ U. again with tenor δ = U − T, i.e., the rate on a loan that starts at time T and has to be repaid at time U. The term rate F(T, U) is based on the same curve as the backward-looking benchmark and can be thought of as LIBOR without the risk premia for credit and liquidity mentioned in the introduction.
In order to account for the risk premia, the actual LIBOR fallback used when converting existing LIBOR contracts will beR(T, U) + c, where c is determined by the historical relationship between LIBOR andR(T, U), see ISDA (2019). For ease of presentation, we shall ignore the spread, and we focus on the implications of changing to a backward-looking benchmark from a classical forward-looking rate. The first thing to notice is that This change, while subtle, has substantial implications for the pricing of interest rate derivatives as we show in the next section.

Swaplet
Let us first consider the expectation of the discounted rate or what could be called a zero-strike swaplet, or one element in the floating leg of a swap. For t ≤ T ≤ U, we have The latter equality reveals that, for t ≤ T ≤ U, backward-looking rates are priced in the same way as the forward-looking equivalent, and by implication it follows that the valuation of linear derivatives remains essentially unaffected by the backward-looking nature of R(T, U). However, for t ∈ [T, U], one has instead Thus when one moves to pricing the derivative after its initiation at t > T the payoff is still undetermined as R(T, U) is only measurable at time U. Relation (2) thus becomes explicitly dependent on the value of the bank account B T , which is not the case in a classic forward-looking setting as given in the Expression (3).

Forward Rate Agreement
The forward rate agreement (FRA) based on a forward-looking benchmark rate, generates a cash flow at time T. In a single-curve setting, this can be considered a linear derivative by noting that for t ≤ T The price of the FRA thus becomes equivalent to the value of a contract with a cash flow δ(F(T, U) − K) at time U. It is not possible to readily substitute F(T, U) with R(T, U) in Equation (4), because the latter is not measurable at time T (aka 'measurability problem'). This is likely to make transitioning FRA contracts from the forward-to the backward-looking benchmarks challenging as discussed in Henrard (2019). The closest backward-looking counterpart to the FRA is a swaplet paying δ(R(T, U) − K) at time U. This would be congruent with the standard FRA for t ≤ T, since However, for t ∈ [T, U] we have The fact that the payoff is random beyond the initiation point T may appear inconsequential, but in practice it matters a great deal, for many risk management systems now have to be reconfigured to address correctly the randomness within the interval [T, U]. See again Henrard (2019) for further discussion.

Caplet
Similarly, a caplet written on the backward-looking rate would presumably pay (R(T, U) − K) + at time U. We can similarly price it for ∀t ≤ U. LetK = 1 + δK, and we have which we can compare with the price of a standard caplet with payoff δ(F(T, While the linear derivatives are priced the same in the backward-and forward-looking setting, at least for t ≤ T, it is clear from the above equation that the caplets are not. In fact, as also shown in Lyashenko and Mercurio (2019) an application of Jensen's inequality to Equation (33) reveals that In other words a caplet written on the backward-looking benchmark is always worth more than or equal to its forward-looking counterpart.

Swaps and Swaptions
An RFR swap can be defined as a contract that has a floating leg paying {R(T i−1 , T i )} {i=a+1,...,b} at times T a+1 , . . . , T b and a fixed leg paying K at the same dates. We can value the swap at time t ≤ T a as We notice that, for t ≤ T a , V a,b t is in fact equivalent to the value of a classical swap paying the discrete forward-looking rate F(T i−1 , T i ) at time T i instead of R(T i−1 , T i ) at time T i . This is not the case for t ∈ [T k−1 , T k ] for some k ∈ (a + 1, . . . , b) where the value of future payments becomes It again demonstrates that the underlying challenge with the backward-looking regime is that one gets explicit dependence on the savings account process once the valuation continues beyond the initiation date of the contract.
A swaption is an option to enter into a swap at time T a , and we can of course price it for t ≤ T a as Since the valuation of the swap is unaffected by whether or not the underlying rate is backwardor forward-looking if priced before the starting time T a , it follows that swaption pricing is unaffected by the shift from forward to backward-looking rates.

Considerations
These preliminary examples illustrate that changing the benchmark from forward-looking to backward-looking rates will alter the pricing and risk management of linear and nonlinear derivatives. In particular (2) and (9) suggest that a modelling paradigm that provides an explicit expression for the bank account (B t ) or more precisely its reciprocal (B −1 t ) is advantageous. This is possible by using an affine short-rate model, via the relation B t = exp( t 0 r s ds), but, as is well-known, this choice does not in general lead to a closed-form solution for swaptions-approximations such as the one in Singleton and Umantsev (2002) have to be applied. Lyashenko and Mercurio (2019) extend the LMM framework to include explicit expressions for the bank account in the time interval [T, U], but the framework suffers from the general tractability issues in the standard LMM, which involve non-trivial issues with Monte-Carlo simulation, and a lack of closed-form pricing for swaptions as well as futures.
A framework that prices swaptions in closed-form is one that results in discounted ZCB prices being linear in the spacial variable. The pricing kernel models by Flesaker and Hughston (1996) and Rutkowski (1997) are established examples thereof. This approach has been further studied and extended in Filipović et al. (2017), who demonstrate that it has excellent empirical properties. Its calibration properties in a multi-curve context are studied in Crépey et al. (2016). Such models of rational form could naturally be applied also to price swaptions on backward-looking rates. However, in the case of derivatives with payoffs requiring direct modelling of the savings account, it turns out that it is not as straightforward as that. Consider the well-known relation where Q is the risk-neutral measure associated with the savings account (numeraire) process (B t ) t≥0 . If the aim is to obtain a rational form for the savings account process (B t ), see Assumption 1 below, then modelling the pricing kernel (π t ) directly may be a tedious path to take, which not necessarily will lead to the desired goal. What we are really interested in is an explicit expression for the process (B t ) t≥0 . Thus, we propose that the savings account process be of rational form, as postulated in Assumption 1 below. We will show that this allows for closed-form pricing under a risk-neutral measure of all derivatives written on a backward-looking interest rate benchmark studied in this section. For the theory underpinning the distinction between an (implied) savings account and a so-called classical savings account (one that may be written in terms of a short rate process), we refer to Döberlein and Schweizer (2001). The discontinuation of interbank offered rate benchmarks has been described by some as the biggest experiment ever in financial markets. These are momentous times while market participants scramble to prepare for the transition to RFR benchmarks. The discontinuation of term-based benchmarks, such as LIBOR, reduces a multiple-curve interest rate system to a single-curve one, akin to the ones in place before the great financial crisis more than a decade ago. As discussed in Backwell et al. (2019), the lack of a term-rate system produces increased exposure to roll-over risk. As a consequence, one can expect a refinancing strategy to become costlier due to the increased roll-over frequency of debt that is rolled-over on an overnight basis. One ought to expect such effects to also impact any so-called 'LIBOR fallback' provisions and in particular the calculation of any adjustment (spread) to legacy contracts written on term-rates, see Mercurio (2018) and Henrard (2019) for example. The adoption of a consistent pricing and hedging approach that models and explains the nature of term rates and the impact of roll-over risk on borrowing strategies appears the right way forward. This would allow for a transparent treatment of term rates and overnight interest rate benchmarks so that a proposed transition to the latter has a chance to be appreciated in its entirety. A lack of term-based benchmark rates may produce a level of opacity in lending markets. Market makers may offer rates, which cannot be 'independently assessed' in the sense that a mark-to-market comparison to a FRA market is not going to be available. For first advancements on consistent approaches to term-and overnight interest rate benchmarks, we here refer to Alfeus et al. (2018), Backwell et al. (2019), and to an article by Macrina & Mahomed (in preparation) with a focus on term risk in interest rate markets. We conclude our reflections by noting that, at the time of writing, there is no global consensus regarding the transition to RFRs. Depending on the jurisdiction, term-based benchmarks are to be superseded by RFRs (e.g., UK), term-based and RFR benchmarks shall coexist until both are accepted as market representatives (e.g., USA, Eurozone). Such ambiguities are perhaps a reflection of the differences between market structures, but also derive from the different natures of the benchmarks themselves that are being used.

Rational Savings Account Models
In the case where one deals with backward-looking rates, we have now shown instances where the payoff itself will explicitly depend on the savings account. It is therefore a clear advantage to directly model the process (B t ) t≥0 . Building on Döberlein and Schweizer (2001), we make the following assumption throughout this paper: Assumption 1. The savings account model (B t ) t≥0 is given by for all t ≥ 0, and where (Y t ) t≥0 is a stochastic, nonnegative process of finite variation. The function a(t) is deterministic, strictly positive, of bounded variation while b(t) is deterministic, nonnegative and of bounded variation. It is assumed they satisfy the initial condition, The financial or economic interpretation of the savings account is no different from that of a classical bank account. The savings account (B t ) t≥0 discounts/accrues cash flows, and, as the name suggests, quantifies the savings made over a period of time. It is also the numeraire associated with an equivalent risk-neutral measure. The difference here is that the savings account model, given by Definition (16), is more general than bank account processes defined in terms of a short-rate model. The rational form of Relation (16) allows for an increased level of tractability, as we will see next, when deriving the price dynamics of interest rate derivatives under a risk-neutral measure Q. Lemma 1. The savings account (B t ) t≥0 , defined by Equation (16), is a strictly positive and finite-variation process.
Proof. Let the sample path for a given ω ∈ Ω be h : t → Y t (ω) and define the mapping g : t → a(t) + b(t)h(t). The functions a, b and h have bounded variation. It follows from the chain rule of bounded variation, see Ambrosio and Dal Masio (1990), that the product (and trivially the sum) of two functions of bounded variation are of bounded variation. Hence, the mapping g has bounded variation. Let f : x → 1/x. Then the savings account is given by the composite function B t (ω) = f • g(t, ω). The function f is Lipschitz-continuous and strictly positive on R >0 . By another application of the chain rule result, it turns out that the composite function B t (ω) is of bounded variation for all ω ∈ Ω and thus that the stochastic process (B t ) t≥0 is strictly positive and has finite variation for all t ≥ 0.
Lemma 1 establishes that the specification in Equation (16) is a savings account in the sense of Döberlein and Schweizer (2001). The probability measure Q, associated with the numeraire process (B t ) t≥0 , is hence a risk-neutral measure.
Remark 1. If any of the conditions underpinning Lemma 1 are relaxed, resulting in (B t ) t≥0 being a positive semi-martingale that is not of finite variation, then (B t ) t≥ is no longer a savings account. It would however still qualify as a numeraire and its associated measure would still be a martingale measure, though not a risk-neutral measure. This means that the true savings account process would no longer be directly specified, thus making the modelling advantage of our approach void.
We are now in the position to write down the price process of a discount bond: where we recall that the conditional expectation is taken with respect to the risk-neutral measure Q.
Assuming that the bond price function P(t, T) is differentiable in T, the instantaneous bond return process (η t ) t≥0 has the form where the dot-notation stands for differentiation with respect to t. It is worth emphasising that the instantaneous bond return may not be given by a short rate, as used in Equation (1) T] is not of finite variation, in general. It follows that the ZCB price process and the instantaneous bond return process may also be of infinite variation. The above expression also reveals that, even in the case of a strictly increasing process (Y t ) t≥0 , the support of (η t ) t≥0 is not limited to a positive domain. Its sign is controlled by the time-derivativesȧ(t) andḃ(t), which may have either sign, and the model can produce positive as well as negative rates. Furthermore, the Form (17) of the discount bond allows for the calibration of the price model to the discount bond initial term-structure (P 0t ) 0≤t≤T . We have where 0 ≤ t ≤ T. Solving for a(t) and inserting the result into Equation (17), one obtains Similar to an HJM-style model, the deterministic function a(t) 0≤t≤T may be automatically calibrated to the initial term structure while b(t) 0≤t≤T and the specific parameters related to the process (Y t ) t≥0 can be used to calibrate to non-linear derivatives.
Assumption 2. The stochastic process (Y t ) t≥0 is taken to be Electronic copy available at: https://ssrn.com/abstract=3523004 where w ∈ R d−1 ≥0 and (Z t ) t≥0 is an affine process, as characterised in the appendix, defined on R d−1 ≥0 and such that w · Z 0 > 0. Döberlein and Schweizer (2001), for a process (Y t ) 0≤t of the form (21)-not necessarily defined in terms of an affine process-such that E[w · Z t ] < ∞ for all t ≥ 0, it follows that the dynamics of the savings account, given by Relation (16), may be written as dB t = r t B t dt, where

Remark 2. In line with Proposition 12 in
is the short rate of interest. This is obtained by applying Ito's quotient rule. It then also follows that the instantaneous return, with Form (18), on the discount bond satisfies η t = r t , 0 ≤ t ≤ T.
It follows from standard theory that (Y t ) t≥0 is strictly positive and of finite variation. Theorem 4.1 in Keller-Ressel (2008) states that if one defines By definition of an affine process, there exist functions φ : for all t ∈ [0, T]. Here, u ∈ S(I ) where and where S(I ) := {z ∈ C d : z ∈ I} is the the strip of existence. We assume that 0 ∈ int(I ) where int denotes the interior of the set I. Furthermore, we will often need the derivatives with respect to the vector u = [u 1 , . . . , u d ] , and so Next, we derive the pricing formulas of the main fixed-income derivative instruments under Assumptions 1 and 2.

Linear Derivatives
We next make use of the savings account process introduced in Assumption 1, and in particular consider the setting in Assumption 2, to derive the no-arbitrage price process of linear derivatives written on a backward-looking benchmark. In the following proposition, we begin with a discount bond.
Electronic copy available at: https://ssrn.com/abstract=3523004 Proposition 1. The discount bond price process is given by Proof. We apply Lemma A4 to obtain where e 1 is a d-dimensional unit-vector with one in the first entry and zero elsewhere. Inserting the above in Equation (17) yields the Equality (27).

Caplets and Swaptions
Proposition 2. The price of a caplet paying (R(T, U) − K) + at time U for any t ≤ U is, where δ = U − T,K = 1 + δK, c 1 = e 1 B T , c 2 = e 1K B U and q(k) : = E t e k(c 0 +c 1 ·X T +c 2 ·X U (32) for any k ∈ C such that kc 1 , kc 2 ∈ S(I ) and where R is any real constant such thatq(R) < ∞.
Proof. Relation (30) follows by direct insertion of the definition of (B t ) and R(T, U) and from rewriting the expression as a function of X T and X U . Result (31) follows directly from Theorem 4 in Filipović et al. (2017), and Equation (33) follows from Lemma A5 in Appendix B.
The following proposition shows that a closed-form formula for swaptions similar to those derived in, e.g., Flesaker and Hughston (1996) and other works on rational models, can be derived in our setup.

Proposition 3. The price of a swaption is given as
Electronic copy available at: https://ssrn.com/abstract=3523004 For all k ∈ C such that kb a,b ∈ S(I ). R > 0 is any real constant such thatq(R) < ∞. (27) into Relation (13) and from rewriting the expression as a function of X T a . Equation (35) follows by an application of Theorem 4 in Filipović et al. (2017).

Futures
Finally, we consider the contract of a future. A standard result, presented in, e.g., Hunt and Kennedy (2004), is that the price of a future contract with a random payoff X, measurable at time U, is priced as the risk-neutral expectation thereof, without discounting, that is, E t [X]. The convention for three-month SOFR/SONIA futures is that their payoff is given by the difference between a notional and the discretely compounded backward-looking rate. Assuming a unit notional, the payoff is 1 −R(T, U), where δ = U − T ≈ three months. As in Mercurio (2018) and define the futures rate for a three-month contract by In a LIBOR market model, this expectation cannot be solved in closed-form; an approximation is developed in Lyashenko and Mercurio (2019). In the following proposition, we show that in our setting no approximation is needed.
Proposition 4. The futures rate for a three-month SOFR futures contract expiring at time U with δ = U − T ≈ three months is given by Then, from which Equation (38) follows. In the third equality, Fubini's theorem is applied. This is valid subject to ψ U−T (−se 1 b(U)) ∈ I for all s ≥ 0. This condition holds because −se 1 b(U) ∈ I for all s ≥ 0.
The convention for one-month SOFR futures is that their payoff is an arithmetic average rate instead of a geometric compounded rate. Following Mercurio (2018), we define the rate of one-month futures expiring at time U, where 1m ≈ δ = U − T, by where R 1m (T, U) is the arithmetic average of the overnight rate over one month. Then, Here, T = T 1 < · · · < T n 1m +1 = U spans the course of one month. This case is not covered in the LMM context of Lyashenko and Mercurio (2019), and it is not clear how one would proceed without further approximations. It is investigated in a short rate model in Mercurio (2018), where the arithmetic average is approximated by an integral over the short rate of interest. This approximation simplifies the expression for Solution (46) in an affine short rate model, but it is in fact not necessary-even for the short rate model they present. As shown below, tractability depends on obtaining an analytical result for the expression E t [1/P TU ], which exists by definition in an exponential affine short rate model. In our savings account setting the resulting expression for the conditional expectation is a little more involved, but nevertheless it only amounts to performing a one-dimensional integration.
Electronic copy available at: https://ssrn.com/abstract=3523004 Similarly, by applying Lemmas A3 and A4 to the second term gives With this result at hand, we can now express the one-month futures rate: Proposition 5.

Example: Multidimensional CIR Process
We conclude with a sketch of an example.
and where Z i t is given by and κ i , θ i , σ i > 0 and (W i t ) is standard Brownian motion independent of (W j t ) for i = j. The conditional moment generating function E t [e X T ] = e φ T−t (u)+ψ T−t (u)·X t for (X t ) t≥0 defined above is found by first noticing that Using the independence property of the processes (Z i t ) t≥0 and applying Lemma A6 gives It follows from Lemma A6 that the result holds for the set {u ∈ R d |u 1 w i < The derivatives needed in the pricing equations can be calculated as follows: The explicit expressions are lengthy and tedious to calculate by hand, but can be obtained with ease by using a software such as Mathematica offering symbolic calculations. In Appendices A and B, we include some results that assist with the mathematical developments in this section.

Conclusions and Outlook
In this paper we develop a new class of rational term structure models that is tailored to the new backward-looking interest rate benchmarks, which are meant to supersede the forward-looking LIBOR. Currently, only the three-month SOFR and SONIA futures can be said to be traded on the new overnight benchmarks, and trading volumes are still rather limited, but nevertheless growing as the LIBOR discontinuation deadline nears. As trade volumes in SOFR/SONIA-based contracts rise, we expect that swaps, caps and swaptions enter the market in a fashion similar to what is described in this paper. A maturing market will generate the needed data also to perform empirical investigations on the validity of the proposed models. At present, not even linear derivatives are traded liquidly, meaning that an initial term structure (P 0t ) t≥0 cannot be implied in a straightforward way. An implied initial term structure would allow for the construction of market-implied, forward-looking term rates via the formula F(t, t + δ) = (1/P t,t+δ − 1)/δ; an important benchmark, much coveted by market participants. The most liquid market instrument referencing RFR rates are the one-month and three-month SOFR futures contracts trading on the Chicago Mercantile Exchange. The non-linearity of the futures contract implies that futures prices do not determine the initial curve (P 0t ) t≥0 in a unique fashion and vice versa, of course. However, model-free attempts to use futures prices to infer the SOFR-curve have been pursued by applying an approximation in Heitfield and Park (2019). The authors envisage to work on a similar procedure, but performed in an arbitrage-and approximation-free manner using the model structure presented in this paper.
A much easier and model-free approach to infer the initial curve would be possible should a liquid market for truly linear derivatives, such as swaps, arise. In this case the initial curve could be constructed using standard single-curve bootstrapping techniques.
It is unclear whether cash market participants will be satisfied with solely using a benchmark free of credit and liquidity risk premia such as SOFR. Alternative indices, such as the ICE Bank Yield Index or AMERIBOR, which include credit risk premia, might soon become increasingly relevant as an alternative to SOFR once LIBOR is discontinued. If this is what the future will bring, then the single-curve framework of this paper would need to be extended to a multi-curve setting in order to price the differences in risk premia across indices. This could be done along the lines of Crépey et al. (2016) and Macrina and Mahomed (2018), properly adjusting for the difference in model structure.

Appendix A. General Results
Lemma A1. Let U be a random variable and assume there exists δ > 0 such that M U (v) := E[e vU ] < ∞ for all v ∈ (−δ, δ). Then Lemma A2. If U > 0 is a random variable then it satisfies Proof. One has that, for x > 0, ∞ 0 e −sx ds = 1/x. The result is produced by applying Fubini's Theorem: Lemma A3. Let U > 0 and W be random variables, and there exists a

Proof. Observe that
∞ 0 e −sx ds = 1/x. The result is obtained by applying Fubini's Theorem:

Appendix B. Affine Processes
Let (Ω, F , Q) be a filtered probability space, where (F t ) 0≤t≤T denotes the associated filtration and 0 <T < ∞. We set We denote by S(I ) the strip S(I ) = {z ∈ C d : z ∈ I}. We define D := {(t, u) ∈ (0,T) × C, u ∈ S(I )}, and we let (X t ) t≥0 be an affine process in the sense of Definition 3.1 in Cuchiero et al. (2011), meaning that there exist functions φ : D → C and ψ : D → C d such that For all (t, u) ∈ D. Since affine processes are also time-homogeneous Markov process, it follows that Lemma A4. For any v 2 ∈ int(I ) and for 0 ≤ t ≤ T ≤T and ∀c 1 , c 2 ∈ R and ∀v 1 ∈ R d . It follows that ) Electronic copy available at: https://ssrn.com/abstract=3523004 Proof. Let X j T and u j be an arbitrary component in X T and u ∈ int(I ), respectively. Then ∂ ∂u j E t e u·X T = E t ∂ ∂u j e u·X T = E t X j T e u·X T .
Differentiating within the expectation hinges, by way of the dominated convergence theorem, on the condition E t [|X j T e u·X T |] < ∞ which can be established by Hölder's inequality. We have E t |X j T e u·X T | ≤ E t |X j T | p 1/p E t e qu·X T 1/q , p > 1, q = p/(p − 1).
Since 0 is in the interior of I, it follows that E[|X j T | p ] < ∞. As p can be chosen arbitrarily large it follows that, since u ∈ int(I ), q can be made sufficiently close to 1 such that qu ∈ I and thus E t [e qu·X T ] < ∞. From Inequality (A8) one then has E t [|X j T e u·X T |] < ∞. Applying the same arguments for all j = 1, . . . , d and using Equations (A5) and (A7) we get from which Equation (A6) follows.
Electronic copy available at: https://ssrn.com/abstract=3523004 Proof. The expressions for A t (u 1 , u 2 ) and B t (u 1 , u 2 ) are stated to hold for u 1 , u 2 ≤ 0 in Proposition 6.3.4.1 in Jeanblanc et al. (2009). The functions A t (u 1 , u 2 ) and B t (u 1 , u 2 ) are verified to be analytical also for u 1 < k 2 /(2σ 2 ) and u 2 ≤ κ/σ 2 . To verify that the analytical extension remains the unique moment generating function in the extended domain, Keller-Ressel and Mayerhofer (2015) show that a sufficient and necessary condition is that the functions A t (u 1 , u 2 ) and B t (u 1 , u 2 ) are solutions to the generalized Riccati equations associated with the affine process ( t 0 Z s ds, Z t ) t≥0 . By standard theory, see for example Duffie et al. (2000), these are given by dA t = κθB t , A 0 (0) = 0.
Inserting the expressions A t (u 1 , u 2 ) and B t (u 1 , u 2 ) above and verifying that they are indeed solutions completes the proof.