# Pricing Spread Options with Stochastic Interest Rates

Authors: Yunguo Jin; Shouming Zhong
Journal: Mathematical Problems in Engineering
 (2014)

---

## Abstract

Although spread options have been extensively studied in the literature, few papers deal with the problem
of pricing spread options with stochastic interest rates. This study presents three novel spread option
pricing models that permit the interest rates to be random. The paper not only presents a good approach
to formulate spread option pricing models with stochastic interest rates but also offers a new test bed to
understand the dynamics of option pricing with interest rates in a variety of asset pricing models. We
discuss the merits of the models and techniques presented by us in some asset pricing models. Finally, we
use regular grid method to the calculation of the formula when underlying stock returns are continuous and
a mixture of both the regular grid method and a Monte Carlo method to the one when underlying stock
returns are discontinuous, and sensitivity analyses are presented.


---

## Body

## 1. Introduction

Spread options are financial instruments whose payoffs depend on the price difference of two underlying assets and are widely traded in a range of financial markets. For example, the price spread between heating oil and crude oil (crack spread) represents the value of production (including profit) for a refinery firm. In particular, to energy trading, the derivative contracts are a risk mitigation tool of crucial importance. Carmona and Durrleman [1] presented a general overview of the common features of all spread options by discussing in detail their roles as speculation devices and risk management tools. Their review concentrated on the mathematics of the pricing and hedging of spread option. Hurd and Zhou [2] showed a numerical integration method for computing spread options in two or higher dimensions using the FFT. Tt is based on square integrable integral formulas for the payoff function and, like those methods, it is applicable to a variety of spread option payoff in any model for which the characteristic function of the joint return process is given analytically. Caldana and Fusai [3] proposed a new accurate method for pricing European spread options and obtained a lower bound, as in Bjerksund and Stensland [4], but for general processes. Lo [5] applied the idea of WKB method to provide a simple derivation of Kirks approximation and discussed its validity.Although spread options [1–9] have come into the limelight in recent years, little research has been published on the pricing of spread options with stochastic interest rates. For options on bonds and other interest-rate-dependent instruments, movements in the interest rate are critical and the assumption of a constant interest rate is inappropriate. In this case, stochastic interest rate may significantly affect the prices.The pricing of options in the presence of stochastic interest rates can generally be difficult. However, as was discussed in Shreve [10], the choice of another numéraire (e.g., a stochastic bond price) may be more convenient for option pricing in a model with stochastic interest rate. As was pointed out in Kwok [11], a model can be complicated or simple, depending on the choice of the numéraire for the model. Merton [12] used a bond maturing on the option expiration date to hedge European options. Jamshidian [13] observed that the forward price of an asset is its price when denominated in the numéraire of the zero-coupon bond maturing at the delivery date. Geman et al. [14] presented a generalized Black-Scholes-Merton option pricing formula that permits the interest rate to be random.However, sometimes it is inappropriate to assume that the volatility of the discounted bond price is a constant. There is no reason why interest rates should behave like stock prices. Interest rates certainly do not exhibit the long-term exponential growth seen in the equity markets.On the other hand, mean-reverting models are used for modeling a random variable that is not going anywhere. That is why they are often used for interest rates. The Vasicek model takes the following form:(1)
d
R
(
t
)
=
(
a
-
b
R
(
t
)
)
d
t
+
σ
d
W
(
t
)
,
where a, b, and σ are positive constants and W
(
t
) is a one-dimensional Brownian motion. The model is mean-reverting to a constant level, which is a good property. It is easy to obtain the solution so that there are explicit formulae for many interest rate derivatives. Thus, it is desirable to have a model that combines the strengths of the existing methods.Furthermore, the technique of measure changes is useful in a variety of asset pricing models. However, few papers focus on the discussion of the merits for it and its applications.Credit risk arises from the possibility that borrowers and counterparties in derivatives transactions may default. Most financial institutes devote considerable resources to the measurement and management of credit risk. As was pointed out in Hull [15], the credit exposure on a derivatives transaction is more complicated than that on a loan. This is because the claim that will be made in the event of a default is more uncertain. Klein [16] presented an improved method of pricing vulnerable Black-Scholes options under assumptions which were appropriate in many business situations. An analytic pricing formula is derived which allowed not only for correlation between the option’s underlying asset and the credit risk of the counterparty but also for the option writer to have other liabilities. Yang et al. [17] extended the model of Klein [16] and considered the pricing of vulnerable options when the underlying asset follows a stochastic volatility model. They used multiscale asymptotic analysis to derive an analytic approximation formula for the price of the vulnerable options and studied the stochastic volatility effect on the option price.This study aims to fill the gaps and presents three novel spread options pricing models that permit the interest rates to be random. We assume that a discounted zero-coupon bond follows a stochastic process. The volatility of the discounted bond price is a function of timet rather than a constant. The study presents the technique of measure changes to the computations of the spread options. We discuss the excellent feature for the models and techniques presented by us in some asset pricing.The organization of the remaining part is as follows. Section2 outlines a spread option pricing model with stochastic interest rates and presents a technique of measure changes to the price of the spread option. Section 3 extends this analysis to jump diffusion model. Section 4 extends this analysis to a vulnerable spread option. Section 5 presents numerical analysis. The final section gives concluding remarks.


## 2. Spread Option with Stochastic Interest Rates

This section constructs a spread option pricing model that permits the interest rate to be random. A key point in the model is to assume that the volatility of the discounted bond price is a function of timet rather than a constant. Using the technique of measure changes, we analyze the price of the spread option when underlying stock returns are continuous. We discuss the excellent features for the model and the technique in some asset pricing models.


### 2.1. Problem Formulation

Let(
Ω
,
F
,
P
~
) be a probability space on which is defined a three-dimensional Brownian motion W
~
(
t
)
=
(
W
~
a
(
t
)
,
W
~
b
(
t
)
,
W
~
c
(
t
)
), 0
≤
t
≤
T, and let P
~ denote the risk neutral probability measure. In particular, W
~
a
(
t
), W
~
b
(
t
), and W
~
c
(
t
) are independent Brownian motions. Let F
(
t
), 0
≤
t
≤
T, be the filtration generated by the Brownian motion. We assume that the interest rate process R
(
t
) is adapted and define a discount process D
(
t
)
=
e
-
∫
0
t
‍
R
(
u
)
d
u. Consider a zero-coupon bond that pays 1 unit of currency at maturity T. According to the risk-neutral pricing formula, the value of this bond at time t
∈
[
0
,
T
] is B
(
t
,
T
)
=
1
/
D
(
t
)
E
~
D
(
T
)
∣
F
(
t
). Let S
1
(
t
) and S
2
(
t
) be two stock price processes. The evolutions of the discounted prices will be modeled by the following stochastic differential equations:
(2)
d
(
D
(
t
)
S
i
(
t
)
)
=
D
(
t
)
S
i
(
t
)
v
i
·
d
W
~
i
(
t
)
,
i
=
1,2
,
(3)
d
D
(
t
)
B
(
t
,
T
)
=
D
(
t
)
B
(
t
,
T
)
v
3
(
t
,
T
)
·
d
W
~
3
(
t
)
,
where W
~
1
(
t
)
=
(
W
~
a
(
t
)
,
0,0
), W
~
2
(
t
)
=
(
W
~
a
(
t
)
,
W
~
b
(
t
)
,
0
), and W
~
3
(
t
)
=
(
W
~
a
(
t
)
,
W
~
b
(
t
)
,
W
~
c
(
t
)
). v
1
=
(
v
11
,
0,0
), v
2
=
(
v
21
,
v
22
,
0
), and v
3
(
t
,
T
)
=
(
v
31
(
t
,
T
)
,
v
32
(
t
,
T
)
,
v
33
(
t
,
T
)
) are three vector volatility processes. The value at time t
∈
[
0
,
T
] of a spread option, expiring at time T with strike price K, is
(4)
V
(
t
)
=
E
~
e
-
∫
t
T
‍
R
(
u
)
d
u
S
1
T
-
S
2
T
-
K
+
∣
F
(
t
)
.


### 2.2. Pricing the Spread Option

