Next Article in Journal
Screen-Printed Metamaterial Absorber Using Fractal Metal Mesh for Optical Transparency and Flexibility
Next Article in Special Issue
Application of the Fractal Brownian Motion to the Athens Stock Exchange
Previous Article in Journal
On Numerical Simulations of Variable-Order Fractional Cable Equation Arising in Neuronal Dynamics
Previous Article in Special Issue
The Importance of Non-Systemically Important Banks—A Network-Based Analysis for China’s Banking System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Forward Starting Option Pricing under Double Fractional Stochastic Volatilities and Jumps

1
School of Science, Xi’an University of Posts and Telecommunications, Xi’an 710121, China
2
School of Computer Science and Technology, Xi’an University of Posts and Telecommunications, Xi’an 710121, China
*
Author to whom correspondence should be addressed.
Fractal Fract. 2024, 8(5), 283; https://doi.org/10.3390/fractalfract8050283
Submission received: 6 March 2024 / Revised: 27 April 2024 / Accepted: 4 May 2024 / Published: 8 May 2024

Abstract

:
This paper aims to provide an effective method for pricing forward starting options under the double fractional stochastic volatilities mixed-exponential jump-diffusion model. The value of a forward starting option is expressed in terms of the expectation of the forward characteristic function of log return. To obtain the forward characteristic function, we approximate the pricing model with a semimartingale by introducing two small perturbed parameters. Then, we rewrite the forward characteristic function as a conditional expectation of the proportion characteristic function which is expressed in terms of the solution to a classic PDE. With the affine structure of the approximate model, we obtain the solution to the PDE. Based on the derived forward characteristic function and the Fourier transform technique, we develop a pricing algorithm for forward starting options. For comparison, we also develop a simulation scheme for evaluating forward starting options. The numerical results demonstrate that the proposed pricing algorithm is effective. Exhaustive comparative experiments on eight models show that the effects of fractional Brownian motion, mixed-exponential jump, and the second volatility component on forward starting option prices are significant, and especially, the second fractional volatility is necessary to price accurately forward starting options under the framework of fractional Brownian motion.

1. Introduction

Forward starting options belong to European-style exotic options which start at the determination time of the strike and expirate further in the future. They are the building blocks of cliquets which are popular in the world of equity derivatives because of their retaining of the upside potential with protection against downside risk. However, the valuation of forward starting options is not easy, because these options not only depend on the terminal value of the underlying price, but also on the asset price at the determination time of the strike.
Rubinstein [1] first evaluated forward starting options based on the assumption of geometric Brownian motion. Via the changing of numeraire, Kruse and Nögel [2] considered the counterpart by deriving a semi-analytical pricing formula with two integrals which have to be evaluated under Heston stochastic volatility. With Feynman-Kac theorem, Ahlip and Rutkowski [3] and Lin and He [4] extend the method in [2] to stochastic volatility and stochastic interest rate, and regime-switching stochastic volatility models, respectively. Using a similar approach, Guo and Hung [5] and Ahlip et al. [6] considered forward starting options pricing under stochastic volatility, jumps, and stochastic interest rates, respectively. Lucié [7] and Haastrecht and Pelsser [8] reduced the formula in [2] to one in the form of a single integral under the Heston model and a stochastic volatility and stochastic interest rate model, respectively, and improved the efficiency of pricing using the fast Fourier transform technique (FFT). Given that the accuracy of FFT depends on a dampening factor, Zhang and Geng [9] and Zhang and Sun [10] expressed forward starting options prices by using a Fourier cosine series expansion under a stochastic volatility with jump model and a stochastic volatilities, stochastic interest rate, and double jumps model, respectively. Recently, Hata et al. [11] proposed a Taylor expansion approach for evaluating forward starting options under the Hull-White stochastic volatility model.
All of the above studies are conducted under the models driven by standard Brownian motion which are Markovian or memoryless. However, many studies show that asset price fluctuations exhibit “long memory” [12,13,14,15,16,17,18] or “short memory” [19,20,21,22] which can be captured by stochastic volatility models driven by fractional Brownian motion with the Hurst index H ( 1 / 2 ,   1 ) or H ( 0 ,   1 / 2 ) , respectively. In addition, jumps in the asset price were observed by Coqueret and Tavin [23], Jin and Hong [24], Bates [25], and Wang and Xia [26]. Recently, by introducing two fractional volatilities into a mixed-exponential jump-diffusion model [27], Zhang et al. [28] propose a double fractional Heston mixed-exponential jump-diffusion model and demonstrate that the model fits the market better than the double Heston mixed-exponential jump-diffusion model and the double Heston model driven by standard Brownian motion.
However, because the model is not Markovian, all pricing methods based on the Feynman-Kac theorem are not easily applied, which poses a huge challenge for forward starting options pricing. Zhang et al. [28] proposed a model approximation method by approximating the pricing model with a semimartingale under which the Feynman-Kac theorem can be activated. Motivated by the performance of the model and the pricing method in [28], this paper extends the model and pricing method in [28] to the case of forward starting option. Recently, Wang and Guo [29] evaluated variance and volatility swaps under a similar model. Different from Wang and Guo [29] in which only one fractional volatility was used, we consider two fractional volatilities and combine model approximation and the Fourier transform technique to evaluate forward starting options.
The main goal of this paper is to provide a fast and efficient method for evaluating forward starting options under two fractional volatilities and jumps. The main contributions of this paper are twofold. Firstly, this paper derives the forward characteristic function of log return proposed in [28]. Secondly, this paper proposes a pricing method for forward starting options under double fractional volatilities and jumps. Compared to existing option pricing methods for fractional Brownian motion setting, our method can be used for both cases of H ( 1 / 2 ,   1 ) and H ( 0 ,   1 / 2 ) . The rest of the paper is organized as follows. Section 2 presents the model and pricing problem. Section 3 provides an approximate model and the derivation of the forward characteristic function. Section 4 dilates the pricing method for forward starting options. Section 5 presents some numerical experiments. Section 6 concludes.

2. The Model and Pricing Problem

2.1. The Model

Assume that Ω , F , F t 0 t T , P is a complete probability space with the right continuous filtration { F t } 0 t T and the risk-neutral probability measure P . In the double fractional Heston mixed-exponential jump-diffusion model, the dynamics of the asset price S t are given by the following system:
d S t S t = ( r λ δ ) d t + υ 1 t d W 1 t S + υ 2 t d W 2 t S + d k = 1 N t ( ζ k 1 ) , d υ 1 t = κ 1 θ 1 υ 1 t d t + σ 1 υ 1 t d B 1 t H 1 , d υ 2 t = κ 2 θ 2 υ 2 t d t + σ 2 υ 2 t d B 2 t H 2 ,
where W 1 t S and W 2 t S are both standard Brownian motion, B 1 t H 1 and B 2 t H 2 are both fractional Brownian motion with a Hurst index H 1 and H 2 , respectively, and N t is a Poisson process with rate λ . Let ζ = ( ζ k ) k 1 be a sequence of independent identically distributed nonnegative random variables, such that Y = ln ζ has a mixed-exponential distribution with the density as follows:
f Y ( y ) = p k = 1 m p k η k e η k y I y 0 + q j = 1 n q j θ ^ j e θ ^ j y I y < 0 ,
where η k > 1 , θ ^ j > 0 , p 0 , q = 1 p 0 ,   p k ,   q j ( ,   ) ,   k = 1 m p k = 1 ,   j = 1 n q j = 1 . To make f Y ( y ) a density, a simple sufficient condition is k = 1 M p k η k 0 ( M = 1 , , m ) ,   j = 1 L q j θ ^ j 0   ( L = 1 , , n ) and a necessary condition is p 1 > 0 , q 1 > 0 ,   k = 1 m p k η k 0 ,   j = 1 n q j θ ^ j 0 . Assume that δ = E [ ζ 1 ] . With the density (2), basic calculus gives the following:
δ = p k = 1 m p k η k 1 η k 1 + q j = 1 n q j θ ^ j 1 θ ^ j + 1 1 .
For j = 1 ,   2 , κ j ,   θ j ,   σ j satisfying 2 κ j θ j σ j 2 denote the mean reversion rate, long-run mean, and volatilities of variance υ j t , respectively. According to Thao [30], fractional Brownian motion B j t H j has the following form:
B j t H j = 1 Γ ( H j + 1 / 2 ) 0 t ( t s ) H j 1 / 2 d W j s v ( j = 1 ,   2 ) ,
where Γ ( ) is the gamma function, and W j s v   are both standard Brownian motion. Suppose processes W j t S and W j t v are independent of ζ and N t , respectively, while W j t S and W j t v are correlated by setting C o v d W 1 t S ,   d W 1 t v = ρ 1 d t , C o v d W 2 t S ,   d W 2 t v = ρ 2 d t to mimic the asset-volatility leverage effect. Suppose S 0 = S ,   υ 10 = υ 1 ,   υ 20 = υ 2 . According to Zhang et al. [28], model (1) can be written as the following stochastic partial integro-differential system:
d S t S t = ( r λ δ ) d t + j = 1 2 υ j t d W j t S + d k = 1 N t ( ζ k 1 ) , d υ 1 t = 1 Γ ( H + 1 / 2 ) 0 t ( t s ) H 1 / 2 κ 1 θ 1 υ 1 s d s + 1 Γ ( H + 1 / 2 ) 0 t ( t s ) H 1 / 2 σ 1 υ 1 s d W 1 s v , d υ 2 t = 1 Γ ( H + 1 / 2 ) 0 t ( t s ) H 1 / 2 κ 2 θ 2 υ 2 s d s + 1 Γ ( H + 1 / 2 ) 0 t ( t s ) H 1 / 2 σ 2 υ 2 s d W 2 s v .

2.2. The Pricing Problem

Under the risk-neutral measure P , the value of a forward starting put option can be expressed on the forward log return as follows:
V F ( t 0 , T ) = e r T E p e ln K e ln S T S t 0 + F 0
with determination time t 0 , maturity T , constant interest rate r , and strike price K . According to [28], if the characteristic function of ln S T S t 0 is given, one can evaluate V F ( t 0 , T ) effectively by using the Fourier transform technology.
Let φ F ( u , τ ) = E p exp i u ln S T S t 0 F t 0   ( τ = T t 0 ) be the characteristic function of ln S T S t 0 , which we call the forward characteristic function. Because B j t H j are not semimartingales, φ F ( u , τ ) is not a semimartingale. Therefore, we cannot obtain φ F ( u , τ ) by solving a PDE like in the standard Brownian motion environment. However, according to El Euch and Gatheral [31], the forward characteristic function φ F ( u , τ ) has a similar structure to that of the classic Heston model with jumps as follows:
φ F ( u , τ ) = E exp i u λ τ δ + i u j = 1 N τ Y j exp j = 1 2 h j ( u , τ ) + g j ( u , τ ) υ j ,
where h j ( u , τ ) = κ j θ j 0 τ g j ( u , t ) d t ,   g j ( u , τ ) = 1 Γ ( H j + 1 / 2 ) 0 τ ( τ t ) H j 1 / 2 L j ( u , t ) d t ,   L j ( u , τ ) are the solutions to the following fractional Riccati equations:
D H 1 + 1 / 2 L 1 ( u , τ ) = 1 / 2 ( u 2 i u ) + ( i u ρ 1 σ 1 κ 1 ) L 1 ( u , τ ) + 1 / 2 σ 1 2 L 1 2 ( u , τ ) , D H 2 + 1 / 2 L 2 ( u , τ ) = 1 / 2 ( u 2 i u ) + ( i u ρ 2 σ 2 κ 2 ) L 2 ( u , τ ) + 1 / 2 σ 2 2 L 2 2 ( u , τ ) , I 1 / 2 H 1 L 1 ( u , 0 ) = I 1 / 2 H 2 L 2 ( u , 0 ) = 0 ,
where
D H j + 1 / 2 L j ( u , τ ) = 1 Γ ( 1 / 2 H j ) d d τ 0 τ ( τ t ) H j 1 / 2 L j ( u , t ) d t ,
I 1 / 2 H j L j ( u , τ ) = 1 Γ ( 1 / 2 H j ) 0 τ ( τ t ) H j 1 / 2 L j ( u , t ) d t .

