Abstract

The pricing problem of a kind of European vulnerable option was studied. The mixed fractional Brownian motion and the jump process were used to characterize the evolution of stock prices. The closed-form solution to European option pricing was obtained by applying martingale measure transformation method. At the end of this paper, some numerical experiments were adopted to compare the new pricing formula introduced in this paper with the classical Black-Scholes pricing formula. The result showed that the new pricing formula conformed to the actual financial market. In fact, the option value is positively correlated with the underlying asset price and the company’s asset price and the jump process has significant influence on the value of option.

1. Introduction

Vulnerable European option is a kind of option with credit risk, which refers to the possibility of one party suffering losses as a result of the other party’s default on the contract. The pricing of credit risk has been concerned by people very early. Merton (1974) [1] introduced option pricing theory into zero coupon bond credit risk pricing at first time. He supposed that the capital structure of company consists of asset and debt. A default occurred when the company becomes insolvent at the expiration of the option. In this model, a baseline of structured-form model ideas was built. Later, the irrational assumption was improved in Merton’s model and the first passage time model was proposed by Black and Cox (1976)[2]. Johnson and Stulz (1987) [3] firstly led default risk to option pricing and put forward a new conception called vulnerable option. Klein (1996) [4] obtained the pricing formula of vulnerable option by the method of martingale. The default risk was decomposed into timing and recovery risks by Madan (1998) [5] who generated estimation strategies evaluating the recovery risk and then constructed implicit prices of contingent securities reflecting the timing risk. Lobo (1999) [6] and Zhou (2001) [7] considered the pricing problem of European option with the credit risk that the underlying asset price carried jump, but there were no analytical solutions in the end. Klein’s option pricing model was expanded by Ammann (2001)[8] by applying structured method. They derived the closed-form solution to vulnerable option price in the environment where interest rate and default intensity complying the stochastic differential equation. In addition, Chang and Huang (2006) [9] also dealt with the problem. In the case of the underlying interest rates following a mean reverting square-root process, the analytical solutions to defaultable bonds were derived by Cathcart (2003) [10]. The reduced-form model was dated back to Jarrow and Turnbull (1995) [11]. Assuming that the default event obeyed the independent homogeneous Poisson process, the discrete methods were used to establish a simplified credit risk model by Jarrow and Turnbull. Cathcart (2006) [12] derived the price of the defaultable bond and proposed a valuation model that combined features of both the structured-form and reduced-form approaches for modelling default risk. Later, the problem that a new model of credit risk was proposed combining structural information with the reduced-form approach also was studied by Ballestra (2014) [13] and a closed-form approximate solution was derived by perturbation approach and the Laplace transform. A quasianalytical approximation of the survival probability was derived based on Madan’s credit risk model, and a quasianalytical approximate expression for CDS par spreads was obtained by Ballestra (2017) [14]. All the discussions above are based on the underlying asset price obeying the geometric Brownian motion, which implies that the changes of underlying asset price are mutually independent random variables and the return on assets obeys the normal distribution. But in recent years many financial empirical studies showed that the distribution of return on assets had the feature of high peaks and heavy tails. Changes in stock price showed different degrees of long-term correlation and autocorrelation at different time. These characteristics are different from standard Brownian motion but similar to fractional Brownian motion. All of these features make fractional Brownian motion more widely used in assets pricing. Sottinen and Valkeila (2001) [15] proposed that there is no equivalent martingale measure in the fractional Black-Scholes market. Yaozhou Hu and øksendal (2003) [16] studied the structural model of credit risk under fractional Brownian motion and proved that the Black-Scholes market driven by fractional Brownian motion is a nonarbitrage market. More studies about fractional Brownian motion are introduced by øksendal (2004) [17]. However, there were some shortcomings when Fractional Brownian motion is utilized to characterize the volatility of asset price by Björk and Hult (2005) [18]. In addition, there is a main problem that fractional Brownian motion is not a semi-martingale when it is applied in financial market. In order to avoid these problems and consider the long memory characteristics of the financial asset price process, it is reasonable to use mixed fractional Brownian motion to characterize the volatility of financial assets. Mixed fractional Brownian motion is a Gaussian process and it is a linear combination of Brownian motion and fractional Brownian motion. He is Cheridito (2001) [19] who initially put mixed fractional Brownian motion of economics into use. Androshchuk (2006) [20] proved that the mixed Brownian-fractional-Brownian model had no arbitrage opportunity in the class of self-financing Markov-type strategies. A European call option pricing formula was derived in mixed Brownian-fractional-Brownian environment with a mean self-financing delta hedging argument in a discrete-time setting and the minimal pricing of European option under transaction costs was obtained by Wang, Zhu and Tang (2010) [21]. Ballestra (2016) [22] studied the pricing barriier options by solving an initial-boundary value partial differential problem, in which an ad-hoc numerical procedure was employed based on product integration. Besides, Necula (2002)[23] introduced the Black-Scholes pricing formula by using quasi-martingale method in the framework of risk neutral measure, which is a generalization of the classical Black-Scholes formula. Wenli Huang (2015) [24] considered the general European option pricing model with stochastic interest rate driven by fractional Brownian motion and partial differential equation. Xiaonan Su (2012) [25] supposed that the intensity of default is driven by a jump diffusion process and got vulnerable option pricing formula in a reduced-form by using martingale method. Based on Xiaonan Su, Jinzhi Li (2013) [26] used the change of measure to obtain the closed-form solution of vulnerable option pricing. Assuming that the stock price obeyed the jump-diffusion model, Chao Wang (2015) [27] got analytic solution to vulnerable European option based on Xiaonan Su. Adopting the Ito’s formula of mixed fractional Brownian motion, Zhiguang Li (2016) [28] obtained the European option pricing formula with short-term interest rate obeying Vasicek model.

