Basket Credit Default Swap Pricing with Two Defaultable Counterparties

In this paper, we study the basket CDS pricing with two defaultable counterparties based on the reduced-form model. ,e default jump intensities of the reference firms and counterparties are all assumed to follow the mean-reverting constant elasticity of variance (CEV) processes. Taking the Vasicek process which is a special case of CEV process as an example, the approximate analytic solutions of the joint survival probability density, the probability densities of the first default and the first two defaults can be solved by using PDE method. In addition, we also extend the Vasciek process to the Vasciek process with cojumps and obtain the approximate closed-form solutions of the relevant default probability densities. ,en with the expressions of the probability densities, we can get the formula of the basket CDS price with two defaultable counterparties. In the numerical analysis, we do sensitivity analysis and compare the basket CDS prices under our model with that with only one defaultable counterparty. ,e numerical results show that our model can be applied into practice.


Introduction
Credit derivatives have been widely used by market participants to manage and hedge credit risks. One of the most popular credit derivatives is the credit default swap (CDS).
ere are mainly two kinds of CDS contracts: single-name CDS and basket CDS. e difference between the two contracts is the number of the reference entities. Since the bankruptcy of the Worldcom Company and the Enron Corporation, people have paid more and more attention to the default probability of large companies. e subprime crisis has made people realize that the correlated default risk plays an important role in the pricing of a basket CDS. With the development of global economic integration, enterprises are more closely related, so the study of pricing the basket CDS has attracted more and more researchers. ere are two common models to study the pricing of credit derivatives in the literature, namely the structural models introduced by Black and Scholes [1] and Merton [2], and the reduced-form intensity-based models pioneered by Jarrow and Turnbull [3].
In the structural model, the default of a firm is deemed to be triggered when the firm value falls below the liability level. Black and Cox [4] proposed the first passage model in which the default time was assumed to be the first time that the firm value broke down the constant barrier. Based on Merton [2] model and  model, Gökgöz et al. [5] studied the evaluation of a single-name CDS via the discounted cash flow method. Chen and He [6] proposed the multi-scale stochastic volatility (SV) model to price a CDS. Based on the structural model and introducing the concept of fuzziness, Wu et al. [7] proposed a new double exponential jump diffusion model with fuzziness for CDS pricing. e reduced-form intensity-based model was introduced by Jarrow and Turnbull [3], in which the Poisson process was used to describe the exogenous default occurrence. Lando [8] proposed a Cox process to describe the default intensity and assumed that the risk-free interest rate satisfied the Vasicek model. Malherbe [9] applied a Poission process to describe the default intensity. He assumed that the default intensities were constant between defaults, but could jump at the times of defaults. Herbertsson and Rootzén [10] used the matrixanalytic method to derive a closed-form expression for a basket CDS. Zheng and Jiang [11] used the total hazard construction method to derive an analytic formula for the joint distribution of default times.
e key of basket CDS pricing is to obtain the relevant default probability of multiple assets.
ere are mainly three methods in the literature to model the default risk of multiple assets: copula function, conditional independence model and contagion model. Using the copula approach, one can derive a joint distribution of default times by combining the marginal distributions of default times. Li [12] studied the default correlations among companies using CreditMetrics model with copula functions. Crépey and Jeanblanc [13] studied CDS pricing with counterparty risk under a Markov chain copula model. Harb and Louhichi [14] used the mixture copula to price the basket CDS with counterparty risk. In the conditional independence model, Finger [15] assumed that the common macro factor affected the default times of all the assets in the portfolio, while the default intensities were independent with each other. Based on the reduced-form default intensity model, Kijima and Muromachi [16] studied the pricing of basket CDS of the firstdefault type and obtained an analytic formula for the basket CDS price. Kijima and Muromachi [17] considered k th -to-default basket CDS pricing, but they did not give the explicit solution. White [7] presented a new model for valuing a CDS contract which was affected by multiple credit risks. ey showed that the default dependency had a significant impact on CDS pricing. Davis and Lo [18] firstly employed the contagion model to describe the default risk of the stock portfolio. Based on the contagion model, Jarrow and Yu [19] introduced the concept of the counterparty risk and illustrated the effect of the counterparty risk on CDS price. Hao and Wang [20] studied the CDS pricing with the contagious risk under the fractional Vasicek interest rate model. Dong and Wang [21] assumed the intensities of the default times were driven by macro-economy described by a Markov regimeswitching model. Gu and Liu [22] established the attenuation model for the contagious risk and derived the pricing formula of CDS in the fractional dimension environment. Huang and Song [23] priced the basket CDS with counterparty risk under a multi-name contagion model.
In this paper, we study the basket CDS pricing with two defaultable counterparties based on the reduced-form model. When the reference asset is a basket of assets and if there are positive correlations among assets, the default probability may be high. In this case, the CDS sellers are willing to sign the CDS contract together to share the default risk. e investors wonder whether more counterparties can reduce the default risk. e main contributions of this paper are as follows. (1) To the best of our knowledge, there is no basket CDS contract with two defaultable counterparties traded in the market yet. However, this kind of contracts can be applied into practice when the time is ripe. At present, it is meaningful to carry out the theoretical research. (2) e default jump intensities of the reference entities and counterparties are all assumed to follow the mean-reverting constant elasticity of variance (CEV) processes. e CEV process is a general process that contains the Vasicek process, CIR process and geometric Brownian motion. Using PDE method, we obtain three PDEs for the joint survival probability density, the probability density of the first default and the probability density of the first two defaults among the reference entities and counterparties. (3) Taking the Vasicek process which is a special case of CEV process as an example, the approximate analytic solutions of the relevant default probability densities can be solved from PDEs. In addition, we also extend the Vasciek process to the Vasciek process with cojumps and obtain the approximate closed-form solutions of the relevant default probability densities.
en with the expressions of the probability densities, we can get the formula of the basket CDS price with two defaultable counterparties. In the numerical analysis, we find that the CDS buyer pay more for the basket CDS contract with two defautable counterparties and there will be almost no price difference if the number of reference assets is large enough. e numerical results show that our model can be applied into practice. It is worthy of note that our model can be extended to the jump model with stochastic volatility. For the introduction of this model, readers can refer to He and Lin [24,25], He and Chen [26,27]. Stochastic volatility model can describe the phenomenon of volatility clustering of default intensity, but the derivation of basket CDS price is quite difficult. When the volatilities of the default intensity of the reference assets and two counterparties are stochastic, the number of state variables will increase which makes the calculation more difficult. e solution for the PDEs satisfied by the relevant default probability densities will not necessarily have analytical solutions. If so, the CDS price can be solved by Monte Carlo simulation and other numerical methods. e article is organized as follows. In Section 2, we assume the default jump intensities of the reference firms and counterparties follow the mean-reverting CEV processes. We obtain three PDEs for the joint survival probability density, the probability density of the first default and the probability density of the first two defaults. In Section 3, we determine the approximate closed-form solutions for relevant default probability densities under the Vasicek processes. In Section 4, we obtain the approximate closedform solutions for relevant default probability densities under the Vasicek processes with cojumps. In Section 5, we derive the formula of the basket CDS price with two defaultable counterparties. In Section 6, we do sensitivity analysis under our model. We compare the price differences of the basket CDS with two defaultable counterparties and that with only one defaultable counterparty. Finally, we offer concluding remarks in Section 7.

Default Probability Density under Reducedform Intensity Model
e traditional basket CDS contract is usually signed with only one credit protection seller. e disadvantage of this kind of contract is that when the credit protection seller defaults, the credit protection buyer is likely to lose the CDS fee or not be compensated. erefore, we consider the basket CDS pricing with two credit protection sellers. If the reference assets in the basket do not default, the credit protection buyer will pay the CDS fee to the two credit protection sellers at the same time. If the default of the reference assets occurs, the sellers will compensate the buyer. In addition, in order to make our model more general and closer to the real market, we assume that two credit protection sellers may also default.
Let T > 0 be a finite time horizon and fix a probability space (Ω, F, P), the probability measure P is the risk-neutral measure. e canonical filtration generated by the underlying stochastic structure is denoted by F t , which defines the information available at each time. e conditional probability measure given F t is denoted by P t and the associated conditional expectation operator is E t . Let default time τ be a stopping time associated with the filtration F t . For sufficiently small Δt ≥ 0, λ(t) is an intensity process for τ if satisfied Suppose λ(t) is a (n + 2)-dimensional Markov process of c� adl� ag state variables drawn from space V ⊂ R n+2 . We assume the reference asset is a basket of bonds F B i (t), i � 1, 2, . . . , n. which are issued by different companies. Each bond may default with the default intensity λ i (t), i � 1, 2, . . . , n . ere are two sellers named F C and F D who will both compensate if one of the assets in the basket defaults. We assume the seller F C has a stochastic default intensity of λ n+1 (t) and seller F D has a stochastic default intensity of λ n+2 (t). All the default intensities λ i (t), i � 1, 2, . . . , n, n + 1, n + 2 are assumed to follow the mean-reverting constant elasticity of variance (CEV) processes where a i , b i , σ i , β are positive constants. Respectively, a i is the mean-reverting rate. b i represents the long-term level of jump intensity. σ i is the volatility of the jump intensity. β can be interpreted as the elasticity. Each W i (t) is a standard Brownian motion. dW i (t)dW j (t) � ρ ij dt for i ≠ j and ρ ij � 1 for i � j. e CEV process in (2) can include some existing processes as special cases. (1) If β � 0, process λ i (t) is the Vasicek process. (2) If β � 1/2, process λ i (t) is the CIR process.
(3) If β � 1 and b i � 0, process λ i (t) is the geometric Brownian motion. Let τ 1 , . . . τ n denote the default times of reference assets F B i (i � 1, 2, . . . n). τ n+1 and τ n+2 represent the default time of counterparty F C and F D respectively throughout this article. Given F T , the default times τ j are conditionally independent. e initial time and expiration date are represented by t and T (0 ≤ t ≤ T). Before we price the basket CDS, we need some conclusions about the default probability densities as shown in eorem 2.1-2.3. eorem 2.1. (Joint survival probability density under the mean-reverting CEV model) Denote λ � λ i , i � 1, 2, . . . , n + 2 , if all the reference assets F B i (i � 1, 2, . . . , n) and two counterparties do not default until time s(t ≤ s ≤ T), the joint survival probability density P(t, λ; s) satisfies the following PDE: Proof. If no default events happen, the CDS buyer will pay the CDS fee continuously until the expiration date T. e conditional independence means the joint survival probability at time s(t ≤ s ≤ T) can be given by . We denote the probability density P(t, λ; s) as With Feynman-Kac theorem, we can get the PDE (3) that P(t, λ; s) satisfies. eorem 2.2 ( e probability density of the first default for the i th company under the meanreverting CEV model). Among the reference assets F B i (i � 1, 2, . . . , n) issued by different companies and two counterparties, if the i th company firstly defaults at time τ i (s ≤ τ i ≤ s + ds), we denote this default probability density at time s to be q i (t, λ; s). en the probability density q i (t, λ; s) satisfies the following PDE □ Discrete Dynamics in Nature and Society Proof. Among n + 2 companies, the probability of the first default for the i th company at time τ i (s ≤ τ i ≤ s + ds) can be given by Because , the probability density of the first default for the i th company at time s is With Feynman-Kac theorem, we can get the PDE (6) that the probability density q i (t, λ; s) satisfies. eorem 2.3. ( e probability density of the first two defaults for the i th and j th companies under the mean-reverting CEV model) Among the reference assets F B i (i � 1, 2, . . . , n) issued by different companies and two counterparties, if the i th company is the first to default at time τ i (t ≤ τ i ≤ s) and the j th company is the second to default at time . . , n + 2 and then the probability density q ij of the default event has a solution form as follows where Q(t, λ; s) satisfies the following PDE □ Proof. Among the n + 2 defaultable companies (including the n reference assets and two counterparties), the i th company is the first to default at time τ i (u ≤ τ i ≤ u + du) and the j th company is the second to default at time e cumulative default process is denoted by en the probability of the default event is So 4 Discrete Dynamics in Nature and Society Denote According to (15), we have With Feynman-Kac theorem, we can get the PDE (10) that Q(t, λ; s) satisfies.