Definition 1.
One defines three new measures.(1)
Suppose one takes the asset priceS
i
(
t
), i
=
1,2, to be the numéraire. One defines the risk-neutral measure for this numéraire by
(5)
P
~
S
i
(
A
)
=
1
S
i
(
0
)
∫
A
‍
D
(
T
)
S
i
(
T
)
d
P
~
,
for all A
∈
F
(
T
).
(2)
LetT be a fixed maturity date. One defines the T-forward measure P
~
T by
(6)
P
~
T
(
A
)
=
1
B
(
0
,
T
)
∫
A
‍
D
(
T
)
d
P
~
,
for all A
∈
F
(
T
).Let(
Ω
2
,
F
2
,
P
2
) be a probability space on which is defined a multidimensional Brownian motion W
(
t
)
=
(
W
1
(
t
)
,
W
2
(
t
)
,
…
,
W
d
(
t
)
), 0
≤
t
≤
T. Let F
2
(
t
), 0
≤
t
≤
T, be the filtration generated by the Brownian motion. The proof of Lemma 2 is like that of the one-dimensional Girsanov Theorem  5.2
.
3 in [10], except it uses a d-dimensional version of Lévy’s Theorem.Lemma 2 (Girsanov, multiple dimensions, [10]).
LetΘ
(
t
)
=
Θ
1
(
t
)
,
Θ
2
(
t
)
,
…
,
Θ
d
(
t
) be a d-dimensional adapted process. Define
(7)
Z
(
t
)
=
exp
⁡
-
∫
0
t
‍
Θ
(
u
)
·
d
W
(
u
)
-
1
2
∫
0
t
‍
Θ
(
u
)
2
d
u
,
W
¯
(
t
)
=
W
(
t
)
+
∫
0
t
‍
Θ
(
u
)
d
u
,
and assume that
(8)
E
∫
0
t
‍
Θ
(
u
)
2
Z
2
(
u
)
d
u
<
∞
.
Then, under the probability measure P
¯ given by
(9)
P
¯
(
A
)
=
∫
A
‍
Z
(
ω
)
d
P
(
ω
)
,
∀
A
∈
F
2
(
T
)
,
the process W
¯
(
t
) is a d-dimensional Brownian motion.For the following Lemma3, we refer the reader to Theorem  4.3
.
1 in [10] for its discussion.Lemma 3 (see [10]).
LetW
0
(
t
) be a Brownian motion and let Δ
(
t
) be an adapted stochastic process that satisfies E
∫
0
T
‍
Δ
2
(
t
)
d
t
<
∞. Then, ∫
0
t
‍
Δ
(
u
)
d
W
0
(
u
) is a martingale.Proposition 4.
Under a new measure, the following price dynamics holds.(1)
Under theT-forward measure P
~
T, the process S
i
(
t
)
/
B
(
t
,
T
), i
=
1,2, is a martingale. Moreover,
(10)
d
S
i
(
t
)
B
(
t
,
T
)
=
S
i
(
t
)
B
(
t
,
T
)
(
v
i
-
v
3
(
t
,
T
)
)
·
d
W
~
T
(
t
)
,
where W
~
T
(
t
)
=
(
W
~
a
(
t
)
,
W
~
b
(
t
)
,
W
~
c
(
t
)
)
-
∫
0
t
‍
v
3
(
u
,
T
)
d
u is a three-dimensional Brownian motion under P
~
T.
(2)
Under the measuresP
~
S
i, i
=
1,2, the processes B
(
t
,
T
)
/
S
i
(
t
) and S
j
(
t
)
/
S
i
(
t
), j
=
1,2, i
≠
j, are a martingale. Moreover,
(11)
d
B
(
t
,
T
)
S
i
(
t
)
=
B
(
t
,
T
)
S
i
(
t
)
(
v
3
(
t
,
T
)
-
v
i
)
·
d
W
~
S
i
(
t
)
,
d
S
j
(
t
)
S
i
(
t
)
=
S
j
(
t
)
S
i
(
t
)
(
v
j
-
v
i
)
·
d
W
~
S
i
(
t
)
,
where W
~
S
i
(
t
)
=
(
W
~
a
(
t
)
,
W
~
b
(
t
)
,
W
~
c
(
t
)
)
-
v
i
t is a three-dimensional Brownian motion under P
~
S
i.The proof of Proposition4 will be given in Appendix A.We need some notations and denote by· the Euclidean norm in R
3. Let
(12)
A
i
3
(
t
)
=
S
i
(
t
)
B
(
t
,
T
)
exp
⁡
-
1
2
∫
t
T
‍
(
v
i
-
v
3
(
u
,
T
)
)
2
d
u
,
B
i
3
(
t
)
=
∫
t
T
‍
v
i
-
v
3
(
u
,
T
)
2
d
u
,
A
3
i
(
t
)
=
B
(
t
,
T
)
S
i
(
t
)
exp
⁡
-
1
2
∫
t
T
‍
(
v
3
(
u
,
T
)
-
v
i
)
2
d
u
,
B
3
i
(
t
)
=
∫
t
T
‍
v
3
(
u
,
T
)
-
v
i
2
d
u
,
A
i
j
t
=
S
i
(
t
)
S
j
(
t
)
exp
⁡
-
1
2
∫
t
T
‍
v
i
-
v
j
2
d
u
,
B
i
j
t
=
v
i
-
v
j
T
-
t
,
h
h
h
h
h
h
h
h
h
h
h
h
h
i
h
h
h
h
h
h
h
i
,
j
=
1,2
,
i
≠
j
,
ρ
1
=
(
v
2
-
v
1
)
·
∫
t
T
‍
v
3
u
,
T
-
v
1
2
d
u
v
2
-
v
1
∫
t
T
‍
v
3
u
,
T
-
v
1
2
d
u
,
ρ
2
=
(
v
1
-
v
2
)
·
∫
t
T
‍
v
3
u
,
T
-
v
2
2
d
u
v
1
-
v
2
∫
t
T
‍
v
3
u
,
T
-
v
2
2
d
u
,
ρ
3
=
∫
t
T
‍
v
1
-
v
3
u
,
T
2
d
u
·
∫
t
T
‍
v
2
-
v
3
u
,
T
2
d
u
∫
t
T
‍
v
1
-
v
3
u
,
T
2
d
u
∫
t
T
‍
v
2
-
v
3
u
,
T
2
d
u
.For spread options, there is a long history of approximation methods for the computation of the prices. The challenge in pricing spread options stems from the fact that there is no explicit solution. Furthermore, the pricing of options with stochastic interest rates can generally be difficult. Therefore, pricing spread options with stochastic interest rates is a challenging task in finance and is important for both theory and applications.Theorem 5.
The price of the spread option, assuming the risk-neutral dynamics given in (2) and (3), is V
(
t
)
=
V
1
(
t
)
+
V
2
(
t
)
+
V
3
(
t
) with
(13)
V
1
t
=
S
1
t
×
Pr
⁡
1
-
A
21
t
e
B
21
t
ξ
11
-
K
A
31
t
e
B
31
(
t
)
ξ
12
>
0
,
V
2
t
=
-
S
2
(
t
)
×
Pr
⁡
A
12
t
e
B
12
t
ξ
21
-
1
-
K
A
32
t
e
B
32
(
t
)
ξ
22
>
0
,
V
3
t
=
-
K
B
(
t
,
T
)
×
Pr
⁡
A
13
t
e
B
13
(
t
)
ξ
31
-
A
23
(
t
)
e
B
23
(
t
)
ξ
32
-
K
>
0
,
where the random variables ξ
i
1 and ξ
i
2, i
=
1,2
,
3, have the bivariate standard normal distribution with correlation ρ
i. For event Λ, a number Pr
⁡
(
Λ
) indicates the probability that Λ will occur.The proof of Theorem5 will be given in Appendix B. The delta of an option or a portfolio of options is the sensitivity of the option or portfolio to the underlying. The sensitivities of the spread option price to initial stock prices S
i, i
=
1,2, and spread K are given by Δ
1
=
∂
V
/
∂
S
1, Δ
2
=
∂
V
/
∂
S
2, and κ
=
∂
V
/
∂
K. The simplest way to calculate the delta, the gamma, and the theta of an option using Monte Carlo simulation is to estimate the option value twice.Remark 6.
According to our knowledge, the paper is the first one investigating spread option pricing with stochastic interest rates. For stochastic interest rates, the study presents a more generalized spread option pricing model. Some models of stochastic interest rates are just a special case of the model presented by us. For Vasicek model (1), interest rates can easily go negative, which is a very bad property. Geman et al. [14] assumed that the volatility of the discounted bond price is a constant. For a random variable that is not going anywhere, sometimes the assumption is inappropriate. We assume that the volatility of the discounted bond price is a function of time t rather than a constant.Remark 7.
The spread option pricing formula can be used in the standard call option. We present also a generalized Black-Scholes-Merton option pricing formula with stochastic interest rates. In fact, the price of the European call, assuming the risk-neutral dynamics given in (2) and (3) (i
=
1), is V
(
t
)
=
V
1
(
t
)
+
V
3
(
t
), with
(14)
V
1
(
t
)
=
S
1
(
t
)
Pr
⁡
1
-
K
A
31
(
t
)
e
B
31
(
t
)
ξ
1
>
0
,
V
3
(
t
)
=
-
K
B
(
t
,
T
)
Pr
⁡
A
13
(
t
)
e
B
13
(
t
)
ξ
3
-
K
>
0
,
where the random variables ξ
1 and ξ
3 have the standard normal distribution. If the interest rate is a constant r, then B
(
t
,
T
)
=
e
-
r
(
T
-
t
); the result reduces to the usual Black-Shcholes-Merton formula. The price of an exchange option (spread options with K
=
0), assuming the risk-neutral dynamics given in (2) and (3), is V
(
t
)
=
V
1
(
t
)
+
V
2
(
t
)
+
V
3
(
t
) with
(15)
V
1
(
t
)
=
S
1
(
t
)
Pr
⁡
1
-
A
21
(
t
)
e
B
21
(
t
)
ξ
1
>
0
,
(16)
V
2
(
t
)
=
-
S
2
(
t
)
Pr
⁡
A
12
(
t
)
e
B
12
(
t
)
ξ
2
-
1
>
0
,
(17)
V
3
t
=
-
K
B
(
t
,
T
)
×
Pr
⁡
A
13
(
t
)
e
B
13
(
t
)
ξ
31
-
A
23
(
t
)
e
B
23
(
t
)
ξ
32
>
0
,
where the random variables ξ
1 and ξ
2 have the standard normal distribution and the random variables ξ
31 and ξ
32 have the bivariate standard normal distribution with correlation ρ
3.Remark 8.
The technique of measure changes is useful in a variety of asset pricing models. The risk-neutral valuation of the spread option price involves a two-dimensional integration rather than a three-dimensional integration. Monte Carlo approximations always contain some randomness. If we use some methods rather than Monte Carlo approximations, for the two-dimensional integration, we can increase the accuracy of the price or improve the computational efficiency. In particular, for the generalized Black-Scholes-Merton option pricing formula with stochastic interest rates and the exchange option pricing formula in Remark7, it is of great significance. In the cases, the risk-neutral valuations involve a one-dimensional integration rather than a two-dimensional integration. Furthermore, by the technique of measure changes, the integrands become simple. As was pointed out in Carmona and Durrleman [1], even the most sophisticated models are very crude approximations of reality, and traders prefer having simple and reliable closed form formulae whose behavior can be well understood.


