Next Article in Journal
Threshold Dynamics and the Density Function of the Stochastic Coronavirus Epidemic Model
Previous Article in Journal
Linearly Implicit High-Order Exponential Integrators Conservative Runge–Kutta Schemes for the Fractional Schrödinger Equation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Barrier Option Pricing in the Sub-Mixed Fractional Brownian Motion with Jump Environment

Department of Mathematics, Zhejiang University of Science and Technology, Hangzhou 310023, China
*
Author to whom correspondence should be addressed.
Fractal Fract. 2022, 6(5), 244; https://doi.org/10.3390/fractalfract6050244
Submission received: 24 March 2022 / Revised: 21 April 2022 / Accepted: 22 April 2022 / Published: 29 April 2022

Abstract

:
This paper investigates the pricing formula for barrier options where the underlying asset is driven by the sub-mixed fractional Brownian motion with jump. By applying the corresponding I t o ^ ’s formula, the B-S type PDE is derived by a self-financing strategy. Furthermore, the explicit pricing formula for barrier options is obtained through converting the PDE to the Cauchy problem. Numerical experiments are conducted to test the impact of the barrier price, the Hurst index, the jump intensity and the volatility on the value of barrier option, respectively.

1. Introduction

Barrier option is a path-dependent exotic option, whose value depends not only on the price of the underlying asset, but also on whether the price of the underlying asset touches the preset barrier price within the effective execution period of the option. For its cheaper premiums against the corresponding vanilla options, barrier options can be seen everywhere in global exchanges and over-the-counter markets. Many companies use various barrier options to hedge risks. In addition, the studies on barrier option pricing can also promote the research of many structured financial products, such as convertible bonds, bank-triggered financial products and so on. For the above reasons, the pricing of barrier options has always been a topical issue [1,2,3,4]. If the option right terminates (starts) when the underlying asset price touches the given barrier price, it is called a knock-out (in) kind; it is called a down (up) option, if the initial underlying asset price is above (below) the barrier price [5]. Therefore, single-barrier options which this paper discusses include eight types: down (up)-and-out (in) call (put) options.
In 1973, Merton [6] gave the closed solution of down-and-out European call options. Later, Reiner and Rubinstein [7] extended the pricing formulas of other European barrier options in 1991. However, these studies are under the Black–Scholes model [8] (the B-S model) which assumes the underlying asset price follows the logarithmic normal distribution. However, in recent years the self-similarity and long-range dependence has been found in the financial asset through numbers of the financial empirical studies [9,10], which is inconsistent with the B-S model. Then, Necula [11] studied the extended B-S model, where the assets price is driven by the fractional Brownian motion (fBm) instead of the Brownian motion. The fBm was first proposed by Kolmogorov [12], which exhibits self-similarity and long-range dependence. Since then, a volume of research on option-pricing models with fBm have been conducted, such as [13,14,15].
However, the fBm is neither a Markov process nor a semi-martingale, except degenerating into the Brownian motion. Although we can use Wick-self-financing strategies to analyze the fBm [16,17], Björk and Hult [18] found the application of the fBm has little economic sense, which limited its applicability in financial market. Therefore, other processes are proposed to describe the fluctuation of financial assets, such as the sub-fractional Brownian motion (sub-fBm) [19] and the sub-mixed fractional Brownian motion (sub-mixed fBm) [20].
The sub-fBm preserves most properties of the fBm, but it has the characteristics of non-stationary second-order moment increment and faster convergence [21]. Moreover, the sub-mixed fBm is a combination of the Brownian motion and the sub-fBm. When the Hurst index H 0.75 , 1 , the sub-mixed fBm becomes a semi martingale, which is equivalent to the Brownian motion [22]. At the same time, inspired by Merton [23] and other recent research [24,25,26], this paper introduces the jump diffusion process to describe the jump points of asset price caused by unsystematic risk factors, which is usually ignored in the pricing of barrier options. The purpose of this paper is to obtain the pricing formula of barrier options where the underlying asset is driven by the sub-mixed fBm and the compensated Poisson process.
The remainder of this paper is organized as follows: In Section 2, some necessary preliminary knowledge about the sub-fBm will be presented. In Section 3, we obtain the corresponding I t o ^ ’s formula of the asset price driven by the sub-mixed fBm with jump, and give the expressions for underlying asset price. In Section 4, the Black–Scholes PDE and the closed-form solution for barrier options are obtained. In Section 5, numerical experiments are carried out to study the influences of several parameters on barrier options. Section 6 gives a summary.

2. Preliminaries

Let { Ω , F t , P } be a complete probability space with a filtration { F t } t 0 satisfying the usual conditions.
Definition 1. 
The sub-mixed fBm ξ t H = ξ t H ( α , β ) t 0 is a linear combination of the Brownian motion B t t 0 and the sub-fBm B t H t 0 , which can be expressed as:
ξ t H ( α , β ) = α B t + β B t H , t 0 ,
where H is the Hurst index, α a n d β are positive constants, B t t 0 and B t H t 0 are independent of each other.
Lemma 1. 
The sub-mixed fBm ξ t H = ξ t H ( α , β ) t 0 has the following properties [20]:
1. 
ξ t H ( α , β ) t 0 is a central Gaussian process.
2. 
When t = 0 , ξ 0 H ( α , β ) = α B 0 + β B 0 H = 0 .
3. 
t , s 0 , the covariance of ξ t H ( α , β ) and ξ s H ( α , β ) is
Cov ξ t H ( α , β ) , ξ s H ( α , β ) = α 2 ( t s ) + β 2 2 t 2 H + s 2 H | t s | 2 H ,
where t s = 1 2 ( t + s | t s | ) .
4. 
t 0 , E ξ t H ( α , β ) 2 = α 2 t + β 2 2 2 2 H 1 t 2 H .