Basket CDS Pricing with the Vasicek Processes
For PDEs (3) (6) (10), there are no analytical solutions generally. e analytic solutions exist only when β is a special value. In this paper, we aim to take the Vasicek process (i e. β � 0) as an example to derive the analytic price for the basket CDS. How to get the general solution under CEV process, CIR process and geometric Brownian motion are the problems we need to solve in the future. Although the Vasicek process may take negative values which is economically unacceptable, the probability of such a pathology to arise is very small (please refer to Rogers [28], Hull and White [29]). Consequently, the inconvenience linked to the Vasicek process and its extensions can be neglected facing to the benefit they could bring. eorem 3.1. (Joint survival probability density under the Vasicek model) Denote λ � λ i , i � 1, 2, . . . , n + 2 , if all the reference assets F B i (i � 1, 2, . . . , n) and two counterparties do not default until time s(t ≤ s ≤ T), the joint survival probability density P(t, λ; s) has a closed-form solution as follows where Proof. According to eorem 2.1, According to ∅ ksendal (2003), P(t, λ; s) has a solution with the following form Substitute the above formula into Equation (21) to get two ODEs Solve the above ODEs and obtain (19) and (20). eorem 3.2. ( e probability density of the first default for the i th company under the Vasicek model) Among the reference assets F B i (i � 1, 2, . . . , n) issued by different companies and two counterparties, if the i th company firstly defaults at time τ i (s ≤ τ i ≤ s + ds), we denote this default probability density at time s to be q i (t, λ; s). q i (t, λ; s) has a closed-form solution as follows Discrete Dynamics in Nature and Society A(t; s) and B k (t; s) are expressed as in (19) and (20). Proof. According to eorem 2.2, q i (t, λ; s) satisfies the following PDE if λ i is a Vasicek process (i e. β � 0) According to ∅ ksendal (2003), the probability density has a solution with the following form A(t; s) and B k (t; s) have been solved by (23). Substitute the above formula into (27) to get two ODEs as follows Solve the above ODEs and obtain (25) and (26). eorem 3.3. ( e probability density of the first two defaults for the i th and j th companies) Among the reference assets F B i (i � 1, 2, . . . , n) issued by different companies and two counterparties, if the i th asset is the first to default at time τ i (t ≤ τ i ≤ s) and the j th asset is the second to default at time . . n + 2 and then the probability density q ij of the default event has a closedform solution as follows Q(t, λ; s) has a form like where 6 Discrete Dynamics in Nature and Society According to ∅ ksendal (2003), Q(t, λ; s) has a solution with the following form where A(t; s) and B k (t; s) satisfy the following ODEs which are similar to (23).

