Abstract

The empirical research shows that the log-return of stock price in finance market rejects the normal distribution and admits a subclass of the asymmetric distribution. Hence, the pricing problem of stock loan is investigated under the assumption that the log-return of stock price follows the CGMY process in this work. Under this framework, the pricing model of stock loan can be described by a free boundary condition problem of space-fractional partial differential equation (FPDE). First of all, in order to change the original model defined in a fixed domain, a penalty term is introduced, and then a first order fully implicit difference schemes is developed. Secondly, based on the numerical scheme, we prove the stock loan value generated by our method does not fall below the value obtained when the contract of stock loan is exercised early. Finally, the numerical experiments are implemented and the impacts of key parameters in the CGMY model on the value and optimal redemption price of stock loan are analyzed, and some reasonable explanation should be given.

1. Introduction

A stock loan can be treated as a contract between two parties: the bank or other financial institution (the lender), and a client (the borrower). The borrower obtains a loan from the lender with their stocks as collateral, and the contract offers the borrower a right to redeem the stock at any valid time. But he/she should pay off the loan as well as the cumulative loan. As described in Ref. [1], the risk aversion investors can use stock loan to transfer the risk of holding stock to the financial institutions, and this kind of derivatives also can establish market liquidity. Therefore, the stock loan is one of the most important financial derivatives in the financial market and currently, both Florida Mortgage Corporation and Shelly Bay Capital specialize in providing stock loan services.

The pricing problem about the stock loan has been a popular topic in the academic since the first research paper [1] about the stock loan published in 2007. In this literature, the risk asset was assumed to follow the Geometric Brownian motion (GBM) and the contract of stock loan is finite maturity. Following this publication, more and more researchers paid attentions to the academical topic of stock loan valuation with other conditions. For instance, Lu and Putri [2] considered the pricing problem of stock loan with finite maturity margin. Under the framework of hyper-exponential jump diffusion model, Cai and Sun [3] studied the value and optimal redemption price of stock loan with infinite and finite maturity. The stock loan with finite maturity was also investigated in Ref. [4] under the case that risk-free interest rate follows the Rednleman–Bartter model without drift term.

It is well known that the classical Black–Scholes (B–S) model was established under a lot of strict assumptions, such as the log-return of stock price follows the nonmoral distribution, frictionless, and so on. However, these strict assumptions do not accord with the dynamic process of stock price in real financial market. The empirical research shows that discontinuities or jumps are believed to be an indispensable element of financial risk-asset price (see e.g., [57]) and the log-return of risk asset appear to be “asymmetric distribution” and “leptokurtic distribution” (see e.g., [812]). For this reason, for capturing these characters of stock price, many scholars use other complex stochastic process to drive the stock price. Prominent examples including the CEV model [13], GEV model [14], KoBoL model [15], and so on. A frequently used stochastic process is called CGMY process and it was presented by Carr et al. [16] with the aim to provide a model for the dynamic of equity log-returns. This model is rich enough to accommodate jumps of finite or infinite activity, and finite or infinite variation. In fact, this stochastic process is a particular type of pure jump Lévy process with four key parameters (e.g., , , , and ) which control its essential characteristics.The parameter may be viewed as a measure of the overall level of activity. The aggregate activity level should be calibrated through movements under the case of keeping the other three parameters constant and integrating over all that moves exceeding a small level. Parameters and control the rate of exponential decay on the right and left of the Lévy density, respectively, leading to skewed distributions when they are unequal. Parameter is particularly useful in characterizing the fine structure of the stochastic process.The CGMY process has infinite variation and finite quadratic variation if . At present, many papers discussed pricing problem of financial derivatives based on the CGMY process. Ballotta and Kyriakou [17] obtained the option price by using the Monte Carlo method in this case. Zhang et al. [18] presented a fast numerical method to double barrier option under framework of CGMY. Chen et al. [19] investigated European option and they obtained the explicit closed-form analytical based on the model.

A better refernce value to the holder of stock loan can be provided in this case of CGMY model. The pricing model is free-boundary problem of partial differential equation with tempered fractional derivatives. The tempered fractional derivatives and nonlinearity associated with the free boundary add the difficulty of solving the pricing model. Aiming at this problem, we consider a penalty method in which the free boundary is removed by adding a small and continuous penalty term to the governing equation. Therefore, the mathematical model should be solved on a fixed domain. Even if on the fixed domain, the analytical solution is seldom available so that a full-implicit scheme is employed. Based on our numerical scheme, we prove that the stock loan value generated by the penalty method cannot fall below the value obtained when the stock loan is exercised early, and we also do simulation to verify the result.

