Propagation dynamics of nonlocal dispersal equations with inhomogeneous bistable nonlinearity

  • Received: 01 November 2019 Revised: 01 July 2020 Published: 24 November 2020
  • Primary: 35K57, 35R20; Secondary: 92D25

  • This paper is concerned with the nonlocal dispersal equations with inhomogeneous bistable nonlinearity in one dimension. The varying nonlinearity consists of two spatially independent bistable nonlinearities, which are connected by a compact transition region. We establish the existence of a unique entire solution connecting two traveling wave solutions pertaining to the different nonlinearities. In particular, we use a "squeezing" technique to show that the traveling wave of the equation with one nonlinearity approaching from infinity, after going through the transition region, converges to the other traveling wave prescribed by the nonlinearity on the other side. Furthermore, we also prove that such an entire solution is Lyapunov stable.

    Citation: Shao-Xia Qiao, Li-Jun Du. Propagation dynamics of nonlocal dispersal equations with inhomogeneous bistable nonlinearity[J]. Electronic Research Archive, 2021, 29(3): 2269-2291. doi: 10.3934/era.2020116

    Related Papers:

    [1] Shao-Xia Qiao, Li-Jun Du . Propagation dynamics of nonlocal dispersal equations with inhomogeneous bistable nonlinearity. Electronic Research Archive, 2021, 29(3): 2269-2291. doi: 10.3934/era.2020116
    [2] Yang Yang, Yun-Rui Yang, Xin-Jun Jiao . Traveling waves for a nonlocal dispersal SIR model equipped delay and generalized incidence. Electronic Research Archive, 2020, 28(1): 1-13. doi: 10.3934/era.2020001
    [3] Shuxia Pan . Asymptotic spreading in a delayed dispersal predator-prey system without comparison principle. Electronic Research Archive, 2019, 27(0): 89-99. doi: 10.3934/era.2019011
    [4] Meng-Xue Chang, Bang-Sheng Han, Xiao-Ming Fan . Global dynamics of the solution for a bistable reaction diffusion equation with nonlocal effect. Electronic Research Archive, 2021, 29(5): 3017-3030. doi: 10.3934/era.2021024
    [5] Minzhi Wei . Existence of traveling waves in a delayed convecting shallow water fluid model. Electronic Research Archive, 2023, 31(11): 6803-6819. doi: 10.3934/era.2023343
    [6] Cui-Ping Cheng, Ruo-Fan An . Global stability of traveling wave fronts in a two-dimensional lattice dynamical system with global interaction. Electronic Research Archive, 2021, 29(5): 3535-3550. doi: 10.3934/era.2021051
    [7] Jianlong Han, Seth Armstrong, Sarah Duffin . An unconditionally stable numerical scheme for competing species undergoing nonlocal dispersion. Electronic Research Archive, 2024, 32(4): 2478-2490. doi: 10.3934/era.2024114
    [8] Léo Girardin, Danielle Hilhorst . Spatial segregation limit of traveling wave solutions for a fully nonlinear strongly coupled competitive system. Electronic Research Archive, 2022, 30(5): 1748-1773. doi: 10.3934/era.2022088
    [9] Yujia Xiang, Yuqi Jiao, Xin Wang, Ruizhi Yang . Dynamics of a delayed diffusive predator-prey model with Allee effect and nonlocal competition in prey and hunting cooperation in predator. Electronic Research Archive, 2023, 31(4): 2120-2138. doi: 10.3934/era.2023109
    [10] Xiao Su, Hongwei Zhang . On the global existence and blow-up for the double dispersion equation with exponential term. Electronic Research Archive, 2023, 31(1): 467-491. doi: 10.3934/era.2023023
  • This paper is concerned with the nonlocal dispersal equations with inhomogeneous bistable nonlinearity in one dimension. The varying nonlinearity consists of two spatially independent bistable nonlinearities, which are connected by a compact transition region. We establish the existence of a unique entire solution connecting two traveling wave solutions pertaining to the different nonlinearities. In particular, we use a "squeezing" technique to show that the traveling wave of the equation with one nonlinearity approaching from infinity, after going through the transition region, converges to the other traveling wave prescribed by the nonlinearity on the other side. Furthermore, we also prove that such an entire solution is Lyapunov stable.



    In the present paper, we focus on the nonlocal dispersal equation with spatially dependent bistable nonlinearity

    ut(x,t)=Ju(x,t)u(x,t)+f(x,u),x,tR, (1)

    where the nonlinearity f(x,u) and kernel function J(x) satisfy the following assumptions.

    (F) The nonlinearity f satisfies that

    {f(x,u)=f1(u)forx0,u[0,1],f(x,u)=f2(u)forxL,u[0,1],f2(u)f(x,u)f1(u),fu(x,u)<1forx[L,0],u[0,1], (2)

    where L>0, f1 and f2 are two given bistable nonlinearities satisfying that

    fiC1,1([0,1]),fi(0)=fi(1)=0,fi(0),fi(1)<0,fi<0in(0,θi),fi>0in(θi,1),0<θ1<θ2<1,and10fi(s)ds>0,i=1,2.

    (J) The kernel function JC1(R) satisfies that

    {J(x)=J(x),J(x)0,RJ(y)dy=1,R|J(y)|dy<,RJ(y)eλydy<for allλ>0. (3)

    A typical example is f(x,u)=u(1u)(a(x)u).

    The nonlocal dispersal equation has attracted so much attention because of its extensive applications to account for diffusion phenomena involving jumps in biology, physics and chemistry [1,12]. Traveling waves as one kind of special solutions with invariant profile and fixed speed as well as entire solutions have been adequately investigated. It is well-known that the importance of the study of entire solutions of reaction-diffusion equations (nonlocal dispersal equations) is frequently recalled in the literature. Since the pioneering works of Hamel and Nadirashvili [13,14], there have been tremendous advances in studying the existence of entire solutions for various models. In particular, when the nonlinearity is homogeneous (i.e. L=0 and f1=f2), the theories of traveling waves and entire solutions for nonlocal dispersal equation 1 with various types of nonlinearities have been well established, the related results refer to [2,6,7,8,9,15,18,23,24,26,27] and references therein. Specifically, when the kernel J(x) is compactly supported, Sun et al. [23] constructed a two-dimensional manifold of entire solutions which behave as two traveling wave solutions coming from both directions for bistable nonlocal dispersal equation.

    However, for the inhomogeneous nonlinearity, several works are devoted to transition fronts (see [3,19]) of nonlocal dispersal equations [17,20] and forced waves in shifting habitats [25]. Traveling waves and spreading speed of monostable nonlocal dispersal equations with space periodic nonlinearity were studied in [21]. Li et al. [16] further obtained the existence of entire solutions for space periodic nonlinearity. Particularly, Eberle [10,11] constructed a heteroclinic orbit connecting two traveling waves for bistable local dispersal equation 1 in cylinders. Meanwhile, Berestycki and Rodríguez [5] considered a bistable nonlocal dispersal equation with a gap in one dimension. As far as we know, there is no result on entire solutions for the nonlocal dispersal equation 1 with inhomogeneous nonlinearity f(x,u) satisfying (F).

    In this paper, we aim to construct an entire solution connecting traveling waves with the two different nonlinearities, motivated by [10,11]. By constructing suitable sub- and super-solutions, we establish the existence and uniqueness of the entire solution behaving as the traveling wave coming from one side and eventually going to the other side, which is different with the one constructed in [23] for the case L=0 and f2=f1. It should be pointed out that there are many differences to [4], though the similar method is applied. More precisely, since the nonlocal operator is not a compact operator as the Laplacian operator in [4], we have to establish the Lipschitz continuity of entire solutions in order to prove its uniqueness as in [23]. Compared to [23], we drop the assumption the kernel J(x) is compactly supported.

    The crucial part of this paper is to figure out the long time asymptotic behavior of the entire solution. Since the lack of compactness of the nonlocal operator, we can not use the Lyapunov function argument as [10,11] to show that the entire solution converges to a translation of the other traveling wave as time goes to positive infinity. However, inspired by the idea in [7], we use a "squeezing" technique to address this issue. Furthermore, we also apply sub- and super-solutions method with comparison principle to establish Lyapunov stability of the entire solution. This can be done because that we can obtain a positive estimate on the derivative of the entire solution with respect to t when it is not so close to 0 as well as 1 and t is large enough.

    Now we state the main results of this paper as follows.

    Theorem 1.1. Let assumptions (F) and (J) hold. Then there exists a unique entire solution u(x,t) of 1 with 0<u(x,t)<1, ut(x,t)>0 for all tR and xR such that

    u(x,t)ϕ1(x+c1t)0astuniformly inxR

    and

    u(x,t)ϕ2(x+c2t+β)0ast+uniformly inxR

    for some βR, where (ϕi,ci)(i=1,2) are the traveling wave solutions solving

    {Jϕiϕi+fi(ϕi)ciϕi=0,ϕi()=0,ϕi(+)=1,0<ϕi<1. (4)

    Theorem 1.2. The entire solution u(x,t) constructed in Theorem 1.1 is Lyapunov stable in the following sense: For any given ϵ>0, there exists δ>0 such that for any uniformly continuous v0(x)[0,1] with supxR|v0(x)u(x+a,t0)|<δ, the solution v(x,t,v0(x)) of 1 with initial value v0(x) satisfies

    |v(x,t,v0)u(x+a,t+t0)|<ϵ

    for xR and t0, where a,t0R are two constants.

    The rest of this paper is organized as follows. In Section 2, we recall some results of the bistable traveling waves for homogeneous nonlinearities and prove the comparison principle for 1. Section 3 is devoted to constructing the unique entire solution. In Section 4, we study the asymptotic profile as time goes to positive infinity. Finally, we establish Lyapunov stability of the entire solution in Section 5.

    In this section, some known results on the traveling waves of 4 are outlined and the comparison principle is established.

    It follows from Theorem 2.7 in [2] and Theorem 2.7 in [23] that 4 admits a solution ϕ1(z) satisfying

    {β0eλ0zϕ1(z)α0eλ0z,z0,β1eλ1z1ϕ1(z)α1eλ1z,z>0, (5)

    where α0,α1,β0 and β1 are positive constants, λ0 and λ1 are the positive roots of

    cλ0=RJ1(y)eλ0ydy1+f1(0)andcλ1=RJ1(y)eλ1ydy1+f1(1),

    respectively. Moreover, we have

    {˜β0eλ0zϕ1(z)˜α0eλ0z,z0,˜β1eλ1zϕ1(z)˜α1eλ1z,z>0 (6)

    for some constants ˜α0,˜α1,˜β0 and ˜β1>0. At last, note that f1C1,1([0,1]), there exists some Lf>0 such that

    |f1(u+v)f1(u)f1(v)|Lfuvfor0u,v1.

    We show that the following comparison principle holds by a contradiction argument.

    Proposition 1. Suppose that assumptions (F) and (J) hold. Furthermore, let u(x,t), v(x,t) be continuous and bounded functions on R×[0,T] for some T>0 satisfy

    {u(x,t)t(RJ(xy)[u(y,t)u(x,t)]dy)f(x,u(x,t))0,(x,t)R×(0,T],u(x,0)0,xR,
    {v(x,t)t(RJ(xy)[v(y,t)v(x,t)]dy)f(x,v(x,t))0,(x,t)R×(0,T],v(x,0)0,xR,

    respectively. Then,

    u(x,t)v(x,t)inR×[0,T].

    Furthermore, if u(x,0)v(x,0) for xR, then u(x,t)>v(x,t) for xR,t[0,T].

    Proof of Proposition 1. Let ˉw(x,t)=u(x,t)v(x,t). In fact, it is sufficient to show ˉw(x,t)0 for (x,t)R×[0,ϵ0T] with ϵ0(0,1). Otherwise, suppose that infR×[0,ϵ0T]ˉw(x,t)<0. Denote ˇw(x,t)=eZtˉw(x,t), where Z=fu(x,u)+1. It follows that infR×[0,ϵ0T]ˇw(x,t)<0. Then there exists a sequence (xn,tn)R×(0,ϵ0T] such that

    limn+ˇw(xn,tn)=infR×[0,ϵ0T]ˇw(x,t)<0.

    Observe that

    ˇwt(x,t)=Zˇw(x,t)+eZtˉwt(x,t)Zˇw(x,t)+eZtRJ(xy)[ˉw(y,t)ˉw(x,t)]dy+eZt[f(x,u(x,t))f(x,v(x,t))]=Zˇw(x,t)+RJ(xy)[ˇw(y,t)ˇw(x,t)]dy+fu(x,uθ(x,t))ˇw(x,t).

    where uθ(x,t) is between v(x,t) and u(x,t). This implies that

    ˇw(xn,tn)ˇw(xn,0)tn0[Jˇw(xn,s)ˇw(xn,s)+Zˇw(xn,s)+fu(x,uθ(x,s))ˇw(xn,s)]dstn0[Jˇw(xn,s)+(fu(x,u)fu(x,uθ(x,s)))infR×[0,ϵ0T]ˇw(x,t)]ds.

    Letting n converge to infinity, we have

    infR×[0,ϵ0T]ˇw(x,t)(1+fu(x,u)minR×[0,1]|fu(x,uθ(x,s))|)ϵ0TinfR×[0,ϵ0T]ˇw(x,t).

    Since infR×[0,ϵ0T]ˇw(x,t)<0, we can choose ϵ0 sufficiently small such that

    infR×[0,ϵ0T]ˇw(x,t)(1+fu(x,u)minR×[0,1]|fu(x,uθ(x,s))|)ϵ0TinfR×[0,ϵ0T]ˇw(x,t)>infR×[0,ϵ0T]ˇw(x,t).

    Thus, we get a contradiction. Therefore, we obtain u(x,t)v(x,t) for xR and t[0,T]. The remaining part of this proposition can be proved similarly by replacing the auxiliary function ˇw(x,t)=eZtˉw(x,t) with ˇw(x,t)=eZtˉw(x,t)ϵt for small ϵ>0. In fact, we can similarly obtain

    infR×[0,ϵ0T]ˇw(x,t)(1+fu(x,u)minR×[0,1]|fu(x,uθ(x,s))|ϵinfR×[0,ϵ0T]ˇw(x,t))ϵ0TinfR×[0,ϵ0T]ˇw(x,t).

    Here, we can also choose ϵ0 such that

    (1+fu(x,u)minR×[0,1]|fu(x,uθ(x,s))|ϵinfR×[0,ϵ0T]ˇw(x,t))ϵ0T<1.

    This finishes the proof.

    In this section, we focus on the construction of the entire solution which behaves like a traveling wave approaching from infinity. The main idea is to establish suitable sub- and super-solutions, which are defined as follows.

    W(x,t)={ϕ1(x+c1tξ(t))ϕ1(x+c1tξ(t)),x0,0,x<0,

    and

    W+(x,t)={ϕ1(x+c1t+ξ(t))+ϕ1(x+c1t+ξ(t)),x0,2ϕ1(c1t+ξ(t)),x<0,

    here ξ(t) is the solution of the following equation

    ˙ξ(t)=Meλ(c1t+ξ),t<T,ξ()=0,

    where M,λ and T are positive constants to be specified later. A direct calculation yields that

    ξ(t)=1λln11c11Meλc1t.

    For the function ξ(t) to be defined, one must have 1c11Meλc1t>0. Besides, there is

    c1t+ξ(t)0for<tT,

    where T:=1λc1lnc1+Mc1.

    Now we verify that W are sub- and super-solutions of 1 for tT. Define

    Lu=ut(Juu)f(x,u).

    We first deal with the sub-solution W(x,t). Since it is obvious for x<0, we only consider the case x0. A straightforward calculation yields that

    LW(x,t)=(c1˙ξ(t))[ϕ1(x+c1tξ(t))ϕ1(x+c1tξ(t))]+0J(xy)[ϕ1(y+c1tξ(t))ϕ1(y+c1tξ(t))]dy+[ϕ1(x+c1tξ(t))ϕ1(x+c1tξ(t))]f(x,W)=˙ξ(t)[ϕ1(x+c1tξ(t))ϕ1(x+c1tξ(t))]+0J(xy)[ϕ1(y+c1tξ(t))ϕ1(y+c1tξ(t))]dy+f1(ϕ1(x+c1tξ(t)))f1(ϕ1(x+c1tξ(t)))f(x,W).

    Recall that f(x,u)=f1(u) for x0 and ϕ1>0, it follows that

    LW(x,t)˙ξ(t)[ϕ1(x+c1tξ(t))ϕ1(x+c1tξ(t))]+f1(ϕ1(x+c1tξ(t)))f1(ϕ1(x+c1tξ(t)))f1(ϕ1(x+c1tξ(t))ϕ1(x+c1tξ(t))).

    Then we continue to show LW(x,t)0 in two cases 0<x<c1t+ξ(t) and xc1t+ξ(t).

    Case 1. For 0<x<c1t+ξ(t), similar to [4], we have

    ϕ1(x+c1tξ(t))ϕ1(x+c1tξ(t))>m[ϕ1(x+c1tξ(t))ϕ1(x+c1tξ(t))]

    for some m>0. Consequently,

    LW(x,t)˙ξ(t)[ϕ1(x+c1tξ(t))ϕ1(x+c1tξ(t))]+Lfϕ1(x+c1tξ(t)))[ϕ1(x+c1tξ(t))ϕ1(x+c1tξ(t))][ϕ1(x+c1tξ(t))ϕ1(x+c1tξ(t))][mMeλ(c1t+ξ(t))+Lfα0eλ0(x+c1tξ(t))]=[ϕ1(x+c1tξ(t))ϕ1(x+c1tξ(t))]eλ(c1t+ξ(t))[mM+Lfα0eλ0(x2ξ(t))].

    Therefore, if we choose MLfα0m then LW(x,t)0.

    Case 2. For xc1t+ξ(t), if λ0λ1, then we obtain

    LW(x,t)Meλ(c1t+ξ(t))[~β1eλ1(x+c1tξ(t))˜αeλ0(x+c1tξ(t))]+Lfα0eλ0(x+c1tξ(t))[α1eλ1(x+c1tξ(t))β0eλ0(x+c1tξ(t))]eλ0(x+c1tξ(t))[M~β1e(λ0λ1)xe(λ0+λ1λ)c1t+(λ+λ0+λ1)ξ(t))M˜αeλ(c1t+ξ(t))Lfα0].

    Since c1t+ξ(t)<0 for t<T and λ<min{λ0,λ1}, LW(x,t)0 provided M2Lfα0~β1.

    For λ0<λ1, which means f1(0)>f1(1), we have

    LW(x,t)Meλ(c1t+ξ(t))ϕ1(x+c1tξ(t))+f1(ϕ1(x+c1tξ(t)))f1(ϕ1(x+c1tξ(t)))f1(ϕ1(x+c1tξ(t))ϕ1(x+c1tξ(t))).

    Moreover, if x+c1tξ(t)>L for some L1 such that ϕ1(x+c1tξ(t))[1σ,1] and ϕ1(x+c1tξ(t))[0,σ]. It follows that

    LW(x,t)Meλ(c1t+ξ(t))ϕ1(x+c1tξ(t))+[f1(1)f(0)]ϕ1(x+c1tξ(t))+o(ϕ1(x+c1tξ(t)))eλ0(x+c1tξ(t))[M˜α0eλ(c1t+ξ(t))+(f1(1)f(0))o(1)].

    Since t<T such that c1t+ξ(t)1, we have LW(x,t)0. Finally, if 0<x+c1tξ(t)L, then, from the above case λ0λ1, there holds

    LW(x,t)eλ0(x+c1tξ(t))[M~β1e(λ0λ1)xe(λ0+λ1λ)c1t+(λ+λ0+λ1)ξ(t))M˜αeλ(c1t+ξ(t))Lfα0]eλ0(x+c1tξ(t))[M~β1e(λ0λ1)Le(λ2λ0)c1t+(λ+2λ0)ξ(t))M˜αeλ(c1t+ξ(t))Lfα0].

    Thanks to t<T, λ<min{λ0,λ1} and M>2Lfα0~β1e(λ1λ0)L, we obtain

    LW(x,t)0.

    We intend to testify the super-solution W+(x,t) in two steps.

    Step 1. W+(x,t)=2ϕ1(c1t+ξ(t)) for x<0. Following a direct calculation, we have

    LW+(x,t)=2(c1+˙ξ(t))ϕ1(c1t+ξ(t))RJ(xy)W+(y,t)dy+2ϕ1(c1t+ξ(t))f(x,2ϕ1(c1t+ξ(t)))2(c1+˙ξ(t))ϕ1(c1t+ξ(t))f1(2ϕ1(c1t+ξ(t)))+0J(xy)[ϕ1(y+c1t+ξ(t))+ϕ1(y+c1t+ξ(t))]dy.

    Denote

    I=+0J(xy)[ϕ1(y+c1t+ξ(t))+ϕ1(y+c1t+ξ(t))]dy.

    Then, it follows that

    I=c1tξ(t)0J(xy)[ϕ1(y+c1t+ξ(t))+ϕ1(y+c1t+ξ(t))]dy++c1tξ(t)J(xy)[ϕ1(y+c1t+ξ(t))+ϕ1(y+c1t+ξ(t))]dyα0eλ0(c1t+ξ(t))+0J(xy)[eλ0y+eλ0y]dy+2+c1tξ(t)J(xy)dy.

    Let

    Jλ0=+0J(y)[eλ0y+eλ0y]dy.

    It follows from assumption (J) that there exists KJ>0 such that

    J(x)KJe2λ0xforx0.

    Furthermore, there hold

    IJλ0α0eλ0(c1t+ξ(t))+2KJ+c1tξ(t)e2λ0(xy)dyJλ0α0eλ0(c1t+ξ(t))+2KJλ0e2λ0(c1t+ξ(t)).

    Therefore, according to t<T such that ϕ1(c1t+ξ(t))1, we get

    LW+(x,t)2(c1+˙ξ(t))ϕ1(c1t+ξ(t))Jλ0α0eλ0(c1t+ξ(t))2KJλ0e2λ0(c1t+ξ(t))f1(2ϕ1(c1t+ξ(t)))(2c1~β0α0Jλ02f1(0)β0o(1))eλ0(c1t+ξ(t))+(2M~β02KJλ0)e2λ0(c1t+ξ(t)).

    Since

    M>KJ~β0λ0andJλ0<2c1~β0+2|f1(0)|β0α0,

    consequently, we have LW+(x,t)0 for t<T.

    Step 2. For x0, W+(x,t)=ϕ1(x+c1t+ξ(t))+ϕ1(x+c1t+ξ(t)). Then we can obtain

    LW+(x,t)=(c1+˙ξ(t))[ϕ1(x+c1t+ξ(t))+ϕ1(x+c1t+ξ(t))]RJ(xy)W+(y,t)dy+[ϕ1(x+c1t+ξ(t))+ϕ1(x+c1t+ξ(t))]f1(W+)=˙ξ(t)[ϕ1(x+c1t+ξ(t))+ϕ1(x+c1t+ξ(t))]+0J(xy)[ϕ1(y+c1t+ξ(t))+ϕ1(y+c1t+ξ(t))2ϕ1(c1t+ξ(t))]dy+f1(ϕ1(x+c1t+ξ(t)))+f1(ϕ1(x+c1t+ξ(t)))f1(ϕ1(x+c1t+ξ(t))+ϕ1(x+c1t+ξ(t))).

    Denote

    II=0J(xy)[ϕ1(y+c1t+ξ(t))+ϕ1(y+c1t+ξ(t))2ϕ1(c1t+ξ(t))]dy.

    We consider two cases.

    Case 1. 0xc1tξ(t). Without loss of generality, let ϕ1(0)=θ1. Since ϕ1>0 and t<T such that ϕ1(c1t+ξ(t))<12θ1, we have

    II=0c1t+ξ(t)J(xy)[ϕ1(y+c1t+ξ(t))+ϕ1(y+c1t+ξ(t))2ϕ1(c1t+ξ(t))]dy+c1t+ξ(t)J(xy)[ϕ1(y+c1t+ξ(t))+ϕ1(y+c1t+ξ(t))2ϕ1(c1t+ξ(t))]dy0c1t+ξ(t)J(xy)Cλ0eλ0(c1t+ξ(t))[eλ0y+eλ0y2(2+e(λ0+η)y+e(λ0+η)y)Cηeη(c1t+ξ(t))]dy(2+Jη)Cλ0Cηe(λ0+η)(c1t+ξ(t)).

    The second inequality follows from

    |ϕ1(x)Cλ0eλ0x|Cηe(λ0+η)xforx0and some0<η<λ0,

    which can be easily obtained by 5 and 6. As a consequence,

    LW+(x,t)Me(λ+λ0)(c1t+ξ(t))~β0(eλ0x+eλ0x)(2+Jη)Cλ0Cηe(λ0+η)(c1t+ξ(t))Lfα20e2λ0(c1t+ξ(t)))e(λ+λ0)(c1t+ξ(t))[2M~β0(2+Jη)Cλ0Cηe(ηλ)(c1t+ξ(t))Lfα20e(λ0λ)(c1t+ξ(t))]e(λ+λ0)(c1t+ξ(t))[2M~β0(2+Jη)Cλ0CηLfα20].

    The last inequality holds since λ<min{λ0,λ1,η}. Thus LW+(x,t)0, provided that

    M(2+Jη)Cλ0Cη+Lfα202~β0.

    Case 2. Here x>c1tξ(t). From Case 1, we know

    II0c1t+ξ(t)J(xy)Cλ0eλ0(c1t+ξ(t))(2+e(λ0+η)y+e(λ0+η)y)Cηeη(c1t+ξ(t))dyCλ0CηKJe(λ0+η)(c1t+ξ(t))0c1t+ξ(t)e2λ0(xy)(3+e(λ0+η)y)dyCλ0CηKJeλ0(x+c1t+ξ(t))λ0x(3λ0eη(c1t+ξ(t))+1η).

    Moreover, if λ1λ0, then

    LW+(x,t)Meλ(c1t+ξ(t))(~β1eλ1(x+c1t+ξ(t))+~β0eλ0(x+c1t+ξ(t)))Lfα0eλ0(x+c1t+ξ(t))Cλ0CηKJeλ0(x+c1t+ξ(t))(3λ0eη(c1t+ξ(t))+1η)eλ(c1t+ξ(t))+λ0(x+c1t+ξ(t))[M~β1e(λ0λ1)x(λ0+λ1)(c1t+ξ(t))3Cλ0CηKJλ0e(ηλ)(c1t+ξ(t))(Lfα0+Cλ0CηKJη)eλ(c1t+ξ(t))].

    Remember that λ<min{λ0,λ1,η} and t<T. As a matter of fact, take M sufficiently large such that

    M~β13Cλ0CηKJλ0(Lfα0+Cλ0CηKJη)0. (7)

    This yields that LW+(x,t)0.

    If λ0<λ1, we have f1(0)>f1(1). In addition, for x>c1tξ(t)+M with M1 such that ϕ1(x+c1t+ξ(t))[1σ,1], we have

    f1(ϕ1(x+c1t+ξ(t)))+f1(ϕ1(x+c1t+ξ(t)))f1(W+)=[f1(0)f1(1)]ϕ1(x+c1t+ξ(t))+o(ϕ1(x+c1t+ξ(t))).

    Meanwhile,

    LW+(x,t)eλ(c1t+ξ(t))+λ0(x+c1t+ξ(t))(M~β1e(λ0λ1)x(λ0+λ1)(c1t+ξ(t))3Cλ0CηKJλ0e(ηλ)(c1t+ξ(t)))+([f1(0)f1(1)]β0Cλ0CηKJηeλ0x)eλ0(x+c1t+ξ(t)).

    Then LW+(x,t)0 by 7 and x<c1t+ξ(t)M1 such that

    [f1(0)f1(1)]β0Cλ0CηKJηeλ0x0.

    When c1tξ(t)<xc1tξ(t)+M, it follows from the case λ0λ1 that

    LW+(x,t)eλ(c1t+ξ(t))+λ0(x+c1t+ξ(t))[M~β1e(λ0λ1)x(λ0+λ1)(c1t+ξ(t))3Cλ0CηKJλ0e(ηλ)(c1t+ξ(t))(Lfα0+Cλ0CηKJη)eλ(c1t+ξ(t))]eλ(c1t+ξ(t))+λ0(x+c1t+ξ(t))[M~β1e(λ0λ1)M2λ0(c1t+ξ(t))3Cλ0CηKJλ0e(ηλ)(c1t+ξ(t))(Lfα0+Cλ0CηKJη)eλ(c1t+ξ(t))].

    Similar to the case λ1λ0, we can obtain LW+(x,t)0 for t<T by letting M be sufficiently large.

    Proof of Theorem 1.1. Let un(x,t) be the unique solution of the Cauchy problem

    {(un)t(x,t)=Jun(x,t)un(x,t)+f(x,un),xR,t>n,un(x,n)=W(x,n),xR.

    Since W(x,n)=un(x,n)W+(x,n), the comparison principle yields for n>T

    W(x,t)un(x,t)W+(x,t)forxR,t[n,T].

    Moreover, since m:=minR×[0,1](1fu(x,u))>0 by the assumption fu(x,u)<1, it is not difficult to see that δun(x,t)=un(x+η,t)un(x,t) with ηR is a sub-solution of

    {v(t)=Lηmv(t),t>n,v(n)=Mη,

    where LR|J(y)|dy and M is some positive constant. It is then similar to the proof of Proposition 2.4 in [23] that un(x,t) and (un)t(x,t) are Lipschiz continuous with respect to x. Besides, a direct calculation gives that |un(x,t)|<C and |(un)t(x,t)|<C for some positive constant C. Then applying Arzela-Ascoli Theorem, there exists a subsequence of {un}+n=1 and {(un)t}+n=1, still denoted by {un}+n=1 and {(un)t}+n=1, such that {un}+n=1 and {(un)t}+n=1 converge to a function u(x,t) and ut(x,t). Now we have obtained the entire solution. In fact, the regularity of this entire solution is not so well. But we can show the entire solution is Lipschitz continuous in xR, which is important to show the uniqueness of the entire solution.

    Proposition 2. Let u(x,t) be the entire solution in Theorem 1.1. Then u(x,t) satisfies

    |u(x+η,t)u(x,t)|Mη,

    and

    |u(x+η,t)tu(x,t)t|M"η

    for M,M">0.

    The proof is similar to that of Proposition 3.1 in [23] by virtue of the fact fu(x,u)<1 for xR and u[0,1], so we omit it here.

    As for the uniqueness, since it is easy to see u(x,t)ϕ1(x+c1t) as t by the sub- and super-solutions established above, we can similarly show that the conclusion of Lemma 3.1 in [4] also holds true here. Then referring to the process of Section 3 in [4], we can obtain the uniqueness of the entire solution by constructing similar sub- and super-solutions. Thus we have shown the first part of Theorem 1.1.

    In this section, we are going to show the entire solution established in previous section converges to a shift of ϕ2(x+c2t) as t+. The main idea here is first to construct suitable sub-and super-solutions taking advantage of the traveling wave ϕ2(x+c2t) to get proper lower and upper bounds, which is a classical and regular method in study of traveling waves. Then we apply a "squeezing" technique which is used to prove the asymptotic stability of the traveling waves for nonlocal dispersal equation in [7] to obtain the convergence of the entire solution. Here, we assume that for some X>0,ω,σ>0 and

    ω<min{|f1(0)|4,|f2(0)|4,|f1(1)|4,|f2(1)|4},

    there hold

    ϕ1(x),ϕ2(x)σ2forxXandϕ1(x),ϕ2(x)1σ2forxX,

    besides,

    f1(s),f2(s)ωfors[0,σ][1σ,1].

    Theorem 4.1. Let assumptions (F) and (J) hold, and u(x,t) be the unique entire solution of 1 with 0<u(x,t)<1, ut(x,t)>0 for all tR and xR such that

    u(x,t)ϕ1(x+c1t)0astuniformly inxR. (8)

    Then

    u(x,t)ϕ2(x+c2t+β)0ast+uniformly inxR

    for some βR.

    Lemma 4.2. Suppose assumptions (F) and (J) hold. Then there exist T, β, δ,ω>0 such that

    u(x,t)max{ϕ2(x+c2tβ)δeωt,0}

    for xR and tT.

    Proof of Lemma 4.2. In fact, it follows from 8 that there exists tR such that

    |u(x,t)ϕ1(x+c1t)|ϵ2

    for any ϵ>0 and xR. Denote

    ˜u(x,t)=u(x,t+t)andu_(x,t)=max{ϕ2(ξ(x,t))Cveωt,0},

    where

    ξ(x,t)=x+c2(t+t)β0+CveωtCvfort0,xR.

    We first show ˜u(x,0)u_(x,0) for xR. Then there exists x1>0 such that

    ϕ1(x+c1t)1ϵ2andϕ2(x+c2t)1ϵ2forxx1.

    As ut>0 and ϕi>0 (i=1,2), we have

    1ϵ2u(x,t)1and1ϵ2ϕi1,i=1,2

    for tt,xx1. Then,

    |ϕiu(x,t)|ϵ,i=1,2fortt,xx1.

    Meanwhile, there exists x2R such that

    0<ϕi(x+cit)ϵ2,i=1,2forxx2.

    Particularly, there holds

    |u(x,t)ϕ2(x+c2t)|ϵforxR(x1,x2).

    In addition, minx[x1,x2]u(x,t)>0 since 0<u(x,t)<1. Then we can choose β0>0 such that

    ϕ2(x+c2tβ0)u(x,t)forx[x1,x2].

    From all the discussion above we obtain ˜u(x,0)u_(x,0). Next, we prove

    Lu_(x,t)=u_t(x,t)(Ju_(x,t)u_(x,t))f(x,u_)0fort0,xR.

    In fact, we only need to consider u_(x,t)=ϕ2(ξ(x,t))Cveωt because 0 is obviously a sub-solution of 1. Hence,

    Lu_(x,t)=(c2Cvωeωt)ϕ2(ξ(x,t))+CvωeωtRJ(xy)ϕ2(ξ(y,t))dy+ϕ2(ξ(x,t))f(x,u_(x,t))Cvωeωtϕ2(ξ(x,t))+f2(ϕ2(ξ(x,t)))f2(ϕ2(ξ(x,t))Cveωt)+Cvωeωt.

    In the following, we continue the proof in three cases.

    Case ξ(x,t)X. In this case,

    ϕ2(ξ)Cveωtϕ2(ξ)σ.

    Then, by ω<|f2(0)|4,

    Lu_(x,t)(f2(0)+2ω)Cveωt0.

    Case ξ(x,t)X. For this case, by 2Cv<σ,

    1σ1σ2Cvϕ2(ξ)Cveωtϕ2(ξ)<1.

    Then, since ω<|f2(1)|4,

    Lu_(x,t)(f2(1)+2ω)Cveωt0.

    Case ξ[X,X]. Let τ0=minξ[X,X]ϕ2(ξ). Then,

    Lu_(x,t)τ0Cvωeωt+f2Cveωt+Cvωeωt(f2Cv++Cvωτ0Cvω)eωt.

    Choose Cv=τ10ω10(f2Cv++Cvω). This gives Lu_(x,t)0. We finish the proof by letting β=β0+Cv and δ=Cv.

    Similarly, we can prove the following lemma.

    Lemma 4.3. Under the assumptions of Lemma 4.2, there exist T+>T, β+, C+,ω>0 such that

    min{ϕ1(x+c1t+β+)+C+eωt,1}u(x,t)

    for xR,tT+, where ω,T are defined in Lemma 4.2.

    Now we are in a position to establish an important lemma.

    Lemma 4.4. Under the assumptions of Lemma 4.2, there exist T>max{T+,T}, β,C,ω>0 such that

    min{ϕ2(x+c2t+β)+Ceωt,1}u(x,t)

    for xR,tT, where ω,T,T+ are defined in Lemmas 4.2 and 4.3.

    Proof of Lemma 4.4. Define ˆu(x,t)=u(x,t+T) and ¯u(x,t)=min{ϕ2(ξ+(x,t))+v(t),1}, where

    ξ+(x,t)=x+c2(t+T)+β++V(t),V(t)=KVt0v(s)ds,

    and

    v(t)=(ϵ1+C3)eωtC3eλ1c2t

    for 0<ϵ1,C3<σ4 and t0. Here, λ1 is the positive root of

    c2λ1=RJ(xy)eλ1ydy1+f2(1).

    Then

    v(0)=ϵ1,0<v(t)<σ2fort0

    and

    ˙v(t)=ωv(t)(wλ1c2)C3eλ1c2tfort0.

    Now we intend to show ˆu(x,0)¯u(x,0). Recall the fact that |1ϕ2(x)|C2eλ1x for xR,C2>0. Meanwhile, we assume that 0<ω<λ1c2. Additionally, by the assumptions of f(x,u), there exists C0>0 such that

    |f1(s)f2(s)|C0|1s|fors[0,1].

    Since ϕi()=0 and ϕi(+)=1, there exists x3<0 such that

    ϕ1(x+c1T+β+)<ϵ12andϕ2(x+c2T)<ϵ1forx<x3

    as well as x4>0 such that

    ϕ2(x+c2Tβ)1ϵ1andϕ2(x+c2T)1ϵ1

    for x>x4. Note that we can take T>0 sufficiently large such that

    max{C+,δ}eωT<ϵ12andXc2Tx0.

    Therefore, it is obvious that

    u(x,T)min{ϕ2(x+c2T)+ϵ1,1}¯u(x,0)forxR(x3,x4).

    Furthermore, since maxx[x3,x4]u(x,T)<1, one can choose β+>0 such that

    ϕ2(x+c2T+β+)u(x,T)forx[x3,x4].

    This yields that ˆu(x,0)¯u(x,0).

    In the following, we are going to prove

    L¯u(x,t)=¯ut(x,t)(J¯u(x,t)¯u(x,t))f(x,¯u)0fort0,xR.

    Here we only need to consider ¯u(x,t)=ϕ2(ξ+(x,t))+v(t) since 1 is obviously a super-solution of 1. A direct calculation gives that

    L¯u(x,t)=˙ξ+(x,t)ϕ2(ξ+(x,t))+˙v(t)RJ(xy)ϕ2(ξ+(y,t))dy+ϕ2(ξ+(x,t))f(x,ϕ2(ξ+(x,t))+v(t))
    =˙V(t)ϕ2(ξ+(x,t))+˙v(t)+f2(ϕ2(ξ+(x,t)))f(x,ϕ2(ξ+(x,t))+v(t)){˙V(t)ϕ2(ξ+(x,t))+˙v(t)+f2(ϕ2(ξ+(x,t)))f2(ϕ2(ξ+(x,t))+v(t))forx<L,˙V(t)ϕ2(ξ+(x,t))+˙v(t)+f2(ϕ2(ξ+(x,t)))f1(ϕ2(ξ+(x,t))+v(t))forxL.

    We continue to prove L¯u(x,t)0 in three cases.

    Case 1. If ξ+>X and xL, then

    L¯u(x,t)˙V(t)ϕ2(ξ+(x,t))+˙v(t)+f1(ϕ2(ξ+(x,t)))f1(ϕ2(ξ+(x,t))+v(t))C0|1ϕ2(ξ+(x,t))|˙V(t)ϕ2(ξ+(x,t))+˙v(t)(f1(1)+ω)˙v(t)C2C0eλ1(x+c2(t+T)+β++V(t))˙V(t)ϕ2(ξ+(x,t))+˙v(t)(f1(1)+ω)˙v(t)C2C0eλ1c2teλ1(Lc2T).

    Particularly, choose C3 such that

    ˙v(t)=ωv(t)+C2C0eλ1c2teλ1(Lc2T),

    which means

    C3=C2C0eλ1(Lc2T)λ1c2ω.

    As a consequence,

    L¯u(x,t)(f1(1)+2ω)v(t)0.

    In addition, if ξ+>X and x<L, then

    L¯u(x,t)ωv(t)(wλ1c2)C3eλ1c2t(f2(1)+ω)v(t)(f2(1)+2ω)v(t)0.

    Case 2. If ξ+<X, then x<L. Similar to Case 1 with x<L, we can obtain L¯u(x,t)0.

    Case 3. For ξ+[X,X], we know x<L. Therefore,

    L¯u(x,t)˙V(t)ϕ2(ξ+(x,t))+˙v(t)f2v(t)˙V(t)τ0˙V(t)ωv(t)f2v(t)(KVτ0ωf2)v(t),

    where τ0 is defined as in Lemma 4.2. It follows from

    KVω+f2τ0

    that L¯u(x,t)0. This finishes the proof by letting β=β++ϵ1+C3 and C=C3+ϵ1.

    Now we are ready to show that the entire solution converges to a shift of ϕ2(x+c2t) using a "squeezing" technique [7]. We should declare that assumption (A) (see (B2) in [7]) holds.

    (A) There exists a positive nonincreasing function ζ(n) defined on [1,+) such that for any u1(x,t),u2(x,t) satisfying 1u1(x,t),u2(x,t)2, Lu1(x,t)0, Lu2(x,t)0 and u1(x,0)u2(x,0), there holds

    minx[n,n]{u1(x,1)u2(x,1)}ζ(n)10[u1(y,0)u2(y,0)]dyform1.

    In fact, it follows from the comparison principle that u(x,t)v(x,t) for xR and t0. Furthermore, denote ˘w(x,t)=eKt[u(x,t)v(x,t)]0 with K1+maxs[1,2]|fs(x,s)|. Then

    ˘wt(x,t)K˘w(x,t)+J˘w(x,t)˘w(x,t)+f(x,u(x,t))f(x,v(x,t))J˘w(x,t)+[K(1+maxs[1,2]|fs(x,s)|)]˘w(x,t)J˘w(x,t),

    which implies that ˘w(x,t)˘w(x,0) for t0. Therefore, ˘wt(x,t)J˘w(x,0) and ˘w(x,t)tJ˘w(x,0). Repeating the same progress on [t,2t],...,[(N1),Nt], we have ˘w(x,Nt)tNJ˘w(x,0) for N1. In addition, since J(x)0 and RJ(y)dy=1, there exists N0=N0(M)1 such that by letting t=TN0 we have

    ˘w(x,T)(TN0)NJ˘w(x,0)(TN0)Nc(M)10˘w(y,0)dy,

    where

    c(M)=minx[M1,M+1]J(x)>0.

    Thus, there exists a positive function η(x,t)C([0,),[0,)) such that

    u(x,t)v(x,t)η(|x|,t)10[u(y,0)v(y,0)]dyforxR,t>0,

    which implies the assumption (A) holds.

    Now, we start to prove the following lemma, which plays an important role in the proof of Theorem 4.1.

    Lemma 4.5. Suppose that assumptions (F) and (J) hold. Then there exists a small ϵ0 such that if for some τ0,ξR,δ(0,σ2], and h>0, there holds

    ϕ2(x+c2τ+ξ)δu(x,t)ϕ2(x+c2τ+ξ+h)+δforxR, (9)

    then for every t>τ+1, there exist ˜ξ(t),˜δ(t) and ˜h(t) satisfying

    ˜ξ(t)[ξδγCm,ξ+h+δγCm],˜δ(t)eω(tτ1)[δ+ϵ0min{h,1}],˜h(t)τhγϵ0Cmmin{h,1}+2δγCm

    such that 9 holds with (τ,ξ,δ,h) replaced by (t,˜ξ(t),˜δ(t),˜h(t)). Here, Cm:=max{δ,C} and γ:=max{Cv,KVω}, where the parameters are defined in Lemmas 4.2 and 4.4.

    Proof of Lemma 4.5. In view of Lemmas 4.2 and 4.4, it is easy to see that

    ϕ2(x+c2tβ)δeωtu(x,t)ϕ2(x+c2t+β)+Ceωt.

    Furthermore, as in the proofs of Lemmas 4.2 and 4.4,

    ϕ2(x+c2tβ0+CveωtCv)δeωtu(x,t)

    and

    ϕ2(x+c2t+β++KVωKVωeωt)+Ceωtu(x,t).

    Denote

    κ=max{β0γ},h=β++KVωκ.

    In addition, by the definition of Cm, γ and letting ˇu(x,t)=u(xκ,t), we have

    ϕ2(x+c2t+γeωt)Cmeωtˇu(x,t)ϕ2(x+c2t+h+γeωt)+Cmeωt.

    Let ¯h=min{h,1} and ϑ0=12min[0,2]ϕ2(x). Then

    10[ϕ2(y+¯h)ϕ2(y)]dy2ϑ0¯h.

    Therefore, at least one of the following two inequalities is true

    (i)10[ˇu(y,0)ϕ2(y)]dyϑ0¯h,(ii)10[ϕ2(y+¯h)ˇu(y,0)]dyϑ0¯h.

    Next, we consider the case (i) since the case (ii) is similar. According to the assumption (A), for ζ=ζ(M+c2+2) with M1 such that γCmϕ2(x)<1 for |x|M, and for every x[Mc2γ,M+c2+γ], there holds

    ˇu(x,1)[ϕ2(x+γeω)Cmeω]ζ10[ˇu(y,0)(ϕ2(y)Cm)]dyζϑ0¯h.

    Now define

    ϵ0=min{σ2,γ2Cm,minx[Mc2γ,M+c2+γ]ζϑ0γ2Cmϕ2(x)}.

    Accordingly, there exists ˜θ(0,1) such that

    ϕ2(x+γeω+2ϵ0¯hγCm)ϕ2(x+γeω)=ϕ2(x+γeω+2ϵ0¯h˜θγCm)2ϵ0¯hγCmζϑ0¯h

    for all x[Mc2(γ1),M+c2+(γ1)]. Hence,

    ˇu(x,1)ϕ2(x+γeω+2ϵ0¯hγCm)Cmeω

    for x[Mc2(γ1),M+c2+(γ1)].

    For |x|M+c2+(γ1), by the choice of M, we know that

    ϕ2(x+γeω)ϕ2(x+γeω+2ϵ0¯hγCm)ϵ0¯h.

    Then, it follows that

    u(x,1)ϕ2(x+κ+γeω+2ϵ0¯hγCm)ϵ0¯hCmeωforxR.

    Note that p:=δeω+ϵ0¯hσ, repeat the operation above with 1 replaced by 1+t, which implies that

    u(x,1+t)ϕ2(x+c2t+κ+γeω+2ϵ0¯hγCm+γeωt)peωtϕ2(x+c2(t+1)+κ+ϵ0¯hγCmδγCm)(δ+ϵ0¯h)eωt.

    Thus we finish the proof by setting

    t=1+t,˜ξ(t)=γϵ0¯hCm,˜δ=(δ+ϵ0¯h)eω(t1)˜h=[h+δγCmeωt]˜ξ(t)=hϵ0γCm¯h+δγCm[2eωt].

    Now we shall prove Theorem 4.1.

    Proof of Theorem 4.1. We shall divide the proof in three steps.

    Step 1. Following from Lemmas 4.2 and 4.4, there exist T,M>0 such that

    ϕ2(x+c2TM)Cmu(x,T)ϕ2(x+c2T+M)+CmforxR. (10)

    Here, Cm is defined as in the proof of Lemma 4.4. Define

    ϵ=min{σ2,ϵ04}andk0=ϵ0γCm2ϵγCmγϵ02Cm>0.

    Meanwhile, fix t2 such that

    eω(t1)[1+ϵ0ϵ]1k.

    Then, replace Cm with ϵ and denote M,T by h02,T0. Assume that h01, otherwise, we directly go to Step 2.

    By 10, applying Lemma 4.5 with τ=T0, ξ=h02, h=h0, δ=ϵ, 9 holds with τ=T0+t, some ξ[h02γϵCm,h02+γϵCm], δ=ϵ and h=h0k, by the definition of k,t, it follows that

    ˆδ(T0+t)[ϵ+ϵ0]eωtϵandˆh(T0+t)h0γCmϵ0+2ϵγCmh0k.

    Repeat the same process, it yields that 9 holds for τ=T0+Nt, δ=ϵ and h=h0Nk with N such that h0(N1)t1. Thus, there exists T1>T0 such that 9 holds for τ=T1, δ=ϵ and h=1 and some ξR.

    Step 2. In this step, we use a mathematical induction to show that for every nonnegative integer k, 9 holds for ξ=ξkR and

    τ=Tk:=T1+kt,δ=δk:=(1k)kϵ,h=hk:=(1k)k.

    It is obvious that the assertion holds for k=0 by Step 2. Then suppose that the assertion is true for k=l0. We show that it is true for k=l+1. In fact, as can be seen in Lemma 4.5 with τ=Tl and t=Tl+1, one can obtain that 9 holds with (τ,ξ,δ,h) replaced by (ˆτ,ˆξ,ˆδ,ˆh) satisfying

    ˆξ[ξlδlγCm,ξl+δlγCm],ˆδ(δl+ϵ0hl)=[1k]lϵ(1+ϵ0ϵ)eω(t1)(1k)l+1ϵ,ˆhhlhlϵ0γCm+2δlγCm=[1k]l[1ϵ0γCm+2ϵγCm]=[1k]l+1

    by the definition of ϵ,k and t. This means that 9 holds for τ=Tl+1, some ξ=ξl+1[ξlδlγCm,ξl+δlγCm], δ=[1k]l+1ϵ and h=[1k]l+1. Now we finish the mathematical induction.

    Step 3. So far, we have known that 9 holds for (τ,ξ,δ,h)=(Tk,ξk,δk,hk) for all k=0,1.... Furthermore, 9 holds with τ[Tk,), δ=δk, h=hk+2δkγCm and ξ=ξkδkγCm,k=0,1....

    Now we define

    δ(t)=δk,ξ(t)=ξkδkγCm,h(t)=hk+2δkγCm

    for t[Tk,Tk+1],k=0,1.... Then,

    ϕ2(x+c2t+ξ(t))δ(t)u(x,t)ϕ2(x+c2t+ξ(t+h(t))+δ(t)fortT1,xR.

    It follows from the definition of δ(t) and h(t) that

    δ(t)=δk=[1k]kϵϵexp[(tT1t1)ln(1k)]fortT1,h(t)=hk+2δkγCm[1+2ϵγCm]exp[(tT1t1)ln(1k)]fortT1.

    Moreover, since for any tτT1,

    ξ(t)[ξ(τ)δ(τ)γCm,ξ(τ)+h(τ)+δ(τ)γCm],

    there holds

    |ξ(t)ξ(τ)|h(τ)+2δ(τ)γCm,

    which implies that ξ():=limt+ξ(t) exists and

    |ξ()ξ(τ)|h(τ)+2δ(τ)γCm[1+4ϵ2δ(τ)γCm]e[(tT1t1)ln(1k)]fortT1.

    Therefore, we have that

    |u(x,t)ϕ2(x+c2t+ξ())|0ast+.

    Furthermore, the convergence is exponential. Then we complete the proof.

    We investigate the Lyapunov stability of the entire solution obtained in Theorem 1.1 in this section. That is, the aim here is to prove Theorem 1.2. The following lemma plays an important role in proving Theorem 1.2.

    Lemma 5.1. Let u(x,t) be the unique entire solution in Theorem 1.1. Then for any φ(0,12], there exist constants Tφ=Tφ(φ)>1 and Kφ=Kφ(φ)>0 such that

    ut(x,t)Kφfor anytTφandxΩφ(t),

    where

    Ωφ(t)={xR:φu(x,t)1φ}.

    Proof of Lemma 5.1. It is easy to choose Tφ and Mφ such that

    Ωφ(t){xR:|x+c2t|Mφ}{xR:x1}.

    Now suppose there exist sequences tk[Tφ,+) and xkΩφ(t) such that

    ut(tk,xk)0ask+.

    Here only two cases happen, tk+ or tkt for some t[Tφ,+) as k+.

    For the former case, denote

    uk(x,t)=u(x+xk,t+tk).

    By Proposition 2, {uk(x,t)}k=1 is equicontinuous in xR and tR. Furthermore, there exists a subsequence still denoted by {uk(x,t)}k=1 such that

    ukuask+

    for some function u satisfying 1 with u(0,0)t=0 and (u(x,t))t0. By the comparison principle, we further have

    u(x,t)t0fort0.

    However, this is impossible because by Theorem 4.1

    u(x,t)=ϕ2(x+c2t+β+a)for somea[Mη,Mη].

    For the second case, xk remains bounded by the definition of Ωφ(t). Therefore, we assume that xkx as k+ and let

    uk(x,t):=u(x+xk,t).

    Then, each uk(x,t) is defined for all (x,t){xRx1}×[Tφ,+) by the definition of Ωφ(t). Similarly, there exists a subsequence, again denoted by {uk}k=1, such that

    ukuask+

    for some function u satisfying 1 with ut(0,t)=0,(u(x,t))t0. Then by the comparison principle, we obtain ut(x,t)0 for tt, but this is impossible since by Theorem 4.1,

    u(x,t)ϕ2(x+β+x+c2t)0ast+.

    This ends the proof.

    Proof of Theorem 1.2. We first define a pair of sub- and super-solutions as follows.

    U±(x,t)=u(x+a,t+t0±˜δϖ(1eωt))±ϖeωt,

    where ω is defined as in Section 4 and t0Tφ and ˜δ,ϖ>0 are constants. Besides, we claim that U± are super- and sub-solutions of 1. Then, by the comparison principle, we have

    U(x,t)v(x,t)U+(x,t)forxR,t0. (11)

    In view of that for all |τ|ϖ with ϖ=:ϖ(ϵ)=ϵ2˜M and ˜M=2+maxu[0,1](fL(x,u))u, we have

    supxR,tR|u(x,t)u(x,t+τ)|supxR,tR|ut(x,t)|ϖϵ2. (12)

    It then follows from 11 and 12 that

    |v(x,t,v0)u(x+a,t+t0)|<ϵ

    for xR and t0.

    Now we prove the claim. We show U+(x,t) is a super-solution of 1 for t0. Then it can be similarly shown that U(x,t) is a sub-solution of 1 for t0. Since it is easy to see that U+(x,0)v0(x)forxRandδ<ϖ2, we only need to show that

    LU+(x,t)=U+t(x,t)(JU+(x,t)U+(x,t))f(x,U+)=˜δϖωeωtϖωeωt+f(x,u)f(x,U+)0.

    We go further to show LU+(x,t)0 in two cases.

    Case 1. For any u=u(x+a,t+γ+˜δϖ(1eωt))[0,σ][1σ,1] with σ defined as in Section 4,

    LU+(x,t)min{|fu(x,0)|2,|(fu(x,1)|2}ϖeωtωϖeωtϖeωt(2ωω)0.

    Case 2. For u[σ,1σ], it follows from Lemma 5.1 that there exists Kσ such that utKσ. Therefore,

    LU+(x,t)˜δϖωKσeωtϖωeωtfu(x,u)ϖeωtϖeωt(˜δωKσωfu(x,u))).

    Let ˜δ be sufficiently large such that

    ˜δωKσω+fu(x,u).

    Thus, we have LU+(x,t)0. Now, we obtain the claim. The proof of Theorem 1.2 is completed.

    Remark 1. In this paper, we have considered the existence, uniqueness, asymptotic behavior and Lyapunov stability of entire solutions of the nonlocal dispersal equation 1 under the assumption J(x)=J(x), xR. It is well-known from [8,22,28,29] that the asymmetry of J has a great influence on the profile of the traveling waves and the sign of the wave speeds, which further makes the properties of the entire solution more diverse. Naturally, an interesting problem is to consider entire solutions of 1 under asymmetric conditions.

    Remark 2. The method used here can be also applied to consider the bistable lattice differential equations with fi satisfying (F)

    ˙ui(t)=ui+1(t)+ui1(t)2ui(t)+fi(ui(t)),iZ,tR. (13)

    The existence, uniqueness, asymptotic behavior and the Lyapunov stability of entire solutions to 13 can be similarly obtained.

    While for bistable random diffusion equations with f(x,u) satisfying (F)

    ut=uxx+f(x,u),xR,tR,

    the construction of the entire solution behaving as the traveling wave pertaining to f1 coming from infinity is similar. In fact, since regularity of the entire solution is known, the Lyapunov function method can be used to show the approach of the entire solution to the traveling wave pertaining to f2, one can refer to [10,11].



    [1] P. W. Bates, On some nonlocal evolution equations arising in materials science, Nonlinear Dynamics and Evolution Equations, in: Fields Inst. Commun., Amer. Math. Soc., Providence, RI, 48 (2006), 13–52.
    [2] Traveling waves in a convolution model for phase transitions. Arch. Rational Mech. Anal. (1997) 138: 105-136.
    [3] Generalized transition waves and their properties. Comm. Pure Appl. Math. (2012) 65: 592-648.
    [4] Bistable traveling waves around an obstacle. Comm. Pure Appl. Math. (2009) 62: 729-788.
    [5] A non-local bistable reaction-diffusion equaiton with a gap. Discrete Contin. Dyn. Syst. (2017) 37: 685-723.
    [6] Uniqueness of travelling waves for nonlocal monostable equations. Proc. Amer. Math. Soc. (2004) 132: 2433-2439.
    [7] Existence, uniqueness and asymptotic stablility of traveling waves in nonlocal evolution equations. Adv. Differential Equations (1997) 2: 125-160.
    [8] Nonlocal anisotropic dispersal with monostable nonlinearity. J. Differential Equations (2008) 244: 3080-3118.
    [9] On a nonlocal reaction diffusion equation arising in population dynamics. Proc. Roy. Soc. Edinburgh Sect. (2007) 137: 727-755.
    [10] A heteroclinic orbit connecting traveling waves pertaining to different nonlinearities. J. Differential Equations (2018) 265: 804-829.
    [11] A heteroclinic orbit connecting traveling waves pertaining to different nonlinearities in a channel with decreasing cross section. Nonlinear Anal. (2018) 172: 99-114.
    [12] P. Fife, Some nonclassical trends in parabolic and parabolic–like evolutions, in: Trends in Nonlinear Analysis, Springer, Berlin (2003), 153–191.
    [13] Entire solutions of the KPP equation. Comm. Pure Appl. Math. (1999) 52: 1255-1276.
    [14] Travelling fronts and entire solutions of the Fisher-KPP equation in RN. Arch. Ration. Mech. Anal. (2001) 157: 91-163.
    [15] Entire solutions in the Fisher-KPP equation with nonlocal dispersal. Nonlinear Anal. Real World Appl. (2010) 11: 2302-2313.
    [16] Entire solutions of nonlocal dispersal equations with monostable nonlinearity in space periodic habitats. J. Differential Equations (2016) 261: 2472-2501.
    [17] Transition fronts for inhomogeneous Fisher-KPP reactions and non-local diffusion. Trans. Amer. Math. Soc. (2016) 368: 8615-8631.
    [18] Existence and stability of traveling wavefronts in a nonlocal diffusion equation with delay. Nonlinear Anal. (2010) 72: 3150-3158.
    [19] Traveling waves in diffusive random media. J. Dynam. Differential Equations (2004) 16: 1011-1060.
    [20] Stability of transition waves and positive entire solutions of Fisher-KPP equations with time and space dependence. Nonlinearity (2017) 30: 3466-3491.
    [21] Spreading speeds for monostable equations with nonlocal dispersal in space periodic habitats. J. Differential Equations (2010) 249: 747-795.
    [22] Traveling waves for a nonlocal anisotropic dispersal equation with monostable nonnlinearity. Nonlinear Anal. (2011) 74: 814-826.
    [23] Entire solutions for nonlocal dispersal equations with bistable nonlineartiy. J. Differential Equations (2011) 251: 551-581.
    [24] Entire solutions of a diffusion and competitive Lotka-Volterra type system with nonlocal delayes. Nonlinearity (2010) 23: 1609-1630.
    [25] J.-B. Wang and C. Wu, Forced waves and gap formations for a Lotka-Volterra competition model with nonlocal dispersal and shifting habitats, Nonlinear Anal. Real World Appl., 58 (2021), 103208. doi: 10.1016/j.nonrwa.2020.103208
    [26] Entire solutions for nonlocal dispersal equations with spatio-temporal delay: Monostable case. J. Differential Equations (2015) 258: 2435-2470.
    [27] Existence of traveling wave solutions for a nonlocal bistable equation: An abstract approach. Publ. Res. Inst. Math. Sci. (2009) 45: 955-979.
    [28] Entire solution in an ignition nonlocal dispersal equation: Asymmetric kernel. Sci. China Math. (2017) 60: 1791-1804.
    [29] Entire solutions in nonlocal bistable equations: Asymmetric case. Acta Math. Sin. (Engl. Ser.) (2019) 35: 1771-1794.
  • This article has been cited by:

    1. Xing-Xing Yang, Guo-Bao Zhang, Entire Solutions for an Inhomogeneous Bistable Discrete Diffusive Equation, 2023, 46, 0126-6705, 10.1007/s40840-022-01451-5
  • Reader Comments
  • © 2021 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(2666) PDF downloads(324) Cited by(1)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog