Loading [MathJax]/jax/output/SVG/jax.js
Research article Special Issues

Studies on invariant measures of fractional stochastic delay Ginzburg-Landau equations on Rn

  • Received: 15 December 2023 Revised: 26 February 2024 Accepted: 29 February 2024 Published: 15 March 2024
  • This paper is concerned with invariant measures of fractional stochastic delay Ginzburg-Landau equations on the entire space Rn. We first derive the uniform estimates and the mean-square uniform smallness of the tails of solutions in corresponding space. Then we deduce the weak compactness of a set of probability distributions of the solutions applying the Ascoli-Arzelˊa. We finally prove the existence of invariant measures by applying Krylov-Bogolyubov's method.

    Citation: Hong Lu, Linlin Wang, Mingji Zhang. Studies on invariant measures of fractional stochastic delay Ginzburg-Landau equations on Rn[J]. Mathematical Biosciences and Engineering, 2024, 21(4): 5456-5498. doi: 10.3934/mbe.2024241

    Related Papers:

    [1] Conrad Bertrand Tabi, Alidou Mohamadou, Timoleon Crepin Kofane . Soliton-like excitation in a nonlinear model of DNA dynamics with viscosity. Mathematical Biosciences and Engineering, 2008, 5(1): 205-216. doi: 10.3934/mbe.2008.5.205
    [2] Ting Kang, Yanyan Du, Ming Ye, Qimin Zhang . Approximation of invariant measure for a stochastic population model with Markov chain and diffusion in a polluted environment. Mathematical Biosciences and Engineering, 2020, 17(6): 6702-6719. doi: 10.3934/mbe.2020349
    [3] Xuan Zhang, Huiqin Jin, Zhuoqin Yang, Jinzhi Lei . Effects of elongation delay in transcription dynamics. Mathematical Biosciences and Engineering, 2014, 11(6): 1431-1448. doi: 10.3934/mbe.2014.11.1431
    [4] Heping Ma, Hui Jian, Yu Shi . A sufficient maximum principle for backward stochastic systems with mixed delays. Mathematical Biosciences and Engineering, 2023, 20(12): 21211-21228. doi: 10.3934/mbe.2023938
    [5] Edward J. Allen . Derivation and computation of discrete-delayand continuous-delay SDEs in mathematical biology. Mathematical Biosciences and Engineering, 2014, 11(3): 403-425. doi: 10.3934/mbe.2014.11.403
    [6] Tuan Anh Phan, Jianjun Paul Tian . Basic stochastic model for tumor virotherapy. Mathematical Biosciences and Engineering, 2020, 17(4): 4271-4294. doi: 10.3934/mbe.2020236
    [7] Tingting Yu, Sanling Yuan . Dynamics of a stochastic turbidostat model with sampled and delayed measurements. Mathematical Biosciences and Engineering, 2023, 20(4): 6215-6236. doi: 10.3934/mbe.2023268
    [8] H.T. Banks, S. Dediu, H.K. Nguyen . Sensitivity of dynamical systems to parameters in a convex subset of a topological vector space. Mathematical Biosciences and Engineering, 2007, 4(3): 403-430. doi: 10.3934/mbe.2007.4.403
    [9] Xinran Zhou, Long Zhang, Tao Zheng, Hong-li Li, Zhidong Teng . Global stability for a class of HIV virus-to-cell dynamical model with Beddington-DeAngelis functional response and distributed time delay. Mathematical Biosciences and Engineering, 2020, 17(5): 4527-4543. doi: 10.3934/mbe.2020250
    [10] Chunhua Shan, Hongjun Gao, Huaiping Zhu . Dynamics of a delay Schistosomiasis model in snail infections. Mathematical Biosciences and Engineering, 2011, 8(4): 1099-1115. doi: 10.3934/mbe.2011.8.1099
  • This paper is concerned with invariant measures of fractional stochastic delay Ginzburg-Landau equations on the entire space Rn. We first derive the uniform estimates and the mean-square uniform smallness of the tails of solutions in corresponding space. Then we deduce the weak compactness of a set of probability distributions of the solutions applying the Ascoli-Arzelˊa. We finally prove the existence of invariant measures by applying Krylov-Bogolyubov's method.



    The nonlinear Ginzburg-Landau equation plays an important role in the studies of physics, which describes many interesting phenomena and has been studied extensively (see [1] for a more detailed description). The fractional Ginzburg-Landau equation [2,3,4] is employed to describe processes in media with fractional dispersion or long-range interaction. It becomes very popular because the fractional derivative and fractional integral have broad applications in different fields of science [5,6,7,8,9,10].

    Our work focuses on the existence of invariant measures of the autonomous fractional stochastic delay Ginzburg-Landau equations on Rn:

    du(t)+(1+iν)(Δ)αu(t)dt+(1+iμ)|u(t)|2βu(t)dt+λu(t)dt=G(x,u(tρ))dt,+k=1(σ1,k(x)+κ(x)σ2,k(u(t)))dWk(t), t>0, (1.1)

    with initial condition

    u(s)=φ(s),s[ρ,0], (1.2)

    where u(x,t) is a complex-valued function on Rn×[0,+). In (1.1), i is the imaginary unit, α,β,μ,ν and λ are real constants with β>0,λ>0 and ρ>0. (Δ)α with 0<α<1 is the fractional Laplace operator, σ1,k(x)L2(Rn) and σ2,k(u):CR are nonlinear functions, κ(x)L2(Rn)L(Rn) and {Wk}k=1 is a sequence of independent standard real-valued Wiener process on a complete filtered probability space (Ω, F, {Ft}tR,P), where {Ft}tR is an increasing right continuous family of sub-σ-algebras of F that contains all P-null sets.

    The Ginzburg-Landau equation with fractional derivative was first introduced in [2]. There is a large amount of literature which was used for investigating fractional deterministic Ginzburg-Landau equations such as [1] and stochastic equations such as [11,12,13,14,15,16,17]. These papers had respectively researched the long-time deterministic as well as random dynamical systems of fractional equations with autonomous forms and non-autonomous forms. However, in spite of quite a lot of contribution of the works, no result is provided for the existence of pathwise pullback random attractors and invariant measures for the delay stochastic Ginzburg-Landau equations.

    The delay differential equations [18] was described the dynamical systems that rely on current and past historical states. For the past few years, researchers had made great progress in the study of linear and nonlinear delay differential equations, see [20,21]. Delay differential equations are widely used in many fields, so investigating the solutions of equations has profound significance. Therefore, it's necessary that we establish the dynamics of delay stochastic Ginzburg-Landau equations.

    The goal of this paper is to prove the existence of invariant measures of the stochastic Eqs (1.1) and (1.2) in L2(Ω;C([ρ,0],L2(Rn))) by applying Krylov-Bogolyubov's method. The main difficulty of this paper is that deducing the uniform estimates of solutions (because of the nonlinear term (1+iμ)|u(t)|2βu(t) and complex-valued solutions), proving the weak compactness of a set distribution laws of the segments of solutions in L2(Ω;C([ρ,0],L2(Rn))) (because the standard Sobolev embeddings are not compact on unbounded domains Rn), and establishing the equicontinuity of solutions in L2(Ω;C([ρ,0],L2(Rn))) (because the uniform estimates in L2(Ω;C([ρ,0],L2(Rn))) are not sufficient, and the uniform estimates in L2(Ω;C([ρ,0],H1(Rn))) are needed).

    For the estimates of the nonlinear term (1+iμ)|u(t)|2βu(t), we apply integrating by parts and nonnegative definite quadratic form. There are Several methods to handle the noncompact on unbounded domain, including weighted spaces [22,23,24], weak Feller approach [25,26] and uniform tail-estimates [23,27]. We first obtain the uniform estimates of the tail of the solution as well as the technique of dyadic division, then establish the weak compactness of a set of probability distribution of solutions in C([ρ,0],L2(Rn)) applying the Ascoli-Arzelˊa theorem.

    Let S be the Schwartz space of rapidly decaying C functions on Rn. The fractional Laplace operator (Δ)α for 0<α<1 is defined by, for uS,

    (Δ)αu(x)=12C(n,α)Rnu(x+y)+u(xy)2u(x)|y|n+2αdy,    xRn,

    where C(n,α) is a positive constant given by

    C(n,α)=α4αΓ(n+2α2)πn2Γ(1α).

    By [28], the inner product ((Δ)α2u,(Δ)α2v) in the complex field is defined by

    ((Δ)α2u,(Δ)α2v)=C(n,α)2RnRn(u(x)u(y))(ˉv(x)ˉv(y))|xy|n+2αdxdy,

    for uHα(Rn). The fractional Sobolev space Hα(Rn) is endowed with the norm

    u2Hα(Rn)=u2L2(Rn)+2C(n,α)(Δ)α2u2L2(Rn).

    About the fractional derivative of fractional Ginzburg-Landau equations, there is another statement in [29].

    We organize the article as follows. In Section 2, we establish the well-posedness of (1.1) and (1.2) in L2(Ω;C([ρ,0],H)). In Sections 3 and 4, we derive the uniform estimates of solutions in L2(Ω;C([ρ,0],H)) and L2(Ω;C([ρ,0],V)), respectively. In Section 5, the existence of invariant measures is obtained.

    In this section, we show the nonlinear drift term and the diffusion term in (1.1) which are needed for the well-posedness of the stochastic delay Ginzburg-Landau Eqs (1.1) and (1.2) defined on Rn.

    We assume that G:Rn×CC is continuous and satisfies

    |G(x,u)||h(x)|+a|u|, xRn, uC (2.1)

    and

    |G(x,u)||ˆh(x)|+ˆa|u|, xRn, uC, (2.2)

    where a and ˆa>0 are constants and h(x),ˆh(x)L2(Rn). Moreover, G(x,u) is Lipschitz continuous in uC uniformly with respect to xRn. More precisely, there exists a constant CG>0 such that

    |G(x,u1)G(x,u2)|CG|u1u2|, xRn, u1,u2C. (2.3)

    For the diffusion coefficients of noise, we suppose that for each kN+

    k=1σ1,k2<, (2.4)

    and that σ2,k(u):CR is globally Lipschitz continuous; namely, for every kN+, there exists a positive number αk such that for all s1,s2C,

    |σ2,k(s1)σ2,k(s2)|αk|s1s2|. (2.5)

    We further assume that for each kN+, there exist positive numbers βk, ˆβk, γk and ˆγk such that

    |σ2,k(s)|βk+γk|s|, sC, (2.6)

    and

    |σ2,k(s)|ˆβk+ˆγk|s|, sC, (2.7)

    where k=1(α2k+β2k+γ2k+ˆβ2k+ˆγ2k)<+. In this paper, we deal with the stochastic Eqs (1.1) and (1.2) in the space C([ρ,0],L2(Rn)). In the following discussion, we denote by H=L2(Rn), V=H1(Rn).

    A solution of problems (1.1) and (1.2) will be understood in the following sense.

    Definition 2.1. We suppose that φ(s)L2(Ω,C([ρ,0],H)) is F0-measurable. Then, a continuous H-valued Ft-adapted stochastic process u(x,t) is named a solution of problems (1.1) and (1.2), if

    1) u is pathwise continuous on [0,+), and Ft-adapted for all t0,

    uL2(Ω,C([0,T],H))L2(Ω,L2([0,T],V))

    for all T>0,

    2) u(s)=φ(s) for ρs0,

    3) For all t0 and ξV,

    (u(t),ξ)+(1+iν)t0((Δ)α2u(s),(Δ)α2ξ)ds+t0Rn(1+iμ)|u(s)|2βu(s)ξ(x)dxds+λt0(u(s),ξ)ds=(φ(0),ξ)+t0(G(s,u(sρ)),ξ)ds+k=1t0(σ1,k(x)+κ(x)σ2,k(u(s)),ξ)dWk(s), (2.8)

    for almost all ωΩ.

    By the Galerkin method and the argument of Theorem 3.1 in [30], one can verify that if (2.1)–(2.7) hold true, then, for every F0-measurable function φ(s)L2(Ω,C([ρ,0],H)), the problems (1.1) and (1.2) has a unique solution u(x,t) in the sense of Definition 2.1.

    Now, we establish the Lipschitz continuity of the solutions of the problems (1.1) and (1.2) with respect to the initial data in L2(Ω,C([ρ,0],H)).

    Theorem 2.2. Suppose (2.1)–(2.6) hold, and F0-measurable function φ1,φ2L2(Ω,C([ρ,0],H)). If u1=u(t,φ1) and u2=u(t,φ2) are the solutions of the problems (1.1) and (1.2) with initial data φ1 and φ2, respectively, then, for any t0,

    E[supρstu(s,φ1)u(s,φ2)2]+E[t0u(s,φ1)u(s,φ2)2Vds]
    C1e˜C1tE[supρs0φ1(s)φ2(s)]2],

    where C1 and ˜C1 are positive constants independent of φ1 and φ2.

    Proof. Since both u1 and u2 are the solutions of the problems (1.1) and (1.2), we have, for all t0,

    u1u2+(1+iν)t0(Δ)α(u1u2)ds+(1+iμ)t0(|u1|2βu1|u2|2βu2)ds+λt0(u1u2)ds=φ1(0)φ2(0)+t0(G(x,u1(sρ))G(x,u2(sρ)))ds+k=1t0κ(x)(σ2,k(u1)σ2,k(u2))dWk. (2.9)

    By (2.9), the integration by parts of Ito's formula and taking the real parts, we get, for all t0,

    u1u22+2t0(Δ)α2(u1u2)2ds+2Ret0Rn(ˉu1ˉu2)[|u1|2βu1|u2|2βu2]dxds+2λt0u1u22ds=φ1(0)φ2(0)2+2Ret0(u1u2,G(x,u1(sρ))G(x,u2(sρ)))ds+k=1t0κ(x)(σ2,k(u1)σ2,k(u2))2ds+2Ret0(u1u2,k=1κ(x)(σ2,k(u1)σ2,k(u2)))dWk(s). (2.10)

    For the third term in the first row of (2.10), one has

        2Ret0Rn(ˉu1ˉu2)[|u1|2βu1|u2|2βu2]dxds=t0Rn2|u1|2β+2+2|u2|2β+22Re(u1ˉu2)(|u1|2β+|u2|2β)dxdst0Rn2|u1|2β+2+2|u2|2β+22|u1||u2|(|u1|2β+|u2|2β)dxdst0Rn2|u1|2β+2+2|u2|2β+2(|u1|2+|u2|2)(|u1|2β+|u2|2β)dxds=t0Rn|u1|2β+2+|u2|2β+2|u1|2β|u2|2|u2|2β|u1|2dxds=t0Rn(|u1|2β|u2|2β)(|u1|2|u2|2)dxds0.

    By (2.10), we deduce that for t0,

        E[sup0rtu1(r)u2(r)2]E[supρs0φ1(s)φ2(s)2]+2E[t0u1u2G(x,u1(sρ))G(x,u2(sρ))ds]   +k=1E[t0κ(σ2,k(u1)σ2,k(u2))2ds]   +2E[sup0rt|k=1r0(u1u2,κ(x)(σ2,k(u1)σ2,k(u2))dWk(s))|]. (2.11)

    For the second term on the right-hand side of (2.11), by (2.3), one has

        2E[t0u1u2G(x,u1(sρ))G(x,u2(sρ))ds]E[t0u1u22ds]+E[t0G(x,u1(sρ))G(x,u2(sρ))2ds]E[t0u1u22ds]+C2GE[t0u1(sρ)u2(sρ)2ds]=E[t0u1u22ds]+C2GE[tρρu1u22ds](1+C2G)E[t0u1u22ds]+C2GE[0ρφ1(s)φ2(s)2ds](1+C2G)t0E[sup0rsu1u22]ds+ρC2GE[supρs0φ1(s)φ2(s)2].

    For the third term on the right-hand side of (2.11), by (2.5), we have

      k=1E[t0κ(x)(σ2,k(u1)σ2,k(u2))2ds]κ(x)2Lk=1α2kE[t0u1u22ds]κ(x)2Lk=1α2kt0E[sup0rsu1u22]ds. (2.12)

    For the forth term on the right-hand side of (2.11), by Burkholder-Davis-Gundy's inequality, one has

        2E[sup0rt|k=1r0(u1u2,κ(x)(σ2,k(u1)σ2,k(u2))dWk(s))|]B1E[(t0k=1|(u1u2,κ(x)(σ2,k(u1)σ2,k(u2)))|2ds)12]B1E[(t0k=1u1u22κ2Lσ2,k(u1)σ2,k(u2)2ds)12]B1E[sup0stu1u2κL(k=1α2k)12(t0u1u22ds)12]12E[sup0stu1u22]+12B21κ2Lk=1α2kE[t0sup0rsu1u22ds], (2.13)

    where B1 is a constant produced by Burkholder-Davis-Gundy's inequality.

    It follows from (2.11)–(2.13) that for all t0,

    E[sup0rtu1(r)u2(r)2]2(1+ρC2G)E[supρs0φ1(s)φ2(s)2]+2[1+C2G+(1+12B21)κ2Lk=1α2k]t0E[sup0rsu1(r)u2(r)2]ds. (2.14)

    Applying Gronwall inequality to (2.14), we obtain that for all t0,

    E[sup0rtu1(r)u2(r)2]2(1+ρC2G)ec1tE[supρs0φ1(s)φ2(s)2], (2.15)

    where c1=2[1+C2G+(1+12B21)κ2Lk=1α2k]. By (2.10), there exists c2 such that for all t0,

    E[t0u1u22Vds]˜c2ec2tE[supρs0φ1(s)φ2(s)2].

    We assume that a, αk and γk are small enough in the sense, there exists a constant p2 such that

    2112p(2p1)2p12pa+2p(2p1)κ2Lk=1(α2k+γ2k)<pλ. (3.1)

    By (3.1), one has

    2κ2Lk=1γ2k<λ, (3.2)

    and

    2a+2κ2Lk=1γ2k<λ. (3.3)

    The inequalities (3.1)–(3.3) are used to establish the uniform tail-estimate of the solution of (1.1) and (1.2).

    Lemma 3.1. Suppose (2.1)–(2.6) and (3.2) hold. If φ(s)L2(Ω;C([ρ,0],H)), then, for all t0, there exists a positive constant μ1 such that the solution u of (1.1) and (1.2) satisfies

    E[u(t)2]+t0eμ1(st)E(u(s)2V)ds+t0eμ1(st)E(u(s)2β+2L2β+2)ds
    M1E[supρs0φ(s)2]+~M1, (3.4)

    and

    t+ρ0E[u(s)2V]ds(M1(t+ρ)+1+2aρC(n,α))E[supρs0φ(s)2]+2(t+ρ)aC(n,α)h(x)2
    +2(t+ρ)C(n,α)k=1(σ1,k2+2β2kκ(x)2)+˜M1(t+ρ),

    where ˜M1 is a positive constant independent of φ.

    Proof. By (1.1) and the integration by parts of Ito's formula, we have for all t0,

        u(t)2+2t0(Δ)α2u(s)2ds+2t0u(s)2β+2L2β+2ds+2λt0u(s)2ds=2Ret0(u(s),G(x,u(sρ)))ds+φ(0)2+k=1t0σ1,k(x)+κ(x)σ2,k(u(s))2ds+2Ret0(u(s),k=1σ1,k(x)+κ(x)σ2,k(u(s)))dWk(s). (3.5)

    The system (3.5) can be rewritten as

        d(u(t)2)+2(Δ)α2u(t)2dt+2u(t)2β+2L2β+2dt+2λu(t)2dt=2Re(u(t),G(x,u(tρ)))dt+k=1σ1,k(x)+κ(x)σ2,k(u(t))2dt+2Re(u(t),k=1σ1,k(x)+κ(x)σ2,k(u(t)))dWk(t). (3.6)

    Assume that μ1 is a positive constant, one has

        eμ1tu(t)2+2t0eμ1s(Δ)α2u(s)2ds+2t0eμ1su(s)2β+2L2β+2ds=(μ12λ)t0eμ1su(s)2ds+φ(0)2+2Ret0eμ1s(u(s),G(x,u(sρ)))ds+k=1t0eμ1sσ1,k+κσ2,k(u(s))2ds+2Ret0eμ1s(u(s),k=1σ1,k(x)+κ(x)σ2,k(u(s)))dWk(s).

    Taking the expectation, we have for all t0,

        eμ1tE(u(t)2)+2E[t0eμ1s(Δ)α2u(s)2ds]+2E[t0eμ1su(s)pLpds]=E(φ(0)2)+(μ12λ)E[t0eμ1su(s)2ds]+2E[t0eμ1sRe(u(s),G(x,u(sρ)))ds]+k=1E[t0eμ1sσ1,k(x)+κ(x)σ2,k(u(s))2ds]. (3.7)

    For the third term on the right-hand side (3.7), by (2.1), we have

        2E[t0eμ1sRe(u(s),G(x,u(sρ)))ds]2t0eμ1sE[u(s)G(x,u(sρ))]ds2at0eμ1sE(u(s)2)ds+22at0eμ1sE[G(x,u(sρ))2]ds2at0eμ1sE(u(s)2)ds+2at0eμ1sh(x)2ds+2at0eμ1sE[u(sρ)2]ds2a(1+eμ1ρ)t0eμ1sE[u(s)2]ds+2ah(x)2t0eμ1sds+2aeμ1ρ0ρeμ1sE[φ(s)2]ds2a(1+eμ1ρ)t0eμ1sE[u(s)2]ds+2eμ1taμ1h(x)2+2aρeμ1ρE[supρs0φ(s)2]. (3.8)

    For the forth term on the right-hand side (3.7), by (2.6), we have

        k=1E[t0eμ1sσ1,k+κσ2,k(u(s))2ds]k=1E[t0eμ1s(2σ1,k2+2κσ2,k(u(s))2)ds]2μ1k=1σ1,k2eμ1t+4k=1t0eμ1sE[β2kκ2+γ2kκ2Lu(s)2]ds2μ1k=1(σ1,k2+2β2kκ(x)2)eμ1t+4k=1γ2kκ(x)2Lt0eμ1sE(u(s)2)ds. (3.9)

    By (3.7)–(3.9), we obtain for all t0,

        eμ1tE(u(t)2)+2E[t0eμ1s(Δ)α2u(s)2ds]+2E[t0eμ1su(s)2β+2L2β+2ds](1+2aρeμ1ρ)E[supρs0φ(s)2]+[μ12λ+2a(1+eμ1ρ)+4k=1γ2kκ2L]t0eμ1sE[u2]ds+2aμ1eμ1th(x)2+2μ1k=1(σ1,k2+2β2kκ(x)2)eμ1t. (3.10)

    By (3.2), there exists a positive constant μ1 sufficiently small such that

    2μ1+2a+2aeμ1ρ+4k=1γ2kκ(x)2L2λ.

    Then, we have, for all t0,

        E(u(t)2)+2t0eμ1(st)E((Δ)α2u(s)2)ds     +μ1t0eμ1(st)E(u(s)2)ds+2t0eμ1(st)E(u(s)2β+2L2β+2)ds(1+2aρeμ1ρ)E(supρs0φ(s)2)+1μ1(2ah(x)2+2k=1(σ1,k2+2β2kκ(x)2)),

    which completes the proof of (3.4).

    Integrating (3.6) on [0,t+ρ] and taking the expectation, one has

      E[u(t+ρ)2]+2E[t+ρ0(Δ)α2u(s)2ds]+2E[t+ρ0u(s)2β+2L2β+2ds]+2λE[t+ρ0u(s)2ds]=E[φ(0)2]+2E[t+ρ0Re(u(s),G(x,u(sρ)))ds]+k=1E[t+ρ0σ1,k+κ(x)σ2,k(u(s))ds]. (3.11)

    For the second term on the right-hand side of (3.11), by (2.1), we have

        2E[t+ρ0Re(u,G(x,u(sρ)))ds]22aE[t+ρ0u2ds]+2aρE[supρs0φ(s)2]+2(t+ρ)ah2. (3.12)

    For the third term on the right-hand side of (3.11), one has

        k=1E[t+ρ0σ1,k+κσ2,k(u(s))ds]2(t+ρ)k=1(σ1,k2+2β2kκ2)+4k=1γ2kκ(x)2LE[t+ρ0u2ds]. (3.13)

    Then, by (3.2) and (3.11)–(3.13), for all t0, we obtain,

    2E[t+ρ0(Δ)α2u(s)2ds](1+2aρ)E[supρs0φ(s)2]
    +2(t+ρ)k=1(σ1,k2+2β2kκ(x)2)+2(t+ρ)ah(x)2.

    The result then follows from (3.4).

    The next lemma is used to obtain the uniform estimates of the segments of solutions in C([ρ,0],H).

    Lemma 3.2. Suppose (2.1)–(2.6) and (3.2) hold. Then, for any φ(s)L2(Ω,F0;C([ρ,0],H)), the solution of (1.1) satisfies that, for all tρ,

    E(suptρrtu(r)2)M2E[supρs0φ(s)2]+˜M2,

    where M2 and ˜M2 are positive constants independent of φ.

    Proof. By (1.1) and integration by parts of Ito's formula and taking the real part, we get for all tρ and tρrt,

        u(r)2+2rtρ(Δ)α2u(s)2ds+2rtρu(s)2β+2L2β+2ds+2λrtρu(s)2ds=u(tρ)2+2Rertρ(u(s),G(x,u(sρ)))ds+k=1rtρσ1,k(x)+κ(x)σ2,k(u(s)))2ds+2Rek=1rtρ(u(s),(σ1,k(x)+κ(x)σ2,k(u(s))dWk(s)). (3.14)

    For the second term on the right-hand side of (3.14), by (2.1) we have, for all tρ and tρrt,

        2Rertρ(u(s),G(x,u(sρ)))ds2rtρu(s)G(x,u(sρ))dsrtρu(s)2ds+rtρG(x,u(sρ))2dsrtρu(s)2ds+2rtρh2ds+2a2rtρu(sρ)2dsrtρu(s)2ds+2ρh2+2a2tρt2ρu(s)2ds. (3.15)

    For the third term on the right-hand side of of (3.14), for all tρ and tρrt, by (2.6), we have

        k=1rtρσ1,k(x)+κ(x)σ2,k(u(s))2ds2ρk=1σ1,k2+4ρκ2k=1β2k+4κ2Lk=1γ2krtρu(s)2ds. (3.16)

    By (3.14)–(3.16), we obtain for all tρ and tρrt,

    u(r)2c3+u(tρ)2+c4rt2ρu(s)2ds
    +2Rek=1rtρ(u(s),(σ1,k(x)+κ(x)σ2,k(u(s))dWk(s)), (3.17)

    where c3=2ρh2+2ρk=1σ1,k2+4ρκ2k=1β2k and c4=1+2a2+4κ2Lk=1γ2k. By (3.17), we find that for all tρ,

    E[suptρrtu(r)2]c3+E[u(tρ)2]+c4tt2ρE[u(s)2]ds+2E[suptρrt|k=1rtρ(u(s),(σ1,k(x)+κ(x)σ2,k(u(s))dWk(s))|]. (3.18)

    For the second term and the third term on the right-hand side of (3.18), by Lemma 3.1, we deduce for all tρ,

    E[u(tρ)2]sups0E[u(s)2]M1E[supρs0φ2]+˜M1 (3.19)

    and

    c4tt2ρE[u(s)2]ds2ρc4supsρE[u(s)2]c5E[supρs0φ2]+c5. (3.20)

    For the last term on the right-hand side of (3.18), by Burkholder-Davis-Gundy's inequality and Lemma 3.1, we obtain for all tρ,

        2E[suptρrt|k=1rtρ(u(s),σ1,k(x)+κ(x)σ2,k(u(s))dWk(s))|]2B2E[(k=1ttρ|(u(s),σ1,k+κσ2,k(u(s)))|2ds)12]12E[suptρstu(s)2]+2B22E[k=1ttρσ1,k+κσ2,k(u(s))2ds]12E[suptρstu(s)2]+2B22(2ρk=1σ1,k2+4ρκ2k=1β2k)+8B22ρκ2Lk=1γ2ksups0E[u(s)2]. (3.21)

    By Lemma 3.1 and (3.18)–(3.21), we deduce that for all tρ,

    E[suptρrtu(r)2]M2E[supρs0φ(s)2]+˜M2.

    This completes the proof.

    To establish the tightness of a family of distributions of solutions, we now derive uniform estimates on the tails of solutions to the problems (1.1) and (1.2).

    Lemma 3.3. Suppose (2.1)–(2.6) and (3.2) hold. If φ(s)L2(Ω,C([ρ,0],H)). Then, for all t0, the solution u of (1.1) and (1.2) satisfies

    lim supmsuptρ|x|mE[|u(t,x)|2]dx=0.

    Proof. We suppose that θ(x):RnR is a smooth function with 0θ(x)1, for all xRn defined by

    θ(x)={0if |x|1,1if |x|2.

    For fixed mN, we denote that θm(x)=θ(xm). By (1.1), we have

    d(θmu)+(1+iν)θm(Δ)αudt+(1+iμ)θm|u|2βudt+λθmudt=θmG(x,u(tρ))dt
    +k=1θm(σ1,k+κσ2,k)dWk(t). (3.22)

    By (3.2), We can find μ2 sufficiently small such that

    μ2+22a+4κ2Lk=1γ2k2λ<0. (3.23)

    By (3.22) and integration by parts of Ito's formula and taking the expectation, we obtain

        E[θmu2]+2t0eμ2(st)E[Rnθ2m|u|2β+2dx]ds=eμ2tE[θmφ(0)2]2t0eμ2(st)E[Re(1+iν)((Δ)α2u,(Δ)α2(θ2mu))]ds+(μ22λ)t0eμ2(st)E[θmu2]ds+2t0eμ2(st)E[Re(θmu,θmG(x,u(sρ)))]ds+k=1t0eμ2(st)E[θm(σ1,k+κ(x)σ2,k(u(s)))2]ds. (3.24)

    For the first term in the second row of (3.24), since φ(s)L2(Ω,C([ρ,0],H)), we have for all s[ρ,0], E[φ(0)2]<. It follows that for any ε>0, there exists a positive N1=N1(ε,φ)1, for all mN1, one has |x|mE[φ2(0,x)]dx<ε. Consequently,

        E[θmφ(0)2]=E[Rn|θ(xm)φ(0,x)|2dx]=E[|x|m|θ(xm)φ(0,x)|2dx]|x|mE[|φ(0,x)|2]dx<ε,   mN1. (3.25)

    Now we consider the second term on the right-hand side of (3.24). We first have

        2E[Re(1+iν)((Δ)α2u(s),(Δ)α2(θ2mu(s)))]=C(n,α)E[Re(1+iν)RnRn[u(x)u(y)][θ2m(x)ˉu(x)θ2m(y)ˉu(y)]|xy|n+2α]dxdy=C(n,α)E[Re(1+iν)RnRn[u(x)u(y)][θ2m(x)(ˉu(x)ˉu(y))+ˉu(y)(θ2m(x)θ2m(y))]|xy|n+2α]dxdy=C(n,α)E[Re(1+iν)RnRnθ2m(x)|u(x)u(y)|2|xy|n+2αdxdy]C(n,α)E[Re(1+iν)RnRn(u(x)u(y))(θ2m(x)θ2m(y))ˉu(y)|xy|n+2αdxdy]C(n,α)E[Re(1+iν)RnRn(u(x)u(y))(θ2m(x)θ2m(y))ˉu(y)|xy|n+2αdxdy]C(n,α)1+ν2E[|RnRn(u(x)u(y))(θ2m(x)θ2m(y))ˉu(y)|xy|n+2αdxdy|]2C(n,α)1+ν2E[Rn|ˉu(y)|(Rn|(u(x)u(y))(θm(x)θm(y))||xy|n+2αdx)dy]2C(n,α)1+ν2E[u(s)(Rn(Rn|(u(x)u(y))(θm(x)θm(y))||xy|n+2αdx)2dy)12]2C(n,α)1+ν2E[u(s)(Rn(Rn|u(x)u(y)|2|xy|n+2αdxRn|(θm(x)θm(y))|2|xy|n+2αdx)dy)12]. (3.26)

    We now prove the following inequality:

    Rn|(θm(x)θm(y))|2|xy|n+2αdxc6m2α. (3.27)

    Let xy=h and hm=z, then, we obtain,

        Rn|(θm(x)θm(y))|2|xy|n+2αdx=Rn|θ(y+hm)θ(ym)|2|h|n+2αdh=Rn|θ(ym+z)θ(ym)|2mn+2α|z|n+2αmndz=1m2αRn|θ(ym+z)θ(ym)|2|z|n+2αdz=1m2α|z|1|θ(ym+z)θ(ym)|2|z|n+2αdz+1m2α|z|>1|θ(ym+z)θ(ym)|2|z|n+2αdzc6m2α|z|1|z|2|z|n+2αdz+4m2α|z|>11|z|n+2αdzc6m2α|z|11|z|n+2α2dz+4m2α|z|>11|z|n+2αdzc6ˉc6m2α+4˜c6m2α=c6ˉc6+4˜c6m2α. (3.28)

    This proves (3.27) with c6:=c6ˉc6+4˜c6. By (3.26) and (3.27), we obtain,

        2E[Re(1+iν)((Δ)α2u(s),(Δ)α2θ2mu(s))]2c6(1+ν2)C(n,α)mαE[u(s)RnRn|u(x)u(y)|2|xy|n+2αdxdy]c6(1+ν2)C(n,α)mα(E(u(s)2)+E(RnRn|u(x)u(y)|2|xy|n+2αdxdy))c6(1+ν2)C(n,α)mαE(u(s)2)+2c6(1+ν2)mαE((Δ)α2u(s)2). (3.29)

    By (3.29), for the second term on the right-hand side of (3.24), we get

        2t0eμ2sE[Re(1+iν)((Δ)α2u(s),(Δ)α2θ2mu(s))]dsc6(1+ν2)C(n,α)mαt0eμ2sE[u(s)2]ds+2c6(1+ν2)mαt0eμ2sE[(Δ)α2u(s)2]ds. (3.30)

    By Lemma 3.1, we have

        c6(1+ν2)C(n,α)mαt0eμ2(st)E[u(s)2]dsc6(1+ν2)C(n,α)mα[M1E[supρs0φ(s)2]+˜M1]t0eμ2(st)dsc6(1+ν2)C(n,α)mα1μ2[M1E[supρs0φ(s)2]+˜M1]. (3.31)

    By (3.31), we deduce that there exists N2(ε,φ)N1, for all t0 and mN2,

    c6(1+ν2)C(n,α)mαt0eμ2(st)E[u(s)2]ds<ε.

    By Lemma 3.1, there exists N3(ε,φ)N2 such that for all t0 and mN3,

        2c6(1+ν2)mαt0eμ2(st)E[(Δ)α2u2]ds2c6(1+ν2)mα[M1E[supρs0φ(s)2]+˜M1]<ε.

    For the forth term on the right-hand side of (3.24), we obtain that there exists N4(ε,φ)N3, for all t0 and mN4,

        2t0eμ2(st)E[Re(θmu,θmG(x,u(sρ)))]ds2at0eμ2(st)E[θmu(s)2]ds+12at0eμ2(st)E[θmG(x,u(sρ))2]ds2aμ2|x|mh2(x)dx+2a0ρeμ2(st)E[θmφ(s)2]ds+22at0eμ2(st)E[θmu(s)2]ds2aμ2ε+2a0ρeμ2(st)E[θmφ(s)2]ds+22at0eμ2(st)E[θmu(s)2]ds.

    Since {φ(s)L2(Ω,H)|s[ρ,0]} is compact, it has a open cover of balls with radius ε2 which denoted by {B(φi,ε2)}li=1. Since φi=φ(si)L2(Ω;C([ρ,0],H)) for i=1,2,,l, we obtain that for given ε>0,

    {φ(s)L2(Ω;C([ρ,0],H))}li=1{XL2(Ω,H)|XφiL2(Ω,H)<ε2}.

    Since φiL2(Ω,H), there exists a positive constant N5=N5(ε,φ)N4, for mN5, we have

    supi=1,2,,l|x|mE[|φ(si,x)|2]dx<ε4.

    Then,

    sups[ρ,0]|x|mE[|φ(s,x)|2]dx<ε2,mN5.

    Consequently, one has

        2t0eμ2(st)E[Re(θmu,θmG(x,u(sρ)))]ds2aμ2ε+2aρε2+22at0eμ2(st)E[θmu(s)2]ds. (3.32)

    For the fifth term on the right-hand side of (3.24), by (2.6), we obtain

        k=1t0eμ2(st)E[θm(σ1,k+κ(x)σ2,k(u(s)))2]ds2k=1t0eμ2(st)θmσ1,k2ds+2k=1t0eμ2(st)E[θmκ(x)σ2,k(u(s))2]ds2μ2k=1|x|m|σ1,k(x)|2dx+4μ2k=1β2k|x|mκ2(x)dx+4κ(x)2Lk=1γ2kt0eμ2(st)E[θmu(s)2]ds.

    Since k=1σ1,k2< and κ(x)L2(Rn)L(Rn), there exists N6=N6(ε,φ)N5, for all t0 and mN6, we have

    k=1|x|m|σ1,k(x)|2dx+|x|mκ2(x)dx<ε.

    Consequently, for the fifth term on the right-hand side of (3.24), we get for all t0 and mN6,

    k=1t0eμ2(st)E[θm(σ1,k+κσ2,k)2]ds2μ2(1+2k=1β2k)ε
    +4κ2Lk=1γ2kt0eμ2(st)E[θmu(s)2]ds.

    Therefore, for all t0 and mN6,

    E[θmu(t)2][2+eμ2t+2aμ2+22aρ+2μ2(1+2k=1β2k)]ε
    +(μ22λ+22a+4κ2Lk=1γ2k)t0eμ2(st)E[θmu(s)2]ds.

    Taking the limit in the above equation and by (3.23), we have

    lim supmsuptρ|x|mE[|u(t,x)|2]dx=0,

    which completes the proof.

    Lemma 3.4. Suppose (2.1)–(2.6) and (3.2) hold. If φ(s)L2(Ω,C([ρ,0],H)), then the solution u of (1.1) and (1.2) satisfies

    lim supmsupt0E[supr[tρ,t]|x|m|u(r,x)|2dx]=0.

    Proof. By (3.22) and integration by parts of Ito's formula and taking the real part, for all tρ and r[tρ,t], we have

        eμ2rθmu(r)2+2rtρeμ2sRnθ2m|u|2β+2dxds=eμ2(tρ)θmu(tρ)22rtρeμ2sRe(1+iν)((Δ)α2u(s),(Δ)α2θ2mu(s))ds+(μ22λ)rtρeμ2sθmu(s)2ds+2Rertρeμ2s(θmu(s),θmG(x,u(sρ)))ds+k=1rtρeμ2sθm(σ1,k+κ(x)σ2,k(u(s)))2ds+2Rek=1rtρeμ2s(θmu(s),θm(σ1,k+κσ2,k(u(s))))dWk(s). (3.33)

    By (3.33), we deduce,

        E[suptρrtθmu(r)2]E[θmu(tρ)2]2E[suptρrtrtρeμ2(sr)Re(1+iν)((Δ)α2u,(Δ)α2θ2mu)ds]+|μ22λ|E[suptρrtrtρθmu2eμ2(sr)ds]+2E[suptρrtrtρeμ2(sr)θmuθmG(x,u(sρ))ds]+k=1E[suptρrtrtρeμ2(sr)θm(σ1,k+κσ2,k(u(s)))2ds]+2E[suptρrt|k=1rtρeμ2(sr)(θmu(s),θm(σ1,k+κσ2,k(u(s))))dWk(s)|]. (3.34)

    For the first term on the right-hand side of (3.34), by Lemma 3.3, one has for any ε>0, there exists ˜N1(ε,φ)1 such that for all m˜N1 and tρ,

    E[θmu(tρ)2]|x|mE[|u(tρ,x)|2]dx<ε. (3.35)

    For the second term on the right-hand side of (3.34), by (3.29), we have

        2E[suptρrtrtρeμ2(sr)Re(1+iν)((Δ)α2u(s),(Δ)α2θ2mu(s))ds]2c6(1+ν2)C(n,α)mαE[suptρrt(rtρeμ2(sr)u(s)(Δ)α2u(s)ds)]2c6(1+ν2)C(n,α)mαeμ2ρE[(ttρeμ2(st)u(s)(Δ)α2u(s)ds)]c6(1+ν2)C(n,α)mαeμ2ρ{ttρeμ2(st)E[u2]ds+E[ttρeμ2(st)(Δ)α2u2ds]}c6(1+ν2)C(n,α)mαeμ2ρ{ρsups[tρ,t]E[u(s)2]+E[ttρeμ2(st)(Δ)α2u2ds]}. (3.36)

    By Lemma 3.1 and (3.36), we deduce that there exists ˜N2(ε,φ)˜N1 such that for all m˜N2 and tρ,

    2E[suptρrtrtρeμ2(sr)Re(1+iν)((Δ)α2u(s),(Δ)α2θ2mu(s))ds]<ε. (3.37)

    For the third term on the right-hand side of (3.34), by Lemma 3.3, we obtain that for all m˜N2 and tρ,

    |μ22λ|E[suptρrtrtρθmu(s)2eμ2(sr)ds]|μ22λ|E[ttρθmu(s)2ds]
    |μ22λ|ρsuptρstE[θmu(s)2]<|μ22λ|ρε. (3.38)

    For the forth term on the right-hand side of (3.34), by (2.1), we obtain

    2E[suptρrtrtρeμ2(sr)θmu(s)θmG(x,u(sρ))ds]ttρE[θmu(s)2]ds+2ρθmh2+2a2tρt2ρE[θmu(s)2]dsρsuptρstE[θmu(s)2]+2ρθmh2+2a2ρsupt2ρstρE[θmu(s)2],

    which along with Lemma 3.3, we deduce that there exists ˜N3(ε,φ)˜N2 such that for all m˜N3 and tρ,

    2E[suptρrtrtρeμ2(sr)θmu(s)θmG(x,u(sρ))ds]<(3+2a2)ρε. (3.39)

    For the fifth term on the right-hand side of (3.34), by (2.6), we have

        k=1E[suptρrtrtρeμ2(sr)θm(σ1,k+κσ2,k(u(s)))2ds]2ρk=1θmσ1,k2+2ρk=1suptρstE[θmκ(x)σ2,k(u(s))2]2ρk=1|x|m|σ1,k(x)|2dx+4ρk=1β2k|x|m|κ(x)|2dx+4ρκ(x)2Lk=1γ2ksuptρstE[θmu(s)2].

    By the condition κ(x)L2(Rn)L(Rn), (2.4) and Lemma 3.3, we deduce that there exists ˜N4(ε,φ)˜N3 such that for all m˜N4 and tρ,

    k=1E[suptρrtrtρeμ2(sr)θm(σ1,k+κσ2,k(u(s)))2ds]<2ρ(1+λ+2k=1β2k)ε. (3.40)

    For the sixth term on the right-hand side of (3.34), by (2.6), (3.40) and Burkholder-Davis-Gundy's inequality, we have,

        2E[suptρrt|k=1rtρeμ2(sr)(θmu(s),θm(σ1,k+κσ2,k(u(s))))dWk(s)|]2eμ2(tρ)E[suptρrt|k=1rtρeμ2s(θmu(s),θmσ1,k+θmκ(x)σ2,k(u(s)))dWk(s)|]2˜B2eμ2(tρ)E[(ttρe2μ2sk=1|(θmu(s),θmσ1,k+θmκ(x)σ2,k(u(s)))|2ds)12]2˜B2eμ2(tρ)E[suptρstθmu(s)(ttρe2μ2sk=1θmσ1,k+θmκσ2,k(u(s))2ds)12]12E[suptρstθmu(s)2]+2˜B22E[e2μ2ρttρe2μ2(st)k=1θmσ1,k+θmκ(x)σ2,k(u(s))2ds]12E[suptρstθmu(s)2]+2˜B22e2μ2ρk=1E[suptρrtrtρeμ2(sr)θmσ1,k+θmκ(x)σ2,k(u(s))2ds]12E[suptρstθmu(s)2]+4ρ(1+λ+2k=1β2k)˜B22e2μ2ρε.

    Above all, for all m˜N4 and tρ, we obtain,

    E[suptρrtθmu(r)2][4+2|μ22λ|ρ+(6+4a2)ρ+4ρ(1+2˜B22e2μ2ρ)(1+λ+2k=1β2k)]ε.

    Therefore, we conclude

    lim supmsupt0E[suptρrt|x|m|u(r,x)|2dx]=0.

    Lemma 3.5. Suppose (2.1)–(2.6) and (3.1) hold. If φ(s)L2(Ω,C([ρ,0],H)), then there exists a positive constant μ3 such that the solution u of (1.1) and (1.2) satisfies

     suptρE[u(t)2p]+supt0E[t0eμ3(st)u(s)2p2(Δ)α2u(s)2ds]        (1+aρeμρ2p(4p2)2p12p)E[φ2pCH]+M3, (3.41)

    where M3 is a positive constant independent of φ.

    Proof. By (3.1), there exist positive constants μ and ϵ1 such that

     μ+aeμρ2p2112p(2p1)2p12p+4(p1)(2p1)ϵ2p2p21k=1(σ1,k2+κ2β2k)         +4θ(2p1)κ2Lk=1γ2k2pλ. (3.42)

    Given nN, let τn be a stopping time as defined by

    τn=inf{t0:u(t)>n},

    and as usual, we set τn=+ if {t0:u(t)>n}=. By the continuity of solutions, we have

    limnτn=+.

    Applying Ito's formula, we obtain

        d(u(t)2p)=d((u(t)2)p)=pu(t)2(p1)d(u(t)2)+2p(p1)u(t)2(p2)    ×k=1|(u(t),σ1,k+κσ2,k(u(t)))|2dt. (3.43)

    Substituting (3.6) into (3.43), we infer

    d(u(t)2p)=2pu(t)2(p1)(Δ)α2u(t)2dt2pu(t)2(p1)u(t)2β+2L2β+2dt2pλu(t)2pdt     +2pu(t)2(p1)Re(u(t),G(x,u(tρ)))dt     +pu(t)2(p1)k=1σ1,k(x)+κ(x)σ2,k(u(t))2dt     +2pu(t)2(p1)Re(u(t),k=1σ1,k(x)+κ(x)σ2,k(u(t)))dWk(t)     +2p(p1)u(t)2(p2)k=1|(u(t),σ1,k+κσ2,k(u(t)))|2dt. (3.44)

    We also get the formula

    d(eμtu(t)2p)=μeμtu(t)2pdt+eμtd(u(t)2p). (3.45)

    Substituting (3.44) into (3.45) and integrating on (0,tτn) with t0, we deduce

        eμ(tτn)u(tτn)2p+2ptτn0eμsu(s)2(p1)(Δ)α2u(s)2ds=2ptτn0eμsu(s)2(p1)u(s)2β+2L2β+2ds+φ(0)2p+(μ2pλ)tτn0eμsu(s)2pds      +2ptτn0eμsu(s)2(p1)Re(u(s),G(x,u(sρ)))ds      +pk=1tτn0eμsu(s)2(p1)σ1,k+κσ2,k(u(s))2ds      +2pk=1tτn0eμsu(t)2(p1)Re(u(s),σ1,k+κσ2,k(u(s)))dWk(s)      +2p(p1)k=1tτn0eμsu(s)2(p2)|(u(s),σ1,k+κσ2,k(u(s)))|2ds. (3.46)

    Taking the expectation, we obtain for t0,

        E[eμ(tτn)u(tτn)2p]+2pE[tτn0eμsu(s)2(p1)(Δ)α2u(s)2ds]=2pE[tτn0eμsu(s)2(p1)u(s)2β+2L2β+2ds]+E[φ(0)2p]+(μ2pλ)E[tτn0eμsu(s)2pds]      +2pE[tτn0eμsu(s)2(p1)Re(u(s),G(x,u(sρ)))ds]      +pk=1E[tτn0eμsu(s)2(p1)σ1,k+κσ2,k(u(s))2ds]      +2p(p1)k=1E[tτn0eμsu(s)2(p2)|(u(s),σ1,k+κσ2,k(u(s)))|2ds]E[φ(0)2p]+(μ2pλ)E[tτn0eμsu(s)2pds]      +2pE[tτn0eμsu(s)2(p1)Re(u(s),G(x,u(sρ)))ds]      +pk=1E[tτn0eμsu(s)2(p1)σ1,k+κσ2,k(u(s))2ds]      +2p(p1)k=1E[tτn0eμsu(s)2(p2)|(u(s),σ1,k+κσ2,k(u(s)))|2ds]. (3.47)

    Next, we estimate the terms on the right-hand side of (3.47).

    For the third term on the right-hand side of (3.47), by Young's inequality and (2.1), we infer

        2θE[tτn0eμsu(s)2(p1)Re(u(s),G(x,u(sρ)))ds]2θE[tτn0eμsu(s)2p1G(x,u(sρ))2ds]aeμρ2p2112p(2p1)2p12pE[tτn0eμsu(s)2pds]        +(2p122p1a2peμρ)2p12pE[tτn0eμsG(x,u(sρ))2ds]aeμρ2p2112p(2p1)2p12pE[tτn0eμsu(s)2pds]        +22p1(2p122p1a2peμρ)2p12pE[tτn0eμs(h2p+a2pu(sρ)2p)ds]aeμρ2p2112p(2p1)2p12pE[tτn0eμsu(s)2pds]        +1μ(4p2a2peμρ)2p12ph2peμt+aρeμρ2p(4p2)2p12pE[φ2pCH]. (3.48)

    For the forth term on the right-hand side of (3.47), we infer

        pk=1E[tτn0eμsu(s)2(p1)σ1,k+κσ2,k(u(s))2ds]2pk=1E[tτn0eμsu(s)2(p1)σ1,k2ds]      +2pk=1E[tτn0eμsu(s)2(p1)κσ2,k(u(s))2ds]. (3.49)

    For the first term on the right-hand side of (3.49), we have

        2pk=1E[tτn0eμsu(s)2(p1)σ1,k2ds]2(p1)ϵ2p2p21k=1σ1,k2E[tτn0eμsu(s)2pds]+2μϵp1k=1σ1,k2eμt. (3.50)

    For the second term on the right-hand side of (3.49), we have

        2pk=1E[tτn0eμsu(s)2(p1)κσ2,k(u(s))2ds]4pκ2k=1β2kE[tτn0eμsu(s)2(p1)ds]+4pκ2Lk=1γ2kE[tτn0eμsu(s)2pds]4(p1)ϵ2p2p21κ2k=1β2kE[tτn0eμsu(s)2pds]    +4μϵp1κ2k=1β2keμt+4pκ2Lk=1γ2kE[tτn0eμsu(s)2pds]. (3.51)

    By (3.49)–(3.51), we obtain

        pk=1E[tτn0eμsu(s)2(p1)σ1,k+κσ2,k(u(s))2ds][4(p1)ϵ2p2p21k=1(σ1,k2+κ2β2k)+4pκ2Lk=1γ2k]E[tτn0eμsu(s)2pds]    +2μϵp1k=1(σ1,k2+2κ2β2k))eμt. (3.52)

    For the fifth term on the right-hand side of (3.47), applying (3.52), we have

        2p(p1)k=1E[tτn0eμsu(s)2(p2)|(u(s),σ1,k+κσ2,k(u(s)))|2ds]2p(p1)k=1E[tτn0eμsu(s)2p2σ1,k+κσ2,k(u(s))2ds][8(p1)2ϵ2p2p21k=1(σ1,k2+κ2β2k)+8p(p1)κ2Lk=1γ2k]E[tτn0eμsu(s)2pads]    +4(p1)μϵp1k=1(σ1,k2+2κ2β2k))eμt. (3.53)

    From (3.47), (3.48), (3.52) and (3.53), we obtain that for t0,

        E[eμ(tτn)u(tτn)2p]+2pE[tτn0eμsu(s)2(p1)(Δ)α2u(s)2ds](1+aρeμρ2p(4p2)2p12p)E[φ2pCH]   +(μ2pλ+aeμρ2p2112p(2p1)2p12p+4(p1)(2p1)ϵ2p2p21    ×k=1(σ1,k2+κ2β2k)+4p(2p1)κ2Lk=1γ2k)E[tτn0eμsu(s)2pds]      +1μ(4p2a2peμρ)2p12ph2peμt+4(p1)μϵp1k=1(σ1,k2+2κ2β2k))eμt. (3.54)

    Then by (3.42) and (3.54), we obtain that for t0,

        E[eμ(tτn)u(tτn)2p]+2pE[tτn0eμsu(s)2(p1)(Δ)α2u(s)2ds](1+aρeμρ2p(4p2)2p12p)E[φ2pCH]+1μ(4p2a2paeμρ)2p12ph2peμt   +4(p1)μϵp1k=1(σ1,k2+2κ2β2k))eμt. (3.55)

    Letting n, by Fatou's Lemma, we deduce that for t0,

        E[eμtu(t)2p]+2pE[t0eμsu(s)2(p1)(Δ)α2u(s)2ds](1+aρeμρ2p(4θ2)2θ12p)E[φ2pCH]+1μ(4p2a2peμρ)2p12ph2peμt   +4(p1)μϵp1k=1(σ1,k2+2κ2β2k))eμt.

    Hence, we have for t0,

        E[u(t)2p]+2pE[t0eμ(st)u(s)2(p1)(Δ)α2u(s)2ds](1+aρeμρ2p(4p2)2p12p)E[φ2pCH]+1μ(4p2a2peμρ)2p12ph2p   +4(p1)μϵp1k=1(σ1,k2+2κ2β2k)).

    This implies the desired estimate.

    In this section, we establish the uniform estimates of solutions of problems (1.1) and (1.2) with initial data in C([ρ,0],V). To the end, we assume that for each kN, the function σ1,kV and

    k=1σ1,k2V<. (4.1)

    Furthermore, we assume that the function κV and there exists a constant C>0 such that

    |κ(x)|C. (4.2)

    In the sequel, we further assume that the constant a, ˆγk in (2.7) are sufficiently small in the sense that there exists a constant p2 such that

    ˆa2112p(2p1)2p12p+2p(2p1)κ2Lk=1(β2k+ˆβ2k+γ2k+ˆγ2k)<pλ2. (4.3)

    By (4.3), we can find

    2ˆa+2κ2Lk=1ˆγ2k<λ2. (4.4)

    Lemma 4.1. Suppose (2.1)–(2.7) and (4.4) hold. If φ(s)L2(Ω;C([ρ,0],V)), then, for all t0, there exists a positive constant μ4 such that the solution u of (1.1) and (1.2) satisfies

    supsρE[u(t)2]+sups0E[t0eμ4(st)(Δ)α+12u(s)2ds]M4(E[φ2CV]+1), (4.5)

    where M4 is a positive constant independent of φ.

    Proof. By (4.4), there exists a positive constant μ1 such that

    μ12λ+8κ2Lk=1ˆγ2k<0. (4.6)

    By (1.1) and applying Ito's formula to eμ1tu(t)2, we have for t0,

        eμ1tu(t)2+2t0eμ1s(Δ)α+12u(s)2ds+2t0eμ1sRe((1+iμ)|u(s)|2βu(s),Δu(s))ds=(μ12λ)t0eμ1su(s)2ds+φ(0)2+2Ret0eμ1s(G(x,u(sρ)),Δu(s))ds    +k=1t0eμ1s(σ1,k+κσ2,k(u(s)))2ds     +2k=1Ret0eμ1s(σ1,k(x)+κ(x)σ2,k(u(s)),Δu(s))dWk(s).

    Taking the expectation, we have for all t0,

        eμ1tE[u(t)2]+2E[t0eμ1s(Δ)α+12u(s)2ds]+2E[t0eμ1sRe((1+iμ)|u(s)|2βu(s),Δu(s))ds]=(μ12λ)E[t0eμ1su(s)2ds]+E[φ(0)2]+2E[Ret0eμ1s(G(x,u(sρ)),Δu(s))ds]    +k=1E[t0eμ1s(σ1,k+κσ2,k(u(s)))2ds]. (4.7)

    First, we estimate the third term on the left-hand side of (4.7). Applying integrating by parts, we have

        Re((1+iμ)|u|2βu,Δu)=Re(1+iμ)Rn((β+1)|u|2β|u|2+β|u|2(β1)(u¯u)2)dx=Rn|u|2(β1)((β+1)|u|2|u|2+β(1+iμ)2(u¯u)2+β(1iμ)2(¯uu)2)dx=Rn|u|2(β1)trace(YMYH), (4.8)

    where

    Y=(¯uuu¯u)H,M=(β+12β(1+iμ)2β(1iμ)2β+12),

    and YH is the conjugate transpose of the matrix Y. We observe that the condition β11+μ21 implies that the matrix M is nonpositive definite. Hence, we obtain

    2E[t0eμ1sRe((1+iμ)|u(s)|2βu(s),Δu(s))ds]0. (4.9)

    Next, we estimate the terms on the right-hand side of (4.7). For the third term on the right-hand side of (4.7), applying (2.2) and Gagliardo-Nirenberg inequality, we have

        2E[Ret0eμ1s(G(x,u(sρ)),Δu(s))ds]2E[t0eμ1su(s)G(x,u(sρ))ds]E[t0eμ1su(s)2ds]+E[t0eμ1sG(x,u(sρ))2ds]E[t0eμ1su(s)2ds]+2E[t0eμ1sˆh(x)2ds]+2ˆa2E[t0eμ1su(sρ)2ds]E[t0eμ1s(Δ)α+12u(s)2ds]+2μ1ˆh(x)2eμ1t    +cμ1sups0E[u(s)2]eμ1t+2ˆa2μ1supρs0E[φ(s)2]eμ1t, (4.10)

    where c is a positive constant from Gagliardo-Nirenberg inequality. For the forth term on the right-hand side of (4.7), applying (2.6) and (2.7), we have

        k=1E[t0eμ1s(σ1,k+κσ2,k(u(s)))2ds]2k=1E[t0eμ1s(σ1,k2+(κσ2,k(u(s)))2)ds]2μ1k=1σ1,k2eμ1t+8k=1E[t0eμ1s(β2kκ2+ˆβ2kκ2+γ2kC2u(s)2+ˆγ2kκ2Lu(s)2)ds]2μ1k=1(σ1,k2+4β2kκ2+4ˆβ2kκ2+4C2γ2ksups0E[u(s)2])eμ1t     +8k=1ˆγ2kκ(x)2LE[t0eμ1su(s)2ds]. (4.11)

    By (4.7), (4.10) and (4.11), we obtain

        E[u(t)2]+E[t0eμ1(st)(Δ)α+12u(s)2ds] E[φ(0)2]eμ1t+2μ1ˆh(x)2+(μ12λ+8κ2Lk=1ˆγ2k)E[t0eμ1su(s)2ds]    +2μ1(c2+4(C2k=1γ2k+cκ2Lk=1ˆγ2k))supsρE[u(s)2]    +2μ1k=1(σ1,k2+4(β2k+ˆβ2k)κ2V)+2ˆa2μ1supρs0E[φ(s)2]. (4.12)

    Then by (4.6) and (4.12), we obtain that for all t0,

        E[u(t)2]+E[t0eμ1(st)(Δ)α+12u(s)2ds] E[φ(0)2]eμ1t+2μ1ˆh(x)2+2μ1(c2+4(C2k=1γ2k+cκ2Lk=1ˆγ2k))supsρE[u(s)2]    +2μ1k=1(σ1,k2+4(β2k+ˆβ2k)κ2V)+2ˆa2μ1supρs0E[φ(s)2]. (4.13)

    Then by (4.13) and Lemma 3.1, we obtain the estimates (4.5).

    Lemma 4.2. Suppose (2.1)–(2.7) and (4.4) hold. If φ(s)L2(Ω;C([ρ,0],V)), then the solution u of (1.1) and (1.2) satisfies

    suptρ{E[suptρrtu(r)2]}M5(E[φ2CV]+1), (4.14)

    where M5 is a positive constant independent of φ.

    Proof. By (1.1) and Ito's formula, we get for all tρ and tρrt,

        u(r)2+2rtρ(Δ)α+12u(s)2ds    +2rtρRe((1+iμ)|u(s)|2βu(s),Δu(s))ds+2λrtρu(s)2ds=u(tρ)2+2Rertρ(G(x,u(sρ)),Δu(s))ds    +k=1rtρ(σ1,k+κσ2,k(u(s)))2ds     +2k=1Rertρ(σ1,k(x)+κ(x)σ2,k(u(s)),Δu(s))dWk(s). (4.15)

    For the third term on the left-hand side of (4.15), applying (4.8), we have

    2rtρRe((1+iμ)|u(s)|2βu(s),Δu(s))ds0. (4.16)

    For the second term on the right-hand side of (4.15), applying (2.2) and Gagliardo-Nirenberg inequality, we have

        2Rertρ(G(x,u(sρ)),Δu(s))ds2rtρu(s)G(x,u(sρ))dsrtρu(s)2ds+rtρG(x,u(sρ))2dsrtρu(s)2ds+2rtρˆh(x)2ds+2ˆa2rtρu(sρ)2dsrtρ(Δ)α+12u(s)2ds+2ρˆh(x)2+2ˆa2rρt2ρu(s)2ds+crtρu(s)2ds. (4.17)

    For the third term on the right-hand side of (4.15), applying (2.6) and (2.7), we have

        k=1rtρ(σ1,k+κσ2,k(u(s)))2ds2k=1rtρ(σ1,k2+(κσ2,k(u(s)))2)ds2ρk=1σ1,k2+8ρ(κ2k=1β2k+κ2k=1ˆβ2k)     +8C2k=1γ2krtρu(s)2ds+8κ2Lk=1ˆγ2krtρu(s)2ds. (4.18)

    By (4.4) and (4.15)–(4.18), we infer that for all tρ and tρrt,

        u(r)2c1+u(tρ)2+c2rt2ρu(s)2ds+2ˆa2rt2ρu(s)2ds     +2k=1Rertρ(σ1,k(x)+κ(x)σ2,k(u(s)),Δu(s))dWk(s), (4.19)

    where c1 and c2 are positive constants. By (4.19), we deduce that for all tρ,

        E[suptρrtu(r)2]c1+E[u(tρ)2]+c2rt2ρE[u(s)2+u(s)2]ds     +2E[suptρrt|k=1rtρ(σ1,k(x)+κ(x)σ2,k(u(s)),Δu(s))dWk(s)|]. (4.20)

    For the second term on the right-hand side of (4.20), by Lemma 4.1 we infer that for all tρ,

        E[u(tρ)2]supsρE[u(s)2]c3E[φ2CV]+c3. (4.21)

    For the third term on the right-hand side of (4.20), by Lemmas 3.1 and 4.1 we infer that for all tρ,

        c2rt2ρE[u(s)2+u(s)2]ds2ρc2supsρE[u(s)2+u(s)2]c4E[φ2CV]+c4. (4.22)

    For the last term on the right-hand side of (4.20), by BDG inequality, (4.18), Lemmas 3.1 and 4.1, we deduce that for all tρ,

        2E[suptρrt|k=1rtρ(σ1,k(x)+κ(x)σ2,k(u(s)),Δu(s))dWk(s)|]2c5E[(k=1ttρ|(σ1,k(x)+κ(x)σ2,k(u(s)),Δu(s))|2ds)12]2c5E[(k=1ttρu(s)2(σ1,k(x)+κ(x)σ2,k(u(s)))2ds)12]2c5E[suptρstu(s)(k=1ttρ(σ1,k(x)+κ(x)σ2,k(u(s)))2ds)12]12E[suptρstu(s)2]+2c25E[k=1ttρ(σ1,k(x)+κ(x)σ2,k(u(s)))2ds]12E[suptρstu(s)2]+c6+c6ttρE[u(s)2+u(s)2]ds12E[suptρstu(s)2]+c6+ρc6(sups0Eu(s)2+sups0u(s)2)12E[suptρstu(s)2]+c7E[φ2CV]+c7. (4.23)

    By (4.20)–(4.23), we obtain that for all tρ,

    E[suptρrtu(r)2]c8E[φ2CV]+c9,

    which completes the proof.

    Lemma 4.3. Suppose (2.1)–(2.7) and (3.1) hold. If φ(s)L2p(Ω,C([ρ,0],V)), then there exists a positive constant μ5 such that the solution u of (1.1) and (1.2) satisfies

     suptρE[u(t)2p]+supt0E[t0eμ5(st)u(s)2(p1)(Δ)α+12u(s)2ds]        M5(E[φ2pCV]+1), (4.24)

    where M5 is a positive constant independent of φ.

    Proof. By (3.1), there exist positive constants μ and ϵ1 such that

     μ+4(p1)ϵpp11+2pˆa2μϵp1+8C2(p1)(2p1)ϵpp11k=1γ2k+8p(2p1)κ2Lk=1ˆγ2k         +2(p1)(2p1)ϵpp11k=1(σ1,k2+4β2kκ2+4ˆβ2kκ2)2pλ. (4.25)

    By (1.1) and applying Ito's formula to eμtu(t)2p, we get for t0,

        eμtu(t)2p+2pt0eμsu(s)2(p1)(Δ)α+12u(s)2ds      +2pt0eμsu(s)2(p1)Re((1+iμ)|u(s)|2βu(s),Δu(s))ds=φ(0)2p+(μ2pλ)t0eμsu(s)2pds     +2pRet0eμsu(s)2(p1)(G(x,u(sρ)),Δu(s))ds     +pk=1t0eμsu(s)2(p1)(σ1,k+κσ2,k(u(s)))2ds     +2pk=1Ret0eμ1su(s)2(p1)(σ1,k+κσ2,k(u(s)),Δu(s))dWk(s)     +2p(p1)k=1Ret0eμ1su(s)2(p2)|(σ1,k+κσ2,k(u(s)),Δu(s))|2ds.

    Taking the expectation, we have for t0,

        eμtE[u(t)2p]+2pE[t0eμsu(s)2(p1)(Δ)α+12u(s)2ds]      +2pE[t0eμsu(s)2(p1)Re((1+iμ)|u(s)|2βu(s),Δu(s))ds]=E[φ(0)2p]+(μ2pλ)E[t0eμsu(s)2pds]     +2pE[Ret0eμsu(s)2(p1)(G(x,u(sρ)),Δu(s))ds]     +pk=1E[t0eμsu(s)2(p1)(σ1,k+κσ2,k(u(s)))2ds]     +2p(p1)k=1E[Ret0eμ1su(s)2(p2)|(σ1,k+κσ2,k(u(s)),Δu(s))|2ds]. (4.26)

    By (4.8), we get the third term on the left-hand side of (4.26) is nonnegative. Next, we estimate each term on the right-hand side of (4.26). For the third term on the right-hand side of (4.26), applying (2.2), Gagliardo-Nirenberg inequality and Young's inequality, we deduce

        2pE[Ret0eμsu(s)2(p1)(G(x,u(sρ)),Δu(s))ds]2pE[t0eμsu(s)2(p1)u(s)G(x,u(sρ))ds]pE[t0eμsu(s)2(p1)u(s)2ds]+pE[t0eμsu(s)2(p1)G(x,u(sρ))2ds]pE[t0eμsu(s)2(p1)u(s)2ds]    +2pE[t0eμsu(s)2(p1)ˆh2ds]+2pˆa2E[t0eμsu(s)2(p1)u(sρ)2ds]pE[t0eμsu(s)2(p1)((Δ)α+12u(s)2+cu(s)2)ds]+2ϵp1ˆh(x)2pt0eμsds    +4(p1)ϵpp11E[t0eμsu(s)2pds]+2pˆa2ϵp1E[t0eμsu(s)2pds]pE[t0eμsu(s)2(p1)(Δ)α+12u(s)2ds]+(4(p1)ϵpp11+2pˆa2μϵp1)E[t0eμsu(s)2pds]    +1μϵp1ˆh(x)2peμt+2pˆa2μϵp1supρs0E[φ(s)2p]eμt+csups0E[u(s)2p]eμt. (4.27)

    For the forth term on the right-hand side of (4.26), applying (2.7), we infer

        pk=1E[t0eμsu(s)2(p1)(σ1,k+κσ2,k(u(s)))2ds]2pk=1E[t0eμsu(s)2(p1)σ1,k2ds]+2pk=1E[t0eμsu(s)2(p1)(κσ2,k(u(s)))2ds]2pk=1E[t0eμsu(s)2(p1)σ1,k2ds]    +8pk=1E[t0eμsu(s)2(p1)(β2kκ2+ˆβ2kκ2+γ2kC2u(s)2+ˆγ2kκ2Lu(s)2)ds]2pk=1E[t0eμsu(s)2(p1)(σ1,k2+4β2kκ2+4ˆβ2kκ2)ds]    +8C2pk=1γ2kE[t0eμsu(s)2(p1)u(s)2ds]+8pκ2Lk=1ˆγ2kE[t0eμsu(s)2pds]. (4.28)

    Then applying Young's inequality, (4.28) can be estimated by

        pk=1E[t0eμsu(s)2(p1)(σ1,k+κσ2,k(u(s)))2ds]k=1(σ1,k2+4β2kκ2+4ˆβ2kκ2)×E[t0eμs((2p2)ϵpp11u(s)2p+2ϵp1)ds]    +2C2k=1γ2kE[t0eμs((4p4)ϵpp11u(s)2p+4ϵp1u(s)2p)ds]    +8pκ2Lk=1ˆγ2kE[t0eμsu(s)2pds]=[(2p2)ϵpp11k=1(σ1,k2+4β2kκ2+4ˆβ2kκ2)    +2C2(4p4)ϵpp11k=1γ2k+8pκ2Lk=1ˆγ2k]E[t0eμsu(s)2pds]    +k=1(σ1,k2+4β2kκ2+4ˆβ2kκ2)2ϵp1eμt+8C2μϵp1k=1γ2ksups0E[u(s)2p]eμt. (4.29)

    For the fifth term on the right-hand side of (4.26), applying integrating by parts and (4.29), we get

        2p(p1)k=1E[Ret0eμ1su(s)2(p2)|(σ1,k+κσ2,k(u(s)),Δu(s))|2ds]2p(p1)k=1E[t0eμ1su(s)2p2(σ1,k+κσ2,k(u(s)))2ds]2(p1)[(2p2)ϵpp11k=1(σ1,k2+4β2kκ2+4ˆβ2kκ2)    +2C2(4p4)ϵpp11k=1γ2k+8pκ2Lk=1ˆγ2k]E[t0eμsu(s)2pds]     +k=1(σ1,k2+4β2kκ2+4ˆβ2kκ2)4(p1)ϵp1eμt+16C2(p1)μϵp1k=1γ2ksups0E[u(s)2p]eμt. (4.30)

    By (4.26), (4.27), (4.29) and (4.30), we obtain

        E[u(t)2p]+pE[t0eμ(st)u(s)2(p1)(Δ)α+12u(s)2ds]E[φ(0)2p]eμt+[μ2pλ+4(ϱ1)ϵpp11+2pˆa2μϵp1+8C2(p1)(2p1)ϵpp11k=1γ2k    +8p(2p1)κ2Lk=1ˆγ2k+2(p1)(2p1)ϵpp11k=1(σ1,k2+4β2kκ2+4ˆβ2kκ2)]    ×E[t0eμ(st)u(s)2pds]+1μϵp1ˆh(x)2p+2pˆa2μϵp1supρs0E[φ(s)2p]    +csups0E[u(s)2p]+k=1(σ1,k2+4β2kκ2+4ˆβ2kκ2)4p2ϵp1    +8C2(2p1)μϵp1k=1γ2ksups0E[u(s)2p]. (4.31)

    Then by (4.25) and (4.31), we deduce that for all t0,

        E[u(t)2p]+pE[t0eμ(st)u(s)2(p1)(Δ)α+12u(s)2ds]E[φ(0)2p]eμt+1μϵp1ˆh(x)2p+2pˆa2μϵp1supρs0E[φ(s)2p]    +csups0E[u(s)2p]+k=1(σ1,k2+4β2kκ2+4ˆβ2kκ2)4p2ϵp1    +8C2(2p1)μϵp1k=1γ2ksups0E[u(s)2p]. (4.32)

    Therefore, by (4.32) and Lemma 3.5, there exists a constant M5 independent of φ such that

        suptρE[u(t)2p]+supt0E[t0eμ(st)u(s)2(p1)(Δ)α+12u(s)2ds]M5(E[φ2pCV]+1). (4.33)

    For convenience, we write A=(1+iν)()α+λI. Then, similar to Theorem 6.5 in [31], the solution of (1.1) and (1.2) can be expressed as

     u(t)=eAtu(0)t0eA(ts)(1+iμ)|u(s)|2βu(s)ds   +t0eA(ts)G(,u(sρ))ds+k=1t0eA(ts)(σ1,k+κσ2,k(u(s)))dWk(s). (4.34)

    The next lemma is concerned with the H¨older continuity ofsolutions in time which is needed to prove the tightness of distributions of solutions.

    Lemma 4.4. Suppose (2.1)–(2.7) and (3.1) hold. If φ(s)L2p(Ω,C([ρ,0],V)), then the solution u of (1.1) and (1.2) satisfies, for any t>r0,

    E[u(t)u(r)2p]M6(|tr|p+|tr|2p), (4.35)

    where M6 is a positive constant depending on φ, but independent of t and r.

    Proof. By (4.34), we get for t>r0,

     u(t)=eA(tr)u(r)treA(ts)(1+iμ)|u(s)|2βu(s)ds   +treA(ts)G(,u(sρ))ds+k=1treA(ts)(σ1,k+κσ2,k(u(s)))dWk(s). (4.36)

    Then we infer

     u(t)u(r)2p52p4[(eA(tr)I)u(r)2p+treA(ts)(1+iμ)|u(s)|2βu(s)ds2p   +treA(ts)G(,u(sρ))ds2p+k=1treA(ts)(σ1,k+κσ2,k(u(s)))dWk(s)2p]. (4.37)

    Taking the expectation of (4.36), we have for all t>r0,

     E[u(t)u(r)2p]52p4E[(eA(tr)I)u(r)2p]+52p4E[treA(ts)(1+iμ)|u(s)|2βu(s)ds2p]   +52p4E[treA(ts)G(,u(sρ))ds2p]+52p4E[k=1treA(ts)(σ1,k+κσ2,k(u(s)))dWk(s)2p]. (4.38)

    For the first term on the right-hand side of (4.38), by Theorem 1.4.3 in [32], we find that there exists a positive number C0 depending on ϱ such that for all t>r0,

    52p4E[(eA(tr)I)u(r)2p]C0(tr)pE[u(r)2pCV].

    Applying Lemmas 3.5 and 4.3, we obtain for all t>r0,

    52p4E[(eA(tr)I)u(r)2p]C1(tr)p. (4.39)

    For the second term on the right-hand side of (4.38), by the contraction property of eAt, we infer that for all t>r0,

       E[treA(ts)(1+iμ)|u(s)|2βu(s)ds2p](1+μ2)pE[(tr|u(s)|2β+1ds)2p](1+μ2)pE[(tr|u(s)|2β+12pds)](tr)2p1(1+μ2)psups0E[(u(s)2(2β+1)L2(2β+1))p](tr)2p.

    We deduce the estimate sups0E[(u(s)2(2β+1)L2(2β+1))p]M(E[φ2CV]+1) similarly to Lemma 3.5 together with Lemma 3.3 in [1]. Hence, the second term on the right-hand side of (4.38) can be estimated by

    E[treA(ts)(1+iμ)|u(s)|2βu(s)ds2p]C2(tr)2p. (4.40)

    For the third term on the right-hand side of (4.38), by the contraction property of eAt and (2.1) and Lemma 3.5, we deduce that for all t>r0,

        52p4E[treA(ts)G(,u(sρ))ds2p]52p4E[(trG(,u(sρ))ds)2p]52p4E[(tr(h+au(sρ))ds)2p]52p4E[(tr(h+au(sρ))2pds)](tr)2p1102p8(tr)2p1tr(h2p+a2pE[u(sρ)2p])ds102p8(h2p+a2psuptρE[u(s)2p])(tr)2pC3(tr)2p. (4.41)

    For the forth term the right-hand side of (4.38), from the BDG inequality, the contraction property of eAt, (2.6) H¨older's inequality and Lemma 3.5, we deduce

    52p4E[k=1treA(ts)(σ1,k+κσ2,k(u(s)))dWk(s)2p]52p4C4E[(trk=1eA(ts)(σ1,k+κσ2,k(u(s)))2ds)p]52p4C4E[(trk=12(σ1,k2+κσ2,k(u(s))2)ds)p]52p4C4E[(trk=12(σ1,k2+2κ2β2k+2κ2Lγ2k(u(s))2)ds)p]52p2C4E[(2k=1(σ1,k2+2κ2β2k)(tr)+4k=1κ2Lγ2ktru(s)2ds)p]102p8C4(k=1(σ1,k2+2κ2β2k)p(tr)p+102p8C4(2k=1κ2Lγ2k)pE[(tru(s)2ds)p]102p8C4(k=1(σ1,k2+2κ2β2k)p(tr)p    +102p8C4(2k=1κ2Lγ2k)p(tr)p1trE[u(s)2p]ds102p8C4(k=1(σ1,k2+2κ2β2k)p(tr)p    +102p8C4(2k=1κ2Lγ2k)p(tr)psups0E[u(s)2p]C5(tr)p. (4.42)

    Therefore, from (4.38)–(4.42), we obtain there exists C6>0 independent of t and r, such that for all t>r0,

    E[u(t)u(r)2p]C6(|tr|p+|tr|2p).

    The proof is complete.

    In this section, we first recall the definition of invariant measure and transition operator. Then we construct a compact subset of C([ρ,0];H) in order to prove the tightness of the sequence of invariant measure mk on C([ρ,0];H).

    Recall that for any initial time t0 and every Ft0-measurable function φ(s)L2(Ω,C([ρ,0],H)), problems (1.1) and (1.2) has a unique solution u(t;t0,φ) for t[t0ρ,). For convenience, given tt0 and Ft0-measurable function φ(s)L2(Ω,C([ρ,0],H)), the segment of u(t;t0,φ) on [tρ,t] is written as

    ut(t0,φ)(s)=u(t+s;t0,φ) for every s[ρ,0].

    Then ut(t0,φ)L2(Ω,C([ρ,0],H)) for all tt0. We introduce the transition operator for (1.1). If ϕ(s):C([ρ,0],H)R is a bounded Borel function, then for initial time r with 0rt and φ(s)C([ρ,0],H), we write

    (pr,tϕ)(φ)=E[ϕ(ut(r,φ))].

    Particularly, for ΓB(C([ρ,0],H)), 0rt and φC([ρ,0],H), we have

    p(r,φ;t,Γ)=(pr,t1Γ)(φ)=P{ωΩ|ut(r,φ)Γ},

    where 1Γ is the characteristic function of Γ. Then p(r,φ;t,) is the distribution of ut(0,φ) in C([ρ,0],H). In the following context, we will write p0,t as pt.

    Recall that a probability measure M on C([ρ,0],H) is called an invariant measure, if for all t0 and every bounded and continuous function ϕ:C([ρ,0];H)R,

    C([ρ,0];H)(ptϕ)(φ)dM(φ)=C([ρ,0];H)ϕ(φ)dM(φ),  for all t0.

    According to [33], we infer that the transition operator {pr,t}0rt has the following properties.

    Lemma 5.1. Suppose (2.1)–(2.7) and (4.1)–(4.3) hold. One has

    (a) The family {pr,t}0rt is Feller; that is, if ϕ:C([ρ,0],H)R is bounded and continuous, then for any 0rt, the function pr,tϕ:C([ρ,0],H)R is also bounded and continuous.

    (b) The family {pr,t}0rt is homogeneous (in time); that is, for any 0rt,

    p(r,φ;t,)=p(0,φ;tr,),φC([ρ,0],H).

    (c) Given r0 and φC([ρ,0],H), the process {ut(r,φ)}tr is a C([ρ,0],H)-valued Markov process. Consequently, if ϕ:C([ρ,0],H)R is a bounded Borel function, then for any 0srt, P-almost surely,

    (ps,tϕ)(φ)=(ps,r(pr,tϕ))(φ),φC([ρ,0],H),

    and the Chapman-Kolmogorov equation is valid:

    p(s,φ;t,Γ)=C([ρ,0],H)p(s,φ;r,dy)p(r,y;t,Γ),

    for any φC([ρ,0],H) and ΓB(C([ρ,0],H)).

    Now, we establish the existence of invariant measures of problems (1.1) and (1.2).

    Theorem 5.2. Suppose (2.1)–(2.7) and (4.1)–(4.3) hold. Then (1.1) and (1.2) processes an invariant measure on C([ρ,0],H).

    Proof. We employ Krylov-Bogolyubov's method to the solution u(t,0,0) of problems (1.1) and (1.2), where the initial condition φ0 at the initial time 0. Because of this particular φC([ρ,0],V)C([ρ,0],H), we know that all results obtained in the previous Sections 3 and 4 are valid. For simplicity, the solution u(t,0,0) is written as u(t) and the segment ut(0,0) as ut. For kN+, we set

    Mk=1kk+ρρp(0,0;t,)dt. (5.1)

    Step 1. We prove the tightness of {Mk}k=1 in C([ρ,0],H). Applying Lemmas 3.2 and 4.2, we get that there exists C1>0 such that for all tρ,

    E[supρs0ut(s)2V]C1. (5.2)

    By (5.2) and Chebyshev's inequality, we have that for all tρ,

    P({supρs0ut(s)VR})1R2E[supρs0ut(s)2V]C1R20  as  R,

    and hence for every ε>0, there exists R1=R1(ε)>0 such that for all tρ,

    P({supρs0ut(s)VR1})13ε. (5.3)

    By Lemma 4.4, we get that there exists C2>0 such that for all tρ and r,s[ρ,0],

    E[ut(r)ut(s)2p]C2(1+|rs|p)|rs|p,

    and hence for all tρ and r,s[ρ,0],

    E[ut(r)ut(s)2p]C2(1+ρp)|rs|p. (5.4)

    Since p2, applying (5.4) and the usual technique of dyadic division, we obtain that there exists R2=R2(ε)>0 such that for all tρ,

    P({supρsr0ut(r)ut(s)|rs|p14pR2})113ε. (5.5)

    By Lemma 3.4, we get that for given ε>0 and mN, there exists an integer nm=nm(ε,m)1 such that for all tρ,

    E[supρs0|x|nm|u(t+s,x)|2dx]ε22m+2,

    which implies that for all tρ and mN,

    P({supρs0|x|nm|u(t+s,x)|2dx12m})2mE[supρs0|x|nm|u(t+s,x)|2dx]ε2m+2. (5.6)

    By (5.6), we infer that for all tρ,

    P(m=1{supρs0|x|nm|u(t+s,x)|2dx12m})k=1ε2m+214ε,

    and hence for all tρ,

    P({supρs0|x|nm|u(t+s,x)|2dx12m for all mN})114ε. (5.7)

    Let

    M1,ε={ζ:[ρ,0]V,supρs0ζ(s)VR1(ε)}, (5.8)
    M2,ε={ζC([ρ,0],H):supρsr0ζ(r)ζ(s)|rs|ϱ14ϱR2(ε)}, (5.9)
    M3,ε={ζC([ρ,0],H):supρs0|x|nm|ζ(s,x)|2dx12m for all mN}, (5.10)

    and

    Mε=M1,εM2,εM3,ε. (5.11)

    From (5.3), (5.5) and (5.7)–(5.11), we obtain that for all tρ,

    P(utMε)>1ε. (5.12)

    By (5.1) and (5.12), we deduce that for all kN,

    Mk(Mε)>1ε. (5.13)

    Next, we prove the set Mε is precompact in C([ρ,0],H). First, we prove for every s[ρ,0] the set {ζ(s):ζMε} is a precompact subset of H. By (5.8) and (5.11), we obtain that for every s[ρ,0], the set {ζ(s):ζMε} is bounded in V. Let Qm0={xRn:|x|<nm0}. Then we get that the set {ζ(s)|Qm0:ζMε} is bounded in H1(Qm0) and hence precompact in L2(Qm0) due to compactness of the embedding H1(Qm0)L2(Qm0). This implies that the set {ζ(s)|Qm0:ζMε} has a finite open cover of balls with radius 12δ in L2(Qm0). Note that for every δ>0, there exists m0=m0(δ)N such that for all ζMε,

    |x|nm0|ζ(s,x)|2dx12m0<δ28. (5.14)

    Hence, by (5.14), the set {ζ(s):ζMε} has a finite open cover of balls with radius 12δ in L2(Rn). Since δ>0 is arbitrary, we obtain that the set {ζ(s):ζMε} is percompact in H. Then from (5.9) and (5.11), we obtain that Mε is equicontinuous in C([ρ,0],H). Therefore, by the Ascoli-Arzelˊa theorem we deduce that Mε is precompact in C([ρ,0],H), which along with (5.13) shows that {mk}k=1 is tight on C([ρ,0],H).

    Step 2. We prove the existence of invariant measures of problems (1.1) and (1.2). Since the sequence {Mk}k=1 is tight on C([ρ,0];H), there exists a probability measure m on C([ρ,0];H), we take a subsequence of {Mk}(not rebel) such that Mkm,  as  k. In the following, we prove M is an invariant measure of (1.1) and (1.2). Applying (5.1) and the Chapman-Kolmogorov equation, we obtain that for every t0 and every ϕ:C([ρ,0];H)R,

        C([ρ,0];H)ϕ(v)dM(v)=limk1kk+ρρ(C([ρ,0];H)ϕ(v)p(0,0;s,dv))ds=limk1kk+ρtρt(C([ρ,0];H)ϕ(v)p(0,0;s+t,dv))ds=limk1kk+ρρ(C([ρ,0];H)ϕ(v)p(0,0;s+t,dv))ds=limk1kk+ρρ(C([ρ,0];H)(C([ρ,0];H)ϕ(v)p(s,φ;s+t,dv))p(0,0;s,dφ))ds=limk1kk+ρρ(C([ρ,0];H)(C([ρ,0];H)ϕ(v)p(0,φ;t,dv))p(0,0;s,dφ))ds=C([ρ,0];H)(C([ρ,0];H)ϕ(v)p(0,φ;t,dv))dM(φ)=C([ρ,0];H)(p0,tϕ)(φ)dM(φ),

    which completes the proof.

    The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

    The authors are grateful to the anonymous referees whose suggestions have in our opinion, greatly improved the paper. This work is partially supported by the NSF of Shandong Province (No. ZR 2021MA055) and USA Simons Foundation (No. 628308).

    The authors declare there is no conflict of interest.



    [1] H. Lu, P. W. Bates, S. Lü, M. Zhang, Dynamics of the 3-D fractional complex Ginzburg-Landau equation, J. Differ. Equations, 259 (2015), 5276–5301. https://doi.org/10.1016/j.jde.2015.06.028 doi: 10.1016/j.jde.2015.06.028
    [2] V. E. Tarasov, G. M. Zaslavsky, Fractional Ginzburg-Landau equation for fractal media, Physica A, 354 (2005), 249–261. https://doi.org/10.1016/j.physa.2005.02.047 doi: 10.1016/j.physa.2005.02.047
    [3] V. E. Tarasov, G. M. Zaslavsky, Fractional dynamics of coupled oscillators with long-range interaction. Chaos, 16 (2006), 023110. https://doi.org/10.1063/1.2197167 doi: 10.1063/1.2197167
    [4] H. Weitzner, G. M. Zaslavsky, Some applications of fractional derivatives, Commun. Nonlinear Sci. Numer. Simul., 8 (2003), 273–281.
    [5] A. I. Saichev, G. M. Zaslavsky, Fractional kinetic equations: solutions and applications, Chaos, 7 (1997), 753–764. https://doi.org/10.1063/1.166272 doi: 10.1063/1.166272
    [6] M. F. Shlesinger, G. M. Zaslavsky, J. Klafter, Strange kinetics, Nature, 363 (1993), 31–37. https://doi.org/10.1038/363031a0 doi: 10.1038/363031a0
    [7] Y. Sire, E. Valdinoci, Fractional Laplacian phase transitions and boundary reactions: a geometric inequality and a symmetry result, J. Funct. Anal., 256 (2009), 1842–1864. https://doi.org/10.1016/j.jfa.2009.01.020 doi: 10.1016/j.jfa.2009.01.020
    [8] G. M. Zaslavsky, Chaos, fractional kinetics, and anomalous transport, Phys. Rep., 371 (2002), 461–580. https://doi.org/10.1016/S0370-1573(02)00331-9 doi: 10.1016/S0370-1573(02)00331-9
    [9] G. M. Zaslavsky, Hamiltonian Chaos and Fractional Dynamics, Oxford University Press, Oxford, 2005.
    [10] G. M. Zaslavsky, M. Edelman, Weak mixing and anomalous kinetics along filamented surfaces, Chaos, 11 (2001), 295–305. https://doi.org/10.1063/1.1355358 doi: 10.1063/1.1355358
    [11] C. Guo, J. Shu, X. Wang, Fractal dimension of random attractors for non-autonomous fractional stochastic Ginzburg-Landau equations, Acta Math. Sin. Engl. Ser., 36 (2020), 318–336. https://doi.org/10.1007/s10114-020-8407-4 doi: 10.1007/s10114-020-8407-4
    [12] H. Lu, P. W. Bates, J. Xin, M. Zhang, Asymptotic behavior of stochastic fractional power dissipative equations on Rn, Nonlinear Anal., 128 (2015), 176–198. https://doi.org/10.1016/j.na.2015.06.033 doi: 10.1016/j.na.2015.06.033
    [13] H. Lu, J. Li, M. Zhang, Stochastic dynamics of non-autonomous fractional Ginzburg-Landau equations on R3, Discrete Contin. Dyn. Syst. - Ser. B, 27 (2022), 6943–6968. https://doi.org/10.3934/dcdsb.2022028 doi: 10.3934/dcdsb.2022028
    [14] H. Lu, L. Wang, L. Zhang, M. Zhang, The asymptotic behavior of non-autonomous stochastic Ginzburg-Landau equations on thin domains, J. Appl. Anal. Comput., 11 (2021), 2313–2333. https://doi.org/10.11948/20200378 doi: 10.11948/20200378
    [15] H. Lu, L. Wang, M. Zhang, Dynamics of fractional stochastic Ginzburg-Landau equation on unbounded domains driven by nonlinear noise, Mathematics, 10 (2022), 4485. https://doi.org/10.3390/math10234485 doi: 10.3390/math10234485
    [16] H. Lu, M. Zhang, Dynamics of non-autonomous fractional Ginzburg-Landau equations driven by colored noise, Discrete Contin. Dyn. Syst. - Ser. B, 25 (2020), 3553–3576. https://doi.org/10.3934/dcdsb.2020072 doi: 10.3934/dcdsb.2020072
    [17] J. Shu, X. Huang, J. Zhang, Asymptotic behavior for non-autonomous fractional stochastic Ginzburg-Landau equations on unbounded domains, J. Math. Phys., 61 (2020), 072704. https://doi.org/10.1063/1.5143404 doi: 10.1063/1.5143404
    [18] Y. Kuang, Delay Differential Equations: with Applications in Population Dynamics, Academic Press, New York, 1993.
    [19] T. Caraballo, M. Garrido-Atienza, T. Taniguchi, The existence and exponential behavior of solutions to stochastic delay evolution equations with a fractional Brownian motion, Nonlinear Anal. Theory Methods Appl., 74 (2011), 3671–3684. https://doi.org/10.1016/j.na.2011.02.047 doi: 10.1016/j.na.2011.02.047
    [20] X. Wang, K. Lu, B. Wang, Random attractors for delay parabolic equations with additive noise and deterministic nonautonomous forcing, SIAM J. Appl. Dyn. Syst., 14 (2015), 1018–1047. https://doi.org/10.1137/140991819 doi: 10.1137/140991819
    [21] F. Wu, G. Yin, H. Mei, Stochastic functional differential equations with infinite delay: existence and uniqueness of solutions, solution maps, Markov properties, and ergodicity, J. Differ. Equations, 262 (2017), 1226–1252. https://doi.org/10.1016/j.jde.2016.10.006 doi: 10.1016/j.jde.2016.10.006
    [22] J. Kim, On the stochastic Burgers equation with polynomical nonlinearity in the real line, Discrete Contin. Dyn. Syst. - Ser. B, 6 (2006), 835–866. https://doi.org/10.3934/dcdsb.2006.6.835 doi: 10.3934/dcdsb.2006.6.835
    [23] J. Kim, On the stochastic Benjamin-Ono equation, J. Differ. Equations, 228 (2006), 737–768. https://doi.org/10.1016/j.jde.2005.11.005 doi: 10.1016/j.jde.2005.11.005
    [24] O. Misiats, O. Stanzhytskyi, N. Yip, Existence and uniqueness of invariant measures for stochastic reaction-diffusion equations in unbounded domains, J. Theor. Probab., 29 (2016), 996–1026. https://doi.org/10.1007/s10959-015-0606-z doi: 10.1007/s10959-015-0606-z
    [25] Z. Brzezniak, E. Motyl, M. Ondrejat, Invariant measure for the stochastic Navier-Stokes equations in unbounded 2D domains, Ann. Prabab., 45 (2017), 3145–3201. https://doi.org/10.1214/16-AOP1133 doi: 10.1214/16-AOP1133
    [26] Z. Brzezniak, M. Ondrejat, J. Seidler, Invariant measures for stochastic nonlinear beam and wave equations, J. Differ. Equations, 260 (2016), 4157– 4179. https://doi.org/10.1016/j.jde.2015.11.007 doi: 10.1016/j.jde.2015.11.007
    [27] J. Kim, Periodic and invariant measures for stochastic wave equations, Electron. J. Differ. Equations, 5 (2004), 1–30.
    [28] E. D. Nezza, G. Palatucci, E. Valdinoci, Hitchhiker's guide to the fractional Sobolev spaces, Bull. Sci. Math., 136 (2012), 521–573. https://doi.org/10.1016/j.bulsci.2011.12.004 doi: 10.1016/j.bulsci.2011.12.004
    [29] L. Tang, Dynamical behavior and multiple optical solitons for the fractional Ginzburg-Landau equation with βderivative in optical fibers, Opt. Quantum Electron., 56 (2024), 175.
    [30] X. Mao, Stochastic Differential Equations and Applications, Second Edition, Woodhead Publishing Limited, Cambridge, 2011.
    [31] G. D. Prato, J. Zabczyk, Stochastic Equations in Infinite Dimension, Cambridge University Press, Cambridge, 1992.
    [32] D. Henry, Geometric Theory of Semilinear Parabolic Equations, Springer, New York, 1981. https://doi.org/10.1007/BFb0089647
    [33] Z. Chen, B. Wang, Invariant measures of fractional stochastic delay reaction-diffusion equations on unbounded domains, Nonlinearity, 34 (2021), 3969–4016. https://doi.org/10.1088/1361-6544/ac0125 doi: 10.1088/1361-6544/ac0125
  • Reader Comments
  • © 2024 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(1334) PDF downloads(131) Cited by(0)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog