`Abstract and Applied AnalysisVolume 2013 (2013), Article ID 780542, 13 pageshttp://dx.doi.org/10.1155/2013/780542`
Research Article

## Equilibrium Asset and Option Pricing under Jump-Diffusion Model with Stochastic Volatility

1School of Economic Mathematics, Southwestern University of Finance and Economics, Chengdu 611130, China
2Department of Mathematics & Statistics, Curtin University, Perth, WA 6102, Australia

Received 24 July 2013; Revised 30 September 2013; Accepted 10 October 2013

Copyright © 2013 Xinfeng Ruan et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

We study the equity premium and option pricing under jump-diffusion model with stochastic volatility based on the model in Zhang et al. 2012. We obtain the pricing kernel which acts like the physical and risk-neutral densities and the moments in the economy. Moreover, the exact expression of option valuation is derived by the Fourier transformation method. We also discuss the relationship of central moments between the physical measure and the risk-neutral measure. Our numerical results show that our model is more realistic than the previous model.

#### 1. Introduction

Option pricing problem is one of the predominant concerns in the financial market. Since the advent of the Black-Scholes option pricing formula in [1], there has been an increasing amount of literature describing the theory and its practice. Due to drawbacks of the Black-Scholes model which cannot explain numerous empirical facts such as large and sudden movements in prices, heavy tails, volatility clustering, the incompleteness of markets, and the concentration of losses in a few large downward moves, many option valuation models have been proposed and tested to fit those empirical facts. Jump-diffusion models with stochastic volatility could overcome these drawbacks of the Black-Scholes model in [221]. Based on those advantages, in this paper, we focus on studying the jump-diffusion model with stochastic volatility.

Different from the Black-Scholes framework, we use jump diffusion to describe the price dynamics of underlying asset. The market of our model is incomplete; that is, it is not possible to replicate the payoff of every contingent claim by a portfolio, and there are several equivalent martingale measures. How to choose a consistent pricing measure from the set of equivalent martingale measures becomes an important problem. This means that we need to find some criteria to determine one from the set of equivalent martingale measures in some economically or mathematically motivated fashion. A unique martingale measure was found by various researchers via using optimal criteria, for instance, minimal martingale criterion, minimal entropy martingale criterion, and utility maximization criterion [2231].

General equilibrium framework method is also a popular method to deal with the option pricing in an incomplete market. General equilibrium framework is initially introduced by Lucas Jr. (1978) [32], Cox et al. (1985) [33] and developed by Vasanttilak and Lee (1990) [34], Pan (2002) [35], Liu and Pan (2003) [36], Liu et al. (2005) [37], Bates (2008) [38], Santa-Clara and Yan (2004) [12], and Zhang et al. (2012) [6]. They assumed that there is a representative investor who wants to maximize an objective function in a rational expectations economy where there are one risk-free asset and one risky asset. When the market is clear, the representative investor takes all money into the risky asset. In this paper, we build a general equilibrium model which is the same as that due to Santa-Clara and Yan (2004) [12]. Under this model, we obtain an exact expression of the equity premium and the pricing kernel in a general equilibrium economy. This can be regarded as a great contribution to the literature.

The pricing kernel which acts like the physical and risk-neutral densities and moments in the economy is also a vitally important problem in mathematical finance. In some constant volatility models with jump diffusions, Pan (2002), Liu and Pan (2003), and Liu et al. (2005) [3537] derived the pricing kernel with some restrictions of jump sizes in a general equilibrium setting. Recently, Zhang et al. (2012) [6] presented an analytical form for the pricing kernel without any distributional assumption on the jumps. In this paper, we extend the results of Zhang et al. (2012) to the pricing kernel with stochastic volatility.

Duffie et al. (2000) [10] and Chacko and Das (2002) [39] presented a transform analysis to price the valuation of options for affine jump diffusions with stochastic volatility. Lorig and Lozano-Carbasse (2013) [21] studied option pricing in exponential Lévy-type models with stochastic volatility and stochastic jump intensity. Lewis (2008) [8] used Fourier transformation methods to obtain the transform-based solution of option price. In this paper, we employ the Fourier transformation method to get the exact expression of European options.

Finally, we get the relationship of central moments between the physical measure and the risk-neutral measure which can help us to study the negative variance risk premium, the implied volatility smirk, and the prediction of realized skewness. Some relevant work has been done by Bakshi et al. (2003) [13], Carr and Wu (2009) [14], and Neuberger (2012) [15]. However, to the best of our knowledge, except for Zhang et al. (2012) [6], there is no literature studying this relationship. In this paper, we extend it to a stochastic volatility case. This can be regarded as another contribution to the literature.

The rest of the paper is organized as follows. In Section 2, we present our jump-diffusion model with stochastic volatility. In Section 3, we study the equity premium in a general equilibrium economy. The option pricing and the relationship of central moments between the physical measure and the risk-neutral measure are studied in Section 4. Numerical results and conclusions are shown in Sections 5 and 6, respectively.

#### 2. Our Model

In this paper, we consider the financial market with the following two basic assets.(i)A Bond whose price at time is given by (ii)A Stock whose price at time is given by where and ; on the filtered complete space , there are , which both are 1-dimensional Brownian motions with , and is a Poisson process with the constant jump intensity . is the expectation under the physical measure. The jump size is stochastic and , , ,  and   are constants.

The integration form of stock process in (2) is given by

We suppose that the portfolio is which, respectively, means the fraction of wealth invested in the stock and money market; then the wealth process and the initial wealth satisfy the equations as follows: where is the equity premium.

The representative investor maximizes his/her expected utility: where is the conditional expectation and equals .

For tractability, we concentrate our attention on the case of constant relative risk aversion (CRRA) utility function: where the constant is the relative risk aversion coefficient.

The equity premium is very important for option pricing in general equilibrium framework. Following the idea of Santa-Clara and Yan (2004) [12] and Zhang et al. (2012) [6], we obtain the equilibrium equity premium by modeling general equilibrium economy in the following proposition.

Proposition 1. In general equilibrium framework, the equilibrium equity premium is given by where

Proof. From the optimal control problem (5), we get the Bellman equation as follows: where
Equating the derivatives of the Bellman equation (9) with respect to to zero, we have following equation:
In equilibrium, the money market is in zero net supply. Therefore, the representative investor holds all the wealth in the stock market; that is, . Then we can get the expression of from (11): Then Bellman equation (9) can be written as follows: From (6), we conjecture that Then, substituting (11) into (14), we obtain This leads to a system of two ordinary differential equations (ODEs): where
This system (16) and (17) can be solved explicitly. First, we solve the first ODE (16), which is the Riccati differential equation. Making the substitution we obtain the second-order differential equation A general solution has the form where Thus, Then, the solution of the second ODE (17) is Substituting (14) into (12), we get Proposition 1.

Remark 2. In the special case where there is no stochastic volatility and jumps, and , and consequently which is constant in Merton (1976) [2]. In the special case where there is no stochastic volatility, , and consequently which is constant in Zhang et al. (2012) [6].

#### 4. Option Pricing

In this section, we will study the pricing kernel and the option pricing in general equilibrium framework. We first derive the pricing kernel which acts like the physical and risk-neutral densities in the economy and is the key to obtain the PDE of option price as follows.

Proposition 3. In general equilibrium framework, the pricing kernel is given in differential form by and the integration is given by where .
The martingale condition, , requires that the jump size satisfies the following restriction:

Proof. To satisfy the martingale condition, , from (3) and (26), we have Substituting (7) into (28), we have Using the property of Poisson process, , it is easy to obtain Proposition 3.

Remark 4. In this market, there is only one tradable asset, a stock with price , but there are at least two dimensions of risk, diffusive risk, and jump risk. Therefore, the market is incomplete and the pricing kernel is not unique. The nonuniqueness of the pricing kernel can be justified by the fact that the distribution of jump size in the pricing kernel can be arbitrary as long as it satisfies the martingale restriction (27). In a special case, we can choose , as in Liu et al. (2005) [37].

Remark 5. With Proposition 3, we define a new probability measure : since, for any assets at time , we have which means is a risk-neutral probability measure.

Lemma 6. Define a new probability measure, , by the following Radon-Nikodym derivative: then the following relation is true.

