Next Article in Journal
Testing the Homogeneity of Differences between Two Proportions for Stratified Bilateral and Unilateral Data across Strata
Previous Article in Journal
PFA-Nipals: An Unsupervised Principal Feature Selection Based on Nonlinear Estimation by Iterative Partial Least Squares
Previous Article in Special Issue
Pricing Equity-Indexed Annuities under a Stochastic Dividend Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Pricing European Vulnerable Options with Jumps and Stochastic Default Obstacles Barrier under Regime Switching

School of Economics, Jinan University, Guangzhou 510632, China
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(19), 4155; https://doi.org/10.3390/math11194155
Submission received: 5 September 2023 / Revised: 20 September 2023 / Accepted: 29 September 2023 / Published: 3 October 2023

Abstract

:
In this paper, we propose an enhanced model for pricing vulnerable options. Specifically, our model assumes that parameters such as interest rates, jump intensity, and asset value volatility are governed by an observable continuous-time finite-state Markov chain. We take into account European vulnerable options that are exposed to both default risk and rare shocks from underlying and counterparty assets. We also consider stochastic default barriers driven by a regime-switching model and geometric Brownian motion, thus improving upon the assumption of fixed default barriers. The risky assets follow a related jump-diffusion process, whereas the default barriers are influenced by a geometric Brownian motion correlated with the risky assets. Within the framework of our model, we derive an explicit pricing formula for European vulnerable options. Furthermore, we conduct numerical simulations to examine the effects of default barriers and other related parameters on option prices. Our findings indicate that stochastic default barriers increase credit risk, resulting in a decrease in option prices. By considering the aforementioned factors, our research contributes to a better understanding of pricing vulnerable options in the context of counterparty credit risk in over-the-counter trading.

1. Introduction

The over-the-counter (OTC) market represents a significant portion of the financial market, with many financial institutions actively trading options on it. Unlike exchange-traded markets, OTC transactions are unorganized and lack transparency, which increases the risk of counterparty default. This means that one party may fail to fulfill its contractual obligations. Since the global financial crisis of 2007–2008, participants in the OTC market have become increasingly concerned about the credit default risk associated with options. The structural model, also known as the firm value model, is a key tool for pricing options with credit risk. This article explores option pricing within the framework of the firm value model. Johnson and Stulz [1] were the first to define options with credit risk as vulnerable options.
There are several models for vulnerable option pricing in the literature. Klein [2,3] proposed a model that allows for correlation between the underlying asset of the option and the counterparty’s asset and assumes that the default barrier is the fixed liability of the option writer. He also proposed an improved method that includes both the fixed liabilities of the option writer and losses from vulnerable option shorts in the default barrier. Hui [4] proposed a European vulnerable option pricing model with a stochastic default barrier representing the option writer’s liability. Wang [5] assumed that both the counterparty’s assets and related assets were driven by a related jump-diffusion process and derived a pricing formula for vulnerable options. These studies suggest that the default barrier has an impact on the pricing of European vulnerable options. However, they only consider the impact of barriers on options and do not take into account how a country’s economic status at different times affects option prices.
Since the dynamics of risky asset values are not continuous, some literature introduces jumping processes, compound Poisson processes, or self-exciting Hawkes processes to better describe their value models. These studies discuss the pricing formula for European vulnerable options (see XU [6], TIAN [7], Ma [8], and Zhou [9]). Wang [10] proposed decomposing the underlying asset risk and counterparty assets into systemic and non-systemic risks and used the obtained pricing formula to discuss the difference between vulnerable option prices. In addition, many empirical analyses have concluded that the volatility of risky assets is not fixed. Therefore, some literature assumes that the underlying asset is driven by stochastic volatility models such as the Heston model and mean reversion process and discusses the pricing formula for European vulnerable options (see Liu [11], Yang [12], Lee [13], Wang [14]). Yoon [15] uses double Merlin transforms to obtain closed-form pricing formulas for European vulnerable options under constant and Hull–White interest rate models (see Liu [16]). Huang [17] investigated the pricing of European-style vulnerable options in the presence of non-affine stochastic volatility and double exponential jumps in the underlying asset price process.
Most option pricing models are within the framework of Brownian motion, which makes it difficult to capture external factors such as the state of the national economy. The literature also incorporates vulnerable option pricing under the regime-switching model. Elliott [18,19] was the first to consider introducing regime switching into the option pricing model. Wang [20] uses the local risk minimization method to obtain the optimal hedging strategy for European vulnerable call options. When the option writer’s asset value is depressed by other financial institutions due to poor sales, Yang [21] used polynomial change numerical techniques to obtain a semi-analytical pricing formula for European vulnerable options. Assuming that the default barrier in the vulnerable option model is fixed, Wang [22] assumes that risk assets are driven by regime-switching geometric Brownian motion and obtains a European vulnerable option pricing formula. NIU [23] assumes that the value of risk assets is driven by a relevant regime-switching double exponential jump-diffusion model and uses two-dimensional Laplace transforms to obtain a European vulnerable option pricing formula. Han [24] assumes that the value of risk assets is driven by a relevant regime-switching jump-diffusion model and uses the Esscher transform to select an equivalent martingale measure to obtain a European vulnerable option pricing formula and discusses the impact of jump risk on European vulnerable option pricing. Capponi [25] presents an efficient method for pricing vulnerable contingent claims in a regime-switching market using a suitable change of probability measure and a Poisson series representation. Damircheli [26] investigated the modeling of default probability for publicly traded companies using a synchronous-jump regime-switching model and incorporating tempered stable processes. They proposed an efficient meshfree collocation method to solve the partial integro-differential equations associated with the credit risk model.
However, based on the above literature, we can see that the factors affecting the pricing of European vulnerable options include default barriers, risk asset value models, and national economic conditions. We found that there is currently no study on the pricing of European vulnerable options under a regime-switching model with stochastic default barriers. In light of this, we naturally consider proposing a model that is closer to the actual situation. This article assumes that risky assets are driven by a relevant regime-switching jump-diffusion model and that the default barrier follows a regime-awitching geometric Brownian motion. We obtain an explicit solution for the European vulnerable option pricing formula. In addition, we use numerical simulation to explore the impact of default barriers and parameters on option prices and find that stochastic default barriers increase credit risk, leading to lower option prices. We differ from Han [24] in that they assume that the default barrier is constant. We assume that the dynamics of the default barrier are regime-switching geometric Brownian motions related to the underlying assets and counterparty assets.
The structure of this paper is as follows: Section 2 establishes a regime-switching model for risk assets and default barriers. Section 3 uses the Esscher transform to select an equivalent martingale measure. Section 4 uses the risk-neutral pricing principle to obtain a pricing formula for European vulnerable options. Section 5 presents a numerical analysis of the proposed model using Monte Carlo simulation methods. The conclusion is given in Section 6.

2. Model Description