3. Asset Pricing Model

In this paper, we adopt the classical financial stochastic analysis theory and make some extensions for the B-S model. Furthermore, the following assumptions are hold:
  • There are two kinds of assets in the financial market: risk-free assets (bonds) and risky assets (stocks).
  • The stock price S t is driven by the sub-mixed fBm with jump:
    d S t = ( μ q ) S t d t + S t d ξ t H ( α , β ) + γ S t d J t = ( μ q ) S t d t + α S t d B t + β S t d B t H + γ S t d J t ,
    where μ is the instantaneous expected return rate of the stock; q is the stock dividend rate; α , β and γ represent the volatility of stock price; J t t 0 is a compensated Poisson process with intensity λ . B t t 0 , B t H t 0 and J t t 0 are independent of each other.
  • The return of risk-free assets in time period t is
    d M t = r M t d t ,
    where constant r is the risk-free interest rate.
  • All assets can be traded freely and continuously without transaction costs and taxes.
  • There is no arbitrage opportunity in the market.
  • Short selling is not limited.
  • The option can be exercised only at the maturity time.
Remark 1. 
It is worth mentioning that there are some limitations when the B-S type model is applied, which are detailed in references [27,28,29,30]. In this paper, we focus on the classical setting of B-S model and will not elaborate too much here. If possible, further research can be carried out in the future.
Theorem 1. 
Assume that Y t = ξ t H ( α , β ) + γ J t with the initial value zero, and f t , Y t is second-order differentiable. Then, the I t o ^ ’s formula of the sub-mixed fBm with jump can be expressed as follows:
f t , Y t = f ( 0 , 0 ) + 0 t f s λ γ f Y d s + 0 t α 2 2 + 2 2 2 H 1 H β 2 s 2 H 1 2 f Y 2 d s + α 0 t f Y d B s + β 0 t f Y d B s H + s t f s , Y s f s , Y s = f ( 0 , 0 ) + 0 t f s + α 2 2 + λ γ 2 2 + 2 2 2 H 1 H β 2 s 2 H 1 2 f Y 2 d s + α 0 t f Y d B s + β 0 t f Y d B s H + γ 0 t f Y d J s .
Proof. 
According to the sub-mixed fBm I t o ^ ’s formula [20] and the jump process analysis method [31], we have
f t , Y t = f ( 0 , 0 ) + 0 t f s d s + 0 t f S d Y s c + 1 2 0 t 2 f S 2 ( d Y s c ) 2 + s t f s , Y s f s , Y s = f ( 0 , 0 ) + 0 t f s λ γ f Y d s + 0 t α 2 2 + 2 2 2 H 1 H β 2 s 2 H 1 2 f Y 2 d s + α 0 t f Y d B s + β 0 t f Y d B s H + s t f s , Y s f s , Y s .
The following identities are used:
d Y t c = α d B t + β d B t H λ γ d t ,
( d Y t c ) 2 = α 2 + 2 2 2 2 H 1 H β 2 t 2 H 1 d t
and Y t c = α B t + β B t H λ γ t is the continuous part of Y t .
If g ( x ) is second order differentiable. Given that Poisson process N t t 0 with intensity λ has the second-order moment increments < d N t , d N t > = λ d t , by generalized I t o ^ ’s formula we obtain
s t g N s g N s = 0 t g N d N s + λ 2 0 t 2 g N 2 d s .
Combining Y t = ξ t H ( α , β ) + γ J t = α B t + β B t H + γ N t λ γ t , we arrive at
s t f s , Y s f s , Y s = γ 0 t f Y d N s + λ γ 2 2 0 t 2 f Y 2 d s .
Substitute (4) back to (3), which yields
f t , Y t = f ( 0 , 0 ) + 0 t f s + α 2 2 + λ γ 2 2 + 2 2 2 H 1 H β 2 s 2 H 1 2 f Y 2 d s + α 0 t f Y d B s + β 0 t f Y d B s H + γ 0 t f Y d J s .
Theorem 2. 
The stock price satisfying (1) has the following explicit solution:
S t = S 0 exp ( μ q ) t α 2 2 + λ γ 2 2 t + 1 2 2 H 2 β 2 t 2 H + α B t + β B t H + γ J t .
Proof. 
Let f t , Y t = S 0 exp ( μ q ) t α 2 2 + λ γ 2 2 t + 1 2 2 H 2 β 2 t 2 H + Y t .
An application of Theorem 1 yields
d f t , Y t = f t + α 2 2 + λ γ 2 2 + 2 2 2 H 1 H β 2 t 2 H 1 2 f Y 2 d t + f Y d Y t . = ( μ q ) f t , Y t d t + f t , Y t d Y t = ( μ q ) f t , Y t d t + f t , Y t d ξ t H ( α , β ) + γ f t , Y t d J t ,
where
f t = ( μ q ) α 2 2 + λ γ 2 2 + 2 2 2 H 1 H β 2 t 2 H 1 f t , Y t , f Y = f t , Y t a n d 2 f Y 2 = f t , Y t .
Comparing (1) and (5), we can deduce d S t = d f ( t , Y t ) , where the values are the same f ( 0 , Y 0 ) = S 0 . Therefore,
S t = f t , Y t = S 0 exp ( μ q ) t α 2 2 + λ γ 2 2 t + 1 2 2 H 2 β 2 t 2 H + α B t + β B t H + γ J t .

