Adopting Feynman–Kac Formula in Stochastic Differential Equations with (Sub-)Fractional Brownian Motion

: The aim of this work is to establish and generalize a relationship between fractional partial differential equations (fPDEs) and stochastic differential equations (SDEs) to a wider class of stochastic processes, including fractional Brownian motions { B Ht , t ≥ 0 } and sub-fractional Brownian motions { ξ Ht , t ≥ 0 } with Hurst parameter H ∈ ( 12 ,1 ) . We start by establishing the connection between a fPDE and SDE via the Feynman–Kac Theorem, which provides a stochastic representation of a general Cauchy problem. In hindsight, we extend this connection by assuming SDEs with fractional-and sub-fractional Brownian motions and prove the generalized Feynman–Kac formulas under a (sub-)fractional Brownian motion. An application of the theorem demonstrates, as a by-product, the solution of a fractional integral, which has relevance in probability theory.


Introduction
Consider the Cauchy problem [1] of the following parabolic partial differential equation (PDE) on R d ∂ ∂t u(x, t) = κ ∂ 2 ∂x 2 u(x, t) + ηB H (t), where u(x, t) ∈ C 2,1 , u 0 (x) is a bounded measurable function and B H (t) is a fractional Brownian motion (cf. Section 2). Without loss of generality, we assume that the parameter κ is constant. This second-order PDE has a stochastic representation for η = 0, according to the Feynman-Kac formula [2,3]. Indeed, we obtain if x t satisfies Equation (3) and the function σ(x t , t) is sufficiently integrable where B H t is a Brownian motion (BM) if the Hurst parameter is of H = 1 2 [4][5][6]. Additionally, the problem of (1) has an intimate relationship to the fractional partial differential equation (fPDE) [7]: Note that this equation contains a fractional derivative in general or a semi-derivative in respect of time in special [8][9][10][11][12][13].
There is a large amount of the literature devoted to each issue of the Cauchy problem [6,14]. This research closes a gap by considering the linking relationships of (sub-)fractional Brownian motions as well as fPDEs. The Feynman-Kac formula (2) provides a unique weak solution to Equation (1). Different versions of the Feynman-Kac formula have been discovered for a variety of problems [15,16]. Some generalizations of the Feynman-Kac formula are discovered by Querdiane and Silva [17] and Hu et al. [18,19]. A Feynman-Kac formula also exists for Lévy processes by Nualart and Schoutens [20]. Advancements in stochastic differential equations and fractional partial differential equations to analyse complex systems are related to our research [21][22][23][24]. Furthermore, recent developments in fractional calculus contributed to a better understanding and further studies of the relationships between fractional PDEs and stochastic calculus [25][26][27][28][29][30][31]. However, we are concerned about the linkage of the Cauchy problem and the representation by a fPDE, as well as the Feynman-Kac formula. For the Cauchy problem, we generalize the stochastic representation of Feynman-Kac by utilizing fractional Brownian motion (fBM) with Hurst parameter H > 1/2.
In addition, the more recent literature looks at the idea of sub-fractional Brownian motion (sub-fBM). A sub-fBM is an intermediate between a Brownian motion and fractional Brownian motion. The existence and properties, such as long-range dependence, self-similarity and non-stationarity were introduced by Bojdecki et al. [32] and Tudor et al. [33,34]. Since the sub-fractional Brownian motion is not a martingale, methods of stochastic analysis are more sophisticated. However, several authors developed stochastic calculus and integration concepts for an fBM [25] and sub-fBM [35][36][37]. Recently, for a sub-fractional Brownian motion with Hurst parameters H > 1 2 , a maximal inequality was established according to the Burkholder-Davis-Gundy inequality for fractional Brownian motion [38]. It turns out that fBM and sub-fBM are adequate stochastic processes in scientific applications [13,39].
In this paper, our purpose is to construct and prove a general link of the Cauchy problem with the Feynman-Kac equation via Itô's formula for fBM and sub-fBM. Consequently, this paper links the solution of u(x, t) defined by Equation (1) with the stochastic Feynman-Kac representation to a fractional Brownian motion {B H t } and sub-fBM {ξ H t }. We prove the result and show the properties of (sub-)fractional processes in stochastic analysis. Note that, throughout this paper, we frequently assume 1 2 < H < 1. The paper is organized as follows. Section 2 contains preliminaries on fractional calculus, particularly fractional Brownian motion. Thereafter, we examine sub-fractional stochastic processes and integration rules in Section 3. Here, we list the definitions and assumptions for the remainder of the article. In Section 4, we link the Cauchy problem to the Feynman-Kac formula with stochastic differential equations driven by fractional and sub-fractional Brownian motions. We state our theorems and prove our statements. In Section 5, we examine the Cauchy problem and the relationship to fractional partial differential equations (fPDE). Furthermore, we find a new fractional derivative and integral with relevance in probability theory. The conclusion is in Section 6.

Preliminaries
In the following section, we define preliminary concepts on fractional stochastic processes and fractional calculus.

Fractional Calculus
Since we deal with the Hurst parameter H, we need to know fractional calculus. Let a, b ∈ R, a < b. Let f ∈ L 1 (a, b) and α > 0. The left-and right-sided fractional integral of f of order α are defined for x ∈ (a, b), respectively, as This is the fractional integral of Riemann-Liouville type. Similarly, the fractional leftand right-sided derivative, for f ∈ I α a (L p ) and 0 < α < 1, are defined by and for all x ∈ (a, b) and I α a (L p ) is the image of L p (a, b). It is easy to see that if f ∈ I 1 a (L 1 ), Note

Fractional Stochastic Process
Mandelbrot and van Ness defined a fractional Brownian Motion (fBM), B H t , as a Brownian motion, B(t), together with a Hurst parameter (or Hurst index) H ∈ (0, 1) in 1968 [8]. The new feature of fBM's is that the increments are interdependent. The latter property is defined as self-similarity. A self-similar process has invariance with respect to changes in timescale (scaling-invariance). Almost all other stochastic processes, such as the standard Brownian Motion or Lévy processes, likely have independent increments. They create the famous class of Markov processes. Empirically, there is ubiquitous evidence in science that fractional stochastic processes, for instance, spectral densities with a sharp peak, are related to the phenomena of long-range interdependence over time. Indeed, the observable presence of interdependence in many real-world applications calls for fractional stochastic processes.

Definition 1.
Let H be 0 < H < 1 and B 0 an arbitrary real number. We call B H (t, ω) a fractional Brownian Motion (fBM) with Hurst parameter H and starting value B 0 at time 0, such as (1) B H (0, ω) = B 0 , and; (3) [Or equivalently by the Riemann-Liouville fractional integral: We immediately obtain the corollary.
For values of H, such as 0 < H < 1 2 or 1 2 < H < 1 the fBM B H (t, ω) has different properties. If 0 < H < 1 2 , we say that it has the property of short memory. Indeed, Mandelbrot and van Ness [8] shows that this range is associated with negative correlation. If 1 2 < H < 1, then the fBM has the property of long-memory or long-range dependence with time-persistence (Mandelbrot and van Ness [8]). Alternatively, we define a fractional Brownian motion by where H ∈ (0, 1) denotes the Hurst parameter.

Remark 1.
The covariance is trivially derived by starting with a standard Brownian motion and extending it with the Hurst index H, such as where, for H = 1 2 , we obtain the Brownian motion. The covariance is derived by the following steps

Corollary 2. The expectation of non-overlapping increments of an fBM is
Thus, we can see that the expected increments are non-zero. Indeed, the increments are interdependent, contrary to Markov processes. The second part of the variance is Proof. Part (1): For t 1 < t 2 < t 3 < t 4 , the equality of the covariance function implies that Y := B H (t 2 ) − B H (t 1 ) has the same distribution as X : where t 1 < t 2 and t 3 < t 4 with ∆t = t 2 − t 1 = t 4 − t 3 . Hence, the incremental behavior at any point in the future is the same. Thus, we say that it has stationary increments. Part (2) hence, we obtain (B H (at)) 2 = a 2H (B H (t)) 2 and this equal to B H (at) = a H B H (t). The proof of part (3) is already in Corollary 2.

Itô's Formula for Fractional Brownian Motion
A fractional Brownian motion is continuous but almost certainly not differentiable [8]. Hence, it is inconvenient that an fBM does not have a derivative or integral. Furthermore, the fBM is neither a martingale nor a semi-martinagle. Therefore, Itô calculus is not applicable to fractional Brownian Motions if H = 1 2 . However, stochastic calculus was developed with respect to fractional Brownian motion by [40] and the stochastic integral was introduced by [25]. The theory is a fractional extension of Itô-calculus, but limited to a Hurst index H ∈ (1/2, 1). If H > 1/2 the fBM exhibits long-range dependence, which is a fundamental property in physics or finance.
By utilizing Wick calculus that has zero mean and explicit expressions for the second moment, we define the stochastic fractional integral, satisfying the property E[ Further, the functions f and g belong to the Hilbert space with the inner product defined by This machinery leads to an analogue Itô formula for a fractional Brownian process. Already, Alòs et al. [41] proved this result under certain conditions for Itô's formula. Theorem 1. (Alòs et al., 2001). Let f be a function of class C 2 (R), satisfying the growth condition where c and λ are positive constants and λ < 1 4 T −2H . Suppose that B H = {B H t , t ∈ [0, T]} is a zero mean continuous Gaussian process whose covariance function R f BM (t, s) is of the form in Equation (8). Then, the process F (B H t ) belongs to a Hilbert space and, for each t ∈ [0, T], the following Itô's formula holds: However, we utilize a result by Duncan et al. [25], which is more convenient in our case. Here, is the Itô-Duncan theorem for a fractional Brownian motion: Theorem 2. (Duncan et al., 2000, Thm 4.1, p. 596). If f : R → R is a twice continuously differentiable function with bounded derivatives to order two, i.e., f ∈ C 2 , then

Remark 2.
If H = 1 2 , we obtain, from Theorem 2, the usual Itô formula for a Brownian motion Similarly, for a function with two parameters f (t, B H t ), a generalized rule exists according to Duncan et al. [25].
For the proof, we refer to Duncan et al. [25]. If F(s) = a(s) is a deterministic function; then, the rule simplifies. Let η t = t 0 a u dB H u , where a ∈ L 2 φ ; then, we obtain If a s ≡ 1, then we obtain Itô's formula, such as in Theorem 2 and in Equation (13).

Sub-Fractional Stochastic Process
A sub-fractional Brownian motion (sub-fBM) is an intermediate between a Brownian motion and fractional Brownian motion. It is a more general, self-similar Gaussian process or a generalization of a fBM. The sub-fBM has the property of H-self-similarity and longrange dependence, such as the fBM, yet it does not have stationary increments [32].
It is well-established that a stochastic process is uniquely determined by its covariance function Cov(ξ H t , ξ H s ). Thus, we define:

Definition 3. A sub-fractional Brownian motion of Hurst parameter H is a centered mean zero
where ξ H 0 = 0 and E[ξ H t ] = 0.
If H = 1 2 , it coincides with a Brownian motion on R + with covariance Cov(ξ H t , ξ H s ) = s ∧ t := min[s, t]. The process ξ H t has the following integral representation for H > 1 2 (see [41]): Hence, the sub-fractional Brownian motion has a kernel of Note that the kernel has similarities to the fBM, as in Equation (9). Next, we discuss the main properties of a sub-fBM: Lemma 1. Let ξ H t be a sub-fBM for all t. It has the following properties: = ξ H u − ξ H s , i.e., the increments are non-stationary.
Part 2. Given property 1, one immediately obtains Part 3. Let s = 0 and t = h > 0, then The difference in both increments is Finally, we prove two differences of fBM and sub-fBM.

Proposition 2.
Let B H t be a fractional Brownian motion and ξ H t be a sub-fractional Brownian motion. For H ∈ ( 1 2 , 1) the following holds: where, only for s = t = 0 or s = 0, t = 0, we obtain equality.

Itô's Formula for Sub-Fractional Brownian Motion
For a Hurst parameter H > 1 2 , the stochastic integral of a sub-fBM T 0 f (t)dξ H t exists. The following theorem holds and is proven by [42]: According to Yan et al. ([36], Theorem 3.2 on p. 139) Itô's formula under a sub-fBM can be computed as follows: Theorem 5. (Yan et al., 2011) Let f ∈ C 2 (R) and H ∈ ( 1 2 , 1). Then, we have Details of the proof are given in ( [36], pp. 139-140). The authors even extend Itô's formula to d−dimensional sub-fBM and convex functions f : ξ H t → R.

Linking Cauchy via Feynman-Kac to SDEs with fBM and Sub-fBM
Next, we derive the link between the Cauchy problem (1) and the stochastic representation according to Feynman-Kac by Equation (2). Consider a stochastic process x s on the time interval [t, T] as the solution to the SDE in Equation (3). Next, use the Dynkin operator or Fokker-Planck operator A defined by After integration and under the assumption that u(x, t) satisfies Equation (22). The time integrals will vanish. Given x t ∈ C 2 and a deterministic σ, we obtain, after taking the expectation and the property that the stochastic integral vanishes, the stochastic representation as follows: If H = 1 2 , the stochastic representation simplifies to the well-known Feynman-Kac formula Next, we state the Feynman-Kac formula for our Cauchy problem (1), given a subfractional Brownian motion.
if x t ∈ C 2 . Note, for H = 1 2 , we obtain the same as in Theorem 6.
The proof follows an equal argument as above in the proof of Theorem 6.

Cauchy Problem and Fractional-PDE
Next, we demonstrate the direct linkage for the Cauchy-problem (1) to the fPDE in Equation (4). In step one, we compute the Laplace transform of the right-hand side of the heat equation: whereū(x, t) := L[u(x, t)]. Thus, we obtain This is a second-order ordinary differential equation in the x−variable. The solution is u(x, t) = c * e − √ sx for some constant c. Determining the constant by the second-derivativē u xx = c * se − √ sx shows that c = 1. In step two, we compute the first-derivative of the solution x t e − x 2 2t . The final term solves the fPDE in Equation (30). Thus, the fractional integral for α = 1 2 must be equal to the probability density function N (x) in order to satisfy the fPDE in Equation (30).

Conclusions
This article studies the relationships of the Cauchy problem (1) and relates them to fractional partial-differential equations, as well as to the stochastic representations by the Feynman-Kac formula with a generalized fractional and sub-fractional Brownian motion with Hurst parameter H > 1/2. In addition, we find fractional derivatives and integrals in relation to the Gaussian probability function by utilizing the novel insight into the linkage of the Cauchy problem and fPDE. This vantage point is of importance in probability theory, fractional calculus and stochastic theory. In future research, we intend to extend our theorems to Hurst parameters H < 1/2 and the stochastic Cauchy problem under a sub-fBM.