We construct a continuous-time financial market on a probability space Ω , F , P . We assume that all risky assets can be continuously traded in time [0, T]. Let X = X t ; t 0 denote a continuous-time Markov process on probability space Ω , F , P with a finite state space e 1 , e 2 , , e N , where e i = 0 , , 1 , , 0 R N with the i -th element is ‘1’. We assume that the states of the economy are modelled by a continuous-time Markov process X = X t ; t 0 . According to Elliott [27] semi-martingale decomposition for X, we have X t = X 0 + 0 t A X s d s + M t , where A = q i j N × N is Q- matrix of X and M t is a R N -valued martingale with respect to the natural filtration generated by X t ; t 0 .
Assumption 1.
The risk-free rate  r t  is related to the state of the economy. Let  r t = r , X t , where  r = r 1 , r 2 , , r N R N . For any  i = 1 ,   2 ,   3 ,   , N , we have  r i > 0 ;  ,  denotes the inner product of the vector. When  i j ,   e i ,   e j = 0 ; otherwise,  e i ,   e j = 1 .
Assumption 2.
Suppose that  S t  is the market value of the underlying assets at time t, and  V t  denote the value of the counterparty’s company assets at time t.  V t  and  S t  are driven by the following jump-diffusion processes:
d S t S t = α t k 1 λ t + λ 1 t d t + σ 1 t d W 1 t + e Z t 1 1 d N t + e Z t 1 1 d N 1 t
d V t V t = β t k 2 λ t + λ 2 t d t + σ 2 t d W 2 t + e Z t 2 1 d N t + e Z t 2 1 d N 2 t
where W 1 t ,   W 2 t are standard Brownian motions. N 1 t and N 2 t , respectively, indicate the individual shocks to the value of the underlying asset S t and the individual shocks to the value of the counterparty’s company assets V t , where N 1 t and N 2 t are pure scalar Poisson process with intensities λ 1 t and λ 2 t , respectively. N t refers to the common shocks to the corresponding assets, where N t are pure scalar Poisson process with intensity λ t .
For i = 1 ,   2 , we suppose that Z t i denotes the underlying asset’s jump amplitude when the jump happens, following a normal distribution with mean μ i and variance δ i 2 . At any time, s t ,   Z t i and Z s i are independent of each other; k i = E e Z t i 1 = e   μ i + 1 2 δ i 2 1 represents the average jump of the asset. α t and σ 1 t are the expected return rate of the underlying asset S t and the volatility of the underlying asset S t , respectively. β t and σ 2 t are the expected return rate of the value of the counterparty’s company assets V t and the volatility of the value of the counterparty’s company assets V t , respectively. We assume that α t , β t , σ 1 t , σ 2 t , λ 1 t , λ 2 t and λ t all depend on the economic state X. Let α t = α ,   X t , β t = β ,   X t , σ 1 t = σ 1 ,   X t , σ 2 t = σ 2 ,   X t , λ 1 t = λ 1 ,   X t , λ 2 t = λ 2 ,   X t , λ t = λ ,   X t , where α = α 1 , , α i , , α N R N , β = β 1 , , β i , , β N R N , σ 1 = σ 11 , , σ 1 i , , σ 1 N R N , σ 2 = σ 21 , , σ 2 i , , σ 2 N R N , λ = λ 1 , , λ i , , λ N R N , λ 1 = λ 11 , ,   λ 1 i , ,   λ 1 N R N , and λ 2 = λ 21 , ,   λ 2 i , ,   λ 2 N R N . For any j = 1 ,   2 ,   3 ,   , N , and i = 1 ,   2 , there are σ i j > 0 .
Assumption 3.
The value of the counterparty’s company liabilities, also called the default barrier, is represented by  D t , and  D t  is driven by the following geometric Brownian motion:
d D t D t = γ t d t + σ 3 t d W 3 t
where W 3 t is the standard Brownian motion, and σ 3 t is the volatility of the counterparty’s company liabilities. We assume that γ t and σ 3 t are related to the economic state X . Let γ t = γ ,   X t , σ 3 t = σ 3 ,   X t , where γ = γ , , γ i , , γ N R N , σ 3 = σ 31 , , σ 3 i , , σ 3 N R N . For any i = 1 ,   2 ,   3 ,   , N , there are σ 3 i > 0 .
Assumption 4.
We assume that the covariance matrix of the Brownian motion  W 1 t ,   W 2 t ,   W 3 t  on probability space  Ω , F , P  is as follows:
1 ρ 12 ρ 13 ρ 12 1 ρ 23 ρ 13 ρ 23 1 t
We assume that N 1 t , N 2 t , N t , Z t 1 and Z t 2 are independent of each other. Moreover, we assume that N 1 t , N 2 t , N t , Z t 1 and Z t 2 are independent of W 1 t ,   W 2 t ,   W 3 t . In addition, for the convenience of calculation in the following chapters, we express the relevant three-dimensional Brownian motion W 1 t ,   W 2 t , and W 3 t as
W 1 t W 2 t   W 3 t = 1 0 0 ρ 12 1 ρ 12 2 0 ρ 13 ρ 23 ρ 12 ρ 13 1 ρ 12 2 1 ρ 12 2 1 ρ 13 2 ρ 23 ρ 12 ρ 13 2 1 ρ 12 2 W 1 t W ¯ 2 t   W ¯ 3 t
where W 1 t W ¯ 2 t and W ¯ 3 t are independent Brownian motions.
Assumption 5.
At the option maturity T, if the value of the counterparty’s company assets is not greater than the value of the counterparty’s company liabilities, we believe that the option has defaulted. The option holder will receive  1 ω D T V T  times the notional amount from the option writer, where  ω  is the recovery rate. Therefore, at the maturity T, the payment function of the vulnerable European call option with a strike price of  K  is as follows:
S T K + I V T D T + 1 ω V T D T I V T < D T
where I · is represented as an indicator function.

3. Selection of Equivalent Martingale Measures

