# Stochastic Current of Bifractional Brownian Motion

Authors: Jingjun Guo
Journal: Journal of Applied Mathematics
 (2014)

---

## Abstract

We study the regularity of stochastic current defined as Skorohod integral with respect to bifractional Brownian motion through Malliavin calculus. Moreover, we similarly derive some results in the case of multidimensional multiparameter. Finally, we consider stochastic current of bifractional Brownian motion as a distribution in Watanabe spaces.


---

## Body

## 1. Introduction

The fractional Brownian motion was first introduced within a Hilbert space framework by Kolmogorov in [1]. It was further studied by Mandelbrot and Van Ness in [2], who provided a stochastic integral representation of this process in terms of a standard Brownian motion in 1968. In recent years, fractional Brownian motion has become an intense object in stochastic analysis and related fields for the moment, due to its interesting properties, such as self-similarity, and its applications in various scientific areas. However, when Hurst parameters H
≠
1
/
2, fractional Brownian motion is neither a semimartingale nor a Markovian process. The techniques used in Brownian motion cannot be directly applied.Nevertheless, every fractional Brownian motion has its limits in modelling certain phenomena. In order to fit better in concrete situations, several authors have recently introduced some generalized fractional Brownian motions. For instance, we mention subfractional Brownian motion (see [3, 4]) and bifractional Brownian motion (see [5, 6]).The concept of current comes from geometric measure theory. The simplest is the functional(1)
φ
⟶
∫
0
T
‍
〈
φ
(
γ
(
t
)
)
,
γ
(
t
)
′
〉
R
d
d
t
,
where φ
:
R
d
→
R
d and γ
(
t
) is a rectifiable curve. This functional ξ
(
x
) can be defined by
(2)
ξ
(
x
)
=
∫
0
T
‍
δ
(
x
-
γ
(
t
)
)
γ
(
t
)
′
d
t
,
where δ
(
x
) is a Dirac function (see [7]). If we want to simulate this current, we need to replace the deterministic curve γ
(
t
) with stochastic process X
t. At the same time, the stochastic integral must be properly interpreted. Recently, people pay attentions to the research on stochastic current. Give the following map:
(3)
φ
⟶
I
(
φ
)
=
∫
0
T
‍
〈
φ
(
X
t
)
,
d
X
t
〉
,
where φ is a vector function on R
d which belongs to some Banach spaces V, X
t is a stochastic process, and the integral is some version of a stochastic integral defined through regularization. Stochastic current is a continuous version of the mapping; that is, stochastic current is regarded as a stochastic element of the dual space of V in [8].The problem of stochastic current is motivated by the study of fluidodynamical models. In [9], in the study of the energy of a vortex filament naturally appear some stochastic double integrals related to Wiener process
(4)
∫
[
0
,
T
]
2
‍
f
(
X
s
-
X
t
)
d
X
s
d
X
t
,
where f
(
x
)
=
K
α
(
x
) is the kernel of the pseudodifferential operator (
1
-
Δ
)
-
α. In the recent years, some results of stochastic currents of Gaussian processes have been obtained through different stochastic integrals in [7, 8, 10]. For example, Flandoli and Tudor [7] have studied the existence and regularity of stochastic currents through Malliavin calculus, where the integrals are defined as Skorohod integrals with respect to the Brownian motion and fractional Brownian motion, respectively. In [10] authors have shown the Sobolev regularity of the stochastic current, which is associated with the pathwise integral.Recall that the bifractional Brownian motionB
H
,
K is a centered Gaussian process with covariance function
(5)
R
H
,
K
(
t
,
s
)
=
1
2
K
(
(
t
2
H
+
s
2
H
)
K
-
|
t
-
s
|
2
H
K
)
,
where parameters H
∈
(
0,1
) and K
∈
(
0,1
]. It is well known that, when K
=
1, bifractional Brownian motion is a fractional Brownian motion. Since bifractional Brownian motion seems to be more flexible and more complex model than fractional Brownian motion, it seems desirable to extend the stochastic current of fractional Brownian motion to the case of bifractional Brownian motion. For this aim, motivated by [7, 11], we use Malliavin calculus and multiple integrals to discuss the stochastic current defined as divergence integral with respect to bifractional Brownian motion. Let us compare our results with the analogous ones from the case of fractional Brownian stochastic current. Note that the regularity condition of bifractional Brownian current does not depend on parameters H and K, while the situation is different in the case of fractional Brownian motion. On the other hand, because the problems of bifractional Brownian motion are more complex, we need some useful techniques to deal with bifractional Brownian current.The paper is organized as follows. In Section2, we provide some background materials from bifractional Brownian motion. In Section 3, we firstly consider the regularity of stochastic current of bifractional Brownian motion with respect to x. Lastly, we regard stochastic current of bifractional Brownian motion as a distribution in Watanabe spaces.


## 2. Bifractional Brownian Motion

In this section, we briefly recall some notations and facts of bifractional Brownian motion, and for details see [5, 6, 11].A bifractional Brownian motionB
H
,
K is a center Gaussian process with variance
(6)
R
H
,
K
(
t
,
s
)
=
1
2
K
(
(
t
2
H
+
s
2
H
)
K
-
|
t
-
s
|
2
H
K
)
,
where parameters H
∈
(
0,1
) and K
∈
(
0,1
]. In the case K
=
1 we retrieve the fractional Brownian motion, while in the case K
=
1 and H
=
1
/
2 bifractional Brownian motion corresponds to the Brownian motion.Letℋ be a Hilbert space. ℋ is defined as the completion of the linear space generated by the I
[
0
,
t
]
,
t
∈
[
0
,
T
] with respect to the inner product
(7)
〈
I
[
0
,
t
]
,
I
[
0
,
s
]
〉
ℋ
=
R
H
,
K
(
t
,
s
)
=
∬
0
T
I
[
0
,
t
]
(
u
)
I
[
0
,
s
]
(
v
)
∂
2
R
H
,
K
(
u
,
v
)
∂
u
∂
v
d
u
d
v
.Sometimes working with the spaceℋ is not convenient, because this space also contains distributions (see [11]) and the norm in this space is not always tractable. We always use the subspace ℋ
1 of ℋ, which is defined as the set of measurable functions f on [
0
,
T
] with
(8)
∥
f
∥
ℋ
1
2
≡
∬
0
T
|
f
(
u
)
|
|
f
(
v
)
|
|
∂
2
R
H
,
K
(
u
,
v
)
∂
u
∂
v
|
d
u
d
v
<
∞
.
We can prove that ℋ
1 is a Banach space for the norm ∥
·
∥
ℋ
1. At the same time, we have
(9)
L
2
(
[
0
,
T
]
)
⊂
ℋ
1
⊂
ℋ
.Denote multiple stochastic integrals byI
n
(
f
n
) with respect to B
H
,
K, where f
n
∈
ℋ
⊗
n. For each F
∈
L
2
(
[
0
,
T
]
), F has chaos expansion F
=
∑
n
=
0
∞
‍
I
n
(
f
n
). Let L be Ornstein-Uhlenbeck operator
(10)
L
F
=
-
∑
n
=
0
∞
‍
n
I
n
(
f
n
)
.
For each p
∈
(
1
,
∞
) and α
∈
ℝ, define Sobolev-Watanabe D
α
,
p as the closure of the set of polynomial variables with respect to the norm
(11)
∥
F
∥
α
,
p
=
∥
(
I
d
-
L
)
α
/
2
∥
L
p
(
Ω
)
,
where I denotes the identity. Malliavin derivative operator D is defined as follows:
(12)
D
t
(
I
n
(
f
n
)
)
=
n
I
n
-
1
(
f
n
(
·
,
t
)
)
.
It is well known that stochastic variable F belongs to D
α
,
2 if and only if
(13)
∑
n
=
0
∞
‍
(
1
+
n
)
α
∥
I
n
(
f
n
)
∥
L
2
(
Ω
)
2
<
∞
.The adjoint ofD is always called the divergence integral (or Skorohod integral). For adapted integrands, the divergence integral coincides with the classical Ito
^ integral. Hence the divergence integral is called generalized It
o
^ integral. If u is a stochastic process, it has the following chaos expansion:
(14)
u
s
=
∑
n
≥
0
‍
I
n
(
f
n
(
·
,
s
)
)
,
where f
n
(
·
,
s
)
∈
ℋ
⊗
(
n
+
1
). Skorohod integral of u is defined as
(15)
∫
0
T
u
s
d
B
s
H
,
K
=
∑
n
≥
0
‍
I
n
+
1
(
f
n
(
·
,
s
)
(
s
)
)
,
where f
n
(
s
) denotes the symmetrization of f
n with respect to n
+
1 variables.


## 3. Stochastic Current of Bifractional Brownian Motion

### 3.1. Stochastic Current of One-Dimensional Case with respect tox

In this section, we give stochastic current of bifractional Brownian motion as follows:(16)
ξ
(
x
)
=
∫
[
0
,
T
]
N
‍
δ
(
x
-
B
s
H
,
K
)
d
B
s
H
,
K
,
where the integral is a Skorohod integral, x
∈
ℝ, and T
>
0.Put(17)
β
n
x
(
s
)
=
p
s
2
H
K
(
x
)
[
s
2
H
K
]
n
/
2
H
n
(
x
s
H
K
)
=
(
R
H
,
K
(
s
)
)
-
(
n
/
2
)
p
R
H
,
K
(
s
)
(
x
)
H
n
(
x
R
H
,
K
(
s
)
1
/
2
)
,
where p
s
2
H
K
(
x
) is a Gaussian kernel function of variance s
2
H
K and H
n
(
x
) is the Hermite polynomial of degree n.By Lemma 3.1 in [7], the following lemma is obtained. Indeed, the lemma can be regarded as a version in the case of bifractional Brownian motion.Lemma 1.
Useβ
n
x
^
(
s
) to denote the Fourier transform of the function x
→
β
n
x
(
s
); then
(18)
β
n
x
^
(
s
)
=
exp
⁡
{
-
x
2
2
R
H
,
K
(
s
)
}
(
-
i
)
n
x
n
n
!
.
Applying Lemma1 and as in [7], we can obtain the stochastic current of bifractional Brownian motion.Theorem 2.
LetB
H
,
K be a bifractional Brownian motion with Hurst parameters H
∈
(
0,1
), K
∈
(
0,1
] satisfying 2
H
K
>
1 and let ξ
(
x
) be given by (16). Then, for each ω
∈
Ω and when r
>
1
/
2
,
ξ
(
x
) belongs to the negative Sobolev space H
-
r
(
ℝ
;
ℝ
).Proof.
By the chaos expansion ofδ
(
x
-
B
H
,
K
) (see [11] or [5]), we have
(19)
δ
(
x
-
B
H
,
K
)
=
∑
n
≥
0
‍
p
s
2
H
K
(
x
)
s
n
H
K
H
n
(
x
s
2
H
K
)
I
n
B
H
,
K
(
I
[
0
,
s
]
⊗
n
(
·
)
)
=
∑
n
≥
0
‍
β
n
x
(
s
)
I
n
B
H
,
K
(
I
[
0
,
s
]
⊗
n
(
·
)
)
,
where I
n
B
H
,
K denotes multiple stochastic integrals with respect to bifractional Brownian motion B
H
,
K.
Let us consider the Fourier transform ofδ
(
x
-
B
H
,
K
):
(20)
δ
^
(
x
-
B
H
,
K
)
=
∑
n
≥
0
‍
β
n
x
^
(
s
)
I
n
B
H
,
K
(
I
[
0
,
s
]
⊗
n
(
·
)
)
.
By Lemma 1, we obtain
(21)
β
n
x
^
(
s
)
=
exp
⁡
{
-
x
2
2
s
2
H
K
}
(
-
i
x
)
n
n
!
.
Hence
(22)
ξ
^
(
x
)
=
∑
n
≥
0
‍
(
-
i
x
)
n
n
!
I
n
+
1
B
H
,
K
(
(
exp
⁡
{
-
x
2
2
s
2
H
K
}
I
[
0
,
s
]
⊗
n
(
·
)
)
(
s
)
)
,
where (
s
) denotes the symmetrization with respect to n
+
1 variables.
By the definition of∥
·
∥
H
-
r
(
ℝ
;
ℝ
) and taking advantage of (22), we get
(23)
E
∥
ξ
(
x
)
∥
H
-
r
(
ℝ
;
ℝ
)
2
=
E
[
∫
ℝ
‍
(
1
+
x
2
)
-
r
|
ξ
^
(
x
)
|
2
d
x
]
≤
∫
ℝ
‍
(
1
+
x
2
)
-
r
∑
n
≥
0
‍
x
2
n
(
n
!
)
2
(
n
+
1
)
!
×
∥
(
exp
⁡
{
-
x
2
2
s
2
H
K
}
I
[
0
,
s
]
⊗
n
(
·
)
)
(
s
)
∥
ℋ
⊗
(
n
+
1
)
2
d
x
.
From the following fact:
(24)
(
exp
⁡
{
-
x
2
2
s
2
H
K
}
I
[
0
,
s
]
⊗
n
(
·
)
)
(
s
)
(
t
1
,
…
,
t
n
+
1
)
=
1
n
+
1
∑
i
=
1
n
+
1
‍
exp
⁡
{
-
x
2
2
t
i
2
H
K
}
I
[
0
,
t
i
]
⊗
n
(
t
1
,
…
,
t
^
i
,
…
,
t
n
+
1
)
,
we show that
(25)
E
∥
ξ
(
x
)
∥
H
-
r
(
ℝ
;
ℝ
)
2
≤
∫
ℝ
‍
(
1
+
x
2
)
-
r
∑
n
≥
0
‍
x
2
n
(
n
!
)
2
(
n
+
1
)
!
×
∬
[
0
,
T
]
n
+
1
∑
i
,
j
=
1
n
‍
1
(
n
+
1
)
2
exp
⁡
{
-
x
2
2
s
i
2
H
K
}
·
exp
⁡
{
-
x
2
2
t
j
2
H
K
}
I
[
0
,
s
i
]
⊗
n
(
s
1
,
…
,
s
^
i
,
…
,
s
n
+
1
)
×
I
[
0
,
t
j
]
⊗
n
(
t
1
,
…
,
t
^
j
,
…
,
t
n
+
1
)
·
∏
q
=
1
n
+
1
‍
∂
2
∂
s
q
∂
t
q
R
H
,
K
×
(
s
q
,
t
q
)
d
s
1
⋯
d
s
n
+
1
d
t
1
⋯
d
t
n
+
1
d
x
=
∫
ℝ
‍
(
1
+
x
2
)
-
r
∑
n
≥
0
‍
x
2
n
(
n
+
1
)
!
×
∑
i
=
1
n
+
1
‍
∬
[
0
,
T
]
n
+
1
exp
⁡
{
-
x
2
2
s
i
2
H
K
}
·
exp
⁡
{
-
x
2
2
t
i
2
H
K
}
I
[
0
,
s
i
]
⊗
n
(
s
1
,
…
,
s
^
i
,
…
,
s
n
+
1
)
×
I
[
0
,
t
i
]
⊗
n
(
t
1
,
…
,
t
^
i
,
…
,
t
n
+
1
)
·
∏
q
=
1
n
+
1
‍
∂
2
∂
s
q
∂
t
q
R
H
,
K
(
s
q
,
t
q
)
d
s
1
⋯
d
s
n
+
1
d
t
1
⋯
d
t
n
+
1
d
x
+
∫
ℝ
‍
(
1
+
x
2
)
-
r
∑
n
≥
0
‍
x
2
n
(
n
+
1
)
!
×
∑
i
,
j
=
1
,
i
≠
j
n
+
1
‍
∬
[
0
,
T
]
n
+
1
exp
⁡
{
-
x
2
2
s
i
2
H
K
}
·
exp
⁡
⁡
{
-
x
2
2
t
j
2
H
K
}
×
I
[
0
,
s
i
]
⊗
n
(
s
1
,
…
,
s
^
i
,
…
,
s
n
+
1
)
×
I
[
0
,
t
j
]
⊗
n
(
t
1
,
…
,
t
^
j
,
…
,
t
n
+
1
)
·
∏
q
=
1
n
+
1
‍
∂
2
∂
s
q
∂
t
q
×
R
H
,
K
(
s
q
,
t
q
)
d
s
1
⋯
d
s
n
+
1
d
t
1
⋯
d
t
n
+
1
d
x
≡
Δ
1
+
Δ
2
.
Firstly, we turn to estimate Δ
1. Using the similar methods in [7, 11] and the following fact:
(26)
∬
0
T
I
[
0
,
t
]
(
u
)
I
[
0
,
s
]
(
v
)
∂
2
∂
u
∂
v
R
H
,
K
(
u
,
v
)
d
u
d
v
=
R
H
,
K
(
t
,
s
)
,
we find
(27)
Δ
1
=
∫
ℝ
‍
(
1
+
x
2
)
-
r
∑
n
≥
0
‍
x
2
n
(
n
+
1
)
!
×
∑
i
=
1
n
+
1
∬
[
0
,
T
]
n
+
1
exp
⁡
{
-
x
2
2
s
1
2
H
K
}
‍
·
exp
⁡
{
-
x
2
2
t
1
2
H
K
}
I
[
0
,
s
1
]
⊗
n
(
s
2
,
…
,
s
^
i
,
…
,
s
n
+
1
)
×
I
[
0
,
t
1
]
⊗
n
(
t
2
,
…
,
t
^
j
,
…
,
t
n
+
1
)
·
∏
q
=
2
n
+
1
‍
∂
2
∂
s
q
∂
t
q
R
H
,
K
(
s
q
,
t
q
)
∂
2
∂
s
1
∂
t
1
×
R
H
,
K
(
s
1
,
t
1
)
d
s
1
·
·
·
d
s
n
+
1
d
t
1
·
·
·
d
t
n
+
1
d
x
=
∫
ℝ
‍
(
1
+
x
2
)
-
r
×
∑
n
≥
0
‍
x
2
n
n
!
∬
[
0
,
T
]
exp
⁡
⁡
{
-
x
2
2
s
1
2
H
K
}
·
exp
⁡
{
-
x
2
2
t
1
2
H
K
}
R
H
,
K
n
(
s
1
,
t
1
)
∂
2
∂
s
1
∂
t
1
×
R
H
,
K
(
s
1
,
t
1
)
d
s
1
d
t
1
d
x
=
∫
ℝ
‍
(
1
+
x
2
)
-
r
∬
[
0
,
T
]
exp
⁡
{
-
x
2
2
s
1
2
H
K
}
exp
⁡
{
-
x
2
2
t
1
2
H
K
}
·
∑
n
≥
0
‍
x
2
n
n
!
R
H
,
K
n
(
s
1
,
t
1
)
∂
2
∂
s
1
∂
t
1
×
R
H
,
K
(
s
1
,
t
1
)
d
s
1
d
t
1
d
x
=
∫
ℝ
‍
(
1
+
x
2
)
-
r
∬
[
0
,
T
]
exp
⁡
⁡
{
-
x
2
2
s
1
2
H
K
}
exp
⁡
⁡
{
-
x
2
2
t
1
2
H
K
}
·
exp
⁡
{
x
2
R
H
,
K
(
s
1
,
t
1
)
}
×
∂
2
∂
s
1
∂
t
1
R
H
,
K
(
s
1
,
t
1
)
d
s
1
d
t
1
d
x
,
where the last equality is established due to Taylor expansion formula of exponential function.
Since(28)
2
R
H
,
K
(
s
1
,
t
1
)
-
s
1
2
H
K
-
t
1
2
H
K
=
-
E
[
(
B
s
1
H
,
K
-
B
t
1
H
,
K
)
2
]
,
we get
(29)
exp
⁡
{
x
2
2
(
2
R
H
,
K
(
s
1
,
t
1
)
-
s
1
2
H
K
-
t
1
2
H
K
)
}
=
exp
⁡
{
-
x
2
2
E
[
(
B
s
1
H
,
K
-
B
t
1
H
,
K
)
2
]
}
.
By [5], for each s
1
,
t
1
∈
[
0
,
T
], we have
(30)
2
-
K
|
s
1
-
t
1
|
2
H
K
≤
E
[
(
B
s
1
H
,
K
-
B
t
1
H
,
K
)
2
]
≤
2
1
-
K
|
s
1
-
t
1
|
2
H
K
,
which implies that
(31)
E
[
(
B
s
1
H
,
K
-
B
t
1
H
,
K
)
2
]
≥
0
.
Use the change of variablesy
=
x
{
E
[
(
B
s
1
H
,
K
-
B
t
1
H
,
K
)
2
]
}
1
/
2. Furthermore, we have
(32)
d
x
=
(
s
1
2
H
K
+
t
1
2
H
K
-
1
2
K
-
1
×
[
(
s
1
2
H
+
t
1
2
H
)
K
|
s
1
-
t
1
|
2
H
K
]
s
1
2
H
K
+
t
1
2
H
K
-
1
2
K
-
1
)
-
1
/
2
d
y
,
(
1
+
x
2
)
-
r
=
{
×
[
(
s
1
2
H
+
t
1
2
H
)
K

-
|
s
1
-
t
1
|
2
H
K
]
1
2
K
-
1
)
-
1

y
2

1
+
(
s
1
2
H
K
+
t
1
2
H
K
-
1
2
K
-
1
×
[
(
s
1
2
H
+
t
1
2
H
)
K

-
|
s
1
-
t
1
|
2
H
K
]
1
2
K
-
1
)
-
1

y
2

}
-
r

.
On the other hand, by [11], there exists a constant C
3,1
(
H
,
K
) depending on H and K such that
(33)
|
∂
2
∂
s
1
∂
t
1
R
H
,
K
(
s
1
,
t
1
)
|
≤
C
3,1
(
H
,
K
)
(
s
1
t
1
)
H
K
-
1
.
Putting (32)-(33) into (27), calculate
(34)
Δ
1
≤
C
3,1
(
H
,
K
)
∬
0
T
(
s
1
t
1
)
H
K
-
1
×
(
s
1
2
H
K
+
t
1
2
H
K
-
1
2
K
-
1
×
[
(
s
1
2
H
+
t
1
2
H
)
K
-
|
s
1
-
t
1
|
2
H
K
]
1
2
K
-
1
)
r
-
1
/
2
×
∫
ℝ
‍
exp
⁡
{
-
y
2
2
}
·
(
s
1
2
H
K
+
t
1
2
H
K
-
1
2
K
-
1
·
[
(
s
1
2
H
+
t
1
2
H
)
K
-
|
s
1
-
t
1
|
2
H
K
]
+
y
2
s
1
2
H
K
+
t
1
2
H
K
-
1
2
K
-
1
)
-
r
d
y
d
s
1
d
t
1
≤
C
3,1
(
H
,
K
)
×
∬
0
T
(
s
1
t
1
)
H
K
-
1
(
s
1
2
H
K
+
t
1
2
H
K
-
1
2
K
-
1
×
[
(
s
1
2
H
+
t
1
2
H
)
K
-
|
s
1
-
t
1
|
2
H
K
]
1
2
K
-
1
)
r
-
1
/
2
×
∫
ℝ
‍
exp
⁡
{
-
y
2
2
}
y
-
2
r
d
y
d
s
1
d
t
1
.
When H
K
-
1
+
2
H
K
(
r
-
1
/
2
)
>
-
1, that is, r
>
0, (34) is finite.
Secondly, using the similar estimation method in the first part, consider the estimation ofΔ
2 as follows:
(35)
Δ
2
=
∫
ℝ
‍
(
1
+
x
2
)
-
r
∑
n
≥
1
‍
x
2
n
(
n
+
1
)
!
n
(
n
+
1
)
×
∬
[
0
,
T
]
n
+
1
exp
⁡
{
-
x
2
2
s
1
2
H
K
}
·
exp
⁡
{
-
x
2
2
t
2
2
H
K
}
I
[
0
,
s
1
]
⊗
n
(
s
2
,
…
,
s
n
+
1
)
×
I
[
0
,
t
2
]
⊗
n
(
t
1
,
t
3
,
…
,
t
n
+
1
)
·
∏
q
=
1
n
+
1
‍
∂
2
∂
s
q
∂
t
q
×
R
H
,
K
(
s
q
,
t
q
)
d
s
1
·
·
·
d
s
n
+
1
d
t
1
·
·
·
d
t
n
+
1
d
x
=
∫
ℝ
‍
(
1
+
x
2
)
-
r
∑
n
≥
1
‍
x
2
n
(
n
-
1
)
!
×
∭
∫
0
T
exp
⁡
⁡
{
-
x
2
2
s
1
2
H
K
}
exp
⁡
⁡
{
-
x
2
2
t
2
2
H
K
}
·
I
[
0
,
s
1
]
(
s
2
)
I
[
0
,
t
2
]
(
t
1
)
∂
2
∂
s
1
∂
t
1
×
R
H
,
K
(
s
1
,
t
1
)
∂
2
∂
s
2
∂
t
2
R
H
,
K
(
s
2
,
t
2
)
·
R
H
,
K
n
-
1
(
s
1
,
t
2
)
d
s
1
d
s
2
d
t
1
d
t
2
d
x
.
By Taylor expansion formula, the following equality is obvious:
(36)
Δ
2
=
∫
ℝ
‍
(
1
+
x
2
)
-
r
x
2
×
∭
∫
0
T
exp
⁡
{
-
x
2
2
(
s
1
2
H
K
+
t
2
2
H
K
-
R
H
,
K
(
s
1
,
t
2
)
)
}
·
I
[
0
,
s
1
]
(
s
2
)
I
[
0
,
t
2
]
(
t
1
)
∂
2
∂
s
1
∂
t
1
R
H
,
K
(
s
1
,
t
1
)
×
∂
2
∂
s
2
∂
t
2
R
H
,
K
(
s
2
,
t
2
)
·
R
H
,
K
n
-
1
(
s
1
,
t
2
)
d
s
1
d
s
2
d
t
1
d
t
2
d
x
.
Recalling some results in [11], there exist parameters H
,
K and a constant C
3,2
(
H
,
K
) depending on H and K such that
(37)
|
∂
2
∂
s
2
∂
t
2
R
H
,
K
(
s
2
,
t
2
)
|
≤
C
3,2
(
H
,
K
)
(
s
2
t
2
)
H
K
-
1
.
Thus, by (33) and (37), we obtain
(38)
Δ
2
≤
C
3,1
(
H
,
K
)
C
3,2
(
H
,
K
)
×
∫
ℝ
‍
(
1
+
x
2
)
-
r
×
∭
∫
0
T
x
2
·
exp
⁡
{
-
x
2
2
[
s
1
2
H
K
+
t
2
2
H
K
-
1
2
K
-
1
×
(
(
s
1
2
H
+
t
2
2
H
)
K
-
|
s
1
-
t
2
|
2
H
K
)
s
1
2
H
K
+
t
2
2
H
K
-
1
2
K
-
1
]
-
x
2
2
[
s
1
2
H
K
+
t
2
2
H
K
-
1
2
K
-
1
}
·
(
s
1
t
1
)
H
K
-
1
(
s
2
t
2
)
H
K
-
1
I
[
0
,
s
1
]
(
s
2
)
×
I
[
0
,
t
2
]
(
t
1
)
d
s
1
d
s
2
d
t
1
d
t
2
d
x
.
Calculate
(39)
∭
∫
0
T
exp
⁡
{
-
x
2
2
[
s
1
2
H
K
+
t
2
2
H
K
-
1
2
K
-
1
×
(
(
s
1
2
H
+
t
2
2
H
)
K
-
|
s
1
-
t
2
|
2
H
K
)
s
1
2
H
K
+
t
2
2
H
K
-
1
2
K
-
1
]
-
x
2
2
}
·
(
s
1
t
1
)
H
K
-
1
(
s
2
t
2
)
H
K
-
1
I
[
0
,
s
1
]
(
s
2
)
×
I
[
0
,
t
2
]
(
t
1
)
d
s
1
d
s
2
d
t
1
d
t
2
=
∬
0
T
exp
⁡
{
-
x
2
2
[
s
1
2
H
K
+
t
2
2
H
K
-
1
2
K
-
1
×
(
(
s
1
2
H
+
t
2
2
H
)
K
-
|
s
1
-
t
2
|
2
H
K
)
s
1
2
H
K
+
t
2
2
H
K
-
1
2
K
-
1
]
-
x
2
2
}
·
(
s
1
t
2
)
H
K
-
1
(
∫
0
T
‍
s
2
H
K
-
1
I
[
0
,
s
1
]
(
s
2
)
d
s
2
)
×
(
∫
0
T
‍
t
1
H
K
-
1
I
[
0
,
t
2
]
(
t
1
)
d
t
1
)
d
s
1
d
t
2
=
1
(
H
K
)
2
×
∬
0
T
exp
⁡
{
-
x
2
2
[
s
1
2
H
K
+
t
2
2
H
K
-
1
2
K
-
1
×
(
(
s
1
2
H
+
t
2
2
H
)
K
-
|
s
1
-
t
2
|
2
H
K
)
s
1
2
H
K
+
t
2
2
H
K
-
1
2
K
-
1
]
-
x
2
2
}
×
(
s
1
t
2
)
2
H
K
-
1
d
s
1
d
t
2
.
Use the change of variables y
=
[
s
1
2
H
K
+
t
2
2
H
K
-
(
1
/
2
K
-
1
)
(
(
s
1
2
H
+
t
2
2
H
)
K
-
|
s
1
-
t
2
|
2
H
K
)
]
1
/
2
x. Hence
(40)
Δ
2
≤
C
3,3
(
H
,
K
)
×
∬
0
T
(
s
1
t
2
)
2
H
K
-
1
(
s
1
2
H
K
+
t
2
2
H
K
-
1
2
K
-
1
×
[
(
s
1
2
H
+
t
2
2
H
)
K
-
|
s
1
-
t
2
|
2
H
K
]
s
1
2
H
K
+
t
2
2
H
K
-
1
2
K
-
1
)
r
-
3
/
2
×
∫
ℝ
‍
exp
⁡
{
-
y
2
2
}
·
(
s
1
2
H
K
+
t
2
2
H
K
-
1
2
K
-
1
×
[
(
s
1
2
H
+
t
2
2
H
)
K
-
|
s
1
-
t
2
|
2
H
K
]
+
y
2
)
-
r
d
y
d
s
1
d
t
2
≤
C
3,3
(
H
,
K
)
∬
0
T
(
s
1
t
2
)
2
H
K
-
1
×
(
s
1
2
H
K
+
t
2
2
H
K
-
1
2
K
-
1
×
[
(
s
1
2
H
+
t
2
2
H
)
K
-
|
s
1
-
t
2
|
2
H
K
]
s
1
2
H
K
+
t
2
2
H
K
-
1
2
K
-
1
)
r
-
3
/
2
×
∫
ℝ
‍
exp
⁡
{
-
y
2
2
}
y
-
2
r
d
y
d
s
1
d
t
2
,
where C
3,3
(
H
,
K
)
=
C
3,1
(
H
,
K
)
C
3,2
(
H
,
K
).
When2
H
K
-
1
+
2
H
K
(
r
-
3
/
2
)
>
-
1, (40) is finite. In other words, when r
>
1
/
2, (40) is finite.
From what we have said above we can draw a conclusion that whenr
>
1
/
2,
(41)
E
∥
ξ
∥
H
r
(
ℝ
;
ℝ
)
2
<
∞
.From Theorem2 we see that when r
>
1
/
2, the mapping ξ
(
x
)
=
∫
[
0
,
T
]
N
‍
δ
(
x
-
B
s
H
,
K
)
d
B
s
H
,
K belongs to negative Sobolev space H
-
r
(
ℝ
;
ℝ
). Note that the regularity condition in Theorem 2 does not depend on H and K. The condition is interesting, because the condition of ξ
(
x
)
=
∫
[
0
,
T
]
N
‍
δ
(
x
-
B
s
)
d
B
s which belongs to negative Sobolev space is also r
>
1
/
2, where B
s is the Brownian motion (see [7]). In other words, they have the same regularity condition. However, the situation is different in the case of fractional Brownian motion, because the regularity condition of fractional Brownian stochastic current is r
>
1
/
2
H
-
1
/
2, which is dependent on Hurst parameter H.


### 3.2. Stochastic Current of d-Dimensional Case with respect tox

As in [7], we can extend stochastic current of one-dimensional bifractional Brownian motion to the case of d-dimensional bifractional Brownian motion.LetB
H
,
K be the vector valued bifractional Brownian motion; that is, B
H
,
K
=
(
B
H
1
,
K
1
,
…
,
B
H
d
,
K
d
), where B
H
i
,
K
i are independent one-dimensional bifractional Brownian sheet. In this part, we consider ξ
(
x
) as follows:
(42)
ξ
(
x
)
=
(
∫
[
0
,
T
]
N
‍
δ
(
x
-
B
s
H
,
K
)
d
B
s
H
1
,
K
1
,
…
,
∫
[
0
,
T
]
N
‍
δ
(
x
-
B
s
H
,
K
)
d
B
s
H
d
,
K
d
)
,
where the integrals are Skorohod integrals with respect to bifractional Brownian motion.Denote(43)
R
H
,
K
l
(
t
,
s
)
=
∏
i
=
1
N
‍
R
H
,
K
l
(
t
i
,
s
i
)
=
〈
I
[
0
,
t
]
,
I
[
0
,
s
]
〉
ℋ
H
l
.Theorem 3.
LetB
H
,
K be d-dimensional bifractional Brownian motion with parameters H
i and K
i satisfying 2
H
i
K
i
>
1
(
i
=
1
,
…
,
d
) and let ξ
(
x
) be given by (42); then, for each ω
∈
Ω and when r
>
d
/
2
-
1, ξ
(
x
) belongs to Sobolev space H
-
r
(
ℝ
;
ℝ
).Proof.
Denoteξ
l
(
x
) by
(44)
ξ
l
(
x
)
=
∫
[
0
,
T
]
N
‍
δ
(
x
-
B
s
H
,
K
)
d
B
s
H
l
,
K
l
=
∫
[
0
,
T
]
N
‍
δ
l
(
x
-
B
s
H
,
K
)
δ
(
x
l
-
B
s
H
l
,
K
l
)
d
B
s
H
l
,
K
l
=
∫
[
0
,
T
]
N
‍
δ
l
(
x
-
B
s
H
,
K
)
×
∑
n
l
≥
0
‍
β
n
l
x
l
(
s
)
I
n
l
B
H
l
,
K
l
(
I
[
0
,
s
]
⊗
n
l
(
·
)
)
d
B
s
H
l
,
K
l
,
where δ
l
(
x
-
B
s
H
,
K
)
=
∏
q
=
1
,
q
≠
l
d
‍
δ
(
x
q
-
B
s
H
q
,
K
q
).
Calculate the Fourier transform of (44) as follows:
(45)
ξ
^
l
(
x
)
=
∑
n
l
≥
0
‍
I
n
l
+
1
B
H
l
,
K
l
(
(
exp
⁡
{
-
i
∑
r
=
1
,
r
≠
l
d
‍
x
r
B
s
H
r
,
K
r
}
×
(
-
i
)
n
l
x
l
n
l
n
l
!
×
exp
⁡
⁡
{
-
x
2
2
|
s
|
2
H
l
K
l
}
I
[
0
,
s
]
⊗
n
l
(
·
)
)
(
s
)
)
.
According to the definition of the normal ∥
ξ
(
x
)
∥
H
-
r
(
ℝ
;
ℝ
)
2 and using Euler formula, we can prove that
(46)
E
|
ξ
^
l
(
x
)
|
2
=
∑
n
l
≥
0
‍
x
l
2
n
l
(
n
l
+
1
)
!
×
∑
i
,
j
=
1
n
l
+
1
‍
∬
[
0
,
T
]
N
(
n
l
+
1
)
∏
q
=
1
n
l
+
1
‍
∂
2
∂
s
∂
t
R
H
,
K
q
(
s
,
t
)
·
cos
⁡
(
∑
r
=
1
,
r
≠
l
d
‍
x
r
B
s
i
H
r
,
K
r
)
cos
⁡
(
∑
r
=
1
,
r
≠
l
d
‍
x
r
B
t
j
H
r
,
K
r
)
×
exp
⁡
{
-
x
l
2
2
|
s
i
|
2
H
l
K
l
}
·
exp
⁡
⁡
{
-
x
l
2
2
|
t
j
|
2
H
l
K
l
}
×
I
[
0
,
s
i
]
⊗
n
l
(
s
1
,
…
,
s
^
i
,
…
,
s
n
l
+
1
)
×
I
[
0
,
t
j
]
⊗
n
l
(
t
1
,
…
,
t
^
j
,
…
,
t
n
l
+
1
)
d
s
d
t
+
∑
n
l
≥
0
‍
x
l
2
n
l
(
n
l
+
1
)
!
×
∑
i
,
j
=
1
n
l
+
1
∬
[
0
,
T
]
N
(
n
l
+
1
)
∏
q
=
1
n
l
+
1
‍
∂
2
∂
s
∂
t
R
H
,
K
q
(
s
,
t
)
‍
·
sin
(
∑
r
=
1
,
r
≠
l
d
‍
x
r
B
s
i
H
r
,
K
r
)
sin
(
∑
r
=
1
,
r
≠
l
d
‍
x
r
B
t
j
H
r
,
K
r
)
×
exp
⁡
{
-
x
l
2
2
|
s
i
|
2
H
l
K
l
}
·
exp
⁡
{
-
x
l
2
2
|
t
j
|
2
H
l
K
l
}
I
[
0
,
s
i
]
⊗
n
l
(
s
1
,
…
,
s
^
i
,
…
,
s
n
l
+
1
)
×
I
[
0
,
t
j
]
⊗
n
l
(
t
1
,
…
,
t
^
j
,
…
,
t
n
l
+
1
)
d
s
d
t
,
where |
s
i
|
2
H
l
K
l
=
∏
j
=
1
N
‍
(
s
j
i
)
2
H
l
K
l
,
I
[
0
,
s
i
]
=
∏
j
=
1
N
‍
I
[
0
,
s
j
i
] and R
H
,
K
q
(
s
,
t
)
=
R
H
q
,
K
q
(
s
q
,
t
q
).
By the bound ofsin
(
x
) and cos
⁡
(
x
), we can get the following inequality:
(47)
E
|
ξ
^
l
(
x
)
|
2
≤
C
3,4
∑
n
l
≥
0
‍
x
l
2
n
l
(
n
l
+
1
)
!
×
∑
i
,
j
=
1
n
l
+
1
‍
∬
[
0
,
T
]
N
(
n
l
+
1
)
∏
q
=
1
n
l
+
1
‍
∂
2
∂
s
∂
t
R
H
,
K
q
(
s
,
t
)
·
exp
⁡
{
-
x
l
2
2
|
s
i
|
2
H
l
K
l
}
×
exp
⁡
{
-
x
l
2
2
|
t
j
|
2
H
l
K
l
}
·
I
[
0
,
s
i
]
⊗
n
l
(
s
1
,
…
,
s
^
i
,
…
,
s
n
l
+
1
)
×
I
[
0
,
t
j
]
⊗
n
l
(
t
1
,
…
,
t
^
j
,
…
,
t
n
l
+
1
)
d
s
d
t
+
C
3,5
∑
n
l
≥
0
‍
x
l
2
n
l
(
n
l
+
1
)
!
×
∑
i
,
j
=
1
n
l
+
1
‍
∬
[
0
,
T
]
N
(
n
l
+
1
)
∏
q
=
1
n
l
+
1
‍
∂
2
∂
s
∂
t
R
H
,
K
q
(
s
,
t
)
·
exp
⁡
{
-
x
l
2
2
|
s
i
|
2
H
l
K
l
}
exp
⁡
{
-
x
l
2
2
|
t
j
|
2
H
l
K
l
}
·
I
[
0
,
s
i
]
⊗
n
l
(
s
1
,
…
,
s
^
i
,
…
,
s
n
l
+
1
)
×
I
[
0
,
t
j
]
⊗
n
l
(
t
1
,
…
,
t
^
j
,
…
,
t
n
l
+
1
)
d
s
d
t
≡
Δ
3
,
l
+
Δ
4
,
l
,
where C
3,4 and C
3,5 are both constants.
By the same estimation techniques ofΔ
2 in Theorem 2, we can obtain the estimation of Δ
4
,
l. Here we need to discuss the estimation of Δ
3
,
l.
Applying∬
0
T
I
[
0
,
t
]
(
u
)
I
[
0
,
s
]
(
v
)
(
∂
2
/
∂
u
∂
v
)
R
H
,
K
(
u
,
v
)
d
u
d
v
=
R
H
,
K
(
t
,
s
) again, we can write
(48)
Δ
3
,
l
=
C
3,4
∑
n
l
≥
0
‍
x
l
2
n
l
n
l
!
×
∬
[
0
,
T
]
N
(
n
l
+
1
)
∏
q
=
1
n
l
+
1
‍
∂
2
∂
s
1
∂
t
1
R
H
,
K
q
(
s
1
,
t
1
)
·
exp
⁡
{
-
x
l
2
2
|
s
1
|
2
H
l
K
l
}
exp
⁡
{
-
x
l
2
2
|
t
1
|
2
H
l
K
l
}
·
I
[
0
,
s
1
]
⊗
n
l
(
s
2
,
…
,
s
n
l
+
1
)
×
I
[
0
,
t
1
]
⊗
n
l
(
t
2
,
…
,
t
n
l
+
1
)
d
s
d
t
=
C
3,4
∑
n
l
≥
0
‍
x
l
2
n
l
n
l
!
×
∬
[
0
,
T
]
N
(
R
H
,
K
l
(
s
1
,
t
1
)
)
n
l
·
∂
2
∂
s
∂
t
R
H
,
K
l
(
s
1
,
t
1
)
×
exp
⁡
{
-
x
l
2
2
|
s
1
|
2
H
l
K
l
}
exp
⁡
{
-
x
l
2
2
|
t
1
|
2
H
l
K
l
}
d
s
1
d
t
1
=
C
3,4
∬
[
0
,
T
]
N
∂
2
∂
s
1
∂
t
1
R
H
,
K
l
(
s
1
,
t
1
)
·
exp
⁡
{
x
l
2
2
(
|
s
1
|
2
H
l
K
l
-
|
t
1
|
2
H
l
K
l
2
R
H
,
K
l
(
s
1
,
t
1
)
-
|
s
1
|
2
H
l
K
l
-
|
t
1
|
2
H
l
K
l
)
x
l
2
2
}
d
s
1
d
t
1
.
According to [11], there exists a constant C
3,6
(
H
,
K
,
l
) depending upon H
,
K, and l such that
(49)
|
∂
2
∂
s
1
∂
t
1
R
H
,
K
l
(
s
1
,
t
1
)
|
≤
C
3,6
(
H
,
K
,
l
)
(
s
1
t
1
)
H
l
K
l
-
1
.
Therefore
(50)
∫
ℝ
d
‍
(
1
+
|
x
|
2
)
-
r
Δ
3
,
l
d
x
≤
C
3,6
(
H
,
K
,
l
)
×
∫
ℝ
d
‍
∬
[
0
,
T
]
N
(
s
1
t
1
)
H
l
K
l
-
1
·
exp
⁡
{
-
x
2
2
E
[
(
B
s
1
H
l
,
K
l
-
B
t
1
H
l
,
K
l
)
2
]
}
d
s
1
d
t
1
d
x
.
Use the change of variables y
l
=
x
l
{
E
[
(
B
s
1
H
l
,
K
l
-
B
t
1
H
l
,
K
l
)
2
]
}
1
/
2. Thus
(51)
∫
ℝ
d
‍
(
1
+
|
x
|
2
)
-
r
Δ
3
,
l
d
x
≤
C
3,7
(
H
,
K
,
l
)
×
∬
[
0
,
T
]
N
(
s
1
t
1
)
H
l
K
l
-
1
{
E
[
(
B
s
1
H
l
,
K
l
-
B
t
1
H
l
,
K
l
)
2
]
}
r
-
d
/
2
·
∫
ℝ
d
‍
exp
⁡
{
-
y
l
2
2
}
×
(
E
[
(
B
s
1
H
l
,
K
l
-
B
t
1
H
l
,
K
l
)
2
]
+
y
l
2
)
-
r
d
s
1
d
t
1
d
y
,
where y
=
(
y
1
,
…
,
y
d
).
On the other hand, by [6] for arbitrary ɛ
≥
0
,
s
1
,
t
1
∈
[
ɛ
,
1
]
N, there exists a constant C
3,8 such that
(52)
E
[
(
B
s
1
H
l
,
K
l
-
B
t
1
H
l
,
K
l
)
2
]
≤
C
3,8
∑
j
=
1
N
‍
|
s
j
1
-
t
j
1
|
2
H
l
K
l
.
In order to be simple, here we only consider the case of T
=
1.
Comparing (51) with (52), we find that
(53)
∫
ℝ
d
‍
(
1
+
|
x
|
2
)
-
r
Δ
3
,
l
d
x
≤
C
3,9
∬
[
0
,
T
]
N
(
s
1
t
1
)
H
l
K
l
-
1
·
(
∑
j
=
1
N
‍
|
s
j
1
-
t
j
1
|
2
H
j
K
j
)
r
-
d
/
2
×
∫
ℝ
d
‍
exp
⁡
{
-
y
l
2
2
}
y
-
2
r
d
y
d
s
d
t
.
When H
l
K
l
-
1
+
2
H
l
K
l
(
r
-
d
/
2
)
>
-
1, that is, r
>
d
/
2
-
1, (53) is finite. That is to say, when r
>
d
/
2
-
1, there is
(54)
∫
ℝ
d
‍
(
1
+
|
x
|
2
)
-
r
E
[
|
ξ
^
(
x
)
|
2
]
d
x
=
∫
ℝ
d
‍
(
1
+
|
x
|
2
)
-
r
E
[
|
ξ
^
1
(
x
)
|
2
+
·
·
·
+
|
ξ
^
d
(
x
)
|
2
]
d
x
<
∞
.It is interesting to contrast Theorems2 and 3 with Propositions 3 and 4 in [7]. Parameters H and K of bifractional Brownian motion do not affect the regularity condition of bifractional Brownian current. However, the regularity condition of stochastic current of one-dimensional fractional Brownian motion is different from the case of d-dimensional setting (see [7]). In other words, Hurst parameters H of fractional Brownian motion have influence on fractional Brownian currents in the case of different dimension.


### 3.3. Stochastic Current of d-Dimensional Bifractional Brownian Motion with respect toω

LetB
H
,
K be vector valued bifractional Brownian motion where vectors H
=
(
H
1
,
…
,
H
d
)
∈
(
0,1
)
d and K
=
(
K
1
,
…
,
K
d
)
∈
(
0,1
]
d; that is, B
H
,
K
=
(
B
H
1
,
K
1
,
…
,
B
H
d
,
K
d
). In this part, ξ
(
x
) is given by
(55)
ξ
(
x
)
=
∫
0
T
‍
δ
(
x
-
B
s
H
,
K
)
d
B
s
H
,
K
,
where
(56)
δ
(
x
-
B
s
H
,
K
)
=
∏
i
=
1
d
‍
δ
(
x
i
-
B
s
H
i
,
K
i
)
=
∑
n
=
(
n
1
,
…
,
n
d
)
‍
β
n
(
s
,
x
)
I
n
B
H
i
,
K
i
(
I
[
0
,
s
]
⊗
|
n
|
(
·
)
)
,
β
n
(
s
,
x
)
=
∏
i
=
1
d
‍
1
R
H
,
K
i
(
s
)
n
i
/
2
p
R
H
,
K
i
(
x
i
)
H
n
i
(
x
i
R
H
,
K
i
(
s
)
1
/
2
)
,
R
H
,
K
i
(
s
)
=
R
H
,
K
i
(
s
,
s
)
=
s
2
H
i
K
i
,
I
n
B
H
,
K
(
I
[
0
,
s
]
⊗
|
n
|
(
·
)
)
=
∏
i
=
1
d
‍
I
n
i
B
H
i
,
K
i
(
I
[
0
,
s
]
⊗
n
i
(
·
)
)
.Using the chaos expansion of divergent integral with respect to bifractional Brownian motion, we can obtain the following expression:(57)
ξ
i
(
x
)
=
∫
0
T
‍
δ
(
x
-
B
s
H
,
K
)
d
B
s
H
i
,
K
i
=
∑
n
=
(
n
1
,
…
,
n
d
)
‍
I
n
i
+
1
B
H
i
,
K
i
[
∏
j
=
1
,
j
≠
i
d
‍
I
n
j
B
H
j
,
K
j
(
I
[
0
,
s
]
⊗
n
j
(
·
)
)
)
(
s
)
(
β
n
(
s
,
x
)
I
[
0
,
s
]
⊗
n
i
(
s
1
,
…
,
s
n
i
)
I
[
0
,
T
]
(
s
)
·
∏
j
=
1
,
j
≠
i
d
‍
I
n
j
B
H
j
,
K
j
(
I
[
0
,
s
]
⊗
n
j
(
·
)
)
)
(
s
)
]
.Theorem 4.
LetB
H
,
K be a bifractional Brownian motion with parameters H
=
(
H
1
,
…
,
H
d
)
∈
(
0,1
)
d and K
=
(
K
1
,
…
,
K
d
)
∈
(
0,1
]
d satisfying 2
H
i
K
i
>
1
(
i
=
1
,
…
,
d
). If ξ
(
x
) is given by (55), then, for every x
∈
ℝ
d and α
<
1
/
2
(
H
K
)
*
-
d
/
2, ξ
(
x
) is one member of Sobolev-Watanabe spaces D
α
-
1,2, where (
H
K
)
*
=
max
⁡
{
H
1
K
1
,
…
,
H
d
K
d
}.Proof.
In order to be convenient, we always replace the normal∥
·
∥
ℋ
⊗
(
n
+
1
) with the normal ∥
·
∥
ℋ
1
⊗
(
n
+
1
). Using the chaos expansion of δ
(
x
-
B
H
,
K
) and by the definition of the normal of ∥
ξ
(
x
)
∥
2
,
α
-
1
2, we verify that
(58)
∥
ξ
i
(
x
)
∥
2
,
α
-
1
2
≤
∑
m
≥
1
‍
(
m
+
1
)
α
-
1
∑
|
n
|
=
n
1
+
·
·
·
+
n
d
=
m
-
1
‍
(
n
i
+
1
)
!
·
∬
[
0
,
T
]
n
i
+
1
∑
l
,
k
=
1
n
i
+
1
‍
β
n
(
s
l
)
β
n
(
t
k
)
I
[
0
,
T
]
(
s
l
)
I
[
0
,
T
]
(
t
k
)
·
I
[
0
,
s
l
]
⊗
n
i
(
s
1
,
…
,
s
^
l
,
…
,
s
n
i
+
1
)
×
I
[
0
,
t
k
]
⊗
n
i
(
t
1
,
…
,
t
^
k
,
…
,
t
n
i
+
1
)
·
∏
j
=
1
,
j
≠
i
d
‍
n
j
!
R
H
,
K
j
(
s
l
,
t
k
)
n
j
×

∏
q
=
1
n
i
+
1
‍

∂
2
∂
s
q
∂
t
q
R
H
,
K
i

(
s
q
,
t
q
)

d
t
1
·
·
·
d
t
n
i
+
1
d
s
1
·
·
·
d
s
n
i
+
1
.
Because α
<
0, the following inequalities: (
m
+
2
)
α
-
1
≤
(
m
+
1
)
α
-
1 and n
i
+
1
≤
m
+
1, are obvious. Similarly to what is performed in [11] (or see [7]), we have
(59)
∥
ξ
i
(
x
)
∥
2
,
α
-
1
2
≤
∑
m
≥
0
‍
(
m
+
1
)
α
∑
|
n
|
=
n
1
+
·
·
·
+
n
d
=
m
‍
n
i
!
×
[
(
1
-
1
n
i
+
1
)
×
∬
[
0
,
T
]
2
β
n
(
s
1
)
β
n
(
t
2
)
·
I
[
0
,
T
]
(
s
1
)
I
[
0
,
T
]
(
t
2
)
×
R
H
,
K
i
(
s
1
,
t
2
)
n
i
-
1
(
1
-
1
n
i
+
1
)
]
×
∏
j
=
1
,
j
≠
i
d
‍
n
j
!
R
H
,
K
j
(
s
1
,
t
2
)
n
j
·
I
[
0
,
s
1
]
(
s
2
)
×
I
[
0
,
t
2
]
(
t
1
)
∂
2
∂
s
1
∂
t
1
R
H
,
K
i
(
s
1
,
t
1
)
×
∂
2
∂
s
2
∂
t
2
R
H
,
K
i
(
s
2
,
t
2
)
d
s
1
d
s
2
d
t
1
d
t
2
+
1
n
i
+
1
×
∬
[
0
,
T
]
2
β
n
(
s
1
)
β
n
(
t
1
)
I
[
0
,
T
]
(
s
1
)
×
I
[
0
,
T
]
(
t
1
)
R
H
,
K
i
(
s
1
,
t
1
)
n
i
-
1
·
∏
j
=
1
,
j
≠
i
d
‍
n
j
!
R
H
,
K
j
(
s
1
,
t
1
)
n
j
I
[
0
,
s
1
]
(
s
2
)
×
I
[
0
,
t
1
]
(
t
2
)
·
∂
2
∂
s
1
∂
t
1
R
H
,
K
i
(
s
1
,
t
1
)
∂
2
∂
s
2
∂
t
2
×
R
H
,
K
i
(
s
2
,
t
2
)
d
s
1
d
s
2
d
t
1
d
t
2
≤
∑
m
≥
0
‍
(
m
+
1
)
α
∑
|
n
|
=
n
1
+
·
·
·
+
n
d
=
m
‍
n
i
!
×
[
(
1
-
1
n
i
+
1
)
×
∫
[
0
,
T
]
2
‍
β
n
(
s
1
)
β
n
(
t
2
)
·
R
H
,
K
i
(
s
1
,
t
2
)
n
i
-
1
×
∏
j
=
1
,
j
≠
i
d
‍
n
j
!
R
H
,
K
j
(
s
1
,
t
2
)
n
j
-
1
·
∂
∂
s
1
R
H
,
K
i
(
s
1
,
t
2
)
×
∂
∂
t
2
R
H
,
K
i
(
s
1
,
t
2
)
d
s
1
d
t
2
+
1
n
i
+
1
×
∫
[
0
,
T
]
2
‍
β
n
(
s
1
)
β
n
(
t
1
)
R
H
,
K
i
(
s
1
,
t
1
)
n
i
-
1
·
∏
j
=
1
,
j
≠
i
d
‍
n
j
!
R
H
,
K
j
(
s
1
,
t
1
)
n
j
×
∂
2
∂
s
1
∂
t
1
R
H
,
K
i
(
s
1
,
t
1
)
d
s
1
d
t
1
]
.
According to (4.37) in [6], for β
∈
[
1
/
4
,
1
/
2
), it follows that
(60)
∏
j
=
1
d
‍
β
n
(
u
)
≤
C
3,10
∏
j
=
1
d
‍
1
n
j
!
(
n
j
∨
1
)
(
(
8
β
-
1
)
/
12
)
,
where C
3,10 is a constant.
On the other hand, as in [11] and using the inequality a
2
+
b
2
≤
2
a
b, for a
,
b
∈
ℝ
+, there exists a constant C
3,11
(
H
,
K
) depending on H and K such that
(61)
|
∂
R
H
,
K
i
(
s
,
t
)
∂
t
∂
R
H
,
K
i
(
s
,
t
)
∂
s
|
≤
C
3,11
(
H
,
K
)
(
s
t
)
2
H
i
K
i
-
1
,
|
∂
2
R
H
,
K
i
(
s
,
t
)
∂
s
∂
t
|
≤
C
3,11
(
H
,
K
)
(
s
t
)
H
i
K
i
-
1
,
|
R
H
,
K
i
(
s
,
t
)
(
s
t
)
H
i
K
i
|
≤
C
3,11
(
H
,
K
)
.
Putting (60) and (61) into (59) and applying the self-similarity of the covariance R
(
s
,
t
), we have
(62)
∥
ξ
i
(
x
)
∥
2
,
α
-
1
2
≤
C
3,12
∑
m
≥
0
‍
(
1
+
m
)
α
×
∑
|
n
|
=
n
1
+
·
·
·
+
n
d
=
m
‍
∏
j
=
1
d
‍
1
(
n
j
∨
1
)
(
8
β
-
1
)
/
6
·
∫
[
0
,
T
]
2
‍
R
H
,
K
i
(
s
,
t
)
n
i
-
1
(
s
t
)
(
n
i
-
1
)
H
i
K
i
(
s
t
)
H
i
K
i
-
1
×
∏
j
=
1
,
j
≠
i
d
‍
R
H
,
K
j
(
s
,
t
)
n
j
(
s
t
)
n
j
H
j
K
j
d
s
d
t
≤
C
3,12
∑
m
≥
0
‍
(
1
+
m
)
α
×
∑
|
n
|
=
n
1
+
·
·
·
+
n
d
=
m
‍
∏
j
=
1
d
‍
1
(
n
j
∨
1
)
(
8
β
-
1
)
/
6
·
∫
0
T
‍
t
2
H
i
K
i
-
1
d
t
∫
0
1
‍
z
H
i
K
i
-
1
(
R
H
i
,
K
i
i
(
1
,
z
)
z
H
i
K
i
)
n
i
-
1
×
∏
j
=
1
,
j
≠
i
d
‍
R
H
,
K
j
(
1
,
z
)
n
j
z
H
j
K
j
d
z
,
where we use the change of variables t
=
t and z
=
s
/
r (see [6, 11]). According to (4.9) in [6], the bounded condition of integral part in (62) is proved.
Indeed,(63)
∫
0
1
‍
∏
j
=
1
d
‍
(
R
H
,
K
j
(
1
,
z
)
z
H
j
K
j
)
n
j
d
z
≤
C
3,13
m
-
1
/
2
(
H
K
)
*
,
where (
H
K
)
*
=
max
⁡
{
H
1
K
1
,
…
,
H
d
K
d
}.
Therefore, whenα
<
1
/
2
(
H
K
)
*
-
d
(
1
-
(
8
β
-
1
)
/
6
),
(64)
∥
ξ
i
(
x
)
∥
2
,
α
-
1
2
≤
C
3,14
∑
m
≥
1
‍
(
1
+
m
)
α
m
-
1
/
2
(
H
K
)
*
-
1
+
d
(
1
-
(
8
β
-
1
)
/
6
)
<
∞
.
Since β
∈
[
1
/
4
,
1
/
2
), we may choose β to tend to 1
/
2. Hence, when α
<
1
/
2
(
H
K
)
*
-
d
/
2, (64) is finite.


## 3.1. Stochastic Current of One-Dimensional Case with respect tox

In this section, we give stochastic current of bifractional Brownian motion as follows:(16)
ξ
(
x
)
=
∫
[
0
,
T
]
N
‍
δ
(
x
-
B
s
H
,
K
)
d
B
s
H
,
K
,
where the integral is a Skorohod integral, x
∈
ℝ, and T
>
0.Put(17)
β
n
x
(
s
)
=
p
s
2
H
K
(
x
)
[
s
2
H
K
]
n
/
2
H
n
(
x
s
H
K
)
=
(
R
H
,
K
(
s
)
)
-
(
n
/
2
)
p
R
H
,
K
(
s
)
(
x
)
H
n
(
x
R
H
,
K
(
s
)
1
/
2
)
,
where p
s
2
H
K
(
x
) is a Gaussian kernel function of variance s
2
H
K and H
n
(
x
) is the Hermite polynomial of degree n.By Lemma 3.1 in [7], the following lemma is obtained. Indeed, the lemma can be regarded as a version in the case of bifractional Brownian motion.Lemma 1.
Useβ
n
x
^
(
s
) to denote the Fourier transform of the function x
→
β
n
x
(
s
); then
(18)
β
n
x
^
(
s
)
=
exp
⁡
{
-
x
2
2
R
H
,
K
(
s
)
}
(
-
i
)
n
x
n
n
!
.
Applying Lemma1 and as in [7], we can obtain the stochastic current of bifractional Brownian motion.Theorem 2.
LetB
H
,
K be a bifractional Brownian motion with Hurst parameters H
∈
(
0,1
), K
∈
(
0,1
] satisfying 2
H
K
>
1 and let ξ
(
x
) be given by (16). Then, for each ω
∈
Ω and when r
>
1
/
2
,
ξ
(
x
) belongs to the negative Sobolev space H
-
r
(
ℝ
;
ℝ
).Proof.
By the chaos expansion ofδ
(
x
-
B
H
,
K
) (see [11] or [5]), we have
(19)
δ
(
x
-
B
H
,
K
)
=
∑
n
≥
0
‍
p
s
2
H
K
(
x
)
s
n
H
K
H
n
(
x
s
2
H
K
)
I
n
B
H
,
K
(
I
[
0
,
s
]
⊗
n
(
·
)
)
=
∑
n
≥
0
‍
β
n
x
(
s
)
I
n
B
H
,
K
(
I
[
0
,
s
]
⊗
n
(
·
)
)
,
where I
n
B
H
,
K denotes multiple stochastic integrals with respect to bifractional Brownian motion B
H
,
K.
Let us consider the Fourier transform ofδ
(
x
-
B
H
,
K
):
(20)
δ
^
(
x
-
B
H
,
K
)
=
∑
n
≥
0
‍
β
n
x
^
(
s
)
I
n
B
H
,
K
(
I
[
0
,
s
]
⊗
n
(
·
)
)
.
By Lemma 1, we obtain
(21)
β
n
x
^
(
s
)
=
exp
⁡
{
-
x
2
2
s
2
H
K
}
(
-
i
x
)
n
n
!
.
Hence
(22)
ξ
^
(
x
)
=
∑
n
≥
0
‍
(
-
i
x
)
n
n
!
I
n
+
1
B
H
,
K
(
(
exp
⁡
{
-
x
2
2
s
2
H
K
}
I
[
0
,
s
]
⊗
n
(
·
)
)
(
s
)
)
,
where (
s
) denotes the symmetrization with respect to n
+
1 variables.
By the definition of∥
·
∥
H
-
r
(
ℝ
;
ℝ
) and taking advantage of (22), we get
(23)
E
∥
ξ
(
x
)
∥
H
-
r
(
ℝ
;
ℝ
)
2
=
E
[
∫
ℝ
‍
(
1
+
x
2
)
-
r
|
ξ
^
(
x
)
|
2
d
x
]
≤
∫
ℝ
‍
(
1
+
x
2
)
-
r
∑
n
≥
0
‍
x
2
n
(
n
!
)
2
(
n
+
1
)
!
×
∥
(
exp
⁡
{
-
x
2
2
s
2
H
K
}
I
[
0
,
s
]
⊗
n
(
·
)
)
(
s
)
∥
ℋ
⊗
(
n
+
1
)
2
d
x
.
From the following fact:
(24)
(
exp
⁡
{
-
x
2
2
s
2
H
K
}
I
[
0
,
s
]
⊗
n
(
·
)
)
(
s
)
(
t
1
,
…
,
t
n
+
1
)
=
1
n
+
1
∑
i
=
1
n
+
1
‍
exp
⁡
{
-
x
2
2
t
i
2
H
K
}
I
[
0
,
t
i
]
⊗
n
(
t
1
,
…
,
t
^
i
,
…
,
t
n
+
1
)
,
we show that
(25)
E
∥
ξ
(
x
)
∥
H
-
r
(
ℝ
;
ℝ
)
2
≤
∫
ℝ
‍
(
1
+
x
2
)
-
r
∑
n
≥
0
‍
x
2
n
(
n
!
)
2
(
n
+
1
)
!
×
∬
[
0
,
T
]
n
+
1
∑
i
,
j
=
1
n
‍
1
(
n
+
1
)
2
exp
⁡
{
-
x
2
2
s
i
2
H
K
}
·
exp
⁡
{
-
x
2
2
t
j
2
H
K
}
I
[
0
,
s
i
]
⊗
n
(
s
1
,
…
,
s
^
i
,
…
,
s
n
+
1
)
×
I
[
0
,
t
j
]
⊗
n
(
t
1
,
…
,
t
^
j
,
…
,
t
n
+
1
)
·
∏
q
=
1
n
+
1
‍
∂
2
∂
s
q
∂
t
q
R
H
,
K
×
(
s
q
,
t
q
)
d
s
1
⋯
d
s
n
+
1
d
t
1
⋯
d
t
n
+
1
d
x
=
∫
ℝ
‍
(
1
+
x
2
)
-
r
∑
n
≥
0
‍
x
2
n
(
n
+
1
)
!
×
∑
i
=
1
n
+
1
‍
∬
[
0
,
T
]
n
+
1
exp
⁡
{
-
x
2
2
s
i
2
H
K
}
·
exp
⁡
{
-
x
2
2
t
i
2
H
K
}
I
[
0
,
s
i
]
⊗
n
(
s
1
,
…
,
s
^
i
,
…
,
s
n
+
1
)
×
I
[
0
,
t
i
]
⊗
n
(
t
1
,
…
,
t
^
i
,
…
,
t
n
+
1
)
·
∏
q
=
1
n
+
1
‍
∂
2
∂
s
q
∂
t
q
R
H
,
K
(
s
q
,
t
q
)
d
s
1
⋯
d
s
n
+
1
d
t
1
⋯
d
t
n
+
1
d
x
+
∫
ℝ
‍
(
1
+
x
2
)
-
r
∑
n
≥
0
‍
x
2
n
(
n
+
1
)
!
×
∑
i
,
j
=
1
,
i
≠
j
n
+
1
‍
∬
[
0
,
T
]
n
+
1
exp
⁡
{
-
x
2
2
s
i
2
H
K
}
·
exp
⁡
⁡
{
-
x
2
2
t
j
2
H
K
}
×
I
[
0
,
s
i
]
⊗
n
(
s
1
,
…
,
s
^
i
,
…
,
s
n
+
1
)
×
I
[
0
,
t
j
]
⊗
n
(
t
1
,
…
,
t
^
j
,
…
,
t
n
+
1
)
·
∏
q
=
1
n
+
1
‍
∂
2
∂
s
q
∂
t
q
×
R
H
,
K
(
s
q
,
t
q
)
d
s
1
⋯
d
s
n
+
1
d
t
1
⋯
d
t
n
+
1
d
x
≡
Δ
1
+
Δ
2
.
Firstly, we turn to estimate Δ
1. Using the similar methods in [7, 11] and the following fact:
(26)
∬
0
T
I
[
0
,
t
]
(
u
)
I
[
0
,
s
]
(
v
)
∂
2
∂
u
∂
v
R
H
,
K
(
u
,
v
)
d
u
d
v
=
R
H
,
K
(
t
,
s
)
,
we find
(27)
Δ
1
=
∫
ℝ
‍
(
1
+
x
2
)
-
r
∑
n
≥
0
‍
x
2
n
(
n
+
1
)
!
×
∑
i
=
1
n
+
1
∬
[
0
,
T
]
n
+
1
exp
⁡
{
-
x
2
2
s
1
2
H
K
}
‍
·
exp
⁡
{
-
x
2
2
t
1
2
H
K
}
I
[
0
,
s
1
]
⊗
n
(
s
2
,
…
,
s
^
i
,
…
,
s
n
+
1
)
×
I
[
0
,
t
1
]
⊗
n
(
t
2
,
…
,
t
^
j
,
…
,
t
n
+
1
)
·
∏
q
=
2
n
+
1
‍
∂
2
∂
s
q
∂
t
q
R
H
,
K
(
s
q
,
t
q
)
∂
2
∂
s
1
∂
t
1
×
R
H
,
K
(
s
1
,
t
1
)
d
s
1
·
·
·
d
s
n
+
1
d
t
1
·
·
·
d
t
n
+
1
d
x
=
∫
ℝ
‍
(
1
+
x
2
)
-
r
×
∑
n
≥
0
‍
x
2
n
n
!
∬
[
0
,
T
]
exp
⁡
⁡
{
-
x
2
2
s
1
2
H
K
}
·
exp
⁡
{
-
x
2
2
t
1
2
H
K
}
R
H
,
K
n
(
s
1
,
t
1
)
∂
2
∂
s
1
∂
t
1
×
R
H
,
K
(
s
1
,
t
1
)
d
s
1
d
t
1
d
x
=
∫
ℝ
‍
(
1
+
x
2
)
-
r
∬
[
0
,
T
]
exp
⁡
{
-
x
2
2
s
1
2
H
K
}
exp
⁡
{
-
x
2
2
t
1
2
H
K
}
·
∑
n
≥
0
‍
x
2
n
n
!
R
H
,
K
n
(
s
1
,
t
1
)
∂
2
∂
s
1
∂
t
1
×
R
H
,
K
(
s
1
,
t
1
)
d
s
1
d
t
1
d
x
=
∫
ℝ
‍
(
1
+
x
2
)
-
r
∬
[
0
,
T
]
exp
⁡
⁡
{
-
x
2
2
s
1
2
H
K
}
exp
⁡
⁡
{
-
x
2
2
t
1
2
H
K
}
·
exp
⁡
{
x
2
R
H
,
K
(
s
1
,
t
1
)
}
×
∂
2
∂
s
1
∂
t
1
R
H
,
K
(
s
1
,
t
1
)
d
s
1
d
t
1
d
x
,
where the last equality is established due to Taylor expansion formula of exponential function.
Since(28)
2
R
H
,
K
(
s
1
,
t
1
)
-
s
1
2
H
K
-
t
1
2
H
K
=
-
E
[
(
B
s
1
H
,
K
-
B
t
1
H
,
K
)
2
]
,
we get
(29)
exp
⁡
{
x
2
2
(
2
R
H
,
K
(
s
1
,
t
1
)
-
s
1
2
H
K
-
t
1
2
H
K
)
}
=
exp
⁡
{
-
x
2
2
E
[
(
B
s
1
H
,
K
-
B
t
1
H
,
K
)
2
]
}
.
By [5], for each s
1
,
t
1
∈
[
0
,
T
], we have
(30)
2
-
K
|
s
1
-
t
1
|
2
H
K
≤
E
[
(
B
s
1
H
,
K
-
B
t
1
H
,
K
)
2
]
≤
2
1
-
K
|
s
1
-
t
1
|
2
H
K
,
which implies that
(31)
E
[
(
B
s
1
H
,
K
-
B
t
1
H
,
K
)
2
]
≥
0
.
Use the change of variablesy
=
x
{
E
[
(
B
s
1
H
,
K
-
B
t
1
H
,
K
)
2
]
}
1
/
2. Furthermore, we have
(32)
d
x
=
(
s
1
2
H
K
+
t
1
2
H
K
-
1
2
K
-
1
×
[
(
s
1
2
H
+
t
1
2
H
)
K
|
s
1
-
t
1
|
2
H
K
]
s
1
2
H
K
+
t
1
2
H
K
-
1
2
K
-
1
)
-
1
/
2
d
y
,
(
1
+
x
2
)
-
r
=
{
×
[
(
s
1
2
H
+
t
1
2
H
)
K

-
|
s
1
-
t
1
|
2
H
K
]
1
2
K
-
1
)
-
1

y
2

1
+
(
s
1
2
H
K
+
t
1
2
H
K
-
1
2
K
-
1
×
[
(
s
1
2
H
+
t
1
2
H
)
K

-
|
s
1
-
t
1
|
2
H
K
]
1
2
K
-
1
)
-
1

y
2

}
-
r

.
On the other hand, by [11], there exists a constant C
3,1
(
H
,
K
) depending on H and K such that
(33)
|
∂
2
∂
s
1
∂
t
1
R
H
,
K
(
s
1
,
t
1
)
|
≤
C
3,1
(
H
,
K
)
(
s
1
t
1
)
H
K
-
1
.
Putting (32)-(33) into (27), calculate
(34)
Δ
1
≤
C
3,1
(
H
,
K
)
∬
0
T
(
s
1
t
1
)
H
K
-
1
×
(
s
1
2
H
K
+
t
1
2
H
K
-
1
2
K
-
1
×
[
(
s
1
2
H
+
t
1
2
H
)
K
-
|
s
1
-
t
1
|
2
H
K
]
1
2
K
-
1
)
r
-
1
/
2
×
∫
ℝ
‍
exp
⁡
{
-
y
2
2
}
·
(
s
1
2
H
K
+
t
1
2
H
K
-
1
2
K
-
1
·
[
(
s
1
2
H
+
t
1
2
H
)
K
-
|
s
1
-
t
1
|
2
H
K
]
+
y
2
s
1
2
H
K
+
t
1
2
H
K
-
1
2
K
-
1
)
-
r
d
y
d
s
1
d
t
1
≤
C
3,1
(
H
,
K
)
×
∬
0
T
(
s
1
t
1
)
H
K
-
1
(
s
1
2
H
K
+
t
1
2
H
K
-
1
2
K
-
1
×
[
(
s
1
2
H
+
t
1
2
H
)
K
-
|
s
1
-
t
1
|
2
H
K
]
1
2
K
-
1
)
r
-
1
/
2
×
∫
ℝ
‍
exp
⁡
{
-
y
2
2
}
y
-
2
r
d
y
d
s
1
d
t
1
.
When H
K
-
1
+
2
H
K
(
r
-
1
/
2
)
>
-
1, that is, r
>
0, (34) is finite.
Secondly, using the similar estimation method in the first part, consider the estimation ofΔ
2 as follows:
(35)
Δ
2
=
∫
ℝ
‍
(
1
+
x
2
)
-
r
∑
n
≥
1
‍
x
2
n
(
n
+
1
)
!
n
(
n
+
1
)
×
∬
[
0
,
T
]
n
+
1
exp
⁡
{
-
x
2
2
s
1
2
H
K
}
·
exp
⁡
{
-
x
2
2
t
2
2
H
K
}
I
[
0
,
s
1
]
⊗
n
(
s
2
,
…
,
s
n
+
1
)
×
I
[
0
,
t
2
]
⊗
n
(
t
1
,
t
3
,
…
,
t
n
+
1
)
·
∏
q
=
1
n
+
1
‍
∂
2
∂
s
q
∂
t
q
×
R
H
,
K
(
s
q
,
t
q
)
d
s
1
·
·
·
d
s
n
+
1
d
t
1
·
·
·
d
t
n
+
1
d
x
=
∫
ℝ
‍
(
1
+
x
2
)
-
r
∑
n
≥
1
‍
x
2
n
(
n
-
1
)
!
×
∭
∫
0
T
exp
⁡
⁡
{
-
x
2
2
s
1
2
H
K
}
exp
⁡
⁡
{
-
x
2
2
t
2
2
H
K
}
·
I
[
0
,
s
1
]
(
s
2
)
I
[
0
,
t
2
]
(
t
1
)
∂
2
∂
s
1
∂
t
1
×
R
H
,
K
(
s
1
,
t
1
)
∂
2
∂
s
2
∂
t
2
R
H
,
K
(
s
2
,
t
2
)
·
R
H
,
K
n
-
1
(
s
1
,
t
2
)
d
s
1
d
s
2
d
t
1
d
t
2
d
x
.
By Taylor expansion formula, the following equality is obvious:
(36)
Δ
2
=
∫
ℝ
‍
(
1
+
x
2
)
-
r
x
2
×
∭
∫
0
T
exp
⁡
{
-
x
2
2
(
s
1
2
H
K
+
t
2
2
H
K
-
R
H
,
K
(
s
1
,
t
2
)
)
}
·
I
[
0
,
s
1
]
(
s
2
)
I
[
0
,
t
2
]
(
t
1
)
∂
2
∂
s
1
∂
t
1
R
H
,
K
(
s
1
,
t
1
)
×
∂
2
∂
s
2
∂
t
2
R
H
,
K
(
s
2
,
t
2
)
·
R
H
,
K
n
-
1
(
s
1
,
t
2
)
d
s
1
d
s
2
d
t
1
d
t
2
d
x
.
Recalling some results in [11], there exist parameters H
,
K and a constant C
3,2
(
H
,
K
) depending on H and K such that
(37)
|
∂
2
∂
s
2
∂
t
2
R
H
,
K
(
s
2
,
t
2
)
|
≤
C
3,2
(
H
,
K
)
(
s
2
t
2
)
H
K
-
1
.
Thus, by (33) and (37), we obtain
(38)
Δ
2
≤
C
3,1
(
H
,
K
)
C
3,2
(
H
,
K
)
×
∫
ℝ
‍
(
1
+
x
2
)
-
r
×
∭
∫
0
T
x
2
·
exp
⁡
{
-
x
2
2
[
s
1
2
H
K
+
t
2
2
H
K
-
1
2
K
-
1
×
(
(
s
1
2
H
+
t
2
2
H
)
K
-
|
s
1
-
t
2
|
2
H
K
)
s
1
2
H
K
+
t
2
2
H
K
-
1
2
K
-
1
]
-
x
2
2
[
s
1
2
H
K
+
t
2
2
H
K
-
1
2
K
-
1
}
·
(
s
1
t
1
)
H
K
-
1
(
s
2
t
2
)
H
K
-
1
I
[
0
,
s
1
]
(
s
2
)
×
I
[
0
,
t
2
]
(
t
1
)
d
s
1
d
s
2
d
t
1
d
t
2
d
x
.
Calculate
(39)
∭
∫
0
T
exp
⁡
{
-
x
2
2
[
s
1
2
H
K
+
t
2
2
H
K
-
1
2
K
-
1
×
(
(
s
1
2
H
+
t
2
2
H
)
K
-
|
s
1
-
t
2
|
2
H
K
)
s
1
2
H
K
+
t
2
2
H
K
-
1
2
K
-
1
]
-
x
2
2
}
·
(
s
1
t
1
)
H
K
-
1
(
s
2
t
2
)
H
K
-
1
I
[
0
,
s
1
]
(
s
2
)
×
I
[
0
,
t
2
]
(
t
1
)
d
s
1
d
s
2
d
t
1
d
t
2
=
∬
0
T
exp
⁡
{
-
x
2
2
[
s
1
2
H
K
+
t
2
2
H
K
-
1
2
K
-
1
×
(
(
s
1
2
H
+
t
2
2
H
)
K
-
|
s
1
-
t
2
|
2
H
K
)
s
1
2
H
K
+
t
2
2
H
K
-
1
2
K
-
1
]
-
x
2
2
}
·
(
s
1
t
2
)
H
K
-
1
(
∫
0
T
‍
s
2
H
K
-
1
I
[
0
,
s
1
]
(
s
2
)
d
s
2
)
×
(
∫
0
T
‍
t
1
H
K
-
1
I
[
0
,
t
2
]
(
t
1
)
d
t
1
)
d
s
1
d
t
2
=
1
(
H
K
)
2
×
∬
0
T
exp
⁡
{
-
x
2
2
[
s
1
2
H
K
+
t
2
2
H
K
-
1
2
K
-
1
×
(
(
s
1
2
H
+
t
2
2
H
)
K
-
|
s
1
-
t
2
|
2
H
K
)
s
1
2
H
K
+
t
2
2
H
K
-
1
2
K
-
1
]
-
x
2
2
}
×
(
s
1
t
2
)
2
H
K
-
1
d
s
1
d
t
2
.
Use the change of variables y
=
[
s
1
2
H
K
+
t
2
2
H
K
-
(
1
/
2
K
-
1
)
(
(
s
1
2
H
+
t
2
2
H
)
K
-
|
s
1
-
t
2
|
2
H
K
)
]
1
/
2
x. Hence
(40)
Δ
2
≤
C
3,3
(
H
,
K
)
×
∬
0
T
(
s
1
t
2
)
2
H
K
-
1
(
s
1
2
H
K
+
t
2
2
H
K
-
1
2
K
-
1
×
[
(
s
1
2
H
+
t
2
2
H
)
K
-
|
s
1
-
t
2
|
2
H
K
]
s
1
2
H
K
+
t
2
2
H
K
-
1
2
K
-
1
)
r
-
3
/
2
×
∫
ℝ
‍
exp
⁡
{
-
y
2
2
}
·
(
s
1
2
H
K
+
t
2
2
H
K
-
1
2
K
-
1
×
[
(
s
1
2
H
+
t
2
2
H
)
K
-
|
s
1
-
t
2
|
2
H
K
]
+
y
2
)
-
r
d
y
d
s
1
d
t
2
≤
C
3,3
(
H
,
K
)
∬
0
T
(
s
1
t
2
)
2
H
K
-
1
×
(
s
1
2
H
K
+
t
2
2
H
K
-
1
2
K
-
1
×
[
(
s
1
2
H
+
t
2
2
H
)
K
-
|
s
1
-
t
2
|
2
H
K
]
s
1
2
H
K
+
t
2
2
H
K
-
1
2
K
-
1
)
r
-
3
/
2
×
∫
ℝ
‍
exp
⁡
{
-
y
2
2
}
y
-
2
r
d
y
d
s
1
d
t
2
,
where C
3,3
(
H
,
K
)
=
C
3,1
(
H
,
K
)
C
3,2
(
H
,
K
).
When2
H
K
-
1
+
2
H
K
(
r
-
3
/
2
)
>
-
1, (40) is finite. In other words, when r
>
1
/
2, (40) is finite.
From what we have said above we can draw a conclusion that whenr
>
1
/
2,
(41)
E
∥
ξ
∥
H
r
(
ℝ
;
ℝ
)
2
<
∞
.From Theorem2 we see that when r
>
1
/
2, the mapping ξ
(
x
)
=
∫
[
0
,
T
]
N
‍
δ
(
x
-
B
s
H
,
K
)
d
B
s
H
,
K belongs to negative Sobolev space H
-
r
(
ℝ
;
ℝ
). Note that the regularity condition in Theorem 2 does not depend on H and K. The condition is interesting, because the condition of ξ
(
x
)
=
∫
[
0
,
T
]
N
‍
δ
(
x
-
B
s
)
d
B
s which belongs to negative Sobolev space is also r
>
1
/
2, where B
s is the Brownian motion (see [7]). In other words, they have the same regularity condition. However, the situation is different in the case of fractional Brownian motion, because the regularity condition of fractional Brownian stochastic current is r
>
1
/
2
H
-
1
/
2, which is dependent on Hurst parameter H.


## 3.2. Stochastic Current of d-Dimensional Case with respect tox

As in [7], we can extend stochastic current of one-dimensional bifractional Brownian motion to the case of d-dimensional bifractional Brownian motion.LetB
H
,
K be the vector valued bifractional Brownian motion; that is, B
H
,
K
=
(
B
H
1
,
K
1
,
…
,
B
H
d
,
K
d
), where B
H
i
,
K
i are independent one-dimensional bifractional Brownian sheet. In this part, we consider ξ
(
x
) as follows:
(42)
ξ
(
x
)
=
(
∫
[
0
,
T
]
N
‍
δ
(
x
-
B
s
H
,
K
)
d
B
s
H
1
,
K
1
,
…
,
∫
[
0
,
T
]
N
‍
δ
(
x
-
B
s
H
,
K
)
d
B
s
H
d
,
K
d
)
,
where the integrals are Skorohod integrals with respect to bifractional Brownian motion.Denote(43)
R
H
,
K
l
(
t
,
s
)
=
∏
i
=
1
N
‍
R
H
,
K
l
(
t
i
,
s
i
)
=
〈
I
[
0
,
t
]
,
I
[
0
,
s
]
〉
ℋ
H
l
.Theorem 3.
LetB
H
,
K be d-dimensional bifractional Brownian motion with parameters H
i and K
i satisfying 2
H
i
K
i
>
1
(
i
=
1
,
…
,
d
) and let ξ
(
x
) be given by (42); then, for each ω
∈
Ω and when r
>
d
/
2
-
1, ξ
(
x
) belongs to Sobolev space H
-
r
(
ℝ
;
ℝ
).Proof.
Denoteξ
l
(
x
) by
(44)
ξ
l
(
x
)
=
∫
[
0
,
T
]
N
‍
δ
(
x
-
B
s
H
,
K
)
d
B
s
H
l
,
K
l
=
∫
[
0
,
T
]
N
‍
δ
l
(
x
-
B
s
H
,
K
)
δ
(
x
l
-
B
s
H
l
,
K
l
)
d
B
s
H
l
,
K
l
=
∫
[
0
,
T
]
N
‍
δ
l
(
x
-
B
s
H
,
K
)
×
∑
n
l
≥
0
‍
β
n
l
x
l
(
s
)
I
n
l
B
H
l
,
K
l
(
I
[
0
,
s
]
⊗
n
l
(
·
)
)
d
B
s
H
l
,
K
l
,
where δ
l
(
x
-
B
s
H
,
K
)
=
∏
q
=
1
,
q
≠
l
d
‍
δ
(
x
q
-
B
s
H
q
,
K
q
).
Calculate the Fourier transform of (44) as follows:
(45)
ξ
^
l
(
x
)
=
∑
n
l
≥
0
‍
I
n
l
+
1
B
H
l
,
K
l
(
(
exp
⁡
{
-
i
∑
r
=
1
,
r
≠
l
d
‍
x
r
B
s
H
r
,
K
r
}
×
(
-
i
)
n
l
x
l
n
l
n
l
!
×
exp
⁡
⁡
{
-
x
2
2
|
s
|
2
H
l
K
l
}
I
[
0
,
s
]
⊗
n
l
(
·
)
)
(
s
)
)
.
According to the definition of the normal ∥
ξ
(
x
)
∥
H
-
r
(
ℝ
;
ℝ
)
2 and using Euler formula, we can prove that
(46)
E
|
ξ
^
l
(
x
)
|
2
=
∑
n
l
≥
0
‍
x
l
2
n
l
(
n
l
+
1
)
!
×
∑
i
,
j
=
1
n
l
+
1
‍
∬
[
0
,
T
]
N
(
n
l
+
1
)
∏
q
=
1
n
l
+
1
‍
∂
2
∂
s
∂
t
R
H
,
K
q
(
s
,
t
)
·
cos
⁡
(
∑
r
=
1
,
r
≠
l
d
‍
x
r
B
s
i
H
r
,
K
r
)
cos
⁡
(
∑
r
=
1
,
r
≠
l
d
‍
x
r
B
t
j
H
r
,
K
r
)
×
exp
⁡
{
-
x
l
2
2
|
s
i
|
2
H
l
K
l
}
·
exp
⁡
⁡
{
-
x
l
2
2
|
t
j
|
2
H
l
K
l
}
×
I
[
0
,
s
i
]
⊗
n
l
(
s
1
,
…
,
s
^
i
,
…
,
s
n
l
+
1
)
×
I
[
0
,
t
j
]
⊗
n
l
(
t
1
,
…
,
t
^
j
,
…
,
t
n
l
+
1
)
d
s
d
t
+
∑
n
l
≥
0
‍
x
l
2
n
l
(
n
l
+
1
)
!
×
∑
i
,
j
=
1
n
l
+
1
∬
[
0
,
T
]
N
(
n
l
+
1
)
∏
q
=
1
n
l
+
1
‍
∂
2
∂
s
∂
t
R
H
,
K
q
(
s
,
t
)
‍
·
sin
(
∑
r
=
1
,
r
≠
l
d
‍
x
r
B
s
i
H
r
,
K
r
)
sin
(
∑
r
=
1
,
r
≠
l
d
‍
x
r
B
t
j
H
r
,
K
r
)
×
exp
⁡
{
-
x
l
2
2
|
s
i
|
2
H
l
K
l
}
·
exp
⁡
{
-
x
l
2
2
|
t
j
|
2
H
l
K
l
}
I
[
0
,
s
i
]
⊗
n
l
(
s
1
,
…
,
s
^
i
,
…
,
s
n
l
+
1
)
×
I
[
0
,
t
j
]
⊗
n
l
(
t
1
,
…
,
t
^
j
,
…
,
t
n
l
+
1
)
d
s
d
t
,
where |
s
i
|
2
H
l
K
l
=
∏
j
=
1
N
‍
(
s
j
i
)
2
H
l
K
l
,
I
[
0
,
s
i
]
=
∏
j
=
1
N
‍
I
[
0
,
s
j
i
] and R
H
,
K
q
(
s
,
t
)
=
R
H
q
,
K
q
(
s
q
,
t
q
).
By the bound ofsin
(
x
) and cos
⁡
(
x
), we can get the following inequality:
(47)
E
|
ξ
^
l
(
x
)
|
2
≤
C
3,4
∑
n
l
≥
0
‍
x
l
2
n
l
(
n
l
+
1
)
!
×
∑
i
,
j
=
1
n
l
+
1
‍
∬
[
0
,
T
]
N
(
n
l
+
1
)
∏
q
=
1
n
l
+
1
‍
∂
2
∂
s
∂
t
R
H
,
K
q
(
s
,
t
)
·
exp
⁡
{
-
x
l
2
2
|
s
i
|
2
H
l
K
l
}
×
exp
⁡
{
-
x
l
2
2
|
t
j
|
2
H
l
K
l
}
·
I
[
0
,
s
i
]
⊗
n
l
(
s
1
,
…
,
s
^
i
,
…
,
s
n
l
+
1
)
×
I
[
0
,
t
j
]
⊗
n
l
(
t
1
,
…
,
t
^
j
,
…
,
t
n
l
+
1
)
d
s
d
t
+
C
3,5
∑
n
l
≥
0
‍
x
l
2
n
l
(
n
l
+
1
)
!
×
∑
i
,
j
=
1
n
l
+
1
‍
∬
[
0
,
T
]
N
(
n
l
+
1
)
∏
q
=
1
n
l
+
1
‍
∂
2
∂
s
∂
t
R
H
,
K
q
(
s
,
t
)
·
exp
⁡
{
-
x
l
2
2
|
s
i
|
2
H
l
K
l
}
exp
⁡
{
-
x
l
2
2
|
t
j
|
2
H
l
K
l
}
·
I
[
0
,
s
i
]
⊗
n
l
(
s
1
,
…
,
s
^
i
,
…
,
s
n
l
+
1
)
×
I
[
0
,
t
j
]
⊗
n
l
(
t
1
,
…
,
t
^
j
,
…
,
t
n
l
+
1
)
d
s
d
t
≡
Δ
3
,
l
+
Δ
4
,
l
,
where C
3,4 and C
3,5 are both constants.
By the same estimation techniques ofΔ
2 in Theorem 2, we can obtain the estimation of Δ
4
,
l. Here we need to discuss the estimation of Δ
3
,
l.
Applying∬
0
T
I
[
0
,
t
]
(
u
)
I
[
0
,
s
]
(
v
)
(
∂
2
/
∂
u
∂
v
)
R
H
,
K
(
u
,
v
)
d
u
d
v
=
R
H
,
K
(
t
,
s
) again, we can write
(48)
Δ
3
,
l
=
C
3,4
∑
n
l
≥
0
‍
x
l
2
n
l
n
l
!
×
∬
[
0
,
T
]
N
(
n
l
+
1
)
∏
q
=
1
n
l
+
1
‍
∂
2
∂
s
1
∂
t
1
R
H
,
K
q
(
s
1
,
t
1
)
·
exp
⁡
{
-
x
l
2
2
|
s
1
|
2
H
l
K
l
}
exp
⁡
{
-
x
l
2
2
|
t
1
|
2
H
l
K
l
}
·
I
[
0
,
s
1
]
⊗
n
l
(
s
2
,
…
,
s
n
l
+
1
)
×
I
[
0
,
t
1
]
⊗
n
l
(
t
2
,
…
,
t
n
l
+
1
)
d
s
d
t
=
C
3,4
∑
n
l
≥
0
‍
x
l
2
n
l
n
l
!
×
∬
[
0
,
T
]
N
(
R
H
,
K
l
(
s
1
,
t
1
)
)
n
l
·
∂
2
∂
s
∂
t
R
H
,
K
l
(
s
1
,
t
1
)
×
exp
⁡
{
-
x
l
2
2
|
s
1
|
2
H
l
K
l
}
exp
⁡
{
-
x
l
2
2
|
t
1
|
2
H
l
K
l
}
d
s
1
d
t
1
=
C
3,4
∬
[
0
,
T
]
N
∂
2
∂
s
1
∂
t
1
R
H
,
K
l
(
s
1
,
t
1
)
·
exp
⁡
{
x
l
2
2
(
|
s
1
|
2
H
l
K
l
-
|
t
1
|
2
H
l
K
l
2
R
H
,
K
l
(
s
1
,
t
1
)
-
|
s
1
|
2
H
l
K
l
-
|
t
1
|
2
H
l
K
l
)
x
l
2
2
}
d
s
1
d
t
1
.
According to [11], there exists a constant C
3,6
(
H
,
K
,
l
) depending upon H
,
K, and l such that
(49)
|
∂
2
∂
s
1
∂
t
1
R
H
,
K
l
(
s
1
,
t
1
)
|
≤
C
3,6
(
H
,
K
,
l
)
(
s
1
t
1
)
H
l
K
l
-
1
.
Therefore
(50)
∫
ℝ
d
‍
(
1
+
|
x
|
2
)
-
r
Δ
3
,
l
d
x
≤
C
3,6
(
H
,
K
,
l
)
×
∫
ℝ
d
‍
∬
[
0
,
T
]
N
(
s
1
t
1
)
H
l
K
l
-
1
·
exp
⁡
{
-
x
2
2
E
[
(
B
s
1
H
l
,
K
l
-
B
t
1
H
l
,
K
l
)
2
]
}
d
s
1
d
t
1
d
x
.
Use the change of variables y
l
=
x
l
{
E
[
(
B
s
1
H
l
,
K
l
-
B
t
1
H
l
,
K
l
)
2
]
}
1
/
2. Thus
(51)
∫
ℝ
d
‍
(
1
+
|
x
|
2
)
-
r
Δ
3
,
l
d
x
≤
C
3,7
(
H
,
K
,
l
)
×
∬
[
0
,
T
]
N
(
s
1
t
1
)
H
l
K
l
-
1
{
E
[
(
B
s
1
H
l
,
K
l
-
B
t
1
H
l
,
K
l
)
2
]
}
r
-
d
/
2
·
∫
ℝ
d
‍
exp
⁡
{
-
y
l
2
2
}
×
(
E
[
(
B
s
1
H
l
,
K
l
-
B
t
1
H
l
,
K
l
)
2
]
+
y
l
2
)
-
r
d
s
1
d
t
1
d
y
,
where y
=
(
y
1
,
…
,
y
d
).
On the other hand, by [6] for arbitrary ɛ
≥
0
,
s
1
,
t
1
∈
[
ɛ
,
1
]
N, there exists a constant C
3,8 such that
(52)
E
[
(
B
s
1
H
l
,
K
l
-
B
t
1
H
l
,
K
l
)
2
]
≤
C
3,8
∑
j
=
1
N
‍
|
s
j
1
-
t
j
1
|
2
H
l
K
l
.
In order to be simple, here we only consider the case of T
=
1.
Comparing (51) with (52), we find that
(53)
∫
ℝ
d
‍
(
1
+
|
x
|
2
)
-
r
Δ
3
,
l
d
x
≤
C
3,9
∬
[
0
,
T
]
N
(
s
1
t
1
)
H
l
K
l
-
1
·
(
∑
j
=
1
N
‍
|
s
j
1
-
t
j
1
|
2
H
j
K
j
)
r
-
d
/
2
×
∫
ℝ
d
‍
exp
⁡
{
-
y
l
2
2
}
y
-
2
r
d
y
d
s
d
t
.
When H
l
K
l
-
1
+
2
H
l
K
l
(
r
-
d
/
2
)
>
-
1, that is, r
>
d
/
2
-
1, (53) is finite. That is to say, when r
>
d
/
2
-
1, there is
(54)
∫
ℝ
d
‍
(
1
+
|
x
|
2
)
-
r
E
[
|
ξ
^
(
x
)
|
2
]
d
x
=
∫
ℝ
d
‍
(
1
+
|
x
|
2
)
-
r
E
[
|
ξ
^
1
(
x
)
|
2
+
·
·
·
+
|
ξ
^
d
(
x
)
|
2
]
d
x
<
∞
.It is interesting to contrast Theorems2 and 3 with Propositions 3 and 4 in [7]. Parameters H and K of bifractional Brownian motion do not affect the regularity condition of bifractional Brownian current. However, the regularity condition of stochastic current of one-dimensional fractional Brownian motion is different from the case of d-dimensional setting (see [7]). In other words, Hurst parameters H of fractional Brownian motion have influence on fractional Brownian currents in the case of different dimension.


## 3.3. Stochastic Current of d-Dimensional Bifractional Brownian Motion with respect toω

LetB
H
,
K be vector valued bifractional Brownian motion where vectors H
=
(
H
1
,
…
,
H
d
)
∈
(
0,1
)
d and K
=
(
K
1
,
…
,
K
d
)
∈
(
0,1
]
d; that is, B
H
,
K
=
(
B
H
1
,
K
1
,
…
,
B
H
d
,
K
d
). In this part, ξ
(
x
) is given by
(55)
ξ
(
x
)
=
∫
0
T
‍
δ
(
x
-
B
s
H
,
K
)
d
B
s
H
,
K
,
where
(56)
δ
(
x
-
B
s
H
,
K
)
=
∏
i
=
1
d
‍
δ
(
x
i
-
B
s
H
i
,
K
i
)
=
∑
n
=
(
n
1
,
…
,
n
d
)
‍
β
n
(
s
,
x
)
I
n
B
H
i
,
K
i
(
I
[
0
,
s
]
⊗
|
n
|
(
·
)
)
,
β
n
(
s
,
x
)
=
∏
i
=
1
d
‍
1
R
H
,
K
i
(
s
)
n
i
/
2
p
R
H
,
K
i
(
x
i
)
H
n
i
(
x
i
R
H
,
K
i
(
s
)
1
/
2
)
,
R
H
,
K
i
(
s
)
=
R
H
,
K
i
(
s
,
s
)
=
s
2
H
i
K
i
,
I
n
B
H
,
K
(
I
[
0
,
s
]
⊗
|
n
|
(
·
)
)
=
∏
i
=
1
d
‍
I
n
i
B
H
i
,
K
i
(
I
[
0
,
s
]
⊗
n
i
(
·
)
)
.Using the chaos expansion of divergent integral with respect to bifractional Brownian motion, we can obtain the following expression:(57)
ξ
i
(
x
)
=
∫
0
T
‍
δ
(
x
-
B
s
H
,
K
)
d
B
s
H
i
,
K
i
=
∑
n
=
(
n
1
,
…
,
n
d
)
‍
I
n
i
+
1
B
H
i
,
K
i
[
∏
j
=
1
,
j
≠
i
d
‍
I
n
j
B
H
j
,
K
j
(
I
[
0
,
s
]
⊗
n
j
(
·
)
)
)
(
s
)
(
β
n
(
s
,
x
)
I
[
0
,
s
]
⊗
n
i
(
s
1
,
…
,
s
n
i
)
I
[
0
,
T
]
(
s
)
·
∏
j
=
1
,
j
≠
i
d
‍
I
n
j
B
H
j
,
K
j
(
I
[
0
,
s
]
⊗
n
j
(
·
)
)
)
(
s
)
]
.Theorem 4.
LetB
H
,
K be a bifractional Brownian motion with parameters H
=
(
H
1
,
…
,
H
d
)
∈
(
0,1
)
d and K
=
(
K
1
,
…
,
K
d
)
∈
(
0,1
]
d satisfying 2
H
i
K
i
>
1
(
i
=
1
,
…
,
d
). If ξ
(
x
) is given by (55), then, for every x
∈
ℝ
d and α
<
1
/
2
(
H
K
)
*
-
d
/
2, ξ
(
x
) is one member of Sobolev-Watanabe spaces D
α
-
1,2, where (
H
K
)
*
=
max
⁡
{
H
1
K
1
,
…
,
H
d
K
d
}.Proof.
In order to be convenient, we always replace the normal∥
·
∥
ℋ
⊗
(
n
+
1
) with the normal ∥
·
∥
ℋ
1
⊗
(
n
+
1
). Using the chaos expansion of δ
(
x
-
B
H
,
K
) and by the definition of the normal of ∥
ξ
(
x
)
∥
2
,
α
-
1
2, we verify that
(58)
∥
ξ
i
(
x
)
∥
2
,
α
-
1
2
≤
∑
m
≥
1
‍
(
m
+
1
)
α
-
1
∑
|
n
|
=
n
1
+
·
·
·
+
n
d
=
m
-
1
‍
(
n
i
+
1
)
!
·
∬
[
0
,
T
]
n
i
+
1
∑
l
,
k
=
1
n
i
+
1
‍
β
n
(
s
l
)
β
n
(
t
k
)
I
[
0
,
T
]
(
s
l
)
I
[
0
,
T
]
(
t
k
)
·
I
[
0
,
s
l
]
⊗
n
i
(
s
1
,
…
,
s
^
l
,
…
,
s
n
i
+
1
)
×
I
[
0
,
t
k
]
⊗
n
i
(
t
1
,
…
,
t
^
k
,
…
,
t
n
i
+
1
)
·
∏
j
=
1
,
j
≠
i
d
‍
n
j
!
R
H
,
K
j
(
s
l
,
t
k
)
n
j
×

∏
q
=
1
n
i
+
1
‍

∂
2
∂
s
q
∂
t
q
R
H
,
K
i

(
s
q
,
t
q
)

d
t
1
·
·
·
d
t
n
i
+
1
d
s
1
·
·
·
d
s
n
i
+
1
.
Because α
<
0, the following inequalities: (
m
+
2
)
α
-
1
≤
(
m
+
1
)
α
-
1 and n
i
+
1
≤
m
+
1, are obvious. Similarly to what is performed in [11] (or see [7]), we have
(59)
∥
ξ
i
(
x
)
∥
2
,
α
-
1
2
≤
∑
m
≥
0
‍
(
m
+
1
)
α
∑
|
n
|
=
n
1
+
·
·
·
+
n
d
=
m
‍
n
i
!
×
[
(
1
-
1
n
i
+
1
)
×
∬
[
0
,
T
]
2
β
n
(
s
1
)
β
n
(
t
2
)
·
I
[
0
,
T
]
(
s
1
)
I
[
0
,
T
]
(
t
2
)
×
R
H
,
K
i
(
s
1
,
t
2
)
n
i
-
1
(
1
-
1
n
i
+
1
)
]
×
∏
j
=
1
,
j
≠
i
d
‍
n
j
!
R
H
,
K
j
(
s
1
,
t
2
)
n
j
·
I
[
0
,
s
1
]
(
s
2
)
×
I
[
0
,
t
2
]
(
t
1
)
∂
2
∂
s
1
∂
t
1
R
H
,
K
i
(
s
1
,
t
1
)
×
∂
2
∂
s
2
∂
t
2
R
H
,
K
i
(
s
2
,
t
2
)
d
s
1
d
s
2
d
t
1
d
t
2
+
1
n
i
+
1
×
∬
[
0
,
T
]
2
β
n
(
s
1
)
β
n
(
t
1
)
I
[
0
,
T
]
(
s
1
)
×
I
[
0
,
T
]
(
t
1
)
R
H
,
K
i
(
s
1
,
t
1
)
n
i
-
1
·
∏
j
=
1
,
j
≠
i
d
‍
n
j
!
R
H
,
K
j
(
s
1
,
t
1
)
n
j
I
[
0
,
s
1
]
(
s
2
)
×
I
[
0
,
t
1
]
(
t
2
)
·
∂
2
∂
s
1
∂
t
1
R
H
,
K
i
(
s
1
,
t
1
)
∂
2
∂
s
2
∂
t
2
×
R
H
,
K
i
(
s
2
,
t
2
)
d
s
1
d
s
2
d
t
1
d
t
2
≤
∑
m
≥
0
‍
(
m
+
1
)
α
∑
|
n
|
=
n
1
+
·
·
·
+
n
d
=
m
‍
n
i
!
×
[
(
1
-
1
n
i
+
1
)
×
∫
[
0
,
T
]
2
‍
β
n
(
s
1
)
β
n
(
t
2
)
·
R
H
,
K
i
(
s
1
,
t
2
)
n
i
-
1
×
∏
j
=
1
,
j
≠
i
d
‍
n
j
!
R
H
,
K
j
(
s
1
,
t
2
)
n
j
-
1
·
∂
∂
s
1
R
H
,
K
i
(
s
1
,
t
2
)
×
∂
∂
t
2
R
H
,
K
i
(
s
1
,
t
2
)
d
s
1
d
t
2
+
1
n
i
+
1
×
∫
[
0
,
T
]
2
‍
β
n
(
s
1
)
β
n
(
t
1
)
R
H
,
K
i
(
s
1
,
t
1
)
n
i
-
1
·
∏
j
=
1
,
j
≠
i
d
‍
n
j
!
R
H
,
K
j
(
s
1
,
t
1
)
n
j
×
∂
2
∂
s
1
∂
t
1
R
H
,
K
i
(
s
1
,
t
1
)
d
s
1
d
t
1
]
.
According to (4.37) in [6], for β
∈
[
1
/
4
,
1
/
2
), it follows that
(60)
∏
j
=
1
d
‍
β
n
(
u
)
≤
C
3,10
∏
j
=
1
d
‍
1
n
j
!
(
n
j
∨
1
)
(
(
8
β
-
1
)
/
12
)
,
where C
3,10 is a constant.
On the other hand, as in [11] and using the inequality a
2
+
b
2
≤
2
a
b, for a
,
b
∈
ℝ
+, there exists a constant C
3,11
(
H
,
K
) depending on H and K such that
(61)
|
∂
R
H
,
K
i
(
s
,
t
)
∂
t
∂
R
H
,
K
i
(
s
,
t
)
∂
s
|
≤
C
3,11
(
H
,
K
)
(
s
t
)
2
H
i
K
i
-
1
,
|
∂
2
R
H
,
K
i
(
s
,
t
)
∂
s
∂
t
|
≤
C
3,11
(
H
,
K
)
(
s
t
)
H
i
K
i
-
1
,
|
R
H
,
K
i
(
s
,
t
)
(
s
t
)
H
i
K
i
|
≤
C
3,11
(
H
,
K
)
.
Putting (60) and (61) into (59) and applying the self-similarity of the covariance R
(
s
,
t
), we have
(62)
∥
ξ
i
(
x
)
∥
2
,
α
-
1
2
≤
C
3,12
∑
m
≥
0
‍
(
1
+
m
)
α
×
∑
|
n
|
=
n
1
+
·
·
·
+
n
d
=
m
‍
∏
j
=
1
d
‍
1
(
n
j
∨
1
)
(
8
β
-
1
)
/
6
·
∫
[
0
,
T
]
2
‍
R
H
,
K
i
(
s
,
t
)
n
i
-
1
(
s
t
)
(
n
i
-
1
)
H
i
K
i
(
s
t
)
H
i
K
i
-
1
×
∏
j
=
1
,
j
≠
i
d
‍
R
H
,
K
j
(
s
,
t
)
n
j
(
s
t
)
n
j
H
j
K
j
d
s
d
t
≤
C
3,12
∑
m
≥
0
‍
(
1
+
m
)
α
×
∑
|
n
|
=
n
1
+
·
·
·
+
n
d
=
m
‍
∏
j
=
1
d
‍
1
(
n
j
∨
1
)
(
8
β
-
1
)
/
6
·
∫
0
T
‍
t
2
H
i
K
i
-
1
d
t
∫
0
1
‍
z
H
i
K
i
-
1
(
R
H
i
,
K
i
i
(
1
,
z
)
z
H
i
K
i
)
n
i
-
1
×
∏
j
=
1
,
j
≠
i
d
‍
R
H
,
K
j
(
1
,
z
)
n
j
z
H
j
K
j
d
z
,
where we use the change of variables t
=
t and z
=
s
/
r (see [6, 11]). According to (4.9) in [6], the bounded condition of integral part in (62) is proved.
Indeed,(63)
∫
0
1
‍
∏
j
=
1
d
‍
(
R
H
,
K
j
(
1
,
z
)
z
H
j
K
j
)
n
j
d
z
≤
C
3,13
m
-
1
/
2
(
H
K
)
*
,
where (
H
K
)
*
=
max
⁡
{
H
1
K
1
,
…
,
H
d
K
d
}.
Therefore, whenα
<
1
/
2
(
H
K
)
*
-
d
(
1
-
(
8
β
-
1
)
/
6
),
(64)
∥
ξ
i
(
x
)
∥
2
,
α
-
1
2
≤
C
3,14
∑
m
≥
1
‍
(
1
+
m
)
α
m
-
1
/
2
(
H
K
)
*
-
1
+
d
(
1
-
(
8
β
-
1
)
/
6
)
<
∞
.
Since β
∈
[
1
/
4
,
1
/
2
), we may choose β to tend to 1
/
2. Hence, when α
<
1
/
2
(
H
K
)
*
-
d
/
2, (64) is finite.



---

*Source: 762484-2014-04-02.xml*