Proof. The change of probability measure formula gives
Since , , is i.i.d. and and are correlated, this means that only one of the is correlated with . Without loss of generality, we assume that is correlated with and other are independent of . Then we have

Remark 7. These results are also true in measure, because the difference between and is the Brownian motion that is independent of the jumps.

Now, we consider a European call written on the stock price at time . The option has a payoff function at time . Its price is denoted as at time . We derive a PDE which has to satisfy in the following proposition.

Proposition 8. In general equilibrium framework, the price of European call option satisfies the following PDE: where

Proof. First, we rewrite the stock with continue part and jump part: where Similarly, where where From (40) and (42), we get where The martingale condition requires Using Lemma 6, we have Denoting and using the restriction condition , in Proposition 1 and the terminal payoff function , we obtain If we denote and , then (48) can be written as

Remark 9. The stock process (2) in a risk-neutral measure can be written as where The proof is very easy. Substituting all equations in Remark 9 to (50) and using the equilibrium equity premium (7) and the restriction (27), we will get (2).
Furthermore, we also can understand by Employing the Feynman-Kac theorem to (38), we also can obtain PDE (36) in Proposition 8.

Remark 10. and are similar forms in Heston (1993) [40].

Following the approach of Lewis (2000) [8], we will derive the pricing formula in a risk-neutral measure from the PDE.

Proposition 11. In general equilibrium framework, the pricing formula is given by where

Proof. We denote and ; then, the PDE (36) can be rewritten as Let be Fourier transform of : with .
Denote ; the inverse Fourier transform is given by Then, the PDE (55) can be rewritten as
Denote ; the PDE (58) can be rewritten as
According to the study of Lewis (2000) [8], to solve the PDE (59) with the initial condition (60), it is enough to solve the same equation with the initial value being equal to one. We call this solution fundamental transform and denote it which satisfies following equation by: and the option price satisfies
Denoting , then, the PDE (61) can be rewritten as
We guess that a solution of (63) is as the following form: with and .
This leads to a system of two ODEs: where .
This system can be solved explicitly. First, we solve the first ODE (65), which is the Riccati differential equation. Making the substitution we obtain the second-order differential equation
A general solution of (68) has the form where Thus,
Then, the solution of the second ODE (66) is

Note. Although the pricing formula (40) contains a complex integral, the result is real.

With the European option price being well understood, we study the relations between the physical and risk-neutral densities and moments which can help us to study the negative variance risk premium, the implied volatility smirk, and the prediction of realized skewness. In this paper, we extend the results about those relations in [6] to a stochastic volatility case.

Proposition 12. The risk-neutral density of the continuously compounded return within , , given in (38), is given by The first moment, second, third, and fourth central moments of the continuously compounded return in the risk-neutral measure are given by where ,   is the first moment, , , and are the second, third, and fourth central moments in the risk-neutral measure of random number .
The physical density is given by where The first moment, second, third, and fourth central moments of the continuously compounded return in the physical measure are given by where ,   is the first moment, , , and are second, third, and fourth central moments in the physical measure of random number .

Proof. We denote the corresponding probability function , and also satisfies the PDE (36) with a different boundary condition . Under the Fourier transform, this boundary condition becomes . Solving the PDE (36) under this new boundary condition, we have
We denote the stock return between time and by so that ; then we can differentiate with respect to to get the density function . We find that the parameters of the physical stock price process in (2) are the same those of risk-neutral stock price process in (50); then the parameters of are just removed by the superscript .
Now we compute the first moment, second, third, and fourth central moments of the continuously compounded return in the physical measure. There are at least two methods to solve the problem of how to compute central moments.
One method is using the ’s characteristic function to calculate the th moment of . Since then
Although this method is very straightforward, it is too complicated to calculate the integration. We introduce another way to obtain the central moments by a direct calculation.
From (3), we have
Then, the first moment of the continuously compounded return in the physical measure can be easy to be obtained:
Define the first central moment, second, third, and fourth central moments in the physical measure of random number as , ,  , and .
From (2) we have The integration is given by Then Using (81), we have The results of risk-neutral moments can be obtained with the same procedure in the risk-neutral probability measure.