4. Pricing Formula for Barrier Options

With the explicit solution of the stock price S t in hand, in this section the pricing formula for battier options can be derived.
Theorem 3. 
Assuming that the underlying asset price S t follows (1), then the value of contingent claims V t = V ( t , S t ) satisfies the following PDE:
V t + ( r q ) S t V S + α 2 2 + λ γ 2 2 + 2 2 2 H 1 H β 2 t 2 H 1 S t 2 2 V S 2 r V t = 0 .
Proof. 
Using the self-financing strategy θ t = θ t 1 , θ t 2 , we hold a number of θ t 1 bonds and θ t 2 stocks to build the wealth process, whose value at time t is
V t = θ t 1 M t + θ t 2 S t .
Combining (1) and (2), we obtain
d V t = θ t 1 d M t + θ t 2 d S t + θ t 2 q S t d t = r θ t 1 M t + μ θ t 2 S t d t + θ t 2 S t α d B t + β d B t H + γ d J t .
At the same time, by applying Theorems 1 and 2, we have
d V t = V t d t + V S d S t + 1 2 2 V S 2 d S t 2 = V t + ( μ q ) S t V S + α 2 2 + λ γ 2 2 + 2 2 2 H 1 H β 2 t 2 H 1 S t 2 2 V S 2 d t + S t V S α d B t + β d B t H + γ d J t ,
where d S t 2 = S t 2 α 2 + λ γ 2 + 2 2 2 2 H 1 H β 2 t 2 H 1 d t .
Comparing (7) and (8), θ t 1 and θ t 2 are given
θ t 1 = r M t 1 V t q S t V S + α 2 2 + λ γ 2 2 + 2 2 2 H 1 H β 2 t 2 H 1 S t 2 2 V S 2 , θ t 2 = V S .
From (6), we obtain
θ t 1 = V t θ t 2 S t M t .
Combining (9) and (10), Theorem 3 is proved. □
Theorem 4. 
Suppose that the underlying asset price S t satisfies (1), then at time t the value of the down-and-out call option C d o t , S t with the fixed strike price K, the fixed barrier L and the maturity time T is given
C d o t , S t = S t e q ( T t ) N d 1 K e r ( T t ) N d 2 S t L κ ( t ) L 2 S t e q ( T t ) N d 3 K e r ( T t ) N d 4 .
The following identities are used: N ( x ) = 1 2 π x e t 2 2 d t , which denotes the cumulative probability of standard normal distribution;
d 1 = ln S t K + r q + α 2 2 + λ γ 2 2 ( T t ) + 1 2 2 H 2 β 2 T 2 H t 2 H α 2 + λ γ 2 ( T t ) + 2 2 2 H 1 β 2 T 2 H t 2 H ; d 2 = d 1 α 2 + λ γ 2 ( T t ) + 2 2 2 H 1 β 2 T 2 H t 2 H ; d 3 = ln L 2 K S t + r q + α 2 2 + λ γ 2 2 ( T t ) + 1 2 2 H 2 β 2 T 2 H t 2 H α 2 + λ γ 2 ( T t ) + 2 2 2 H 1 β 2 T 2 H t 2 H ; d 4 = d 3 α 2 + λ γ 2 ( T t ) + 2 2 2 H 1 β 2 T 2 H t 2 H ; κ ( t ) = 1 2 ( r q ) ( T t ) α 2 + λ γ 2 ( T t ) + 2 2 2 H 1 β 2 T 2 H t 2 H .
Proof. 
Let V t ( t , S t ) = C d o ( t , S t ) = C d o . Then according to Theorem 3, the value of the down-and-out call option C d o t , S t is given by
C d o t + ( r q ) S t C d o S + α 2 2 + λ γ 2 2 + 2 2 2 H 1 H β 2 t 2 H 1 S t 2 2 C d o S 2 r C d o = 0 ,
with the initial condition C d o T , S T = S T K + , L < S t < + , and the boundary condition C d o ( t , L ) = 0 , 0 t T .
Let
x = ln S t L , C d o t , S t   =   L C ^ ( t , x ) .
Then,
C d o t = L C ^ t , C d o S = L C ^ x x S = L S t C ^ x a n d 2 C d o S 2 = L S t 2 2 C ^ x 2 C ^ x .
Therefore, we can deduce that
C ^ t + ( r q ) C ^ x + α 2 2 + λ γ 2 2 + 2 2 2 H 1 H β 2 t 2 H 1 2 C ^ x 2 C ^ x r C ^ = 0 ,
with the initial condition C ^ T , ln S T L = e x K L + , 0 < x < + , and the boundary condition C ^ ( 0 , t ) = 0 , 0 t T .
Next, let
ω ( τ , η ) = C ^ ( x , t ) e b ( t ) , τ = c ( t ) , η = x + a ( t ) ,
where a ( t ) , b ( t ) and c ( t ) are undetermined functions about t, which are first order differentiable. Then, we derive
C ^ t = e b ( t ) a ( t ) ω η + c ( t ) ω τ b ( t ) ω , C ^ x = e b ( t ) ω η , 2 C ^ x 2 = e b ( t ) 2 ω η 2
and
c ( t ) ω τ + σ ( t ) 2 ω η 2 + r q + a ( t ) σ ( t ) C ^ x r + b ( t ) ω = 0 ,
where σ ( t ) = α 2 2 + λ γ 2 2 + ( 2 2 2 H 1 ) H β 2 t 2 H 1 .
In order for the solution, let
c ( t ) + σ ( t ) = 0 , r q + a ( t ) σ ( t ) = 0 , r + b ( t ) = 0 , a ( T ) = b ( T ) = c ( T ) = 0 ,
to convert (13) into the heat equation.
From (14), a ( t ) , b ( t ) and c ( t ) are given
a ( t ) = t T r q σ ( s ) d s = r q α 2 2 λ γ 2 2 ( T t ) 1 2 2 H 2 β 2 T 2 H t 2 H , b ( t ) = t T r d s = r ( T t ) , c ( t ) = α 2 2 + λ γ 2 2 ( T t ) + 1 2 2 H 2 β 2 T 2 H t 2 H .
Substitute (15) into (13), then the value of the down-and-out call option C d o t , S t is given by
ω τ = 2 ω η 2 ,
with the initial condition ω ( 0 , η ) = ( e η K ) + , 0 < η < + , and the boundary condition ω ( τ , a ( t ) ) = 0 , 0 t T .
Notice that if we just consider PDE with its initial condition, the solution can be obtained by the Poisson formula
ω ( τ , η ) = 1 2 π τ + φ ( z ) e ( η z ) 2 4 τ d z .
To deal with the boundary condition, let G ( z ) = φ ( z ) e [ a ( t ) z ] 2 4 τ when z > 0 . Then we extend G ( z ) to become an odd function in the whole real number field
G ( z ) = φ ( z ) e [ a ( t ) z ] 2 4 τ , z > 0 , φ ( z ) e [ a ( t ) + z ] 2 4 τ , z 0 .
Comparing the above equation and the original initial condition in (16), we obtain extended initial condition which contains the boundary condition
φ ( z ) = e z K L + , z > 0 , e z K L + e a ( t ) z τ , z 0 .
Therefore, (16) will be transformed into a Cauchy problem
ω τ = 2 ω η 2 ,
with the initial condition ω ( 0 , η ) = φ ( η ) , < η < + .
According to (17),
ω ( τ , η ) = 1 2 π τ + φ ( z ) e ( η z ) 2 4 τ d z = 1 2 π τ l n K L + e z K L e ( η z ) 2 4 τ d z 1 2 π τ l n K L e z K L e ( η z ) 2 + 4 a ( t ) z 4 τ d z = 1 2 π τ l n K L + e z ( η z ) 2 4 τ d z 1 2 π τ K L ln K L + e ( η z ) 2 4 τ d z 1 2 π τ l n K L + e z ( η + z ) 2 4 a ( t ) z 4 τ d z + 1 2 π τ K L l n K L + e ( η + z ) 2 4 a ( t ) z 4 τ d z = I 1 + I 2 + I 3 + I 4 .
For convenience, define N ( x ) as follows: N ( x ) = 1 2 π x e t 2 2 d t , which represents the cumulative probability of standard normal distribution.
For I 1 ,
I 1 = 1 2 π τ l n K L + e z ( η z ) 2 4 τ d z = e τ + η 1 2 π τ l n K L + e ( z η 2 τ ) 2 4 τ d z .
Let t = z η 2 τ 2 τ , then
I 1 = e τ + η 1 2 π l n K L η 2 τ 2 τ + e t 2 2 d t = e τ + η N ( d 1 ) ,
where d 1 = η + 2 τ ln K L 2 τ .
Similarly, let t = z η 2 τ and obtain
I 2 = 1 2 π τ K L l n K L + e ( η z ) 2 4 τ d z = K L 1 2 π l n K L η 2 τ + e t 2 2 d t = K L N d 2 ,
where d 2 = η ln K L 2 τ = d 1 2 τ .
For I 3 ,
I 3 = 1 2 π τ l n K L + e z ( η + z ) 2 4 a ( t ) z 4 τ d z = e [ τ + a ( t ) ] [ τ + a ( t ) η ] τ 1 2 π τ l n K L + e [ z + η 2 a ( t ) 2 τ ] 2 4 τ d z .
Let t = z + η 2 a ( t ) 2 τ 2 τ , then
I 3 = e [ τ + a ( t ) ] [ τ + a ( t ) η ] τ 1 2 π l n K L + η 2 a ( t ) 2 τ 2 τ + e t 2 2 d t = e [ τ + a ( t ) ] [ τ + a ( t ) η ] τ N d 3 ,
where d 3 = 2 a ( t ) + 2 τ η ln K L 2 τ .
For I 4 , let t = z + η 2 a ( t ) 2 τ , we derive
I 4 = 1 2 π τ K L l n K L + e ( η + z ) 2 4 a ( t ) z 4 τ d z = K L e a ( t ) [ a ( t ) η ] τ 1 2 π τ l n K L + e [ z + η 2 a ( t ) ] 2 4 τ d z = e a ( t ) [ a ( t ) η ] τ 1 2 π l n K L + η 2 a ( t ) 2 τ + K L e t 2 2 d t = K L e a ( t ) [ a ( t ) η ] τ N d 4 ,
where, d 4 = 2 a ( t ) η ln K L 2 τ = d 3 2 τ .
Substitute (11) and (12) in them, then
I 1 = S t L e ( r q ) ( T t ) N d 1 ; I 2 = K L N d 2 ; I 3 = e ( r q ) ( T t ) τ l n S t L τ N d 3 = e ( r q ) ( T t ) + 1 ( r q ) ( T t ) τ l n S t L l n S t L N d 3 = e ( r q ) ( T t ) S t L 1 ( r q ) ( T t ) τ L S t N d 3 ; I 4 = K L e [ τ ( r q ) ( T t ) ] l n S t L τ N d 4 = K L S t L 1 ( r q ) ( T t ) τ N d 4 .
By combining I 1 , I 2 , I 3 , a n d I 4 , we have
C d o t , S t = L C ^ ( t , x ) = L e r ( T t ) ω ( τ , η ) = L e r ( T t ) ( I 1 + I 2 + I 3 + I 4 ) = S t e q ( T t ) N d 1 K e r ( T t ) N d 2 S t L κ ( t ) L 2 S t e q ( T t ) N d 3 K e r ( T t ) N d 4 ,
where κ ( t ) = 1 2 ( r q ) ( T t ) α 2 + λ γ 2 ( T t ) + 2 2 2 H 1 β 2 T 2 H t 2 H . □
Corollary 1. 
Suppose that the underlying asset price S t satisfies (1), then at time t the value of the vanilla call option C v a n i l l a ( t , S t ) with the fixed strike price K and the maturity time T is given
C v a n i l l a ( t , S t ) = S t e q ( T t ) N d 1 K e r ( T t ) N d 2 ,
where N ( x ) , d 1 a n d d 2 are shown in Theorem 4.
Proof. 
The proof process is similar to that of Theorem 4.
Let
  • x ˜ = ln S t L , C v a n i l l a t , S t = L C ˜ ( t , x ˜ ) ;
  • ω ˜ ( τ ˜ , η ˜ ) = C ˜ ( t , x ˜ ) e b ( t ) ; τ ˜ = c ( t ) ; η ˜ = x ˜ + a ( t ) ,
