A Fractional Heston-Type Model as a Singular Stochastic Equation Driven by Fractional Brownian Motion

: This paper introduces the fractional Heston-type ( fHt ) model as a stochastic system comprising the stock price process modeled by a geometric Brownian motion. In this model, the infinitesimal return volatility is characterized by the square of a singular stochastic equation driven by a fractional Brownian motion with a Hurst parameter H ∈ ( 0,1 ) . We establish the Malliavin differentiability of the fHt model and derive an expression for the expected payoff function, revealing potential discontinuities. Simulation experiments are conducted to illustrate the dynamics of the stock price process and option prices.


Introduction
Allowing volatility to be stochastic in a financial market model was one of the great achievements in the history of quantitative finance.This innovation led to stochastic volatility modeling, as previously discussed by Heston [1] and several other researchers, addressing shortcomings in the standard Black-Scholes model (see, for example, Alòs et al. [2] for a summary).In the context of Heston [1], the stock price process is described by a geometric Brownian motion of the form: where η and (Y t ) t≥0 represent a constant drift and the stochastic variance of the instantaneous rate of return dX t := dS t /S t , respectively.The stochastic process (Y t ) t≥0 takes the form of the standard Cox-Ingersoll-Ross process, satisfying the following stochastic differential equation: The parameter θ represents the speed of reversion of the stochastic process (Y t ) t≥0 towards its long-run mean µ, and the parameter ν represents the volatility of the stochastic variance (Y t ) t≥0 .The Brownian motions (B t ) t≥0 and ( Bt ) t≥0 are assumed to be correlated.This model is well known in the literature as the Heston model.
It has recently been demonstrated that volatility and the volatility of volatility exhibit rough behaviors.This implies that the paths tend to be rougher, showing short-range dependency and can be effectively modelled using fractional Brownian motion with a Hurst parameter H < 1/2.For further insights, refer to studies such as Alos et al. [3], Fukasawa [4], Gatheral et al. [5], Livieri et al. [6], Takaishi [7], Fukasawa [8], Brandi and Di Matteo [9], and related findings.
On the other hand, volatility persistence is also associated with long-memory properties, indicated by the slow decay of the autocorrelation function.In this regard, Comte and Renault [10] demonstrated that long-memory stochastic volatility models are better suited to reproduce the gradual flattening of implied volatility skews and smiles observed in financial market data.This finding has been corroborated and tested by several other researchers, including Chronopoulou and Viens [11], Chronopoulou and Viens [12], Tripathy [13], and subsequent studies.In contrast to rough volatility, this modeling approach involves using fractional Brownian motion with a Hurst parameter H > 1/2.
From the above, we may notice a contradiction regarding whether volatility is rough or exhibits long-range dependency, a subject of debate in the literature.However, Alòs and Lorite [14] observed that both properties are not mutually exclusive.A process can exhibit both long and short dependency properties, with each dominating at different scales, and consequently, at different maturities in the implied volatility surface.This idea is supported by Funahashi and Kijima [15], who demonstrated that if the volatility where Y H t and Y H ′ t are the fractional Ornstein-Uhlenbeck process driven by fractional Brownian motions with Hurst parameters H > 1/2 and H ′ < 1/2, respectively, then Y H t does not have an impact on the ATM short-time limit skew.
To incorporate both roughness and long-range dependency properties into the Heston model, the standard Brownian motion is replaced by a fractional Brownian motion (W H t ) t≥0 with a Hurst parameter H > 1/2 or H < 1/2, resulting in the fractional Heston model.For further details, refer to studies such as Alòs and Yang [16], Mishura and Yurchenko-Tytarenko [17], Mehrdoust and Fallah [18], Tong and Liu [19], Richard et al. [20], along with the references therein.In this context, the volatility process can be represented by: The stochastic process (Y t ) t≥0 is well-known as the fractional Cox-Ingersoll-Ross (fCIR) process.The Equation ( 1) is well defined only when H ∈ (2/3, 1) as shown by Mishura et al. [21].This limitation was overcome by defining fCIR process as a square of a stochastic process with additive fBm.In other words, the stochastic process (Y t ) t≥0 can be written as the square of a stochastic process (Z t ) t≥0 that verifies The stochastic volatility process (Y t ) t≥0 described above gives rise to complex models in option pricing or risk analysis, which are not easily manageable, particulary when the volatility drift is not constant.In this paper, we propose a general form of the fractional Heston-type model in a simple and natural manner.As before, we assume that the stock price process is driven by a geometric Brownian motion (S t ) t≥0 , satisfying the following stochastic differential equation: where σ(Y t ) represents the volatility of the infinitesimal log-return dX t := dS t /S t with (Y t ) t≥0 a fractional Cox-Ingersoll-Ross (fCIR) process that captures both long and short-range dependency.We adopt the definition provided by Mishura and Yurchenko-Tytarenko [22] or Mpanda et al. [23] and describe the stochastic process (Y t ) t≥0 as where the stochastic process (Z t ) t≥0 is referred to as a general form of the fCIR process that satisfies the following differential equation: and τ is the first time the process (Z t ) t≥0 hits zero, defined by It was shown in Mpanda et al. [23] that the stochastic process (Y t ) t≥0 satisfies This implies that the function f (t, z) can be defined as the drift of the volatility process (Y t ) t≥0 .Additionally, the stochastic process (W H t ) t≥0, H∈(0,1) is well-known as fractional Brownian motion (fBm) with Hurst parameter H, defined as a centered Gaussian process with a covariance function: There exist several representations of fBm.Nourdin [24] summarized them.In financial volatility modelling, the Volterra representation is widely used, particulary due to its simplicity.In this representation, the fBm is written in terms of a standard Brownian motion in the time interval [0, t] as follows: where (V t ) t∈[0,T] is a standard Brownian motion and where κ H (s, t) is a square integrable kernel that may take different forms.One effective expression is through the Euler hypergeometric integral with κ H (s, t) given by with Γ(•) and 2 F 1 (a, b, c; d) being the gamma and Gaussian hypergeometric functions, respectively.A truncated expression of the kernel (9) was suggested by Decreusefond and Üstünel [25] where This representation is referred to as a Type II fBm or Riemann-Liouville fBm.Here, κ H (s, t) is defined by This representation was used by Gatheral et al. [5] to model rough volatility.The standard Brownian motions (B t ) t∈[0,T] and (V t ) t∈[0,T] are assumed to be correlated, meaning there exists ρ ∈ [−1, 1] such that E B t V t = ρt.This implies that there exists a Brownian motion ( Ṽt ) t∈[0,T] independent to (V t ) t∈[0,T] , that is E V t , Ṽt = 0, such that A new and natural way of defining a fractional Heston-type (fHt) model as singular stochastic differential equation driven by a fractional Brownian motion is through the following stochastic system: The existence of the stochastic process (Z t ) t≥0 in (5) was previously discussed by Nualart and Ouknine [26].They proposed that for H < 1/2, the drift function g(t, z) := f (t, z)z −1 must satisfy the linear growth condition, and for H > 1/2, g(t, z) must verify the Hölder continuity condition.Additionally, particular cases of fHt model (12) has been previously investigate by Alòs and Yang [16], Mishura and Yurchenko-Tytarenko [17] show that the fHt model is free of arbitrage.The above stochastic model is also a generalisation of the rough volatility model previously discussed by Gatheral et al. [5] with σ(y) = e y , f (t, y) = −θy and the Type II fBm of the form (10) with small Hurst parameters.
The remainder of this paper is structured as follows: Section 2 constructs approximating sequences of stock prices and fCIR processes.Malliavin differentiability within the fHt model is discussed in Section 3. Finally, Section 4 derives the expected payoff function and performs simulations of option prices.

Approximating Sequences in fHt Model
The main purpose of introducing approximating sequences of both fractional volatility and stock price processes relies on their positiveness.The following theorems discuss the positiveness of (Z t ) t≥0 and before this, we consider the following assumption below.
Under this assumption, the following theorems follow.Theorem 1.Let (Z t ) t≥0 be a stochastic process that verifies (5) Proof.Here, we highlight the proof of this theorem by contradiction, and we refer the reader to Mpanda et al. [23] for a complete and comprehensive proof.Let τ(ω) = inf{t ≥ 0 : Z t (ω) = 0} be the first time that the process (Z t ) t∈[0,T] hits zero and τ ε (ω) = sup{t ∈ (0, τ(ω)) : Z t (ω) = ε} be the last time (Z t ) hits ε before reaching zero.In addition, define where g(t, z) = f (t, z)z −1 and z 0 is any point chosen between zero and the initial value of (Z t ) t∈[0,T] , that is, 0 < z 0 < Z 0 and ζ = inf{t ≥ 0 : Z t (ω) = z 0 }.Then, from (5), we have: and by Hölder continuity of fBm, we have which yields the following inequality On the other hand, the critical point ẑ of the function F ε (z) is given by One may notice that there exists 0 This goes in contradiction with (13).
Theorem 2. Consider for each k > 0, the stochastic process (Z Then, for any T > 0 and H < 1/2, Proof.The proof can be carried out as previously shown by defining τ For more details, refer to Mpanda et al. [23].

Approximating Sequences of (Z t ) t≥0
Inspired by Alos and Ewald [28], we construct an approximating sequence (Z ε t ) t≥0, ε>0 of the fCIR process that satisfies the following differential equation: where the function Λ ε (z) in ( 14) is defined by It is easy to verify that Λ ε (z) > 0 for all ε > 0. As a straight consequence, the drift of (Z ε t ) t≥0, ε>0 is also positive.In addition, The next step is to show that for every t ≥ 0, the sequence Z ε t converges to Z t in L p as ε → 0. Proposition 3. The sequence of estimated random variables Z ε t converges to Z t in L p (Ω) for all p ≥ 1.

Proof.
Case 1. H = 1/2.This case was discussed previously by Alos and Ewald [28] (Proposition 2.1) and can be easily extended to the case where Λ ε (z) is defined by (15).
Case 2. For H > 1/2, the dominated convergence theorem shall be applied.Firstly, we need to show the pointwise convergence of the approximated stochastic process (Z ε t ) t≥0 towards (Z t ) t≥0 , that is lim ε→0 Z ε t = Z t .For this, let τ ε (ω) = inf{t ≥ 0 : Z t (ω) ≤ ε} be the first time the process (Z t ) t≥0 hits ε .Since the sample paths of the stochastic process (Z t ) t≥0 are positive everywhere almost surely as in Theorem 1, then the stochastic process (Z t ) t≥0 up to stopping time τ ε .Then, for all t ∈ [0, τ ε ] and using the definition of Λ ε (z) given by ( 16), Z τ ε t = Z ε t almost surely when ε → 0 since the drift function f (t, z) is monotonic.Again, the positiveness of (Z t ) t≥0 means that lim ε→0 Z τ ε t = Z t a.s.We may conclude that lim ε→0 Z τ ε t = lim ε→0 Z ε t = Z t almost surely and for all t ≥ 0.
On the other hand, the result from Hu et al. [29] (Theorem 3.1) shows that for a fixed T > 0 and for all p ≥ 1, are nonrandom constants depending on parameters γ, β, T, and This result also implies that It follows that sup t∈[0,T] Z ε t (ω) } ∈ L p (Ω) which yields the desired L p convergence.Case 3.For H < 1/2, we consider a sequence of an increasing drift function ) t≥0 as follows: where Λ(z) is defined by (15) and ) t≥0 hits ε, then from Theorem 2, for any fixed T > 0, P(ω ) t≥0 remains positive up to time T which is not necessary equal to infinity unlike the previous case.After using similar arguments for Case 2, one may conclude that lim ε→0 To achieve this, we borrow some ideas from Mishura and Yurchenko-Tytarenko [30].Firstly, let Z0 be a small positive value less than the initial value Z 0 such that 0 < Z0 < Z 0 and let τ 1 = τ 1 (ε, ω) be the last time the stochastic process (Z ε t ) t≥0, ε>0 hits (or before hits) Z0 , that is, Technically, there exists a constant M ≥ 2 such that Z0 = Z 0 M .Now we can consider two cases: t ∈ [0, τ 1 ] and t ∈ (τ 1 , T].
By applying the Callebaut's inequality theorem, it will be easy to show that for all p ≥ 1, From ( 17), we may deduce that Since the drift function satisfies the linear growth condition, this means there exists a positive constant k such that f (t, z) ≤ k(1 + |z|).It follows that Inequalities ( 19)-(21) yield the following: On the other hand, recall that g., Nourdin [24]) and since From the Grönwall-Bellman inequality theorem, we obtain which can be shortly written as |Z ε t | p ≤ C, where C = C(r, k, T, Z 0 , ν, H) is a nonrandom constant in parameters r, k, T, Z 0 , ν, and H taking the following form where and Then we have: As previously stated, the integral in the last inequality of ( 24) can be expressed as follows On the other hand, we may observe that It follows that, From this expression, we may also conclude that |Z ε t | p ≤ C, where C = C(C 1 , C 2 ) and where C 1 and C 2 are a non-random constants defined by ( 22) and (23), respectively.This shows that E Z ε t p < ∞ and consequently, E sup t∈[0,T] Z t p < ∞.
This concludes the proof of the proposition.
Assumption 2. The volatility function σ(y) is strictly positive and Lipschitz continuous.
Corollary 4.Under Assumption 2, and for any p ≥ 1, Proof.This follows immediately from the previous proposition.
Remark 1.One may use similar arguments to Mishura and Yurchenko-Tytarenko [30] to show that the stochastic process (Z ε t ) t≥0, ε>0 is strictly positive almost surely for all H ∈ (0, 1).Consequently, it is also well suitable for rough volatility processes, that is, a fractional volatility process with H < 1/2.