Basket CDS Pricing with Cojumps
In this section, we will extend the Vasicek processes to the jump processes. We assume there exist simultaneous jumps called cojumps among all the companies when an extreme event occurs. All the default intensities λ i (t), i � 1, 2, . . . , n, n + 1, n + 2 are assumed to follow the Vasicek processes with cojumps where N(t) is a Poisson counter with constant intensity λ J . prob(dN(t) � 1) � λ J dt and prob(dN(t) � 0) � 1 − λ J dt.
ε i is the percentage jump size (conditional on a jump occurring). Before we price the basket CDS, we will give some results as shown in eorem 4.1-4.3 about the default probability densities under the jump models (40). eorem 4.1. (Joint survival probability density) If all the reference assets F B i (i � 1, 2, . . . , n) and two counterparties do not default until time s(t ≤ s ≤ T), the joint survival probability density P ′ (t, λ; s) will have a closed-form solution as follows: where, Discrete Dynamics in Nature and Society Proof. Similar with eorem 3.1, according to Feynman-Kac theorem, the joint survival probability P ′ (t, λ; s) at time s(t ≤ s ≤ T) satisfies the following PDE According to ∅ ksendal (2003), P ′ (t, λ; s) has a solution with the following form Substitute the above formula into (44) to obtain With the approximate formula We substitute (47) into (46) to obtain two ODEs Solve the above ODEs and obtain (42) and (43). eorem 4.2. ( e probability density of the first default for the i th company) Among the reference assets F B i (i � 1, 2, . . . , n) issued by different companies and two counterparties, if the i th asset firstly defaults at time τ i (s ≤ τ i ≤ s + ds), the default probability density q i ′ (t, λ; s) at time s has a closed-form solution as follows where And A ′ (t;s) and B k (t;s) are expressed as in (42) According to ∅ ksendal(2003), q i ′ (t, λ; s) has a solution with the following form Substitute the above formula into (52), we have With (42) and (43), we can obtain two ODEs Solve the above ODEs and obtain (50) and (51). eorem 4.3 (the probability density of the first two defaults for the i th and j th companies). Among the reference assets F B i (i � 1, 2, . . . , n) issued by different companies and two counterparties, if the i th asset is the first to default at time τ i (t ≤ τ i ≤ s) and the j th asset is the second to default at time τ j (s ≤ τ j ≤ s + ds) with the constrain that τ i ≤ T, τ i ≤ τ j ≤ T, then the probability density q ij ′ of the default event has a closed-form solution Q ′ (t, λ; s) has a form like where Discrete Dynamics in Nature and Society □ Proof. Similar with eorem 3.3, the probability density q ij ′ (t, λ; s) has a form as follows where Denote ε � ε j , j � 1, 2, . . . , i − 1, i + 1, . . . , n + 2 and with Feynman-Kac theorem, Q ′ (t, λ; s) satisfies the following PDE According to ∅ ksendal(2003), Q ′ (t, λ; s) has a solution with the following form Solve the above ODEs and obtain (58) and (59).

