Large Time Behavior on the Linear Self-Interacting Diffusion Driven by Sub-Fractional Brownian Motion II: Self-Attracting Case

In this study, as a continuation to the studies of the self-interaction diffusion driven by subfractional Brownian motion S H , we analyze the asymptotic behavior of the linear self-attracting diffusion: d X t H = d S t H − θ ∫ 0 t ( X t H − X s H ) d s d t + ν d t , X 0 H = 0 , where θ > 0 and ν ∈ R are two parameters. When θ < 0, the solution of this equation is called self-repelling. Our main aim is to show the solution X H converges to a normal random variable X ∞ H with mean zero as t tends to infinity and obtain the speed at which the process X H converges to X ∞ H as t tends to infinity.


INTRODUCTION
In a previous study (I) (see [12]), as an extension to classical result, we considered the linear selfinteracting diffusion as follows: In the present study, we consider the case θ > 0 and study its large time behaviors.
Let us recall the main results concerning the system (Eq. 1). When H 1 2 , as a special case of path-dependent stochastic differential equations, in 1995, Cranston and Le Jan [8] introduced a linear self-attracting diffusion (Eq. 1) with θ > 0. They showed that the process X t converges in L 2 and almost surely as t tends infinity. This path-dependent stochastic differential equation was first developed by Durrett and Rogers [10] introduced in 1992 as a model for the shape of a growing polymer (Brownian polymer). The general form of this kind of model can be expressed as follows: where B is a d-dimensional standard Brownian motion and f is Lipschitz continuity. X t corresponds to the location of the end of the polymer at time t. Under some conditions, they established asymptotic behavior of the solution of the stochastic differential equation. The model is a continuous analog of the notion of edge (respectively, vertex) self-interacting random walk (see, e.g., Pemantle [22]). By using the local time of the solution process X, we can make it clear how the process X interacts with its own occupation density. In general, Eq. 2 defines a self-interacting diffusion without any assumption on f. We call it self-repelling (respectively, self-attracting) if, for all x ∈ R d , x · f(x) ≥ 0 (respectively, ≤ 0). More examples can be found in Benaïm et al. [2,3], Cranston and Mountford [9], Gan and Yan [11], Gauthier [13], Herrmann and Roynette [14], Herrmann and Scheutzow [15], Mountford and Tarr [20], Sun and Yan [26,27], Yan et al [34], and the references therein.
In this present study, our main aim is to expound and prove the following statements: (I) For θ > 0 and 1 2 < H < 1, the random variable exists as an element in L 2 , where the function is defined as follows: (II) For θ > 0 and 1 2 < H < 1, we have X H t → X H ∞ in L 2 and almost surely as t → ∞.
(III) For θ > 0 and 1 (IV) For θ > 0 and 1 Then the convergence holds in L 2 as T tends to infinity. This article is organized as follows. In Section 2, we present some preliminaries for sub-fBm and Malliavin calculus. In Section 3, we obtain some lemmas. In Section 4, we prove the main results given as before. In Section 5, we give some numerical results.

PRELIMINARIES
In this section, we briefly recall the definition and properties of stochastic integral with respect to sub-fBm. We refer to Alós et al [1], Nualart [21], and Tudor [31] for a complete description of stochastic calculus with respect to Gaussian processes.
As we pointed out in the previous study (I) (see [12]), the sub-fBm S H is a rather special class of self-similar Gaussian processes such that S H 0 0 and for all s, t ≥ 0. For H 1/2, S H coincides with the standard Brownian motion B. S H is neither a semimartingale nor a Markov process unless H 1/2, so many of the powerful techniques from stochastic analysis are not available when dealing with S H . As a Gaussian process, it is possible to construct a stochastic calculus of variations with respect to S H . The sub-fBm appeared in Bojdecki et al [4] in a limit of occupation time fluctuations of a system of independent particles moving in R d according a symmetric α-stable Lévy process. More examples for sub-fBm and related processes can be found in Bojdecki et al. [4][5][6][7], Li [16][17][18][19], Shen and Yan [23,24], Sun and Yan [25], C. A. Tudor [32], Tudor [28][29][30][31], C. A. Tudor [33], Yan et al [33,35,36], and the references therein.
The normality and Hölder continuity of the sub-fBm S H imply that t1S H t admits a bounded p H variation on any finite interval with p H > 1 H . As an immediate result, one can define the Young integral of a process u {u t , t ≥ 0} with respect to a sub-fBm S H φ(t)φ(s)ψ H (t, s)dsdt.
Let now D and δ be the (Malliavin) derivative and divergence operators associated with the sub-fBm S H . And let D 1,2 denote the Hilbert space with respect to the norm as follows: holds for any F ∈ D 1,2 and D 1,2 ⊂ Dom(δ). Moreover, for any u ∈ D 1,2 , we have  for an adapted process u, and it is called the Skorohod integral. By using Alós et al [1], we can obtain the relationship between the Skorohod and the Young integral as follows: provided u has a bounded q variation with 1 ≤ q < 1 H and u ∈ D 1,2 such that

SOME BASIC ESTIMATES
For simplicity, we throughout let C stand for a positive constant which depends only on its superscripts, and its value may be different in different appearances, and this assumption is also suitable to c. Recall that the linear self-attracting diffusion with sub-fBm S H is defined by the following stochastic differential equation:  with θ > 0. The kernel (t, s)1h θ (t, s) is defined as follows: for s, t ≥ 0. By the variation of constants method (see, Cranston and Le Jan [8]) or Itô's formula, we may introduce the following representation: for t ≥ 0. The kernel function (t, s)1h θ (t, s) with θ > 0 admits the following properties (these properties are proved partly in Cranston and Le Jan [8]): • For all s ≥ 0, the limit exists.
• For all t ≥ s ≥ 0, we have h θ (s) ≤ h θ (t, s), and • For all t ≥ s, r ≥ 0 and θ ≠ 0, we have • For all t > 0, we have Lemma 3.1. Let 1 2 < H < 1 and θ > 0. Then the random variable exists as an element in L 2 .
Proof. This is a simple calculus exercise. In fact, we have for all θ > 0 and 1 2 < H < 1. Clearly, Eq. 10 implies that and for all θ > 0 and 1 2 < H < 1. These show that the random variable X H ∞ exists as an element in L 2 .
Lemma 3.2. Let θ > 0. We then have Proof. This is a simple calculus exercise. In fact, we have for all t ≥ 0 and θ > 0. Noting that lim t→∞ t e we see that by L'Hopital's rule.
for all t ≥ 0.
where we have used the equation for all 0 ≤ s < t ≤ T, where C and c are two positive constants depending only on H, θ, ] and T.
Proof. The lemma is similar to Lemma 3.5 in the previous study (I).
Lemma 3.6. Let θ > 0 and 1 2 ≤ H < 1. Then the convergence holds in L 2 and almost surely as t tends to infinity.
Proof. Convergence (18) in L 2 follows from Lemma (3.1). In fact, by Eq. 10, we have as t tends to infinity.
Frontiers in Physics | www.frontiersin.org January 2022 | Volume 9 | Article 791858 8 On the other hand, by Lemma (3.5), 3.3 and the equation S H t t → 0 almost surely as t tends to infinity, we find that as t tends to infinity. It follows from the integration by parts that almost surely as t tends to infinity.

SOME LARGE TIME BEHAVIORS
In this section, we consider the long time behaviors for X H with 1 2 < H < 1 and θ > and our objects are to prove the statements given in Section 1.
holds in L 2 and almost surely as t tends to infinity.
Proof. When H 1 2 , the convergence is obtained in Cranston-Le Jan [8]. Consider the decomposition for all t ≥ 0. We first check that Eq. 19 holds in L 2 . By Lemma 3.6 and Lemma 3.2, we only need to prove Y H t converges to zero in L 2 . It follows from the equation for all θ > 0 as t tends to infinity and Lemma 3.4 that for all θ > 0 and 1 2 < H < 1 as t tends to infinity, which implies that Eq. 19 holds in L 2 .
We now check that Eq. 19 holds almost surely as t tends to infinity. By Lemma 3.6, we only need check that Y H t converges to zero almost surely as t tends to infinity. We have for every integer n ≥ 1 and 0 ≤ k ≤ n, which implies that for any ε > 0, every integer n ≥ 1 and 0 ≤ k ≤ n.
On the other hand, for every s ∈ (0, 1), we denote Then {R n,k s , 0 ≤ s ≤ 1} also is Gaussian for every integer n ≥ 1 and 0 ≤ k ≤ n. It follows that as T tends to infinity.
in L 2 as T tends to infinity.
Proof. Given 1 2 < H < 1 and θ > 0, for all t ≥ 0. We now prove the lemma in three steps.
Step I. We claim that as t tends to infinity. Clearly, we have Thus, 29 is equivalent to Frontiers in Physics | www.frontiersin.org January 2022 | Volume 9 | Article 791858 for all 1 2 < H < 1.
Step II. We claim that as T tends to infinity. We have that for all t > s > 0. An elementary calculation may show that