Approximating Sequences of Stock Price Process (S t ) t≥0
With (Z ε t ) t≥0, ε>0 , let us construct the approximating sequence (S ε t ) t≥0, ε>0 of the stock price process (S t ) t≥0 defined by the following geometric Brownian motion: where 2 , with (Z ε t ) t≥0, ε>0 the approximating sequence that satisfies (14).The next step is to show that S ε t converges to S t in L p , p ≥ 1.
Proposition 5. Set X t := log S t and X ε t := log S ε t .Then, the sequence X ε t converges to X t in L p (Ω) for all p ≥ 1.
Proof.Firstly, we have from Itô formula that where X 0 := log S 0 .Then, for some non-random constant C > 0, one may have: Then it follows firstly that T 1 → 0 from Corollary 4. To analyse convergence of T 2 , the Burkholder-Davis-Gundy inequality can be used and one may deduce that which also converges to zero from Corollary 4. It follows that lim ε→0 sup t≥0 X ε t − X t p = 0, ∀p > 0 that implies the desired L p convergence of X ε t to X t and S ε t to S t .
1.The stochastic process (S t ) t≥0 is a unique solution of a geometric Brownion motion of the form that can be found using the standard Itô formula, yielding: 2. The approximated stochastic volatility and stock price processes will be compulsory for H ≤ 1/2 and optional for H > 1/2.However, for the sake of consistency, we shall use the approximated sequences ( 14) with ε = 0 for H > 1/2 and with ε > 0 for H ≤ 1/2.