Due to the presence of jumps and the uncertainty introduced by regime switching, the market is incomplete and there are infinitely many equivalent martingale measures. Furthermore, it is not easy to obtain an equivalent martingale measure for two jump-diffusion processes with correlated jumps. It is important to note that while systemic risks can be hedged, non-systematic risks are specific to individual risk assets and cannot be hedged, so their premiums need not be considered. In this paper, we assume that the jump component represents systemic risk and its risk premium must be taken into account. We should measure changes in the jump component. In this section, we apply the regime-switching Esscher transform to select equivalent martingale measures for three stochastic processes with associated jump risks under regime switching.
Using Ito’s theorem, the solutions of Equations (1)–(3) are given by
S t S 0 = e L 1 t = e x p 0 t α s k 1 λ s + λ 1 s 1 2 σ 1 s 2 d s + 0 t σ 1 s d W 1 s + 0 t Z s 1 d N 1 s + 0 t Z s 1 d N s
V t V 0 = e L 2 t = e x p 0 t β s k 2 λ s + λ 2 s 1 2 σ 2 s 2 d s + 0 t σ 2 s d W 2 s + 0 t Z s 2 d N 2 s + 0 t Z s 2 d N s
D t D 0 = e L 3 t = e x p 0 t γ s 1 2 σ 3 s 2 d s + 0 t σ 3 s d W 3 s
where L 1 t is the return process of the underlying asset S t , L 2 t is the return process of the underlying asset V t , and L 3 t is the return process of the underlying asset D t . According to the above Equations (7)–(9), we have the three process L 1 t , L 2 t , and L 3 t as follows:
L 1 t = 0 t α s k 1 λ s + λ 1 s 1 2 σ 1 s 2 d s + 0 t σ 1 s d W 1 s + 0 t Z s 1 d N 1 s + 0 t Z s 1 d N s = C 1 t + J 1 t + J t
L 2 t = 0 t β s k 2 λ s + λ 2 s 1 2 σ 2 s 2 d s + 0 t σ 2 s d W 2 s + 0 t Z s 2 d N 2 s + 0 t Z s 2 d N s = C 2 t + J 2 t + J t
L 3 t = 0 t γ s 1 2 σ 3 s 2 d s + 0 t σ 3 s d W 3 s = C 3 t
Let F t X t 0 , T and F t i t 0 , T   i = 1 , 2 , 3 be the natural σ -filtrations generated by X t t 0 , T and L i t t 0 , T   i = 1 , 2 , 3 , respectively. For any i = 1 , 2 , 3 and t 0 , T , define a filtration H t i t 0 , T such that H t i is the σ - algebra with H t i =   F t X F t i . For any t 0 , T , let R t = H t 1 H t 2 H t 3 . Let V 1 t c 1 , V 2 t c 2 , V 3 t c 3 , V 1 t j 1 , V 2 t j 2 , V t j and V t j be the regime-switching Esscher transform parameters, and they are all related to the economic state of X . For any m c 1 , j 1 , l j , j , and n c 2 , j 2 , we assume that V 1 t m = V 1 m ,   X t , V 2 t n = V 2 n ,   X t , and V t l = V l ,   X t , where V 1 m = V 11 m , , V 1 i m , , V 1 N m R N , V 2 n = V 21 n , , V 2 i n , , V 2 N n R N , and V l = V 1 l , , V i l , , V N l R N . For the regime-switching Esscher transform parameter V 3 t c 3 , we assume that V 3 t c 3 = V 3 c 3 ,   X t , where V 3 c 3 = V 31 c 3 , , V 3 i c 3 , , V 3 N c 3 R N .
Under the given condition F t X , we use the regime-switching Esscher transform proposed by Elliott [18] to select an equivalent martingale measure. Define the equivalent measure of P ¯ as P , the Radon-Nikodym derivative of the regime-switching Esscher transform can be written as follows:
d P ¯ d P = e x p 0 T V 1 t c 1 d C 1 t + 0 T V 2 t c 2 d C 2 t + 0 T V 3 t c 3 d C 3 t E P e x p 0 T V 1 t c 1 d C 1 t + 0 T V 2 t c 2 d C 2 t + 0 T V 3 t c 3 d C 3 t | F T X × e x p 0 T V 1 t j 1 d J 1 t + 0 T V 2 t j 2 d J 2 t + 0 T V t j d J t + 0 T V t j d J t E P e x p 0 T V 1 t j 1 d J 1 t + 0 T V 2 t j 2 d J 2 t + 0 T V t j d J t + 0 T V t j d J t | F T X
By Equation (5) and calculation, we have
d P ¯ d P = e x p 0 T V 1 t c 1 σ 1 t d W 1 t + 0 T V 2 t c 2 σ 2 t d W 2 t + 0 T V 3 t c 3 σ 3 t d W 3 t × e x p 1 2 0 T V 1 t c 1 2 σ 1 t 2 + V 2 t c 2 2 σ 2 t 2 + V 3 t c 3 2 σ 3 t 2 d t + 0 T V t j Z t 1 + V t j Z t 2 d N t × e x p 0 T ρ 12 V 1 t c 1 σ 1 t V 2 t c 2 σ 2 t + ρ 13 V 1 t c 1 σ 1 t V 3 t c 3 σ 3 t + ρ 23 V 3 t c 3 σ 3 t V 2 t c 2 σ 2 t d t × e x p 0 T λ 1 t e x p μ 1 V 1 t j 1 + 1 2 δ 1 V 1 t j 1 2 1 d t + 0 T V 1 t j 1 Z t 1 d N 1 t × e x p 0 T λ 2 t e x p μ 2 V 2 t j 2 + 1 2 δ 2 V 2 t j 2 2 1 d t + 0 T V 2 t j 2 Z t 2 d N 2 t × e x p 0 T λ t e x p μ 1 V t j + μ 2 V t j + 1 2 δ 1 V t j 2 + 1 2 δ 2 V t j 2 1 d t
Delbaen [28] showed that the absence of arbitrage is essentially equivalent to the existence of equivalent martingale measures. Under an equivalent martingale measure, the price process of a risky asset is a martingale, which is known as the martingale condition. Due to the uncertainty introduced by the economic state X, the discount process must be defined on a larger filtration R. From the principles of asset pricing and no-arbitrage, it follows that under the equivalent martingale measures P ¯ , the discount process of risky assets is a martingale.
Theorem 1.
The discounting process of risk assets satisfies the following formula:
E P ¯ e x p 0 t r s d s S t R u = e x p 0 u r s d s S u  
E P ¯ e x p 0 t r s d s V t R u = e x p 0 u r s d s V u  
E P ¯ e x p 0 t r s d s D t R u = e x p 0 u r s d s D u  
if and only if the regime-switching Esscher transform parameters satisfy the following conditions:
α s k 1 λ s + λ 1 s r s + V 1 s c 1 σ 1 s 2 + ρ 12 σ 1 s V 2 s c 2 σ 2 s + ρ 13 V 3 s c 3 σ 1 s σ 3 s + λ ˜ 1 s k 1 j 1 + λ ˜ s k 1 j = 0
β s k 2 λ s + λ 2 s r s + V 2 s c 2 σ 2 s 2 + ρ 12 σ 1 s V 1 s c 1 σ 2 s + ρ 23 V 3 s c 3 σ 3 s + λ ˜ 2 s k 2 j 2 + λ ˜ s k 2 j = 0
γ s r s + V 3 s c 3 σ 3 s 2 + ρ 13 σ 1 s V 1 s c 1 σ 3 s + ρ 23 V 2 s c 2 σ 2 s σ 3 s = 0
where  λ ˜ 1 s = λ 1 s   e x p μ 1 V 1 s j 1 + 1 2 δ 1 V 1 s j 1 2 ,  λ ˜ 2 s = λ 2 s   e x p μ 2 V 2 s j 2 + 1 2 δ 2 V 2 s j 2 2 ,  λ ˜ s = λ s e x p μ 1 V s j + μ 2 V s j + 1 2 δ 1 V s j 2 + 1 2 δ 2 V s j 2 ,  k 1 j 1 = k 1 + 1 e δ 1 2 V 1 s j 1 1 ,  k 2 j 2 = k 2 + 1 e δ 2 2 V 2 s j 2 1 ,  k 1 j = k 1 + 1 e δ 1 2 V 1 s j 1 ,  k 2 j = k 2 + 1 e δ 2 2 V 2 s j 1 .
Proof. 
By Bayes’ rule, we have
  E P ¯ e 0 t r s d s S t R u = E P e 0 t r s d s S t d P ¯ d P R u = E P e 0 t r s d s S t d P ¯ d P R u E P d P ¯ d P R u = e x p 0 u r s d s S u E P e x p M S A + M S B + M S C + M S J = e x p 0 u r s d s S u
So, if and only if the following formula holds, the risk asset discounting process is martingale
α s k 1 λ s + λ 1 s r s + V 1 s c 1 σ 1 s 2 + ρ 12 σ 1 s V 2 s c 2 σ 2 s + ρ 13 V 3 s c 3 σ 1 s σ 3 s + λ ˜ 1 s k 1 j 1 + λ ˜ s k 1 j = 0
Similarly, for the risk asset V t , the following conditions can be proved:
β s k 2 λ s + λ 2 s r s + V 2 s c 2 σ 2 s 2 + ρ 12 σ 1 s V 1 s c 1 σ 2 s + ρ 23 V 3 s c 3 σ 2 s σ 3 s + λ ˜ 2 s k 2 j 2 + λ ˜ s k 2 j = 0
The following proves the conditions for the equivalent martingale measures of the risk asset D t . According to Baye’s rule, we also know
E P ¯ e 0 t r s d s D t R u = E P e 0 t r s d s D t d P ¯ d P R u = E P e 0 t r s d s D t d P ¯ d P R u E P d P ¯ d P R u = E P e x p u t γ s r s + V 3 s c 3 σ 3 s 2 + ρ 13 σ 1 s V 1 s c 1 σ 3 s + ρ 23 V 2 s c 2 σ 2 s σ 3 s d s e x p 0 u r s d s D u = e x p 0 u r s d s D u
We have γ s r s + V 3 s c 3 σ 3 s 2 + ρ 13 σ 1 s V 1 s c 1 σ 3 s + ρ 23 V 2 s c 2 σ 2 s σ 3 s = 0 .
The proof of Theorem 1 is complete. □
Theorem 2.
Under the equivalent martingale measure  P ¯ ,  S t ,  V t  and  D t  are driven by the following formula:
log S t S 0 = 0 t r s λ ˜ s k 1 j λ ˜ 1 s k 1 j 1 1 2 σ 1 s 2 d s + 0 t σ 1 s d W ˜ 1 s + 0 t Z ˜ s 1 d N ˜ 1 s + 0 t Z ^ s 1 d N ˜ s
log V t V 0 = 0 t r s λ ˜ s k 2 j λ ˜ 2 s k 2 j 2 1 2 σ 2 s 2 d s + 0 t σ 2 s d W ˜ 2 s + 0 t Z ˜ s 2 d N ˜ 2 s + 0 t Z ^ s 2 d N ˜ s
log D t D 0 = 0 t r s 1 2 σ 3 s 2 d s + 0 t σ 3 s d W ˜ 3 s
where  W ˜ 1 s ,  W ˜ 2 s  and  W ˜ 3 s  are the standard Brownian motion under the equivalent martingale measure  P ¯ , and satisfyes the following relationship:
W ˜ 1 t = W 1 t 0 t V 1 s c 1 σ 1 s d s 0 t ρ 12 V 2 s c 2 σ 2 s + ρ 13 V 3 s c 3 σ 3 s d s
W ˜ 2 t = W 2 t 0 t V 2 s c 2 σ 2 s d s 0 t ρ 12 V 1 s c 1 σ 1 s + ρ 23 V 3 s c 3 σ 3 s d s
W ˜ 3 t = W 3 t 0 t V 3 s c 3 σ 3 s d s 0 t ρ 13 V 1 s c 1 σ 1 s + ρ 23 V 2 s c 2 σ 2 s d s
l o g D t D 0 = 0 t r s 1 2 σ 3 s 2 d s + 0 t σ 3 s d W ˜ 3 s
The covariance matrix of the Brownian motion  W ˜ 1 t ,   W ˜ 2 t ,   W ˜ 3 t  is the same as that of  W 1 t ,   W 2 t ,   W 3 t .  N ˜ 1 t ,  N ˜ 2 t  and  N ˜ t  are Poisson processes with respective intensities  λ ˜ 1 t = λ 1 t   e x p μ 1 V 1 t j 1 + 1 2 δ 1 V 1 t j 1 2 ,  λ ˜ 2 t = λ 2 t   e x p μ 2 V 2 t j 2 + 1 2 δ 2 V 2 t j 2 2  and  λ ˜ t = λ t e x p μ 1 V t j + μ 2 V t j + 1 2 δ 1 V t j 2 + 1 2 δ 2 V t j 2 , respectively.  Z ˜ t 1 ,  Z ˜ t 2 ,  Z ^ t 1  and  Z ^ t 2  are normal random variables, and  Z ˜ t 1 ~ N μ 1 + V 1 t j 1 δ 1 2 ,   δ 1 2 ,  Z ˜ t 2 ~ N μ 2 + V 2 t j 2 δ 2 2 ,   δ 2 2 ,  Z ^ t 1 ~ N μ 1 + V 1 t j δ 1 2 ,   δ 1 2  and  Z ^ t 2 ~ N μ 2 + V 2 t j δ 2 2 ,   δ 2 2 .  k 1 j 1 = k 1 + 1 e δ 1 2 V 1 t j 1 1 ,  k 2 j 2 = k 2 + 1 e δ 2 2 V 2 t j 2 1 ,  k 1 j = k 1 + 1 e δ 1 2 V 1 t j 1  and  k 2 j = k 2 + 1 e δ 2 2 V 2 t j 1  are expressed as the average jumps of  Z ˜ t 1 ,  Z ˜ t 2 ,  Z ^ t 1  and  Z ^ t 2 , respectively.
Proof. 
According to the multi-dimensional Girsanov theorem, under the given conditions F t X and the equivalent martingale measure P ¯ , we know that W ˜ 1 t ,   W ˜ 2 t and W ˜ 3 t are Brownian motion such that Equations (29)–(31) are established.
Using the theorem T10 [29] (p. 241), we observe that Z ˜ t 1 , Z ˜ t 2 , Z ^ t 1 and Z ^ t 2 are normal random variables, i.e., Z ˜ t 1 ~ N μ 1 + V 1 t j 1 δ 1 2 ,   δ 1 2 , Z ˜ t 2 ~ N μ 2 + V 2 t j 2 δ 2 2 ,   δ 2 2 , Z ^ t 1 ~ N μ 1 + V 1 t j δ 1 2 ,   δ 1 2 and Z ^ t 2 ~ N μ 2 + V 2 t j δ 2 2 ,   δ 2 2 . We can calculate k 1 j 1 = E e Z ˜ t 1 1 = k 1 + 1 e δ 1 2 V 1 t j 1 1 , k 2 j 2 = E e Z ˜ t 2 1 = k 2 + 1 e δ 2 2 V 2 t j 2 1 , k 1 j = E e Z ^ t 1 1 = k 1 + 1 e δ 1 2 V 1 t j 1 and k 2 j = E e Z ^ t 2 1 = k 2 + 1 e δ 2 2 V 2 t j 1 .
The proof of Theorem 2 is complete. □
According to Theorem 1, we obtain the martingale condition, but from algebraic knowledge, we know that there are infinite solutions to Equations (18)–(20). In order to facilitate the numerical simulation in Section 5, we provide a set of special solutions. Let
α s k 1 λ s + λ 1 s r s + V 1 s c 1 σ 1 s 2 + ρ 12 σ 1 s V 2 s c 2 σ 2 s + ρ 13 V 3 s c 3 σ 1 s σ 3 s = 0 β s k 2 λ s + λ 1 s r s + V 2 s c 2 σ 2 s 2 + ρ 12 σ 1 s V 1 s c 1 σ 2 s + ρ 23 V 3 s c 3 σ 2 s σ 3 s = 0 γ s r s + V 3 s c 3 σ 3 s 2 + ρ 13 σ 1 s V 1 s c 1 σ 3 s + ρ 23 V 2 s c 2 σ 2 s σ 3 s = 0 λ 1 s   e x p μ 1 V 1 s j 1 + 1 2 δ 1 V 1 s j 1 2 k 1 + 1 e x p δ 1 2 V 1 s j 1 1 = 0 λ 2 s   e x p μ 2 V 2 s j 2 + 1 2 δ 2 V 2 s j 2 2 k 2 + 1 e x p δ 2 2 V 2 s j 2 1 = 0 λ s e x p μ 1 V s j + μ 2 V s j + 1 2 δ 1 V s j 2 + 1 2 δ 2 V s j 2 k 1 + 1 e δ 1 2 V s j 1 = 0 λ s e x p μ 1 V s j + μ 2 V s j + 1 2 δ 1 V s j 2 + 1 2 δ 2 V s j 2 k 1 + 1 e δ 2 2 V s j 1 = 0
Solving the above equation can obtain the regime-switching Esscher transform parameters as
V 1 s c 1 = π 1 ρ 23 2 1 σ 3 s σ 2 s + π 2 ρ 12 ρ 23 ρ 13 σ 1 s σ 3 s + π 3 ρ 13 ρ 23 ρ 12 σ 1 s σ 2 s σ 3 s σ 2 s σ 1 s 2 1 ρ 12 2 1 ρ 13 2 ρ 23 ρ 13 ρ 12 2
V 2 s c 2 = π 1 ρ 12 ρ 23 ρ 13 σ 3 s σ 2 s + π 2 ρ 13 2 1 σ 1 s σ 3 s + π 3 ρ 23 ρ 13 ρ 12 σ 1 s σ 2 s σ 3 s σ 1 s σ 2 s 2 1 ρ 12 2 1 ρ 13 2 ρ 23 ρ 13 ρ 12 2
V 3 s c 3 = π 1 ρ 13 ρ 23 ρ 12 σ 3 s σ 2 s + π 2 ρ 23 ρ 13 ρ 12 σ 1 s σ 3 s + π 3   ρ 12 2 1 σ 1 s σ 2 s σ 1 s σ 2 s σ 3 s 2 1 ρ 12 2 1 ρ 13 2 ρ 23 ρ 13 ρ 12 2
V 1 s j 1 = V s j = μ 1 + 1 2 δ 1 2 δ 1 2 , V s j = V 2 s j 2 = μ 2 + 1 2 δ 2 2 δ 2 2 , λ ˜ s = λ s e x p μ 1 2 2 δ 1 2 + δ 1 2 8 μ 2 2 2 δ 2 2 + δ 2 2 8 , λ ˜ 1 s = λ 1 s   e x p μ 1 2 2 δ 1 2 + δ 1 2 8 , λ ˜ 2 s = λ 2 s   e x p μ 2 2 2 δ 2 2 + δ 2 2 8 , k 1 j 1 = k 2 j 2 = k 1 j = k 2 j = 0 , where π 1 = α s k 1 λ s + λ 1 s r s , π 2 = β s k 2 λ s + λ 2 s r s and π 3 = γ s r s .

4. Pricing Vulnerable Options