Basket CDS Price
In this section, we will discuss the pricing of the basket CDS with two defaultable counterparties under the Vasicek model and the jump model in (40) respectively. Firstly, we consider the basket CDS price under the Vasicek model using eorems 3.1-3.3 as an example. We assume the credit protection buyer F A holds a basket of n reference assets F B i (i � 1, 2, . . . , n). At initial time t, F A buys a basket CDS contact with the credit protection sellers F C and F D . e maturity of the basket CDS is T. During time t to T, the CDS buyer F A will pay the CDS fee continuously to the CDS sellers F C and F D until T if no defaults happen. Once one of the assets F B i (i � 1, 2, . . . n) defaults first, the CDS sellers F C and F D will both compensate to the CDS buyer F A . However, if the CDS seller F C or F D defaults first, the CDS buyer F A will stop paying premiums to F C or F D and not receive any compensations from F C or F D . Note that when one counterparty defaults, the CDS buyer will still pay the premium to the other counterparty that does not default. e default events include three types of defaults, namely, the default of reference asset F B i , the default of counterparty F C and the default of counterparty F D . at is, we need to consider the default order of F B i , F C and F D . Now we analysis all the possible default events once the basket CDS contract becomes effective from time t : For any i(i � 1, 2, . . . n), reference asset F B i is the first to default at time τ i . Next we will discuss how to compute the basket CDS price that the credit protection buyer F A pay to counterparty F C . We assume the CDS costs are continuously paid by F A and denote the cost rate to be W 1 . Firstly, we need to analyze the present value at time t of the basket CDS costs received by counterparty F C under different situations. Situation 1. According to eorem 3.2, the default probability density under this situation is q i (t, λ; s). e present value of the CDS costs received by coun- Situation 2. Correspondingly, according to the notation stated in eorem 3.3, (1) For any i(i � 1, 2, . . . n), the probability density of the default event is q n+1,i (t, λ; s) (2) e probability density of the default event is q n+1,n+2 (t, λ; s) (3) e probability density of the default event is q n+1 (t, λ; s) − n i�1 q n+1,i (t, λ; s) − q n+1,n+2 (t, λ; s) e present value of the total CDS costs received by counterparty  1, 2, . . . n), the probability density of the default event is q n+2,i (t, λ; s) (2) e probability density of the default event is q n+2,n+1 (t, λ; s) (3) e probability density of the default event is q n+2 (t, λ; s) − n i�1 q n+2,i (t, λ; s) − q n+2,n+1 (t, λ; s) e present value of the total CDS costs received by counterparty Situation 4. According to eorem 3.1, the joint survival probability density is P(t, λ; s).
e present value of the CDS costs received by counterparty F C is W 1 T t e − r(s− t) P(t, λ; s)ds. en we will analyze the present values of compensations paid by counterparty F C under different situations. Denote R to be the recovery rate and L i to be the face value of the reference assets F B i (i � 1, 2 . . . n). We assume that the counterparty F C and counterparty F D prefer to share default risk. Under 1, counterparty F C will compensate (1/2)L i (1 − R) T t e − r(s− t) q i (t, λ; s)ds to the basket CDS buyer F A . Under 2, counterparty F C is the first to default, so there will be no compensations paid by counterparty F C . Under 3, only when the counterparty F D defaults firstly and F B i defaults secondly, counterparty F C will compensate. e present values of the compensations paid by counterparty T t e − r(s− t) q n+2,i (t, λ; s)ds. Otherwise, there will be no compensations paid by counterparty F C . Under Situation 4, there will be no compensations paid by counterparty F C if no defaults happen from time t to T.
Finally, according to the no-arbitrage pricing principal, the present value of the total CDS costs received by counterparty F C should be equal to the present value of the total compensations paid by counterparty F C . us we have So the CDS price W 1 the buyer F A paid to counterparty F C is Discrete Dynamics in Nature and Society T t e −r(s−t) q i (t, λ; s)ds + T t e −r(s−t) P(t, λ; s)ds .