In this paper, we will mainly study the pricing problem of a class European vulnerable option under a mixed fractional Brownian motion environment. The paper is organized as follows: In Section 2, we will give pricing formula of structured-form and reduced-form model respectively. Then we will obtain the pricing formula of vulnerable option by combing the two models into Hybrid model when the stock price obeys mixed fractional Brownian motion under the jump-diffusion model. In Section 3, the analytic solution is obtained by the pricing formula. In Section 4, in order to investigate the impact of various parameters on option pricing, we give the figures of option value under different parameters with respect to the underlying asset price changes. In addition, for the reader’s convenience, the last part of this paper is an appendix to the concrete proof process of vulnerable option pricing.

2. Hybrid Model

Merton (1974) [1] first introduced the structured-form model, which is a default model based on specific information of the company. He thought that default was related to the company’s own business condition. The structured-form model considered the company’s asset as underlying asset, the value of creditor’s rights as the strike price, and used the option pricing method to consider default. In Merton’s model, he assumed that there would be no default when the option expired, but this obviously did not match the actual financial market. Later, Black and Cox (1976) [2] improved the unreasonable assumption in Merton’s model and proposed the first passage time model. The model set a boundary for company’s asset and when the value of the asset fell below this boundary, a default happened. Under risk-neutral probability measure, it is assumed that the firm’s value process satisfies the stochastic differential equationwhere is risk-free interest rate, is standard Brownian motion on probability space , and is fractional Brownian motion whose parameter is on this space. is the stock price volatility driven by standard Brownian motion and is volatility driven by fractional Brownian motion.

It is supposed that there exists a default threshold . When , the company will default, and on the contrary, there will be no default. Hence the default time is given by

Klein (1996)[4] assumed that if a default happened, the option seller will only pay part of the value of the equity on maturity and we denoted this ratio by where the percentage of the cost of the bankruptcy to the option seller’s asset is .

The reduced-form model can trace back to Jarrow and Turnbull (1995) [11]. Unlike the structured-form model, it did not consider the relationship between the default and asset value. It regarded the process of default as a jump process. The model treated the default event as a kind of jump and often portrayed the default as a Poisson process with intensity . The probability of nondefault at the moment is and is defined as the default density. Later Lando (1998) [29] introduced Cox process and supposed that the default process was characterized by a Cox process. The default time was the first time that Cox process jumped. This model assumed that is a Cox process and its random intensity process satisfiedAccording to the reduced-form model, we know that the value of the European vulnerable call option at time 0 is as follows:where the moment of default is , expiration dare is , strike price is , and the asset recovery rate is .

Applying the method of structured-form model, let and the value of European vulnerable call option at time 0 is given by

We consider to combine the reduced-form model with structured-form model and to define the moment of occurrence of default as the moment of jump or the moment when the default boundary is reached for the first time. We assume that the asset price obeys the mixed fraction Brownian motion and the stock price is subject to the mixed fraction Brownian motion under the jump diffusion model. Then the vulnerable option price model is given bywhere is constant and denotes the default boundary and follows the assumption that .

By law of total expectation and Fubini’s theorem, we havesince the full path of at time is known, applying (4), we have [30]soTherefore, (7) can be represented as

We assume that stock price and company’s asset price under risk-neutral martingale measure obey the next two processeswhere is Poisson process and satisfies and denotes the jump range of J’s jump . denotes Poisson process with a strength of and is independent of , . , , , , is a list of independent and identically distributed random variables and is independent of , , . Assumed that obeys lognormal distribution, that is, . By lemma, we have

Supposing that default intensity obeys model under risk-neutral measure where , , , and are constants and denotes the intensity of the Cox process characterizing the default process.

The covariance coefficient of , , and isand similarly, the covariance coefficient of , , and is

