Research Article  Open Access
Chaoqun Ma, Shengjie Yue, Yishuai Ren, "Pricing Vulnerable European Options under Lévy Process with Stochastic Volatility", Discrete Dynamics in Nature and Society, vol. 2018, Article ID 3402703, 16 pages, 2018. https://doi.org/10.1155/2018/3402703
Pricing Vulnerable European Options under Lévy Process with Stochastic Volatility
Abstract
This paper considers the pricing issue of vulnerable European option when the dynamics of the underlying asset value and counterparty’s asset value follow two correlated exponential Lévy processes with stochastic volatility, and the stochastic volatility is divided into the longterm and shortterm volatility. A meanreverting process is introduced to describe the common longterm volatility risk in underlying asset price and counterparty’s asset value. The shortterm fluctuation of stochastic volatility is governed by a meanreverting process. Based on the proposed model, the joint moment generating function of underlying logasset price and counterparty’s logasset value is explicitly derived. We derive a closedform solution for the vulnerable European option price by using the Fourier inversion formula for distribution functions. Finally, numerical simulations are provided to illustrate the effects of stochastic volatility, jump risk, and counterparty credit risk on the vulnerable option price.
1. Introduction
In general, options traded in the overthecounter (OTC) market are associated with counterparty credit risk. Since the global financial crisis of 20072008, regulators and investors have widely recognized the effect of counterparty credit risk on OTC options (see, e.g., Crépey [1], Sayah [2], and Breton and Marzouk [3]). Unlike the options traded in regulated exchanges, there is no organized exchange to ensure that option writers will fulfill the contractual obligations in the OTC markets. Once the option writer defaults, the option holder will get only a fraction of necessary payments required by the contract. The holders of OTC options are likely to suffer a large loss if they ignore the counterparty credit risk. Thus, the participants in the OTC markets need to take into account counterparty credit risk when pricing OTC options. Those options with counterparty credit risk are referred to as vulnerable options.
The valuation of vulnerable options has gained substantial attention in the finance literature. For modeling credit risk, the structural approach has been used to study vulnerable options by a lot of scholars, in which default events happen when the firm’s asset value is lower than a specific boundary. Based on this approach, the vulnerable European option pricing model is first presented by Johnson and Stulz [4]. Klein [5] extends the work of Johnson and Stulz [4] by assuming option writer to have other liabilities, in which the prices of underlying asset and the asset of counterparty are correlated. Subsequently, many extensions of the results in Klein [5] have been considered. For example, Hui et al. [6] and Hui et al. [7] discuss the pricing of vulnerable European option with stochastic default barrier. In contrast to constant interest rate in Klein [5], Klein and Inglis [8] incorporate stochastic interest rate into the vulnerable option pricing model. Liao and Huang [9] consider the possible counterparty default prior to the maturity of option; they derive explicit formulas for vulnerable options with stochastic interest rates. Taking jump into account, Tian et al. [10] investigate the vulnerable European option pricing where the dynamics of asset prices are governed by two correlated jumpdiffusion processes. More broadly, Niu and Wang [11] discuss the pricing of vulnerable European option under regime switching setting.
However, the above mentioned literature generally assume that the volatilities of the underlying asset and counterparty’s asset are constant, which is contrary to the actual situation, because it cannot explain many wellknown empirical evidences such as the volatility clustering and the volatility smile. It is obvious that models with stochastic volatility will be more realistic. Some scholars have suggested alternative models for the price of vulnerable option. Lee et al. [12] derive the explicit expression for the vulnerable option price where the price processes of underlying asset and counterparty’s asset satisfy the Heston stochastic volatility model. Wang et al. [13] extend the model of Lee et al. [12] and assume that the price processes of underlying asset and counterparty’s asset follow two correlated twofactor stochastic volatility diffusion processes; they obtain a pricing formula of vulnerable option in a special case. However, Lee et al. [12] and Wang et al. [13] do not consider the effects of jump risk of underlying asset and counterparty’s asset on vulnerable option price. A lot of empirical evidences show that the arrival of new information and sudden disastrous events may lead to discontinuous changes in asset prices and the jump components of models are understood as the effects of market shocks on the asset prices. Therefore, many researchers have proposed pricing models that incorporate both stochastic volatility and jumps components. For example, Bates [14] studies pricing problem of option where the dynamics of underlying asset follows a jumpdiffusion process with stochastic volatility, in which the jump components of asset price follow the compound Poisson process. However, Carr and Wu [15] and Huang and Wu [16] find the fact that the compound Poisson process cannot describe the many small jumps of financial assets. Moreover, many researches have considered more general jump structures. A more general process is Lévy process, which allows the jump components to have infinite activity and admits nearly an arbitrary distribution. Li et al. [17] apply Markov chain Monte Carlo (MCMC) methods to study infiniteactivity Lévy processes with stochastic volatility; their empirical results show that infiniteactivity jumps are essential for modeling the S&P 500 index returns. Subsequently, Zaevski et al. [18] derive the general formula for a European call option under stochastic volatility and tempered stable Lévy jumps. Liang and Li [19] study pricing of the European option under Normal Tempered Stable process with stochastic volatility. Li et al. [20] discuss the equity premium and option pricing under the general equilibrium under Lévy process and stochastic volatility.
Along with the above research lines, based on the structural approach for default, we study the valuation of vulnerable European option when the price processes of the underlying asset and counterparty’s asset follow two correlated exponential Lévy processes with stochastic volatility. The volatilities of underlying asset and counterparty’s asset are both decomposed into longterm and shortterm volatility. The shortterm variance of asset price follows a meanreverting process. A meanreverting process is introduced to allow for the common longterm volatility risk in the underlying asset price and counterparty’s asset value, which differs from assumption that longterm volatility is constant in Wang et al. [13]. This assumption is consistent with the models in Christoffersen et al. [21] and Wong and Zhao [22]. The market is fluctuating in reality, and the market factor affects the longterm volatilities of underlying asset price and counterparty’s asset value. Thus, the longterm volatility of asset price is also changing over time. Based on the proposed model, the joint moment generating function of underlying logasset price and counter’s logasset value is derived. Moreover, we obtain an explicit formula for vulnerable European call option price by applying the Fourier inversion transform. Finally, sensitivity analysis and numerical simulation are provided for our results. By the numerical analysis, we discover that (i) the longrun mean of common longterm variance of underlying asset and counterparty’s asset exerts positive effect on the vulnerable option price; (ii) vulnerable option price is an increasing function with respect to the longrun mean of shortterm variance of underlying asset price; (iii) a higher longrun mean of shortterm variance of counterparty’s asset value induces a lower vulnerable option price; (iv) when the jump components of underlying asset price and counterparty’ asset value follow the CGMY model, the aggregate activity rate of jumps of underlying asset price has positive effect on the value of vulnerable option price. The higher aggregate activity rate of jumps of counterparty’s asset value, the lower vulnerable option price.
Our contributions to the literature are summarized below: First, we extend the model of Wang et al. [13] by considering a more general financial market. In particular, compared to the constant longterm volatilities of underlying asset price and counterparty’s asset value in Wang et al. [13], we assume that the common longterm volatility risk in the underlying asset price and counterparty’s asset value is governed by a meanreverting process. Besides, we simultaneously incorporate infiniteactivity Lévy jumps and stochastic volatility into the dynamics of underlying asset and counterparty’s asset, which differs from the jumpfree model in Wang [13]. Second, we derive an explicit expression for vulnerable option price by the Fourier inversion formula for distribution functions, which is different from the methods in Lee et al. [12] and Wang et al. [13].
The rest of this paper is organized as follows. Section 2 introduces the basic model setup and formulates the pricing problem of vulnerable European call option problem. In Section 3, the joint characteristic function of the underlying logasset price and counterparty’s logasset value is derived. By applying Fourier inversion transform, we derive an analytic solution for vulnerable European option price. Section 4 demonstrates our results with numerical examples, and Section 5 concludes this paper.
2. Basic Setting of the Model
In this section, we describe theoretical framework for valuing of vulnerable European call option. The underlying asset price and counterparty’s asset value follow exponential Lévy processes with stochastic volatility. The volatility risk is divided into shortterm risk and longterm risk. We incorporate a meanreverting process to describe the common longterm volatility risk in underlying asset price and counterparty’s asset value. The shortterm fluctuation of stochastic volatility also follows a meanreverting process.
Time is indexed by , where T is a positive finite constant representing the time horizon. To model uncertainties, suppose that we have a given complete probability space with filtration , and the filtration satisfies the usual condition (increase, rightcontinuous and augmented). It represents the flow of information available to the investor. is a riskneutral probability. Suppose , , , , and are all standard Brownian motions which are adapted. Moreover, cov, cov, cov, cov, cov, and any other Brownian motions are pairwise independent. Under a risk neutral measure , the underlying asset price is modeled as the exponential of the lévy process with stochastic volatility:where stands for the value of before a possible jump occurs, ; is the interest rate, and the parameters , , and are the meanreverting rate, the longterm mean; and volatility of the variance process , respectively. We require to ensure that the process remains strictly positive. represents the longterm variance of underlying asset price, which reflect the effect of corporate performance and market states on the underlying asset price. Here we assume that follow a meanreverting process, which represents the common systematic risk caused by market factors. is shortterm variance of underlying asset value, which depends on investors’ trading and reflects company’s idiosyncratic risk. In addition, is a compensated jump measure, where is the jump measure and Lévy kernel (or density) satisfies ; the underlying asset price in (1) contains two orthogonal martingales: two purely continuous martingale and a purely discontinuous (jump) martingale.
Next, let denote the market value of counterparty’s asset, and the dynamic of is given by the SDE under risk neutral measure ,where stands for the value of before a possible jump occurs, ; represents the longterm variance of counterparty’s asset value, and describes the systematic risk caused by market shocks. is shortterm variance of counterparty’s asset price, which depends on the investor’s trading activities. The parameters , , and are the meanreverting rate, the longterm mean, and volatility of the variance process , respectively. We require to ensure that the process remains strictly positive. In addition, is a compensated jump measure, where is the jump measure and Lévy kernel (or density) satisfies .
In the proposed framework, we consider stochastic volatility and jumps in the dynamics of both the underlying asset and the counterparty’s asset, and the longterm stochastic volatility and the shortterm stochastic volatility are not constants. The underlying asset price and the counterparty’s asset value are correlated, and the correlation coefficient is . Furthermore, correlation between stochastic volatility and asset prices is also considered in the proposed model. The shortterm stochastic volatilities of the underlying asset and the counterparty reflect company’s idiosyncratic risk while the longterm stochastic volatilities of the underlying asset and the counterparty describe systematic risk caused by common market shocks. In addition, we also focus on the discontinuous changes in the underlying asset and counterparty’s asset; the jump components in (1) and (4) reflect the abnormal vibrations in the asset price due to the arrival of important new information and sudden disastrous events.
3. Vulnerable European Call Option Pricing
In this section, we derive the joint moment generating function of the underlying logasset price and counterparty’ logasset value. In this framework, by using the Fourier inversion transform, we obtain the explicit expression for vulnerable European option price.
Let denote the option maturity; for vulnerable European call options, the boundary condition is determined by a payoff function depending upon the financial distress circumstance at the time . As in Klein (1996) and its subsequent results, then the payoff of a vulnerable European call option is given bywhere is the strike price of the option. is the deadweight cost related with the bankruptcy or reorganization process of the firm expressed as a percentage of the counterparty’ asset value. is the total amount of claims of the counterparty, and is an exogenous default barrier. may be less than due to the possibility of a counterparty continuing its operation. A credit loss happens if the counterparty’s asset value is less than the amount . If is higher than or equal to the default barrier , the claim is paid out in full. Otherwise, default event occurs and only the proportion of the nominal claim at is paid out by the option writer.
By the standard riskneutral arguments, under a riskneutral measure , the price of a vulnerable European call option at time is given aswhere denotes risk neutral expectations. denotes the stochastic discount factor. Expression (7) shows that the value of a vulnerable European call option contains two parts which are conditional on whether is higher or less than the default barrier under riskneutral measure . If is higher than , the option writer can afford to pay all the claim. If is lower than , only a proportion of the nominal claim is paid.
In the following, we derive the explicit expression for the value of vulnerable option price by using Fourier inversion formula for the probability distribution functions. In the process of calculation, we need to use the joint moment generating function of and . Thus, it is necessary for us to give this definition of the joint moment generating function. Given the dynamics of the underlying asset price and the counterparty’s asset value under the probability measure , it is possible to obtain the joint moment generating function of and . Let , , and the joint moment generating function of and under probability measure is defined aswhere , . Then, the following lemma holds.
Lemma 1. Suppose that underlying asset price and counterpart’s asset value follow the dynamics in (1) and (4), respectively. Let , , , , , and the joint moment generating function of and in (8) is given bywherewhere
Proof. See Appendix A.
Now, we turn to derive the pricing formula for vulnerable option pricing. We need to rewrite pricing formula (7) of vulnerable European call option as follows:where are given by
In the following, from the standard probability theory, we can obtain the probability distribution function by using the joint moment generating function. In Theorem 2, we can derive the explicit expression of vulnerable European option by applying Fourier inversion transform.
Theorem 2. The price of vulnerable option pricing at time is formally given bywhere is the moment generating function in (8), and , , , and are given in (B.11), (B.12), (B.17), and (B.22), respectively.
Proof. See Appendix B.
4. Numerical Analysis
In this section, we present numerical results of the vulnerable European call option price under Lévy process with stochastic volatility. Numerical examples are presented to investigate the effects of the parameters on the vulnerable European call option. The pricing models of vulnerable option in Klein [5] and Wang et al. [13] and pricing models of nonvulnerable option in Sun [23] and Li et. al [20] are chosen as reference models. Our choice of the four reference models is made to study the effects of stochastic volatility, jump risks, and credit risk on vulnerable option prices. In the following numerical analysis, we change one of the parameter values to illustrate its effect on the price with other parameters taking on the values in Tables 1–3. Our model allows a general distribution for jump components of underlying asset price and counterparty’s asset value and thus can be easily deduced to the simple cases such that the jump components follow compound Poisson process in Merton [24] and Kou [25], the CGMY model of Carr et al. [26], etc. By assuming different types of Lévy kernel, our model has different forms. In the following, we will first study a simple case with jump components following compound Poisson processes and then a more general case with the jump components following a CGMY model.