In this section, we derive an explicit formula for pricing vulnerable options. Due to the similarity between call and put options, we only provide detailed calculations for vulnerable European call options.
Theorem 3.
Under the proposed model framework, we obtain the pricing formula of European vulnerable call options as
C   S ,   V ,   D ,   X = 0 , T N C * ψ J 1 , J 2 , , J N d J 1 d J 2 d J N
where  ψ J 1 , J 2 , , J N  is the joint distribution of accumulated stay time  J = J 1 , J 2 , , J N  given  X 0 = X  under the equivalent martingale measure  P ¯ . Under the equivalent martingale measure  P ¯ , the characteristic function of  J = J 1 , J 2 , , J N  satisfies the following formula:
E P ¯ e x p i ξ , J X 0 = X = X e x p A T + i G T , I
where  i 2 = 1 ,  ξ = ξ 1 , ξ 2 , , ξ N R N ,  I = 1 , 1 , , 1 R N  and  G  is a N order diagonal matrix generated by  ξ 1 , ξ 2 , , ξ N ;
C * = E P ¯ e 0 T r s d s S T K + I H T 1 + 1 ω H T I H T < 1 F T X = n = 0 m 1 = n m 2 = n Q N ˜ s = n , N ˜ 1 s = n 1 , N ˜ s = n 2 e T   R ¯ T A 1 m 1 ,   m 2 + A 2 m 1 ,   m 2 + A 3 m 1 ,   m 2 + A 4 m 1 ,   m 2 = n = 0 m 1 = n m 2 = n λ ˜ ¯ s T n n ! λ ˜ ¯ 1 s T m 1 n m 1 n ! λ ˜ ¯ 2 s T m 2 n m 2 n ! e T λ ˜ ¯ s + λ ˜ ¯ 1 s + λ ˜ ¯ 2 s + R ¯ T × A 1 m 1 ,   m 2 + A 2 m 1 ,   m 2 + A 3 m 1 ,   m 2 + A 4 m 1 ,   m 2
Proof. 
According to the principle of risk-neutral pricing, under the equivalent martingale measures P ¯ , S 0 = S , D 0 = D , V 0 = V and the initial state of the national economy X 0 = X , the European vulnerable options C   S ,   V ,   D ,   X can be written as follows:
C   S ,   V ,   D ,   X = E P ¯ e x p 0 T r s d s S T K + I V T D T + 1 ω V T D T I V T < D T = E P ¯ e x p 0 T r s d s S T K + I H T 1 + 1 ω H T I H T < 1
where H T = V T D T . Under the equivalent martingale measure P ¯ , the expression of V T and D T are expressed by Equations (22)–(24), and the expression of H T is
log H T H 0 = 0 T 1 2 σ 3 s 2 λ ˜ s k 2 j λ ˜ 2 s k 2 j 2 1 2 σ 2 s 2 d s + 0 T σ 2 s d W ˜ 2 s 0 T σ 3 s d W ˜ 3 s + 0 T Z ˜ s 2 d N ˜ 2 s + 0 T Z ^ s 2 d N ˜ s
The model parameters we propose include the interest rate, jump intensity, and volatility of asset value, which are governed by an observable continuous-time finite-state Markov chain. To facilitate calculations, we need to determine the residence time of each state in the finite-state Markov chain (which can be interpreted as the duration of “good” and “bad” states in the market economy). We assume that J i represents the staying time of the market economy state X t t 0 , T in state i within [0, T] time, we have
R ¯ T = 1 T 0 T r s d s = 1 T i = 0 N r i J i ,   U ¯ 1 T = 1 T 0 T σ 1 s 2 d s = 1 T i = 0 N σ 1 i 2 J i ,
U ¯ 2 T = 1 T 0 T σ 2 s 2 d s = 1 T i = 0 N σ 2 i 2 J i ,   U ¯ 3 T = 1 T 0 T σ 3 s 2 d s = 1 T i = 0 N σ 3 i 2 J i ,
U ¯ 4 T = 1 T 0 T σ 1 s σ 2 s d s = 1 T i = 0 N σ 1 i σ 2 i J i ,   U ¯ 5 T = 1 T 0 T σ 1 s σ 3 s d s = 1 T i = 0 N σ 1 i σ 3 i J i ,
U ¯ 6 T = 1 T 0 T σ 2 s σ 3 s d s = 1 T i = 0 N σ 2 i σ 3 i J i ,   λ ˜ ¯ s = 1 T 0 T λ ˜ s d s = 1 T i = 0 N λ ˜ i J i ,
λ ˜ ¯ 1 s = 1 T 0 T λ ˜ 1 s d s = 1 T i = 0 N λ ˜ 1 i J i ,   λ ˜ ¯ 2 s = 1 T 0 T λ ˜ 2 s d s = 1 T i = 0 N λ ˜ 2 i J i .
We need to calculate the European call vulnerable options C   S ,   V ,   D ,   X under condition S 0 = S , D 0 = D , V 0 = V and X 0 = X , and first calculate the European call vulnerable option under the given conditions F T X , S 0 = S , D 0 = D , V 0 = V and X 0 = X , i.e.,
C * = C S , V , D ,   R ¯ T ,   U ¯ 1 T , U ¯ 2 T , U ¯ 3 T , U ¯ 4 T , U ¯ 5 T , U ¯ 6 T , λ ˜ ¯ s , λ ˜ ¯ 1 s , λ ˜ ¯ 2 s = E P ¯ e 0 T r s d s S T K + I H T 1 + 1 ω H T I H T < 1 F T X
In order to calculate the European call vulnerable option under conditions S 0 = S , D 0 = D , V 0 = V and X 0 = X , we need to know the joint conditional distribution of R ¯ T ,   U ¯ 1 T , U ¯ 2 T , U ¯ 3 T , U ¯ 4 T , U ¯ 5 T ,   U ¯ 6 T , λ ˜ ¯ s , λ ˜ ¯ 1 s and λ ˜ ¯ 2 s under the equivalent martingale measure P ¯ . But R ¯ T ,   U ¯ 1 T , U ¯ 2 T , U ¯ 3 T , U ¯ 4 T , U ¯ 5 T ,   U ¯ 6 T , λ ˜ ¯ s , λ ˜ ¯ 1 s and λ ˜ ¯ 2 s are regulated by the observable continuous-time finite state Markov chain, so the joint conditional distribution of R ¯ T ,   U ¯ 1 T , U ¯ 2 T , U ¯ 3 T , U ¯ 4 T , U ¯ 5 T ,   U ¯ 6 T , λ ˜ ¯ s , λ ˜ ¯ 1 s and λ ˜ ¯ 2 s under the equivalent martingale measure P ¯ is determined by the conditional joint distribution of the cumulative stay time J = J 1 , J 2 , , J N under conditions S 0 = S , D 0 = D , V 0 = V and X 0 = X under the equivalent martingale measure P ¯ . Thus, we have
C   S ,   V ,   D ,   X = 0 , T N C * ψ J 1 , J 2 , , J N d J 1 d J 2 d J N
where ψ J 1 , J 2 , , J N is the joint distribution of cumulative stay time J = J 1 , J 2 , , J N under the equivalent martingale measure P ¯ and X 0 = X . Under the equivalent martingale measure P ¯ , the characteristic function of J = J 1 , J 2 , , J N satisfies the following formula:
E P ¯ e x p i ξ , J X 0 = X = X e x p A T + i G T , I
where i 2 = 1 , ξ = ξ 1 , ξ 2 , , ξ N R N , I = 1 , 1 , , 1 R N and G is an N order diagonal matrix generated by ξ 1 , ξ 2 , , ξ N . See Buffmgton [26] for proof of theoretical properties.
We next calculate C * , given the condition G T n ,   n 1 ,   n 2 = N ˜ s = n , N ˜ 1 s = n 1 , N ˜ 2 s = n 2 and F T X , S T all jumps are m 1 = n 1 + n and H T all jumps are m 2 = n 2 + n , we have
C * = E P ¯ e 0 T r s d s S T K + I H T 1 + 1 ω H T I H T < 1 F T X = n = 0 m 1 = n m 2 = n Q N ˜ s = n , N ˜ 1 s = n 1 , N ˜ s = n 2 E P ¯ { e 0 T r s d s ( S T , m 1 K ) + I H T , m 2 1 + 1 ω H T , m 2 I H T , m 2 < 1 }
By calculation, we have
E P ¯ e 0 T r s d s S T , m 1 K + I H T , m 2 1 + 1 ω H T , m 2 I H T , m 2 < 1 = e T   R ¯ T A 1 m 1 ,   m 2 + A 2 m 1 ,   m 2 + A 3 m 1 ,   m 2 + A 4 m 1 ,   m 2
By observing A 1 A 4 , we find that log S T , m 1 S 0 , log H T , m 2 H 0 is a bivariate normal distribution under the equivalent martingale measure P ¯ . In order to calculate A 1 A 4 , we first obtain the properties of the bivariate normal distribution log S T , m 1 S 0 , log H T , m 2 H 0 as follows:
M 1 m 1 = E P ¯ log S T , m 1 S 0 = 0 T r s 1 2 σ 1 s 2 d s 1 2 m 1 δ 1 2 = T R ¯ T 1 2   U ¯ 1 T 1 2 m 1 δ 1 2
M 2 m 2 = E P ¯ log H T , m 2 H 0 = 0 T 1 2 σ 3 s 2 σ 2 s 2 d s 1 2 m 2 δ 2 2 = T 2   U ¯ 3 T U ¯ 2 T 1 2 m 2 δ 2 2
H 2 m 1 = V a r log S T , m 1 S 0 = 0 T σ 1 s 2 d s + m 1 δ 1 2 = T U ¯ 1 T + m 1 δ 1 2
H 2 m 2 = V a r log H T , m 2 H 0 = T U ¯ 2 T + U ¯ 3 T 2 U ¯ 6 T ρ 23 + m 2 δ 2 2
Cov log S T , m 1 S 0 , log H T , m 2 H 0 = 0 T σ 1 s σ 2 s ρ 12 σ 1 s σ 3 s ρ 13 d s = T U ¯ 4 T ρ 12 + U ¯ 5 T ρ 13
ρ ¯ 1 m 1 ,   m 2 = c o r r log S T , m 1 S 0 , log H T , m 2 H 0 = T U ¯ 4 T ρ 12 + U ¯ 5 T ρ 13 H 2 m 1 H 2 m 2
In order to simplify the calculation, we can rewrite log S T , m 1 S 0 , log H T , m 2 H 0 as
log S T , m 1 S 0 = M 1 m 1 + H 2 m 1 ε 1
log H T , m 2 H 0 = M 2 m 2 + H 2 m 2 ε 2
where ε 1 , ε 2 is a two-dimensional standard normal distribution, and the correlation coefficient is ρ ¯ 1 m 1 ,   m 2 . Next, we have
A 1 m 1 ,   m 2 = S 0 e x p M 1 m 1 + 1 2 H 2 m 1 N 2 a 1 m 1 , a 2 m 1 , m 2 , ρ ¯ 1 m 1 ,   m 2
A 2 m 1 ,   m 2 = K N 2 b 1 m 1 , b 2 m 2 , ρ ¯ 1 m 1 ,   m 2
A 3 m 1 ,   m 2 = S 0 1 ω H 0 e x p { M 1 m 1 + M 2 m 2 + 1 2 H 2 m 1 + H 2 m 2   + ρ ¯ 1 m 1 ,   m 2 H m 1 H m 2 } N 2 c 1 m 1 ,   m 2 , c 2 m 1 ,   m 2 , ρ ¯ 1 m 1 ,   m 2
A 4 m 1 ,   m 2 = K 1 ω H 0 e M 2 m 2 + 1 2 H 2 m 2 N 2 d 1 m 1 ,   m 2 , d 2 m 1 ,   m 2 , ρ ¯ 1 m 1 ,   m 2
Therefore, given the condition F T X , we can obtain:
C * = E P ¯ e 0 T r s d s S T K + I H T 1 + 1 ω H T I H T < 1 F T X = n = 0 m 1 = n m 2 = n Q N ˜ s = n , N ˜ 1 s = n 1 , N ˜ s = n 2 e T   R ¯ T A 1 m 1 ,   m 2 + A 2 m 1 ,   m 2 + A 3 m 1 ,   m 2 + A 4 m 1 ,   m 2 = n = 0 m 1 = n m 2 = n λ ˜ ¯ s T n n ! λ ˜ ¯ 1 s T m 1 n m 1 n ! λ ˜ ¯ 2 s T m 2 n m 2 n ! e T λ ˜ ¯ s + λ ˜ ¯ 1 s + λ ˜ ¯ 2 s + R ¯ T   × A 1 m 1 ,   m 2 + A 2 m 1 ,   m 2 + A 3 m 1 ,   m 2 + A 4 m 1 ,   m 2
The proof of Theorem 3 is complete. □