3. Vulnerable Option Pricing Model

Theorem 1. Let stock price and asset price satisfy (12) and (13) and the default intensity satisfies (16); then the value of the European vulnerable call option at time iswhere is strike price, is the riskless interest rate, is maturity, is default boundary, is default recovery rate, and , and are constant, obey two-dimensional normal distribution, , , where , , , , and are constant and , , and .

Proof. For convenience, (11) can be written aswhereBy calculating, we haveWe show the derivation process of . For and , the derivation process is included in the Appendix A and B which are similar to . We definewhereStep one, we calculate . At first, we introduce a new measureSubstituting (15) into , we haveFor satisfying theorem, we havethus under measure, (14) can be written asforand when , we havethusthen (30) is written aswhereStep two, we calculate . At first we introduce a new measureBy multidimensional theorem, we getunder measure, (14) can be written asBy introducing derivative, we haveand similarly, when , there existsand hence we havewhereTherefore, we haveThe proof is finished.

4. Numerical Test

The images of option value under different parameters with respect to the underlying asset price changes are given respectively in order to investigate the influence of various parameters on the value of the vulnerable option. Assuming that each parameter has the following values: the constant interest rate is , the strike price is , the maturity is , the default threshold is , the expected value of jump range is , the percentage of the cost of the bankruptcy to the option seller’s asset is , the stock price volatility under standard Brownian motion are , , , the stock price volatility under fractional Brownian motion are , , , the expectation and standard deviation of are , , the correlation coefficients between and are .

The value of vulnerable option is obtained from (19). By using software, we show some results of numerical tests.

Figure 1 shows that the value of the option changes with both underlying asset and the price of the company’s asset. As we can see from the figure, the option value increases with the increase of the underlying asset price and the company asset price. The underlying asset price is positively correlated with the company asset price and the option value.

From Figure 2, when there is a jump process that the jump range expectation is positive, the option value is higher than the normal one and this is in line with the actual situation. This situation shows that the jump process has a significant influence on option pricing. In actual market, introducing jump process can make option pricing more accurate. Since in this article the situation of default is considered, the value of the option is slightly lower than the value of the classical Black-Scholes option and this is also consistent with the actual situation. The introduction of default in this article can make option pricing more perfect.

Figures 3 and 4 depict the relationship between option value and underlying asset under different parameters respectively. Figure 3 depicts the value of options under different default intensity and the figure shows that as the default intensity increases, the option value is also growing. Figure 4 depicts the change of option value under different boundaries. From the figure, as the default boundary increases, the option value gradually decreases due to the increased default risk. This situation agrees with the actual market.

Figure 5 depicts the effect of the correlation coefficients and on the price of the option. As we can see from Figure 5, as the correlation coefficient increases, the option value continues to decrease. However, the correlation coefficient in (a) has a greater impact on option price than that in (b).

5. Conclusion

The traditional definition of default probability is divided into two types: reduced-form model and structured-form model. The reduced-form model supposes that default occurs only when jump occurs and structured-form model supposes that default occurs when the company’s asset value is less than the default boundary. The probability of default is determined by the relationship between the company’s asset value and default boundary. We employ the combination of reduced-form model and structured-form model to pricing European vulnerable option. In this paper, it is assumed that the default can be caused by jump process or by the insolvent of company. The expression of vulnerable option value is given under the condition of risk-free interest rate and the analytical solution is obtained. Finally, we conduct numerical experiments on the analytical solutions and analyze the results.

Appendix

A. Derivation Process of II

Derived of . Definewhere

Step one, we calculus . The equivalent martingale measure of defined by derivative definition of as follows

Solving the solution of Equation (16), we derivewe integrate two sides from exclusivelyapplying Fubini’s theorem, we obtainletwhere is given asletusing multidimensional fractional lemmaso

Using fractional Girsanov’s theoremaccording to Equation (14), therefore, we have the expression for and under measure

Using derivativewhen , we bring and into sowhere

Then we calculate , let measure inletunder the measure , we bring Equation (A.14) into Equation (A.21)hence

Using fractional Girsanov’s theorem, we deriveunder measure , we have

Using derivative, we havewhen , we bring and into sowhereso

B. Derivation Process of III

Deduced . We definewhere

Firstly we calculate . Let new measure inletunder measure , we havethen we haveso

Using fractional Girsanov’s theoremsince under measure using derivative, we obtainwhen , we bring and intosowhere

Then we will calculate . Let measure inletunder measure we haveso

Using fractional Girsanov’s theorem, we havewe have the expression for and under measure using derivative, we obtainwhen , we bring and intohencewhereso

The proof of the theorem is completed.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work is supported by the National Natural Science Foundation of China (no. 61304088) and the Fundamental Research Funds for the Central Universities (no. 2013QNA37).