Pricing FX Options in the Heston / CIR Jump-Diffusion Model with Log-Normal and Log-Uniform Jump Amplitudes

We examine foreign exchange options in the jump-diffusion version of the Heston stochastic volatility model for the exchange rate with log-normal jump amplitudes and the volatility model with log-uniformly distributed jump amplitudes. We assume that the domestic and foreign stochastic interest rates are governed by the CIR dynamics. The instantaneous volatility is correlated with the dynamics of the exchange rate return, whereas the domestic and foreign short-term rates are assumed to be independent of the dynamics of the exchange rate and its volatility. The main result furnishes a semianalytical formula for the price of the foreign exchange European call option.


Introduction
We extend the results from Ahlip and Rutkowski [1] by deriving a closed-form pricing formula for the foreign exchange (FX) options in a model where the spot exchange rate and its volatility are jump-diffusions, with log-normal and loguniform jump amplitudes, respectively, whereas the domestic and foreign interest rates are governed by the Cox-Ingersoll-Ross (CIR) dynamics postulated in [2].In particular, our model allows for correlation between the exchange rate process and its instantaneous volatility.The interest rate processes are independent of one another, and they are also independent of the foreign exchange rate and its volatility.
In the seminal paper by Heston [3], the author noted that increasing the volatility of volatility only increases the kurtosis of spot returns and does not capture skewness.In order to capture the skewness, it is crucial to include also the properly specified correlation between the volatility and the spot exchange rate returns.In papers by Bakshi et al. [4], Bates [5], and Duffie et al. [6], the authors showed that stochastic volatility models do not offer reliable prices for close to expiration derivatives.This motivated Bates [5] and Bakshi et al. [4] to introduce jumps to the dynamics of the exchange rate.However, as observed by Andersen and Andreasen [7] and Alizadeh et al. [8], the addition of jumps to the dynamics of the exchange rate is not sufficient to capture the sudden increase in volatility due to market turbulence.Since the overall volatility in financial markets consists of a highly persistent slow moving and rapid moving components, Eraker et al. [9] proposed to introduce jump process to the dynamics of the volatility process in order to enhance the cross-sectional impact on option prices.
More recently, D'Ippoliti et al. [10] obtained closed-form solutions, in the spirit of Heston, in a model with jumps in both spot returns of the underlying asset and its volatility.Yan and Hanson [11] consider a model in which the stock prices follow a jump-diffusion process with log-uniformly distributed jump amplitudes under the Heston volatility model.In the above-mentioned papers, the authors assume constant interest rates.Although the assumption of constant interest rates leads to highly tractable FX models, empirical results have confirmed that such models do not reflect the market reality, especially for long-dated hybrid FX products.For these products, the fluctuations of both the exchange rate 2 International Journal of Stochastic Analysis and the interest rates are critical, so that the constant interest rates assumption is clearly inappropriate for reliable valuation and hedging.
Let us comment briefly on the existing literature in the same vein.van Haastrecht et al. [12] have extended the stochastic volatility model of Schobel and Zhu [13] to equity/currency derivatives by including stochastic interest rates and assuming all driving model factors to be instantaneously correlated.Since their model is based on the Gaussian processes, it enjoys analytical tractability even in the most general case of a full correlation structure.By contrast, when the squared volatility is driven by the CIR process and the interest rate is driven either by the Vasicek [14] or the Cox et al. [2] process, a full correlation structure leads to intractability of equity options even under a partial correlation of the driving factors, as have been documented by, among others, van Haastrecht and Pelsser [15] and Grzelak and Oosterlee [16,17], who examined, in particular, the Heston/Vasicek and Heston/CIR hybrid models (see also Grzelak et al. [18], where the Schobel-Zhu/Hull-White and Heston/Hull-White models for equity derivatives are studied).
Our goal is to derive semianalytical solutions for prices of plain-vanilla FX options in a model in which the instantaneous volatility component is specified by the extended Heston model with log-normally and log-uniformly distributed jump amplitudes for the exchange rate and the volatility process, respectively, whereas the short-term interest rates for the domestic and foreign economies are governed by the independent CIR processes.The model thus incorporates important empirical characteristics of exchange rate return variability: (a) the correlation between the exchange rate and its stochastic volatility, (b) the presence of jumps in the exchange rate and volatility processes, and (c) the random character of interest rates.The practical importance of this feature of newly developed FX models is rather clear in view of the existence of complex FX products that have a long lifetime and are sensitive to smiles or skews in the market.The results obtained in this paper extend results obtained by Guoqing et al.In their model only the stock price process is subject to jumps, but the volatility of volatility is modeled by the Heston dynamics.
The paper is organized as follows.In Section 2, we set the foreign exchange model examined in this work.The options pricing problem is introduced in Section 3. The main result, Theorem 3 of Section 4, furnishes the pricing formula for FX options.It is worth stressing that the independence of volatility and interest rates appears to be a crucial assumption from the point of view of analytical tractability and thus it cannot be relaxed.Numerical illustrations of our method are provided in Section 5 where the diffusion and jump-diffusion models are compared.