where a ( t ) , b ( t ) and c ( t ) are given in (15).
Then, the value of vanilla call option C v a n i l l a ( t , S t ) can be obtained by solving the following Cauchy problem
ω ˜ τ ˜ = 2 ω ˜ η ˜ 2 ,
with the initial condition ω ˜ ( 0 , η ˜ ) = ( e η ˜ K ) + , 0 < η ˜ < + .
The remaining calculation process can be obtained by referring to the solution process of (18). □
Corollary 2. 
Suppose that the underlying asset price S t satisfies (1), then at time t the value of the vanilla put option P v a n i l l a ( t , S t ) at time t with the fixed strike price K and the maturity time T is
P v a n i l l a ( t , S t ) = K e r ( T t ) N d 2 S t e q ( T t ) N d 1 ,
where N ( x ) , d 1 , d 2 are given in Theorem 4.
Proof. 
We just need change the condition to ( K S T ) + and the rest of prove process are similar to Corollary 1. □
Theorem 5. 
Suppose that the underlying asset price S t satisfies (1), then at time t there is the following parity formula between the value of the down-and-out call option C d o ( t , S t ) and the value of the down-and-out put option P d o ( t , S t ) , if the options have the same fixed strike price K, the same fixed barrier L and the same maturity time T:
C d o t , S t + K e r ( T t ) N d 6 S t L κ ( t ) N d 8 = P d o t , S t + S t e q ( T t ) N d 5 S t L κ ( t ) 2 N d 7 ,
where N ( x ) = 1 2 π x e t 2 2 d t , which denotes the cumulative probability of standard normal distribution;
d 5 = ln S t L + r q + α 2 2 + λ γ 2 2 ( T t ) + 1 2 2 H 2 β 2 T 2 H t 2 H α 2 + λ γ 2 ( T t ) + 2 2 2 H 1 β 2 T 2 H t 2 H ; d 6 = d 5 α 2 + λ γ 2 ( T t ) + 2 2 2 H 1 β 2 T 2 H t 2 H ; d 7 = ln L S t + r q + α 2 2 + λ γ 2 2 ( T t ) + 1 2 2 H 2 β 2 T 2 H t 2 H α 2 + λ γ 2 ( T t ) + 2 2 2 H 1 β 2 T 2 H t 2 H ; d 8 = d 7 α 2 + λ γ 2 ( T t ) + 2 2 2 H 1 β 2 T 2 H t 2 H ; κ ( t ) = 1 2 ( r q ) ( T t ) α 2 + λ γ 2 ( T t ) + 2 2 2 H 1 β 2 T 2 H t 2 H .
Proof. 
Let
W d o ( t , S t ) = C d o ( t , S t ) P d o ( t , S t ) ,
which is the difference between the value of the down-and-out call option C d o ( t , S t ) and the down-and-out put option P d o ( t , S t ) at the moment of t. Notice that W d o ( t , S t ) satisfies the following PDE
W d o t + ( r q ) S t W d o S + α 2 2 + λ γ 2 2 + 2 2 2 H 1 H β 2 t 2 H 1 S t 2 2 W d o S 2 r W d o = 0 ,
with the initial condition W d o T , S T = S T K , L < S t < + , and the boundary condition W d o ( t , L ) = 0 , 0 t T .
By analogy with the solution procedure of (16), it can be obtained
W d o t , S t = S t e q ( T t ) N d 5 K e r ( T t ) N d 6 S t L κ ( t ) L 2 S t e q ( T t ) N d 7 K e r ( T t ) N d 8 .
Combining the above equation and (19), Theorem 5 is proved. □
Theorem 6. 
Suppose that the underlying asset price S t satisfies (1), then at time t the value of the down-and-out put option P d o t , S t with the fixed strike price K, the fixed barrier L and the maturity time T is
P d o t , S t = S t e q ( T t ) N d 1 N d 5 K e r ( T t ) N d 2 N d 6 S t L κ ( t ) L 2 S t e q ( T t ) N d 3 N d 7 K e r ( T t ) N d 4 N d 8 ,
where N ( x ) , d 1 d 8 and κ ( t ) are given in Theorems 4 and 5.
Proof. 
By combining Theorems 4 and 5, Theorem 6 is easily proved. □
Theorem 7. 
Suppose that the underlying asset price S t follows Equation (1), the maturity date is T, the fixed strike price is K and the fixed barrier is B, and then the value of the down-and-in call option C d i ( t , S t ) and the value of the down-and-in put option P d i ( t , S t ) at time t are given, respectively:
C d i t , S t = S t L κ ( t ) L 2 S t e q ( T t ) N d 3 K e r ( T t ) N d 4 , P d i t , S t = K e r ( T t ) N d 6 S t e q ( T t ) N d 5 + S t L K ( t ) L 2 S t e q ( T t ) N d 3 N d 7 K e r ( T t ) N d 4 N d 8 ,
where N ( x ) ,   d 3 d 8 and κ ( t ) are detailed in Theorems 4 and 5.
Proof. 
When other conditions are the same, a portfolio with a out option and the corresponding in option will always be able to exercise one of their option right, which is equivalent to a vanilla option
V v a n i l l a ( t , S t ) = V d o ( t , S t ) + V d i ( t , S t ) = V u o ( t , S t ) + V u i ( t , S t ) ,
where V v a n i l l a ( t , S t ) is the European option, V d o ( t , S t ) , V d i ( t , S t ) , V u o ( t , S t ) and V u i ( t , S t ) , respectively, denote the corresponding value of the down-and-out option, the down-and-in option, the up-and-out option and the up-and-in option.
Therefore, C d i ( t , S t ) = C v a n i l l a ( t , S t ) C d o ( t , S t ) ; P d i ( t , S t ) = P v a n i l l a ( t , S t ) P d o ( t , S t ) .
Using Corollary 1, Corollary 2, Theorems 4 and 6, Theorem 7 is proved. □
Above all, the pricing formulas of all four types of downward barrier options have been given. Similarly, the pricing formulas corresponding to four types of upward barrier options can be deduced.

5. Numerical Experiment

In this section, numerical experiments are conducted to discuss the effects of the barrier price L, the Hurst index H, the jump intensity λ and volatility α , β , γ on barrier options by MATLAB and R language software. In this section, we just take the down-and-out call option as an example for space constraints.
Firstly, parameters are assumed as follows:
t = 0 , T = 0.5 , K = 100 , H = 0.75 , α = β = γ = 0.4 , λ = 1 .
According to Theorem 4, the value of down-and-out call option C d o ( t , S t ) under different barrier prices and stock prices can be obtained, which are given in Table 1. At the same time, Figure 1 is drawn to describe the trend of option value affected by barrier price under different stock prices.
Observing Table 1 and Figure 1, it can be seen that when the stock price is fixed, the value of down-and-out call option decreases with the growth of barrier price. When other conditions remain unchanged, with the rising barrier price, the possibility of down-and-out call option termination is increasing, so the option value will continue to decline. In particular, when the barrier price increases to the initial stock price, the option will be knocked out at once, which means it has no value any more.
Then, in order to discuss the impact of the Hurst index H and the jump intensity λ on the option price, a new hypothesis is proposed as follows:
t = 0 , T = 0.5 , K = 100 , L = 70 , α = β = γ = 0.4 .
Take the different H, λ , and other assumptions remain unchanged to obtain the option value under various conditions, as shown in Table 2.
Figure 2 is the variation diagram of the value of down-and-out call option with the different Hurst index and jump intensity, when S 0 is fixed at 100. The relationships between the value of down-and-out option and the Hurst index is positive. The larger the Hurst index is, the more stable the underlying asset price is. This means the price fluctuation will be smaller, which denotes the corresponding option value will be smaller.
At the same time, the value of down-and-out option and the jump intensity change in the same direction. The jump intensity represents the unsystematic risk. When it increases, the underlying asset will has more intense fluctuations, which means higher upper limit and invariant lower bound. Therefore, the option value will rise.
Finally, for rigorousness, we verify the positive correlation between the volatility and the value of down-and-out call option, where α , β , a n d γ are different. Assume that the parameter selection is as follows:
t = 0 , T = 0.5 , K = 100 , L = 70 , H = 0.75 , λ = 1 .
Let σ ^ = ( α , β , γ ) , and make
σ ^ 1 = ( 0.1 , 0.15 , 0.2 ) ; σ ^ 2 = ( 0.2 , 0.25 , 0.3 ) ; σ ^ 3 = ( 0.3 , 0.35 , 0.4 ) ; σ ^ 4 = ( 0.4 , 0.45 , 0.5 ) .
According to Theorem 4, the value of down-and-out call option under different volatility can be obtained, which are shown in Table 3. The value of down-and-out call option increases with the rise of the volatility, which is consistent with the fact.

6. Conclusions

This paper investigated the barrier option pricing model in the environment of the sub-mixed fractional Brownian motion with jump intensity. Through the self-financing strategy, we derive the B-S type PDE of the derivatives. Then, the value of the down-and-out option is obtained by applying transformation techniques. Meanwhile, the parity formula between barrier call option and barrier put option can be given by a similar method. Next, using the linear relationship between the knock-out option and the knock-in option, the value of the knock-in option can be deduced. In Section 5, the numerical experiment is carried out where we take the down-and-out call option as an example. According to the results shown in Figure 1, Figure 2, and Table 3, the following relationships can be found: The barrier price and Hurst index are inversely related to the value of the down-and-out call option, while the jump intensity and volatility are positively correlated with it. In the above numerical experiment, the parameter values are averages. It may affect the application of the model according to reference [32], and we will try our best to overcome this limitation in future research. Meanwhile, the Asian barrier options can also be considered to extend the model used in this paper.

Author Contributions

Writing—original draft preparation, B.J.; writing—review and editing, X.T. and Y.J. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (No. 11771399).

Data Availability Statement

Not applicable.

Acknowledgments

We are very grateful to the reviewers for their valuable comments.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dassios, A.; Lim, J.W. Recursive formula for the double-barrier Parisian stopping time. J. Appl. Probab. 2018, 55, 282–301. [Google Scholar] [CrossRef]
  2. Funahashi, H.; Higuchi, T. An analytical approximation for single barrier options under stochastic volatility models. Ann. Oper. Res. 2018, 266, 129–157. [Google Scholar] [CrossRef]
  3. Guillaume, T. Closed form valuation of barrier options with stochastic barriers. Ann. Oper. Res. 2021, 1–30. [Google Scholar] [CrossRef]
  4. Gao, Y.; Jia, L. Pricing formulas of barrier-lookback option in uncertain financial markets. Chaos Solitons Fractals 2021, 147, 110986. [Google Scholar] [CrossRef]
  5. Shreve, S.E. Stochastic Calculus for Finance II: Continuous-Time Models; Springer: New York, NY, USA, 2004. [Google Scholar]
  6. Merton, R.C. Theory of rational option pricing. Bell Econ. Manag. Sci. 1973, 4, 141–183. [Google Scholar] [CrossRef] [Green Version]
  7. Rubinstein, M. Breaking down the barriers. Risk 1991, 4, 28–35. [Google Scholar]
  8. Black, F.; Scholes, M. The Pricing of Options and Corporate Liabilities. J. Political Econ. 1973, 81, 637–654. [Google Scholar] [CrossRef] [Green Version]
  9. Ding, Z.; Granger, C.W.; Engle, R.F. A long memory property of stock market returns and a new model. J. Empir. Financ. 1993, 1, 83–106. [Google Scholar] [CrossRef]
  10. Shiryaev, A.N. Essentials of Stochastic Finance: Facts, Models, Theory; World Scientific: Singapore, 1999. [Google Scholar]
  11. Necula, C. Option pricing in a fractional Brownian motion environment. Adv. Econ. Financ. Res.-Dofin Work. Pap. Ser. 2008, 2, 259–273. [Google Scholar] [CrossRef]
  12. Kolmogorov, A.N. Wienersche spiralen und einige andere interessante kurven in hilbertscen raum, cr (doklady). Acad. Sci. URSS (NS) 1940, 26, 115–118. [Google Scholar]
  13. Chen, Q.; Zhang, Q.; Liu, C. The pricing and numerical analysis of lookback options for mixed fractional Brownian motion. Chaos Solitons Fractals 2019, 128, 123–128. [Google Scholar] [CrossRef]
  14. Bian, L.; Li, Z. Fuzzy simulation of European option pricing using sub-fractional Brownian motion. Chaos Solitons Fractals 2021, 153, 111442. [Google Scholar] [CrossRef]
  15. Wang, J.; Yan, Y.; Chen, W.; Shao, W.; Tang, W. Equity-linked securities option pricing by fractional Brownian motion. Chaos Solitons Fractals 2021, 144, 110716. [Google Scholar] [CrossRef]
  16. Cheridito, P. Arbitrage in fractional Brownian motion models. Financ. Stochastics 2003, 7, 533–553. [Google Scholar] [CrossRef]
  17. Bender, C.; Elliott, R.J. Arbitrage in a discrete version of the Wick-fractional Black-Scholes market. Math. Oper. Res. 2004, 29, 935–945. [Google Scholar] [CrossRef]
  18. Björk, T.; Hult, H. A note on Wick products and the fractional Black-Scholes model. Financ. Stochastics 2005, 9, 197–209. [Google Scholar] [CrossRef]
  19. Bojdecki, T.; Gorostiza, L.G.; Talarczyk, A. Sub-fractional Brownian motion and its relation to occupation times. Stat. Probab. Lett. 2004, 69, 405–419. [Google Scholar] [CrossRef]
  20. Charles, E.N.; Mounir, Z. On the sub-mixed fractional Brownian motion. Appl.-Math.-J. Chin. Univ. 2015, 30, 27–43. [Google Scholar] [CrossRef] [Green Version]
  21. Tudor, C. Some properties of the sub-fractional Brownian motion. Stochastics Int. J. Probab. Stoch. Process. 2007, 79, 431–448. [Google Scholar] [CrossRef]
  22. Xu, F.; Zhou, S. Pricing of perpetual American put option with sub-mixed fractional Brownian motion. Fract. Calc. Appl. Anal. 2019, 22, 1145–1154. [Google Scholar] [CrossRef]
  23. Merton, R.C. Option pricing when underlying stock returns are discontinuous. J. Financ. Econ. 1976, 3, 125–144. [Google Scholar] [CrossRef] [Green Version]
  24. Zhou, Q.; Yang, J.J.; Wu, W.X. Pricing vulnerable options with correlated credit risk under jump-diffusion processes when corporate liabilities are random. Acta Math. Appl. Sin. Engl. Ser. 2019, 35, 305–318. [Google Scholar] [CrossRef]
  25. Sun, W.; Zhao, Y.; MacLean, L. Real Options in a Duopoly with Jump Diffusion Prices. Asia-Pac. J. Oper. Res. 2021, 38, 2150009. [Google Scholar] [CrossRef]
  26. Zhang, W.G.; Li, Z.; Liu, Y.J.; Zhang, Y. Pricing European option under fuzzy mixed fractional Brownian motion model with jumps. Comput. Econ. 2021, 58, 483–515. [Google Scholar] [CrossRef]
  27. Liu, M. Two possible types of superfluidity in crystals. Phys. Rev. B 1978, 18, 1165. [Google Scholar] [CrossRef]
  28. Callen, H.B. Thermodynamics and an Introduction to Thermostatistics; John Wiley & Sons: New York, NY, USA, 1985. [Google Scholar]
  29. Appel, D.; Grabinski, M. The origin of financial crisis: A wrong definition of value. Port. J. Quant. Methods 2011, 2, 33. [Google Scholar]
  30. Klinkova, G.; Grabinski, M. Conservation laws derived from systemic approach and symmetry. Int. J. Latest Trends Fin. Ecol. Sci. Vol. 2017, 7, 1307. [Google Scholar]
  31. Tankov, P. Financial Modelling with Jump Processes; Chapman and Hall/CRC: London, UK, 2003. [Google Scholar]
  32. Grabinski, M.; Klinkova, G. Wrong use of averages implies wrong results from many heuristic models. Appl. Math. 2019, 10, 605. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The change curve of the value of down-and-out option for different barrier prices.
Figure 1. The change curve of the value of down-and-out option for different barrier prices.
Fractalfract 06 00244 g001
Figure 2. Plot of down-and-out option value against different Hurst index and jump intensity values.
Figure 2. Plot of down-and-out option value against different Hurst index and jump intensity values.
Fractalfract 06 00244 g002
Table 1. The value of down-and-out options for different barrier prices and stock prices.
Table 1. The value of down-and-out options for different barrier prices and stock prices.
L S 0 = 120 S 0 = 110 S 0 = 100 S 0 = 90 S 0 = 80
1156.350
11012.1120.000
10022.44511.5910.000
9526.96916.5415.534
9031.04220.92610.5510.000
8534.65424.73115.0955.086
8037.80327.95619.1499.5100.000
7540.49530.61322.70213.2644.507
7042.74832.73025.75716.3518.275
6544.58734.35228.32418.79811.321
6046.04935.53330.43020.65313.680
Table 2. The value of down-and-out option for different Hurst index and jump intensity.
Table 2. The value of down-and-out option for different Hurst index and jump intensity.
S 0 H = 0.75 H = 0.85 H = 0.95
λ = 0 λ = 2 λ = 4 λ = 0 λ = 2 λ = 4 λ = 0 λ = 2 λ = 4
12031.11936.62140.01729.76135.81239.50128.47235.04539.018
11527.25732.65235.83925.86031.87935.36024.50031.14134.910
11023.55428.76031.70822.14428.03331.27020.73727.33430.857
10520.02424.94927.62618.63224.27827.23217.21223.62826.859
10016.67921.21723.59015.34420.61323.24413.95320.02522.916
9513.52417.56319.60012.29117.03819.30510.98216.52219.026
9010.55913.98215.6509.47913.54515.4118.31113.11515.183
857.77010.46211.7346.89510.12511.5535.9359.79111.380
805.1296.9887.8424.5106.7597.7223.8216.5317.606
752.5823.5343.9592.2593.4183.9001.8963.3023.843
Table 3. The value of down-and-out option against the volatility of the underlying asset.
Table 3. The value of down-and-out option against the volatility of the underlying asset.
S 0 σ ^ 1 =
(0.1, 0.15, 0.2)
σ ^ 2 =
(0.2, 0.25, 0.3)
σ ^ 3 =
(0.3, 0.35, 0.4)
σ ^ 4 =
(0.4, 0.45, 0.5)
12025.32328.44132.36036.149
11520.91824.46728.50632.202
11016.77020.70224.79028.338
10512.96317.17621.22224.559
1009.58113.91717.80720.867
956.70410.94814.54917.259
904.3858.28111.44313.729
852.6375.9108.47710.267
801.4163.8035.6236.856
750.6081.8862.8383.467
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ji, B.; Tao, X.; Ji, Y. Barrier Option Pricing in the Sub-Mixed Fractional Brownian Motion with Jump Environment. Fractal Fract. 2022, 6, 244. https://doi.org/10.3390/fractalfract6050244

AMA Style

Ji B, Tao X, Ji Y. Barrier Option Pricing in the Sub-Mixed Fractional Brownian Motion with Jump Environment. Fractal and Fractional. 2022; 6(5):244. https://doi.org/10.3390/fractalfract6050244

Chicago/Turabian Style

Ji, Binxin, Xiangxing Tao, and Yanting Ji. 2022. "Barrier Option Pricing in the Sub-Mixed Fractional Brownian Motion with Jump Environment" Fractal and Fractional 6, no. 5: 244. https://doi.org/10.3390/fractalfract6050244

APA Style

Ji, B., Tao, X., & Ji, Y. (2022). Barrier Option Pricing in the Sub-Mixed Fractional Brownian Motion with Jump Environment. Fractal and Fractional, 6(5), 244. https://doi.org/10.3390/fractalfract6050244

Article Metrics

Back to TopTop