Remark 13. If , where is constant, our results will degenerate into Proposition 3.19 by Zhang et al. (2012) [6].
We denote , as the skewness and kurtosis in risk-neutral measure, respectively. Skewness and kurtosis are very important for asset pricing. For example, Bakshi et al. (2003) [13] concluded that variations in the risk-neutral skew were instrumental in explaining the differential pricing of individual equity options and found that less negatively skewed stocks have flatter smiles. From Proposition 12, we easily obtain the following three corollaries about skewness and kurtosis.

Corollary 14. The skewness and kurtosis in physical measure are given by And the skewness and kurtosis in risk-neutral measure are given by

Corollary 15. For jump size , the first moment and the second, third, and fourth central moments in the physical measure are given by The first moment and the second, third, and fourth central moments in the risk-neutral are given by where ,

Proof. To compute central moments in the physical measure, we have , , , , and .

Then, from Proposition 12, we have In order to compute central moments in the physical measure, from Proposition 8, we have Then, from Proposition 12, we have

Corollary 16. For small risk aversion coefficient, one has the following relationship between the third central moments in the neutral-risk measure and in the physical measure:

Proof. From Corollary 15, for small risk aversion coefficient , we have

#### 5. Numerical Results

In this section, we use our model to solve the equity premium puzzle. First, we get the long equilibrium risk premium for , where is sufficiently large. In this case, ; ; then

To compare with the results from the model in [6], we take nonrandom constant jump size and ,  , , . Then, we obtain the relative risk aversion coefficient which is the same as the numerical value in [6]. However, the risk premium is stochastic. We get the expectation of as follows:

If we assume that other parameter values in our model are , , , and , we will find that the expectation of tends to 0.6 from Figure 1, which is consistent with above analysis.

Figure 1: The expectation of .

As we know that the risk premium is stochastic and it has the same stochastic characteristic with volatility , the risk premium can be generated by a discrete scheme of system as follows: where is the time interval and is a sample from the standard normal distribution.

One path of is given in Figure 2.

Figure 2: One path of .

In realistic world, the risk premium is stochastic. For example, Figure and Figure 5 in Graham and Harvey (2012) show that risk premium follows mean-reverting stochastic process in [41]. Thus, our model is more realistic than that in [6].

#### 6. Conclusions

We successfully extend the model in [6] and obtain more general results. First, the equilibrium equity premium in general equilibrium framework links not only the jump risk but also the stochastic volatility risk. Moreover, the pricing kernel with stochastic volatility is shown in this paper. We employ the Fourier transformation method to obtain the exact expression of option price. Next, we get the relationship of central moments between the physical measure and risk-neutral measure. It is a vitally crucial problem in financial area. Our numerical results show that our model is more realistic than the one in [6].

#### Acknowledgment

This work is supported by the Fundamental Research Funds for the Central Universities (JBK130401).

#### References