## 2.1. Problem Formulation

Let(
Ω
,
F
,
P
~
) be a probability space on which is defined a three-dimensional Brownian motion W
~
(
t
)
=
(
W
~
a
(
t
)
,
W
~
b
(
t
)
,
W
~
c
(
t
)
), 0
≤
t
≤
T, and let P
~ denote the risk neutral probability measure. In particular, W
~
a
(
t
), W
~
b
(
t
), and W
~
c
(
t
) are independent Brownian motions. Let F
(
t
), 0
≤
t
≤
T, be the filtration generated by the Brownian motion. We assume that the interest rate process R
(
t
) is adapted and define a discount process D
(
t
)
=
e
-
∫
0
t
‍
R
(
u
)
d
u. Consider a zero-coupon bond that pays 1 unit of currency at maturity T. According to the risk-neutral pricing formula, the value of this bond at time t
∈
[
0
,
T
] is B
(
t
,
T
)
=
1
/
D
(
t
)
E
~
D
(
T
)
∣
F
(
t
). Let S
1
(
t
) and S
2
(
t
) be two stock price processes. The evolutions of the discounted prices will be modeled by the following stochastic differential equations:
(2)
d
(
D
(
t
)
S
i
(
t
)
)
=
D
(
t
)
S
i
(
t
)
v
i
·
d
W
~
i
(
t
)
,
i
=
1,2
,
(3)
d
D
(
t
)
B
(
t
,
T
)
=
D
(
t
)
B
(
t
,
T
)
v
3
(
t
,
T
)
·
d
W
~
3
(
t
)
,
where W
~
1
(
t
)
=
(
W
~
a
(
t
)
,
0,0
), W
~
2
(
t
)
=
(
W
~
a
(
t
)
,
W
~
b
(
t
)
,
0
), and W
~
3
(
t
)
=
(
W
~
a
(
t
)
,
W
~
b
(
t
)
,
W
~
c
(
t
)
). v
1
=
(
v
11
,
0,0
), v
2
=
(
v
21
,
v
22
,
0
), and v
3
(
t
,
T
)
=
(
v
31
(
t
,
T
)
,
v
32
(
t
,
T
)
,
v
33
(
t
,
T
)
) are three vector volatility processes. The value at time t
∈
[
0
,
T
] of a spread option, expiring at time T with strike price K, is
(4)
V
(
t
)
=
E
~
e
-
∫
t
T
‍
R
(
u
)
d
u
S
1
T
-
S
2
T
-
K
+
∣
F
(
t
)
.


## 2.2. Pricing the Spread Option

Definition 1.
One defines three new measures.(1)
Suppose one takes the asset priceS
i
(
t
), i
=
1,2, to be the numéraire. One defines the risk-neutral measure for this numéraire by
(5)
P
~
S
i
(
A
)
=
1
S
i
(
0
)
∫
A
‍
D
(
T
)
S
i
(
T
)
d
P
~
,
for all A
∈
F
(
T
).
(2)
LetT be a fixed maturity date. One defines the T-forward measure P
~
T by
(6)
P
~
T
(
A
)
=
1
B
(
0
,
T
)
∫
A
‍
D
(
T
)
d
P
~
,
for all A
∈
F
(
T
).Let(
Ω
2
,
F
2
,
P
2
) be a probability space on which is defined a multidimensional Brownian motion W
(
t
)
=
(
W
1
(
t
)
,
W
2
(
t
)
,
…
,
W
d
(
t
)
), 0
≤
t
≤
T. Let F
2
(
t
), 0
≤
t
≤
T, be the filtration generated by the Brownian motion. The proof of Lemma 2 is like that of the one-dimensional Girsanov Theorem  5.2
.
3 in [10], except it uses a d-dimensional version of Lévy’s Theorem.Lemma 2 (Girsanov, multiple dimensions, [10]).
LetΘ
(
t
)
=
Θ
1
(
t
)
,
Θ
2
(
t
)
,
…
,
Θ
d
(
t
) be a d-dimensional adapted process. Define
(7)
Z
(
t
)
=
exp
⁡
-
∫
0
t
‍
Θ
(
u
)
·
d
W
(
u
)
-
1
2
∫
0
t
‍
Θ
(
u
)
2
d
u
,
W
¯
(
t
)
=
W
(
t
)
+
∫
0
t
‍
Θ
(
u
)
d
u
,
and assume that
(8)
E
∫
0
t
‍
Θ
(
u
)
2
Z
2
(
u
)
d
u
<
∞
.
Then, under the probability measure P
¯ given by
(9)
P
¯
(
A
)
=
∫
A
‍
Z
(
ω
)
d
P
(
ω
)
,
∀
A
∈
F
2
(
T
)
,
the process W
¯
(
t
) is a d-dimensional Brownian motion.For the following Lemma3, we refer the reader to Theorem  4.3
.
1 in [10] for its discussion.Lemma 3 (see [10]).
LetW
0
(
t
) be a Brownian motion and let Δ
(
t
) be an adapted stochastic process that satisfies E
∫
0
T
‍
Δ
2
(
t
)
d
t
<
∞. Then, ∫
0
t
‍
Δ
(
u
)
d
W
0
(
u
) is a martingale.Proposition 4.
Under a new measure, the following price dynamics holds.(1)
Under theT-forward measure P
~
T, the process S
i
(
t
)
/
B
(
t
,
T
), i
=
1,2, is a martingale. Moreover,
(10)
d
S
i
(
t
)
B
(
t
,
T
)
=
S
i
(
t
)
B
(
t
,
T
)
(
v
i
-
v
3
(
t
,
T
)
)
·
d
W
~
T
(
t
)
,
where W
~
T
(
t
)
=
(
W
~
a
(
t
)
,
W
~
b
(
t
)
,
W
~
c
(
t
)
)
-
∫
0
t
‍
v
3
(
u
,
T
)
d
u is a three-dimensional Brownian motion under P
~
T.
(2)
Under the measuresP
~
S
i, i
=
1,2, the processes B
(
t
,
T
)
/
S
i
(
t
) and S
j
(
t
)
/
S
i
(
t
), j
=
1,2, i
≠
j, are a martingale. Moreover,
(11)
d
B
(
t
,
T
)
S
i
(
t
)
=
B
(
t
,
T
)
S
i
(
t
)
(
v
3
(
t
,
T
)
-
v
i
)
·
d
W
~
S
i
(
t
)
,
d
S
j
(
t
)
S
i
(
t
)
=
S
j
(
t
)
S
i
(
t
)
(
v
j
-
v
i
)
·
d
W
~
S
i
(
t
)
,
where W
~
S
i
(
t
)
=
(
W
~
a
(
t
)
,
W
~
b
(
t
)
,
W
~
c
(
t
)
)
-
v
i
t is a three-dimensional Brownian motion under P
~
S
i.The proof of Proposition4 will be given in Appendix A.We need some notations and denote by· the Euclidean norm in R
3. Let
(12)
A
i
3
(
t
)
=
S
i
(
t
)
B
(
t
,
T
)
exp
⁡
-
1
2
∫
t
T
‍
(
v
i
-
v
3
(
u
,
T
)
)
2
d
u
,
B
i
3
(
t
)
=
∫
t
T
‍
v
i
-
v
3
(
u
,
T
)
2
d
u
,
A
3
i
(
t
)
=
B
(
t
,
T
)
S
i
(
t
)
exp
⁡
-
1
2
∫
t
T
‍
(
v
3
(
u
,
T
)
-
v
i
)
2
d
u
,
B
3
i
(
t
)
=
∫
t
T
‍
v
3
(
u
,
T
)
-
v
i
2
d
u
,
A
i
j
t
=
S
i
(
t
)
S
j
(
t
)
exp
⁡
-
1
2
∫
t
T
‍
v
i
-
v
j
2
d
u
,
B
i
j
t
=
v
i
-
v
j
T
-
t
,
h
h
h
h
h
h
h
h
h
h
h
h
h
i
h
h
h
h
h
h
h
i
,
j
=
1,2
,
i
≠
j
,
ρ
1
=
(
v
2
-
v
1
)
·
∫
t
T
‍
v
3
u
,
T
-
v
1
2
d
u
v
2
-
v
1
∫
t
T
‍
v
3
u
,
T
-
v
1
2
d
u
,
ρ
2
=
(
v
1
-
v
2
)
·
∫
t
T
‍
v
3
u
,
T
-
v
2
2
d
u
v
1
-
v
2
∫
t
T
‍
v
3
u
,
T
-
v
2
2
d
u
,
ρ
3
=
∫
t
T
‍
v
1
-
v
3
u
,
T
2
d
u
·
∫
t
T
‍
v
2
-
v
3
u
,
T
2
d
u
∫
t
T
‍
v
1
-
v
3
u
,
T
2
d
u
∫
t
T
‍
v
2
-
v
3
u
,
T
2
d
u
.For spread options, there is a long history of approximation methods for the computation of the prices. The challenge in pricing spread options stems from the fact that there is no explicit solution. Furthermore, the pricing of options with stochastic interest rates can generally be difficult. Therefore, pricing spread options with stochastic interest rates is a challenging task in finance and is important for both theory and applications.Theorem 5.
The price of the spread option, assuming the risk-neutral dynamics given in (2) and (3), is V
(
t
)
=
V
1
(
t
)
+
V
2
(
t
)
+
V
3
(
t
) with
(13)
V
1
t
=
S
1
t
×
Pr
⁡
1
-
A
21
t
e
B
21
t
ξ
11
-
K
A
31
t
e
B
31
(
t
)
ξ
12
>
0
,
V
2
t
=
-
S
2
(
t
)
×
Pr
⁡
A
12
t
e
B
12
t
ξ
21
-
1
-
K
A
32
t
e
B
32
(
t
)
ξ
22
>
0
,
V
3
t
=
-
K
B
(
t
,
T
)
×
Pr
⁡
A
13
t
e
B
13
(
t
)
ξ
31
-
A
23
(
t
)
e
B
23
(
t
)
ξ
32
-
K
>
0
,
where the random variables ξ
i
1 and ξ
i
2, i
=
1,2
,
3, have the bivariate standard normal distribution with correlation ρ
i. For event Λ, a number Pr
⁡
(
Λ
) indicates the probability that Λ will occur.The proof of Theorem5 will be given in Appendix B. The delta of an option or a portfolio of options is the sensitivity of the option or portfolio to the underlying. The sensitivities of the spread option price to initial stock prices S
i, i
=
1,2, and spread K are given by Δ
1
=
∂
V
/
∂
S
1, Δ
2
=
∂
V
/
∂
S
2, and κ
=
∂
V
/
∂
K. The simplest way to calculate the delta, the gamma, and the theta of an option using Monte Carlo simulation is to estimate the option value twice.Remark 6.
According to our knowledge, the paper is the first one investigating spread option pricing with stochastic interest rates. For stochastic interest rates, the study presents a more generalized spread option pricing model. Some models of stochastic interest rates are just a special case of the model presented by us. For Vasicek model (1), interest rates can easily go negative, which is a very bad property. Geman et al. [14] assumed that the volatility of the discounted bond price is a constant. For a random variable that is not going anywhere, sometimes the assumption is inappropriate. We assume that the volatility of the discounted bond price is a function of time t rather than a constant.Remark 7.
The spread option pricing formula can be used in the standard call option. We present also a generalized Black-Scholes-Merton option pricing formula with stochastic interest rates. In fact, the price of the European call, assuming the risk-neutral dynamics given in (2) and (3) (i
=
1), is V
(
t
)
=
V
1
(
t
)
+
V
3
(
t
), with
(14)
V
1
(
t
)
=
S
1
(
t
)
Pr
⁡
1
-
K
A
31
(
t
)
e
B
31
(
t
)
ξ
1
>
0
,
V
3
(
t
)
=
-
K
B
(
t
,
T
)
Pr
⁡
A
13
(
t
)
e
B
13
(
t
)
ξ
3
-
K
>
0
,
where the random variables ξ
1 and ξ
3 have the standard normal distribution. If the interest rate is a constant r, then B
(
t
,
T
)
=
e
-
r
(
T
-
t
); the result reduces to the usual Black-Shcholes-Merton formula. The price of an exchange option (spread options with K
=
0), assuming the risk-neutral dynamics given in (2) and (3), is V
(
t
)
=
V
1
(
t
)
+
V
2
(
t
)
+
V
3
(
t
) with
(15)
V
1
(
t
)
=
S
1
(
t
)
Pr
⁡
1
-
A
21
(
t
)
e
B
21
(
t
)
ξ
1
>
0
,
(16)
V
2
(
t
)
=
-
S
2
(
t
)
Pr
⁡
A
12
(
t
)
e
B
12
(
t
)
ξ
2
-
1
>
0
,
(17)
V
3
t
=
-
K
B
(
t
,
T
)
×
Pr
⁡
A
13
(
t
)
e
B
13
(
t
)
ξ
31
-
A
23
(
t
)
e
B
23
(
t
)
ξ
32
>
0
,
where the random variables ξ
1 and ξ
2 have the standard normal distribution and the random variables ξ
31 and ξ
32 have the bivariate standard normal distribution with correlation ρ
3.Remark 8.
The technique of measure changes is useful in a variety of asset pricing models. The risk-neutral valuation of the spread option price involves a two-dimensional integration rather than a three-dimensional integration. Monte Carlo approximations always contain some randomness. If we use some methods rather than Monte Carlo approximations, for the two-dimensional integration, we can increase the accuracy of the price or improve the computational efficiency. In particular, for the generalized Black-Scholes-Merton option pricing formula with stochastic interest rates and the exchange option pricing formula in Remark7, it is of great significance. In the cases, the risk-neutral valuations involve a one-dimensional integration rather than a two-dimensional integration. Furthermore, by the technique of measure changes, the integrands become simple. As was pointed out in Carmona and Durrleman [1], even the most sophisticated models are very crude approximations of reality, and traders prefer having simple and reliable closed form formulae whose behavior can be well understood.


## 3. Extension to Jump Diffusion Model

This section constructs a spread option pricing model when underlying stock returns are discontinuous. As was pointed out in Merton (Ch. 9 in [18]), the abnormal vibrations in price are due to the arrival of important new information about the stock that has more than a marginal effect on price. This component is modeled by a jump process reflecting the nonmarginal impact of the information. In particular, a Poisson process is used to describe the arrival rate of catastrophe events and is applied for the pricing of CAT insurance products.


### 3.1. The Model

Let(
Ω
,
F
,
P
~
) be a probability space on which are defined a three-dimensional Brownian motion W
~
(
t
)
=
(
W
~
1
(
t
)
,
W
~
2
(
t
)
,
W
~
3
(
t
)
), 0
≤
t
≤
T, and two Poisson processes N
1
(
t
) and N
2
(
t
), 0
≤
t
≤
T, and P
~ denote the risk neutral probability measure. In particular, W
~
a
(
t
), W
~
b
(
t
), and W
~
c
(
t
) are independent Brownian motions. Let F
(
t
), 0
≤
t
≤
T, be the filtration generated by the Brownian motion and the Poisson processes. Let λ
1 and λ
2 be the intensities of the Poisson processes N
1
(
t
) and N
2
(
t
), respectively. We define two compound Poisson processes as follows:
(18)
Q
i
t
=
∑
j
=
1
N
i
t
‍
Y
i
j
,
-
1
<
X
i
j
,
i
=
1,2
,
where Y
i
1
,
Y
i
2
,
…, random variables take values in the set {
y
i
1
,
y
i
2
,
…
,
y
i
M
i
2
}. The random variables Y
i
1
,
Y
i
2
,
…, are independent and identically distributed, with P
{
Y
i
j
=
y
i
m
i
}
=
p
(
y
i
m
i
), m
i
=
1,2
,
…
,
M
i. We assume E
Y
i
j
=
∑
m
i
=
1
M
i
‍
y
i
m
i
p
(
y
i
m
i
)
=
β
i. The evolutions of the discounted stock prices will be modeled by the stochastic differential equations:
(19)
d
D
t
S
i
t
=
D
(
t
)
S
i
(
t
)
v
i
·
d
W
~
i
(
t
)
+
D
(
t
)
S
i
(
t
-
)
d
(
Q
i
(
t
)
-
β
i
λ
i
t
)
,
i
=
1,2
.


### 3.2. Pricing the Spread Option

Lemma 9.
Under the measuresP
~
S
i, i
=
1,2, the process W
~
S
i
(
t
)
=
(
W
~
a
(
t
)
,
W
~
b
(
t
)
,
W
~
c
(
t
)
)
-
v
i
t is a three-dimensional Brownian motion. Q
i
(
t
) is a compound Poisson process with intensity (
1
+
β
i
)
λ
i and independent, identically distributed jump sizes satisfying
(20)
P
~
S
i
Y
i
k
=
y
i
m
i
=
1
+
y
i
m
i
p
y
i
m
i
1
+
β
i
,
k
=
1,2
,
…
,
m
i
=
1,2
,
…
,
M
i
,
and the processes W
~
S
i
(
t
) and Q
i
(
t
) are independent.The proof of Lemma9 will be given in Appendix C.Proposition 10.
Under the following measures, the following price dynamics holds.(1)
Under theT-forward measure P
~
T, the process S
i
(
t
)
/
B
(
t
,
T
), i
=
1,2, is a martingale. Moreover,
(21)
d
S
i
(
t
)
B
(
t
,
T
)
=
S
i
t
-
B
t
,
T
v
i
-
v
3
t
,
T
·
d
W
~
T
t
+
d
Q
i
t
-
β
i
λ
i
t
,
where W
~
T
(
t
)
=
(
W
~
a
(
t
)
,
W
~
b
(
t
)
,
W
~
c
(
t
)
)
-
∫
0
t
‍
v
3
(
u
,
T
)
d
u is a three-dimensional Brownian motion under P
~
T. The jump measure is unaffected by the change of measure.
(2)
Under the measuresP
~
S
i, i
=
1,2, the processes B
(
t
,
T
)
/
S
i
(
t
) and S
j
(
t
)
/
S
i
(
t
), j
=
1,2, i
≠
j, are a martingale. Moreover,
(22)
d
B
(
t
,
T
)
S
i
(
t
)
=
B
(
t
,
T
)
S
i
(
t
-
)
∏
k
=
1
N
i
(
t
)
‍
Y
i
k
1
+
Y
i
k
-
β
i
λ
i
t
(
v
3
(
t
,
T
)
-
v
i
)
·
d
W
~
S
i
(
t
)
hhhhhhh
-
d
∏
k
=
1
N
i
(
t
)
‍
Y
i
k
1
+
Y
i
k
-
β
i
λ
i
t
,
d
S
j
t
S
i
t
=
S
j
t
-
S
i
t
-
∏
k
=
1
N
i
t
‍
Y
i
k
1
+
Y
i
k
-
β
i
λ
i
t
v
j
-
v
i
·
d
W
~
S
i
t
hhhhhhh
+
d
(
Q
j
(
t
)
-
β
j
λ
j
t
)
hhhhhhh
-
d
∏
k
=
1
N
i
(
t
)
‍
Y
i
k
1
+
Y
i
k
-
β
i
λ
i
t
,
where W
~
S
i
(
t
)
=
(
W
~
a
(
t
)
,
W
~
b
(
t
)
,
W
~
c
(
t
)
)
-
v
i
t is a three-dimensional Brownian motion under P
~
S
i.The proof of Proposition10 will be given in Appendix D.By Proposition10 and a similar discussion as the one in Theorem 5, we may obtain the following result.Theorem 11.
The price of the spread option, assuming the risk-neutral dynamics given in (3) and (19), is
(23)
V
t
=
∑
i
=
0
∞
‍
∑
j
=
0
∞
∑
L
=
1
2
‍
P
~
S
L
N
1
(
T
)
-
N
1
t
=
i
,
‍
hhhhhhhhhhhh
N
2
T
-
N
2
(
t
)
=
j
S
L
(
t
)
V
L
(
t
,
i
,
j
)
hhhhhhh
+
P
~
T
N
1
T
-
N
1
t
=
i
,
hhhhhhhhhhhh
N
2
(
T
)
-
N
2
t
=
j
hhhhhhh
×
K
B
t
,
T
V
3
t
,
i
,
j
∑
L
=
1
2
‍
with
(24)
V
1
t
,
i
,
j
=
P
~
S
1
1
-
∏
k
=
1
i
‍
∏
L
=
1
j
‍
1
+
Y
2
L
1
+
Y
1
k
A
21
(
t
)
e
B
21
(
t
)
ξ
11
hhhhh
-
∏
k
=
1
i
‍
1
1
+
Y
1
k
K
A
31
(
t
)
e
B
31
(
t
)
ξ
12
>
0
,
V
2
t
,
i
,
j
=
-
P
~
S
2
∐
k
=
1
i
‍
∐
L
=
1
j
‍
1
+
Y
1
k
1
+
Y
2
L
A
12
(
t
)
e
B
12
(
t
)
ξ
21
-
1
hhhhhh
-
∏
L
=
1
j
‍
1
1
+
Y
2
L
K
A
32
(
t
)
e
B
32
(
t
)
ξ
22
>
0
,
V
3
t
,
i
,
j
=
-
P
~
T
∏
k
=
1
i
‍
1
+
Y
1
k
A
13
t
e
B
13
t
ξ
31
hhhhhh
-
∏
L
=
1
j
‍
(
1
+
Y
2
L
)
A
23
hhhhihhhh
×
(
t
)
e
B
23
(
t
)
ξ
31
-
K
>
0
∏
L
=
1
j
‍
,
where
(25)
A
i
3
t
=
S
i
(
t
)
B
(
t
,
T
)
exp
⁡
-
1
2
∫
t
T
‍
(
v
i
-
v
3
(
u
,
T
)
)
2
d
u
h
h
h
h
h
h
h
h
h
h
-
β
i
λ
i
(
T
-
t
)
1
2
∫
t
T
‍
,
A
3
i
t
=
B
(
t
,
T
)
S
i
(
t
)
exp
⁡
-
1
2
∫
t
T
‍
(
v
3
(
u
,
T
)
-
v
i
)
2
d
u
hhhhhihhhhi
+
β
i
λ
i
(
T
-
t
)
1
2
∫
t
T
‍
,
A
i
j
t
=
S
i
(
t
)
S
j
(
t
)
exp
⁡
-
1
2
∫
t
T
‍
(
v
i
-
v
j
)
2
d
u
hhhhhhhii
+
β
i
λ
i
T
-
t
-
β
j
λ
j
T
-
t
1
2
∫
t
T
‍
,
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
i
,
j
=
1,2
,
i
≠
j
.Remark 12.
We can explore further applications of the formula beyond the case that the jump distribution is finitely many jump sizes. Suppose the jump sizesy
i
j, i
=
1,2, j
=
1,2
,
…, have a density f
i
(
y
i
) rather than a probability mass function P
(
y
i
1
)
,
P
(
y
i
2
)
,
…
,
P
(
y
i
M
i
). These densities are strictly positive on a set B
⊂
(
-
1
,
∞
) and zero elsewhere. In the case β
i
=
∫
-
1
∞
‍
f
i
(
x
)
d
x. By Lemma 9 and the definition of multiple integrals, we have that, under the measures P
~
S
i, i
=
1,2, the jumps in Q
i
(
t
) are independent and identically distributed with density (
1
+
y
i
)
f
i
(
y
i
)
/
(
1
+
β
i
).Remark 13.
We can also apply the formula to a catastrophe option or European call in a jump model. The price of the European call, assuming the risk-neutral dynamics given in (3) and (19) (i
=
1), is
(26)
V
t
=
∑
i
=
0
∞
‍
P
~
S
1
N
1
T
-
N
1
t
=
i
S
1
t
V
1
t
,
i
hhhh
+
P
~
T
(
N
1
(
T
)
-
N
1
(
t
)
=
i
)
hhhh
×
K
B
(
t
,
T
)
V
3
(
t
,
i
,
j
)
P
~
S
1
N
1
T
-
N
1
t
=
i
S
1
t
V
1
t
,
i
,
where
(27)
V
1
(
t
,
i
)
=
P
~
S
1
1
-
∏
k
=
1
i
‍
1
1
+
Y
1
k
K
A
31
(
t
)
e
B
31
(
t
)
ξ
12
>
0
,
V
3
(
t
,
i
,
j
)
=
-
P
~
T
∏
k
=
1
i
‍
(
1
+
Y
1
k
)
A
13
(
t
)
e
B
13
(
t
)
ξ
31
-
K
>
0
.Remark 14.
The technique of measure changes plays also a significant role for a variety of option pricing models when underlying stock returns are discontinuous. However, for the assumption of a constant interest rate, sometimes it is also useful. The European call in Remark13 is a good example.


## 3.1. The Model

Let(
Ω
,
F
,
P
~
) be a probability space on which are defined a three-dimensional Brownian motion W
~
(
t
)
=
(
W
~
1
(
t
)
,
W
~
2
(
t
)
,
W
~
3
(
t
)
), 0
≤
t
≤
T, and two Poisson processes N
1
(
t
) and N
2
(
t
), 0
≤
t
≤
T, and P
~ denote the risk neutral probability measure. In particular, W
~
a
(
t
), W
~
b
(
t
), and W
~
c
(
t
) are independent Brownian motions. Let F
(
t
), 0
≤
t
≤
T, be the filtration generated by the Brownian motion and the Poisson processes. Let λ
1 and λ
2 be the intensities of the Poisson processes N
1
(
t
) and N
2
(
t
), respectively. We define two compound Poisson processes as follows:
(18)
Q
i
t
=
∑
j
=
1
N
i
t
‍
Y
i
j
,
-
1
<
X
i
j
,
i
=
1,2
,
where Y
i
1
,
Y
i
2
,
…, random variables take values in the set {
y
i
1
,
y
i
2
,
…
,
y
i
M
i
2
}. The random variables Y
i
1
,
Y
i
2
,
…, are independent and identically distributed, with P
{
Y
i
j
=
y
i
m
i
}
=
p
(
y
i
m
i
), m
i
=
1,2
,
…
,
M
i. We assume E
Y
i
j
=
∑
m
i
=
1
M
i
‍
y
i
m
i
p
(
y
i
m
i
)
=
β
i. The evolutions of the discounted stock prices will be modeled by the stochastic differential equations:
(19)
d
D
t
S
i
t
=
D
(
t
)
S
i
(
t
)
v
i
·
d
W
~
i
(
t
)
+
D
(
t
)
S
i
(
t
-
)
d
(
Q
i
(
t
)
-
β
i
λ
i
t
)
,
i
=
1,2
.