In the numerical analysis, we use the function quadgk available in MATLAB to deal with the single integrals for the Fourier inversion transform. The double integrals for the Fourier inversion transform are computed by using the function integral2 available in MATLAB.
4.1. The Merton Model
When jump components of underlying asset and counterparty’s asset follow compound Poisson processes, the jump magnitude of asset price is assumed to be drawn from a lognormal distribution. Furthermore, the Lévy kernels are given bywhere is the jump intensity of the asset process, is the mean of jump size, and is the variance of jump size. Following the Merton [24] and Sun [23], the two integrals in (13) are given by
In Table 2, we will give the specimen values for the model parameters of assets’ jump components. Figures 1–5 illustrate the effects of basic parameters on vulnerable option prices where the parameters include the default barrier, longrun mean of stochastic volatility, and jump intensity. Denote prices of nonvulnerable European option in Sun [23] by , and let be values of vulnerable option in proposed model. We denote the values of vulnerable option in Wang et al. [13] and Klein [5] prices by and , respectively.
Figure 1 illustrates that the vulnerable option price is a decreasing function of default barrier . When the default barrier is relatively small, the value of vulnerable European option decreases slowly with ; while the default barrier is close to the initial value of counterparty’s asset, the vulnerable option price decreases more dramatically with . Large default barrier will increase the default probability; the counterparty credit risk has a more significant impact on the vulnerable option price if the default barrier increases. In addition, the value in proposed model is higher than those of models proposed by Klein [5] and Wang et al. [13]. The phenomenon is consistent with the results presented in existing literature, because the jump parts of underlying asset increase the chance that the underlying asset price will become very low or very high at maturity; the option holder will gain benefit from the increase of underlying asset price and suffers limited downside risk when underlying asset price decreases. Therefore jump risk enhances option value. Moreover, the value of nonvulnerable option price is higher than those of vulnerable European options. Once the credit events happen, the option holder will only get a fraction of promised payoff; thus, the vulnerable option is less valuable.
Figure 2 displays that the value of option price with stochastic volatility increases with the longrun mean of the process . The reason for this is that a higher value yields a higher value of . Large will increase probability that the underlying asset price and counterparty’s asset value will become very low or very high at maturity; when the underlying asset price and counterparty’s asset value increase, the option holder will gain benefit; when underlying asset price and counterparty’s asset value decrease, the option holder faces a limited loss. Therefore the option becomes valuable if increases.
Figure 3 plots the effect of the longterm mean of stochastic volatility on vulnerable option price. The solid and dotted lines correspond to the vulnerable option prices obtained by the proposed model against longrun means of shortterm variances and , respectively. For parameter , we can observe that exerts positive effect on the vulnerable option price. The reason is the same as parameter . Large longrun mean of variance will increase the probability that underlying asset price turns into very high or very low at maturity. Consequently, the vulnerable option price is valuable under large . For parameter , the large longrun mean of variance will lead to counterparty’s asset value being very high or very low at maturity. Consequently, the default probability will increase. Thus, the vulnerable option price decreases with .
As shown in Figure 4, the value of option with jump risk increases with the jump intensity . The reason is the same as that of . The jump intensity will result in underlying asset price doing very poorly and very well at maturity, but the option holder has a limited loss. Thus, the option with jump risk is valuable under large jump intensity . Figure 5 shows that the jump intensity has a negative effect on the value of vulnerable option with jump risk. Counterparty credit risk is related to the jump intensity of the counterparty’s asset. A stronger jump intensity of counterparty’ asset value corresponds to a larger default probability. Accordingly, the vulnerable option price decreases with jump intensity . In addition, the vulnerable option price is a decreasing function of , the deadweight cost coefficient affects the recovery rate when default event happens, and a higher value of gives rise to a smaller recovery rate. Thus, the option holder suffer a greater loss under large .
4.2. The Double Exponential Jump Model
When jump components of underlying asset and counterparty’s asset follow compound Poisson processes, we assume that the jump magnitude of asset price has asymmetric double exponential distribution. The Lévy kernels are given bywhere is the jump intensity of the asset process and and are probability of upmove jump and probability of downmove jump, respectively. and are mean of positive jumps and mean of negative jumps, respectively. Following Kou [25], the two integrals in (13) are given byIn Table 3, we will give the specimen values for the model parameters of assets’ jump components. Figures 6 and 7 illustrate the effects of jump intensity and jump intensity on vulnerable option price.
Figures 6 and 7 illustrate the sensitivities of the vulnerable option price with respect to the jump intensity and the jump intensity . From Figure 6, we observe that high jump intensity leads to a high vulnerable option. Large jump intensity will increase the chance that the underlying asset price will become very low or very high; the option holder will benefit from the price increase. Figure 7 shows that the vulnerable option price is a decreasing function of . Large jump intensity will increase default probability. Therefore the vulnerable option price becomes less valuable if increases.
4.3. The CGMY Model
A more general case can be considered if jump components of underlying asset and counterparty’s asset follow CGMY model developed by Carr et al. [26], and the Lévy kernels are given bywhere , , , and . determines the aggregate activity rate of the CGMY model while the parameters and determine the height of the upward and downward jumps, respectively. controls both tails of the spot return distribution. Based on this, following Carr et al. [26], the two integrals in (13) are given by
In Table 4, we will give the specimen values for the model parameters of jump parts of underlying asset price and counterparty’s asset value. Figures 8–11 illustrate the effects of basic parameters on vulnerable option prices. We denote prices of nonvulnerable European option in Li et al. [20] by .