The Heston/CIR Jump-Diffusion Foreign Exchange Model
Let (Ω, F, P) be an underlying probability space.Let the exchange rate  = (  ) ∈[0,] , its instantaneous squared volatility V = (V  ) ∈[0,] , the domestic short-term interest rates  = (  ) ∈[0,] , and the foreign short-term interest rate r = (r  ) ∈[0,] be governed by the following system of SDEs: We work under the following standing assumptions: (A.1) Processes   = (   ) ∈[0,] and  V = ( V  ) ∈[0,] are correlated Brownian motions with a constant correlation coefficient, so that the quadratic covariation between the processes   and ) ∈[0,] are independent Brownian motions and they are also independent of the Brownian motions   and  V ; hence the processes , , and r are independent.(A.5)The Poisson processes   ,  V and sequences of random variables (   ) ∞ =1 and ( V  ) ∞ =1 are mutually independent, as well as independent of the Brownian motions   ,  V ,   ,   .
Note that we postulate that the instantaneous squared volatility process V, the domestic short-term interest rate , and the foreign interest rate r are independent stochastic processes.We will argue in what follows that this assumption is indeed crucial for analytical tractability.For brevity, we refer to the foreign exchange model given by SDEs (1) under Assumptions (A.1)-(A.6)as the Heston/CIR jump-diffusion FX model.

Foreign Exchange Call Option
We will first establish the general representation for the value of the foreign exchange (i.e., currency) European call option with maturity  > 0 and a constant strike level  > 0. The probability measure P is interpreted as the domestic spot martingale measure (i.e., the domestic risk-neutral probability).We denote by F = (F  ) ∈[0,] the filtration generated by the Brownian motions   ,  V ,   ,   and the compound Poisson processes   and  V .We write E P  (⋅) and P  (⋅) to denote the conditional expectation and the conditional probability under P with respect to the field F  , respectively.In our computations, we will adopt the "domestic" point of view, which will frequently be represented by the subscript .Similarly, we will use the subscript  when referring to a foreign denominated variable.Hence the arbitrage price   (, ) of the foreign exchange call option at time  ∈ [0, ] is given as the conditional expectation with respect to the -field F  of the option's payoff at expiration discounted by the domestic money market account; that is, or, equivalently, Similarly, the arbitrage price of the domestic discount bond maturing at time  equals, for every  ∈ [0, ], and an analogous formula holds for the price process   (, ) of the foreign discount bond under the foreign spot martingale measure (see, e.g., Chapter 14 in Musiela and Rutkowski [20]).As a preliminary step towards the general valuation result presented in Section 4, we state the following well-known proposition (see, e.g., Cox et al. [2] or Chapter 10 in Musiela and Rutkowski [20]).It is worth stressing that we use here, in particular, the postulated independence of the foreign interest rate r and the exchange rate process .Under this standing assumption, the dynamics of the foreign bond price   (, ) under the domestic spot martingale measure P can be seen as an immediate consequence of formula (14.3) in Musiela and Rutkowski [20].The simple form of the dynamics of   (, ) under P is a consequence of the postulated independence of   and   (see Assumption (A.2)).This crucial feature underpins our further calculations and thus it cannot be easily relaxed.
The following result is also well known (see, e.g., Section 14.1.1 in Musiela and Rutkowski [20]).

Pricing Formula for the FX Call Option
We are in a position to state the main result of the paper, which furnishes a semianalytical formula for the arbitrage price of the FX call option of European style under the Heston stochastic volatility for the exchange rate combined with the independent CIR models for the domestic and foreign short-term rates.Since the proof of Theorem 3 relies on the derivation of the conditional characteristic function of the logarithm of the exchange rate, any suitable version of the Fourier inversion technique or simulation technique can be International Journal of Stochastic Analysis applied to obtain the option price.The interested reader is referred to, for instance, Carr and Madan [21,22] or Lord and Kahl [23,24] and the references therein, as well as the recent papers by Bernard et al. [25] and Levendorskii [26] who developed and examined in detail methods with essential improvements in accuracy and/or efficiency.

Lemma 4. (i) Under Assumptions (A.3) and (A.5), the following equalities are valid:
(ii) Under Assumptions (A.4) and (A.5), the following equalities are valid for  =  +  with  ≤ 0: The next result extends Lemma 6.1 in Ahlip and Rutkowski [28] (see also Duffie et al. [6]) where the model without the jump component in the dynamics of V was examined.
Lemma 5. Let the dynamics of processes V, , and r be given by SDEs (1) with independent Brownian motions  V ,   , and   .For any complex numbers , , μ, λ, μ, λ, we set Then where where one denotes Proof.For the reader's convenience, we sketch the proof of the lemma.Let us set, for  ∈ [0, ], International Journal of Stochastic Analysis Then the process  = (  ) ∈[0,] satisfies and thus it is an F-martingale under P.By applying the Itô formula to the right-hand side in (25) and by setting the drift term in the dynamics of  to be zero, we deduce that the function (, V, , r) satisfies the following partial integrodifferential equation (PIDE): with the initial condition (0, V, , r) = exp(−V − λ − λr).
We search for a solution to this PIDE in the form with By substituting this expression in the PIDE and using part (ii) in Lemma 4, we obtain the following system of ODEs for the functions  1 ,  2 ,  3 ,  1 ,  2 ,  3 (for brevity, we suppress the last three arguments): By solving these equations, we obtain the stated formulae.
Under the assumptions of Lemma 5, it is possible to factorize  as a product of two conditional expectations.This means that the functions  1 ( 1 ),  2 ( 2 ), and  3 ( 3 ) are of the same form, except that they correspond to different sets of parameters, , ,  V for  1 ,  1 ,   ,   ,   for  2 ,  2 , and   ,   ,   for  3 ,  3 .Note, however, that the roles played by the processes V, , and r in our model are clearly different.
It should also be stressed that no closed-form analytical expression for (, V  ,   , r ) is available in the case of correlated Brownian motions  V ,   ,   .Brigo and Alfonsi [29], who deal with this issue in a different context, propose to use a simple Gaussian approximation, instead of the exact solution.More recently, Grzelak and Oosterlee [16] proposed more sophisticated approximations in the framework of the Heston/CIR hybrid model.We do not follow this path here, however, and we focus instead on finding a semianalytical solution, since this goal can be achieved under Assumptions (A.1)-(A.6).
Let us now introduce a convenient change of the underlying probability measure, from the domestic spot martingale measure P to the domestic forward martingale measure P  .Definition 6.The domestic forward martingale measure P  , equivalent to P on (Ω, F  ), is defined by the Radon-Nikodým derivative process  = (  ) ∈[0,] , where An application of the Girsanov theorem shows that the process   = (   ) ∈[0,] , which is given by the equality ( The following auxiliary result is easy to establish and thus its proof is omitted.Recall that   (, ) is given by equality (18) This shows that P is a martingale measure associated with the choice of the price process     (, ) as a numéraire asset.
We are now in a position to state the following lemma.

Lemma 9. The price of the FX call option satisfies
or, equivalently, To complete the proof Theorem 3, it remains to evaluate the conditional probabilities arising in formula (43).By another application of the Girsanov theorem, one can check that the process (, V, , r) has the Markov property under the probability measures P  and P .In view of Proposition 1 and Lemma 2, the random variable   is a function of   ,   , and r .We thus conclude that where we denote  1 (,   , V  ,   , r , ) = P (  > ln  |   , V  ,   , r ) ,  2 (,   , V  ,   , r , ) = P  (  > ln  |   , V  ,   , r ) .

(45)
To obtain explicit formulae for the conditional probabilities above, it suffices to derive the corresponding conditional characteristic functions: The idea is to use the Radon-Nikodým derivatives in order to obtain convenient expressions for the characteristic functions in terms of conditional expectations under the domestic spot martingale measure P. The following lemma will allow us to achieve this goal.
Using (32), we now obtain which is the desired expression.
In view of the formula established in Lemma 10 and the abstract Bayes formula, to compute  1 () =  1 (, ,   , V  ,   , r ), it suffices to focus on the following conditional expectation under P: Similarly, in view of formula (31), we obtain for  2 () =  2 (, ,   , V  ,   , r ) To proceed, we will need the following result, which is an immediate consequence of Lemma 7.

Corollary 11. Under Assumptions (A.1)-(A.4), the process
= ln (, ) admits the following representation under the domestic forward martingale measure P  : For the sake of conciseness, we denote  = 1 + ,  = , and   = exp(  ).After simplifications and rearrangement, the formula above becomes In view of Assumptions (A.1)-(A.6),we may use the following representation for the Brownian motion   : where  = (  ) ∈[0,] is a Brownian motion under P independent of the Brownian motions  V ,   , and   .Consequently, the conditional characteristic function  1 () can be represented in the following way: By combining Proposition 1 with Definition 6, we obtain the following auxiliary result, which will be helpful in the proof of Theorem 3. Lemma 12.Given the dynamics (1) of processes V, , and r and formula (32), we obtain the following equalities: Proof.The first asserted formula is an immediate consequence of (1).For the second, we recall that the function   (, ) is known to satisfy the following differential equation, for any fixed  > 0: International Journal of Stochastic Analysis with the terminal condition   (, ) = 0. Therefore, using the Itô formula and equality (32), we obtain This yields the second asserted formula, upon integration between  and .The derivation of the last one is based on the same arguments and thus it is omitted.

Proof of Theorem 3.
We split the proof of Theorem 3 into two steps in which we deal with  1 () and  2 ().
Step 1.We will first compute  1 ().By combining (57) with the equalities derived in Lemma 12, we obtain the following representation for  1 (): Recall the well-known property that if  has the standard normal distribution then E(  ) =   2 /2 for any complex number  ∈ C.
Consequently, by conditioning first on the sample path of the process (V, , r) and using the independence of the processes (V, , r) and  under P and Lemma 4, we obtain where we denote  = 1 − .This in turn implies that the following equality holds: where the constants  1 ,  2 ,  3 ,  4 ,  5 ,  6 are given by (16).A direct application of Lemma 5 furnishes an explicit formula for  1 (), as reported in the statement of Theorem 3.
Step 2. In order to compute the conditional characteristic function (67) Similarly as in the case of  1 (), we condition on the sample path of the process (V, , r) and we use the postulated independence of the processes (V, , r) and  under P.By invoking also Lemma 4, we obtain Using Lemma 12, we conclude that with the coefficients  1 ,  2 ,  3 ,  4 ,  5 ,  6 reported in formula (17).Another straightforward application of Lemma 5 yields the closed-form expression (14) for the conditional characteristic function  2 ().
To complete the proof of Theorem 3, it suffices to combine formula (44) with the standard inversion formula (12) providing integral representations for the conditional probabilities: This ends the derivation of the pricing formula for the foreign exchange call option.The price of the corresponding put option is readily available as well, due to the put-call parity relationship for FX options (see formula (72) in Section 5).

Numerical Results
The goal of the final section is to illustrate our approach by means of numerical examples in which we apply our FX market model, that is, the Heston/CIR jump-diffusion model, and we compare this approach with other related models that were proposed in Moretto et al. [30] and Ahlip and Rutkowski [1] to deal with the exchange rate derivatives.
Let us start by noting that the foreign exchange market differs from equity markets in that quotes for options are not made in terms of strikes.Indeed, the FX option prices are quoted in terms of the associated implied volatilities for a fixed forward delta Δ  and a fixed time to expiry  =  − .
For more information about the market conventions, the interested reader is referred to, for instance, Moretto et al. [30] or Reiswich and Uwe [31].
For a quoted volatility , the corresponding strike price  is obtained using the following conversion formula, which is based on the classic Garman-Kohlhagen lognormal model for the exchange rate: where  −1 is the inverse of the standard normal cumulative distribution function and the auxiliary parameter  satisfies  = 1 ( = −1, resp.) for the call (put, resp.)option.Formula (71) makes it clear that market quotations prices based on the implied volatility for fixed deltas are in fact equivalent to quoting prices for fixed strikes.For more details, the interested reader is referred to Hakala and Wystup [32] or Reiswich and Uwe [31].Another relevant feature is that currency derivatives are based on the notion of at-the-money forward (ATMF) rate, that is, the forward exchange rate (, ) obtained by exploiting the interest rate parity implicit in (8).Recall that the universal put-call parity formula for plain-vanilla foreign exchange options reads where   (, ) and   (, ) are prices of currency call and put options, respectively.In particular, the prices of ATMF call and put options are equal in any arbitrage-free market model.5.1.Market Data.In the numerical results presented in Tables 1, 2, and 3, we make use (with the kind permission of the authors) of the data for the USD/EUR exchange rate derivatives and interest rates from the paper by Moretto et al. [30] (see page 469 therein).

Comparison of Model Prices.
The dynamics of the exchange rate and volatility, as given by (1), involve the parameters:   ,   ,   ,  V , , , and  V .In addition, there are three parameters for each of the interest rates.In our numerical examples, the values of parameters , , and  V are borrowed from Moretto et al. [30], who proposed an extension of the Heston model for the exchange rate under the

Option value
Figure 2: Graphs for 25% options prices given in Table 6.
ATM and 25% and 15% volatilities (refer to Figures 3, 2, and 5), respectively.For each maturity, the corresponding strike level was obtained using (71).Prices for ATM, 25%, and 15% USD/EUR currency call options are computed in the Heston, HCIR, and HCIR-LN-LU models using data of June 13, 2005, and parameter values given in Table 4.  9.
Paul Glasserman, and participants at the IMS-FPS-2014 Workshop, July 2-6, 2014, UTS, Sydney, for valuable comments on an earlier version of the paper.All remaining errors are theirs.

International Journal of Stochastic Analysis 7 is
the Brownian motion under the domestic forward martingale measure P  .Using the standard change of a numéraire technique, one can check that the price of the European foreign exchange call option admits the following representation under the probability measure P  :   (, ) =   (, ) E P   ( (, ) 1 {(,)>} ) −   (, ) E P   (1 {(,)>} ) .
[27]proof of Theorem 3 hinges on a number of lemmas.We start by stating the well-known result, which can be easily obtained from Proposition 8.6.3.4 in Jeanblanc et al.[27].Let us denote  =  −  and let us set, for all 0 ≤  < , International Journal of Stochastic Analysis 5 4.1.Auxiliary Results.

Table 5 :
Prices of ATM USD/EUR European exchange rate call options using data of June 13, 2005.Figure1: Graphs for ATM options prices given in Table5.

Table 7 :
Prices for ATM USD/EUR European exchange rate call options using data of June 13, 2005.Values in brackets are strike levels obtained using formula (71).

Table 8 :
Prices for 25% USD/EUR European exchange rate call options using data of June 13, 2005.Values in brackets are strike levels obtained using formula (71).

Table 9 :
Prices for 15% USD/EUR European exchange rate call options using data of June 13, 2005.Values in brackets are strike levels obtained using formula (71).