It is well known in the standard Black-Scholes  model that the log-normal stock diffusion with constant volatility is not consistent with the market price movement. More importantly, there are evidences to indicate that the B-S model cannot describe real stock price behavior. Sometimes sudden changes of prices could happen at a random time and these changes cannot be captured by the log-normal distribution characteristic of the stock price in the Black Scholes model. In order to overcome these shortcomings of the Black-Scholes model, a variety of alternative models are proposed in the financial literature. Among these, the jump diffusion models proposed by Merton  and Kou  are widely used models. Furthermore, the SVJ model of Bates can not only make up for the shortcomings of the B-S model but also can describe the financial market in a more suitable way.
In order to reflect the effect of stochastic volatility on the market volatility, we have two ways, one assumes that the volatility is determined, and the volatility function is then determined by calibration to market price. The other assume that the stochastic volatility approach   . In this model the volatility of the stock price is considered to be a mean reverting diffusion process, which is usually related to the stock process itself. A general approach, originally suggested by Bates, can be much more preferable to describe the financial market.
The valuation of option under jump diffusion process satisfies a partial integro-differential equation with boundary conditions. There are several numerical methods available to approximately solve the above equation. For example, in  , Almendral and Osterlee presented an implicit second order accurate time discretization with finite difference. Recently, Patidar  developed an efficient method for pricing Merton jump diffusion option. The scheme proposed by Halluin  required to use an iterative procedure to solve discrete equations. In particular, an approach based on implicit-explicit schemes in which integral term is treated explicitly was proposed by YongHoon Kwon  . This article aims to solve the pricing financial derivatives with jumps which based on time semi-discretization.
The paper is organized as follows. In Section 2, mathematical models for pricing option with jump diffusion process are given in terms of partial integro-differential equations and provide a brief review of both the Merton and Kou jump diffusion models, then the SVJ model of Bates was presented. Section 3 deals with the construction of time level implicit explicit scheme to discretize the jump diffusion model. The discretization method is called IMEX-BDF2 method with three time levels. Section 4 we consider the penalized nonlinear equation to approximate the LCP system. In Section 5, we will give a numerical simulation example. Finally the paper ends with some discussion and conclusive remarks in Section 6.
2. The Mathematical Model
In this part, we briefly discuss the mathematical model for pricing option with jump diffusion process. Consider an asset which the asset price is S, and then the movement of stock price is modeled by the following stochastic differential equation.
where μ is drift rate, t as the time to maturity, σ represents the constant volatility, is generally satisfied with the Gauss process. The constant λ is the intensity of the independent Poisson process, κ used to express the relative size of the expected jump .
The price of risk asset can be described by the Brown movement. However, the Brown movement is subject to normal distribution, and its value can be negative, which does not conform to the nature of the price. According to the geometric Brown motion obeying the lognormal distribution, which avoids the defect that Brown motion may take negative value. The geometric Brown motion has been widely applied when we establish the financial asset price model.
The distribution of the risk asset which obeys the lognormal distribution is symmetric. However, the empirical analysis shows that the price is asymmetric. The log-double-exponential distribution can make up for the shortage of the lognormal distribution. It can simulate the financial market more vividly.
If we define that is probability density function of the jump with amplitude η, under Merton’s model is given by the log-normal density, and under Kou’s jump-diffusion model is the following log-double-exponential density. The function can be written as:
In Merton’s model,
In Kou’s model, where , and is the Heaviside function, we can show that
although the model is well, there are still some shortcomings, there is no consideration that the volatility is also random.
In order to fit the financial market well, it should be noted that the volatility rate σ is a random fluctuation. Then the model of Merton is popularized to get the SVJ (stochastic volatility-jump) model of Bates.
where the are same as the above, represents the instantaneous variance or volatility of the rate of return on assets, which is described by a square root process. It can guarantee the nonnegativity of . ρ is the correlation coefficient of asset price shock and volatility impact on the , reflects the leverage effect of asset price. According to Ito lemma, the SDE equation of the logarithmic asset price satisfaction is:
Let represent the value of a contingent claim that depends on the underlying asset price S with current time t, then satisfy following backward partial integration differential equation.
where and the price of a European put option which the initial condition and asymptotic behavior are described by
In this section, we will discuss the stability of some implicit-explicit (IMEX) time semi-discrete methods. For non-linear partial differential equations, a class of such IMEX methods have already been discussed in  . But now it is entirely different in the sense of analysis and applications. We shall construct an implicit-explicit backward difference method of order two time semi-discretization, which can fit the original equation much more properly.
Let be a partition of the interval [0, T]. We will consider the following three time semi-discretization.
IMEX-BDF2 (Implicit-explicit backward difference method of order two)
where the , , and , In order to prove the stability of the time semi-discretization, we will need the following lemmas.
Lemma 1: (Discrete Gronwall’s inequality) Let and be four non-negative sequences, such that for non-negative integers & N
then, , where the summation the is assumed to be zero if . This lemma can be easily proved by mathematical induction.
The Discrete Gronwall’s inequality illustrates that the function which satisfies the integro-differential equation, and then there is a corresponding inequality.
Lemma 2: For IMEX-BDF2, there exists an such that & , we have
Proof: From relation IMEX-BDF2 we see that, ,
which implies that
By taking and using symmetry of , we have
We can get
Now by using Lemma1, we have
for some generic constants .
After summing from to , , we get
Hence there exists an such that above relation with the use of Discrete Gronwall’s inequality gives
Further, we have
By taking and using the relation:
Again by using Lemma 1, we can get
for some generic constants .
After summing from to , ,
Thus there exists an such that , above relation with the use of Lemma 1 gives
Finally with the help of this above relation, we have
Now Lemma1, Lemma2, imply the following result: under the assumptions made in above, all the three time semi-discrete IMEX methods are stable.
4. Unconditional Stability
The types of initial and boundary conditions can be properly defined as different types of options. It is different from that the American option can be exercised at any time up to the maturity date. The types of initial and boundary conditions can not be defined easily, but it can be formulated as the linear complementarity problem (LCP) of the form.
In combination with the above equation and the linear complementarity problem, we can make a transformation to simplify the operation.
where K is the strike price of the options, for all , and are the lower and upper truncated boundary respectively.
Since fractional derivative is nonlocal operator, the truncation error cannot be ignored for the nonhomogeneous Dirichlet boundary condition. The boundary condition of FPDEs can be transformed to homogeneous Dirichlet boundary condition.
Let is defined by
using the new variable , then we have
where , the boundary and terminal conditions will become
We have discussed that American options are described by the linear complementarity problem. Then we consider the following penalized nonlinear equation to approximate the LCP system (4.4),
where , and (4.5) converges to (4.2) as ρ goes to infinity. A semi-implicit finite difference scheme is proposed to discretize the nonlinear FPDEs system (4.5).
Now divide the intervals and into N + 1 subintervals and M subintervals, respectively. The spatial and temporal meshes are defined as follows
where , are positive integers and ρ is chosen such that .
Denoting and , we obtain the following finite difference scheme for the nonlinear FPDE system
Lemma 3: The nonlinear scheme above is unconditionally stable.
Proof: Let be the solution of above in the m-th time step. Let be the i0-th entry of such that
thus, . Then, the proof will be completed by discussing the following two cases.
Case 1: for , in view of above, we have
Case 2: For suﬃciently small such that ，we have ,
and C is generic constant depending on the parameter and T.
This case can be proved (see Theorem1 in  ).
5. A Numerical Simulation Example
In this section, we present several numerical experiments to illustrate the efficiency and accuracy of proposed method.
Example: We will simulate the HS300 to fit the SVJ model. The sample number is 10,000. For the convenience of simulation, we will eliminate half of the abnormal data. The stability and effectiveness of time-discretization are verified by the simulation of the remaining data. All numerical experiments result at strike price with parameters as provided in Table 1.
The s-d and MC-error indicates the error of the simulation. In Table 1, we can find that the data are so small that it can be ignored.
Finally, through the above data simulation, we proved that the proposed method is much more efficient and accurate.
6. Discussion and Conclusion
In this article, it has proposed and analyzed three implicit-explicit (IMEX) time
Table 1. Numerical simulation analysis.
semi-discretization namely, IMEX-BDF2, for solving partial integro-differential equations which arise in option pricing theory when the underlying asset follows a jump diffusion process. All the IMEX time semi-discretization is shown to be stable. The American option whose types of initial and boundary conditions cannot be defined easily, but it can be formulated as the linear complementarity problem (LCP). The numerical simulation with European put/call under SVJ model has been carried out.
 Heston, S.L. (1993) A Closed-Form Solution for Options with Stochastic Volatility with Applications to Bond and Currency Options. Review of Financial & Studies, 6, 327-343.
 Ngounda, E., Patidar, K.C. and Pindza, E. (2013) Contour Integral Method for European Options with Jumps. Communications in Nonlinear Science & Numerical Simulation, 18, 478-492.
 Kwon, Y. and Lee, Y. (2011) A Second-Order Finite Difference Method for Option Pricing under Jump-Diffusion Models. SIAM Journal on Numerical Analysis, 49, 2598-2617.
 Ascher, U.M., Ruuth, S.J. and Spiteri, R.J. (1997) Implicit-Explicit Runge-Kutta Methods for Time-Dependent Partial Differential Equations. Elsevier Science Publishers B. V.
 Kadalbajoo, M.K., Kumar, A. and Tripathi, L.P. (2016) A Radial Basis Function Based Implicit-Explicit Method for Option Pricing under Jump-Diffusion Models. Applied Numerical Mathematics, 110, 159-173.