## 3.2. Pricing the Spread Option

Lemma 9.
Under the measuresP
~
S
i, i
=
1,2, the process W
~
S
i
(
t
)
=
(
W
~
a
(
t
)
,
W
~
b
(
t
)
,
W
~
c
(
t
)
)
-
v
i
t is a three-dimensional Brownian motion. Q
i
(
t
) is a compound Poisson process with intensity (
1
+
β
i
)
λ
i and independent, identically distributed jump sizes satisfying
(20)
P
~
S
i
Y
i
k
=
y
i
m
i
=
1
+
y
i
m
i
p
y
i
m
i
1
+
β
i
,
k
=
1,2
,
…
,
m
i
=
1,2
,
…
,
M
i
,
and the processes W
~
S
i
(
t
) and Q
i
(
t
) are independent.The proof of Lemma9 will be given in Appendix C.Proposition 10.
Under the following measures, the following price dynamics holds.(1)
Under theT-forward measure P
~
T, the process S
i
(
t
)
/
B
(
t
,
T
), i
=
1,2, is a martingale. Moreover,
(21)
d
S
i
(
t
)
B
(
t
,
T
)
=
S
i
t
-
B
t
,
T
v
i
-
v
3
t
,
T
·
d
W
~
T
t
+
d
Q
i
t
-
β
i
λ
i
t
,
where W
~
T
(
t
)
=
(
W
~
a
(
t
)
,
W
~
b
(
t
)
,
W
~
c
(
t
)
)
-
∫
0
t
‍
v
3
(
u
,
T
)
d
u is a three-dimensional Brownian motion under P
~
T. The jump measure is unaffected by the change of measure.
(2)
Under the measuresP
~
S
i, i
=
1,2, the processes B
(
t
,
T
)
/
S
i
(
t
) and S
j
(
t
)
/
S
i
(
t
), j
=
1,2, i
≠
j, are a martingale. Moreover,
(22)
d
B
(
t
,
T
)
S
i
(
t
)
=
B
(
t
,
T
)
S
i
(
t
-
)
∏
k
=
1
N
i
(
t
)
‍
Y
i
k
1
+
Y
i
k
-
β
i
λ
i
t
(
v
3
(
t
,
T
)
-
v
i
)
·
d
W
~
S
i
(
t
)
hhhhhhh
-
d
∏
k
=
1
N
i
(
t
)
‍
Y
i
k
1
+
Y
i
k
-
β
i
λ
i
t
,
d
S
j
t
S
i
t
=
S
j
t
-
S
i
t
-
∏
k
=
1
N
i
t
‍
Y
i
k
1
+
Y
i
k
-
β
i
λ
i
t
v
j
-
v
i
·
d
W
~
S
i
t
hhhhhhh
+
d
(
Q
j
(
t
)
-
β
j
λ
j
t
)
hhhhhhh
-
d
∏
k
=
1
N
i
(
t
)
‍
Y
i
k
1
+
Y
i
k
-
β
i
λ
i
t
,
where W
~
S
i
(
t
)
=
(
W
~
a
(
t
)
,
W
~
b
(
t
)
,
W
~
c
(
t
)
)
-
v
i
t is a three-dimensional Brownian motion under P
~
S
i.The proof of Proposition10 will be given in Appendix D.By Proposition10 and a similar discussion as the one in Theorem 5, we may obtain the following result.Theorem 11.
The price of the spread option, assuming the risk-neutral dynamics given in (3) and (19), is
(23)
V
t
=
∑
i
=
0
∞
‍
∑
j
=
0
∞
∑
L
=
1
2
‍
P
~
S
L
N
1
(
T
)
-
N
1
t
=
i
,
‍
hhhhhhhhhhhh
N
2
T
-
N
2
(
t
)
=
j
S
L
(
t
)
V
L
(
t
,
i
,
j
)
hhhhhhh
+
P
~
T
N
1
T
-
N
1
t
=
i
,
hhhhhhhhhhhh
N
2
(
T
)
-
N
2
t
=
j
hhhhhhh
×
K
B
t
,
T
V
3
t
,
i
,
j
∑
L
=
1
2
‍
with
(24)
V
1
t
,
i
,
j
=
P
~
S
1
1
-
∏
k
=
1
i
‍
∏
L
=
1
j
‍
1
+
Y
2
L
1
+
Y
1
k
A
21
(
t
)
e
B
21
(
t
)
ξ
11
hhhhh
-
∏
k
=
1
i
‍
1
1
+
Y
1
k
K
A
31
(
t
)
e
B
31
(
t
)
ξ
12
>
0
,
V
2
t
,
i
,
j
=
-
P
~
S
2
∐
k
=
1
i
‍
∐
L
=
1
j
‍
1
+
Y
1
k
1
+
Y
2
L
A
12
(
t
)
e
B
12
(
t
)
ξ
21
-
1
hhhhhh
-
∏
L
=
1
j
‍
1
1
+
Y
2
L
K
A
32
(
t
)
e
B
32
(
t
)
ξ
22
>
0
,
V
3
t
,
i
,
j
=
-
P
~
T
∏
k
=
1
i
‍
1
+
Y
1
k
A
13
t
e
B
13
t
ξ
31
hhhhhh
-
∏
L
=
1
j
‍
(
1
+
Y
2
L
)
A
23
hhhhihhhh
×
(
t
)
e
B
23
(
t
)
ξ
31
-
K
>
0
∏
L
=
1
j
‍
,
where
(25)
A
i
3
t
=
S
i
(
t
)
B
(
t
,
T
)
exp
⁡
-
1
2
∫
t
T
‍
(
v
i
-
v
3
(
u
,
T
)
)
2
d
u
h
h
h
h
h
h
h
h
h
h
-
β
i
λ
i
(
T
-
t
)
1
2
∫
t
T
‍
,
A
3
i
t
=
B
(
t
,
T
)
S
i
(
t
)
exp
⁡
-
1
2
∫
t
T
‍
(
v
3
(
u
,
T
)
-
v
i
)
2
d
u
hhhhhihhhhi
+
β
i
λ
i
(
T
-
t
)
1
2
∫
t
T
‍
,
A
i
j
t
=
S
i
(
t
)
S
j
(
t
)
exp
⁡
-
1
2
∫
t
T
‍
(
v
i
-
v
j
)
2
d
u
hhhhhhhii
+
β
i
λ
i
T
-
t
-
β
j
λ
j
T
-
t
1
2
∫
t
T
‍
,
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
i
,
j
=
1,2
,
i
≠
j
.Remark 12.
We can explore further applications of the formula beyond the case that the jump distribution is finitely many jump sizes. Suppose the jump sizesy
i
j, i
=
1,2, j
=
1,2
,
…, have a density f
i
(
y
i
) rather than a probability mass function P
(
y
i
1
)
,
P
(
y
i
2
)
,
…
,
P
(
y
i
M
i
). These densities are strictly positive on a set B
⊂
(
-
1
,
∞
) and zero elsewhere. In the case β
i
=
∫
-
1
∞
‍
f
i
(
x
)
d
x. By Lemma 9 and the definition of multiple integrals, we have that, under the measures P
~
S
i, i
=
1,2, the jumps in Q
i
(
t
) are independent and identically distributed with density (
1
+
y
i
)
f
i
(
y
i
)
/
(
1
+
β
i
).Remark 13.
We can also apply the formula to a catastrophe option or European call in a jump model. The price of the European call, assuming the risk-neutral dynamics given in (3) and (19) (i
=
1), is
(26)
V
t
=
∑
i
=
0
∞
‍
P
~
S
1
N
1
T
-
N
1
t
=
i
S
1
t
V
1
t
,
i
hhhh
+
P
~
T
(
N
1
(
T
)
-
N
1
(
t
)
=
i
)
hhhh
×
K
B
(
t
,
T
)
V
3
(
t
,
i
,
j
)
P
~
S
1
N
1
T
-
N
1
t
=
i
S
1
t
V
1
t
,
i
,
where
(27)
V
1
(
t
,
i
)
=
P
~
S
1
1
-
∏
k
=
1
i
‍
1
1
+
Y
1
k
K
A
31
(
t
)
e
B
31
(
t
)
ξ
12
>
0
,
V
3
(
t
,
i
,
j
)
=
-
P
~
T
∏
k
=
1
i
‍
(
1
+
Y
1
k
)
A
13
(
t
)
e
B
13
(
t
)
ξ
31
-
K
>
0
.Remark 14.
The technique of measure changes plays also a significant role for a variety of option pricing models when underlying stock returns are discontinuous. However, for the assumption of a constant interest rate, sometimes it is also useful. The European call in Remark13 is a good example.


## 4. Spread Option with Correlated Credit Risk and Stochastic Interest Rates