5. Numerical Simulation

In this section, we present the numerical results of European vulnerable call options prices given by Equations (1)–(3). To facilitate our analysis, we use the set of baseline parameters given in Table 1 and Table 2, and we assume that there are 252 trading days in a year. To facilitate our simulation, we assume that for X t there are only two regimes, i.e., N = 2. The first regime e 1 = 1 , 0 can be interpreted as a good economic state. The second plan e 2 = 0 , 1 can be interpreted as a bad economic state. We refer to relevant literature (such as Han [24]) to choose the initial state of the economy X = X 0 = 1 , 0 , and the generation matrix of the Markov chain process X t
A = p 11 p 12 p 21 p 22 = 0.2 0.8 0.7 0.3
From Figure 1, it can be observed that as the initial price of the underlying asset, S_0, increases, the price of the European fragile option also increases. This is because as the price of the underlying asset rises, the potential profit for the option holder increases, thereby increasing the value of the option.
Additionally, we can observe that under State 1 (good economy), the price of the call option is relatively low, whereas under State 2 (bad economy), the price of the call option is relatively high. This is mainly because in a bad economy, the volatility of the risky asset increases, leading to an increase in the risk value of the option. Therefore, investors need to pay a higher price to purchase options under adverse economic conditions.
Furthermore, when considering a mixed state of the economy, i.e., during a transitional period, the value of the option lies between State 1 and State 2. This is because in a mixed state, the economic conditions are uncertain, and the value of the option is influenced by different economic states. As a result, the price of the option falls between the prices observed in a good economy and a bad economy.
In light of these observations, we propose the switching model to better explain the impact of economic conditions on option valuation. In the switching model, the parameters of the risky asset value model evolve under two different economic states, allowing for a more accurate reflection of the changing economic conditions on option values. By introducing the switching model, our model can more flexibly capture changes in economic conditions, thereby better reflecting the pricing of fragile options in practice.
Now, let us turn to Figure 2. Figure 2 shows that as the initial value of the default barrier (counterparty’s debt), D_0, increases, the price of the European fragile option exhibits a decreasing trend. This is because a larger D_0 increases the risk of credit default, reducing the demand for options from investors and leading to a decline in option prices.
Furthermore, by examining Figure 3, we compare the prices of the European fragile option under four different models with the initial price of the underlying asset, S_0. Compared to the Han [24] model, our proposed model has lower prices. This indicates that our model, which considers the dynamic changes in the default barrier, more accurately reflects the increase in default risk, resulting in lower option values. Compared to the Hui [4] model, our model has higher option prices. This is because our model incorporates a jump process of the risky asset, increasing the prices of options and better capturing the impact of economic conditions on options. Compared to the Wang [22] model, when S_0 is between 8 and 9, our model has similar prices to the Wang [22] model. However, when S_0 is between 9 and 12, our model has higher prices.
Figure 4 shows the price of European Vulnerable options plotted against the initial value of the default barrier (D_0) for four different models. As D_0 increases, the price of European Vulnerable options decreases. Notably, our proposed model exhibits lower prices compared to the Han [24] model, indicating that the uncertainty of the default barrier increases over time, resulting in higher default risk and lower option prices. To provide a more comprehensive explanation of these findings, let us delve into the observed models, starting with the Han [24] and Wang [22] models. In these models, as the fixed default barrier (representing the initial value of the counterparty’s liabilities) increases, the option holder assumes higher credit risk. This suggests that if the default barrier is significantly smaller than the counterparty’s assets, the likelihood of default diminishes, leading to higher option prices. Next, let us consider our proposed model and the Hui [4] model. Both models exhibit a similar rate of decline in option prices, indicating that they share a comparable sensitivity to changes in the default barrier. As D_0 increases, the prices of options in both models decline proportionally. Lastly, we compare our proposed model with the Wang [22] model. When D_0 is small, we find that the Wang [22] model yields lower prices compared to our proposed model. This discrepancy arises because, in such instances, the impact of the jump process on credit risk outweighs the effect of the default barrier. These observations provide a detailed explanation of the relationship between the price of European Vulnerable options and the initial value of the default barrier across different models. As D_0 increases, default risk amplifies, leading to a decrease in option prices. Moreover, the influence of the default barrier and the jump process on option prices and credit risk varies among the different models. These findings contribute to a better understanding and assessment of pricing and risk management for vulnerable options.
Figure 5, Figure 6, Figure 7 and Figure 8 show the impact of basic parameters on vulnerable option prices, including time to maturity, recovery rate, the correlation coefficient between the underlying asset and the counterparty’s assets, and common jump intensity. These figures display the option prices calculated using our proposed model plotted against the initial price of the underlying asset S_0 from 8 to 12 with alternative basic parameters.
As shown in Figure 5, as the option’s time to maturity increases, the option price also increases. This indicates that a longer time to maturity increases the risk of counterparty default and, therefore, increases the value of the option. Figure 6 shows that the option price increases as the recovery rate ω increases. This is because a higher recovery rate results in a higher return at maturity. In Figure 7, we can see that the option price decreases as the correlation coefficient ρ_13 increases. The correlation coefficient reflects the relationship between two assets that move together. This means that as the initial value of the underlying asset S_0 increases, the distance between the three lines also increases. A stronger correlation between the underlying assets and the default barrier means that default risk has a greater impact on vulnerable option prices and results in lower option prices. Figure 8 shows that the option price increases as common jump intensity increases. This result is consistent with our intuitive understanding.