Malliavin Differentiability
Nowadays, the application of Malliavin calculus in stochastic volatility modelling has increased, particularly due to the introduction of fBm, which exacerbates the complexity of derivative pricing models.In this section, we discuss the Malliavin differentiability of both the stock price process and its stochastic volatility as given in the stochastic system (12).This analysis will pave the way for the first application, which involves deriving the expected payoff function.

Preliminaries on Malliavin Calculus for fBm
Malliavin calculus, initially introduced by Paul Malliavin in the 1970s, is a powerful mathematical tool used to analyze stochastic processes and their associated functionals.It provides a systematic framework for differentiating stochastic processes with respect to underlying Brownian motions, enabling the analysis of complex stochastic systems.Malliavin calculus finds wide application in quantitative finance for pricing and hedging financial derivatives, risk management, and portfolio optimization.In this section, we provide some preliminaries on Malliavin calculus for fBm.for a complete background in Malliavin calculus, we refer the reader to Nualart [31] and Da Prato [32].For applications in quantitative finance, we refer to Nunno et al. [33] and Alòs and Lorite [14].
On the time interval [0, T], consider the Hilbert space H constructed with the closure of the set of real-valued step functions on [0, T] denoted by E with respect to the scalar product ⟨1 [0,t] , 1 [0,s] ⟩ := ψ(s, t).If the fBm (W H t ) t≥0 takes the Volterra representation (8), then its covariance function is given by where κ H (•, r) is a Kernel taking the form ( 9) or (10).The covariance function (30) can be further developed to reach the following expression: where α H = H(2H − 1).For any step function ϕ 1 , ϕ 2 ∈ E , one may generalise the above as The mapping 1 [0,t] → W H t can be extended to an isometry ϕ → W H (ϕ) between the Hilbert space H and the Gaussian space denoted by H 1 spanned by W H . Now, consider the operator κ * H : E → L 2 ([0, T]) that provides the previous isometry between H and H 1 := L 2 ([0, T]) defined by then, for any where κ H is the operator defined by: The space H H represents the fractional version of the Cameron-Martin space.Additionally, let ψ H : H → H H be the operator defined as: Moreover, it is important to note that ψ H is Holder continuous of order H.This leads to the definition of the Malliavin derivative as presented by Nualart [31].
Definition 1.Consider a space S of smooth random variables of the form The domain of D (denoted by D 1,2 ) is a Sobolev space defined as the closure of the space of smooth random variables S, with respect to the norm: The directional Malliavin derivative is defined as the scalar product In other words, the directional Malliavin derivative ⟨DF, ϕ⟩ H is the derivative at ε = 0 of the smooth random variable F composed with the shifted process (W H (ϕ ′ ) + ε⟨ϕ ′ , ϕ⟩ H ).
The Malliavin derivatives have several important properties.Bouleau and Hirsch [34] showed that if F ∈ D 1,2 and ||DF|| H > 0 a.s., then the law of F has a density with respect to the Lebesgues measure on [0, T].In addition, the following properties were also proved in Nualart [31]: (1) Integration by parts, in the sense that for all ϕ ∈ H, (2) Chain rule, that is, for , then the smooth function (3) The future Malliavin derivative of an adapted process is zero, that is, for all r > t, D r F t = 0 a.s.
As an example, the computation of Malliavin derivative of fBm (represented by ( 8)) with respect to Brownian motion (V t ) t≥0 is given by D V s W H t = κ H (t, s), and for s > t (that is, the future derivative), D V s W H t = 0. We may write this shortly as When the kernel κ H (t, s) is defined as in (10), the computation of the Malliavin derivative is straightforward and is given by: Since fBm can be represented in terms of standard Brownian motion, it is worth mentioning how the Malliavin derivative of stochastic differential equations is computed.There are several approaches in the literature; one of them was previously presented by Detemple et al. [35] based on a transformation to a volterra integral equation, and was proven to be more efficient in the numerical estimation of Malliavin derivatives.In this approach, consider a stochastic differential equation of the form: and set For more details see Detemple et al. [35], Mishura and Yurchenko-Tytarenko [17] and Alòs and Lorite [14].The computation of the Malliavin derivative of the standard Heston model was presented by Alos and Ewald [28], where this approach was utilised.With the above background, we can now discuss the Malliavin derivatives of both stochastic processes (S t ) t≥0 and (Z t ) t≥0 of the stochastic system (12).

Differentiability of Stochastic
Processes (Z t ) t≥0 and (S t ) t≥0 Proposition 6.In the stochastic system (12), the law of stochastic processes (Z t ) t≥0 and (S t ) t≥0 is absolutely continuous with respect to the Lebesgue measure over any finite interval [0, t].In addition, set Then, S t , Z t ∈ D 1,2 , and and where D B Y s = 2 ρ Z s D V Z s , with D V Z s given by (45).
Proof.First, let us show the expression (45).We have As in Detemple et al. [35], we set U t = D V u Z t , then we retrieve the following Volterra integral equation of the second kind with the kernel function F(s, Z s ) and unknown function U t .For any u < t, this equation takes the form: to which the solution is given by which yields (45).Similarly, to find the expression (47), we may compute U t := D W H u Z t through the Volterra integral equation of the form: The derivation of (47) is straightforward; we just have to apply the integration by parts formula and chain rules to (X t ), which takes the following form: We obtain To demonstrate the absolute continuity of (Z t ) t≥0 with respect to the Lebesgue measure over any finite interval [0, t], we first note that the solution to the stochastic differential Equation ( 14) takes the following form: where the function ϱ t (with ϱ 0 = 0) belongs to the class C λ ([0, t]), λ < H.This expression arises from the intricate nature of the set {t ≥ 0|Z t = 0}, which is associated with the level sets of fBm particularly for small Hurst parameters.For further insights, refer to Mukeru [36] and Mishura and Yurchenko-Tytarenko [30].We have: The Taylor expansion of (51) yields: where for some θ ε s ∈ [0, 1].The solution to Equation (52) is obtained using Expression (31).We have: which yields the following:

H
This holds almost surely in L 2 (Ω), and consequently: as previously, and It follows that Z t ∈ D 1,2 , and consequently, according to Bouleau and Hirsch [34], the law of stochastic process (Z t ) t≥0 is absolutely continuous with respect to the Lebesgue measure over any finite interval [0, t].Similar reasoning can also be applied to the stock price process (S t ) t≥0 .

Expected Payoff Function
The aim of this section is to derive the expected payoff function E[h(S T )] by using some results from Malliavin calculus.We follow Altmayer and Neuenkirch [37] closely.

Differentiability of Expected Payoff Function
Let h : R → R be the payoff function that satisfies the following assumption.Proof.Firstly, it is straightforward to check that E[L 2 (S T )] < ∞ since L(x) also verifies the linear growth condition and the law of stock price process (S t ) t∈[0,T] are bounded almost surely.On the other hand, since L verifies Assumption 3 and the sample paths of the stock price process (S t ) t∈[0,T] is absolutely continuous with respect to the Lebesgue measure on R (See Proposition 6), then from the chain rule formula for Malliavin derivatives, we may deduce The first inequality is due to Holder inequality and the finiteness of the last expression makes sense since S t ∈ D 1,2 as discussed previously.It follows that ||L|| 1,2 < ∞, which concludes the proof.
Lemma 8. Let h(x), x ∈ R be a payoff function that satisfies Assumption 3 and denote h(e x ) := g(x) with its antiderivative G(x) that also satisfies the Lipschitz condition.Set and where X T := log S T and Proof.We follow the idea of Altmayer and Neuenkirch [37].To establish the equality (54), we rewrite E[g(X T )] as From Proposition 7, we may deduce that G(X T ) ∈ D 1,2 and We now obtain In addition, from Proposition 6, and since the integral du is well defined from Assumption 2, then we have: and defining I T by (53), we obtain (54).To establish (55), we rewrite the function G(x) (which is the antiderivative of g(x)) as follows where C is a constant taking the form C = 1 0 h(u)du and by using the standard integration by part formula, one may obtain With this setting, we have

Some Simulations 4.2.1. Simulations of Stock Price Process
For simulating the stock price process, one can use the Euler-Maruyama approximation scheme.This involves dividing the time interval [0, T] into N sub-intervals of equal length, where 0 = t 0 , t 1 , • • • , t N = T with t i = iT/N and the lag ∆t = T/N.The estimated stock price at time t i denoted by ( Ŝt i ) i=1,••• ,N and the estimated volatility ( Ŷt i ) i=1,••• ,N are, respectively, given by: where , respectively, the increment of Brownian motions V t∈[0,10] , Ṽt∈[0,T] , and fBm W H t∈[0,T] .In addition, fBm is represented by the Volterra stochastic integral (7), which can be discretized as follows: for all j = 1, • • • , N; i = 0, • • • , j and where δV i = V i − V i−1 is the increment of standard Brownian motion with W H t 0 = 0. Here, κ H (t j , s) is a discretised square integrable kernel (9) given by As an illustrative example, the following figures represent 10 sample paths of the stock price process on the interval [0, T] with N = 1000, ρ = 0.6, X 0 = 100, η = r = 0.05, ν = 0.1, σ( Ŷt i ) = 0.8 Ŷt i + 0.1.The drift of the fractional volatility process is defined by with κ = 1, c = 2. Referring to remark (2), we will choose ε = 0 when H > 1/2, and for H ≤ 1/2, we set ε = 0.01 as shown in Figure 1 below.To the scheme (57), we associate the discrete approximation of the integral I T provided below, which will be used in the computation of the expression (55).

ÎT =
1 Firstly, we consider the stochastic process (Z t ) t≥0 defined as a fractional Ornstein-Uhlenbeck process, specifically with f (t, z) = −θz 2 , where θ is a positive parameter, ν = 2 and H > 1/2.Under these settings, one may recover the model discussed by Bezborodov et al. [27] with Y t = Z 2 t instead.In this case, the volatility process may not necessarily be positive almost surely, as it violates the Assumption 1, rendering Theorems 1 and 2 inapplicable.To address this, the volatility function σ(y) is chosen to be strictly positive.In addition, we define the payoff function h(x) as a combination of European and binary options with the same strike price K and time to maturity T, that is, The expression for L(S T ) can be readily deduced from (56) as follows: The payoff function h(S T ) and L(S T ) can be visualised in the Figure 2 below with strike price K = 0.5.4) and ( 5) cannot be used as the drift function f (t, z) = −θz 2 < 0 for all θ > 0, we consider the direct form of the stochastic volatility (Y t ) t≥0 that takes the form of the Ornstein-Uhlenbeck process satisfying the following differential equation: In this case, we observe that the values of option prices are not significantly different for ρ = 0 and H ≥ 1/2 from Bezborodov et al. [27].The option prices increase or decrease when ρ is positive or negative, respectively.
Next, consider the fractional volatility process described by a standard fCIR process, that is, with f (t, z) = µ − θz 2 and correlation ρ between infinitesimal returns and volatility, the option prices are simulated with ρ = 0.5 and µ = 0.1.
The following tables present the mean prices and their corresponding coefficients of variation for a European-Binary option with the payoff function defined by Equation (62) under different Hurst parameters.The mean values are obtained from an average of 10 3 trials, along with their respective coefficients of variation and option prices are calculated using the expected payoff function discounted by the net present value.We employ two approaches to compute the expectation of the payoff function: 1.The direct expectation E[h(S T )] = E[(S T − K) + + 1 S T >K ], for a fixed strike price K = 0.5.2. The formula obtained via Malliavin derivative, as given by Equation (55), with L(S T ) defined by (63).
In Table 1, the payoff values at the maturity date T are obtained by performing 10 4 simulations.We observe that the direct estimation of expected values tends to stabilize starting from N = 8000 (where N represents the number of steps between 0 and T).For example, with H = 0.1 and σ(Y t ) = √ Y t + 0.1, the expected values for different values of N are represented in Figure 3.In Table 2 below, we perform simulations again, this time using the expressions (55), (62), and (63).We may observe that the expected option values are slightly different from those in the previous table.However, the following table demonstrates that the values of the expected payoff function stabilize starting from N = 4000.One of the main reasons attributed to these satisfactory observations is that the expression of the expectation (55) includes a continuous functional of the stock price process (S t ) t∈[0,T] along with a weight term (1 + I T ) that is independent of the functional.This property is even more efficient for discontinuous payoff functions h(S T ).As in the previous example, choose the Hurst H = 0.   We shall use Z 0 = 1, ν = 0.4, c = 0.02, θ = 1.To keep positiveness of the stochastic process (Z t ) t≥0 for all H ∈ (0, 1), we shall rather use its approximated stochastic process (Z ε t ) t, ε≥0 defined by (14), that is where the function Λ ε (z) is defined by with ε = 0.01 for H ≤ 1/2 and ε = 0 for H > 1/2.As previously stated, the fBm is simulated by using the formula (58) and (59).We perform again 10 3 trials for 10 4 simulations and various time-steps on the time interval [0, 1].We get the mean of option prices with their corresponding coefficient of variations for different volatility functions σ(y) under the European-Binary option as given in Table 3 for direct estimations and in Table 4 by using (55).Note that observations from the previous sections also apply to this one.

Conclusions
In this paper, we have constructed the fractional Heston-type model as a stochastic system comprising the stock price process (S t ) t≥0 modeled by a geometric Brownian motion.The volatility of this process is represented as a strictly positive and Lipschitz continuous function σ(Y t ) of fractional Cox-Ingersoll-Ross process (Y t ) t≥0 , which is characterized by the square of a stochastic process (Z t ) t≥0 that satisfies a stochastic differential equation with additive fractional Brownian motion.
To ensure the positivity of the stochastic process (Z t ) t≥0 for all Hurst parameters H ∈ (0, 1), we have considered an approximating sequence (Z ε t ) t≥0,,ε>0 converging to (Z t ) t≥0 in L p (Ω) for all p ≥ 1.This construction also enables us to demonstrate that (Z t ) t≥0 , (S t ) t≥0 , and the payoff function h(S t ) are Malliavin differentiable.Furthermore, we establish that the law of the stochastic processes (Z t ) t≥0 and (S t ) t≥0 is absolutely continuous with respect to the Lebesgue measure over any finite interval [0, t].
To support our findings, we conducted simulations.Firstly, we modeled volatility using the Ornstein-Uhlenbeck process, corroborating the results found in Bezborodov et al. [27].Secondly, we explored the fractional Cox-Ingersoll-Ross process with time-varying parameters.Our observations indicate that option prices exhibit greater stability under the expected value of option prices obtained through Malliavin calculus.

Assumption 3 .
The payoff function h : R → R and its antiderivative denoted by L(x) (such that L(x) = h(x)) are bounded and verify the Lipschitz condition.Proposition 7. L(S T ) ∈ D 1,2 .

Figure 2 .
Figure 2. The payoff function h(S T ) and L(S T ) with K = 0.5.Now, to find expected values of the payoff function, we use the same parameters (η = r = 0.2, θ = 0.6, T = 1) with different forms of volatility process σ(Y t ) of the infinitesimal log-return process dS t /S t as in Bezborodov et al. [27].Since the fCIR process of the form (4) and (5) cannot be used as the drift function f (t, z) = −θz 2 < 0 for all θ > 0, we consider the direct form of the stochastic volatility (Y t ) t≥0 that takes the form of the Ornstein-Uhlenbeck process satisfying the following differential equation:

Figure 3 .
Figure 3. Values of mean option prices for different N under direct estimations.
1 and the volatility function σ(Y t ) = √ Y t + 0.1, the expected values for different values of N under the formula (55) are represented in Figure 4.

Figure 4 .
Figure 4. Values of mean option prices for different N through Formula (55).
, Bezborodov et al. [27] for H > 1/2.One can use the same idea of Bezborodov et al. [27] (Theorem 3.6) to Then, the Malliavin derivative of F denoted by DF is a H − valued random variable given by DF

Table 1 .
Option prices under direct estimations.

Table 3 .
Option prices using direct estimations.