1. F. Black and M. Scholes, “The pricing of options and corporate liabilities,” Journal of Political Economy, vol. 81, no. 3, pp. 637–659, 1973.
2. R. C. Merton, “Option pricing when underlying stock returns are discontinuous,” Journal of Financial Economics, vol. 3, no. 1-2, pp. 125–144, 1976.
3. G. Bakshi, C. Charles, and Z. Chen, “Empirical performance of alternative option pricing models,” Journal of Finance, vol. 52, no. 5, pp. 2003–2049, 1997.
4. R. Cont and P. Tankov, Financial Modelling with Jump Processes, Chapman & Hall, Boca Raton, Fla, USA, 2004.
5. S. G. Kou, “A jump-diffusion model for option pricing,” Management Science, vol. 48, no. 8, pp. 1086–1101, 2002.
6. J. E. Zhang, H. Zhao, and E. C. Chang, “Equilibrium asset and option pricing under jump diffusion,” Mathematical Finance, vol. 22, no. 3, pp. 538–568, 2012.
7. P. Santa-Clara and S. C. Yan, “Volatility, and the equity premium: lessons from S&P 500 options,” The Review of Economics and Statistics, vol. 92, no. 2, pp. 435–451, 2010.
8. A. L. Lewis, Option Valuation under Stochastic Volatility, Finance Press, Newport Beach, Calif, USA, 2000.
9. A. Sepp, “Fourier transform for option pricing under affine jump-diffusions: an overview,” In press.
10. D. Duffie, J. Pan, and K. Singleton, “Transform analysis and asset pricing for affine jump-diffusions,” Econometrica, vol. 68, no. 6, pp. 1343–1376, 2000.
11. G. R. Duffee, “Term premia and interest rate forecasts in affine models,” Journal of Finance, vol. 57, no. 1, pp. 405–443, 2002.
12. P. Santa-Clara and S. Yan, “Jump and volatility risk and risk premia: a new model and lessons from S&P 500 options,” NBER Working Paper w10912, National Bureau of Economic Research, USA, 2004.
13. G. Bakshi, N. Kapadia, and D. Madan, “Stock return characteristics, skew laws, and the differential pricing of individual equity options,” Review of Financial Studies, vol. 16, no. 1, pp. 101–143, 2003.
14. P. Carr and L. Wu, “Variance risk premiums,” Review of Financial Studies, vol. 22, no. 3, pp. 1311–1341, 2009.
15. A. Neuberger, “Realized skewness,” Review of Financial Studies, vol. 25, no. 11, pp. 3423–3455, 2012.
16. H. Föllmer, “Hedging of contingent claims under incomplete information,” in Applied Stochastic Analysis, M. Davis and R. Elliott, Eds., vol. 5 of Stochastics Monographs, pp. 389–414, Gordon and Breach, New York, NY, USA, 1991.
17. X. Ruan, W. Zhu, S. Li, and J. Huang, “Option pricing under risk-minimization criterion in an incomplete market with the finite difference method,” Mathematical Problems in Engineering, vol. 2013, Article ID 165727, 9 pages, 2013.
18. S. I. Boyarchenko and S. Z. Levendorskiǐ, “Option pricing for truncated Lévy processes,” International Journal of Theoretical and Applied Finance, vol. 3, no. 3, pp. 549–552, 2000.
19. M. Montero, “Renewal equations for option pricing,” European Physical Journal B, vol. 65, no. 2, pp. 295–306, 2008.
20. P. Carr, H. Geman, D. B. Madan, and M. Yor, “Pricing options on realized variance,” Finance and Stochastics, vol. 9, no. 4, pp. 453–475, 2005.
21. M. Lorig and O. Lozano-Carbasse, “Exponential Lévy-type models with stochastic volatility and stochastic jump-intensity,” http://arxiv.org/abs/1205.2398.
22. Y. Miyahara, “Canonical martingale measures of incomplete assets markets,” in Proceedings of the Seventh Japan Russia Symposium, Probability Theory and Mathematical Statistics, pp. 343–352, Tokyo, Japan, 1995.
23. Y. Miyahara, “Minimal entropy martingale measures of jump type price processes in incomplete assets markets,” Asia-Pacific Financial Markets, vol. 6, no. 2, pp. 97–113, 1999.
24. T. Chan, “Pricing contingent claims on stocks diven by Lévy processes,” Annals of Applied Probability, vol. 9, no. 2, pp. 504–528, 1999.
25. H. U. Gerber and E. S. W. Shiu, “Actuarial bridges to dynamic hedging and option pricing,” Insurance, vol. 18, no. 3, pp. 183–218, 1996.
26. S. Müller, Arbitrage Pricing of Contingent Claims, Lecture Notes in Economics and Mathematical Systems, Springer, Berlin, Germany, 1985.
27. M. Mania and R. Tevzadze, “A semimartingale backward equation and the variance-optimal martingale measure under general information flow,” SIAM Journal on Control and Optimization, vol. 42, no. 5, pp. 1703–1726, 2003.
28. O. Bobrovnytska and M. Schweizer, “Mean-variance hedging and stochastic control: beyond the Brownian setting,” IEEE, Transactions on Automatic Control, vol. 49, no. 3, pp. 396–408, 2004.
29. S. Hodges and A. Neuberger, “Optimal replication of contingent claims under transactions costs,” Review of Futures Markets, vol. 8, no. 2, pp. 222–239, 1989.
30. I. Karatzas, J. P. Lehoczky, S. E. Shreve, and G. L. Xu, “Martingal