This study has five main sections. Section 1 provides the introduction. The pricing model of stock loan is set in Section 2 and the numerical method is introduced in Section 3. The simulation and discussions are presented in Section 4, and conclusion is displayed in Section 5.

2. Mathematical Model

In this section, the governing equation of the stock loan with finite maturity, which is a partial differential equation with fractional derivative, is presented under the CGMY framework. And then, financially, the corresponding free moving boundary and terminal conditions will be given to complete the pricing model.

The CGMY process assumes the log value of stock price follows a stochastic differential equation under the risk-neutral measure [16]

where , , and are the risk free interest rate, the dividend and the current time, respectively. is a stochastic variable and its characteristic function is

In addition,

where the parameters and are constants, and , .

2.1. Governing Equation

In fact, as described in Ref. [1] a stock loan problem can be regarded as an American call option with a time-varying strike price , where is the principal and is the continuously compounded loan interest rate and in general. Therefore, the payoff function of stock loan at maturity can be written as

In addition, typically as same as the American call, at each valid time there is a particular value of the underlying, which yields a boundary between two regions: in one side(called exercising region) the investor exercise the contract of stock loan and in other side (called holding region) one should hold this contract. The particular value of the asset is called optimal redemption price and denoted by . Then, based on the principle of non-arbitrage pricing, in the holding region, the value of stock loan at time satisfies

where is the conditional mean operator under the measure and denotes the information flow at time . Famous authors, Cartea [20] and Del-Castillo-Negrete obtained the governing equation of through the Fourier transformation as follows:

where , and

are left and right sided Riemann–Liouville fractional derivatives.

In fact, the -order left and right fractional derivatives are closely related to CGMY process. The fractional differentiation is nonlocal, which relates to the stock loan value in the stopping region . This is where the stock loan value is equal to its intrinsic value. The nonlocalness of this fractional operator means that over a time step , the price of risk-asset should diffuse to value far away from , provides a way in efficiently simulating the existence of jumps of the risk-asset.

2.2. The Boundary Conditions

In this subsection, to complete the pricing system, a set of appropriate boundary conditions will be given. First of all, as a rational investor, she/he does not redeem stock if the level of stock price is very low and here we assume the stock price is close to zero, i.e., , then one can obtain a boundary condition as

On the other hand, the continuity of stock loan value at optimal redemption price should still be retained to ensure the smooth pasty of the stock loan value across the free boundary. For these reasons, we still impose the two boundary conditions as follows:

And the terminal condition is the payoff function

To sum up, the complete pricing model for stock loan with finite maturity under CGMY process can be written as

It must be noted that, a rational investor cannot redeem the stock if the value of redemption is less than the holding value; therefore, the value of the stock loan should satisfy the following inequality

for all and .

2.3. Model Normalization

It must be noted that there is a time-factor in the boundary conditions (14) and (17). And we find the time-factor should influence numerical results, hence, we introduce a new variable system to the model (14)–(17) for avoiding the effect of . Take

then

let , therefore we obtain

And

let , so one obtains

Now, substituting Eq. (20), (21), (23), (25) into Eq. (13) and after transforming the boundary conditions (14)–(17) through Eq. (19), the PDE system of function can be obtained as follows:

where ,

And the equation (18) should be changed as

Mathematically, the function in the model (26)–(30) can be viewed as an American call option with strike price and free-boundary . According to the property of American call, the values of are not less than the payoff . And the free-boundary is monotonically increasing with respect to time to expiry of option contract. In the next section, a numerical scheme based on the penalty method will be proposed to solve this model.

3. Numerical Method

In this section, we first consider a penalty approach in which the free moving boundary is removed by adding a small and continuous penalty term, so that the stock loan pricing problem can be solved on a fixed domain. Then, a different scheme is proposed and an effective nonlinear Newton-type iteration strategy is employed. Furthermore, due to the coefficient matrixes of the finally linear system, which contains the full matrix with Toeplitz structure, the fast biconjugate gradient stabilized method (FBi-CGSTAB) is used to solve our system.

3.1. Model Transformation

In this paper the penalty function is defined as

where is a regularization constant and , is a constant,

If adding the penalty term (34) to the FPDE in (26), one can obtain a nonlinear FPDE system defined on a fixed domain as follows:

where denotes the maximum stock price, . And according to the results in Ref. [21], the maximum value of stock is equal to 3 or 4 times of strike price, hence we take in this paper. In addition, we will omit the subscript of function for convenience.

3.2. Difference Scheme

For the FPDE of (37) is nonlinear, it results in the analytical solution of the mathematical model (37)–(39) that is laborious and even impossible to achieve even if the model is defined on the fixed domain . Therefore, an effective numerical algorithm should be preferred.

Specifically, taking as spatial step such that , where is a positive integer, and placing uniform grids in the time direction, namely , that is

and . For the first order spatial and time, we use the following difference scheme

And the fractional derivative can be approximated by the first-order Grnwald–Letnikov formula as [22]

where is the value of function at grid point , are fractional difference coefficients

and the parameters satisfy the following system

If taking as the free variable, the solution of system (46) can be obtained as

So, the fully implicit difference scheme without the truncation error for (36) is got as follows:

with the boundary and terminal conditions

where

For our proposed difference scheme, we have a discrete analogue form of the important property inherited by the stock loan model. Prior to presenting the proof, we need the following two lemmas.

Lemma 1. If the parameter , then in the Eq. (50) are more than zero.

The proof process of Lemma 1 is provided in Appendix A.1

Lemma 2 (see [23]). For , the coefficients satisfywhere .

In fact, according to Lemma 2, we obtain the following equation:

where .

Theorem 3. If the time step andthen the approximate stock loan values generated by the scheme Eq. (38) satisfyfor all , , whereFor details of proof process, please refer to Appendix A.2.

4. Numerical Examples and Discussions

In this part, the numerical simulation should be executed for verifying the theoretical consideration of our proposed difference scheme. The system (48) is nonlinear, hence a practical Newton-type iteration strategy is given to realize the nonlinear difference scheme. The special structure preserved by the coefficient matrix resulted by difference scheme is carefully exploited. Furthermore, the fractional derivatives usually result in a dense coefficient matrix in the system. It has significant computational and storage requirements. In terms of computational cost; it is very important to use an effective and efficient method to solve our linear system. Therefore, the fast biconjugate gradient stabilized method (FBi-CGSTAB) [18] is employed. Finally, the impacts of important parameters to the optimal redemption price should be analyzed.

Practically, to implement the numerical scheme (48) in computer, the semi-infinite domain must be truncated into a finite domain:

here we take , and now . In this case, we should redefine the space step , and

If we let

then the matrix of numerical scheme (48) can be written as:

where

and denotes the identity matrix, means the transposition of matrix . Both and are Toeplitz matrix

In fact, the system (60) is not solved directly for the penalty function is nonlinear with respect to , so this nonlinear system can be solved through Newton-type iteration approach, which is provided in Appendix A.3.

4.1. Numerical Testing

To ensure that both the theoretical model and method are feasible, many results in our work must be verified before quantitative analysis; therefore, in this part we will carry out some numerical testing. First of all, the fact that our proposed numerical algorithm satisfies the discrete analogue of the positive constraint, in Theorem 3 will be verified. As shown in Figures 1 and 2, the mesh surface of for the different time and stock price under various parameter setting is presented. And the mesh surface shows that the present difference scheme conserves the inequality for all .

As described in Section 2.3, the function in the model (36)–(39) can be recognized as an American call option with strike and optimal exercise boundary . Then, mathematically, the values of are not less than payoff function and is increasing with respect to time to expiry. Figures 3 and 4 show the two facts, respectively.

4.2. Impact of Parameters

As described in Section 2, the essential characteristics of CGMY model are controlled by four key parameters, namely, and and. In other words these parameters should affect the value and optimal redemption price of stock loan. Hence, in this subsection we employ the proposed numerical method to investigate the impacts of the four parameters on the stock loan value and capture the optimal redemption price for different parameters setting with some reasonable explanation.

Figures 5 and 6 display the behaviours of stock loan valuation and optimal redemption price for different , respectively. As shown in Figure 5, the mesh surf of is higher with the increase of . Financially, this phenomenon could be explained from the fact that except keeping the other parameters constant, the aggregate activity level of CGMY process may be calibrated through movements in and the level is increasing with respect to . In addition, the stock loan and American call have a similar property that the higher activity level of underlying asset results in the higher value. Hence, it is reasonable to suggest that the stock loan value mesh surface should move upwards as becomes larger. Mathematically, in the exercising region where the payoff must be equal to the stock loan value and the optimal redemption price is the boundary of exercising region, the higher stock loan values should yield the bigger optimal redemption price, namely, a bigger shall yield a higher optimal redemption boundary shown in Figure 6.

Theoretically, the two parameters and control the rate of exponential decay on the right and left of the density of CGMY process,respectively. As for , the left (right) tail of the distribution for this process is heavier than the right (left) tail, which is consistent with the neutral distribution typically implied from stock loan values [16]. Also a smaller value of results in a heavier left (right) tail. Moreover, according to the results in Ref. [24], it is suggested that the stock loan should have a positive value only if there is a large decrease in the risk-asset. Hence, the stock loan value relies on the left (right) tail of the risk-neutral distribution of stock. In other words, the fatter the left (right) tail is, the more valuable the stock loan value. Therefore, as shown in Figures 7 and 8, the mesh surface of is decreasing with respect to and during with other parameters given, as depicted in the complete time dimension.

According to the analysis in explained earlier, it is not difficult to obtain that a smaller value of or would yield a higher stock loan value, hence, as a rational investor shehe should raise the redemption price at any valid time when there are smaller value of or . And in other words, the higher optimal redemption boundary is increasing with respect to and at any valid time, as shown in Figures 9 and 10.

Finally, we should focus on the impact of parameter . Carr [16] claims that parameter decides whether the up jumps and down jumps of CGMY process have a completely monotone Lévy density. Moreover, both activity and variation of this process will become larger as becomes bigger. So, the investors must be willing to obtain a higher value of contract under the case of “out-the-money”, then the mesh-surface is increasing with respect to with other parameters are kept fixed as shown in Figure 11. However, when stock prices are less than principal, our numerical results are consistent with the conclusion of call option calculated by the fat tails distribution under the “in-the-money” case [14]. As for the stock loan contract, the stock prices must be larger than principal in the exercising region where the value of stock loan decreases monotonically with respect to . This leads to a fact that the bigger the value of is, the higher the optimal redemption boundary is, as depicted in Figure 12.

5. Conclusions

In this paper, the stock loan valuation under the CGMY framework is investigated. In this case, the pricing mathematical model is a FPDF free-boundary problem. A nonlinear penalty term is introduced to transform the free-boundary model into one with fixed domain. Then, a full nonlinear implicit numerical scheme is proposed and we also proved that the numerical solution of transformed model satisfied the inequality . In addition, to improve efficiency of computation, the FBi-CGSTAB method was employed.

In numerical simulation, the impacts of key parameters , and on the stock loan value and optimal redemption price are analyzed. Parameter decides activity level of underlying asset so that both stock loan value and optimal redemption price are increasing with respect to . For the two parameters and , they control the rate of exponential decay on the right and left of the Lévy density, respectively; in other words, the bigger value of results in the fatter left (right) tail. For this reason, both mesh surface of and optimal redemption boundary decrease with respect to and . For parameter , in the case of “out-of-the-moneyness”, the bigger can yield the higher mesh surface and optimal redemption price. However, this phenomenon would disappear if stock price became less than the strike price.

Appendix

A.1. Proof Process of Lemma 1

Proof. We prove . Substituting (47) into the formula of , one can obtainDue to , so . Moreover,and , therefore if , then . Similarly, we also prove under the condition of .

A.2. Proof Process of Theorem 3

Proof. Without loss of generality, we prove this theorem in two steps. The scheme (48) can be rewritten asIn order to prove , for all , we introduceand it is straightforward to obtain . Hence, by substituting into (A.3), it yieldsAfter simplification one can obtainwhereSince,we havewhereIn addition, from the fact , and according to [15], when , , we getAnd according to Eq. (50) one can obtainFinally, the following inequality can be obtainedwhere .
Definetake be an index such that . Noticing the fact that and , thereforeand for , for , it follows from (A.3) thatAfter simple algebra operation, one obtainsDefine a functionwhere . So if we assume that , directly, we getAccording the Eq. (53), we haveIn additionbased on the condition , we can obtain , naturally . Consequentlyfor all .
Next, we should prove that for all . Following the above idea, we defineand let be an index such that . It follows from (48) thatAfter arrangement, it givesIn first step, for all have been proved, then , soNoticing the formula in the square brackets in the above equation is more than 0, and for all , namely . Hence, by mathematical induction, one obtainsfor all , which completes the proof.

A.3. Newton-Type Iteration Approach to the System (48)

where , with the initial value for each time level as the given initial guess and . is Jacobian matrix of column vector and is a damping parameter. We choose , if for some as the stopping criterion, where is a sufficiently small positive control tolerance number. In this paper, we take and . Now the challenging point that should be emphasized is that both matrices and are dense matrix, and the storage requirement and computational efforts are very high, which presents difficulty in capturing the optimal exercise boundary under the CGMY framework. So, the FBi-CGSTAB method should be employed to overcome the challenging point. And finally, the total storage requirement and computational costs have been significantly reduced from to and from to respectively.

Data Availability

There isn’t any data in our manuscript, but for the Matlab code, we can provide to anyone who want.

Conflicts of Interest

The author declare that they have no conflicts of interest.

Acknowledgments

We are especially grateful for anonymous referees for comments and suggestions that greatly improved this paper. The work on this paper was partially supported by the Guizhou science and Technology Department Basic Research Projects ([2019]1175); Science and Technology Bureau Project of Baiyun District ([2017]37); PhD Research Fund of Guizhou University of Commerce ([2018]001).