6. Conclusions

This paper discusses the pricing of vulnerable European options using a regime-switching jump-diffusion model. In the structured model, it is assumed that the counterparty’s asset-liability structure (or default barrier) is driven by a geometric Brownian motion under a regime-switching model, and there exists a correlation among all risk assets. By applying the Esscher transform to select an equivalent martingale measure, a pricing formula for vulnerable options is derived. The results of numerical simulations show that the European vulnerable option prices in the proposed model are lower than those in the model assuming a fixed counterparty’s asset–liability structure. This result is intuitive because the counterparty’s liability level is stochastic, increasing the risk of option default and thus lowering the option prices.
This paper mainly investigates the pricing problem of options with random default barriers and credit risk under a regime-switching model. Several research findings are obtained. However, there are still some gaps between the proposed model and the real financial market. Further research can be conducted in the following aspects:
  • More data on OTC option trading can be collected to discuss the parameter estimation issues of the relevant pricing models for vulnerable options, referring to existing literature.
  • As both models in this paper assume fixed volatility of risk assets and risk-free interest rates, further research directions can consider introducing stochastic volatility and risk-free interest rates into the proposed models to make them adhere to more general models. Additionally, with the advancement of computer performance, alternative methods can be considered to numerically solve the option pricing formulas.

Author Contributions

Methodology, X.L.; Software, Z.Z.; Formal analysis, X.L.; Data curation, Z.Z.; Writing—original draft, Z.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (No.71471075), Fundamental Research Funds for the Central University (No.19JNLH09), Innovation Team Project in Guangdong Province, P.R. China (No.2016WCXTD004) and Industry -University- Research Innovation Fund of Science and Technology Development Center of Ministry of Education, P.R. China (No.2019J01017).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare that they have no competing interest.

References

  1. Johnson, H.; Stulz, R. The pricing of options with default risk. J. Financ. 1987, 42, 267–280. [Google Scholar] [CrossRef]
  2. Klein, P. Pricing Black-Scholes options with correlated credit risk. J. Bank. Financ. 1996, 20, 1211–1229. [Google Scholar] [CrossRef]
  3. Klein, P.; Inglis, M. Pricing vulnerable European options when the option’s payoff can increase the risk of financial distress. J. Bank. Financ. 2001, 25, 993–1012. [Google Scholar] [CrossRef]
  4. Hui, C.H.; Lo, C.F.; Ku, K.C. Pricing vulnerable European options with stochastic default barriers. IMA J. Manag. Math. 2007, 18, 315–329. [Google Scholar] [CrossRef]
  5. Wang, X. Pricing vulnerable options with stochastic default barriers. Financ. Res. Lett. 2016, 19, 305–313. [Google Scholar] [CrossRef]
  6. Xu, W.; Xu, W.; Li, H.; Xiao, W. A jump-diffusion approach to modelling vulnerable option pricing. Financ. Res. Lett. 2012, 9, 48–56. [Google Scholar] [CrossRef]
  7. Tian, L.H.; Wang, G.Y.; Wang, X.C. Pricing vulnerable options with correlated credit risk under jump-diffusion processes. J. Futures Mark. 2014, 34, 957–979. [Google Scholar] [CrossRef]
  8. Ma, Y.; Shrestha, K.; Xu, W. Pricing vulnerable options with jump clustering. J. Futures Mark. 2017, 37, 1155–1178. [Google Scholar] [CrossRef]
  9. Zhou, Q.; Wang, Q.; Wu, W.X. Pricing vulnerable options with variable default boundary under jump-diffusion processes. Adv. Differ. Equ. 2018, 2018, 465. [Google Scholar] [CrossRef]
  10. Wang, X. Differences in the prices of vulnerable options with different counterparties. J. Futures Mark. 2017, 37, 148–163. [Google Scholar] [CrossRef]
  11. Liu, S.I.; Liu, Y.C. Pricing Vulnerable Options under Stochastic Asset and Liability. Prog. Appl. Math. 2011, 1, 12–29. [Google Scholar] [CrossRef]
  12. Yang, S.J.; Lee, M.K.; Kim, J.H. Pricing vulnerable options under a stochastic volatility model. Appl. Math. Lett. 2014, 34, 7–12. [Google Scholar] [CrossRef]
  13. Lee, M.K.; Yang, S.J.; Kim, J.H. A closed form solution for vulnerable options with Heston’s stochastic volatility. Chaos Solitons Fractals 2016, 86, 23–27. [Google Scholar] [CrossRef]
  14. Wang, G.; Wang, X.; Zhou, K. Pricing vulnerable options with stochastic volatility. Phys. A Stat. Mech. Appl. 2017, 485, 91–103. [Google Scholar] [CrossRef]
  15. Yoon, J.H.; Kim, J.H. The pricing of vulnerable options with double Mellin transforms. J. Math. Anal. Appl. 2015, 422, 838–857. [Google Scholar] [CrossRef]
  16. Huang, S.; Guo, X. Valuation of European-style vulnerable options under the non-affine stochastic volatility and double exponential jump. Chaos Solitons Fractals 2022, 158, 112003. [Google Scholar] [CrossRef]
  17. Capponi, A.; Figueroa-López, J.E.; Nisen, J. Pricing and Semimartingale Representations of Vulnerable Contingent Claims in Regime-Switching Markets. Math. Financ. 2014, 24, 250–288. [Google Scholar] [CrossRef]
  18. Elliott, R.J.; Chan, L.; Siu, T.K. Option pricing and Esscher transform under regime switching. Ann. Financ. 2005, 1, 423–432. [Google Scholar] [CrossRef]
  19. Elliott, R.J.; Osakwe, C.J.U. Option pricing for pure jump processes with Markov switching compensators. Financ. Stoch. 2006, 10, 250–275. [Google Scholar] [CrossRef]
  20. Wang, W.; Jin, Z.; Qian, L.; Su, X. Local risk minimization for vulnerable European contingent claims on nontradable assets under regime switching models. Stoch. Anal. Appl. 2016, 34, 662–678. [Google Scholar] [CrossRef]
  21. Yang, Q.Q.; Ching, W.K.; He, W.; Siu, T.K. Pricing vulnerable options under a Markov-modulated jump-diffusion model with fire sales. J. Ind. Manag. Optim. 2019, 15, 293–318. [Google Scholar] [CrossRef]
  22. Wang, W.; Wang, W.S. Pricing vulnerable options under a Markov-modulated regime-switching model. Commun. Stat.-Theory Methods 2010, 39, 3421–3433. [Google Scholar] [CrossRef]
  23. Niu, H.W.; Wang, D.C. Pricing vulnerable options with correlated jump-diffusion processes depending on various states of the economy. Quant. Financ. 2016, 16, 1129–1145. [Google Scholar] [CrossRef]
  24. Han, M.; Song, X.; Niu, H. Pricing Vulnerable Options with Market Prices of Common Jump Risks under Regime-Switching Models. Discret. Dyn. Nat. Soc. 2018, 2018, 8545841. [Google Scholar] [CrossRef]
  25. Buffington, J.; Elliott, R.J. Regime switching and European options. In Stochastic Theory and Control-Proceedings of a Workshop held in Lawrence, Kansas, Volume 280 of Lecture Notes in Control and Information Sciences; Springer: Berlin/Heidelberg, Germany, 2002; pp. 73–82. [Google Scholar]
  26. Damircheli, D.; Razzaghi, M.; Bastani, A.F. A de-singularized meshfree approach to default probability estimation under a regime-switching synchronous-jump tempered stable Lévy model. Eng. Anal. Bound. Elem. 2023, 150, 364–373. [Google Scholar] [CrossRef]
  27. Elliott, R.J.; Aggoun, L.; Moore, J.B. Hidden Markov Models: Estimation and Control; Springer: New York, NY, USA, 1994. [Google Scholar]
  28. Delbaen, F.; Schachermayer, W. A general version of the fundamental theorem of asset pricing. Math. Ann. 1994, 300, 463–520. [Google Scholar] [CrossRef]
  29. Brémaud, P. Point Processes and Queues: Martingale Dynamics; Springer: New York, NY, USA, 1981. [Google Scholar]