From Figure 8, we can see that the value of option with jumps increases with parameter . The parameter determines the aggregate activity rate of jumps of the underlying asset price; the large jump frequency will lead to the underlying asset price doing either very well or very poorly at maturity. While the option holder has limited downside risk, the vulnerable option becomes valuable if increases. As shown in Figure 9, the price of option with jumps is a monotonic increasing function of . The parameter controls both tails of the underlying return distribution. As increases, both tails of underlying return distribution will become fatter, and the underlying asset price becomes very high at maturity; it is not surprising that the option with jumps would become more valuable if increases.
Figure 10 denotes vulnerable option price is a decreasing function of parameter . The parameter controls the aggregate activity rate of jumps of the counterparty’s asset value. Large will increase the probability that counterparty’s asset value becomes very low at maturity, which implies that the default probability will increase. Thus, the vulnerable option price will decrease with . Figure 11 shows that effect of parameter on the vulnerable option price; the parameter controls both tails of the spot return distribution of counterparty’s asset. As increases, default probability will increase. Accordingly, the vulnerable option price is likely to decrease with .
5. Conclusion
This paper investigates the pricing problem of vulnerable European call options where the price processes of underlying asset and counterparty’s asset are governed by two correlated exponential Lévy processes with stochastic volatility. We divide volatility risk of underlying asset and counterparty’s asset into longterm and shortterm risk. A meanreverting process is introduced to allow for the common longterm volatility risk in the underlying asset price and counterparty’s asset value. The shortterm fluctuation of stochastic volatility also follows a meanreverting process. Based on the proposed model, we explicitly obtain the joint moment generating function of underlying logasset price and counter’s logasset value. By using the Fourier inversion transform, an analytic solution of vulnerable European call option is derived. Finally, sensitivity analysis and numerical simulation are provided to illustrate the effects of model parameters on the option price. The main findings are as follows: (i) the longrun means of common longterm variance of underlying asset and counterparty’s asset have positive effect on the vulnerable option price. A large longrun mean of common longterm variance will increase the possibility that underlying asset price and counterparty’s asset value will become very well or very poor at maturity. The option holder will benefit form increase of underlying asset price and faces a limited loss. Therefore the vulnerable option becomes valuable. (ii) The value of vulnerable option price increases with longrun mean of shortterm variance of underlying asset price and decreases with longrun mean of shortterm variance of counterparty’ asset value; (iii) the value in proposed model is higher than those of Wang et al. [13] and Klein [5]. The jump components of underlying asset price will increase the probability that the underlying asset price will turn to very low or very high at maturity, and the vulnerable option with jump risk becomes valuable; (iv) when the jump components of underlying asset price and counterparty’ asset value follow the CGMY model, the value of option price increases with the aggregate activity rate of jumps of underlying asset price and decreases with the aggregate activity rate of jumps of counterparty’s asset value.
However, two further extensions of vulnerable European option pricing would be interesting. One possible extension is to incorporate stochastic interest into pricing problem of vulnerable option. Another extension would be to extend the study on vulnerable option pricing under Markov regimeswitching setting.
Appendix
A. Proof of Lemma 1
To obtain the joint moment generating function for the logasset price and the logassets value , we first take the following transformation:Let , , , , and , and the generalized FeynmanKac theorem implies that solves the following partial integrodifferential equation (PIDE):The boundary condition for (A.3) is given byThis PDE has an exponentialaffine solution of the formwith the boundary conditionsSubstituting (A.5) into PDE in (A.3) yieldsBy matching coefficients, we can derive the following three ordinary differential equations (ODEs):According to the boundary condition (A.6), the above ODEs can be easily solved as below:where
B. Proof of Theorem 2
Proof. Recall the joint moment generating function of and in Lemma 1,To deal with the term in (16), we define a new probability measure on by the following RandonNikodym derivative:Based on the definition of , we have the characteristic function of and under ,where we have used the definition of in Appendix A.
From standard probability theory (see, e.g., Kendall and Stuart [27] and Shephard [28]), one obtains that the distribution function corresponding to the joint characteristic function is given bywhere represents real part of a complex number, and and are the marginal distributions for and ,Here, we have thatRecalling the definition of the probability measure , we have thatRecalling the definition of the probability measure , we have thatwhich in turn implies