3. The Forward Characteristic Function

3.1. Model Approximation

By introducing two small perturbed parameters ε 1 ,   ε 2   ( 0 < ε 1 ,   ε 2 1 ) , the processes B j t H j   ( j = 1 ,   2 ) in Formula (3) can be approximated by the following semimartingales:
B j t ε j , H j = 0 t ( t s + ε j ) H j 1 / 2 d W j s v
in L 2 ( Ω ) as ε 1 0 , ε 2 0 . Differentiating Formula (8), one has the following:
d B j t ε j , H j = ( H j 1 / 2 ) ψ j t d t + ε j H j 1 / 2 d W j t v ,
where ψ j t = 0 t ( t s + ε j ) H j 3 / 2 d W j s v are semimartingales. Putting Formula (9) into model (1), one can approximate model (1) by applying the following classic stochastic partial integro-differential system:
d S t ε 1 , ε 2 S t ε 1 , ε 2 = ( r λ δ ) d t + + υ 1 t d W 1 t S + υ 2 t d W 2 t S + d k = 1 N t ( ζ k 1 ) , d υ 1 t ε 1 = ( H 1 1 / 2 ) ψ 1 t σ 1 υ 1 t ε 1 + κ 1 θ 1 υ 1 t ε 1 d t + ε 1 H 1 1 / 2 σ 1 υ 1 t ε 1 d W 1 t v , d υ 2 t ε 2 = ( H 2 1 / 2 ) ψ 2 t σ 2 υ 2 t ε 2 + κ 2 θ 2 υ 2 t ε 2 d t + ε 2 H 2 1 / 2 σ 2 υ 2 t ε 2 d W 2 t v .
Because ψ j t = 0 t ( t s + ε j ) H j 3 / 2 d W j s v are semimartingales, the forward characteristic function φ F ε 1 , ε 2 ( u , τ ) = E p exp i u ln S T ε 1 , ε 2 S t 0 ε 1 , ε 2 F t 0 is a semimartingale. Based on the tower law of conditional expectations, one can write φ F ε 1 , ε 2 ( u , τ ) as follows:
φ F ε 1 , ε 2 ( u , τ ) = E P E P ( exp i u ln S T ε 1 , ε 2 K i u ln S t 0 ε 1 , ε 2 K | F t 0 ) | F 0 = E P φ ( u , τ ) exp i u ln S t 0 ε 1 , ε 2 K | F 0 ,
where φ ( u , τ ) is the characteristic function of ln ( S T ε 1 , ε 2 / K ) conditional on the filtration F t 0 . If we consider exp i u ln S t 0 ε 1 , ε 2 K as the proportion of φ ( u , τ ) , we only need to calculate the expectation of the proportion characteristic function to obtain the forward characteristic function.

3.2. Calculation of the Expectation of the Proportion Characteristic Function

Lemma 1. 
Under the risk-neutral measure P , the characteristic function of the log asset price ln ( S T ε 1 , ε 2 / K ) conditional on the filtration F t 0 , φ ( u , τ ) is given by the following:
φ ( u , τ ) = exp i u ln S t 0 ε 1 , ε 2 K + M ( u , τ ) + j = 1 2 N j ( u , τ ) υ j t 0 ,
where
M ( u , τ ) = i u r τ + λ τ Λ + j = 1 2 κ j θ j Δ j 2 ( D j + γ j ) τ 2 ln 1 G j e D j τ 1 G j , N j ( u , τ ) = j = 1 2 1 Δ j 2 ( D j + γ j ) 1 e D j τ 1 G j e D j τ , τ = T t 0 , D j = γ j 2 4 a j b j ,   γ j = i u ρ j Δ j κ j ,   G j = γ j + D j γ j - D j , a j = 1 2 Δ j 2 ,   b j = 1 2 u ( i + u ) , Δ j = ε j H j 1 / 2 σ j , Λ = p i = 1 m p i η i 1 η i i u + q j = 1 n q j θ ^ j 1 θ ^ j + i u 1 i u δ .
Proof. 
See Zhang et al. [28]. □
Lemma 2. 
Under the risk-neutral measure P ,
E P j = 1 2 N j ( u , τ ) υ j t 0 ( t ) F 0 = exp M ¯ ( u , τ ) + j = 1 2 N ˜ j ( u , t ) υ j ,
where for j = 1 ,   2 ,
M ¯ ( u , τ ) = j = 1 2 2 κ j θ j Δ j 2 ln 1 N j ( u , τ ) l t j , N ˜ j ( u , τ ) = N j ( u , τ ) e κ t 0 1 N j ( u , τ ) l t j , l t j = 2 κ j Δ j 2 ( 1 e κ j t 0 ) .
Proof. 
Set Π = E P j = 1 2 N j ( u , τ ) υ j t 0 ( t ) F 0 . Π As a semimartingale; the Feynman-Kac PDE for Π is followed as follows:
Π t = j = 1 2 κ j ( θ j υ j t ε j ) + H j 1 2 ψ j t σ j υ j t ε j Π υ j t ε j + 1 2 Δ j 2 υ j t ε j 2 Π υ j t ε j 2 , Π | t = 0 = exp j = 1 2 N j ( u , τ ) υ j .
Since ψ j t   ( j = 1 , 2 ) are semimartingales, one has ψ j 0 = E ( ψ j t ) = 0 . In the light of the affine structure of model (10), we conjuncture that Π has the following form:
Π = exp M ¯ ( u , t ) + j = 1 2 N ˜ j ( u , t ) υ j .
Putting Formula (14) into Formula (13), one has the following ODEs:
M ¯ t = κ 1 θ 1 N ˜ 1 + k 2 θ 2 N ˜ 2 , N ˜ 1 t = κ 1 N ˜ 1 + 1 / 2 Δ 1 2 N ˜ 1 2 , N ˜ 2 t = κ 2 N ˜ 2 + 1 / 2 Δ 2 2 N ˜ 2 2
with M ¯ ( u , 0 ) = 0 ,   N ˜ 1 ( u , 0 ) = N 1 ( u , τ ) ,   N ˜ 2 ( u , 0 ) = N 2 ( u , τ ) . Solving the ODEs (15), Lemma 2 follows. □
Combining Lemma 1 with Lemma 2, one can obtain the expectation of the proportion characteristic function; thus, one has the forward characteristic function φ F ε 1 , ε 2 ( u , τ ) .
Theorem 1. 
Assuming that the dynamics of asset price S t ε 1 , ε 2 are given by model (10), then the forward characteristic function of ln S T ε 1 , ε 2 S t 0 ε 1 , ε 2 ,   φ F ε 1 , ε 2 ( u , τ ) is given by the following:
φ F ε 1 , ε 2 ( u , τ ) = exp M ˜ ( u , τ ) + j = 1 2 N ˜ j ( u , τ ) υ j t 0 ,
where
M ˜ ( u , τ ) = M ( u , τ ) j = 1 2 2 κ j θ j Δ j 2 ln 1 N j ( u , τ ) / l t j , N ˜ j ( u , τ ) = N j ( u , τ ) e κ t 0 1 N j ( u , τ ) l t j , l t j = 2 κ j Δ j 2 ( 1 e κ j t 0 ) .
Based on Theorem 1, differential calculation gives the following outcome.
Theorem 2. 
Assuming that the dynamics of asset price S t ε 1 , ε 2 are given by model (10), and c n = n ln φ F ε 1 , ε 2 ( u , τ ) u n u = 0 denotes the n-th cumulant of ln S T ε 1 , ε 2 S t 0 ε 1 , ε 2 , one has the following:
c 1 = j = 1 2 θ j ( 1 e κ j τ κ j τ ) 2 κ j + e κ j τ 1 2 k j 2 κ j θ j Δ j 2 h t j + e κ j t υ j    + λ τ p k = 1 m p k η k - q l = 1 n q l θ ^ l - p k = 1 m p k η k η k - 1 + q l = 1 n q l θ ^ l θ ^ l + 1 - 1 , c 2 = j = 1 2 κ j θ j γ j τ Δ j 2 5 θ j Δ j 2 8 κ j 3 + 2 θ j ρ j Δ j κ j 2 θ j κ j    + e - κ j τ θ j Δ j 2 + θ j κ j τ Δ j 2 2 κ j 3 + θ j - θ j ρ j τ Δ j κ j - 2 θ j ρ j Δ j κ j 2 + θ j Δ j 2 8 κ j 3 e 2 κ j τ 2 κ j θ j Δ j 2 B j h t j + ( B j ) 2 h t j 2 e κ j t B j h t j + 2 B j h t j υ j + λ τ p k = 1 m 2 p k η k 2 + q l = 1 n 2 q l θ ^ l 2 , c 4 = 0 ,
where
γ j = Δ j 2 κ j ρ j Δ j 3 κ j 2 + Δ j 4 4 κ j 3 , B j = i 1 e κ j τ 2 κ j , B j = Δ j 2 4 κ j 3 e 2 κ j τ + e κ j τ Δ j 2 τ 2 κ j 2 + 1 Δ j ρ j τ κ j Δ j ρ j κ j 2 1 κ j Δ j 2 4 κ j 3 + Δ j ρ j κ j 2 .

4. The Pricing Algorithm

Let ξ = ln S T ε 1 , ε 2 S t 0 ε 1 , ε 2 ; one can express (5) as follows:
V F ( t 0 , T ) = e r T K ( e ξ 1 ) + f F ε 1 , ε 2 ( ξ | x ) d ξ ,
where f F ε 1 , ε 2 ( ξ | x ) is the density of ξ conditional on x = ln ( S / K ) .
Since f F ε 1 , ε 2 ( ξ | x ) decays to zero as ξ ± , f F ε 1 , ε 2 ( ξ | x ) can be approximated by applying a truncated Fourier cosine series expansion (COS) on the domain [ a , b ] as follows:
f F ε 1 , ε 2 ( ξ | x ) 2 b a k = 0 N 1 φ F ε 1 , ε 2 k π b a , τ e i k π b a ( ln K + a ) cos k π b a ( ξ a ) ,
where [ . ] denotes taking the real part of the argument, and N denotes the number of terms in series expansion. Inserting (18) into (17) and interchanging summation and integration, one can approximate forward starting option value by the following:
V F ( t 0 , T ) 2 b a K e r T k = 0 N 1 φ F ε 1 , ε 2 k π b a , τ e i k π a b a U K
where U k = a b ( e ξ 1 ) + cos ( ξ a b a k π ) d ξ denotes the cosine series coefficients. By using a basic calculus calculation, one has the analytical solutions of U k as follows:
U k = 1 1 + k π ( b a ) 2 k π b a sin a k π b a cos a k π b a + e a + a b k π sin a k π b a k 0 , e a a 1 k = 0
.
According to Fang and Oosterrlee [32], for some constant L , the integration domain [ a ,   b ] is chosen as follows:
[ a ,   b ] = c 1 L | c 2 | + | c 4 | ,   c 1 + L | c 2 | + | c 4 |
where c j   ( j = 1 ,   2 ,   4 ) can be obtained by the following:
c j = 1 i j j ( ln φ F ε 1 , ε 2 ( u , τ ) ) u j u = 0   ( j = 1 ,   2 ,   4 ) .
We present the following Algorithm 1 for approximating a forward starting put option price.
Algorithm 1. The COS-based algorithm for pricing a forward starting put option
Step 1: Initialization K ,   S ,   t 0 , T ,   r , κ 1 ,   θ 1 ,   σ 1 , υ 1 ,   ρ 1 ,   κ 2 ,   θ 2 ,   σ 2 ,   υ 2 ,   ρ 2 ,   λ ,   m ,   n ,   p ,   p 1 ,
q 1 ,   θ ^ 1 ,   η 1 ,   θ ^ 2 ,   η 2 ,   H 1 ,   H 2 ,   L ,   N
Step 2: Choose ε 1 ,   ε 2 to approximate the pricing model (1) by applying model (10)
Step 3: Compute the forward characteristic function by applying Formula (16)
Step 4: Compute the cumulant c 1 ,   c 2 ,   c 4 by applying Theorem 2
Step 5: Compute cosine series coefficients U k by applying Formula (20)
Step 6: Compute the truncated domain [ a ,   b ] by applying Formula (21)
Step 7: Approximate the value of a forward starting put option by applying Formula (19)

5. Numerical Experiments

Under the approximate model (10), we used the COS-based algorithm to evaluate forward starting put options. According to [28], the parameters selected were as follows: κ 1 = 12 , θ 1 = 0.05 , σ 1 = 0.9 , ρ 1 = 0.5 , υ 1 = 0.05 ,   κ 2 = 16 ,   θ 2 = 0.03 ,   σ 2 = 0.9 ,   ρ 2 = 0.5 , υ 2 = 0.02 ,   p = 0.4 ,   p 1 = 1.3 ,   q 1 = 1.2 ,   θ ^ 1 = 20 ,   θ ^ 2 = 20 ,   η 1 = 50 , η 2 = 50 ,   H 1 = 0.8 ,   H 2 = 0.7 ,   λ = 1 , m = 2 , n = 2 ,   r = 0.0165 ,   S = 100 ,   K = 100 , t 0 = 1 ,   T = 5 . To test the convergence of the pricing algorithm, we set four groups ( ε 1 ,   ε 2 ) : (0.01, 0.01), (0.001, 0.001), (0.0001, 0.0001), (0.00001, 0.00001), specifying L = 5 , L = 10 , L = 15 , and specifying N = 16 , N = 32 , N = 64 for each L .
To test the effectiveness of the pricing algorithm, we also evaluated forward starting put options by using a Euler-Maruyama scheme for the approximate model (10) as follows:
S t ε 1 , ε 2 = S t ε 1 , ε 2 e ( r λ δ 1 2 j = 1 2 υ j t 1 ε j ) Δ t + j = 1 2 υ j t 1 ε j Δ t W j t S + d k = 1 N t ( ζ k 1 ) j = N t 1 N t e Y j , υ 1 t ε 1 = υ 1 t 1 ε 1 + ( H 1 1 / 2 ) ψ 1 t 1 σ 1 υ 1 t 1 ε 1 + κ 1 θ 1 υ 1 t 1 ε 1 Δ t    + ε 1 H 1 1 / 2 σ 1 υ 1 t 1 ε 1 Δ t ρ 1 W 1 t S + 1 ρ 1 2 Z 1 t , υ 2 t ε 2 = υ 2 t 1 ε 2 + ( H 2 1 / 2 ) ψ 2 t 1 σ 2 υ 2 t 1 ε 2 + κ 2 θ 2 υ 2 t 1 ε 2 Δ t    + ε 2 H 2 1 / 2 σ 2 υ 2 t 1 ε 2 Δ t ρ 2 W 2 t S + 1 ρ 2 2 Z 2 t .
where Y j is a random number generated by two double exponential distributions, W 1 t S ,   W 2 t S ,   Z 1 t ,   Z 2 t are independent standard normal random variables, Δ t = T / N 2 denotes the time step. ψ j t   t N 1 k = 0 N 1 1 t 1 k N 1 + ε j H j 3 2 z j k ( j = 1 , 2 ) , where { z j k } denotes the sequence of independent standard normal random variables and N 1 is the number of subintervals in [ 0 ,   t ] . For Monte Carlo simulation, we used the number of subintervals N 1 = 100 , the number of simulations M = 100 , 000 , and the number of time steps N 2 = 1000 . Table 1 examines the convergence and accuracy of the pricing algorithm for evaluating forward starting puts under approximate model (10).
From Table 1, we can see that the accuracy of the COS-based algorithm depends on the values of ( ε 1 ,   ε 2 ) , L , and the number of terms N in series expansion. As the value of ( ε 1 ,   ε 2 ) decreases, option prices tend to stabilize. The effect of L on option prices depends on the number of terms N . When N = 64 , the effect of L on option prices is almost negligible. Table 1 shows that when ( ε 1 ,   ε 2 ) = ( 0.00001 ,   0.00001 ) and N = 64 , the option values calculated by the two methods are the very closest.
Furthermore, we examined the effects of the fractional Brownian motion, jump, and the second volatility component on forward starting put option prices. To this end, we set ( ε 1 ,   ε 2 ) = ( 0.00001 ,   0.00001 ) , N = 64 , L = 10 and specified two determination time t 0 = 1/4, 1, and two maturities T = 1/2, 5 years. We evaluated forward starting put options under model (10) (FDHestonMEM) applying the model parameters and option parameters used in Table 1. For comparison, we also evaluated forward starting put options under the FDHeston mode by setting λ = 0 , the DHeston model by setting λ = 0 , H 1 = H 2 = 0.5 , the FHestonMEM model by setting κ 2 = θ 2 = σ 2 = υ 2 = ρ 2 = 0 , the FHeston model by setting κ 2 = θ 2 = σ 2 = υ 2 = ρ 2 = λ = 0 , the DHestonMEM model by setting H 1 = H 2 = 0.5 , the HestonMEM model by setting κ 2 = θ 2 = σ 2 = υ 2 = ρ 2 = 0 ,   H 1 = 0.5 , and the Heston model by setting κ 2 = θ 2 = σ 2 = υ 2 = ρ 2 = λ = 0 ,   H 1 = 0.5 . Table 2 compares forward starting put option prices generated from the above eight models under t 0 = 1 / 4 ,   T = 1 / 2 and t 0 = 1 ,   T = 5 , respectively.
From Table 2, we can see that for the above eight models, the mixed-exponential jump and the second volatility component both increase forward starting options prices for all maturities and determination time. The effects of fractional Brownian motion on forward starting option prices depend on the determination time, maturity, and the value of K / S . To comprehensively examine the effects of fractional Brownian motion, jump, and the second volatility component on forward starting put option prices, we computed the relative error between forward starting put options prices under the above models with the same parameters as in Table 2. Figure 1, Figure 2 and Figure 3 report the main outcomes.
To examine the effect of mixed-exponential jump on forward starting put options prices, we investigated four cases including the single volatility model and two volatilities model driven by fractional Brownian motion and standard Brownian motion, respectively. Figure 1 reports the effect of the mixed-exponential jump on the values of short-term ( t 0 = 1 / 4 ,   T = 1 / 2 ) and long-term ( t 0 = 1 ,   T = 5 ) forward starting put options.
From Figure 1, we can see that for the short-term options, when K / S < 1.05, the effect of jump on option prices with single volatility is more significant than that with two volatilities. Moreover, the effects of jump on option prices under the framework of fractional Brownian motion are more significant than those in the standard Brownian motion case. When K / S 1.05 , the effect of jump on option prices with single volatility and with two volatilities are almost the same. For the long-term options, the effects of the mixed-exponential jump on the options prices with single volatility are both more significant than the counterpart with two volatilities under the two frameworks of Brownian motion. However, these effects of jump on the options prices are significantly weaker than those in the short-term case.
To examine the effects of fractional Brownian motion on forward starting put options prices, we investigated four cases including the single volatility model and two volatilities model with jump and without jump, respectively. Figure 2 reports the effects of fractional Brownian motion on the values of short-term ( t 0 = 1 / 4 ,   T = 1 / 2 ) and long-term ( t 0 = 1 ,   T = 5 ) forward starting put options.
From Figure 2 we can see that for the short-term options, when K / S 0.95, fractional Brownian motion both decreases option prices with single volatility and with two volatilities, and the effects of fractional Brownian motion on option prices with single volatility are more significant than those with two volatilities. When K / S > 0.95, fractional Brownian motion both increase the option prices with single volatility and with two volatilities, and the effects of fractional Brownian motion on option prices with single volatility and with two volatilities are almost the same. For the long-term options, when K / S < 0.9 or K / S > 1 , the effects of fractional Brownian motion on the options prices with single volatility are both more significant than the counterpart with two volatilities under two frameworks with jump and without jump. When 0.95 K / S < 1 , the effects of fractional Brownian motion on the options prices with two volatilities are slightly significant than the counterpart with single volatility. However, these effects of fractional Brownian motion on the options prices are significantly weaker than those in the short-term case.
To examine the effects of the second volatility on forward starting put options prices, we investigated four cases including the jump-diffusion model and diffusion model under the framework of fractional Brownian motion and standard Brownian motion, respectively. Figure 3 reports the effects of the second volatility on the values of short-term ( t 0 = 1 / 4 ,   T = 1 / 2 ) and long-term ( t 0 = 1 ,   T = 5 ) forward starting put options.
From Figure 3 we can see that for the short-term options, when K / S < 0.95 , the effects of the second volatility on option prices without jump are more significant than those with jump. As K / S increases, the intensity of the effect gradually weakens. When K / S 0.95 , the effects of the second volatility on option prices with jump and without jump are almost the same. For the long-term options, the effects of the second volatility on the options prices are significantly weaker than those in the short-term case. For all maturities and determination time, the effects of the second volatility on option prices in the framework of fractional Brownian motion are more significant than those in the standard Brownian motion case, which implies that the second volatility is necessary to price accurately forward starting options under the framework of fractional Brownian motion.
Furthermore, based on the model parameters and option parameters used in Table 2, we examined the effects of the main parameters of the second fractional volatility on the short-term ( t 0 = 1 / 4 ,   T = 1 / 2 ) forward starting put options prices. Due to limited space, Figure 4 only reports the effects of Hurst index H 2 and long-run mean θ 2 on the options prices.
From Figure 4 we can see that when H 2 > 0.5 , increasing the value of H 2 rapidly raises the options prices and the intensity of this effect gradually tends to stabilize. When H 2 < 0.5 , the effect of H 2 on the options prices is similar to that in the case of H 2 > 0.5 , and this effect is more significant. Figure 4 displays this compared to the standard Brownian motion ( H 2 = 0.5 ) case, which shows that increasing or decreasing the value of H 2 both significantly change forward starting options prices. Increasing θ 2 raises the options prices and the intensity of this effect gradually increases as the strike prices increase. Figure 4 shows that the effects of H 2 and θ 2 on forward starting options prices are both significant, which implies that the effect of the second fractional volatility on forward starting options prices is significant and the FDHestonMEM model may fit the forward implied volatility better.

6. Conclusions

This paper proposes an effective pricing algorithm for forward starting options under double fractional Heston stochastic volatilities and jumps. We express the value of a forward starting option in terms of the expectation of the forward characteristic function of log return. By introducing two small perturbed parameters, we approximate the pricing model using a semimartingale. Then, we rewrite the forward characteristic function as a conditional expectation of the proportion characteristic function which is expressed in terms of the solution to a classic PDE. With the affine structure of the approximate model, we obtained the solution to the PDE. Based on the obtained forward characteristic function and Fourier transform technique, we propose a pricing algorithm for forward starting options. For comparison, we also developed a simulation scheme for evaluating forward starting options. The numerical results show that the proposed pricing algorithm is effective. Via exhaustive comparative experiments on eight models, we found that the effects of fractional Brownian motion, mixed-exponential jump, and the second volatility component on forward starting option prices are all significant, and especially, the second fractional volatility is necessary to price accurately the forward starting options under the framework of fractional Brownian motion. The pricing algorithm for forward starting options proposed in the paper can be extended to the case of stochastic interest rates.

Author Contributions

Conceptualization, S.Z.; Formal analysis, S.Z. and H.X.; Methodology, S.Z. and H.Y.; Validation, H.Y. and H.X.; Writing—original draft, S.Z.; Writing—review & editing, S.Z. and H.X. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the National Natural Science Foundation (No. 11601420) and the Natural Science Foundation of Shaanxi Province (No. 2024JC-YBMS-008).

Data Availability Statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.

Acknowledgments

The authors would like to thank the editors and anonymous reviewers for their detailed and valuable suggestions and comments, which improved the quality of the original manuscript.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Rubinstein, M. Pay now, choose later. Risk 1991, 4, 44–47. [Google Scholar]
  2. Kruse, S.; Nögel, U. On the pricing of forward starting options in Heston’s model on stochastic volatility. Financ. Stoch. 2005, 9, 233–250. [Google Scholar] [CrossRef]
  3. Ahlip, R.; Rutkowski, M. Forward start options under stochastic volatility and stochastic interest rates. Int. J. Theor. Appl. Financ. 2011, 12, 209–225. [Google Scholar] [CrossRef]
  4. Lin, S.; He, X.J. A closed-form pricing formula for forward start options under a regime-switching stochastic volatility model. Chaos Solitons Fractals 2021, 144, 110644. [Google Scholar] [CrossRef]
  5. Guo, J.H.; Hung, M.W. A generalization of Rubinstein’s “Pay now, choose later”. J. Futures Mark. 2008, 28, 488–515. [Google Scholar] [CrossRef]
  6. Ahlip, R.; Park, L.A.F.; Prodan, A.; Weissenhofer, S. Forward start options under Heston affine jump-diffusions and stochastic interest rate. Int. J. Financ. Eng. 2021, 8, 2150005. [Google Scholar] [CrossRef]
  7. Lucié, V. Forward-start options in stochastic volatility models. Wilmott Mag. 2003, 5, 72–75. [Google Scholar] [CrossRef]
  8. Haastrecht, A.V.; Pelsser, A.A.J. Accounting for stochastic interest rates, stochastic volatility and a general dependency structure in the valuation of forward starting options. J. Futures Markets. 2011, 31, 103–125. [Google Scholar] [CrossRef]
  9. Zhang, S.M.; Geng, J.H. Fourier-cosine method for pricing forward starting options with stochastic volatility and jumps. Commun. Stat.–Theory Methods. 2017, 46, 9995–10004. [Google Scholar] [CrossRef]
  10. Zhang, S.M.; Sun, Y.D. Forward starting options pricing with double stochastic volatility, stochastic interest rates and double jumps. J. Comput. Appl. Math. 2017, 325, 34–41. [Google Scholar] [CrossRef]
  11. Hata, H.; Liu, N.L.; Yasuda, K. Expressions of forward starting option price in Hull–White stochastic volatility model. Decis. Econ. Financ. 2021, 45, 101–135. [Google Scholar] [CrossRef]
  12. Comte, F.; Renault, E. Long memory in continuous-time stochastic volatility models. Math. Financ. 1998, 8, 291–323. [Google Scholar] [CrossRef]
  13. Hassani, H.; Yarmohammadi, M.; Mashhad, L.M. Uncovering hidden insights with long-memory process detection: An in-depth overview. Risks 2023, 11, 113. [Google Scholar] [CrossRef]
  14. Han, Y.; Li, Z.; Liu, C. Option pricing under the fractional stochastic volatility model. Anziam J. 2021, 63, 123–142. [Google Scholar]
  15. Najafi, A.; Mehrdoust, F. Conditional expectation strategy under the long memory Heston stochastic volatility model. Commun. Stat.—Simul. Comput. 2023. [Google Scholar] [CrossRef]
  16. Alhagyan, M.; Yassen, M.F. Incorporating stochastic volatility and long memory into geometric Brownian motion model to forecast performance of standard and Poor’s 500 index. AIMS Math. 2023, 8, 18581–18595. [Google Scholar] [CrossRef]
  17. Alhagyan, M. The effects of incorporating memory and stochastic volatility into GBM to forecast exchange rates of Euro. Alex. Eng. J. 2022, 61, 9601–9608. [Google Scholar] [CrossRef]
  18. Hassler, U.; Pohle, M.O. Forecasting under long memory. J. Financ. Econom. 2023, 21, 742–778. [Google Scholar] [CrossRef]
  19. Gatheral, J.; Jaisson, T.; Rosenbaum, M. Volatility is rough. Quant. Financ. 2018, 18, 933–949. [Google Scholar] [CrossRef]
  20. Livieri, G.; Mouti, S.; Pallavicini, A.; Rosenbaum, M. Rough volatility: Evidence from option prices. IISE Trans. 2018, 50, 767–776. [Google Scholar] [CrossRef]
  21. Fukasawa, M. Volatility has to be rough. Quant. Financ. 2021, 21, 1–8. [Google Scholar] [CrossRef]
  22. Brandi, G.; Di Matteo, T. Multiscaling and rough volatility: An empirical investigation. Int. Rev. Financ. Anal. 2022, 84, 102324. [Google Scholar] [CrossRef]
  23. Coqueret, G.; Tavin, B. An investigation of model risk in a market with jumps and stochastic volatility. Eur. J. Oper. Res. 2016, 253, 648–658. [Google Scholar] [CrossRef]
  24. Jin, X.; Hong, Y. Jump-diffusion volatility models for variance swaps: An empirical performance analysis. Int. Rev. Financ. Anal. 2023, 87, 102606. [Google Scholar] [CrossRef]
  25. Bates, D.S. Empirical option pricing models. Annu. Rev. Financ. Econ. 2022, 14, 369–389. [Google Scholar] [CrossRef]
  26. Wang, L.; Xia, W.X. Power-type derivatives for rough volatility with jumps. J. Futures Mark. 2022, 42, 1369–1406. [Google Scholar] [CrossRef]
  27. Cai, N.; Kou, S.G. Option pricing under a mixed-exponential jump diffusion model. Manag. Sci. 2011, 57, 2067–2081. [Google Scholar] [CrossRef]
  28. Zhang, S.; Yong, H.; Xiao, H. Option pricing with fractional stochastic volatilities and jumps. Fractal Fract. 2023, 7, 680. [Google Scholar] [CrossRef]
  29. Wang, K.; Guo, X.X. Valuations of variance and volatility swaps under double Heston jump-difusion model with approximative fractional stochastic volatility. Comput. Econ. 2024, 63, 1543–1573. [Google Scholar] [CrossRef]
  30. Thao, T.H. An approximate approach to fractional analysis for finance. Nonlinear Anal.–Real World Appl. 2006, 7, 124–132. [Google Scholar] [CrossRef]
  31. Euch, O.E.; Rosenbaum, M. The characteristic function of rough Heston models. Math. Financ. 2019, 29, 3–38. [Google Scholar] [CrossRef]
  32. Fang, F.; Oosterlee, C.W. A novel pricing method for European options based on Fourier cosine series expansions. SIAM J. Sci. Comput. 2008, 31, 826–848. [Google Scholar] [CrossRef]
Figure 1. The effect of the mixed-exponential jump on the values of short-term ( t 0 = 1 / 4 ,   T = 1 / 2 ) and long-term ( t 0 = 1 ,   T = 5 ) forward starting put options.
Figure 1. The effect of the mixed-exponential jump on the values of short-term ( t 0 = 1 / 4 ,   T = 1 / 2 ) and long-term ( t 0 = 1 ,   T = 5 ) forward starting put options.
Fractalfract 08 00283 g001
Figure 2. The effects of fractional Brownian motion on the values of short-term ( t 0 = 1 / 4 ,   T = 1 / 2 ) and long-term ( t 0 = 1 ,   T = 5 ) forward starting put options.
Figure 2. The effects of fractional Brownian motion on the values of short-term ( t 0 = 1 / 4 ,   T = 1 / 2 ) and long-term ( t 0 = 1 ,   T = 5 ) forward starting put options.
Fractalfract 08 00283 g002
Figure 3. The effects of the second volatility on the values of short-term ( t 0 = 1 / 4 ,   T = 1 / 2 ) and long-term ( t 0 = 1 ,   T = 5 ) forward starting put options.
Figure 3. The effects of the second volatility on the values of short-term ( t 0 = 1 / 4 ,   T = 1 / 2 ) and long-term ( t 0 = 1 ,   T = 5 ) forward starting put options.
Fractalfract 08 00283 g003
Figure 4. The effects of Hurst parameter H 2 and long-run mean θ 2 of the second volatility on the values of short-term ( t 0 = 1 / 4 ,   T = 1 / 2 ) forward starting put options.
Figure 4. The effects of Hurst parameter H 2 and long-run mean θ 2 of the second volatility on the values of short-term ( t 0 = 1 / 4 ,   T = 1 / 2 ) forward starting put options.
Fractalfract 08 00283 g004
Table 1. Convergence and accuracy of the COS-based algorithm for evaluating forward starting put options under the approximate model.
Table 1. Convergence and accuracy of the COS-based algorithm for evaluating forward starting put options under the approximate model.
ε 1 ε 2 LCOS(16)COS(32)COS(64)MC ± Std
518.568218.568218.5682
0.010.011018.613518.568218.568218.5686 ± 0.0159
1519.105818.570018.5682
518.583818.583818.5838
0.0010.0011018.627418.583818.583818.5842 ± 0.0162
1519.110918.585518.5838
518.591818.591818.5918
0.00010.00011018.634418.591818.591818.5920 ± 0.0150
1519.113318.593418.5918
518.596018.596018.5960
0.000010.000011018.638118.596018.596018.5960 ± 0.0146
1519.114518.597518.5960
Table 2. Comparison of forward starting put options prices under the FDHestonMEM, FDHeston, FDHestonMEM, FDHeston, DHestonMEM, DHeston, HestonMEM, and Heston model.
Table 2. Comparison of forward starting put options prices under the FDHestonMEM, FDHeston, FDHestonMEM, FDHeston, DHestonMEM, DHeston, HestonMEM, and Heston model.
t 0 TKFDHestonMEMFDHestonFHestonMEMHestonDHestonMEMDHestonHestonMEMFHeston
800.33770.29630.11520.23380.51630.48060.26330.0808
850.83570.77100.37890.52921.03380.97680.58380.3135
901.77201.68601.02011.11991.92131.84071.20840.9212
953.29613.19722.28172.21183.32423.22392.33422.1609
1/41/21005.49705.39674.35014.05215.36895.25994.19084.2280
1058.37948.28857.26966.82698.12028.01726.94997.1656
11011.872311.797610.934910.513711.551911.467610.597510.8588
11515.858515.802215.155014.862515.552315.492014.907515.1063
12020.208720.169519.732119.571119.964619.926519.591619.7044
809.49819.17606.26246.02949.55889.24436.40065.8740
8511.503311.15747.96407.636111.518811.17828.04977.5372
9013.694113.32819.89079.463813.660613.29719.91509.4315
9516.061615.678812.035711.509615.976915.593711.992811.5507
1510018.596018.199914.389713.767818.459118.059614.276813.8855
10521.287220.881116.942116.230321.098520.686016.758816.4249
11024.125223.712119.680818.887323.886123.463819.429019.1565
11527.100126.682722.593221.728026.812826.383622.276922.0673
12030.202329.783225.666624.740529.869529.436325.290925.1438
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhang, S.; Xiao, H.; Yong, H. Forward Starting Option Pricing under Double Fractional Stochastic Volatilities and Jumps. Fractal Fract. 2024, 8, 283. https://doi.org/10.3390/fractalfract8050283

AMA Style

Zhang S, Xiao H, Yong H. Forward Starting Option Pricing under Double Fractional Stochastic Volatilities and Jumps. Fractal and Fractional. 2024; 8(5):283. https://doi.org/10.3390/fractalfract8050283

Chicago/Turabian Style

Zhang, Sumei, Haiyang Xiao, and Hongquan Yong. 2024. "Forward Starting Option Pricing under Double Fractional Stochastic Volatilities and Jumps" Fractal and Fractional 8, no. 5: 283. https://doi.org/10.3390/fractalfract8050283

APA Style

Zhang, S., Xiao, H., & Yong, H. (2024). Forward Starting Option Pricing under Double Fractional Stochastic Volatilities and Jumps. Fractal and Fractional, 8(5), 283. https://doi.org/10.3390/fractalfract8050283

Article Metrics

Back to TopTop