Figure 1. Vulnerable European call option price against the initial price of the underlying assets S. The solid, dashed, and dot-dashed lines correspond to State 1, the proposed model, and State 2, respectively.
Figure 1. Vulnerable European call option price against the initial price of the underlying assets S. The solid, dashed, and dot-dashed lines correspond to State 1, the proposed model, and State 2, respectively.
Mathematics 11 04155 g001
Figure 2. Vulnerable European call option price against the initial price of the default barrier D. The solid, dashed, and dot-dashed lines correspond to State 1, the proposed model, and State 2, respectively.
Figure 2. Vulnerable European call option price against the initial price of the default barrier D. The solid, dashed, and dot-dashed lines correspond to State 1, the proposed model, and State 2, respectively.
Mathematics 11 04155 g002
Figure 3. Vulnerable European call option price against the initial price of the underlying assets S. The solid, dashed, dotted, and dot-dashed lines correspond to Hui’s [4] model, the proposed model, Han’s [22] model, and Wang’s [20] model, respectively.
Figure 3. Vulnerable European call option price against the initial price of the underlying assets S. The solid, dashed, dotted, and dot-dashed lines correspond to Hui’s [4] model, the proposed model, Han’s [22] model, and Wang’s [20] model, respectively.
Mathematics 11 04155 g003
Figure 4. Vulnerable European call option price against the initial price of the default barrier D. The solid, dashed, dotted, and dot-dashed lines correspond to Wang’s [22] model, the proposed model, Han’s [24] model, and Hui’s [4] model, respectively.
Figure 4. Vulnerable European call option price against the initial price of the default barrier D. The solid, dashed, dotted, and dot-dashed lines correspond to Wang’s [22] model, the proposed model, Han’s [24] model, and Hui’s [4] model, respectively.
Mathematics 11 04155 g004
Figure 5. The option price with different maturity T for ω = 0.5 , λ 1 = 1 , λ 2 = 2 and ρ 13 = 0.5 . The solid, dashed, and dot-dashed lines correspond to T = 0.5, T = 1, and T = 1.5, respectively.
Figure 5. The option price with different maturity T for ω = 0.5 , λ 1 = 1 , λ 2 = 2 and ρ 13 = 0.5 . The solid, dashed, and dot-dashed lines correspond to T = 0.5, T = 1, and T = 1.5, respectively.
Mathematics 11 04155 g005
Figure 6. The option price with different recovery rate ω for T = 1, λ 1 = 1 , λ 2 = 2 and ρ 13 = 0.5 . The solid, dashed, dotted, and dot-dashed lines correspond to ω = 0.3 , ω = 0.5 , ω = 1 , and ω = 0.8 , respectively.
Figure 6. The option price with different recovery rate ω for T = 1, λ 1 = 1 , λ 2 = 2 and ρ 13 = 0.5 . The solid, dashed, dotted, and dot-dashed lines correspond to ω = 0.3 , ω = 0.5 , ω = 1 , and ω = 0.8 , respectively.
Mathematics 11 04155 g006
Figure 7. The option price with a correlation coefficient between the underlying asset and the assets of the counterparty for ω = 0.5 , λ 1 = 1 , λ 2 = 2 and T = 1. The solid, dashed, and dot-dashed lines correspond to ρ 13 = 0.9 , ρ 13 = 0.5 , and ρ 13 = 1 , respectively.
Figure 7. The option price with a correlation coefficient between the underlying asset and the assets of the counterparty for ω = 0.5 , λ 1 = 1 , λ 2 = 2 and T = 1. The solid, dashed, and dot-dashed lines correspond to ρ 13 = 0.9 , ρ 13 = 0.5 , and ρ 13 = 1 , respectively.
Mathematics 11 04155 g007
Figure 8. The option price with different common jump intensity λ for T = 1, ω = 0.5 , and ρ 13 = 0.5 . The solid, dashed, and dot-dashed lines correspond to λ 1 ,   λ 2 = 0.5 , 1 , λ 1 ,   λ 2 = 1 , 2 , and λ 1 ,   λ 2 = 2 , 4 , respectively.
Figure 8. The option price with different common jump intensity λ for T = 1, ω = 0.5 , and ρ 13 = 0.5 . The solid, dashed, and dot-dashed lines correspond to λ 1 ,   λ 2 = 0.5 , 1 , λ 1 ,   λ 2 = 1 , 2 , and λ 1 ,   λ 2 = 2 , 4 , respectively.
Mathematics 11 04155 g008
Table 1. Parameter values of vulnerable European call options in the base case.
Table 1. Parameter values of vulnerable European call options in the base case.
ParameterValueParameterValue
Strike priceK = 10Mean jump size of V μ 1 = 0
Time to maturityT = 1Mean jump size of S μ 2 = 0
Initial price of S S 0 = 10 correlation coefficient ρ 12 = 0.7
Initial price of D D 0 = 10 correlation coefficient ρ 23 = 0.6
Initial price of V V 0 = 10 correlation coefficient ρ 13 = 0.5
Volatility of jump size of V δ 2 = 0.1 Volatility of jump size of S δ 1 = 0.1
Table 2. Parameter values of vulnerable European call options in the base case.
Table 2. Parameter values of vulnerable European call options in the base case.
ParameterValue in State e 1 Value in State e 2
r r 1 = 0.02 r 1 = 0.02
S’s Volatility σ S 1 = 0.3 σ S 2 = 0.4
V’s Volatility σ V 1 = 0.3 σ V 2 = 0.4
D’s Volatility σ D 1 = 0.3 σ D 2 = 0.4
common jump intensity λ 1 = 1 λ 2 = 2
S jump intensity λ S 1 = 1 λ S 2 = 2
V jump intensity λ V 1 = 1 λ V 2 = 2
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Liu, X.; Zhang, Z. Pricing European Vulnerable Options with Jumps and Stochastic Default Obstacles Barrier under Regime Switching. Mathematics 2023, 11, 4155. https://doi.org/10.3390/math11194155

AMA Style

Liu X, Zhang Z. Pricing European Vulnerable Options with Jumps and Stochastic Default Obstacles Barrier under Regime Switching. Mathematics. 2023; 11(19):4155. https://doi.org/10.3390/math11194155

Chicago/Turabian Style

Liu, Xiangdong, and Zanbin Zhang. 2023. "Pricing European Vulnerable Options with Jumps and Stochastic Default Obstacles Barrier under Regime Switching" Mathematics 11, no. 19: 4155. https://doi.org/10.3390/math11194155

APA Style

Liu, X., & Zhang, Z. (2023). Pricing European Vulnerable Options with Jumps and Stochastic Default Obstacles Barrier under Regime Switching. Mathematics, 11(19), 4155. https://doi.org/10.3390/math11194155

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop