Abstract
In this study, as a continuation to the studies of the self-interaction diffusion driven by subfractional Brownian motion SH, we analyze the asymptotic behavior of the linear self-attracting diffusion:
where θ > 0 and are two parameters. When θ < 0, the solution of this equation is called self-repelling. Our main aim is to show the solution XH converges to a normal random variable with mean zero as t tends to infinity and obtain the speed at which the process XH converges to as t tends to infinity.
1 Introduction
In a previous study (I) (see [12]), as an extension to classical result, we considered the linear self-interacting diffusion as follows:with θ ≠ 0, where θ and ν are two real numbers, and SH is a sub-fBm with the Hurst parameter . The solution of Eq. 1 is called self-repelling if θ < 0 and is called self-attracting if θ > 0. When θ < 0, in a previous study (I), we showed that the solution XH diverges to infinity as t tends to infinity andandin L2 and almost surely, for all n = 1, 2, …, where ( − 1)!! = 1 and
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 , 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 Xt converges in L2 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. Xt 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 (respectively, ). 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 , the random variable
exists as an element in
L2, where the function is defined as follows:
with
θ> 0.
(II) For θ > 0 and , we have
in
L2and almost surely as
t→
∞.
(III) For θ > 0 and , we have
in distribution as
t→
∞, where
(IV) For θ > 0 and , we have
Then the convergence
holds in
L2as
Ttends 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.
2 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 SH is a rather special class of self-similar Gaussian processes such that andfor all s, t ≥ 0. For H = 1/2, SH coincides with the standard Brownian motion B. SH 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 SH. As a Gaussian process, it is possible to construct a stochastic calculus of variations with respect to SH. The sub-fBm appeared in Bojdecki et al [4] in a limit of occupation time fluctuations of a system of independent particles moving in according a symmetric α-stable Lévy process. More examples for sub-fBm and related processes can be found in Bojdecki et al. [4–7], Li [16–19], Shen and Yan [23, 24], Sun and Yan [25], C. A. Tudor [32], Tudor [28–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 SH imply that admits a bounded pH variation on any finite interval with . As an immediate result, one can define the Young integral of a process u = {ut, t ≥ 0} with respect to a sub-fBm SHas the limit in probability of a Riemann sum. Clearly, when u is of bounded qH variation on any finite interval with qH > 1 and , the integral is well-defined andfor all t ≥ 0.
Let be the completion of the linear space generated by the indicator functions 1[0,t], t ∈ [0, T] with respect to the inner product:for s, t ∈ [0, T]. For every , we can define the Wiener integral with respect to SH, denoted byas a linear (isometric) mapping from onto by using the limit in probability of a Riemann sum, where is the Gaussian Hilbert space generating by andfor any . In particular, when , we can show thatwherefor s, t ∈ [0, T]. Thus, when if for every T > 0, the integral exists in L2 andwe can define the integral as follows:and
Let now D and δ be the (Malliavin) derivative and divergence operators associated with the sub-fBm SH. And let denote the Hilbert space with respect to the norm as follows:Then the duality relationshipholds for any and . Moreover, for any , we havewhere (Du)∗ is the adjoint of Du in the Hilbert space given as follows: . We denotefor 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 and such that
3 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 SH is defined by the following stochastic differential equation:with θ > 0. The kernel (t, s)↦hθ(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)↦
hθ(
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
and
• For all t > 0, we have
Lemma 3.1Letandθ > 0. Then the random variableexists as an element in L2.
ProofThis is a simple calculus exercise. In fact, we havefor all θ > 0 and . Clearly, Eq. 10 implies that and andfor all θ > 0 and . These show that the random variable exists as an element in L2.
Lemma 3.2Letθ > 0. We then have
ProofThis is a simple calculus exercise. In fact, we havefor all t ≥ 0 and θ > 0. Noting thatandwe see thatby L’Hopital’s rule.
Lemma 3.3Letθ > 0. We then havefor all t ≥ 0.
Lemma 3.4Let θ > 0 and . We then have
ProofBy L’Hopital’s rule and the change of variable , it follows thatwhere we have used the equationThis completes the proof.
Lemma 3.5Let θ > 0 and . We then havefor all 0 ≤ s < t ≤ T, where C and c are two positive constants depending only on H, θ, ν and T.
ProofThe lemma is similar to Lemma 3.5 in the previous study (I).
Lemma 3.6Let θ > 0 and . Then the convergenceholds in L2 and almost surely as t tends to infinity.
ProofConvergence (18) in L2 follows from Lemma (3.1). In fact, by Eq. 10, we haveas t tends to infinity.On the other hand, by Lemma (3.5), 3.3 and the equation almost surely as t tends to infinity, we find thatas t tends to infinity. It follows from the integration by parts thatalmost surely as t tends to infinity.
4 Some Large Time Behaviors
In this section, we consider the long time behaviors for XH with and θ > and our objects are to prove the statements given in Section 1.
Let θ > 0 and . Then the convergenceholds in L2 and almost surely as t tends to infinity.
ProofWhen , the convergence is obtained in Cranston-Le Jan [8]. Consider the decompositionfor all t ≥ 0.We first check that Eq. 19 holds in L2. By Lemma 3.6 and Lemma 3.2, we only need to prove converges to zero in L2. It follows from the equationfor all θ > 0 as t tends to infinity and Lemma 3.4 thatfor all θ > 0 and as t tends to infinity, which implies that Eq. 19 holds in L2.We now check that Eq. 19 holds almost surely as t tends to infinity. By Lemma 3.6, we only need check that converges to zero almost surely as t tends to infinity. We havefor all θ > 0 and as t tends to infinity. To obtain the convergence, we define the random sequencefor every integer n ≥ 1. Then {Zn,k, k = 0, 1, 2, …, n} is Gaussian for every integer n ≥ 1. It follows from Lemma 3.4 thatfor every integer n ≥ 1 and 0 ≤ k ≤ n, which implies thatfor any ɛ > 0, every integer n ≥ 1 and 0 ≤ k ≤ n.On the other hand, for every s ∈ (0, 1), we denoteThen also is Gaussian for every integer n ≥ 1 and 0 ≤ k ≤ n. It follows thatfor all s, s′ ∈ [0, 1]. Thus, for any ɛ > 0, by Slepian’s theorem and Markov’s inequality, one can getfor every integer n ≥ 1 and 0 ≤ k ≤ n. Combining this with the Borel–Cantelli lemma and the relationshipwe show that almost surely as t tends to infinity. This completes the proof.
Let θ > 0 and . Then the convergenceholds in distribution, where is a central normal random variable with its variance
ProofWhen , this result also is unknown. We only consider the case and similarly one can prove the convergence for . By Eq. 20, Slutsky’s theorem, and Lemma 3.2, we only need to show thatin probability andin distribution.First, Eq. 22 follows from Eq. 10 andfor all θ > 0 and as t tends to infinity.We now obtain convergence (23). By the equationas t tends to infinity and Lemma 3.4, we getfor all θ > 0 and as t tends to infinity. Thus, convergence (23) follows from the normality of for all and the theorem follows.At the end of this section, we obtain a law of large numbers. Consider the process YH defined byThen the self-attracting diffusion XH satisfiesandby integration by parts. It follows thatfor all and t ≥ 0. By the variation of constant method, we can give the explicit representation of YH as follows:
Lemma 4.1Let and θ > 0. Then we havealmost surely and in L2 as T tends to infinity.
ProofThis lemma follows from Eq. 24 and the estimatesas T tends to infinity.
Let and θ > 0. Then we havein L2 as T tends to infinity.
ProofGiven and θ > 0,for all t ≥ 0. Thenfor all t ≥ 0. We now prove the lemma in three steps.
Step IWe claim thatas t tends to infinity. Clearly, we haveThus, 29 is equivalent toBy L’Hôspital’s rule and Lemma 3.4, it follows thatfor all .
Step IIWe claim thatas T tends to infinity. We have thatfor all t > s > 0. An elementary calculation may show thatfor all t > s > 0. It follows from the equation with x ≥ 0 and β > − 1 thatfor all t > s > 0 and 0 ≤ α ≤ 1. For the term Λ2(H; t, s), by the proof of Lemma 3.4, we find thatfor all . Combining this with the equationand the equation with x > 0 and 0 < ϱ < 1, we getfor all t > s > 0, and 0 ≤ γ ≤ 2 − 2H. Thus, we have showed that the estimateholds for all t > s ≥ 0. In particular, we havefor all t, s ≥ 0. As a corollary, we getas T tends to infinity.
Step IIIWe claim thatas t tends to infinity. By steps I and II, we find that Eq. 37 is equivalent toas t tends to infinity. Noting that the equationfor all t, s > 0, we further find that convergence (38) also is equivalent toas T tends to infinity. We now check that convergence (40) in two cases.
Case 1Let . Clearly, by Eq. 36, we have to
Case 2Let . By Eq. 36, we have thatwith andwith for all T > 1. Similarly, by Eq. 35, we also havefor all T > 1 and since 0 < t2 − s2 < 1 for . Thus, we have shown thatwith andas T tends to infinity. This shows that convergence (40) holds for all . Similarly, we can also show the theorem holds for and the theorem follows.
Remark 1By using the Borel–Cantelli lemma and Theorem 4.3, we can check that convergence (28) holds almost surely.
5 Simulation
We have applied our results to the following linear self-attracting diffusion driven by a sub-fBm
SHwith
as follows:
where
θ> 0 and
are two parameters. We will simulate the process with
ν= 0 in the following cases:
FIGURE 1

Path of XH with θ = 1 and H = 0.7.
FIGURE 2

Path of XH with θ = 10 and H = 0.7.
FIGURE 3

Path of XH with θ = 100 and H = 0.7.
TABLE 1
| t | t | t | |||
|---|---|---|---|---|---|
| 0.000 0 | 0.000 0 | 0.343 8 | −0.121 6 | 0.687 5 | −0.097 9 |
| 0.015 6 | 0.008 7 | 0.359 4 | −0.129 0 | 0.703 1 | −0.096 8 |
| 0.031 3 | 0.011 3 | 0.375 0 | −0.140 6 | 0.718 8 | −0.101 7 |
| 0.046 9 | 0.003 9 | 0.390 6 | −0.146 7 | 0.734 4 | −0.109 0 |
| 0.062 5 | −0.015 3 | 0.406 3 | −0.145 9 | 0.750 0 | −0.108 8 |
| 0.078 1 | −0.023 8 | 0.421 9 | −0.157 9 | 0.765 6 | −0.118 8 |
| 0.093 8 | −0.022 9 | 0.437 5 | −0.162 4 | 0.781 3 | −0.116 3 |
| 0.109 4 | −0.027 0 | 0.453 1 | −0.166 6 | 0.796 9 | −0.112 5 |
| 0.125 0 | −0.033 5 | 0.468 8 | −0.170 1 | 0.812 5 | −0.123 1 |
| 0.140 6 | −0.035 3 | 0.484 4 | −0.171 7 | 0.828 1 | −0.140 0 |
| 0.156 3 | −0.037 0 | 0.500 0 | −0.173 8 | 0.843 8 | −0.146 5 |
| 0.171 9 | −0.049 8 | 0.515 6 | −0.177 4 | 0.859 4 | −0.155 4 |
| 0.187 5 | −0.054 4 | 0.531 3 | −0.176 6 | 0.875 0 | −0.160 4 |
| 0.203 1 | −0.059 3 | 0.546 9 | −0.171 3 | 0.890 6 | −0.170 9 |
| 0.218 8 | −0.076 5 | 0.562 5 | -0.166 7 | 0.906 3 | −0.174 3 |
| 0.234 4 | −0.085 0 | 0.578 1 | −0.166 4 | 0.921 9 | −0.178 1 |
| 0.250 0 | −0.098 1 | 0.593 8 | −0.152 1 | 0.937 5 | −0.179 4 |
| 0.265 6 | −0.106 2 | 0.609 4 | −0.142 2 | 0.953 1 | −0.180 3 |
| 0.281 3 | −0.112 7 | 0.625 0 | −0.139 5 | 0.968 8 | −0.175 8 |
| 0.296 9 | −0.113 2 | 0.640 6 | −0.128 2 | 0.984 4 | −0.193 5 |
| 0.312 5 | −0.114 0 | 0.656 3 | −0.123 4 | 1.000 0 | −0.198 0 |
| 0.328 1 | −0.121 4 | 0.671 9 | −0.105 4 |
Data of with θ = 1 and H = 0.7
TABLE 2
| t | t | t | |||
|---|---|---|---|---|---|
| 0.000 0 | 0.000 0 | 0.343 8 | −0.098 3 | 0.687 5 | −0.110 9 |
| 0.015 6 | −0.006 4 | 0.359 4 | −0.110 4 | 0.703 1 | −0.112 1 |
| 0.031 3 | −0.010 4 | 0.375 0 | −0.110 8 | 0.718 8 | −0.112 6 |
| 0.046 9 | −0.010 1 | 0.390 6 | −0.109 8 | 0.734 4 | −0.103 4 |
| 0.062 5 | −0.017 9 | 0.406 3 | −0.111 9 | 0.750 0 | −0.099 1 |
| 0.078 1 | −0.017 7 | 0.421 9 | −0.110 6 | 0.765 6 | −0.090 1 |
| 0.093 8 | −0.024 2 | 0.437 5 | −0.112 6 | 0.781 3 | −0.089 0 |
| 0.109 4 | −0.031 9 | 0.453 1 | −0.117 0 | 0.796 9 | −0.089 4 |
| 0.125 0 | −0.030 6 | 0.468 8 | −0.118 5 | 0.812 5 | −0.090 9 |
| 0.140 6 | −0.041 6 | 0.484 4 | −0.120 5 | 0.828 1 | −0.085 7 |
| 0.156 3 | −0.052 3 | 0.500 0 | −0.113 1 | 0.843 8 | −0.085 1 |
| 0.171 9 | −0.057 7 | 0.515 6 | −0.106 8 | 0.859 4 | −0.095 1 |
| 0.187 5 | −0.063 7 | 0.531 3 | −0.106 7 | 0.875 0 | −0.090 9 |
| 0.203 1 | −0.069 0 | 0.546 9 | −0.113 7 | 0.890 6 | −0.089 0 |
| 0.218 8 | −0.070 8 | 0.562 5 | −0.110 5 | 0.906 3 | −0.094 0 |
| 0.234 4 | −0.067 0 | 0.578 1 | −0.110 1 | 0.921 9 | −0.097 6 |
| 0.250 0 | −0.063 0 | 0.593 8 | −0.107 8 | 0.937 5 | −0.100 6 |
| 0.265 6 | −0.074 4 | 0.609 4 | −0.107 8 | 0.953 1 | −0.099 8 |
| 0.281 3 | −0.083 1 | 0.625 0 | −0.106 9 | 0.968 8 | −0.094 1 |
| 0.296 9 | −0.086 5 | 0.640 6 | −0.105 9 | 0.984 4 | −0.093 3 |
| 0.312 5 | −0.088 1 | 0.656 3 | −0.108 5 | 1.000 0 | −0.092 8 |
| 0.328 1 | −0.096 2 | 0.671 9 | −0.110 7 |
Data of with θ = 10 and H = 0.7
TABLE 3
| t | t | t | |||
|---|---|---|---|---|---|
| 0.000 0 | 0.000 0 | 0.343 8 | 0.015 3 | 0.687 5 | 0.001 5 |
| 0.015 6 | −0.004 7 | 0.359 4 | 0.013 5 | 0.703 1 | 0.011 6 |
| 0.031 3 | −0.021 0 | 0.375 0 | 0.000 5 | 0.718 8 | 0.014 8 |
| 0.046 9 | −0.024 1 | 0.390 6 | −0.002 0 | 0.734 4 | 0.002 7 |
| 0.062 5 | −0.029 0 | 0.406 3 | 0.002 5 | 0.750 0 | -0.000 8 |
| 0.078 1 | −0.020 0 | 0.421 9 | 0.002 3 | 0.765 6 | -0.008 6 |
| 0.093 8 | −0.014 3 | 0.437 5 | 0.011 6 | 0.781 3 | -0.006 9 |
| 0.109 4 | −0.012 9 | 0.453 1 | 0.003 8 | 0.796 9 | 0.000 1 |
| 0.125 0 | −0.020 6 | 0.468 8 | −0.007 4 | 0.812 5 | 0.006 0 |
| 0.140 6 | −0.015 7 | 0.484 4 | −0.010 5 | 0.828 1 | 0.016 0 |
| 0.156 3 | 0.004 7 | 0.500 0 | −0.012 4 | 0.843 8 | 0.006 2 |
| 0.171 9 | 0.013 6 | 0.515 6 | −0.009 0 | 0.859 4 | 0.011 1 |
| 0.187 5 | 0.011 0 | 0.531 3 | −0.009 4 | 0.875 0 | 0.009 0 |
| 0.203 1 | 0.006 7 | 0.546 9 | −0.016 0 | 0.890 6 | 0.000 3 |
| 0.218 8 | 0.020 0 | 0.562 5 | −0.011 4 | 0.906 3 | -0.003 2 |
| 0.234 4 | 0.016 2 | 0.578 1 | −0.006 7 | 0.921 9 | 0.010 5 |
| 0.250 0 | 0.002 4 | 0.593 8 | −0.002 8 | 0.937 5 | -0.001 1 |
| 0.265 6 | 0.002 5 | 0.609 4 | 0.002 8 | 0.953 1 | 0.001 0 |
| 0.281 3 | 0.008 2 | 0.625 0 | 0.000 9 | 0.968 8 | 0.005 6 |
| 0.296 9 | 0.007 6 | 0.640 6 | −0.006 2 | 0.984 4 | 0.001 9 |
| 0.312 5 | 0.008 3 | 0.656 3 | −0.015 8 | 1.000 0 | -0.004 6 |
| 0.328 1 | 0.008 6 | 0.671 9 | −0.005 1 |
Data of with θ = 100 and H = 0.7
FIGURE 4

Path of XH with θ = 1 and H = 0.5.
FIGURE 5

Path of XH with θ = 10 and H = 0.5.
FIGURE 6

Path of XH with θ = 100 and H = 0.5.
TABLE 4
| t | t | t | |||
|---|---|---|---|---|---|
| 0.000 0 | 0.000 0 | 0.343 8 | 0.939 3 | 0.687 5 | 1.188 3 |
| 0.015 6 | −0.176 1 | 0.359 4 | 0.991 3 | 0.703 1 | 0.992 1 |
| 0.031 3 | 0.009 9 | 0.375 0 | 1.036 3 | 0.718 8 | 0.956 4 |
| 0.046 9 | −0.040 0 | 0.390 6 | 1.218 0 | 0.734 4 | 0.994 3 |
| 0.062 5 | 0.019 0 | 0.406 3 | 1.204 2 | 0.750 0 | 0.885 2 |
| 0.078 1 | 0.088 3 | 0.421 9 | 1.122 9 | 0.765 6 | 0.861 1 |
| 0.093 8 | 0.020 0 | 0.437 5 | 1.111 0 | 0.781 3 | 0.688 6 |
| 0.109 4 | 0.274 4 | 0.453 1 | 1.021 1 | 0.796 9 | 0.653 8 |
| 0.125 0 | 0.231 7 | 0.468 8 | 1.066 0 | 0.812 5 | 0.731 2 |
| 0.140 6 | 0.246 1 | 0.484 4 | 1.007 0 | 0.828 1 | 0.750 8 |
| 0.156 3 | 0.200 4 | 0.500 0 | 1.099 5 | 0.843 8 | 0.866 3 |
| 0.171 9 | 0.172 3 | 0.515 6 | 1.149 7 | 0.859 4 | 0.746 9 |
| 0.187 5 | 0.233 2 | 0.531 3 | 1.162 0 | 0.875 0 | 0.608 0 |
| 0.203 1 | 0.485 9 | 0.546 9 | 1.222 9 | 0.890 6 | 0.618 4 |
| 0.218 8 | 0.697 4 | 0.562 5 | 1.435 0 | 0.906 3 | 0.655 0 |
| 0.234 4 | 0.684 8 | 0.578 1 | 1.447 4 | 0.921 9 | 0.632 1 |
| 0.250 0 | 0.627 5 | 0.593 8 | 1.453 5 | 0.937 5 | 0.610 1 |
| 0.265 6 | 0.777 4 | 0.609 4 | 1.479 4 | 0.953 1 | 0.623 8 |
| 0.281 3 | 0.825 0 | 0.625 0 | 1.276 4 | 0.968 8 | 0.406 6 |
| 0.296 9 | 0.775 4 | 0.640 6 | 1.281 4 | 0.984 4 | 0.389 3 |
| 0.312 5 | 0.878 3 | 0.656 3 | 1.284 8 | 1.000 0 | 0.234 5 |
| 0.328 1 | 0.838 0 | 0.671 9 | 1.168 9 |
Data of with θ = 1 and H = 0.5
TABLE 5
| t | t | t | |||
|---|---|---|---|---|---|
| 0.000 0 | 0.000 0 | 0.343 8 | −0.024 7 | 0.687 5 | −0.492 7 |
| 0.015 6 | −0.054 8 | 0.359 4 | −0.366 6 | 0.703 1 | −0.589 4 |
| 0.031 3 | 0.122 7 | 0.375 0 | −0.452 2 | 0.718 8 | −0.689 0 |
| 0.046 9 | 0.167 9 | 0.390 6 | −0.690 7 | 0.734 4 | −0.507 9 |
| 0.062 5 | 0.151 5 | 0.406 3 | −0.915 4 | 0.750 0 | −0.370 3 |
| 0.078 1 | −0.217 7 | 0.421 9 | −0.954 1 | 0.765 6 | −0.283 2 |
| 0.093 8 | 0.041 1 | 0.437 5 | −1.020 5 | 0.781 3 | −0.445 5 |
| 0.109 4 | −0.061 7 | 0.453 1 | −0.906 9 | 0.796 9 | −0.551 5 |
| 0.125 0 | −0.069 7 | 0.468 8 | −0.855 3 | 0.812 5 | −0.579 9 |
| 0.140 6 | −0.359 2 | 0.484 4 | −0.820 1 | 0.828 1 | −0.509 3 |
| 0.156 3 | −0.348 9 | 0.500 0 | −0.735 7 | 0.843 8 | −0.556 1 |
| 0.171 9 | −0.481 8 | 0.515 6 | −0.822 0 | 0.859 4 | −0.589 2 |
| 0.187 5 | −0.296 6 | 0.531 3 | −0.785 2 | 0.875 0 | −0.501 7 |
| 0.203 1 | −0.471 7 | 0.546 9 | −0.814 6 | 0.890 6 | −0.458 0 |
| 0.218 8 | −0.417 5 | 0.562 5 | −0.823 9 | 0.906 3 | −0.689 5 |
| 0.234 4 | −0.169 3 | 0.578 1 | −0.833 7 | 0.921 9 | −0.784 6 |
| 0.250 0 | −0.126 5 | 0.593 8 | −0.735 3 | 0.937 5 | −0.825 7 |
| 0.265 6 | −0.017 8 | 0.609 4 | −0.539 7 | 0.953 1 | −0.903 4 |
| 0.281 3 | −0.053 6 | 0.625 0 | −0.515 2 | 0.968 8 | −0.736 4 |
| 0.296 9 | −0.071 4 | 0.640 6 | −0.524 5 | 0.984 4 | −0.669 2 |
| 0.312 5 | −0.115 8 | 0.656 3 | −0.489 9 | 1.000 0 | −0.506 1 |
| 0.328 1 | −0.132 2 | 0.671 9 | −0.525 8 |
Data of with θ = 10 and H = 0.5
TABLE 6
| t | t | t | |||
|---|---|---|---|---|---|
| 0.000 0 | 0.000 0 | 0.343 8 | −0.207 4 | 0.687 5 | -0.149 3 |
| 0.015 6 | −0.012 9 | 0.359 4 | −0.373 2 | 0.703 1 | −0.230 8 |
| 0.031 3 | −0.134 8 | 0.375 0 | −0.464 9 | 0.718 8 | 0.164 4 |
| 0.046 9 | 0.069 7 | 0.390 6 | −0.292 5 | 0.734 4 | −0.050 0 |
| 0.062 5 | 0.111 5 | 0.406 3 | −0.244 5 | 0.750 0 | −0.131 7 |
| 0.078 1 | 0.002 9 | 0.421 9 | −0.246 7 | 0.765 6 | −0.218 2 |
| 0.093 8 | −0.058 9 | 0.437 5 | 0.062 8 | 0.781 3 | −0.313 7 |
| 0.109 4 | −0.288 8 | 0.453 1 | −0.091 7 | 0.796 9 | −0.069 1 |
| 0.125 0 | −0.195 6 | 0.468 8 | −0.307 2 | 0.812 5 | −0.239 1 |
| 0.140 6 | −0.046 9 | 0.484 4 | −0.216 2 | 0.828 1 | −0.306 2 |
| 0.156 3 | −0.139 1 | 0.500 0 | −0.241 8 | 0.843 8 | −0.147 8 |
| 0.171 9 | −0.183 3 | 0.515 6 | −0.159 3 | 0.859 4 | −0.203 4 |
| 0.187 5 | −0.117 5 | 0.531 3 | −0.250 9 | 0.875 0 | −0.219 3 |
| 0.203 1 | −0.261 6 | 0.546 9 | −0.344 2 | 0.890 6 | −0.376 9 |
| 0.218 8 | −0.156 8 | 0.562 5 | −0.129 5 | 0.906 3 | 0.051 5 |
| 0.234 4 | −0.221 5 | 0.578 1 | −0.113 0 | 0.921 9 | −0.107 6 |
| 0.250 0 | −0.173 6 | 0.593 8 | −0.191 5 | 0.937 5 | −0.117 3 |
| 0.265 6 | −0.198 5 | 0.609 4 | −0.131 3 | 0.953 1 | −0.274 6 |
| 0.281 3 | 0.067 4 | 0.625 0 | −0.175 8 | 0.968 8 | −0.155 6 |
| 0.296 9 | −0.163 3 | 0.640 6 | −0.100 8 | 0.984 4 | −0.223 2 |
| 0.312 5 | −0.121 9 | 0.656 3 | −0.104 9 | 1.000 0 | −0.232 0 |
| 0.328 1 | −0.161 0 | 0.671 9 | −0.270 3 |
Data of with θ = 100 and H = 0.5
Remark 2From the following numerical results, we can find that it is important to study the estimates of parameters θ and ν.
Statements
Data availability statement
The original contributions presented in the study are included in the article/supplementary material; further inquiries can be directed to the corresponding authors.
Author contributions
All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1.
AlósEMazetONualartDStochastic Calculus with Respect to Gaussian Processes. Ann Prob (2001) 29:766–801. 10.1214/aop/1008956692
2.
BenaïmMCiotirIGauthierC-ESelf-repelling Diffusions via an Infinite Dimensional Approach. Stoch Pde: Anal Comp (2015) 3:506–30. 10.1007/s40072-015-0059-5
3.
BenaïmMLedouxMRaimondOSelf-interacting Diffusions. Probab Theor Relat Fields (2002) 122:1–41. 10.1007/s004400100161
4.
BojdeckiTGorostizaLGTalarczykASome Extensions of Fractional Brownian Motion and Sub-fractional Brownian Motion Related to Particle Systems. Elect Comm Probab (2007) 12:161–72. 10.1214/ecp.v12-1272
5.
BojdeckiTGorostizaLGTalarczykASub-fractional Brownian Motion and its Relation to Occupation Times. Stat Probab Lett (2004) 69:405–19. 10.1016/j.spl.2004.06.035
6.
BojdeckiTGorostizaLGTalarczykAOccupation Time Limits of Inhomogeneous Poisson Systems of Independent Particles. Stochastic Process their Appl (2008) 118:28–52. 10.1016/j.spa.2007.03.008
7.
BojdeckiTGorostizaLGTalarczykASelf-Similar Stable Processes Arising from High-Density Limits of Occupation Times of Particle Systems. Potential Anal (2008) 28:71–103. 10.1007/s11118-007-9067-z
8.
CranstonMLe JanYSelf Attracting Diffusions: Two Case Studies. Math Ann (1995) 303:87–93. 10.1007/bf01460980
9.
CranstonMMountfordTSThe strong Law of Large Numbers for a Brownian Polymer. Ann Probab (1996) 24:1300–23. 10.1214/aop/1065725183
10.
DurrettRTRogersLCGAsymptotic Behavior of Brownian Polymers. Probab Th Rel Fields (1992) 92:337–49. 10.1007/bf01300560
11.
GanYYanLLeast Squares Estimation for the Linear Self-Repelling Diffusion Driven by Fractional Brownian Motion (In Chinese). Sci CHINA Math (2018) 48:1143. 10.1360/scm-2017-0387
12.
GaoHGuoRJinYYanL, Large Time Behavior on the Linear Self-Interacting Diffusion Driven by Sub-fractional Brownian Motion I: Self-Repelling Case. Front Phys (2021) volume 2021. 10.3389/fphy.2021.795210
13.
GauthierC-ESelf Attracting Diffusions on a Sphere and Application to a Periodic Case. Electron Commun Probab (2016) 21(No. 53):1–12. 10.1214/16-ecp4547
14.
HerrmannSRoynetteBBoundedness and Convergence of Some Self-Attracting Diffusions. Mathematische Annalen (2003) 325:81–96. 10.1007/s00208-002-0370-0
15.
HerrmannSScheutzowMRate of Convergence of Some Self-Attracting Diffusions. Stochastic Process their Appl (2004) 111:41–55. 10.1016/j.spa.2003.10.012
16.
LiMModified Multifractional Gaussian Noise and its Application. Phys Scr (2021) 96(1212):125002. 10.1088/1402-4896/ac1cf6
17.
LiMGeneralized Fractional Gaussian Noise and its Application to Traffic Modeling. Physica A (2021) 579:126138. 10.1016/j.physa.2021.126138
18.
LiMMulti-fractional Generalized Cauchy Process and its Application to Teletraffic. Physica A: Stat Mech its Appl (2020) 550(14):123982. 10.1016/j.physa.2019.123982
19.
LiMFractal Time Series-A Tutorial Review. Math Probl Eng (2010) 2010:1–26. 10.1155/2010/157264
20.
MountfordTTarrésPAn Asymptotic Result for Brownian Polymers. Ann Inst H Poincaré Probab Statist (2008) 44:29–46. 10.1214/07-aihp113
21.
NualartDMalliavin Calculus and Related Topics. 2nd ed.New York: Springer (2006).
22.
PemantleRPhase Transition in Reinforced Random Walk and RWRE on Trees. Ann Probab (1988) 16:1229–41. 10.1214/aop/1176991687
23.
ShenGYanLAn Approximation of Subfractional Brownian Motion. Commun Stat - Theor Methods (2014) 43:1873–86. 10.1080/03610926.2013.769598
24.
ShenGYanLEstimators for the Drift of Subfractional Brownian Motion. Commun Stat - Theor Methods (2014) 43:1601–12. 10.1080/03610926.2012.697243
25.
SunXYanLA central Limit Theorem Associated with Sub-fractional Brownian Motion and an Application (In Chinese). Sci Sin Math (2017) 47:1055–76. 10.1360/scm-2016-0748
26.
SunXYanLA Convergence on the Linear Self-Interacting Diffusion Driven by α-stable Motion. Stochastics (2021) 93:1186–208. 10.1080/17442508.2020.1869239
27.
SunXYanLThe Laws of Large Numbers Associated with the Linear Self-Attracting Diffusion Driven by Fractional Brownian Motion and Applications, to Appear in. J Theoret Prob (2021). Online. 10.1007/s10959-021-01126-0
28.
TudorCSome Properties of the Sub-fractional Brownian Motion. Stochastics (2007) 79:431–48. 10.1080/17442500601100331
29.
TudorCInner Product Spaces of Integrands Associated to Subfractional Brownian Motion. Stat Probab Lett (2008) 78:2201–9. 10.1016/j.spl.2008.01.087
30.
TudorCOn the Wiener Integral with Respect to a Sub-fractional Brownian Motion on an Interval. J Math Anal Appl (2009) 351:456–68. 10.1016/j.jmaa.2008.10.041
31.
TudorCSome Aspects of Stochastic Calculus for the Sub-fractional Brownian Motion. Ann Univ Bucuresti Mathematica (2008) 24:199–230.
32.
TudorCA. Analysis of Variations for Self-Similar Processes. Heidelberg, New York: Springer (2013).
33.
YanLHeKChenCThe Generalized Bouleau-Yor Identity for a Sub-fractional Brownian Motion. Sci China Math (2013) 56:2089–116. 10.1007/s11425-013-4604-2
34.
YanLSunYLuYOn the Linear Fractional Self-Attracting Diffusion. J Theor Probab (2008) 21:502–16. 10.1007/s10959-007-0113-y
35.
YanLShenGOn the Collision Local Time of Sub-fractional Brownian Motions. Stat Probab Lett (2010) 80:296–308. 10.1016/j.spl.2009.11.003
36.
YanLShenGHeKItô’s Formula for the Sub-fractional Brownian Motion. Comm Stochastic Anal (2011) 5:135–59. 10.31390/cosa.5.1.09
Summary
Keywords
subfractional Brownian motion, self-attracting diffusion, law of large numbers, Malliavin calculus, asymptotic distribution
Citation
Guo R, Gao H, Jin Y and Yan L (2022) Large Time Behavior on the Linear Self-Interacting Diffusion Driven by Sub-Fractional Brownian Motion II: Self-Attracting Case. Front. Phys. 9:791858. doi: 10.3389/fphy.2021.791858
Received
09 October 2021
Accepted
19 November 2021
Published
25 January 2022
Volume
9 - 2021
Edited by
Ming Li, Zhejiang University, China
Reviewed by
Xiangfeng Yang, Linköping University, Sweden
Yu Sun, Our Lady of the Lake University, United States
Updates
Copyright
© 2022 Guo, Gao, Jin and Yan.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Han Gao, 1061760802@qq.com
This article was submitted to Interdisciplinary Physics, a section of the journal Frontiers in Physics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.