The section constructs a spread option pricing model that allows not only for correlation between the option’s underlying assets and the credit risk of the counterparty, but also for the interest rate to be random. For counterparty risk, we adopt some setting in the model of Klein [16].The value of the asset of the counterpartyS
4
(
t
) will be modeled by the stochastic differential equation:
(28)
d
(
D
(
t
)
S
4
(
t
)
)
=
D
(
t
)
S
i
(
t
)
v
4
·
d
W
~
4
(
t
)
,
where W
~
4
(
t
)
=
(
W
~
a
(
t
)
,
W
~
b
(
t
)
,
W
~
c
(
t
)
,
W
~
d
(
t
)
) and v
4
=
(
v
41
,
v
42
,
v
43
,
v
44
). The value at time t
∈
[
0
,
T
] of a spread option with counterparty risk, expiring at time T with strike price K, is
(29)
V
t
=
E
~
(
1
-
α
)
S
4
(
T
)
γ
e
-
∫
t
T
‍
R
(
u
)
d
u
S
1
T
-
S
2
T
-
K
+
h
h
h
×
(
1
-
α
)
S
4
(
T
)
γ
I
{
S
4
(
T
)
≥
γ
*
}
+
I
{
S
4
(
T
)
<
γ
*
}
h
h
h
h
h
h
×
(
1
-
α
)
S
4
(
T
)
γ
∣
F
(
t
)
e
-
∫
t
T
‍
R
(
u
)
d
u
.
At the exercise date T, a credit loss occurs if S
4
(
T
) is less than some amount γ
*. In the event of a credit loss, only the proportion (
1
-
α
)
S
4
(
T
)
/
γ of the nominal claim is paid out by the counterparty, where α represents the deadweight costs associated with bankruptcy expressed as a percentage of the value of the assets of the counterparty. γ
* may be less than γ due to the possibility of a counterparty continuing in operation even while S
4
(
T
) is less than γ. These deadweight costs include the direct cost of the bankruptcy or reorganization process, as well as the effects of distress on the business operations of the firm. If S
4
(
T
) is greater than γ
* no credit loss occurs and the time T claim receives full payment. We also refer the reader to the model of Klein [16] (a model of expected credit loss) for its relative discussion.Theorem 15.
The price of the spread option with correlated credit risk, assuming the risk-neutral dynamics given in (2), (3), and (28), is V
(
t
)
=
V
1
(
t
)
+
V
2
(
t
)
+
V
3
(
t
) with
(30)
V
1
t
=
S
1
(
t
)
E
I
{
1
-
A
21
(
t
)
e
B
21
(
t
)
ξ
11
-
K
A
31
(
t
)
e
B
31
(
t
)
ξ
12
>
0
}
I
¯
+
A
43
(
t
)
e
B
43
(
t
)
ξ
3
I
_
h
h
h
h
h
h
h
×
I
{
1
-
A
21
(
t
)
e
B
21
(
t
)
ξ
11
-
K
A
31
(
t
)
e
B
31
(
t
)
ξ
12
>
0
}
,
V
2
t
=
-
S
2
(
t
)
E
I
{
A
12
(
t
)
e
B
12
(
t
)
ξ
21
-
1
-
K
A
32
(
t
)
e
B
32
(
t
)
ξ
22
>
0
}
I
¯
+
A
43
(
t
)
e
B
43
(
t
)
ξ
3
I
_
h
h
h
i
h
h
h
h
h
×
I
{
A
12
(
t
)
e
B
12
(
t
)
ξ
21
-
1
-
K
A
32
(
t
)
e
B
32
(
t
)
ξ
22
>
0
}
,
V
3
t
=
-
K
B
(
t
,
T
)
E
I
{
A
13
(
t
)
e
B
13
(
t
)
ξ
31
-
A
23
(
t
)
e
B
23
(
t
)
ξ
32
-
K
>
0
}
I
¯
+
A
43
(
t
)
e
B
43
(
t
)
ξ
3
I
_
h
h
h
h
h
h
h
h
h
h
h
×
I
{
A
13
(
t
)
e
B
13
(
t
)
ξ
31
-
A
23
(
t
)
e
B
23
(
t
)
ξ
32
-
K
>
0
}
,
where
(31)
A
43
t
=
S
4
t
B
t
,
T
exp
⁡
-
1
2
∫
t
T
‍
v
4
2
-
v
3
u
,
T
2
d
u
,
B
43
(
t
)
=
∫
t
T
‍
v
4
-
v
3
u
,
T
2
d
u
.
I
¯
=
I
{
A
43
(
t
)
e
B
43
(
t
)
ξ
3
≥
γ
*
}
,
I
_
=
I
{
A
43
(
t
)
e
B
43
(
t
)
ξ
3
<
γ
*
}
.
The random variables ξ
3 and ξ
i
1 and ξ
3 and ξ
i
2 have the bivariate standard normal distribution with correlation ρ
i
1
~ and ρ
i
2
~, respectively, i
=
1,2
,
3, where
(32)
ρ
~
k
1
=
(
v
j
-
v
k
)
·
∫
t
T
‍
v
4
-
v
3
u
,
T
2
d
u
v
j
-
v
k
∫
t
T
‍
v
4
-
v
3
u
,
T
2
d
u
,
ρ
~
k
2
=
∫
t
T
‍
v
3
u
,
T
-
v
k
2
d
u
·
∫
t
T
‍
v
4
-
v
3
u
,
T
2
d
u
∫
t
T
‍
v
3
u
,
T
-
v
k
2
d
u
∫
t
T
‍
v
4
-
v
3
u
,
T
2
d
u
,
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
k
,
j
=
1,2
,
k
≠
j
,
ρ
~
3
L
=
∫
t
T
‍
v
L
-
v
3
u
,
T
2
d
u
·
∫
t
T
‍
v
4
-
v
3
u
,
T
2
d
u
∫
t
T
‍
v
L
-
v
3
u
,
T
2
d
u
∫
t
T
‍
v
4
-
v
3
u
,
T
2
d
u
,
hhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhh
L
=
1,2
.Remark 16.
Theorem15 also presents a formula for a vulnerable Black-Scholes option with stochastic interest rates and correlated credit risk, which improves the result of the vulnerable Black-Scholes option in a model of Klein [16]. The price of a vulnerable Black-Scholes option, assuming the risk-neutral dynamics given in (2) (i
=
1), (3), and (28), is V
(
t
)
=
V
1
(
t
)
+
V
3
(
t
) with
(33)
V
1
t
=
S
1
(
t
)
E
×
I
¯
+
A
43
(
t
)
e
B
43
(
t
)
ξ
3
I
_
I
{
1
-
K
A
31
(
t
)
e
B
31
(
t
)
ξ
12
>
0
}
,
V
3
t
=
-
K
B
(
t
,
T
)
E
×
I
¯
+
A
43
(
t
)
e
B
43
(
t
)
ξ
3
I
_
I
{
A
13
(
t
)
e
B
13
(
t
)
ξ
31
-
K
>
0
}
.


## 5. Numerical Experiments

In this section, we use regular grid method (the area of each square of the grid is 0.000025) to the calculation of the formula when underlying stock returns are continuous and a mixture of both the regular grid method (the area of each square of the grid is 0.0001) and a Monte Carlo method (a random sample of 50000 forN
(
T
)
-
N
(
t
)
≥
1) to the one when underlying stock returns are discontinuous. We omit the numerical results for the spread option with credit risk to save space. However, they are available from the authors upon request.Our modeling assumption is supplied by the following dynamics of the discounted price processesD
(
t
)
S
1 and D
(
t
)
S
2 and the short interest rate R
(
t
):
(34)
d
(
D
(
t
)
S
1
(
t
)
)
=
σ
1
D
(
t
)
S
1
(
t
)
d
W
~
a
(
t
)
,
d
D
(
t
)
S
2
(
t
)
=
ρ
12
σ
2
D
(
t
)
S
2
(
t
)
d
W
~
a
(
t
)
+
1
-
ρ
12
2
σ
2
D
(
t
)
B
(
t
,
T
)
d
W
~
b
(
t
)
,
d
R
t
=
γ
(
τ
-
R
(
t
)
)
d
t
-
σ
R
ρ
13
d
W
~
a
(
t
)
-
σ
R
η
1
d
W
~
b
(
t
)
-
σ
R
η
2
d
W
~
c
(
t
)
,
where η
1
=
ρ
23
-
ρ
12
ρ
13
/
1
-
ρ
12
2, η
2
=
1
-
ρ
13
2
-
η
1
2. We see that σ
i
2
,
i
=
1,2, is the instantaneous variance of the return on the underlying stock S
i
(
t
), and σ
R is the instantaneous volatility of short-term interest rate. (
W
~
a
(
t
)
,
ρ
13
W
~
a
(
t
)
+
η
1
W
~
b
(
t
)
+
η
2
W
~
c
(
t
)
) and (
ρ
12
W
~
a
(
t
)
+
1
-
ρ
12
2
W
~
b
(
t
)
,
ρ
13
W
~
a
(
t
)
+
η
1
W
~
b
(
t
)
+
η
2
W
~
c
(
t
)
) are a pair of correlated Brownian motions with correlation coefficients ρ
13 and ρ
23, respectively. γ and τ are the respective constant reversion rate and long-term mean in the drift coefficient of the Vasicek dynamics of the interest rate. The price of the T-maturity zero coupon bond B
(
t
,
T
) is given by B
t
,
T
=
exp
⁡
⁡
(
A
t
,
T
-
C
t
,
T
R
t
), where
(35)
A
t
,
T
=
τ
-
σ
R
2
2
γ
2
C
t
,
T
-
T
-
t
-
σ
R
2
4
γ
C
2
t
,
T
,
C
t
,
T
=
1
γ
1
-
exp
⁡
-
γ
T
-
t
,
and the bond price satisfies the following:
(36)
d
B
(
t
,
T
)
=
B
(
t
,
T
)
R
(
t
)
d
t
+
B
(
t
,
T
)
C
(
t
,
T
)
σ
R
ρ
13
d
W
~
a
(
t
)
+
B
(
t
,
T
)
C
(
t
,
T
)
σ
R
η
1
d
W
~
b
(
t
)
+
B
(
t
,
T
)
C
(
t
,
T
)
σ
R
η
2
d
W
~
b
(
t
)
.
Some parameters in Figures 1–4 and Tables 1–3 are fixed (t
=
0, T
=
1, τ
=
5
%, γ
=
0.3).Table 1
Spread option price.
K
λ
1
0
0.04
0.08
0.12
0.16
0.2
2
9.0963
9.3481
9.5591
9.7337
9.8761
9.9900
4
8.1906
8.5078
8.7794
9.0103
9.2046
9.3664
6
7.3595
7.6421
7.8821
8.0840
8.2522
8.3903
8
6.5753
6.8796
7.1398
7.3606
7.5464
7.7009
10
5.8556
6.1613
6.4232
6.6459
6.8336
6.9901Table 2
Spread option price.
ρ
12
λ
1
0
0.04
0.08
0.12
0.16
0.2
−0.8
10.5947
10.8515
11.0598
11.2249
11.3516
11.4443
−0.6
9.7783
10.0964
10.3634
10.5845
10.7648
10.9087
−0.4
8.9031
9.1473
9.3474
9.5083
9.6344
9.7297
−0.2
7.9746
8.2249
8.4324
8.6018
8.7374
8.8429
0
6.9684
7.2559
7.4994
7.7034
7.8722
8.0098
0.2
5.8622
6.1258
6.3488
6.5356
6.6902
6.8161
0.4
4.6314
4.9373
5.2013
5.4276
5.6203
5.7830
0.6
3.2155
3.5425
3.8272
4.0741
4.2872
4.4702
0.8
1.5161
1.8841
2.2092
2.4958
2.7479
2.9691Table 3
Spread option price.
ρ
13
λ
1
0
0.04
0.08
0.12
0.16
0.2
−0.8
6.0963
6.4456
6.7721
7.0777
7.3642
7.6332
−0.6
6.0512
6.3431
6.6131
6.8632
7.0952
7.3110
−0.4
6.0086
6.3172
6.5990
6.8566
7.0924
7.3086
−0.2
5.9608
6.2466
6.5026
6.7319
6.9372
7.1210
0
5.9153
6.2018
6.4529
6.6722
6.8632
7.0287
0.2
5.8629
6.1612
6.4163
6.6325
6.8140
6.9644
0.4
5.8057
6.0576
6.2611
6.4216
6.5440
6.6328
0.6
5.7472
6.0189
6.2281
6.3818
6.4862
6.5469
0.8
5.6825
5.9355
6.1092
6.2126
6.2539
6.3406Figure 1
The price of the spread option varies with the changes in underlying asset pricesS
1
(
t
) and S
2
(
t
).Figure 2
The price of the spread option varies with the changes in the volatilitiesσ
1 and σ
2.The price of the spread option varies with the changes in the correlation coefficientsρ
12 or ρ
13.
(a)
(b)The price of the spread option varies with the changes in the spot interest rateR
(
t
) and the volatility σ
R.
(a)
(b)Figure1 shows how the price of the spread option varies with the changes in underlying asset prices S
1
(
t
) and S
2
(
t
). The remaining model parameters are K
=
10, R
(
t
)
=
1
%, σ
1
=
20
%, σ
2
=
20
%, σ
R
=
15
%, ρ
12
=
0.2, ρ
13
=
0.2, and ρ
23
=
0.01. Figure 2 shows how the price of the spread option varies with the changes in the volatilities σ
1 and σ
2. Larger volatilities induce more uncertainty and the option price increases as the volatility σ
1 or σ
2 increases. The remaining model parameters are K
=
10, R
(
t
)
=
1
%, S
1
=
100, S
2
=
100, σ
R
=
15
%, ρ
12
=
0.2, ρ
13
=
0.2, and ρ
23
=
0.01. Figure 3 shows how the price of the spread option varies with the changes in the correlation coefficients ρ
12 or ρ
13. The changes in ρ
12 or ρ
13 significantly affect the value of the option. If the asset price of S
1 is strongly negatively correlated to interest rates and the asset price increases, then the probability that the option will be deep in the money at maturity date increases. The remaining model parameters are K
=
10, R
(
t
)
=
1
%, S
1
=
100, S
2
=
100, σ
1
=
20
%, σ
2
=
20
%, σ
R
=
15
%, ρ
23
=
0.01, ρ
12
=
0.2 (in Figure 3(b)), and ρ
13
=
0.2 (in Figure 3(a)). Figure 4 shows how the price of the spread option varies with the changes in the spot interest rate R
(
t
) and the volatility σ
R. If the the spot interest rate R
(
t
) increases, then the price of the zero coupon B
(
t
,
T
) decreases, which increases the probability that the option becomes in-the-money. The remaining model parameters in Figures 4(a) and 4(b) are K
=
10, S
1
=
100, S
2
=
100, σ
1
=
20
%, σ
2
=
20
%, σ
R
=
15
%, ρ
12
=
0.2, ρ
13
=
0.5, and ρ
23
=
0.01 and K
=
30, R
(
t
)
=
5
%, S
1
=
100, S
2
=
80, σ
1
=
20
%, σ
2
=
20
%, ρ
12
=
0.1, ρ
13
=
-
0.5, and ρ
23
=
0.5, respectively.Now we show some numerical results for the spread option pricing in a jump model. All the following results are consistent with the behavior of the spread option price when the returns are discontinuous for all the two underlying stocks and, to simplify the discussion, we will focus on the casesλ
1
≠
0 and λ
2
=
0. We only assume (λ
2
=
0) that the dynamics of the discounted price process D
(
t
)
S
1 is
(37)
d
D
t
S
1
t
=
σ
1
D
(
t
)
S
1
(
t
)
d
W
~
a
(
t
)
+
D
(
t
)
S
1
(
t
-
)
d
(
Q
1
(
t
)
-
β
1
λ
1
t
)
,
where Y
11
,
Y
12
,
…, have a density f
1
(
y
1
)
≈
(
1
/
2
π
0.4
)
exp
⁡
(
-
y
1
2
/
0.32
) for -
1
<
y
1
<
1.Table1 shows how the price of the spread option varies with the changes in the exercise price K and the Poisson active rate λ
1. Larger intensity induces more uncertainty and the option price increases as the Poisson active rate increases. The remaining model parameters are R
(
t
)
=
1
%, S
1
=
100, S
2
=
100, σ
1
=
20
%, σ
2
=
20
%, σ
R
=
15
%, ρ
12
=
0.2, ρ
13
=
0.2, and ρ
23
=
0.01. Table 2 shows how the price of the spread option varies with the changes in the correlation coefficient ρ
12 and the Poisson active rate λ
1. The remaining model parameters are K
=
10, R
(
t
)
=
1
%, S
1
=
100, S
2
=
100, σ
2
=
20
%, σ
R
=
15
%, ρ
13
=
0.2, and ρ
23
=
0.01. Table 3 shows how the price of the spread option varies with the changes in the correlation coefficient ρ
13 and the Poisson active rate λ
1. The remaining model parameters are K
=
10, R
(
t
)
=
1
%, S
1
=
100, S
2
=
100, σ
2
=
20
%, σ
R
=
15
%, ρ
12
=
0.2, and ρ
23
=
0.01. For K, ρ
12, and ρ
13, the results in Tables 1–3 are consistent with the behavior of the spread option price when underlying stock returns are continuous.


## 6. Concluding Remarks

In this paper, we have explored further analysis of spread option pricing beyond the case that the interest rate is a constant. The model of stochastic interest rates and the technique of measure changes are useful in a variety of asset pricing models. Through numerical experiments and simulations, we demonstrated that the changes in correlation coefficients, the volatilities of the returns on underlying stocks, and the Poisson active rate significantly affect the price, and stochastic interest rates play a significant role. In general, existing methods for pricing spread options can be divided into two groups: numerical methods and analytical approximations. Analytical approximation methods are generally faster than numerical methods but often lack accuracy and robustness. For the spread option pricing formula (without counterparty risk) when underlying stock returns are continuous, by the technique of measure changes, the integrand becomes simple. In the case, the method is both accurate and fast (by regular grid method, it is easy to derive suitable accurate price). By a mixture of both the regular grid method and a Monte Carlo method to the calculation of the formula when underlying stock returns are discontinuous, we can obtain a suitable approximate price.



---

*Source: 734265-2014-11-13.xml*