(68)
Similarly, the CDS price W 2 paid to counterparty F D should satisfy the following equation (69) When there exist cojumps, using the similar method and the results in eorems 4.1-4.3, the CDS price the buyer F A paid to counterparty F C and F D can be given respectively: With the closed-form solution, the sensitivity of CDS price to the initial default intensity can be measured by partial derivative. Denote f(t, λ; s), g(t, λ; s) as follows, And we can have analytical partial derivatives where And Other derivatives (zW 2 /zλ j ), (zW 1 ′ /zλ j ), (zW 2 ′ /zλ j ) can be obtained in a similar way.

Numerical Analysis
In this section, we will do some numerical analysis to show the impacts of main parameters on CDS prices. In order to verify the correctness of our formulas for CDS prices, we do Monte Carlo simulations. Due to the similar structures of CDS prices, we compute W 1 ′ under the model in Section 4 using formula (71) as an example. In order to obtain W 1 ′ , we need to know the relevant probability densities P ′ (t, λ; s), q i ′ (t, λ; s), q n+2 ′ (t, λ; s). Without the approximate closed-form solutions for P ′ (t, λ; s), q i ′ (t, λ; s), q n+2 ′ (t, λ; s), they can also be numerically solved, which can serve as the benchmark solution to the original problem. us we use the solution for CDS price calculated by Monte Carlo simulation as the benchmark solution. For i � 1, 2, . . . , n, n + 1, n + 2, we use the following formulas in a discrete time form to simulate the paths of the default intensities for reference assets and counterparties.
where ξ i follows the multivariate standard normal distribution with an (n + 2) × (n + 2) correlation matrix. N(t) is a Poisson processes with constant intensities λ J . We assume the jump sizes ε j i are all constants for simplicity. By simulating the default intensities, we can calculate the relevant probability densities, so as to obtain the price of CDS. When the jump intensity λ J is equal to zero, the simulation results can be used to calculate the CDS price of formula (68) under the continuous model in Section 3. In the procedure of the Monte Carlo simulations, we set the number of time step to be 100, and the number of simulations to be 100000. In Table 1, we show the relative percentage differences are all less than 1%. erefore, formula (71) is effective to compute CDS prices. Formula (68) can be also proved to be effective using the same method (the results here are omitted).
Our paper mainly considers two kinds of models, that is, the Vasicek model and the Vasicek model with cojumps (hereinafter referred to as VCJ model). Next, we will use the numerical solution to do some numerical analysis under the two kinds of models. From Figures 1-4, the curve "-" represents the CDS price under the Vasieck model and the curve "--" represents the CDS price under the VCJ model. Some of the parameters values we used in the following numerical analysis refer to Wang and Liang [30], Leung and Kwok [31]. e basic parameters values are T � 5, r � 0.05 and a i � 0.5, b i � 0.1, σ i � 0.05 and L i � 1 for i � 1, 2, . . . , n, n + 1, n + 2. Figure 1 shows the impact of number of assets in the basket on CDS price under the Vasicek model and VCJ model. We set ρ ij � 0.3(i ≠ j), R � 0.4, λ i (t) � 0.1 for all i � 1, 2, . . . , n + 2. Generally the larger the number of assets in the basket is, the default event is more likely to occur. If the default probability increases, the CDS buyer should pay more for the basket CDS contract against lager default risk. Compared with the Vasicek model, the CDS price under the VCJ model is higher. is is because cojumps make the risk of default higher. Figure 2 shows the impact of the initial default intensities on CDS price under the Vasicek model and VCJ model. We set ρ ij � 0.3(i ≠ j), n � 10, R � 0.4. Higher initial default intensities of the companies lead to higher default risks, it is seen from Figure 2 that the basket CDS prices increase with the initial default intensities. e basket CDS price under the VCJ model is higher than that in the Vasicek model. We can also find that the curves in Figure 2 become flatter with larger initial default intensity. Intuitively, if the initial default intensity is too large, the CDS buyers will be less willing to pay more for the CDS contract. erefore, with the increase of the initial jump intensity, the basket CDS prices under the VCJ model and Vasicek model tend to be consistent. Figure 3 shows the impact of the recovery rate on CDS price under the Vasicek model and VCJ model. In Figure 3, we set ρ ij � 0.3(i ≠ j), n � 10, λ i (t) � 0.1 for all i � 1, 2, . . . , n + 2. Recovery rate R reflects the extent of loss. e greater the recovery rate is, the smaller the loss will be. As seen from Figure 3, the basket CDS price decreases when the recovery rate increases. Although the basket CDS price under the VCJ model is higher than that under the Vasicek model, the CDS prices under the two models become close with the increase of recovery rate. Figure 4 shows the impact of the correlation coefficient among the reference assets and counterparties on CDS price under the Vasicek model and VCJ model. We set Discrete Dynamics in Nature and Society R � 0.4, n � 10, λ i (t) � 0.1 for all i � 1, 2, . . . , n + 2. As seen from Figure 4, the basket CDS price with cojumps is higher and the CDS price under both models decreases when the correlation coefficient increases. Larger correlation coefficient means that if one of the companies defaults, the default probabilities of the counterparties will be quite larger. If both of the two counterparties default, the basket CDS buyer will lose a lot and get no compensations, thus the CDS contract will become worthless. Figure 5 and Figure 6 study the sensitivity of CDS price to initial default intensity under the Vasicek model and VCJ model respectively. For example, the partial derivative (zW 1 /zλ j ) means the change speed of W 1 with λ j when other parameters are fixed. In Figure 5, the curve "-"   represents the partial derivative of W 1 about λ j (j � 1, 2, . . . , n). e curve "--" represents the partial derivative of W 1 about λ n+1 . e curve "-." represents the partial derivative of W 1 about λ n+2 . W 1 is the premium paid by the credit protection buyer to the counterparty F C . When the initial default intensity of the reference asset and the counterparty F D become larger, the higher compensation the CDS buyer can obtain from the counterparty F C , so the partial derivative (zW 1 /zλ j )(j � 1, 2, . . . , n) and (zW 1 /zλ n+2 ) are always positive. On the contrary, when the initial default intensity of counterparty F C increases, the lower compensation that CDS buyer can obtain from the counterparty F C , so the partial derivative (zW 1 /zλ n+1 ) is always negative. It is seen from Figure 5 that the sensitivity of CDS price to the initial default intensity of counterparty F D is greater than that of the reference asset. is may be because the default of counterparty F D will increase the fear of credit protection buyer about the default of counterparty F C , so the CDS price is more sensitive to the initial default intensity of counterparty F D . e absolute value of the sensitivity of CDS price to initial default intensity is decreasing and tends to be stable with the increase of initial default intensity. e reason is that when the default risk is larger, the willingness of CDS buyer to pay excess premium for the default risk is decreasing. In Figure 6, the curve "-" represents the partial derivative of W 1 ′ about λ j (j � 1, 2, . . . , n). e curve "--" represents the partial derivative of W 1 ′ about λ n+1 . e curve "-." represents the partial derivative of W 1 ′ about λ n+2 . e sensitivity of CDS price to initial default intensity under the VCJ model is similar to that in Figure 5. Compared with the Vasicek model, it is worth noting that the sensitivity of CDS price to the initial default intensity of counterparty F D under the VCJ model decreases slightly. e default risk of counterparty F D may lead to a higher default risk of counterparty F C under the VCJ model. erefore, this concern will reduce the willingness of CDS buyer to pay premiums to counterparty F C .
From the above sensitivity analysis, we can see that the existence of cojump does have impacts on CDS prices. And then, we investigate the impacts of cojumps on the basket CDS prices. We show the basket CDS prices under different λ J in Figure 7 and 8. Without loss of generality, we set ε i to be constants for all i ∈ 1, 2, . . . , n + 2 { }. From Figure 7, we see the basket CDS price increases when λ J increases if ε i are positive constants. Positive constants ε i mean the default intensities jump upward, which can be caused by bad news. A larger λ J means the bad news is more likely to happen, thus the default probability and the basket CDS price will both increase. On the contrary, negative constants ε i mean the default intensities jump downward, which can be caused by good news. A larger λ J means the good news is more likely to happen, thus the default probability and the basket CDS price will decrease. Figure 9 shows the basket CDS prices under different ε i when λ J is fixed. From the previous analysis in Figures 7 and  8, we can find that the basket CDS prices increase when ε i increase. It can be seen that the curve in Figure 9 becomes flatter with larger ε i . Intuitively, larger ε i means larger default intensity, thus the CDS buyer will be less willing to pay more for the basket CDS contract.
It is generally assumed that there is only one defaultable counterparty in the traditional basket CDS pricing model. We investigate the basket CDS price differences with two defaultable counterparties and with only one defaultable counterparty as shown in Figure 10. e curve "-" represents the result under the Vasieck model and the curve "--" represents the result under the VCJ model. Take the Vasieck model as an example, when there are two defaultable counterparties, we assume the CDS buyer F A spends a total cost of W ′ (W ′ � W 1 + W 2 ) on the basket CDS contract. If there is only one defaultable counterparty under the Vasieck model, the CDS price is assumed to be W. While all the other model assumptions remain unchanged, we compute the price difference W ′ − W under different number of assets in the basket. In Figure 10, we set ρ ij � 0.3(i ≠ j), R � 0.4, λ i (t) � 0.1 for all i � 1, 2, . . . , n + 2. e basket CDS buyer will need to pay more for the basket CDS contract when there are two defautable counterparties. Especially when the number of reference assets in the basket is small (for example, less than 5), the price difference increases with the increase of the number of reference assets. is may be because the default risk is overestimated when considering two counterparties if the number of reference assets in the basket is small. However, when the number of assets in the reference assets is larger (for example, lager than 5), the probability that at least one asset in the basket defaults will increase. Higher default risk makes the CDS buyer do not want to pay more for the basket CDS contract. us no matter how many counterparties there are, there will be almost no price difference if the number of reference assets is large enough. For the VCJ model, we assume the CDS buyer F A spends a total cost of W ″ (W ″ � W 1 ′ + W 2 ′ ) on the basket CDS contract. e price difference of premium is W ″ − W. Compared with the Vasieck model, there are similar conclusions for the VCJ model. Because the default risk under the VCJ model is higher than that under the Vasieck model, the price difference under the VCJ model is larger. From the above analysis, we find that the basket CDS pricing with two defaultable counterparties has the advantage to explain some empirically phenomenons.

Conclusion
In this paper, we mainly study the pricing of the basket CDS with two defaultable counterparties based on the Vasicek processes and the Vasicek processes with cojumps. Using the PDE method, we obtain the approximate closed-form formula of the basket CDS price. In the numerical analysis, we analyze the impacts of main parameters such as the number of reference assets, initial jump intensity, recovery rate and correlations among assets, jump size and jump intensity of cojumps on the basket CDS prices. Comparing the basket CDS prices with two defaultable counterparties and that with only one defaultable counterparty, we find that the price difference will be obvious if the number of assets in the basket is not too large. However, there will be almost no price difference if the number of assets is large enough. e numerical results help us better understand the basket CDS pricing with multiple defaultable counterparties. Moreover, we investigate the basket CDS pricing under the reduced-form