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

Regularization effect of the mixed-type damping in a higher-dimensional logarithmic Keller-Segel system related to crime modeling


  • We study a logarithmic Keller-Segel system proposed by Rodríguez for crime modeling as follows:

    {ut=Δuχ(ulnv)κuv+h1,vt=Δvv+u+h2,

    in a bounded and smooth spatial domain ΩRn with n3, with the parameters χ>0 and κ>0, and with the nonnegative functions h1 and h2. For the case that κ=0, h10 and h20, recent results showed that the corresponding initial-boundary value problem admits a global generalized solution provided that χ<χ0 with some χ0>0.

    In the present work, our first result shows that for the case of κ>0 such problem possesses global generalized solutions provided that χ<χ1 with some χ1>χ0, which seems to confirm that the mixed-type damping κuv has a regularization effect on solutions. Besides the existence result for generalized solutions, a statement on the large-time behavior of such solutions is derived as well.

    Citation: Bin Li, Zhi Wang, Li Xie. Regularization effect of the mixed-type damping in a higher-dimensional logarithmic Keller-Segel system related to crime modeling[J]. Mathematical Biosciences and Engineering, 2023, 20(3): 4532-4559. doi: 10.3934/mbe.2023210

    Related Papers:

    [1] A. Q. Khan, Ibraheem M. Alsulami . Discrete Leslie's model with bifurcations and control. AIMS Mathematics, 2023, 8(10): 22483-22506. doi: 10.3934/math.20231146
    [2] Xiongxiong Du, Xiaoling Han, Ceyu Lei . Dynamics of a nonlinear discrete predator-prey system with fear effect. AIMS Mathematics, 2023, 8(10): 23953-23973. doi: 10.3934/math.20231221
    [3] Fatao Wang, Ruizhi Yang, Yining Xie, Jing Zhao . Hopf bifurcation in a delayed reaction diffusion predator-prey model with weak Allee effect on prey and fear effect on predator. AIMS Mathematics, 2023, 8(8): 17719-17743. doi: 10.3934/math.2023905
    [4] Abdul Qadeer Khan, Ayesha Yaqoob, Ateq Alsaadi . Neimark-Sacker bifurcation, chaos, and local stability of a discrete Hepatitis C virus model. AIMS Mathematics, 2024, 9(11): 31985-32013. doi: 10.3934/math.20241537
    [5] Xiaoming Su, Jiahui Wang, Adiya Bao . Stability analysis and chaos control in a discrete predator-prey system with Allee effect, fear effect, and refuge. AIMS Mathematics, 2024, 9(5): 13462-13491. doi: 10.3934/math.2024656
    [6] Ahmad Suleman, Rizwan Ahmed, Fehaid Salem Alshammari, Nehad Ali Shah . Dynamic complexity of a slow-fast predator-prey model with herd behavior. AIMS Mathematics, 2023, 8(10): 24446-24472. doi: 10.3934/math.20231247
    [7] Weili Kong, Yuanfu Shao . Bifurcations of a Leslie-Gower predator-prey model with fear, strong Allee effect and hunting cooperation. AIMS Mathematics, 2024, 9(11): 31607-31635. doi: 10.3934/math.20241520
    [8] Ali Al Khabyah, Rizwan Ahmed, Muhammad Saeed Akram, Shehraz Akhtar . Stability, bifurcation, and chaos control in a discrete predator-prey model with strong Allee effect. AIMS Mathematics, 2023, 8(4): 8060-8081. doi: 10.3934/math.2023408
    [9] Abdul Qadeer Khan, Ayesha Yaqoob, Ateq Alsaadi . Discrete Hepatitis C virus model with local dynamics, chaos and bifurcations. AIMS Mathematics, 2024, 9(10): 28643-28670. doi: 10.3934/math.20241390
    [10] Fethi Souna, Salih Djilali, Sultan Alyobi, Anwar Zeb, Nadia Gul, Suliman Alsaeed, Kottakkaran Sooppy Nisar . Spatiotemporal dynamics of a diffusive predator-prey system incorporating social behavior. AIMS Mathematics, 2023, 8(7): 15723-15748. doi: 10.3934/math.2023803
  • We study a logarithmic Keller-Segel system proposed by Rodríguez for crime modeling as follows:

    {ut=Δuχ(ulnv)κuv+h1,vt=Δvv+u+h2,

    in a bounded and smooth spatial domain ΩRn with n3, with the parameters χ>0 and κ>0, and with the nonnegative functions h1 and h2. For the case that κ=0, h10 and h20, recent results showed that the corresponding initial-boundary value problem admits a global generalized solution provided that χ<χ0 with some χ0>0.

    In the present work, our first result shows that for the case of κ>0 such problem possesses global generalized solutions provided that χ<χ1 with some χ1>χ0, which seems to confirm that the mixed-type damping κuv has a regularization effect on solutions. Besides the existence result for generalized solutions, a statement on the large-time behavior of such solutions is derived as well.



    The well-known Leslie-Gower model with a Holling Ⅱ functional response is expressed as follows:

    {dxdt=rx(1xK)bxya+x,dydt=sy(1yx), (1.1)

    where the biological interpretation of each parameter is provided in Table 1. It is noteworthy that the Leslie-Gower model incorporates a response function indicating the predator's carrying capacity, which is proportional to the prey population, a feature absent in the Lotka-Volterra model. This distinction has been extensively investigated in many references, such as [1,2,3,4,5].

    Table 1.  Biological meaning of parameters in model (1.1).
    Parameter Interpretation
    x prey population density
    y predator population density
    r intrinsic growth rate of the prey
    K maximum prey carrying capacity of the environment
    b maximum predation rate when prey is abundant
    a prey population at which predator attack capability saturates
    s intrinsic growth rate of the predator

     | Show Table
    DownLoad: CSV

    Previous studies predominantly focused on direct predation effects. In 2016, Wang et al. [6] demonstrated through experiments that predator-induced fear (indirect effects) in prey leads to a reduction in prey birth rates. They introduced the fear factor F(k,y)=11+ky, where k represents the intensity of fear driving antipredator behavior. To further investigate the influence of fear on population dynamics, we enhance the system (1.1) by incorporating the fear factor F(k,y). The modified model is formulated as follows:

    {dxdt=rx1+ky(1xK)bxya+x,dydt=sy(1yx). (1.2)

    For additional significant findings on the impact of fear in predator-prey models, refer to [7,8,9].

    There has been increasing recognition that traditional integer-order differential equations may not sufficiently capture the complexity of biological systems [10,11,12,13,14]. In reality, the behaviors of most organisms in nature are influenced by their historical context. Fractional-order derivatives, which extend the concept of integer-order differentiation, offer a more flexible and accurate framework for modeling memory and hereditary properties in population dynamics. Therefore, the authors of this paper aim to apply the Caputo fractional derivative to (1.2), thereby extending it into a fractional model. The Caputo fractional derivative of a function u(t) of order α(0,1] is given in [15] as follows

    Dαtu(t)=1Γ(1α)t0u(τ)(tτ)αdτ,

    where Γ is the Gamma function. By replacing the integer-order derivative with the Caputo fractional derivative, the following model is obtained:

    {Dαtx=rx1+ky(1xK)bxya+x,Dαty=sy(1yx). (1.3)

    Specifically, when α=1, the fractional-order model (1.3) reduces to the integer-order model (1.2), demonstrating that the fractional-order model serves as a generalization of the integer-order model. Moreover, in fractional calculus, the rate of change at any given moment, expressed by the fractional-order derivative, depends on the population density over a specified time interval. This feature gives the fractional-order model (1.3) a distinct advantage in capturing memory effects within populations.

    When studying populations with nonoverlapping generations or small sizes, mathematical models are often expressed in discrete terms. Although some variables may change continuously in real-world scenarios, the recording of these changes typically happens at specific time intervals during data collection. Thus, employing discrete systems to examine the dynamic behavior of biological populations is highly practical and significant. In references [16,17], the authors employed the piecewise constant approximation method to discretize continuous fractional predator-prey models and explored their dynamical properties. Singh and Sharma [18] examine a discrete prey-predator model with Holling type Ⅱ functional response and prey refuge, identifying bifurcations and controlling chaos through state feedback, pole placement, and hybrid techniques. Berkal and Almatrafi [19] used the exponential piecewise constant argument to discretize continuous fractional activator-inhibitor system. Their analysis includes stability assessments, investigations of Neimark-Sacker and period-doubling bifurcations, and numerical simulations that validate the theoretical findings on the system's dynamics. For a more comprehensive exploration of discrete model studies, readers are advised to consult references [20,21,22,23,24,25] and the related literature cited therein. However, current literature lacks studies on fractional discrete-time predator-prey Leslie-Gower systems that incorporate the fear effect in prey population. Using the same method as in references [16,17] to discretize model (1.3), we can obtain the following discrete model:

    {xn+1=xn+hαΓ(α+1)[rxn1+kyn(1xnK)bxnyna+xn],yn+1=yn+hαΓ(α+1)[syn(1ynxn)]. (1.4)

    Here h>0 represents the time interval of production.

    Furthermore, the key contributions and findings of this study are summarized as follows:

    ● The Caputo fractional derivative of order (0,1] is utilized to incorporate the memory effect into the dynamical behavior of the proposed model.

    ● The existence and stability of fixed points are investigated.

    ● Conditions for the occurrence and direction of period-doubling bifurcation and Neimark–Sacker bifurcation at the positive fixed point are established.

    ● State feedback and hybrid control strategies are employed to manage bifurcations and chaotic behavior in the model.

    ● To validate the accuracy of our theoretical findings, numerical examples for the fractional-order discrete-time Leslie-Gower model with a fear factor are provided.

    The remainder of this paper is structured as follows: In Section 2, we investigate the existence and stability of the equilibrium points of model (1.4). Section 3 analytically demonstrates that model (1.4), under specific parametric conditions, undergoes period-doubling or Neimark-Sacker bifurcation. Section 4 explores the control of chaos toward an unstable equilibrium point using feedback control or hybrid control approaches. Section 5 presents a quantitative analysis of the dynamics of model (1.4) to validate our analytical findings. Finally, Section 6 offers brief conclusions.

    This section explores the existence of equilibrium points in model (1.4) and assesses their stability by evaluating the eigenvalues of the Jacobian matrix at these points. The following definition and lemma are introduced to assist in this stability analysis of equilibrium points.

    Definition 1. [26] Let λ1 and λ2 denote the two roots of the characteristic equation F(λ)=λ2+pλ+q=0 associated with the Jacobian matrix J(x,y). The equilibrium point (x,y) is termed

    (1) sink if |λ1|<1 and |λ2|<1, and the sink is locally asymptotically stable;

    (2) source if |λ1|>1 and |λ2|>1, and the source is locally unstable;

    (3) saddle if |λ1|>1 and |λ2|<1 (or |λ1|<1 and |λ2|>1);

    (4) non-hyperbolic if either |λ1|=1 or |λ2|=1.

    Lemma 1. [26] Let F(1)>0 in F(λ)=λ2Mλ+N, where λ1,λ2 are the two roots of F(λ)=0. Then, the following results hold true:

    (1) |λ1|<1 and |λ2|<1 if, and only if, F(1)>0 and N<1;

    (2) |λ1|<1 and |λ2|>1 (or |λ1|>1 and |λ2|<1) if, and only if, F(1)<0;

    (3) |λ1|>1 and |λ2|>1 if, and only if, F(1)>0 and N>1;

    (4) λ1=1 and |λ2|1 if, and only if, F(1)=0 and M0,2;

    (5) λ1 and λ2 are conjugate complex and |λ1,2|=1 if, and only if, M2<4N and N=1.

    It is evident that the equilibrium points of model (1.4) satisfy the following equations:

    {x=x+hαΓ(α+1)[rx1+ky(1xK)bxya+x],y=y+hαΓ(α+1)[sy(1yx)].

    This algebraic system is satisfied if x=K and y=0, indicating that the model (1.4) has a boundary equilibrium point E0(K,0) for all model parameters. To find the interior equilibrium point E1, we will solve the following system simultaneously:

    r1+ky(1xK)bya+x=0,s(1yx)=0.

    From the second equation, we find y=x. Substituting y=x into the first equation yields:

    (rK+bk)x2+(arK+br)xar=0,

    Given that all model parameters are positive, we obtain the following results through direct calculations.

    Theorem 1. The model (1.4) always has a boundary equilibrium point E0(K,0) and a positive equilibrium point E1=(x,y), where

    x=y=(arK+br)+(arK+br)2+4ar(rK+bk)2(rK+bk).

    The Jacobian matrix of model (1.4) evaluated at any point (x,y) is as follows:

    J(x,y)=[1+hαΓ(α+1)[r1+ky(12xK)aby(a+x)2]hαΓ(α+1)[rkx(1+ky)2(1xK)+bxa+x]hαΓ(α+1)sy2x21+hαΓ(α+1)(s2syx)]. (2.1)

    Theorem 2. The boundary equilibrium point E0(K,0) exhibits the following behaviors:

    (1) It is a source point if r>2Γ(α+1)hα;

    (2) It is a saddle point if 0<r<2Γ(α+1)hα;

    (3) It is a non-hyperbolic point if r=2Γ(α+1)hα.

    Proof. The Jacobian matrix (2.1) evaluated at the boundary equilibrium point E0(K,0) is given by:

    J(E0)=[1hαrΓ(α+1)hαΓ(α+1)bKa+K01+hαsΓ(α+1)].

    The eigenvalues of J(E0) are λ1=1hαrΓ(α+1) and λ2=1+hαsΓ(α+1). Clearly, |λ2|>1 and:

    |λ1|={<1if r>2Γ(α+1)hα,>1if 0<r<2Γ(α+1)hα,=1if r=2Γ(α+1)hα.

    By applying Definition 1, the proof is complete.

    Next, we analyze the local dynamics of model (1.4) at the positive equilibrium point E1(x,y). The Jacobian matrix (2.1) evaluated at E1(x,y) is expressed as:

    J(E1)=(1+Aa11Aa12As1As), (2.2)

    where

    A=hαΓ(α+1),a11=r1+kx(12xK)abx(a+x)2,a12=rkx(1+kx)2(1xK)bxa+x.

    The characteristic equation of J(E1) is given by:

    λ2Mλ+N=0, (2.3)

    where

    M=2+(a11s)A,N=1+(a11s)A(a11+a12)sA2.

    Let F(λ)=λ2Mλ+N, then:

    F(0)=N,F(1)=1+M+N,F(1)=1M+N.

    As a result, we establish the following theorem.

    Theorem 3. Let E1 be the unique positive equilibrium point of model (1.4). Then:

    (1) E1 is a sink point if |M|<1+N<2;

    (2) E1 is a source point if |M|<|1+N| and |N|>1;

    (3) E1 is a saddle point if M2>4N and |M|>|1+N|;

    (4) E1 is non-hyperbolic if |M|=|1+N| or N=1 and |M|<2.

    Proof. (1) According to Lemma 1, E1 is a sink point if and only if F(1)>0, F(1)>0, and N<1. This condition is satisfied when |M|<1+N<2. Similarly, the proofs for Theorem 3 (2)–(4) follow straightforwardly from the definitions and properties of M and N.

    In this section, we will analyze the period-doubling and Neimark-Sacker bifurcation behaviors of model (1.4) at the positive equilibrium point E1(x,y). The reason for not analyzing the boundary equilibrium point is that, at E0(K,0), the predators have become extinct, leaving only the prey.

    Assume M2>4N and

    s=s1:=4+2Aa112A+A2(a11+a12),

    and we can ascertain that the characteristic equation (2.3) satisfies

    F(1)=1+M+N=1+[2+(a11s1)A]+[1+(a11s1)A(a11+a12)s1A2]=4+2a11A[2A+(a11+a12)A2]s1=0.

    Therefore, the eigenvalues of the characteristic equation (2.3) are

    λ1=1,λ2=1+M=3+(a11s1)A.

    Let

    s2=a11+2A,s3=a11+4A.

    The condition |λ2|1 implies that s1s2,s3.

    Based on the above analysis, we can conclude that when the parameters vary within the set

    P.D={(h,α,r,k,K,b,a,s):M2>4N,s=s1andss2,s3}, (3.1)

    model (1.4) will undergo period-doubling bifurcation at E1(x,y).

    Select parameters (h,α,r,k,K,b,a,s)P.D and consider s as a small perturbation of s, i.e., s=ss1, where |s|1. After this perturbation, model (1.4) can be represented as follows:

    {xn+1=xn+hαΓ(α+1)[rxn1+kyn(1xnK)bxnyna+xn],yn+1=yn+hαΓ(α+1)[(s1+s)yn(1ynxn)]. (3.2)

    Let un=xnx and vn=yny, transforming the equilibrium point E1(x,y) into the origin O(0,0). The model (3.2) can thus be expressed as:

    {un+1=un+hαΓ(α+1)[r(un+x)1+k(vn+y)(1un+xK)b(un+x)(vn+y)a+(un+x)],vn+1=vn+hαΓ(α+1)[(s1+s)(vn+y)(1vn+yun+x)]. (3.3)

    The Taylor expansion of model (3.3) around (un,vn)=(0,0) yields the following form:

    [un+1vn+1]=[1+Aa11Aa12As11As1][unvn]+[f(un,vn,s)g(un,vn,s)], (3.4)

    where

    f(un,vn,s)=c13u2n+c14unvn+c15v2n+c16u3n+c17u2nvn+c18unv2n+c19v3n+O((|un|+|vn|+|s|)4),g(un,vn,s)=c23u2n+c24unvn+c25v2n+c26u3n+c27u2nvn+c28unv2n+c29v3n+d1uns+d2vns+d3u2ns+d4unvns+d5v2ns+O((|un|+|vn|+|s|)4),

    and

    c13=ArK(1+kx)+abAx(a+x)3,c14=Ark(1+kx)2(12xK)abA(a+x)2,c15=Ark2x(1+kx)3(1xK),c16=Aabx(a+x)4,c17=ArkK(1+kx)2+abA(a+x)3,c18=Ark2(1+kx)3(12xK),c19=Ark3x(1+kx)4(1xK),c23=c25=As1x,c24=2As1x,c26=c28=As1x2,c27=2As1x2,c29=0,d1=A,d2=A,d3=d5=Ax,d4=2Ax.

    Define

    T=[Aa12Aa122Aa11λ21Aa11],

    then

    T1=1Aa12(1+λ2)[λ21Aa11Aa122+Aa11Aa12].

    Using the transformation:

    [unvn]=T[˜un˜vn],

    model (3.4) is transformed into the following form:

    [˜un+1˜vn+1]=[100λ2][˜un˜vn]+[˜f(˜un,˜vn,s)˜g(˜un,˜vn,s)], (3.5)

    where

    ˜f(˜un,˜vn,s)=(λ21Aa11)c13Aa12c23Aa12(1+λ2)u2n+(λ21Aa11)c14Aa12c24Aa12(1+λ2)unvn+(λ21Aa11)c15Aa12c25Aa12(1+λ2)v2n+(λ21Aa11)c16Aa12c26Aa12(1+λ2)u3n+(λ21Aa11)c17Aa12c27Aa12(1+λ2)u2nvn+(λ21Aa11)c18Aa12c28Aa12(1+λ2)unv2n+(λ21Aa11)c19Aa12c29Aa12(1+λ2)v3nd11+λ2unsd21+λ2vnsd31+λ2u2nsd41+λ2unvnsd51+λ2v2ns+O((|un|+|vn|+|s|)4),
    ˜g(˜un,˜vn,s)=(2+Aa11)c13+Aa12c23Aa12(1+λ2)u2n+(2+Aa11)c14+Aa12c24Aa12(1+λ2)unvn+(2+Aa11)c15+Aa12c25Aa12(1+λ2)v2n+(2+Aa11)c16+Aa12c26Aa12(1+λ2)u3n+(2+Aa11)c17+Aa12c27Aa12(1+λ2)u2nvn+(2+Aa11)c18+Aa12c28Aa12(1+λ2)unv2n+(2+Aa11)c19+Aa12c29Aa12(1+λ2)v3n+d11+λ2uns+d21+λ2vns+d31+λ2u2ns+d41+λ2unvns+d51+λ2v2ns+O((|un|+|vn|+|s|)4),

    and

    un=Aa12˜un+Aa12˜vn,vn=(2+Aa11)˜un+(λ21Aa11)˜vn.

    Next, we apply the center manifold theorem [26] to analyze the dynamics around the equilibrium point (˜un,˜vn)=(0,0) at s=0. According to the theorem, the model (3.5) has a center manifold, which can be represented as:

    Wc(0,0,0)={(˜un,˜vn,s)R3+:˜vn=w(˜un,s),w(0,0)=0,Dw(0,0)=0}.

    Assume that

    w(˜un,s)=η1˜u2n+η2˜uns+η3(s)2+O((|˜un|+|s|)3).

    Then, the center manifold must satisfy

    w(˜un+˜f(˜un,w(˜un,s),s),s)λ2w(˜un,s)˜g(˜un,w(˜un,s),s)=0.

    By comparing the coefficients, it can be obtained that

    η1=Aa12[(2+Aa11)c13+Aa12c23](2+Aa11)[(2+Aa11)c14+Aa12c24]1λ22+(2+Aa11)2[(2+Aa11)c15+Aa12c25]Aa12(1λ22),η2=(2+Aa11)d2Aa12d1(1+λ2)2,η3=0.

    Thus, the model (3.5), when restricted to the center manifold Wc(0,0,0), is given by:

    ˜W:˜un+1=˜un+m1˜u2n+m2˜uns+m3˜u3n+m4˜u2ns+m5˜un(s)2+O((|˜un|+|s|)4) (3.6)

    where

    m1=Aa12[(λ21Aa11)c13Aa12c23]1+λ2(2+Aa11)[(λ21Aa11)c14Aa12c24]1+λ2+(2+Aa11)2[(λ21Aa11)c15Aa12c25]Aa12(1+λ2)m2=d1Aa121+λ2+d2(2+Aa11)1+λ2m3=2Aa12η1[(λ21Aa11)c13Aa12c23]1+λ2+Aa12η1(λ232Aa11)[(λ21Aa11)c14Aa12c24]1+λ22η1(2+Aa11)(λ21Aa11)[(λ21Aa11)c15Aa12c25]Aa12(1+λ2)+(Aa12)2[(λ21Aa11)c16Aa12c26]1+λ2Aa12(2+Aa11)[(λ21Aa11)c17Aa12c27]1+λ2+(2+Aa11)2[(λ21Aa11)c18Aa12c28]1+λ2(2+Aa11)3[(λ21Aa11)c19Aa12c29]Aa12(1+λ2)m4=2Aa12η2[(λ21Aa11)c13Aa12c23]1+λ2+η2(λ232Aa11)[(λ21Aa11)c14Aa12c24]1+λ22η2(2+Aa11)(λ21Aa11)[(λ21Aa11)c15Aa12c25]Aa12(1+λ2)Aa12η1d11+λ2(λ21Aa11)η1d21+λ2(Aa12)2d31+λ2+Aa12(2+Aa11)d41+λ2(2+Aa11)2d51+λ2m5=Aa12η2d11+λ2(λ21Aa11)η2d21+λ2.

    In order for Eq (3.6) to undergo period-doubling bifurcation, it is necessary that the following two quantities possess nonzero values.

    β1=(2˜W˜uns+12˜Ws2˜W˜u2n)|(0,0)=m2β2=[163˜W˜u3n+(122˜W˜u2n)2]|(0,0)=m3+m21.

    Summarize the above analysis into the following theorem.

    Theorem 4. If β1β20, then model (1.4) undergoes a period-doubling bifurcation at the positive equilibrium point E1(x,y) when parameters vary in a small neighborhood of P.D. Additionally, when β2>0 (respectively, β2<0), model (1.4) bifurcates from the equilibrium point E1(x,y) to a stable (respectively, unstable) 2-periodic orbit.

    Assume M2<4N and

    s=s4:=a111+A(a11+a12),

    and we can determine that the eigenvalues of the characteristic equation (2.3) satisfy

    |λ1,2|2=|M2±4NM22i|2=1+(a11s4)A(a11+a12)s4A2=1+a11A[1+(a11+a12)A]As4=1.

    Namely, Eq (2.3) has two complex conjugate roots with unit modulus. It is clear from the above discussion that

    d|λ1,2|ds|s=s4=A2[1+(a11+a12)A]0.

    In addition, it is crucial that when s=s4, λθ1,2(s4)1(θ=1,2,3,4), which is equivalent to M(s4)2,1,0,2. Since M2<4N, we deduce M(s4)2,2. Additionally, we necessitate that M(s4)0,1, which leads to

    s4s5:=a11+3A,s4s6:=a11+2A.

    Based on the preceding analysis, we conclude that when the parameters vary within a small neighborhood of the set

    N.S={(h,α,r,k,K,b,a,s):M2<4N,s=s4,ss5,ss6}, (3.7)

    model (1.4) undergoes a Neimark-Sacker bifurcation at E1(x,y).

    Next, assuming (h,α,r,k,K,b,a,s)N.S. and |s|1 represents a small perturbation of s4, model (1.4) can be described as follows:

    {xn+1=xn+hαΓ(α+1)[rxn1+kyn(1xnK)bxnyna+xn],yn+1=yn+hαΓ(α+1)[(s4+s)yn(1ynxn)]. (3.8)

    Let un=xnx and vn=yny, then we have

    [un+1vn+1]=[1+Aa11Aa12As41As4][unvn]+[F(un,vn)G(un,vn)] (3.9)

    where

    F(un,vn)=c13u2n+c14unvn+c15v2n+c16u3n+c17u2nvn+c18unv2n+c19v3n+O((|un|+|vn|)4),G(un,vn)=c23u2n+c24unvn+c25v2n+c26u3n+c27u2nvn+c28unv2n+c29v3n+O((|un|+|vn|)4),

    and c13, c14, c15, c16, c17, c18, c19, c23, c24, c25, c26, c27, c28, c29 are given in (2.3) by substituting s1 for s4+s.

    In order to obtain the normal form of model (3.9) at s=0, we use the following transformation:

    [unvn]=[ω1+Aa11ρ0As4][ˆunˆvn],

    with

    ρ=M2,ω=4NM22.

    Using this transformation, model (3.9) will transform as follows:

    [ˆun+1ˆvn+1]=[ρωωρ][ˆunˆvn]+[ˆF(ˆun,ˆvn)ˆG(ˆun,ˆvn)], (3.10)

    where

    ˆF(ˆun,ˆvn)=As4c13+(ρ1Aa11)c23As4ωu2n+As4c14+(ρ1Aa11)c24As4ωunvn+As4c15+(ρ1Aa11)c25As4ωv2n+As4c16+(ρ1Aa11)c26As4ωu3n+As4c17+(ρ1Aa11)c27As4ωu2nvn+As4c18+(ρ1Aa11)c28As4ωunv2n+As4c19+(ρ1Aa11)c29As4ωv3n+O((|un|+|vn|)4),ˆG(ˆun,ˆvn)=1As4(c23u2n+c24unvn+c25v2n+c26u3n+c27u2nvn+c28unv2n+c29v3n)+O((|un|+|vn|)4),

    and

    un=ωˆun+(1+Aa11ρ)ˆvn,vn=As4ˆvn.

    Next, a nonzero real number is defined as follows:

    Ω=Re[(12λ1)λ221λ1ξ11ξ12]12|ξ11|2|ξ21|2+Re(λ2ξ22) (3.11)

    where

    ξ11=14[2ˆFˆu2n+2ˆFˆv2n+i(2ˆGˆu2n+2ˆGˆv2n)]|s=0ξ12=18[2ˆFˆu2n2ˆFˆv2n+22ˆGˆunˆvn+i(2ˆGˆu2n2ˆGˆv2n22ˆFˆunˆvn)]|s=0ξ21=18[2ˆFˆu2n2ˆFˆv2n22ˆGˆunˆvn+i(2ˆGˆu2n2ˆGˆv2n+22ˆFˆunˆvn)]|s=0ξ22=116[3ˆFˆu3n+3ˆFˆunˆv2n+3ˆGˆu2nˆvn+3ˆGˆv3n+i(3ˆGˆu3n+3ˆGˆunˆv2n3ˆFˆu2nˆvn3ˆFˆv3n)]|s=0.

    Through some complicated calculations, we get

    2ˆFˆu2n|s=0=2ωAs4[c13As4+c23(ρ1Aa11)],2ˆFˆunˆvn|s=0=1As4[c14(As4)2As4(ρ1Aa11)(2c13c24)3c23(ρ1Aa11)2],2ˆFˆv2n|s=0=1ωAs4[2c15(As4)22(As4)2(ρ1Aa11)(c14c25)2As4(c24c13)(ρ1Aa11)2+2c23(ρ1Aa11)3],3ˆFˆu3n|s=0=6ω2As4[c16As4+c26(ρ1Aa11)],3ˆFˆu2nˆvn|s=0=ωAs4[2c17(As4)2As4(ρ1Aa11)(6c162c27)6c26(ρ1Aa11)2],3ˆFˆunˆv2n|s=0=1As4[2c18(As4)2+2(As4)2(c282c17)(ρ1Aa11)+6As4(c16c27)(ρ1Aa11)2+6c26(ρ1Aa11)3],3ˆFˆv3n|s=0=6(1+Aa11ρ)ωAs4[c18(As4)3+(As4)2(c28c17)(ρ1Aa11)+As4(c16c27)(ρ1Aa11)2+c26(ρ1Aa11)3],

    and

    2ˆGˆu2n|s=0=2c23ω2As4,2ˆGˆunˆvn|s=0=ωAs4[c24As42c23(ρ1Aa11)],2ˆGˆv2n|s=0=2As4[c25(As4)2As4c24(ρ1Aa11)+c23(ρ1Aa11)2],3ˆGˆu3n|s=0=6c26ω3As4,3ˆGˆu2nˆvn|s=0=ω2As4[2c27As46c26(ρ1Aa11)],3ˆGˆunˆv2n|s=0=ωAs4[2c28(As4)24c27As4(ρ1Aa11)+6c26(ρ1Aa11)2],3ˆGˆv3n|s=0=6(1+Aa11ρ)As4[c26(1+Aa11)2+c27As4(1+Aa11)+c28(As4)2(2(1+Aa11)c26+As4c27)ρ+c26ρ2].

    Based on the aforementioned calculations, we establish the theorem regarding the existence and direction of Neimark-Sacker bifurcation.

    Theorem 5. If (h,α,r,k,K,b,a,s)N.S. and Ω0, then model (1.4) undergoes a Neimark-Sacker bifurcation at the equilibrium point E1(x,y) when the parameter s varies in the vicinity of s4. Moreover, if Ω<0 (respectively, Ω>0), then an attracting (respectively, repelling) closed invariant curve bifurcates from the equilibrium point for s>0 (respectively, s<0).

    Chaos often has detrimental effects on biological systems, disrupting the ecological balance of populations and directly influencing long-term population growth projections. Implementing effective control policies not only safeguards the size of ecological populations but also establishes a strong foundation for the sustainable exploitation of ecological resources [27,28]. This section explores two control methods aimed at effectively managing the chaos generated by model (1.4).

    In this subsection, the state feedback control method [22] will be employed to regulate the chaos exhibited by model (1.4). To achieve this, we introduce the following controlled model.

    {xn+1=xn+hαΓ(α+1)[rxn1+kyn(1xnK)bxnyna+xn]+Sn,yn+1=yn+hαΓ(α+1)[syn(1ynxn)], (4.1)

    which corresponds to model (1.4). The feedback controlling force is defined as

    Sn=p1(xnx)p2(yny), (4.2)

    where (x,y) represents the positive equilibrium point of the model (1.4), and p1, p2 stand for the feedback gains. The Jacobian matrix of the controlled model (4.1) evaluated at the positive equilibrium point E1(x,y) is given by

    J1(x,y)=(1+Aa11p1Aa12p2As1As), (4.3)

    where the variables A, a11, and a12 are defined in Eq (2.2). The corresponding characteristic equation of the Jacobian matrix J1(x,y) is

    λ2(2+Aa11Asp1)λ+(1+Aa11p1)(1As)(Aa12p2)As=0. (4.4)

    Let λ1 and λ2 be the roots of the Eq (4.3), then

    λ1λ2=(1+Aa11p1)(1As)(Aa12p2)As. (4.5)

    The lines of marginal stability l1, l2, and l3 are derived by solving λ1λ2=1, λ1=1 and λ2=±1, respectively. These conditions ensure that |λ1,2|=1. Then, we derive the marginal stability lines as follows:

    l1:(1As)p1Asp2=A2s(a11+a12)+A(a11s), (4.6)
    l2:p1+p2=A(a11+a12), (4.7)
    l3:(2As)p1Asp2=A2s(a11+a12)+2A(a11s)+4. (4.8)

    Therefore, l1, l2, and l3 in the (p1,p2)-plane form a triangular region which leads to |λ1,2|<1.

    Theorem 6. If p1 and p2 lie within a triangular region bounded by the lines l1, l2, and l3, it can be concluded that the model (4.1) is stable.

    Next, we apply the hybrid control approach proposed by [23] to control chaos. The controlled model of (1.4) with the hybrid control approach is depicted below:

    {xn+1=ρxn+ρhαΓ(α+1)[rxn1+kyn(1xnK)bxnyna+xn]+(1ρ)xn,yn+1=ρyn+ρhαΓ(α+1)[syn(1ynxn)]+(1ρ)yn. (4.9)

    where 0<ρ<1, and the controlled strategy in (4.9) combines feedback control and parameter perturbation. By appropriately choosing the controlled parameter ρ, the chaotic behaviors of the equilibrium point (x,y) of the controlled model (4.9) can be accelerated (delayed) or even entirely eliminated. The Jacobian matrix of the controlled model (4.9), evaluated at the positive equilibrium point (x,y), is given by

    J2(x,y)=(1+Aρa11Aρa12Aρs1Aρs), (4.10)

    where the variables A, a11, and a12 are defined in Eq (2.2). Then, the positive equilibrium point (x,y) of the controlled model (4.9) is locally asymptotically stable if the roots of the characteristic polynomial of (4.10) lie within the open unit disk. According to the Jury condition, the equilibrium point of the model remains stable if, and only if, the following conditions are met:

    |2+Aρa11Aρs|<1+(1+Aρa11)(1Aρs)(Aρa12)(Aρs)<2.

    Theorem 7. If |2+Aρa11Aρs|<1+(1+Aρa11)(1Aρs)(Aρa12)(Aρs)<2 can hold, it can be concluded that the model (4.9) is stable.

    In this section, by exploring specific cases of model (1.4), we confirm the theoretical analysis above and discover new intriguing complex dynamic behaviors. Additionally, we validate the effectiveness of linear feedback techniques and hybrid control strategies for chaos control through numerical simulations.

    Take

    h=0.69,α=0.85,r=2.27,k=1.64,K=2.16,b=0.24,a=2.02. (5.1)

    After straightforward calculations, we determine the positive equilibrium point E1=(1.7498,1.7498) and the critical bifurcation value s1=2.8472. The Jacobian matrix of model (1.4) evaluated at E1 with s=s1 is given by

    J(E1)=(0.67330.14972.19651.1965),

    whose characteristic equation is

    F(λ)=λ2+0.5232λ0.4768=0. (5.2)

    The roots of (5.2) are λ1=1 and λ2=0.4768. Moreover, we have β1=0.79590 and β2=1.0080>0. According to Theorem 4, model (1.4) undergoes period-doubling bifurcation at E1 as s passes through s1. This behavior, verified through corresponding bifurcation diagrams shown in Figure 1, utilizes initial conditions (x0,y0)=(1.74,1.74) and varies s in the range [2.7,3.7].

    Figure 1.  Bifurcation diagrams of model (1.4) with parameter values as given in (5.1) and initial conditions (1.74,1.74).

    Take

    h=0.75,α=0.8,r=4.2,k=0.5,K=3.5,b=3.6,a=8.1. (5.3)

    After straightforward calculations, we determine the positive equilibrium point E1=(2.1746,2.1746) and the critical bifurcation value s4=1.1872. The Jacobian matrix of model (1.4) evaluated at E1 with s=s4 is given by

    J(E1)=(0.07120.98841.01260.0126),

    whose characteristic equation is

    F(λ)=λ20.0586λ+1=0. (5.4)

    The roots of (5.4) are λ1,2=0.02930±0.9996i. Moreover, we have d=0.39120 and Ω=1392.0464<0. According to Theorem 5, model (1.4) undergoes Neimark-Sacker bifurcation at E1 as s passes through s4. This behavior, verified through corresponding bifurcation diagrams shown in Figure 2, utilizes initial conditions (x0,y0)=(2.2,2.2) and varies s in the range [1.1,1.6].

    Figure 2.  Bifurcation diagrams of model (1.4) with parameter values as given in (5.3) and initial conditions (2.2,2.2).

    We chose the parameter values as follows:

    h=0.75,α=0.8,r=4.2,k=0.5,K=3.5,b=3.6,a=8.1,s=1.58. (5.5)

    By Figure 3, we can get that the variables xn and yn in the model (1.4) are in a chaotic state.

    Figure 3.  Plots for model (1.4) with parameter values as given in (5.5) and initial conditions (2.2,2.2).

    We use the same parameter values as in (5.5). In Figure 4, the triangular region defined by Theorem 6 bounds the parameters p1 and p2. Inside this region, the chaotic behavior produced by model (1.4) is effectively managed, resulting in asymptotic convergence toward the equilibrium point E1=(2.1746,2.1746).

    Figure 4.  Stability region for the controlled model (4.1).

    In this case, with feedback gains set to p1=0.57 and p2=0.38 and the controller activated at the 3000th iteration of the model, Figure 5 demonstrates the control effect. A chaotic trajectory is successfully stabilized at the equilibrium point E1=(2.1746,2.1746). This indicates that the feedback control approach is effective in mitigating bifurcation and chaos.

    Figure 5.  Plots for controlled model (4.1) with parameter values as given in (5.5), p1=0.57, p2=0.38, and initial conditions (2.2,2.2).

    Finally, using the same parameter values as in (5.5), the Jacobian matrix of the controlled model (4.9), evaluated at E1, is given by:

    (10.9288ρ0.988439ρ1.34765ρ11.34765ρ). (5.6)

    The characteristic polynomial of (5.6) is given by

    λ2(22.27645ρ)λ+2.58377ρ22.27645ρ+1=0. (5.7)

    The roots of (5.7) lie within the open unit disk if, and only if, 0<ρ<0.881058. Additionally, the plots for xn and yn of the controlled model (4.9) are shown in Figure 6 with ρ=0.87. From Figure 6, it is clear that the positive equilibrium point E1 is stable. Therefore, it can be concluded that employing the hybrid control approach is effective in mitigating bifurcation and chaos.

    Figure 6.  Plots for controlled model (4.9) with parameter values as given in (5.5), ρ=0.87, and initial conditions (2.2,2.2).

    In this study, we proposed a fractional Leslie-Gower model with a Holling type Ⅱ functional response and antipredator behavior. By employing the piecewise constant approximation method, we derived the discrete model (1.4) and analyzed its dynamical behavior, including the existence and stability of equilibrium points and the possibility of local bifurcations. The following conclusions can be drawn from our research:

    (1) The discrete model (1.4) has two equilibrium points: E0(K,0) and E1(x,y). The only positive coexistence equilibrium point, E1(x,y), reflects the coexistence of predators and prey.

    (2) Our theoretical analysis and numerical simulations of the positive equilibrium point E1(x,y) indicate that the model undergoes period-doubling bifurcation and Neimark-Sacker bifurcation under specific parameter conditions. Figures 1 and 2 illustrate how these bifurcations can lead to chaotic behavior at E1.

    (3) By applying state feedback control and hybrid control methods, we effectively managed the chaotic behavior generated by the model (1.4), as shown in Figures 36. These interventions mitigated the adverse effects of chaos and bifurcations, consequently enhancing ecosystem resilience.

    This study advances our understanding of the complex dynamics in ecological models influenced by fear effects and provides practical techniques for controlling chaotic behavior in such models. Future work could explore different discretization methods for the model, and new parameters could be chosen to study the influence of various ecological effects on population dynamics.

    Yao Shi: Formal analysis, validation, writing-original draft & editing, visualization, methodology; Zhenyu Wang: Conceptualization, investigation, writing-review & editing, software, supervision. All authors have read and approved the final version of the manuscript for publication.

    This work was supported by the Innovation Foundation of Hebei University of Engineering (Grant Nos. SJ2401002097), the National Natural Science Foundation of China (Grant Nos. 12401519), and the Central Guidance on Local Science and Technology Development Fund of Hebei Province (Grant Nos. 246Z1825G).

    The authors declare no conflict of interest.



    [1] N. Rodríguez, On the global well-posedness theory for a class of PDE models for criminal activity, Phys. D Nonlinear Phenom., 260 (2013), 191–200. https://doi.org/10.1016/j.physd.2012.08.003 doi: 10.1016/j.physd.2012.08.003
    [2] M. Short, M. D'Orsogna, V. Pasour, G. Tita, P. Brantingham, A. Bertozzi, et al., A statistical model of criminal behavior, Math. Mod. Meth. Appl. Sci., 18 (2008), 1249–1267. https://doi.org/10.1142/S0218202508003029 doi: 10.1142/S0218202508003029
    [3] M. Short, A. Bertozzi, P. Brantingham, G. Tita, Dissipation and displacement of hotspots in reaction-diffusion model of crime, Proc. Natl. Acad. Sci. USA, 107 (2010), 3961–3965. https://doi.org/0.1073/pnas.0910921107 doi: 10.1073/pnas.0910921107
    [4] H. Berestycki, J. Wei, M. Winter, Existence of symmetric and asymmetric spikes for a crime hotspot model, SIAM J. Math. Anal., 46 (2014), 691–719. https://doi.org/10.1137/130922744 doi: 10.1137/130922744
    [5] R. Cantrell, C. Cosner, R. Manásevich, Global bifurcation of solutions for crime modeling equations, SIAM J. Math. Anal., 44 (2012), 1340–1358. https://doi.org/10.1137/110843356 doi: 10.1137/110843356
    [6] Y. Gu, Q. Wang, G. Yi, Stationary patterns and their selection mechanism of urban crime models with heterogeneous near-repeat victimization effect, Eur. J. Appl. Math., 28 (2017), 141–178. https://doi.org/10.1017/S0956792516000206 doi: 10.1017/S0956792516000206
    [7] T. Kolokolnikov, M. Ward, J. Wei, The stability of steady-state hot-spot patterns for a reaction-diffusion model of urban crime, Discrete Contin. Dyn. Syst. Ser. B, 19 (2014), 1373–1410. https://doi.org/10.3934/dcdsb.2014.19.1373 doi: 10.3934/dcdsb.2014.19.1373
    [8] D. Lloyd, H. O'Farrell, On localised hotspots of an urban crime model, Phys. D Nonlinear Phenom., 253 (2013), 23–39. https://doi.org/10.1016/j.physd.2013.02.005 doi: 10.1016/j.physd.2013.02.005
    [9] D. Lloyd, N. Santitissadeekorn, M. Short, Exploring data assimilation and forecasting issues for an urban crime model, Eur. J. Appl. Math., 27 (2016), 451–478. https://doi.org/10.1017/S0956792515000625 doi: 10.1017/S0956792515000625
    [10] L. Mei, J. Wei, The existence and stability of spike solutions for a chemotax is system modeling crime pattern formation, Math. Models Methods Appl. Sci., 30 (2020), 1727–1764. https://doi.org/10.1142/S0218202520500359 doi: 10.1142/S0218202520500359
    [11] M. Short, A. Bertozzi, P. Brantingham, Nonlinear patterns in urban crime: hotspots, bifurcations, and suppression, SIAM J. Appl. Dyn. Syst., 9 (2010), 462–483. https://doi.org/10.1137/090759069 doi: 10.1137/090759069
    [12] W. Tse, M. Ward, Hotspot formation and dynamics for a continuum model of urban crime, Eur. J. Appl. Math., 27 (2016), 583–624. https://doi.org/10.1017/S0956792515000376 doi: 10.1017/S0956792515000376
    [13] N. Rodríguez, A. Bertozzi, Local existence and uniqueness of solutions to a PDE model for criminal behavior, Math. Models Methods Appl. Sci., 20 (2010), 1425–1457. https://doi.org/10.1142/S0218202510004696 doi: 10.1142/S0218202510004696
    [14] N. Rodríguez, M. Winkler, On the global existence and qualitative behavior of one-dimensional solutions to a model for urban crime, Eur. J. Appl. Math., 33 (2022), 919–959. https://doi.org/10.1017/S0956792521000279 doi: 10.1017/S0956792521000279
    [15] Q. Wang, D. Wang, Y. Feng, Global well-posedness and uniform boundedness of urban crime models: One-dimensional case, J. Differ. Equations, 269 (2020), 6216–6235. https://doi.org/10.1016/j.jde.2020.04.035 doi: 10.1016/j.jde.2020.04.035
    [16] M. Freitag, Global solutions to a higher-dimensional system related to crime modeling, Math. Meth. Appl. Sci., 41 (2018), 6326–6335. https://doi.org/10.1002/mma.5141 doi: 10.1002/mma.5141
    [17] J. Shen, B. Li, Mathematical analysis of a continuous version of statistical models for criminal behavior, Math. Meth. Appl. Sci., 43 (2020), 409–426. https://doi.org/10.1002/mma.5898 doi: 10.1002/mma.5898
    [18] J. Ahn, K. Kang, J. Lee, Global well-posedness of logarithmic Keller-Segel type systems, J. Differ. Equations, 287 (2021), 185–211. https://doi.org/10.1016/j.jde.2021.03.053 doi: 10.1016/j.jde.2021.03.053
    [19] Y. Tao, M. Winkler, Global smooth solutions in a two-dimensional cross-diffusion system modeling propagation of urban crime, Commun. Math. Sci., 19 (2021), 829–849. https://doi.org/10.4310/CMS.2021.v19.n3.a12 doi: 10.4310/CMS.2021.v19.n3.a12
    [20] M. Winkler, Global solvability and stabilization in a two-dimensional cross-diffusion system modeling urban crime propagation, Ann. Inst. Henri Poincaré, Anal. Non Linéaire, 36 (2019), 1747–1790. https://doi.org/10.1016/j.anihpc.2019.02.004 doi: 10.1016/j.anihpc.2019.02.004
    [21] Y. Jiang, L. Yang, Global solvability and stabilization in a three-dimensional cross-diffusion system modeling urban crime propagation, Acta Appl. Math., 178 (2022). https://doi.org/10.1007/s10440-022-00484-z doi: 10.1007/s10440-022-00484-z
    [22] N. Rodríguez, M. Winkler, Relaxation by nonlinear diffusion enhancement in a two-dimensional cross-diffusion model for urban crime propagation, Math. Models Methods Appl. Sci., 30 (2020), 2105–2137. https://doi.org/10.1142/S0218202520500396 doi: 10.1142/S0218202520500396
    [23] F. Heihoff, Generalized solutions for a system of partial differential equations arising from urban crime modeling with a logistic source term, Z. Für Angew. Math. Phys., 71 (2020). https://doi.org/10.1007/s00033-020-01304-w doi: 10.1007/s00033-020-01304-w
    [24] B. Li, L. Xie, Generalized solution to a 2D parabolic-parabolic chemotaxis system for urban crime: Global existence and large time behavior, submitted for publication, 2022.
    [25] P. Jones, P. Brantingham, L. Chayes, Statistical models of criminal behavior: The effects of law enforcement actions, Math. Models Methods Appl. Sci., 20 (2010), 1397–1423. https://doi.org/10.1142/S0218202510004647 doi: 10.1142/S0218202510004647
    [26] A. Pitcher, Adding police to a mathematical model of burglary, Eur. J. Appl. Math., 21 (2010), 401–419. https://doi.org/10.1017/S0956792510000112 doi: 10.1017/S0956792510000112
    [27] J. Zipkin, M. Short, A. Bertozzi, Cops on the dots in a mathematical model of urban crime and police response, Discrete Contin. Dyn. Syst. Ser. B, 19 (2014), no. 5, 1479–1506. https://doi.org/10.3934/dcdsb.2014.19.1479 doi: 10.3934/dcdsb.2014.19.1479
    [28] W. Tse, M. Ward, Asynchronous instabilities of crime hotspots for a 1-D reaction-diffusion model of urban crime with focused police patrol, SIAM J. Appl. Dyn. Syst., 17 (2018), 2018–2075. https://doi.org/10.1137/17M1162585 doi: 10.1137/17M1162585
    [29] A. Buttenschoen, T. Kolokolnikov, M. Ward, J. Wei, Cops-on-the-dots: the linear stability of crime hotspots for a 1-D reaction-diffusion model of urban crime, Eur. J. Appl. Math., 31 (2020), 871–917. https://doi.org/10.1017/S0956792519000305 doi: 10.1017/S0956792519000305
    [30] B. Li, L. Xie, Global large-data generalized solutions to a two-dimensional chemotaxis system stemming from crime modelling, Discrete Contin. Dyn. Syst. Ser. B, 2022 (2022). https://doi.org/10.3934/dcdsb.2022167 doi: 10.3934/dcdsb.2022167
    [31] N. Rodríguez, Q. Wang, L. Zhang, Understanding the effects of on- and off-hotspot policing: Evidence of hotspot, oscillating, and chaotic activities, SIAM J. Appl. Dyn. Syst., 20 (2021), 1882–1916. https://doi.org/10.1137/20M1359572 doi: 10.1137/20M1359572
    [32] E. Keller, L. Segel, Initiation of slime mold aggregation viewed as an instability, J. Theor. Biol., 26 (1970), 399–415. https://doi.org/10.1016/0022-5193(70)90092-5 doi: 10.1016/0022-5193(70)90092-5
    [33] K. Fujie, Boundedness in a fully parabolic chemotaxis system with strongly singular sensitivity, J. Math. Anal. Appl., 424 (2015), 675–684. https://doi.org/10.1016/j.jmaa.2014.11.045 doi: 10.1016/j.jmaa.2014.11.045
    [34] M. Winkler, Global solutions in a fully parabolic chemotaxis system with singular sensitivity, Math. Methods Appl. Sci., 34 (2011), 176–190. https://doi.org/10.1002/mma.1346 doi: 10.1002/mma.1346
    [35] J. Lankeit, A new approach toward boundedness in a two-dimensional parabolic chemotaxis system with singular sensitivity, Math. Meth. Appl. Sci., 39 (2016), 394–404. https://doi.org/10.1002/mma.3489 doi: 10.1002/mma.3489
    [36] C. Stinner, M. Winkler, Global weak solutions in a chemotaxis system with large singular sensitivity, Nonliear Anal. Real Word Appl., 12 (2011), 3727–3740. https://doi.org/10.1016/j.nonrwa.2011.07.006 doi: 10.1016/j.nonrwa.2011.07.006
    [37] J. Lankeit, M. Winkler, A generalized solution concept for the Keller-Segel system with logarithmic sensitivity: global solvability for large nonradial data, NoDEA-Nonlinear Differ. Equ. Appl., 24 (2017). https://doi.org/10.1007/s00030-017-0472-8 doi: 10.1007/s00030-017-0472-8
    [38] A. Zhigun, Generalised supersolutions with mass control for the Keller-Segel system with logarithmic sensitivity, J. Math. Anal. Appl., 467 (2018), 1270–1286. https://doi.org/10.1016/j.jmaa.2018.08.001 doi: 10.1016/j.jmaa.2018.08.001
    [39] M. Winkler, Unlimited growth in logarithmic Keller-Segel systems, J. Differ. Equations, 309 (2022), 74–97. https://doi.org/10.1016/j.jde.2021.11.026 doi: 10.1016/j.jde.2021.11.026
    [40] M. Winkler, T. Yokota, Stabilization in the logarithmic Keller-Segel system, Nonlinear Anal., 170 (2018), 123–141. https://doi.org/10.1016/j.na.2018.01.002 doi: 10.1016/j.na.2018.01.002
    [41] J. Ahn, Global well-posedness and asymptotic stabilization for chemotaxis system with signal-dependent sensitivity, J. Differ. Equations, 266 (2019), 6866–6904. https://doi.org/10.1016/j.jde.2018.11.015 doi: 10.1016/j.jde.2018.11.015
    [42] Q. Hou, Z. Wang, K. Zhao, Boundary layer problem on a hyperbolic system arising from chemotaxis, J. Differ. Equations, 261 (2016), 5035–5070. https://doi.org/10.1016/j.jde.2016.07.018 doi: 10.1016/j.jde.2016.07.018
    [43] H. Jin, J. Li, Z. Wang, Asymptotic stability of traveling waves of a chemotaxis model with singular sensitivity, J. Differ. Equations, 255 (2013), 193–219. doilinkhttps://doi.org/10.1016/j.jde.2013.04.002 doi: 10.1016/j.jde.2013.04.002
    [44] H. Li, K. Zhao, Initial-boundary value problems for a system of hyperbolic balance laws arising from chemotaxis, J. Differ. Equations, 258 (2015), 302–308. https://doi.org/10.1016/j.jde.2014.09.014 doi: 10.1016/j.jde.2014.09.014
    [45] J. Li, T. Li, Z. Wang, Stability of traveling waves of the Keller-Segel system with logarithmic sensitivity, Math. Models Methods Appl. Sci., 24 (2014), 2819–2849. https://doi.org/10.1142/S0218202514500389 doi: 10.1142/S0218202514500389
    [46] Y. Tao, L. Wang, Z. Wang, Large-time behavior of a parabolic-parabolic chemotaxis model with logarithmic sensitivity in one dimension, Discrete Contin. Dyn. Syst. Ser. B, 18 (2013), 821–845. https://doi.org/10.3934/dcdsb.2013.18.821 doi: 10.3934/dcdsb.2013.18.821
    [47] Z. Wang, Z. Xiang, P. Yu, Asymptotic dynamics in a singular chemotaxis system modeling onset of tumor angiogenesis, J. Differ. Equations, 260 (2016), 2225–2258. https://doi.org/10.1016/j.jde.2015.09.063 doi: 10.1016/j.jde.2015.09.063
    [48] M. Winkler, Renormalized radial large-data solutions to the higher-dimensional Keller-Segel system with singular sensitivity and signal absorption, J. Differ. Equations, 264 (2018), 2310–2350. https://doi.org/10.1016/j.jde.2017.10.029 doi: 10.1016/j.jde.2017.10.029
    [49] M. Winkler, The two-dimensional Keller-Segel system with singular sensitivity and signal absorption: global large-data solutions and their relaxation properties, Math. Models Methods Appl. Sci., 26 (2016), 987–1024. https://doi.org/10.1142/S0218202516500238 doi: 10.1142/S0218202516500238
    [50] B. Li, L. Xie, Generalized solution and its eventual smoothness to a logarithmic Keller-Segel system for criminal activities, submitted for publication, 2022.
    [51] M. Winkler, Large-data global generalized solutions in a chemotaxis system with tensor-valued sensitivities, SIAM J. Math. Anal., 47 (2015), 3092–3115. https://doi.org/dx.doi.org/10.1137/140979708 doi: 10.1137/140979708
    [52] M. Aida, K. Osaka, T. Tsujikawa, M. Mimura, Chemotaxis and growth system with sigular sensitivity function, Nonliear Anal. Real Word Appl., 6 (2005), 323–336. https://doi.org/10.1016/j.nonrwa.2004.08.011 doi: 10.1016/j.nonrwa.2004.08.011
    [53] X. Cao, Global bounded solutions of the higher-dimensional Keller-Segel system under smallness conditions in optimal spaces, Discrete Contin. Dyn. Syst., 35 (2015), 1891–1904. https://doi.org/10.3934/dcds.2015.35.1891 doi: 10.3934/dcds.2015.35.1891
    [54] M. Winkler, Aggregation vs. global diffusive behavior in the higher-dimensional Keller-Segel model, J. Differ. Equations, 248 (2010), 2889–2905. https://doi.org/10.1016/j.jde.2010.02.008 doi: 10.1016/j.jde.2010.02.008
    [55] J. Simon, Compact sets in the space Lp(0,T;B), Ann. Mat. Pura Appl., 146 (1996), 65–96. https://doi.org/10.1007/BF01762360 doi: 10.1007/BF01762360
    [56] O. Ladyzhenskaya, N. Ural'tseva, Linear and Quasilinear Elliptic Equations, Translated from the Russian by Scripta Technica, Inc. Translation editor: Leon Ehrenpreis, Academic Press, New York-London, 1968.
  • This article has been cited by:

    1. Yao Shi, Xiaozhen Liu, Zhenyu Wang, Bifurcation Analysis and Chaos Control of a Discrete Fractional-Order Modified Leslie–Gower Model with Nonlinear Harvesting Effects, 2024, 8, 2504-3110, 744, 10.3390/fractalfract8120744
    2. Rutaba Saleem, Ahmad Javid, Interaction solutions and qualitative analysis of the Painlevé integrable (2+1)-dimensional extended Kadomtsev–Petviashvili equation, 2025, 128, 11100168, 611, 10.1016/j.aej.2025.05.013
  • Reader Comments
  • © 2023 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(1947) PDF downloads(98) Cited by(6)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog