Existence and stability results for neutral stochastic  
delay differential equations driven by a  
fractional Brownian motion  
1
2
R. Maheswari *, S. Karunanithi  
1
Department of Mathematics, Sri Eshwar College of Engg., Coimbatore-641 202, Tamilnadu,  
Department of Mathematics, Kongunadu Arts and Science College (Autonomous), Coimbatore-641 029, Tamilnadu, India  
Corresponding author E-mail:mahesenthil12@gmail.com  
2
*
Copyright © 2015R. Maheswari, S. Karunanithi. This is an open access article distributed under the Creative Commons Attribution License, which permits  
unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.  
Abstract  
In this paper we investigate the existence, uniqueness, asymptotic behavior of mild solutions to neutral stochastic  
differential equations with delays driven by a fractional Brownian motion in a Hilbert space. The cases of finite and  
infinite delays are analyzed.  
Keywords: Asymptotic Behaviors; Delays; Fractional Brownian Motion; Mild Solution; Wiener Integral.  
1
. Introduction  
The theory of stochastic differential equations driven by a fractional Brownian motion (fBm) has been studied  
intensively in recent years [1], [2], [3], [4], and [5]. The fBm received much attention because of its huge range of  
potential applications in several fields like telecommunications, networks, finance markets, biology and so on [6], [7],  
[
8], [9]. Moreover, one of the simplest stochastic processes that is Gaussian, self-similar, and has stationary increments  
is fBm [10]. In particular, fBm is a generalization of the classical Brownian motion, which depends on a parameter H ∈  
ꢀ
(
0, 1) called the Hurst index [9]. It should be mentioned that when H = , the stochastic process is a standard Brownian  
2
ꢀ
motion; when H ≠ 2, it behaves completely in a different way than the standard Brownian motion, in particular neither  
is a semimartingale nor a Markov process. It is a self-similar process with stationary increments and has a long-memory  
ꢀ
whenH ≠ 2. These significant properties make fBm a natural candidate as a model for noise in a wide variety of  
physical phenomena such as mathematical finance, communication networks, hydrology and medicine. The existence  
and uniqueness of mild solutions for a class of stochastic differential equations in a Hilbert space with a standard,  
ꢀ
cylindrical fBm with the Hurst parameter in the interval ꢁ , 1ꢂ has been studied [11]. Maslowski and Nualart [12] have  
2
studied the existence and uniqueness of a mild solution for nonlinear stochastic evolution equations in a Hilbert space  
driven by a cylindrical fBm under some regularity and boundedness conditions on the coefficients. Recently, Caraballo  
and et al [13] investigated the existence and uniqueness of mild solutions to stochastic delay equations driven by fBm  
ꢀ
with Hurst parameter H ∈ ꢁ2 , 1ꢂ. An existence and uniqueness result of mild solutions for a class of neutral stochastic  
differential equation with finite delay, driven by an fBm in a Hilbert space has been investigated [14] in Boufoussi and  
Hajji. The asymptotic behavior of solutions for stochastic differential equations with fBm has only been investigated by  
a few authors [13], [14], [15]. Moreover Nguyen [16] has studied the asymptotic behaviors of mild solutions to neutral  
stochastic differential equations driven by an fBm. Motivated by this consideration in this paper we investigate the  
existence and uniqueness and asymptotic behaviors of mild solutions for a neutral stochastic differential equation with  
finite or infinite delays driven by fBm in the following form  
d [x(t) − g ꢁt, xꢃt − r(t)ꢄꢂ] = [Ax(t) + f ꢁt, xꢃt − ρ(t)ꢄꢂ] dt + h ꢁt, xꢃt − δ(t)ꢄꢂ dW(t)  
ꢅ
+
σ(t)dB (t) ; t ≥ 0  
Q
x(t) = ϕ(t), t ∈ (−τ, 0ꢆ(0 < ꢇ ≤ ∞)  
(1.1)  
where A is an infinitesimal generator of an analytic semigroup of bounded linear operators,ꢃS(t)ꢄꢈꢉꢊ in a Hilbert  
ꢀ
ꢅ
space X with norm‖∙‖, B (t)denotes an fBm with H > on a real and separable Hilbert space Y, r, ρ, δ: ꢋ0, ∞) → ꢋ0, τ)  
Q
2
ꢊ
ꢊ
are continuous, f, g: ꢋ0, ∞) × X → X , h: ꢋ0, ∞) × X → ꢌ , σ: ꢋ0, ∞) → ℒ (Y, X) , the initial data ϕ ∈ Cꢃ(−τ, 0ꢆ, Xꢄ the  
2
Q
space of all continuous functions from(−τ, 0ꢆ to X and has finite second moments. Further we assume W and Bꢅ are  
independent. The main tool of this paper is the fixed point theory which was proposed by Burton [17].  
2
. Preliminaries  
In this section we first recall the fBm as well as the Wiener integral with respect to it. We also establish some  
important results which will be needed throughout the paper. Let (Ω, ℱ, P) be a complete probability space and T > 0  
be an arbitrary fixed horizon. Let X and be two real, separable Hilbert spaces and L(Y, X) be the space of bounded  
linear operators from Y to X.  
ꢊ
ꢊ
Let ꢌ (Y, X) be separable Hilbert space with respect to the Hilbert-Schmidt norm‖∙‖ꢍ  
ꢏ
. Let ℒ (Y, X) be the space of all  
2
Q
ꢎ
ꢑ
ꢑ
ψ ∈ ℒ(Y, X) such that ψꢐ  
ꢎ
is a Hilbert Schmidt operator. The norm is given be‖ψ‖2ℒ  
ꢏ
= ꢓψꢐ  
ꢎ
ꢓ = tr(ψꢐψ∗). Then ψ  
ꢒ
ꢊ
is called a ꢐ Hilbert –Schmidt operator from Y to X.Let ℒ (Ω, X) denote the space of ℱꢊ- measurable, X -valued and  
2
square integrable Stochastic processes. Consider a time interval ꢋ0, Tꢆ with arbitrary fixed horizon T and let  
ꢀ
{
βꢅ(t), t ∈ ꢋ0, Tꢆ} be a two sided one dimensional fBm with Hurst parameter H ∈ ꢁ , 1ꢂ with the covariance function  
2
ꢀ
Rꢅ(t, s) = Eꢋβꢅ(t)βꢅ(s)ꢆ = 2 (|t|2ꢅ + |s|2ꢅ − |t − s|2ꢅ), t, s ∈ ℝ.  
ꢀ
It is known that βꢅ(t) with H > 2 admits the following Wiener integral representation  
ꢈ
βꢅ(t) = ∫ Kꢅ(t, s) dW(s),  
ꢊ
where W = {W(t): t ∈ ꢋ0, Tꢆ} a Wiener is process and Kꢅ(t, s) is the kernel given by  
ꢅ
ꢔꢀ⁄2  
du, t > ꢗ.  
Kꢅ(t, s) = cꢅ ∫ (u − s)ꢅꢔ3⁄  
ꢈ
ꢖ
2
ꢁꢕꢂ  
ꢕ
ꢅ
(2ꢅꢔꢀ)  
ꢑ
ꢎ
with cꢅ = √ꢘ  
with β(∙) representing the beta function. Let Kꢅ(t, s) = 0 if t ≤ s .  
ꢁ2ꢔ2ꢅ ,ꢅꢔ ꢂ  
We will denote by ℋthe reproducing kernel Hilbert space of the fBm. Infact ℋis the closure of set of indicator  
functions ꢙ1ꢋꢊ;ꢈꢆ, t ∈ ꢋ0, Tꢆꢚ with respect to the scalar product  
〈
1ꢋꢊ,ꢈꢆ, 1ꢋꢊ,ꢕꢆ〉ℋ = Rꢅ (t, s).  
The mapping 1ꢋꢊ,ꢈꢆ → βꢅ(t) can be extended to an isometry between ℋ and the first Wiener chaos and we will denote  
ꢅ
by β (φ) by the previous isometry. We recall for ψ, φ ∈ ℋ their scalar product in ℋ is given by  
ꢝ ꢝ  
ψ, φ〉ℋ = H(ꢛH − 1) ꢜ ꢜ ψ(s)φ(t) |t − s|2ꢅꢔ2dsdt  
ꢊ ꢊ  
〈
For the deterministic function φ ∈ L2(ꢋ0, Tꢆ), the fractional Weiner integral of φ with respect to βꢅ is defined by  
ꢝ
ꢝ
φ(s)dβꢅ(t) = ∫ K φ(s)dβꢅ(s) ,  
∗
ꢅ
ꢊ
∫
ꢊ
ꢝ
∂k  
∗
where K φ(s) = ∫ φ(r) ∂ꢞ (r, s)dr  
ꢅ
ꢕ
Letꢙβnꢅ(t)ꢚn∈ℕbe a sequence of two-sided one dimensional standard fBm mutually independent on (Ω, ℱ, P). When one  
consider the following series  
∞
ꢅ
n
∑
β (t)en,t ≥ 0  
nꢟꢀ  
where {en} is a orthonormal basis in Y, this series does not necessarily converge in the space Y.Thus we consider a  
ꢅ
Y - valued stochastic process B (t) given formally by the following series  
Q
ꢑ
ꢅ
∞
ꢅ
n
B (t) = ∑  
β (t) ꢐ  
ꢎ
en,, t ≥ 0  
Q
nꢟꢀ  
If ꢐ is a non-negative self-adjoint trace class operator, then this series converges in the space Y, that is, it holds that  
ꢅ
2
ꢅ
Q
B (t) ∈ L (Ω, Y) then we say that the above B (t) is a Y - valued ꢐ - cylindrical fBm with covariance operator ꢐ. For  
Q
example if {λn} is a bounded sequence of non-negative real numbers such that ꢐen = λnen assuming that ꢐ is a nuclear  
operator in Y (ie.,∑n∞ꢟꢀ λn < ∞) then the stochastic process  
ꢑ
ꢎ
ꢅ
ꢅ
∞
ꢅ
n
en = ∑∞nꢟꢀ ꢠλnβn (t) en t ≥ 0  
B (t) = ∑  
β (t) ꢐ  
Q
nꢟꢀ  
is well defined as a Y -valued ꢐ - cylindrical fBm.  
ꢅ
n
where β (t) are real, independent fBm’s. This process is a Y –valued Gaussian, it starts from 0, zero mean and  
covariance  
E〈βꢅ(t), x〉〈βꢅ(s), y〉 = R(s, t)〈ꢐ(x), y〉For all x, y ∈ Y and t, s ∈ ꢋ0, Tꢆ.  
ꢓK φꢐꢀ⁄2en,ꢓ (ꢋꢊ,ꢝꢆ:ꢡ) < ∞  
ꢊ
∞
∗
Definition 2.1: Let φ: ꢋ0, Tꢆ → ℒ (Y, X) such that ∑  
(2.1)  
Q
nꢟꢀ  
ꢅ
ꢍꢎ  
ꢅ
Then its stochastic integral with respect to the fBm B is defined, for t ≥ 0, as follows  
Q
ꢈ
ꢈ
ꢅ
n
ꢈ
ꢅ
∞
∞
∗
ꢅ
∫
φ(s) dB (s) ≔ ∑  
∫ φ(s)ꢠλnend β (s) = ∑  
∫ ꢠλnꢢ (φen) (s)dβn(s)  
Q
nꢟꢀ  
nꢟꢀ  
ꢊ
ꢊ
ꢊ
Notice that if  
∞
∑
ꢣφꢠλnenꢣꢍ  
ꢑ
ꢤ (ꢋꢊ,ꢝꢆ:ꢅ)  
< ∞  
(2.2)  
nꢟꢀ  
ꢈ
ꢊ
2
ℒ
Lemma 2.2: For any φ: ꢋ0, Tꢆ → ℒ (Y, X) satisfies ∫ ‖φ(s)‖  
ꢏ
ds < ∞ such that (2.2) holds, and for any α, β ∈ ꢋ0, Tꢆ  
Q
ꢊ
ꢒ
with α > ꢥ,  
2
2
α
α
ꢀ⁄  
E ꢓ∫β φ(s)dB (s)ꢓ ≤ cH(ꢛH − 1)(α − β)2ꢅꢔꢀ ∑n∞ꢟꢀ ∫ ꢓφ(s)ꢐ 2enꢓ ds  
ꢅ
Q
β
ꢀ⁄  
where c = c(H). If in addition, ∑n∞ꢟꢀ ꢓφ(t)ꢐ 2enꢓ is uniformly convergent for t ∈ ꢋ0, Tꢆ. Then  
2
α
α
E ꢓ∫β φ(s)dB (s)ꢓ ≤ cH(ꢛH − 1)(α − β)2ꢅꢔꢀ ∫ ‖φ(s)‖  
ꢅ
2
ℒ
ꢏ
ds  
Q
β
ꢒ
We now suppose that 0 ∈ ρ(A), where ρ(A) is the resolvent set of A, and the semigroup S(t) is uniformly bounded that  
is to say‖S(t)‖ ≤ M for some constant M ≥ 1 and every t ≥ 0. Then, for 0 < ꢦ ≤ 1, it is possible to define the  
α
α
α
fractional power operator (−A) as a closed linear operator on its domain D(−A) . Furthermore the subspace D(−A) is  
α
α
α
dense in X and the expression ‖X‖α = ‖(−A) x‖ , x ∈ D(−A) defines a norm on Xα = D(−A) . The following  
properties are well known in [18].  
Lemma 2.3:[18] Under the above conditions the following properties hold.  
i)  
Xα is a Banach space for 0 < ꢦ ≤ 1  
ii)  
iii) For every 0 < ꢦ ≤ 1, there exists Mα such that‖(−A) S(t)‖ ≤ Mαt  
If the resolvent operator of A is compact, then embedding Xβ ⊂ Xα is continuous and compact for 0 < ꢦ ≤ ꢥ.  
e , λ > 0, ꢧ ≥ 0.  
α
ꢔα ꢔλꢈ  
ꢊ
Lemma 2.4:[19] For any r ≥1 and for arbitrary L valued predictable process φ(∙),  
2
ꢑ
ꢩ
ꢂ dsꢪ  
ꢞ
2
ꢞ
ꢞ
ꢈ
ꢕꢖp  
ꢕ
2ꢞ  
ꢍꢎ  
Eꢣ∫ꢊ φ(u) dw(u)ꢣ ≤ ꢃr(ꢛr − 1)ꢄ ꢨ∫ꢊ ꢁE‖φ(s)‖  
ꢏ
ꢕ
∈ꢋꢊ,ꢈ ꢆ  
ꢡ
3
. Existence and uniqueness of a solution  
In this section we study the existence and uniqueness of mild solutions for equation (1.1). For this equation we assume  
that the following conditions hold.  
(
Hꢀ)A is the infinitesimal generator of an analytic semigroup,ꢃS(t)ꢄꢈꢉꢊ, of bounded linear  
operators on X. Further we suppose that 0 ∈ ρ(A) and that ‖S(t)‖ ≤ M and  
ꢫꢑꢬβ  
(−A)ꢀꢔβS(t)ꢣ ≤ ꢈꢑꢬβ  
ꢣ
.
For some constants M, Mꢀꢔβ and every t ∈ ꢋ0, Tꢆ  
(
H2) There exists a positive constants Cꢀ, C2 and Lꢀ ,L2 > 0 such that, for all t ≥ 0, x, y ∈ X  
‖f(t, x) − f(t, y)‖2 ≤ Lꢀ ‖x − y‖2  
2 2  
‖f(t, x)‖ ≤ Cꢀ(1 + ‖x‖ )  
2
i)  
ii)  
iii) ‖h(t, x) − h(t, y)‖ ≤ L2 ‖x − y‖2  
2
2
iv) ‖h(t, x)‖ ≤ C2(1 + ‖x‖ )  
(
H3)There exists a constant 0 < ꢥ < 1 , L3 > 0 such that the function g is Xβ - valued and  
satisfies for all t ≥ 0,x, y ∈ X such that  
2
ꢣ(−A)βg(t, x) − (−A)βg(t, y)ꢣ ≤ L3 ‖x − y‖2  
i)  
2
ꢣ(−A)βg(t, x)ꢣ ≤ C3(1 + ‖x‖2)  
iii) The constants L3 and ꢥ, ꢣ(−A)ꢔβꢣL3 < 1  
ii)  
H4) The function (−A)βg is continuous in the quadratic mean sense square. For all functions X  
(
2
limꢈ→∞ꢣ(−A)βgꢃt, x(t)ꢄ − (−A)βgꢃs, x(s)ꢄꢣ = 0  
ꢝ
ꢊ
2
ℒ
(
H5) The function σ: ꢋ0, ∞) → ℒ (Y, X) satisfies∫ꢊ ‖σ(s)‖  
ꢏ
ds < ∞,for all T > 0.  
Q
ꢒ
Moreover we assume that ϕ ∈ ꢭꢃꢋ−τ, 0 ꢆ, ꢌ2(Ω, X)ꢄ.  
Definition 3.1: A X valued stochastic process {x(t), t ∈ (−τ, T)} is called a mild solution of Eqn. (1.1) if  
x(∙) ∈ ꢭꢃꢋ−τ, 0 ꢆ, ꢌ2(Ω, X)ꢄ  
x(t) = ϕ(t) for −τ ≤ t ≤ 0  
i)  
ii)  
iii) For arbitrary t ∈ ꢋ0, Tꢆ we have  
x(t) = S(t) ꢮφ(0) + g ꢁ0, φꢃ−r(0)ꢄꢂꢯ − g ꢁt, xꢃt − r(t)ꢄꢂ  
ꢈ ꢈ  
∫ AS(t − s)g ꢁs, xꢃs − r(s)ꢄꢂ ds + ∫ S(t − s) f ꢁs, xꢃs − ρ(s)ꢄꢂ ds  
ꢊ ꢊ  
−
ꢈ ꢈ  
∫ S(t − s) h ꢁs, xꢃs − δ(s)ꢄꢂ dW(s) + ∫ S(t − s) σ(s)dBꢅ(s)  
ꢊ ꢊ  
+
Theorem 3.1: Suppose(Hꢀ) to (H5) hold. Then for all T > 0, then Eqn. (1.1) has a unique mild solution on ꢋ−τ, Tꢆ.  
Proof: Fix T > 0 and let ℬ ≔ ꢭꢋ−τ, Tꢆ,ꢌ2(Ω, X),be the Banach space of all continuous functions from ꢰ– τ, Tꢱ into  
ꢀ
ꢌ2  
x ∈ ℬꢝ: x(s) = ϕ(s) for s ∈ ꢋ−τ, 0ꢆ}.  
Sꢝ is a closed subset of ℬꢝ provided with the norm ‖∙‖ℬꢲ . Define the operator ψ on Sꢝ by ψ(x)(t) = ϕ(t) for t ∈  
−τ, 0ꢆ and for t ∈ ꢋ0, Tꢆ  
(Ω, X) , equipped with the supremum norm ‖ξ‖ℬꢲ = suꢳꢖ∈ꢋꢔτ,ꢝꢆ(E‖ξ(u)‖2) ⁄2 and let us consider the set Sꢝ  
=
{
ꢋ
ψ(x)(t) = S(t) ꢮϕ(0) + g ꢁ0, ϕꢃ−r(0)ꢄꢂꢯ − g ꢁt, xꢃt − r(t)ꢄꢂ  
ꢈ ꢈ  
∫ AS(t − s)g ꢁs, xꢃs − r(s)ꢄꢂ ds + ∫ S(t − s) f ꢁs, xꢃs − ρ(s)ꢄꢂ ds  
ꢊ ꢊ  
−
+
ꢈ
ꢈ
∫ S(t − s) h ꢁs, xꢃs − δ(s)ꢄꢂ dW(s) + ∫ S(t − s) σ(s)dBꢅ(s)  
ꢊ
ꢊ
Then it is clear to prove the existence of mild solutions to Eqn. (1.1) is equivalent to find a fixed point for the operator ψ.  
Next we will show by Banach fixed point theorem that ψ has a unique fixed point. We divide the proof into two steps.  
2
Step 1. For arbitrary x ∈ Sꢝ , let us prove that t → ψ(x)(t) is continuous on the interval ꢋ0, Tꢆ in the ꢌ (Ω, X) sense .Let  
0
< ꢧ < ꢴ and |h| be sufficiently small. Then for any x ∈ Sꢝ, we have  
‖
ψ(x)(t + h) − ψ(x)(t)‖ ≤ ꢓꢃS(t + h) − S(t)ꢄ ꢮϕ(0) + g ꢁ0, ϕꢃ−r(0)ꢄꢂꢯꢓ  
+
+
+
+
+
ꢓg ꢁt + h, xꢃt + h − r(t + h)ꢄꢂ − g ꢁt, xꢃt − r(t)ꢄꢂꢓ  
ꢈ
ꢵꢶ  
ꢈ
ꢓ∫ꢊ AS(t + h − s)g ꢁs, xꢃs − r(s)ꢄꢂ ds − ∫ AS(t − s)g ꢁs, xꢃs − r(s)ꢄꢂ dsꢓ  
ꢊ
ꢈ
ꢵꢶ  
ꢈ
ꢓ∫ꢊ S(t + h − s) f ꢁs, xꢃs − ρ(s)ꢄꢂ ds − ∫ S(t − s) f ꢁs, xꢃs − ρ(s)ꢄꢂ dsꢓ  
ꢊ
ꢈ
ꢵꢶ  
ꢈ
ꢓ∫ꢊ S(t + h − s) h ꢁs, xꢃs − δ(s)ꢄꢂ dw(s) − ∫ S(t − s) h ꢁs, xꢃs − δ(s)ꢄꢂ dw(s)ꢓ  
ꢊ
ꢈ
ꢵꢶ  
ꢈ
ꢓ∫ꢊ S(t + h − s) σ(s)dBꢅ(s) − ∫ S(t − s) σ(s)dBꢅ(s)ꢓ  
ꢊ
=
∑ꢀꢸꢷꢸ6 Iꢷ(h)  
By the strong continuity of S(t), we have  
limꢶ→ꢊ ꢓꢃS(t + h) − S(t)ꢄ ꢮϕ(0) + g ꢁ0, ϕꢃ−r(0)ꢄꢂꢯꢓ = 0  
The condition (Hꢀ) assure that  
ꢃS(t + h) − S(t)ꢄ ꢮϕ(0) + g ꢁ0, ϕꢃ−r(0)ꢄꢂꢯꢓ ≤ ꢛM ꢓϕ(0) + g ꢁ0, ϕꢃ−r(0)ꢄꢂꢓ ∈ ꢌ2(Ω)  
ꢓ
Then we conclude by the Lebesgue dominated theorem that limꢶ→ꢊ ꢹ|Iꢀ(h)| = 0  
ꢔ
β
By using the fact that the operator (−A) is bounded, we obtained that  
2
2
|I2(h)|2 ≤ ꢣ(−A)ꢔβꢣ ꢹ ꢓ(−A)βg ꢁt + h, xꢃt + h − r(t + h)ꢄꢂ − (−A)βg ꢁt, xꢃt − r(t)ꢄꢂꢓ  
ꢹ
Then we conclude by condition (H4) that limꢶ→ꢊ ꢹ|I2(h)| = 0  
For I3(h), we suppose that h > 0, then we have  
ꢈ
I3(h) ≤ ꢓ∫ꢊ Sꢃ(h) − Iꢄ(−A)ꢀꢔβS(t − s)(−A)βg ꢁs, xꢃs − r(s)ꢄꢂ dsꢓ  
ꢈꢵꢶ  
ꢓ∫ꢈ (−A)ꢀꢔβS(t + h − s)(−A)βg ꢁs, xꢃs − r(s)ꢄꢂ dsꢓ  
+
By Holder’s inequality  
2
ꢈ
ꢓ∫ꢊ Sꢃ(h) − Iꢄ(−A)ꢀꢔβS(t − s)(−A)βg ꢁs, xꢃs − r(s)ꢄꢂ dsꢓ  
ꢹ
ꢈ
2
t ꢹ ꢜ ꢓSꢃ(h) − Iꢄ(−A)ꢀꢔβS(t − s)(−A)βg ꢁs, xꢃs − r(s)ꢄꢂꢓ ds  
≤
ꢊ
By using the strong continuity of S(t), we have for each s ∈ ꢋ0, tꢆ,  
ꢶli→mꢊSꢃ(h) − Iꢄ(−A)ꢀꢔβS(t − s)(−A)βg ꢁs, xꢃs − r(s)ꢄꢂ = 0  
By using condition (Hꢀ), condition (ii) in (H3) and the fact that 0 < ꢥ < 1 , we obtain  
Sꢃ(h) − Iꢄ(−A)ꢀꢔβS(t − s)(−A)βg ꢁs, xꢃs − r(s)ꢄꢂꢓ  
ꢓ
(
ꢫꢵꢀ)ꢫꢑꢬβ  
β
ꢓꢃ– Aꢄ g ꢁs, xꢃs − r(s)ꢄꢂꢓ ∈ ꢌ2(ꢋ0, tꢆ × Ω)  
≤
(
ꢈꢔꢕ)ꢑꢬ  
β
Then we conclude by the Lebesgue dominated theorem that  
2
ꢈ
limꢶ→ꢊ ꢹ ꢓ∫ꢊ Sꢃ(h) − Iꢄ(−A)ꢀꢔβS(t − s)(−A)βg ꢁs, xꢃs − r(s)ꢄꢂ dsꢓ = 0  
By conditions (Hꢀ), condition (ii)in (H3) and Holder’s inequality we get  
2
ꢈꢵꢶ  
ꢓ∫ꢈ (−A)ꢀꢔβS(t + h − s)(−A)βg ꢁs, xꢃs − r(s)ꢄꢂ dsꢓ  
ꢹ
ꢫꢎꢑꢬβ  
≤ 2βꢔꢀ  
ꢝ
h2βꢔꢀ ∫ C (ꢹ‖x(s) − r(s)‖2 + 1) ds  
2
3
ꢊ
2
ꢈ
ꢵꢶ  
Then limꢶ→ꢊ ꢓ∫ꢈ (−A)ꢀꢔβS(t + h − s)(−A)βg ꢁs, xꢃs − r(s)ꢄꢂ dsꢓ = 0  
Hence that limꢶ→ꢊ ꢹ|I3(h)| = 0  
ꢈ
I4(h) ≤ ꢓ∫ꢊ Sꢃ(h) − IꢄS(t − s) f ꢁs, xꢃs − ρ(s)ꢄꢂ dsꢓ  
ꢈ
ꢵꢶ  
+
ꢓ∫ꢈ S(t + h − s) f ꢁs, xꢃs − ρ(s)ꢄꢂ dsꢓ  
2
Now ꢹ ꢓ∫ꢊ Sꢃ(h) − IꢄS(t − s) f ꢁs, xꢃs − ρ(s)ꢄꢂ dsꢓ  
ꢈ
ꢈ
2
2
≤
∫ Sꢃ(h) − IꢄS(t − s) C ꢁꢹꢣxꢃs − ρ(s)ꢄꢣ + 1ꢂ ds  
ꢀ
ꢊ
ꢝ
2
2
M ∫ Sꢃ(h) − Iꢄ C ꢹ ꢁꢣxꢃs − ρ(s)ꢄꢣ + 1ꢂ ds  
ꢀ
ꢊ
≤
2
Since ꢶli→mꢊꢹ ꢓ∫ꢊ Sꢃ(h) − IꢄS(t − s) f ꢁs, xꢃs − ρ(s)ꢄꢂ dsꢓ = 0  
ꢈ
2
ꢈ
ꢵꢶ  
and ꢹ ꢓ∫ꢈ S(t + h − s) f ꢁs, xꢃs − ρ(s)ꢄꢂ dsꢓ  
2
ꢈꢵꢶ  
M2ꢹ ꢓ∫ꢈ f ꢁs, xꢃs − ρ(s)ꢄꢂ dsꢓ  
≤
ꢈꢵꢶ  
2
M2 ∫ C ꢹ ꢁꢣxꢃs − ρ(s)ꢄꢣ + 1ꢂ ds  
2
ꢀ
ꢈ
≤
2
Then ꢶli→mꢊꢹ ꢓ∫ꢈ S(t + h − s) f ꢁs, xꢃs − ρ(s)ꢄꢂ dsꢓ = 0  
ꢈ
ꢵꢶ  
Hence limꢶ→ꢊ ꢹ|I4(h)| = 0  
ꢈ
Now I5(h) ≤ ꢓ∫ꢊ Sꢃ(h) − IꢄS(t − s) h ꢁs, xꢃs − δ(s)ꢄꢂ dW(s)ꢓ  
ꢈ
ꢵꢶ  
+
ꢓ∫ꢈ S(t + h − s) h ꢁs, xꢃs − δ(s)ꢄꢂ dW(s)ꢓ  
2
Now ꢹ ꢓ∫ꢊ Sꢃ(h) − IꢄS(t − s) h ꢁs, xꢃs − δ(s)ꢄꢂ dW(s)ꢓ  
ꢈ
ꢈ
2
2
≤
∫ Sꢃ(h) − IꢄS(t − s) C ꢁꢹꢣxꢃs − δ(s)ꢄꢣ + 1ꢂ ds  
2
ꢊ
ꢈ
2
M2 ∫ Sꢃ(h) − Iꢄ C ꢁꢹꢣxꢃs − δ(s)ꢄꢣ + 1ꢂ ds  
2
2
ꢊ
≤
2
Therefore ꢶli→mꢊꢹ ꢓ∫ꢊ Sꢃ(h) − IꢄS(t − s) h ꢁs, xꢃs − δ(s)ꢄꢂ dW(s)ꢓ = 0  
ꢈ
2
and ꢹ ꢓ∫ꢈ S(t + h − s) h ꢁs, xꢃs − δ(s)ꢄꢂ dW(s)ꢓ  
ꢈ
ꢵꢶ  
ꢈ
ꢵꢶ  
2
M2 ∫ C ꢹ ꢁꢣxꢃs − δ(s)ꢄꢣ + 1ꢂ ds  
2
≤
2
ꢈ
2
Then ꢶli→mꢊꢹ ꢓ∫ꢈ S(t + h − s) h ꢁs, xꢃs − δ(s)ꢄꢂ dW(s)ꢓ = 0  
ꢈ
ꢵꢶ  
Hence limꢶ→ꢊ ꢹ|I5(h)| = 0  
For the term I6(h), we have  
ꢈ
ꢈꢵꢶ  
I6(h) ≤ ꢓ∫ꢊ Sꢃ(h) − IꢄS(t − s) σ(s)dBꢅ(s)ꢓ + ꢓ∫ꢈ S(t + h − s) σ(s)dBꢅ(s)ꢓ  
By condition (Hꢀ) and lemma 2.2 we get  
2
ꢈ
ꢓ∫ꢊ Sꢃ(h) − IꢄS(t − s) σ(s)dBꢅ(s)ꢓ  
ꢹ
ꢈ
2
ℒ
cH(ꢛH − 1)(t)2ꢅꢔꢀ ∫ ꢣSꢃ(h) − IꢄS(t − s)σ(s)ꢣ  
≤
ꢏ
ds  
ꢊ
ꢒ
ꢝ
2
ℒ
cH(ꢛH − 1)(t)2ꢅꢔꢀM ∫ ꢣSꢃ(h) − Iꢄσ(s)ꢣ  
≤
ꢏ
ds  
ꢊ
ꢒ
2
Since ꢶli→mꢊꢣSꢃ(h) − Iꢄσ(s)ꢣ = 0 and  
ℒꢒꢏ  
2
Sꢃ(h) − Iꢄσ(s)ꢣℒ  
≤ ꢺM2‖σ(s)‖2  
∈ ℒꢀ(ꢋ0, Tꢆ, ds)  
ꢣ
ꢏ
ꢏ
ℒꢒ  
ꢒ
we conclude by the dominated convergence theorem that  
2
ꢈ
ꢶli→mꢊꢹ ꢓ∫ꢊ Sꢃ(h) − IꢄS(t − s) σ(s)dBꢅ(s)ꢓ = 0  
Applying lemma 2.2, we get  
ꢈ
ꢵꢶ  
ꢈꢵꢶ  
ꢓ∫ꢈ S(t + h − s) σ(s)dBꢅ(s)ꢓ ≤ cH(h)2ꢅꢔꢀM2 ∫ꢈ ‖σ(s)‖ℒ2  
ꢹ
ꢏ
ds → 0  
ꢒ
Hence that limꢶ→ꢊ ꢹ|I6(h)| = 0  
The above arguments show that ꢶli→mꢊꢹ‖ψ(x)(t + h) − ψ(x)(t)‖2 = 0. Hence we conclude that the function t → ψ(x)(t)  
2
is continuous on ꢋ0, Tꢆ in the ꢌ - sense.  
Step 2. Now we are going to show that ψ is a contraction mapping in Sꢝꢑ with some Tꢀ ≤ T.  
ꢀ
3
3
3
Let x, y ∈ S ꢝ by using the inequality (a + b + c + d)2 ≤ a2 +  
b2 +  
c2 +  
d2  
k
where ꢢ = L3ꢣ(−A)ꢔβꢣ < 1, we obtain any fixed t ∈ ꢋ0, Tꢆ  
ꢀꢔk  
ꢀꢔk  
ꢀꢔk  
2
ꢀ
‖ψ(x)(t) − ψ(y)(t)‖2 ≤ ꢓg ꢁt, xꢃt − r(t)ꢄꢂ − g ꢁt, yꢃt − r(t)ꢄꢂꢓ  
ꢹ
k
2
3
ꢈ
+ ꢀ  
ꢓ∫ꢊ AS(t − s) ꢮg ꢁs, xꢃs − r(s)ꢄꢂ − g ꢁs, yꢃs − r(s)ꢄꢂꢯ dsꢓ  
ꢔk  
2
ꢓ∫ꢊ S(t − s) ꢮf ꢁs, xꢃs − ρ(s)ꢄꢂ − f ꢁs, yꢃs − ρ(s)ꢄꢂꢯ dsꢓ  
3
ꢈ
+ ꢀ  
ꢔk  
2
ꢓ∫ꢊ S(t − s) ꢮh ꢁs, xꢃs − δ(s)ꢄꢂ − h ꢁs, yꢃs − δ(s)ꢄꢂꢯ dW(s)ꢓ  
3
ꢈ
+ ꢀ  
ꢔk  
2
2
ꢣ(−A)ꢔβꢣ ꢓ(−A)βg ꢁt, xꢃt − r(t)ꢄꢂ − (−A)βg ꢁt, yꢃt − r(t)ꢄꢂꢓ  
ꢀ
≤
k
2
3
ꢈ
ꢓ∫ꢊ(−A)ꢀꢔβS(t − s)(−A)β ꢮg ꢁs, xꢃs − r(s)ꢄꢂ − g ꢁs, yꢃs − r(s)ꢄꢂꢯꢓ  
+ ꢀ  
ꢔk  
2
3
ꢈ
ꢓ∫ꢊ S(t − s)(−A)β ꢮf ꢁs, xꢃs − ρ(s)ꢄꢂ − f ꢁs, yꢃs − ρ(s)ꢄꢂꢯ dsꢓ  
+ ꢀ  
+ ꢀ  
ꢔk  
2
3
ꢈ
ꢓ∫ꢊ S(t − s)(−A)β ꢮh ꢁs, xꢃs − δ(s)ꢄꢂ − h ꢁs, yꢃs − δ(s)ꢄꢂꢯ dW(s)ꢓ  
ꢔk  
By Lipchitz property of (−A)βg and f, h combined with Holder’s inequality, we obtain  
‖ψ(x)(t) − ψ(y)(t)‖2 ≤ ꢢꢹ‖x(t − r) − y(t − r)‖2  
ꢹ
3
ꢈꢎβꢬꢑ  
ꢈ
2
2
ꢂ ∫ꢊ ꢹ‖x(s − r) − y(s − r)‖2ds  
+ ꢀ  
C M  
ꢁ
3
ꢀꢔβ  
ꢔk  
2βꢔꢀ  
3
ꢈ
2
tM2C ∫ ꢹꢣxꢃs − ρ(s)ꢄ − yꢃs − ρ(s)ꢄꢣ ds  
2
+ ꢀ  
+ ꢀ  
ꢀ
ꢔk  
ꢊ
3
ꢈ
2
tM2C ∫ ꢹꢣxꢃs − δ(s)ꢄ − yꢃs − δ(s)ꢄꢣ ds  
2
2
ꢔk  
ꢊ
Hencesuꢳꢕ∈ꢋꢔτ,ꢈꢆ ꢹ‖ψ(x)(s) − ψ(y)(s)‖2 ≤ γ(t)ꢕ∈sꢋuꢔꢳτ,ꢈꢆꢹ‖x(s) − y(s)‖2  
ꢎ
ꢎ
3ꢫꢎ  
ꢎ
ꢈꢎ 3ꢫꢎ  
ꢎ
ꢈꢎ  
3
ꢻ
ꢫ
ꢻ
ꢻ
ꢼ
ꢑꢬβ  
ꢑ
ꢎ
t2β+  
+
where γ(t) = ꢢ + (  
ꢀꢔk)(2βꢔꢀ)  
ꢀꢔk  
ꢀꢔk  
By condition (iii)in (H3) we have,γ(0) = ꢢ = ꢣ(−A)ꢔβꢣL3 < 1.Then there exists 0 < ꢴꢀ ≤ T such that  
0
< ꢽ(Tꢀ) < 1 and ψ is a contraction mapping on Sꢝꢑ and therefore has a unique fixed point, which is a mild solution of  
Eqn.(1.1) on ꢰ– τ, Tꢀꢱ.This procedure can be repeated in order to extend the solution to the entire interval ꢰ– τ, Tꢱ in  
finitely many steps. This completes the proof.  
4
. Stability analysis  
In this section we establish the results for the case of finite delays, then the case of infinite delays can be proved. Our  
method is based on the contraction mapping principle.  
In order to prove the required results, we assume the following additional conditions.  
ꢔ
λꢈ  
(
(
H6) The semigroup S(t) satisfies ∃λ > 0, ∃M > 0such that ‖S(t)‖ ≤ Me , t ≥ 0  
∞
ꢊ
γꢕ  
2
H7) The function σ: ꢋ0, ∞) → ℒ satisfies∫ꢊ e ‖σ(s)‖ℒ  
ꢏ
ds < ∞ .  
Q
ꢒ
We first consider the case of finite delays, τ < ∞.  
Theorem 4.1: (Finite delays) Assume that f(t, 0) = g(t, 0) = h(t, 0) = 0, t ≥ 0, the assumptions (Hꢀ) - (H7)holds and  
that  
Γꢎ(β)  
λꢎβ  
2
ꢁL3ꢣ(−A)ꢔβꢣ + Mꢀ2ꢔβL3  
+ M2Lꢀλꢔ2+M2L2(ꢛλ)ꢔ2ꢂ < 1  
ꢺ
(4.1)  
where Γ(∙) is the gamma function,Mꢀꢔβ is the corresponding constant in lemma (2.3). Then the mild solution to (1.1)  
exists uniquely and is exponential decay to zero in mean square. i.e., there exists a pair of positive constants a > 0 and  
∗
∗
M = M (φ, a) such that  
‖x(t)‖2 ≤ M∗eꢔꢾꢈ, ∀ t ≥ 0.  
ꢹ
Proof: Denote by ꢿ the space of all stochastic process x(t, ω): (−τ, ∞) × Ω → X satisfying x(t) = ϕ(t), t ∈ (−τ, 0ꢆ and  
∗
∗
there exists some constants a > 0 and M = M (φ, a) > 0 such that  
‖x(t)‖2 ≤ M∗eꢔꢾꢈ, ∀ t ≥ 0.  
ꢹ
(4.2)  
It is to check that ꢿ is a Banach space endowed with a norm|x| = suꢳꢈꢉꢊ ꢹ|x(t)|2. Without loss of generality, we  
2
ꢿ
assume that a < ꣀ . We denote the operator Φ on ꢿ by  
(Φx)(t) = ϕ(t), t ∈ (−τ, 0ꢆ and  
(Φx)(t) = S(t) ꢮϕ(0) + g ꢁ0, ϕꢃ−r(0)ꢄꢂꢯ − g ꢁt, xꢃt − r(t)ꢄꢂ  
ꢈ
ꢈ
−
∫ AS(t − s)g ꢁs, xꢃs − r(s)ꢄꢂ ds + ∫ S(t − s) f ꢁs, xꢃs − ρ(s)ꢄꢂ ds  
ꢊ
ꢊ
ꢈ
ꢈ
∫ S(t − s) h ꢁs, xꢃs − δ(s)ꢄꢂ dW(s) + ∫ S(t − s) σ(s)dBꢅ(s)  
+
ꢊ
ꢊ
To get desired results, it is enough to show that operator Φ has a unique fixed point in ꢿ. For this purpose, we use the  
contraction mapping principle.  
Step 1. We first verify that Φ(ꢿ) ⊂ ꢿ. For convenience of notation, we denote by M ꢷ∗ , i = 1,ꢛ, … the finite positive  
constants depending on φ, a. By assumption (H6) we have  
2
‖ꢳꢀ(t)‖2 ≤ M2ꢹ ꢓϕ(0) + g ꢁ0, ϕꢃ−r(0)ꢄꢂꢓ eꢔλꢈ ≤ M eꢔλꢈ  
∗
ꢹ
(4.3)  
ꢀ
To estimate ꢳꢷ(t), i = ꢛ, … ,ꣁ. we observe that for x ∈ ꢿ and u(t) = r(t) orρ(t) or δ(t)the following useful estimate  
holds  
2 2  
ꢣxꢃt − u(t)ꢄꢣ ≤ ꢁM∗eꢔꢾꢃꢈꢔꢖ(ꢈ)ꢄ + ꢹꢣφꢃt − u(t)ꢄꢣ ꢂ  
ꢹ
ꢃM∗eꢔꢾꢃꢈꢔꢖ(ꢈ)ꢄ + ꢹE‖φ‖ꣂ2eꢔꢾꢃꢈꢔꢖ(ꢈ)ꢄꢄ  
≤
≤
(M∗ + ꢹ‖ϕ‖ꣂ2)eꢾτeꢔꢾꢈ  
where ‖ϕ‖ꣂ = suꢳꢔτꣃ꣄ꢸꢊ‖φ(s)‖ < ∞.  
Then by assumption (i) in (H3) we have  
2
2
‖ꢳ2(t)‖2 ≤ ꢣ(−A)ꢔβꢣ ꢹ ꢓ(−A)βg ꢁt, xꢃt − r(t)ꢄꢂ − (−A)β g(t, 0)ꢓ  
ꢹ
2
L3ꢣ(−A)ꢔβꢣ ꢹꢣxꢃt − r(t)ꢄꢣ2  
≤
2
L3ꢣ(−A)ꢔβꢣ (M∗ + ꢹ‖ϕ‖ꣂ2)eꢾτeꢔꢾꢈ  
M e  
2
≤
≤
∗
ꢔꢾꢈ  
(4.4)  
Using lemma 2.3, Holder’s inequality and assumption (i) in (H3) obtain that  
2
ꢈ
‖ꢳ3(t)‖2 ≤ ꢹ ꢓ∫ꢊ AS(t − s)g ꢁs, xꢃs − r(s)ꢄꢂ dsꢓ  
ꢹ
ꢈ
2
ꢣ(−A)ꢀꢔβS(t − s)ꢣds ꢜꢣ(−A)ꢀꢔβS(t − s)ꢣꢹ ꢓ(−A)βg ꢁs, xꢃs − r(s)ꢄꢂꢓ ds  
≤
ꢊ
ꢈ
ꢈ
2
2
ꢔβL3 ꢜ(t − s)βꢔꢀ eꢔλ(ꢈꢔꢕ)ds ꢜ(t − s)βꢔꢀ eꢔλ(ꢈꢔꢕ)ꢹꢣxꢃs − r(s)ꢄꢣ ds  
M
ꢀ
≤
ꢊ
ꢊ
Γ(β)  
λβ  
ꢈ
ꢀ2ꢔβL3  
∫ (t − s)βꢔꢀ eꢔλ(ꢈꢔꢕ)(M∗ + ꢹ‖φ‖ꣂ2)eꢾτeꢔꢾꢕds  
≤
≤
M
ꢊ
Γ(β)  
ꢈ
ꢀ2ꢔβL3 λβ (M∗ + ꢹ‖φ‖2ꣂ )eꢾτeꢔꢾꢈ ∫ (t − s)βꢔꢀ e(ꢾꢔλ)(ꢈꢔꢕ)ds  
M
ꢊ
Γꢎ(β)  
(M∗ + ꢹ‖ϕ‖ꣂ2)eꢾτeꢔꢾꢈ  
λβ(λꢔꢾ)β  
≤
M
ꢀ2ꢔβL3  
We therefore have ꢹ‖ꢳ3(t)‖2 ≤ M eꢔꢾꢈ  
∗
(4.5)  
3
Similarly, we obtain by assumption(i) iꣅ (H2) that  
2
ꢈ
‖ꢳ4(t)‖2 = ꢹ ꢓ∫ꢊ S(t − s) f ꢁs, xꢃs − ρ(s)ꢄꢂ dsꢓ  
ꢹ
ꢈ ꢈ  
2
M2Lꢀ ∫ eꢔλ(ꢈꢔꢕ)ds ∫ eꢔλ(ꢈꢔꢕ)ꢹꢣxꢃs − ρ(s)ꢄꢣ ds  
ꢊ ꢊ  
≤
M2Lꢀλꢔꢀ ∫ eꢔλ(ꢈꢔꢕ)(M∗ + ꢹ‖ϕ‖ꣂ2)eꢾτeꢔꢾꢕds  
ꢈ
≤
ꢊ
M2Lꢀλꢔꢀ(λ − a)ꢔꢀ(M∗ + ꢹ‖ϕ‖ꣂ2)eꢾτeꢔꢾꢈ  
≤
≤
M eꢔꢾꢈ  
‖ꢳ5(t)‖2 = ꢹ ꢓ∫ꢊ S(t − s) h ꢁs, xꢃs − δ(s)ꢄꢂ dw(s)ꢓ  
∗
(4.6)  
4
2
ꢈ
ꢹ
ꢈ
2
M2L2 ∫ eꢔ2λ(ꢈꢔꢕ)ꢹꢣxꢃs − δ(s)ꢄꢣ ds  
≤
≤
ꢊ
ꢈ
ꢈ
2
M2L2 ∫ eꢔ2λ(ꢈꢔꢕ)ds ∫ꢊ eꢔ2λ(ꢈꢔꢕ)ꢹꢣxꢃs − δ(s)ꢄꢣ ds  
ꢊ
ꢀ
M2L2λꢔꢀ ∫ eꢔ2λ(ꢈꢔꢕ)(M∗ + ꢹ‖ϕ‖ꣂ2)eꢾτeꢔꢾꢕds  
ꢈ
≤ 2  
ꢊ
ꢀ M2L2λꢔꢀ(ꢛλ − a)ꢔꢀ(M∗ + ꢹ‖ϕ‖ꣂ2)eꢾτeꢔꢾꢈ  
≤ 2  
M eꢔꢾꢈ  
∗
(4.7)  
≤
5
By using lemma 2.2 we get that  
2
ꢈ
‖ꢳ6(t)‖2 = ꢹ ꢓ∫ꢊ S(t − s) σ(s)dBꢅ(s)ꢓ  
ꢹ
ꢈ
M2cH(ꢛH − 1)t2ꢅꢔꢀ ∫ eꢔ2λ(ꢈꢔꢕ)‖σ(s)‖ꢍ2  
≤
ꢏ
ds  
ds  
ꢊ
ꢒ
From this inequality we can infer that  
‖ꢳ6(t)‖2 ≤ M2cH(ꢛH − 1)t2ꢅꢔꢀeꢔ2λ′ꢈ ∫ e2γꢕ‖σ(s)‖2ℒ  
ꢈ
ꢹ
ꢏ
(4.8)  
ꢊ
ꢒ
′
′
′
where λ = λ⋀γ.Indeed, if λ < ꢽ, then λ = λ and we have  
∞
‖ꢳ6(t)‖2 ≤ M2cH(ꢛH − 1)t2ꢅꢔꢀeꢔ2λꢈ ∫ e2λꢕ‖σ(s)‖ꢍ2  
ꢹ
ꢏ
ds  
ꢊ
ꢒ
M2cH(ꢛH − 1)t2ꢅꢔꢀeꢔ2λ′ꢈ ∫ e2γꢕ‖σ(s)‖ℒ2  
∞
≤
ꢏ
ds  
ꢊ
ꢒ
′
If γ < ꣀ then λ = γ and we have  
‖ꢳ6(t)‖2 ≤ M2cH(ꢛH − 1)t2ꢅꢔꢀeꢔ2λ′ꢈ ∫ e2γꢕ‖σ(s)‖2ℒ  
ꢈ
ꢹ
ꢏ
ds  
ꢊ
ꢒ
since suꢳꢈꢉꢊ ꢁt2ꢅꢔꢀ  
e
ꢔλ′ꢈꢂ < ∞, this together with (4.8) we have  
∗
‖ꢳ6(t)‖2 ≤ M eꢔꢾꢈ  
ꢹ
(4.9)  
6
Combining (4.3)-(4.7), (4.9) we see that there exists M̅ ∗ > 0 and a > 0 such that  
2
∗ ꢔꢾꢈ  
ꢹ
‖(Φx)(t)‖ ≤ M e , t ≥ 0.  
̅
Hence we conclude that Φ (ꢿ) ⊂ ꢿ.  
Step 2. We now show that Φ is a contraction mapping. For any x, y ∈ ꢿ, we have  
‖(Φx)(t) − (Φy)(t)‖2 ≤ ꢺ ∑nꢷꢟꢀ Iꢷ  
ꢹ
Since x(t) = y(t) = φ(t), t∈ (−τ, 0ꢆ, this implies that  
2
ꢣxꢃt − r(t)ꢄ − yꢃt − r(t)ꢄꢣ ≤ sꢈuꢉꢊꢳꢹ‖x(t) − y(t)‖2  
ꢹ
Then we have by assumption(i) in (H3)  
2
Iꢀ = ꢹ ꢓg ꢁt, xꢃt − r(t)ꢄꢂ − g ꢁt, yꢃt − r(t)ꢄꢂꢓ  
2
2
L3ꢣ(−A)ꢔβꢣ ꢹꢣxꢃt − r(t)ꢄ − yꢃt − r(t)ꢄꢣ  
≤
≤
2
L3ꢣ(−A)ꢔβꢣ sꢈuꢉꢊꢳꢹ‖x(t) − y(t)‖2  
2
I2 = ꢹ ꢓ∫ꢊ AS(t − s) [g ꢁs, xꢃs − r(s)ꢄꢂ − g ꢁs, yꢃs − r(s)ꢄꢂ] dsꢓ  
ꢈ
ꢈ ꢈ  
2
2
ꢔβL3 ∫ (t − s)βꢔꢀ eꢔλ(ꢈꢔꢕ)ds ∫ꢊ(t − s)βꢔꢀ eꢔλ(ꢈꢔꢕ)ꢹꢣxꢃs − r(s)ꢄ − yꢃs − r(s)ꢄꢣ ds  
M
ꢀ
ꢊ
≤
Γ(β)  
λβ  
ꢈ
2
M
ꢀ2ꢔβL3  
∫ (t − s)βꢔꢀ eꢔλ(ꢈꢔꢕ)ꢹꢣxꢃs − r(s)ꢄ − yꢃs − r(s)ꢄꢣ ds  
≤
ꢊ
Γꢎ(β)  
λꢎβ  
M
ꢀ2ꢔβL3  
sꢈuꢉꢊꢳꢹ‖x(t) − y(t)‖2  
≤
By assumption (i) iꣅ (H2)  
ꢈ
I3 = ꢹ ꢓ∫ꢊ S(t − s) [f ꢁs, xꢃs − ρ(s)ꢄꢂ − f ꢁs, yꢃs − ρ(s)ꢄꢂ] dsꢓ  
ꢈ
ꢈ
2
M2Lꢀ ∫ eꢔλ(ꢈꢔꢕ)ds ∫ eꢔλ(ꢈꢔꢕ)ꢹꢣxꢃs − ρ(s)ꢄ − yꢃs − ρ(s)ꢄꢣ ds  
≤
ꢊ
ꢊ
2
M2Lꢀλꢔꢀ ∫ eꢔλ(ꢈꢔꢕ)ꢹꢣxꢃs − ρ(s)ꢄ − yꢃs − ρ(s)ꢄꢣ ds  
ꢈ
≤
≤
ꢊ
M2Lꢀλꢔ2sꢈuꢉꢊꢳꢹ‖x(t) − y(t)‖2  
By assumption(iii) iꣅ (H2) , we have  
ꢈ
I4 = ꢹ ꢓ∫ꢊ S(t − s) [h ꢁs, xꢃs − δ(s)ꢄꢂ − h ꢁs, xꢃs − δ(s)ꢄꢂ] dw(s)ꢓ  
ꢈ
ꢈ
2
M2L2 ∫ eꢔ2λ(ꢈꢔꢕ)ds ∫ eꢔ2λ(ꢈꢔꢕ)ꢹꢣxꢃs − δ(s)ꢄ − yꢃs − δ(s)ꢄꢣ ds  
≤
ꢊ
ꢊ
ꢈ
2
M2L2(ꢛλ)ꢔꢀ ∫ eꢔ2λ(ꢈꢔꢕ)ꢹꢣxꢃs − δ(s)ꢄ − yꢃs − δ(s)ꢄꢣ ds  
≤
ꢊ
M2L2(ꢛλ)ꢔ2sꢈuꢉꢊꢳꢹ‖x(t) − y(t)‖2  
≤
Thusꢹ‖(Φx)(t) − (Φy)(t)‖2  
L3 Γ2(β)  
2
ꢺ ꢨL3ꢣ(−A)ꢔβꢣ + Mꢀ2ꢔβ  
+ M2Lꢀλꢔ2+M2L2(ꢛλ)ꢔ2ꢪ sꢈuꢉꢊꢳꢹ‖x(t) − y(t)‖2  
≤
λ2β  
By the condition (4.1), we claim that Φ is a contractive. So, applying the Banach fixed point principle, the proof is  
complete.  
Now we consider the case of infinite delays. We assume that t − r(t) → ∞, t − ρ(t) → ∞, t − δ(t) → ∞ as t → ∞.  
Theorem 4.2:[Infinite delays]Under the conditions of theorem 4.1, the mild solution to (1.1) exists uniquely and  
2
converges to zero in mean square i.e., limꢈ→∞ ꢹ‖x(t)‖ = 0  
Proof: Denote by ꢿ′ the space of all stochastic processes x(t, ω): (−∞, ∞) × Ω → X satisfying  
x(t) = ϕ(t), t ∈ (−∞, oꢆ and  
limꢈ→∞ꢹ‖x(t)‖2 = 0  
(4.10)  
We define the operator Ψ on ꢿ′ by (Ψx)(t) = ϕ(t), t ∈ (−∞, 0ꢆand  
(Ψx)(t) = S(t) ꢮϕ(0) + g ꢁ0, ϕꢃ−r(0)ꢄꢂꢯ − g ꢁt, xꢃt − r(t)ꢄꢂ  
ꢈ
ꢈ
−
∫ AS(t − s)g ꢁs, xꢃs − r(s)ꢄꢂ ds + ∫ S(t − s) f ꢁs, xꢃs − ρ(s)ꢄꢂ ds  
ꢊ
ꢊ
ꢈ
ꢈ
∫ S(t − s) h ꢁs, xꢃs − δ(s)ꢄꢂ dw(s) + ∫ S(t − s) σ(s)dBꢅ(s)  
+
ꢊ
ꢊ
Since (Ψx)(t) = (Φx)(t) on ꢋ0, ∞), this implies that Ψ is contractive. Hence it remains to check Φ(ꢿ′) ⊂ ꢿ′. In order to  
2
obtain this claim, we need to show thatlimꢈ→∞ꢹ‖Ψx(t)‖ = 0 for all  
′
x ∈ ꢿ .  
By the definition of ꢿ′, assumption (H7)and the fact t − r(t) → ∞, t → ∞. We get  
lim ꢹꢈ→∞‖Pꢀ(t)‖2 = lꢈi→m∞ꢹ‖P2(t)‖2 = lꢈi→m∞ꢹ‖P6(t)‖ = 0  
2
ꢈ
We further have ꢹ‖ꢳ3(t)‖2 ≤ ꢹ ꢓ∫ꢊ AS(t − s)g ꢁs, xꢃs − r(s)ꢄꢂ dsꢓ  
ꢈ
ꢈ
2
2
ꢔβL3 ∫ (t − s)βꢔꢀ eꢔλ(ꢈꢔꢕ)ds ∫ (t − s)βꢔꢀ eꢔλ(ꢈꢔꢕ)ꢹꢣxꢃs − r(s)ꢄꢣ ds  
≤
≤
M
M
ꢀ
ꢊ
ꢊ
ꢔβL3Γ(β)λꢔβ ∫ (t − s)βꢔꢀ eꢔλ(ꢈꢔꢕ)ꢹꢣxꢃs − r(s)ꢄꢣ ds  
ꢈ
2
2
ꢀ
ꢊ
2
For any x ∈ ꢿ′and ε > 0 it follows from (4.10) that there exists sꢀ > 0 such that ꢹꢣxꢃs − r(s)ꢄꢣ < ꣇ for all  
s ≥ sꢀ.Thus we obtain  
‖ꢳ3(t)‖2 ≤ M L3Γ(β)λꢔβ ∫ (t − s)βꢔꢀ eꢔλ(ꢈꢔꢕ)ꢹꢣxꢃs − r(s)ꢄꢣ ds + M L3Γ2(β)λꢔ2β  
ꢕ
ꢑ
2
2 2  
ϵ
ꢀꢔβ ꢀꢔβ  
ꢊ
ꢹ
which proves that  
limꢈ→∞ ꢹ‖P3(t)‖2 ≤ M L3Γ2(β)λꢔ2β,∀ϵ > 0 and hence, limꢈ→∞ ꢹ‖P3(t)‖2 = 0  
2
ꢀꢔβ  
2
ꢈ
‖ꢳ4(t)‖2 = ꢹ ꢓ∫ꢊ S(t − s) f ꢁs, xꢃs − ρ(s)ꢄꢂ dsꢓ  
ꢹ
ꢈ ꢈ  
2
M2Lꢀ ∫ eꢔλ(ꢈꢔꢕ)ds ∫ eꢔλ(ꢈꢔꢕ)ꢹꢣxꢃs − ρ(s)ꢄꢣ ds  
ꢊ ꢊ  
≤
M2Lꢀλꢔꢀ ∫ eꢔλ(ꢈꢔꢕ)ꢹꢣxꢃs − ρ(s)ꢄꢣ ds  
ꢈ
2
ꢊ
≤
2
For any x ∈ ꢿ′and ε > 0 there exists a sꢀ > 0 such that ꢹꢣxꢃs − r(s)ꢄꢣ < ꣇ for all s ≥ sꢀ.  
2
Thus we obtain ꢹ‖ꢳ4(t)‖2 ≤ M2Lꢀλꢔꢀ  
ꢕ
ꢑ
ꢔλ(ꢈꢔꢕ)  
ꢹꢣxꢃs − ρ(s)ꢄꢣ ds + M2Lꢀλꢔ2  
∫
e
ϵ
ꢊ
which proves that limꢈ→∞ E‖P4(t)‖2 ≤ M2Lꢀλꢔ2ϵ,∀ϵ > 0 and hence,imꢈ→∞ E‖P4(t)‖2 = 0  
2
ꢈ
‖ꢳ5(t)‖2 = ꢹ ꢓ∫ꢊ S(t − s) h ꢁs, xꢃs − δ(s)ꢄꢂ dw(s)ꢓ  
ꢹ
ꢈ
2
M2L2 ∫ eꢔ2λ(ꢈꢔꢕ)ꢹꢣxꢃs − δ(s)ꢄꢣ ds  
ꢊ
≤
2
For any x ∈ ꢿ′and ε > 0 there exists asꢀ > 0 such that ꢹꢣxꢃs − r(s)ꢄꢣ < ꣇for all s ≥ sꢀ.  
ꢕ
ꢑ
2
Thus we obtain ꢹ‖ꢳ5(t)‖2 ≤ M2L2 ∫  
ꢔ2λ(ꢈꢔꢕ)  
ꢹꢣxꢃs − δ(s)ꢄꢣ ds + M2L2(ꢛλ)ꢔꢀϵ  
e
ꢊ
2
which proves that limꢈ→∞ ꢹ‖P5(t)‖ ≤ M2L2(ꢛλ)ꢔꢀϵ, ∀ϵ > 0 and hence, limꢈ→∞ ꢹ‖P5(t)‖2 = 0  
Once again, by applying the Banach fixed point principle, we complete the proof of the theorem.  
To illustrate the obtained theory, let us end this section with an example.  
5
. Example  
Consider the following neutral stochastic partial differential equations with delays driven by a fBm in the following  
form  
꣈2  
꣈꣈t ꢰ  
ꢅ
z(t, ξ) + αꢀGꢃt, z(t − τ, ξ)ꢄꢱ = ꣈ξ2 z(t, ξ) + α2Fꢃt, z(t − τ, ξ)ꢄ ꣈t + α3φꢃt, z(t − τ, ξ)ꢄdw(t) + Θ(t)dB (t)  
Q
z(t, 0) = z(t, π) = 0  
(5.1)  
z(t, ξ) = ϕ(t, ξ), t ∈ (−∞, 0ꢆ, 0 ≤ ξ ≤ π  
where αꢀ, α2, α3 > 0 are constants.  
Let X = L (0, π) with the norm ‖∙‖ and inner product 〈∙,∙〉 .Define A: X → X by Ax = x" with domain  
2
′
"
D(A) = ꢙx ∈ X: x, x are absolutely coꣅtiꣅuous x ∈ X, x(0) = x(π) = 0ꢚ. Then Eqn. (5.1) can be written in the form  
of equation (1.1) with the co-efficient g(t, x) = αꢀGꢃt, z(t − τ, ξ)ꢄ, f(t, x) = α2Fꢃt, z(t − τ, ξ)ꢄ and  
h(t, x) = α3φꢃt, z(t − τ, ξ)ꢄ, σ(t) = Θ(t).  
For the operator A, it is known from Pazy [18] that the following properties hold.  
2
Ax = ∑∞nꢟꢀ A2〈x, en〉 en, x ∈ D(A) , where en(t) = √π siꣅ(ꣅt) , ꣅ = 1,ꢛ,꣉, … is the orthogonal set of eigenvectors.  
∗
A is the infinitesimal generator of an analytic semigroup ꢃS(t)ꢄꢈꢉꢊin X.  
ꢎꢈ  
S(t)x = ∑∞nꢟꢀ eꢔn 〈x, en〉 en, For all x ∈ X and for every t > 0 .  
ꢎ
ꢔ
π
ꢈ
Furthermore‖S(t)‖ ≤ e  
, t ≥ 0.  
The bounded linear operator(−A)3꣊4 is well defined and is given by  
(
−A)3꣊4x = ∑∞nꢟꢀ(ꣅ)3꣊2〈x, en〉 en  
with domain Dꢃ(−A)3꣊4ꢄ = ꢙx ∈ X: ∑∞nꢟꢀ eꢔnꢎꢈ〈x, en〉 en ∈ Xꢚ . Furtheremore ꢣ(−A)3꣊4ꢣ = 1 and ꢣ(−A)3꣊4ꢣ ≤  
ꢀ
∞
ꢀ
t ꢔ(ꢀ꣊4)‖S(t)‖dt < πꢼ⁄ꢎ.Thus (Hꢀ) holds with M = 1, λ = π2 ,(H2) holds with Lꢀ = α22 ,(H3) and (H4) holds  
∫
Γ(3꣊4)  
ꢊ
2
ꢀ
with β = ꣉꣊ꢺ,L2 = ꢣ(−A)3꣊4ꢣ αꢀ2 = αꢀ2 and (H6) holds with γ = .Consequently we conclude by theorem 4.2, that  
2
the stochastic partial equation (5.1) has a unique mild solution and that this solution converges to zero in mean square if  
the parameters αꢀ, α2, α3 satisfy the following conditions.  
αꢑꢎ  
ꢫꢑꢎ  
4
3
αꢑꢎ  
αꢎꢎ  
αꢼꢎ  
ꢀ
Γ2 ꢁ4ꢂ πꢼ + π꣋ + π꣋ < .  
πꢼ +  
4
6
. Conclusion  
In this paper we derived existence, conditions ensuring the exponential decay to zero in mean square of a neutral  
stochastic delay differential equations driven by fBm. In addition we also established the case of infinite delays which  
has not yet been discussed in the context of neutral stochastic delay differential equations driven by fBm.  
References  
[
[
[
1]  
2]  
3]  
L. Coutin, Z. Qian, Stochastic differential equations for fractional Brownian motions. ComptesRendus de I'Académie des Sciences-Series I-  
Mathematics 331 (2000) 75-80. http://dx.doi.org/10.1016/S0764-4442(00)01594-9.  
Y. Hu, D. Nualart, Differential equations driven by H¨older continuous functions of order greater than 1/2, Stochastic analysis and  
applications 2 (2007) 399-413. http://dx.doi.org/10.1007/978-3-540-70847-6_17.  
D. Nualart, Y. Ouknine, Regularization of differential equations by fractional noise, Stochastic Process and their Applications 102 (2002) 103-  
1
16. http://dx.doi.org/10.1016/S0304-4149(02)00155-2.  
[
[
4]  
5]  
D. Nualart, A. R˘a¸scanu, Differential equations driven by fractional Brownian motion, CollectaneaMathematica 53 (2002) 55-81.  
D. Nualart, B. Saussereau, Malliavin calculus for stochastic differential equations driven by a fractional Brownian motion, Stochastic  
Processes and their Applications119 (2009) 391-409. http://dx.doi.org/10.1016/j.spa.2008.02.016.  
[
[
[
6]  
7]  
8]  
T. Caraballo, J. Real, T. Taniguchi, The exponential stability of neutral stochastic delay partial differential equations, Discrete and Continuous  
Dynamical Systems 18 (2007) 295–313. http://dx.doi.org/10.3934/dcds.2007.18.295.  
W. Dai, C. C. Heyde, Ito's formula with respect to fractional Brownian motion and its application, Journal of Applied Mathematics and  
Stochastic Analysis 9 (1996) 439–448. http://dx.doi.org/10.1155/S104895339600038X.  
D. Feyel, A. de la Pradelle, on fractional Brownian processes, Potential Analysis 10 (1999) 273–288.  
http://dx.doi.org/10.1023/A:1008630211913.  
[
[
[
[
[
9]  
E. Hernandez, D. N. Keck, M. A. McKibben, On a class of measure-dependent stochastic evolution equations driven by fBm, Journal of  
Applied Mathematics and Stochastic Analysis 2007 (2007) 1-26. http://dx.doi.org/10.1155/2007/69747.  
10] P. H. Bezandry, Existence of almost periodic solutions for semilinear stochastic evolution equations driven by fractional Brownian motion,  
Electronic Journal of Differential Equations 2012 (2012) 1–21.  
11] T. E. Duncan, B. Maslowski, B. Pasik-Duncan, Fractional Brownian motion and stochastic equations in Hilbert spaces, Stochastics and  
Dynamics 2 (2002) 225–250. http://dx.doi.org/10.1142/S0219493702000340.  
12] B. Maslowski, D. Nualart, Evolution equations driven by a fractional Brownian motion, Journal of Functional Analysis 2 (2003) 277–305.  
http://dx.doi.org/10.1016/S0022-1236(02)00065-4.  
13] T. Caraballo, M. J. Garrido-Atienza, T. Taniguchi, The existence and exponential behavior of solutions to stochastic delay evolution equations  
with  
a
fractional Brownian motion, Nonlinear Analysis: Theory, Methods and Applications 74 (2011) 3671–3684.  
http://dx.doi.org/10.1016/j.na.2011.02.047.  
[
[
[
14] B. Boufoussi, S. Hajji, Neutral stochastic functional differential equations driven by a fractional Brownian motion in a Hilbert space, Statistics  
&
Probability Letters 82 (2012)1549–1558. http://dx.doi.org/10.1016/j.spl.2012.04.013.  
15] T. E. Duncan, B. Maslowski, B. Pasik-Duncan, Stochastic equations in Hilbert space with a multiplicative fractional Gaussian noise,  
Stochastic process and their applications 115(8) (2005) 1357-1383.  
16] Nguyen-Tien-Dung, Neutral stochastic differential equations driven by a fractional Brownian motion with impulsive effects and varying time  
delays, Journal of the Korean Statistical society 43 (2014) 599–608. http://dx.doi.org/10.1016/j.jkss.2014.02.003.  
17] T. A. Burton. Stability by fixed point theory for functional differential equations, Mineola, NY, Dover Publications, Inc. (2006).  
18] A. Pazy, Semigroups of Linear Operators and Applications to Partial Differential Equations, Applied Mathematical Sciences, vol. 44,  
Springer-Verlag, New York. (1983). http://dx.doi.org/10.1007/978-1-4612-5561-1.  
[
[
[
19] G. Da Prato, J. Zabczyk, Stochastic equations in infinite dimensions, Cambridge University Press, Cambridge, (1992).  
http://dx.doi.org/10.1017/CBO9780511666223.