Research article Special Issues

Dynamical properties of a novel one dimensional chaotic map


  • Received: 09 November 2021 Revised: 24 December 2021 Accepted: 31 December 2021 Published: 07 January 2022
  • In this paper, a novel one dimensional chaotic map K(x)=μx(1x)1+x, x[0,1],μ>0 is proposed. Some dynamical properties including fixed points, attracting points, repelling points, stability and chaotic behavior of this map are analyzed. To prove the main result, various dynamical techniques like cobweb representation, bifurcation diagrams, maximal Lyapunov exponent, and time series analysis are adopted. Further, the entropy and probability distribution of this newly introduced map are computed which are compared with traditional one-dimensional chaotic logistic map. Moreover, with the help of bifurcation diagrams, we prove that the range of stability and chaos of this map is larger than that of existing one dimensional logistic map. Therefore, this map might be used to achieve better results in all the fields where logistic map has been used so far.

    Citation: Amit Kumar, Jehad Alzabut, Sudesh Kumari, Mamta Rani, Renu Chugh. Dynamical properties of a novel one dimensional chaotic map[J]. Mathematical Biosciences and Engineering, 2022, 19(3): 2489-2505. doi: 10.3934/mbe.2022115

    Related Papers:

    [1] Rebecca Everett, Kevin B. Flores, Nick Henscheid, John Lagergren, Kamila Larripa, Ding Li, John T. Nardini, Phuong T. T. Nguyen, E. Bruce Pitman, Erica M. Rutter . A tutorial review of mathematical techniques for quantifying tumor heterogeneity. Mathematical Biosciences and Engineering, 2020, 17(4): 3660-3709. doi: 10.3934/mbe.2020207
    [2] Ziyou Zhuang . Optimization of building model based on 5G virtual reality technology in computer vision software. Mathematical Biosciences and Engineering, 2021, 18(6): 7936-7954. doi: 10.3934/mbe.2021393
    [3] Luca Bertelli, Frédéric Gibou . Fast two dimensional to three dimensional registration of fluoroscopy and CT-scans using Octrees on segmentation maps. Mathematical Biosciences and Engineering, 2012, 9(3): 527-537. doi: 10.3934/mbe.2012.9.527
    [4] Christian Engwer, Markus Knappitsch, Christina Surulescu . A multiscale model for glioma spread including cell-tissue interactions and proliferation. Mathematical Biosciences and Engineering, 2016, 13(2): 443-460. doi: 10.3934/mbe.2015011
    [5] Yu Lei, Zhi Su, Xiaotong He, Chao Cheng . Immersive virtual reality application for intelligent manufacturing: Applications and art design. Mathematical Biosciences and Engineering, 2023, 20(3): 4353-4387. doi: 10.3934/mbe.2023202
    [6] Wencong Zhang, Yuxi Tao, Zhanyao Huang, Yue Li, Yingjia Chen, Tengfei Song, Xiangyuan Ma, Yaqin Zhang . Multi-phase features interaction transformer network for liver tumor segmentation and microvascular invasion assessment in contrast-enhanced CT. Mathematical Biosciences and Engineering, 2024, 21(4): 5735-5761. doi: 10.3934/mbe.2024253
    [7] Christophe Prud'homme, Lorenzo Sala, Marcela Szopos . Uncertainty propagation and sensitivity analysis: results from the Ocular Mathematical Virtual Simulator. Mathematical Biosciences and Engineering, 2021, 18(3): 2010-2032. doi: 10.3934/mbe.2021105
    [8] Bethan Morris, Lee Curtin, Andrea Hawkins-Daarud, Matthew E. Hubbard, Ruman Rahman, Stuart J. Smith, Dorothee Auer, Nhan L. Tran, Leland S. Hu, Jennifer M. Eschbacher, Kris A. Smith, Ashley Stokes, Kristin R. Swanson, Markus R. Owen . Identifying the spatial and temporal dynamics of molecularly-distinct glioblastoma sub-populations. Mathematical Biosciences and Engineering, 2020, 17(5): 4905-4941. doi: 10.3934/mbe.2020267
    [9] Duane C. Harris, Giancarlo Mignucci-Jiménez, Yuan Xu, Steffen E. Eikenberry, C. Chad Quarles, Mark C. Preul, Yang Kuang, Eric J. Kostelich . Tracking glioblastoma progression after initial resection with minimal reaction-diffusion models. Mathematical Biosciences and Engineering, 2022, 19(6): 5446-5481. doi: 10.3934/mbe.2022256
    [10] Liang-Sian Lin, Susan C Hu, Yao-San Lin, Der-Chiang Li, Liang-Ren Siao . A new approach to generating virtual samples to enhance classification accuracy with small data—a case of bladder cancer. Mathematical Biosciences and Engineering, 2022, 19(6): 6204-6233. doi: 10.3934/mbe.2022290
  • In this paper, a novel one dimensional chaotic map K(x)=μx(1x)1+x, x[0,1],μ>0 is proposed. Some dynamical properties including fixed points, attracting points, repelling points, stability and chaotic behavior of this map are analyzed. To prove the main result, various dynamical techniques like cobweb representation, bifurcation diagrams, maximal Lyapunov exponent, and time series analysis are adopted. Further, the entropy and probability distribution of this newly introduced map are computed which are compared with traditional one-dimensional chaotic logistic map. Moreover, with the help of bifurcation diagrams, we prove that the range of stability and chaos of this map is larger than that of existing one dimensional logistic map. Therefore, this map might be used to achieve better results in all the fields where logistic map has been used so far.



    Let (M,g) be a complete Riemannian manifold with nonnegative Ricci tensor. We are going to consider bounded "eternal" classical solutions uC2,1(M×R) of the semilinear heat equation

    ut=Δu+|u|p (1.1)

    for p>1 and T+, with particular attention to "triviality", that is, to conditions forcing the solutions to be identically zero.

    A reason for the interest in such eternal (or also "ancient") solutions is that they typically arise as blow-up limits when the solutions of semilinear parabolic equations in time intervals as [0,T) develop a singularity at a certain time TR, i.e., the solution u becomes unbounded as t goes to T.

    In the Euclidean space, it is well known (see [14] and [7, Proposition B]) that noncostant, positive global radial (static, solution of the elliptic problem) solutions on Rn×R exist for the "critical" or any "supercritical" exponent ppS=n+2n2 (this latter is the critical Sobolev exponent), hence providing a counterexample to triviality in such case. In the non-static (parabolic) situation, while triviality of eternal radial (parabolic) positive solutions can be shown in the range of subcritical exponents 1<p<pS, the same expected result for general (not necessarily radial) solutions is known only in the range 1<p<n(n+2)(n1)2 (see Remark 1.2 below). Indeed, such triviality when n(n+2)(n1)2p<pS is a quite long standing open problem (see [10,11]), which might have been finally solved by Quittner in the preprint [12], appeared during the redaction of this work.

    We mention that these triviality issues for eternal (and also ancient) solutions have been recently addressed and partially extended to the cases of compact or bounded geometry Riemannian manifolds, by the first and third author in [4,5].

    Our aim in this paper is to show the following triviality theorem for eternal solutions, "pointwise" monotone in time (mentioned without proof in [11, Remark 4.3 (b)] for M=Rn).

    Theorem 1.1. Let (M,g) be a complete Riemannian manifold of dimension n5 with nonnegative Ricci tensor. Let uC2,1(M×R) be a bounded eternal solution of Eq (1.1) with ut0 and 1<p<pS. Then, u0.

    Remark 1.2. We observe that, essentially with the same proof of [2], we can prove the following result, extending to manifolds the analogous one in Rn.

    Let (M,g) be a complete Riemannian manifold of dimension n3 with nonnegative Ricci tensor. Let uC2,1(M×R) be an eternal solution of Eq (1.1) with

    1<p<n(n+2)(n1)2.

    Then, u0.

    Indeed, as observed in Remark 2.4, Lemma 3.1 in [3] and thus also Lemma 3 and Lemma 4 in [2], hold true on any Riemannian manifolds with nonnegative Ricci tensor (since they are based only on Bochner inequality). Moreover, thanks to [5,Theorem 2.6], we know that the solutions are either strictly positive or identically zero. In particular, using good cutoff functions that can be constructed on manifolds with nonnegative Ricci tensor (see Remark 2.9), from Lemma 4 in [2], given (x,t)M×R, for every ρ>0 and every

    1<p<n(n+2)(n1)2,

    there exists a constant C=C(n,p) such that

    B(x,ρ)×(tρ2,t+ρ2)u2pCρ24pp1μ(B(x,ρ)),

    where B(x,ρ)M is the metric ball with center x and radius ρ and μ the canonical Riemannian measure of (M,g). By Bishop-Gromov inequality (see [6]), we then obtain

    B(x,ρ)×(tρ2,t+ρ2)u2pCρ24pp1μ(B(x,ρ))Cρn+24pp1.

    Since p<pS, letting ρ+ we conclude that u=0 in the whole M×R.

    First of all, using [5, Theorem 2.6], we have that all solutions are either strictly positive or identically zero. The following Hamilton-type gradient estimate can be shown analogously to Lemma 3.1 in [5] and will be crucial in the proof of the integral estimate in Proposition 2.5.

    Lemma 2.1. Let uC2,1(Ω×(0,T)), ΩRn, be a positive bounded solution of Eq (1.1) with 0<uD and let ρ>0 such that ¯B(¯x,2ρ)×[¯t4ρ2,¯t+4ρ2]Ω×(0,T). Then, for n1 and p>1, there exists a constant C=C(n,p)>0, such that

    |u|uC(1ρ+pDp1)(1+logDu)

    in Q=B(¯x,ρ)×(¯tρ2,¯t+ρ2). In particular, for every δ>0, there exists a constant C=C(n,p,δ)>0, such that

    |u|C(1ρ+pDp1)u1δ

    in Q.

    To keep track of the dependencies of the constants, we define

    K=C(1ρ+pDp1) (2.1)

    where C=C(n,p,δ)>0 is given in the previous proposition. In particular, we have the gradient estimate

    |u|Ku1δ. (2.2)

    Notice that K actually depends only on n,p,δ,D and ρ0, if ρ is larger than some fixed constant ρ0>0.

    Remark 2.2. The same result also holds on complete Riemannian manifolds with nonnegative Ricci tensor (see [5]).

    We recall the following technical lemma [3, Lemma 3.1].

    Lemma 2.3. Let U be an open set of Rn. Then, for any positive function wC2(U), any nonnegative ηC2c(U), any real numbers d,mR such that dm+2, the following inequality holds:

    2(nm)d(n1)(m2+d2)4nUηwm2|w|4n1nUηwm(Δw)22(n1)m+(n+2)d2nUηwm1|w|2Δwm+d2Uwm1|w|2w,η+UwmΔww,η+12Uwm|w|2Δη.

    Remark 2.4. The proof is based on the Bochner formula

    12Δ|f|2=|2f|2+f,Δf1n(Δf)2+f,Δf,

    hence, the conclusion of Lemma 2.3 is also true on Riemannian manifolds with nonnegative Ricci tensor, since this inequality holds in such spaces (see [9], for instance).

    From the previous lemma and the pointwise gradient inequality (2.2), we can show the following integral estimate which generalizes the one of Lemma 3 in [2] (we will use the same notation and line of proof).

    Proposition 2.5. Let ΩRn, with n3 and uC2,1(Ω×(0,T)) be a solution of Eq (1.1), with 1<p<pS and such that 0<uD. We assume that for some ρ>0 we have ¯B(¯x,2ρ)×[¯t4ρ2,¯t+4ρ2]Ω×(0,T) and we let Q=B(¯x,ρ)×(¯tρ2,¯t+ρ2), with ρ larger than some ρ0>0. Suppose that ζC2c(Q) takes values in [0,1] and let r>4, δ(0,1/100), then there exist constants m=m(n,p)>2n2 and C=C(n,p,r,δ,D,ρ0)>0 such that

    Qζrumu2t+Qζrum2|u|4+Qζrum+p1|u|2 (2.3)
    CQζr4um+2(1δ)(|ζ|2+|ζt|+(Δζ)2+|ζ|4+ζ2t)+CQζrum+12δ|ut| (2.4)

    where ut=ut0. In particular, the constants m and C are independent of ρ>ρ0.

    Proof. By sake of clarity, we set U=B(¯x,ρ), t1=¯tρ2 and t2=¯t+ρ2.

    Applying Lemma 2.3 to w=u(,t) and η=ζr(,t), for any t[t1,t2] and any reals d,m with dm+2, we get

    2(nm)d(n1)(m2+d2)4nUζrum2|u|4n1nUζrum(Δu)22(n1)m+(n+2)d2nUζrum1|u|2Δum+d2Uum1|u|2u,(ζr)+UumΔuu,(ζr)+12Uum|u|2Δ(ζr).

    Hence, substituting Δu=utup, we obtain

    2(nm)d(n1)(m2+d2)4nUζrum2|u|4n1nUζrumutΔu+n1nUζrum+pΔu2(n1)m+(n+2)d2nUζrum1|u|2ut+2(n1)m+(n+2)d2nUζrum+p1|u|2m+d2Uum1|u|2u,(ζr)+Uumutu,(ζr)Uum+pu,(ζr)+12Uum|u|2Δ(ζr). (2.5)

    Integrating by parts in the two integrals in the second line above, we get

    n1nUζrumutΔu+n1nUζrum+pΔu=n1nUζrumu,ut+m(n1)nUζrum1|u|2ut+n1nUumutu,(ζr)n1nUum+pu,(ζr)(m+p)(n1)nUζrum+p1|u|2.

    Thus, substituting and setting

    X=Uζrumu,ut,Y=Uζrum1|u|2ut,
    Z=Uumutu,(ζr),V=Uum+pu,(ζr),
    W=Uum1|u|2u,(ζr),R=Uum|u|2Δ(ζr),

    the inequality (2.5) becomes

    2(nm)d(n1)(m2+d2)4nUζrum2|u|4+n1nX+m(n1)nY+n1nZn1nV(m+p)(n1)nUζrum+p1|u|22(n1)m+(n+2)d2nY+2(n1)m+(n+2)d2nUζrum+p1|u|2m+d2W+ZV+12R,

    Hence, rearranging and simplifying, we conclude

    αUζrum2|u|4+βUζrum+p1|u|2n1nX+(n+2)d2nY+1nZ1nV+m+d2W+12R, (2.6)

    where

    α=2(nm)d(n1)(m2+d2)4nandβ=(n+2)d2(n1)p2n.

    We now choose

    m=dn1

    obtaining

    α=[2(n1)(n2)d]d4(n1)andβ=(n+2)d2(n1)p2n.

    It is then easy to see that both constants α and β are positive if d satisfies

    2(n1)pn+2<d<2(n1)n2 (2.7)

    which is a meaningful condition, since p<pS=n+2n2. Thus, we set d to be equal to the average mean of the two values, that is

    d=(n1)pn+2+n1n2

    and consequently

    m=m(n,p)=pn+21n2, (2.8)

    which implies dm+2. Indeed, d=m+2 if and only if p=(n+2)(n4)n(n2) but this value is always smaller than 1, for nN. Moreover, from the right inequality (2.7) we have

    m=dn1>2n2.

    To get the thesis we now bound some of the right-hand side terms in inequality (2.6).

    By using Young's inequality, for any ε>0 and r>4, there exists C(ε,r)>0 such that

    |V|εUζrum+p1|u|2+C(ε)Uζr2um+p+1|ζ|2, (2.9)
    WεUζrum2|u|4+C(ε)Uζr4um+2|ζ|4, (2.10)
    R=r(r1)Uζr2um|u|2|ζ|2+rUζr1um|u|2ΔζεUζrum2|u|4+C(ε)Uζr4um+2|ζ|4+C(ε)Uζr2um+2(Δζ)2εUζrum2|u|4+C(ε)Uζr4um+2|ζ|4+C(ε)Uζr4um+2(Δζ)2=εUζrum2|u|4+C(ε)Uζr4um+2(|ζ|4+(Δζ)2) (2.11)

    since ζ1 everywhere. Moreover,

    X=Uζrumu,ut=dfdtPm2Y (2.12)

    with

    f=12Uζrum|u|2andP=12Uum|u|2(ζr)t.

    Then, by equation ut=Δu+up and integrating by parts the Laplacian from the second to the third line as we did before, we get

    Uζrumu2t=Uζrumut(Δu+up)=Uζrum+put+UζrumutΔu=dgdtSZXmY=d(gf)dt+PZSm2Y, (2.13)

    where

    g=1m+p+1Uζrum+p+1andS=1m+p+1Uum+p+1(ζr)t,

    which implies

    Z=d(gf)dt+PSm2YUζrumu2t. (2.14)

    We bound P as follows,

    P=r2Uζr1ζtum|u|2εUζrum2|u|4+C(ε,r)Uζr2um+2ζ2t. (2.15)

    Substituting equalities (2.12), (2.14) and inequalities (2.9), (2.10), (2.11) into estimate (2.6), we obtain

    αUζrum2|u|4+βUζrum+p1|u|2n1n(dfdtPm2Y)+(n+2)d2nY+1n(d(gf)dt+PSm2YUζrumu2t)+εnUζrum+p1|u|2+(m+d+1)ε2Uζrum2|u|4+CUζr4um+2((Δζ)2+|ζ|4)+CUζr2um+p+1|ζ|2dfdt+1ndgdt+((n1)m2n+(n+2)d2nm2n)Y1nUζrumu2t+εnUζrum+p1|u|2+(m+d+3)ε2Uζrum2|u|4+CUζr4um+2((Δζ)2+|ζ|4+ζ2t)+CUζr2um+p+1(|ζ|2+|ζt|)

    for some constant C=C(n,p,r,ε), where we estimated S simply taking the modulus of the integrand and used inequality (2.15) to deal with P. We notice that the coefficient of the term Y is given by

    (n1)m2n+(n+2)d2nm2n=nd2(n1)>0

    as m=dn1.

    Taking ε small enough and then "absorbing" the two integrals in u in the left side of the inequality, we can conclude that there exists a constant C1 depending only on n,p,r, such that for every t(t1,t2) we have

    α2Uζrum2|u|4+β2Uζrum+p1|u|2dfdt+1ndgdt+nd2(n1)Y1nUζrumu2t.+C1Uζr4um+2((Δζ)2+|ζ|4+ζ2t)+C1Uζr2um+p+1(|ζ|2+|ζt|)

    (notice that possibly varying the constant C1, instead of the constants α/2 and β/2 in front of the first two integrals we could have chosen αδ and βδ, for any δ>0).

    Integrating this inequality between t1 and t2 and observing that f(ti)=g(ti)=0, for i=1,2, since ζC2c(Q), we get

    α2Qζrum2|u|4+β2Qζrum+p1|u|2C1Qζr4um+2((Δζ)2+|ζ|4+ζ2t)+C1Qζr2um+p+1(|ζ|2+|ζt|)+nd2(n1)t2t1Y1nQζrumu2t. (2.16)

    Finally, in order to estimate t2t1Y, we make use of the gradient estimate (2.2) from Lemma 2.1. For any δ>0 we have

    Y=Uζrum1|u|2utUζrum1|u|2u+tK2Uζrum+12δu+t,

    where K is defined in formula (2.1) and we set ut=u+t+ut, u+t=ut00 and ut=ut00. It follows

    t2t1YK2Ut2t1ζrum+12δu+t=K2Ut2t1ζrum+12δutK2Ut2t1ζrum+12δut=K2Ut2t1(ζrum+2(1δ))tK2Ut2t1um+2(1δ)(ζr)tK2Ut2t1ζrum+12δutrK2Qζr1um+2(1δ)|ζt|K2Qζrum+12δut (2.17)

    since ζ(x,ti)=0, for i=1,2 and every xU.

    Substituting in inequality (2.16) we get that there exists a constant C2=C2(n,p,r,δ,D,ρ0) (since we have seen after Lemma 2.1 that K depends on n,p,δ,D and ρ0) such that

    1nQζrumu2t+α2Qζrum2|u|4+β2Qζrum+p1|u|2C2(Qζr2um+p+1(|ζ|2+|ζt|)+Qζr4um+2((Δζ)2+|ζ|4+ζ2t))+C2Qζr1um+2(1δ)|ζt|C2Qζrum+12δut.

    Since uD, 0ζ1 and m+p+1>m+2>m+2(1δ)>0, for δ sufficiently small, as m>2n2 and n3, we have um+p+1Cum+2(1δ), then there exists a positive constant C3=C3(n,p,r,δ,D,ρ0) such that

    Qζrumu2t+Qζrum2|u|4+Qζrum+p1|u|2C3Qζr4um+2(1δ)(|ζ|2+|ζt|+(Δζ)2+|ζ|4+ζ2t)C3Qζrum+12δut,

    which is the thesis.

    Remark 2.6. The same proposition also holds for n=1,2, for every p>1, considering m=0, d=1, if n=1 and m=d with any d>2(n1)p/(n+2), when n=2.

    Now we see that the estimate of this proposition implies an interior integral estimate on u2p+m (see [2,Lemma 4]).

    Lemma 2.7. In the same setting and with m as in Proposition 2.5, if r>2(2p+m)p1+δ there exists a constant C=C(n,p,r,δ,D,ρ0)>0 such that

    Qζru2p+mCQζr2(2+m)p1+δ(|ζ|2+|ζt|+(Δζ)2+|ζ|4+ζ2t)2p+m2p2(1δ)+CQζrum+12δ|ut|. (2.18)

    Proof. Multiplying Eq (1.1) by ζrum+p and integrating on U, we obtain

    Uζru2p+m=Uζrum+putUζrum+pΔu=dgdtS+V+(m+p)Uζrum+p1|u|2.

    Estimating S and V as in the proof of Proposition 2.5, then integrating between t1 and t2 and using estimate (2.4), we get

    Qζru2p+mCQζr4um+2(1δ)(|ζ|2+|ζt|+(Δζ)2+|ζ|4+ζ2t)+CQζrum+12δ|ut|

    with a constant C=C(n,p,r,δ,D,ρ0). Observing that

    2p+m>m+2(1δ),

    the conclusion of the lemma follows easily by means of Young's inequality.

    We then have the following integral decay estimate if u is "pointwise" monotone nondecreasing in time, that is ut0 everywhere.

    Proposition 2.8. Let ΩRn, n3, and uC2,1(Ω×(0,T)) be a positive solution of Eq (1.1) with p<pS=n+2n2, such that 0<uD and ut0. Let 0<¯t<T, ¯xΩ, m>2n2 as in (2.8), δ(0,1/100), r>2(2p+m)p1+δ and ρ0>0. Then there exists a constant C=C(n,p,r,δ,D,ρ0)>0 such that, for every ρ>ρ0 with ¯B(¯x,2ρ)×[¯t4ρ2,¯t+4ρ2]Ω×(0,T), letting Q(ρ)=B(¯x,ρ)×(¯tρ2,¯t+ρ2), there holds

    Q(ρ/2)u2p+mCρn+22p+mp1+δ.

    Proof. Since ut=0, we obtain the thesis by applying Lemma 2.7, choosing ζ(x,t)=φ(x)ψ(t) with φC2c(B(¯x,ρ)) and ψC1c(¯tρ2,¯t+ρ2), both taking values in [0,1], such that φ=1 on B(¯x,ρ/2), ψ=1 on [¯tρ2/4,¯t+ρ2/4] and

    |Δφ|+|φ|2+|ψt|C(n)ρ2. (2.19)

    Remark 2.9. Given a complete Riemannian manifold with nonnegative Ricci tensor, it is possible to construct cutoff functions satisfying conditions (2.19) (see [8], Theorem 2.2], for instance). In particular, using Remarks 2.2 and 2.4, it is then clear that the same proof goes through also in this case.

    Proposition 2.10. Let uC2,1(Ω×(0,T)) be a positive bounded solution of Eq (1.1) with p<pS=n+2n2, such that 0<uD and ut0, where Ω is an open subset of a Riemannian manifold (M,g) with nonnegative Ricci tensor of dimension n3. Let 0<¯t<T, ¯xΩ, m>2n2 as in (2.8), δ(0,1/100), r>2(2p+m)p1+δ and ρ0>0. Then there exists a constant C=C(n,p,r,δ,D,ρ0)>0 such that, for every ρ>ρ0 with ¯B(¯x,2ρ)×[¯t4ρ2,¯t+4ρ2]Ω×(0,T), letting Q(ρ)=B(¯x,ρ)×(¯tρ2,¯t+ρ2), there holds

    Q(ρ/2)u2p+mCρ22p+mp1+δμ(B(¯x,ρ)),

    where B(¯x,ρ)M is the metric ball with center ¯x and radius ρ and μ the canonical Riemannian measure of (M,g).

    Remark 2.11. The assumption that u is "pointwise" monotone nondecreasing in time (that is, ut is identically zero, hence the second integral in the right hand side of formula (2.18) in Lemma 2.7 vanishes) can be weakened a little. Indeed, notice that if δ>0 is sufficiently small, we have

    m+12δ>2n2+12δ>0,

    when n>4, hence, using Young's inequality, we get

    CQζrum+12δ|ut|12Qζru2p+m+¯CQζr|ut|2p+m2p1+2δ

    where C is the constant in formula (2.18). Then, we obtain

    12Qζru2p+mCQζr2(2+m)p1+δ(|ζ|2+|ζt|+(Δζ)2+|ζ|4+ζ2t)2p+m2p2(1δ)+¯CQ|ut|2p+m2p1+2δ,

    being 0ζ1.

    Hence, if there exists a constant ˜C=˜C(n,p,r,δ,D,ρ0), such that, for ρ>ρ0 with ¯B(¯x,2ρ)×[¯t4ρ2,¯t+4ρ2]Ω×(0,T), there holds

    Q|ut|2p+m2p1+2δ˜Cρ22p+mp1+δμ(B(¯x,ρ)), (2.20)

    then we have the same conclusion of the above proposition, along the same line of proof.

    We are going to apply the following estimates of Aronson-Serrin in [1], which have been extended to manifolds with Ricci tensor bounded from below by Saloff-Coste in [13]. For ΩRn, let uC2,1(Ω×(0,T)) be a positive bounded solution of Eq (1.1), with p<pS=n+2n2. Following their notation, we set

    d=up1hence, ut=Δu+up=Δu+du.

    Fixed a pair (¯x,¯t)Ω×(0,T) and ρ>0, we consider the parabolic cylinder

    ˜Q(ρ)=B(¯x,ρ)×(¯tρ2,¯t)Q(ρ)

    (where Q(ρ)=B(¯x,ρ)×(¯tρ2,¯t+ρ2) was the set defined in the previous section) and we will use the symbol q,ρ to denote the Lq-norm of a function in ˜Q(ρ). In [1] the following weak Harnack estimate is proved (we mention that the assumption q>n+22 guarantees the validity of the assumption (4) in [1], necessary for Lemma 3).

    Theorem 3.1 ([1,Theorem 2]). For ΩRn, let uC2,1(Ω×(0,T)) be a solution of equation ut=Δu+du. Suppose that ˜Q(3ρ)Ω×(0,T), then, for any q>n+22, there exists a constant C=C(n,q,ρ,d)>0 such that

    |u(x,t)|Cρn+22u2,3ρ

    for every (x,t)˜Q(ρ).

    Inspecting carefully the proof and keeping track of the dependencies of the constants, it follows that the constant C can be made more explicit, that is,

    C=C(n)ρ1n+22qd1/2q,3ρ,

    hence, in the same hypotheses of this theorem, we can conclude

    |u(x,t)|Cρ1n+22qn+22d1/2q,3ρu2,3ρCρ|˜Q(3ρ)|q+12qd1/2q,3ρu2,3ρ,

    for every (x,t)˜Q(ρ), with C=C(n), as |˜Q(3ρ)|=(3ρ)2|B(¯x,3ρ)|=ωn3n+2ρn+2 is the Lebesgue measure of the parabolic cylinder ˜Q(3ρ).

    The same estimate can be extended to Riemannian manifolds (M,g) with nonnegative Ricci tensor, see [13, Theorem 5.6] (with K=0), obtaining

    |u(x,t)|Cρ1qμ(B(¯x,3ρ))q+12qd1/2q,3ρu2,3ρ,

    where C=C(n,M) and B(¯x,ρ)M is the metric ball with center ¯x and radius ρ and μ the canonical Riemannian measure of (M,g).

    Setting then d=up1 in this estimate we get the following corollary.

    Corollary 3.2. Let (M,g) be a complete Riemannian manifold with nonnegative Ricci tensor. Let uC2,1(M×R) be a nonnegative solution of equation ut=Δu+up with p>1, then for any q>n+22 there exists a positive constant C=C(n,M)>0 such that in ˜Q(ρ) we have the estimate,

    u(x,t)Cρ1qμ(B(¯x,3ρ))q+12qup12q(p1),3ρu2,3ρ

    for every ρ>0.

    We are then ready to show Theorem 1.1.

    Proof of Theorem 1.1. Let (M,g) be a complete Riemannian manifold with nonnegative Ricci tensor. Let uC2,1(M×R) be a bounded eternal solution of equation (1.1) such that ut0, with 1<p<pS=n+2n2 and n3. Then, by Theorem 2.6 in [5], the function u is identically zero or positive everywhere. The conclusion in the case 1<p<pBV=n(n+2)(n1)2 is a consequence of the estimates in [2] (see Theorem A in [11]) for the Euclidean setting, which can be extended to Riemannian manifolds with nonnegative Ricci tensor as observed in Remark 1.2, in the introduction. Thus, we assume

    n(n+2)(n1)2p<pS.

    We choose m=dn1>2n2, as in (2.8) and we observe that

    2p+m2n(n+2)(n1)22n2=2n3n22n1(n2)(n1)2=2+2(3n27n+1)(n2)(n1)2>2

    when n3.

    Now let γ>0 such that

    q=(1+γ)2p+mp1>n+22 (3.1)

    and

    1+2(1+γ)1+2γ>p, (3.2)

    (we will discuss the existence of such γ later on).

    It follows that

    limδ0+(1p12(1+γ)(p1+δ)1p1+δ)=1+2γ2(1+γ)1p1<0,

    hence, the quantity

    σ=1p12(1+γ)(p1+δ)1p1+δ (3.3)

    is negative, choosing δ>0 sufficiently small.

    Now, fixing (¯x,¯t)M×R and ρ>0, from Holder's inequality and Proposition 2.10 (notice that ˜Q(3ρ)B(¯x,3ρ)×[¯t9ρ2,¯t+9ρ2]), we obtain

    u2,3ρ=(˜Q(3ρ)u2)1/2(ρ2μ(B(¯x,3ρ)))2p+m22(2p+m)(˜Q(3ρ)u2p+m)12p+mC(ρ2μ(B(¯x,3ρ)))1/2ρ1p1+δ (3.4)

    for every ρ>0.

    Then, using the fact that u is bounded and γ>0, again by Holder's inequality and Proposition 2.10 we have

    up1q(p1),3ρ=(˜Q(3ρ)uq(p1))1/q=(˜Q(3ρ)u(1+γ)(2p+m))1/qC(˜Q(3ρ)u2p+m)1/qC(ρ2μ(B(¯x,3ρ)))1/qρp1(1+γ)(p1+δ).

    By Corollary 3.2, since q>n+22, there exists a positive constant C=C(n,M) such that in ˜Q(ρ) we have

    u(x,t)Cρ1qμ(B(¯x,3ρ))q+12qup12q(p1),3ρu2,3ρ

    hence, by estimates (3.4) and (3.5), we obtain that for every (x,t)˜Q(ρ) there holds

    u(x,t)Cρ1qμ(B(¯x,3ρ))q+12qup12q(p1),3ρu2,3ρCρ1qμ(B(¯x,3ρ))q+12q(ρ2μ(B(¯x,3ρ)))12qρp12(1+γ)(p1+δ)(ρ2μ(B(¯x,3ρ)))12ρ1p1+δ=Cρ1p12(1+γ)(p1+δ)1p1+δ=Cρσ

    where σ<0 is defined by formula (3.3).

    Since u is nonnegative, letting ρ+, we obtain u(x,t)=0 for every (x,t)M×(,¯t) and being ¯tR arbitrary, u=0 everywhere.

    Finally, one can check that, if n5 and p<pS, there exists γ>0 such that conditions (3.1) and (3.2) are satisfied. In fact, in order to satisfy condition (3.1), thanks to the fact that 2p+m2 and p<n+2n2, for any n3 we can choose a γ>0 such that

    γ>n+24(n+2n21)1=4n2.

    On the other hand, to fulfill condition (3.2) it is sufficient to choose γ that satisfies

    11+2γ>p2,

    which, for p<n+2n2 and n6, is true for any γ>0. Finally, for n=5, it is sufficient to choose γ=1.

    The authors declare no conflict of interest.



    [1] J. Gleick, Chaos: Making a New Science, Viking Books, New York, 1997.
    [2] I. Prigogine, I. Stengers, A. Toffler, Order Out of Chaos: Man's New Dialogue with Nature, Bantam Books, New York, 1984.
    [3] E. Ott, Chaos in Dynamical Systems, Cambridge University Press, Cambridge, 2002.
    [4] E. N. Lorenz, Deterministic nonperiodic flow, J. Atmos Sci., 20 (1963), 130–141. https://doi.org/10.1175/1520-0469(1963)020<0130:DNF>2.0.CO;2. doi: 10.1175/1520-0469(1963)020<0130:DNF>2.0.CO;2
    [5] R. M. May, Simple mathematical models with very complicated dynamics, Nature, 261 (1976), 459–465. https://doi.org/10.1038/261459a0. doi: 10.1038/261459a0
    [6] M. F. Barnsley, Fractals Everywhere, 2nd edition, Revised with the assistance of and a foreword by Hawley Rising, Academic Press Professional, Boston, 1993.
    [7] F. M. Atay, J. Jost, A. Wende, Delays, connection topology, and synchronization of coupled chaotic maps, Phys. Rev. Lett., 92 (2004), 1–5. https://doi.org/10.1103/PhysRevLett.92.144101. doi: 10.1103/PhysRevLett.92.144101
    [8] K. M. Cuomo, A. V. Oppenheim, Circuit implementation of synchronized chaos with applications to communications, Phys. Rev. Lett., 71 (1993), 65–68. https://doi.org/10.1103/PhysRevLett.71.65. doi: 10.1103/PhysRevLett.71.65
    [9] I. Klioutchnikova, M. Sigovaa, N. Beizerov, Chaos Theory in Finance, Procedia Comput. Sci., 119 (2017), 368–375. https://doi.org/10.1016/j.procs.2017.11.196. doi: 10.1016/j.procs.2017.11.196
    [10] Y. Sun, G. Qi, Z. Wang, B. J. Y. Wyk, Y. Haman, Chaotic particle swarm optimization, in Proceedings of the 11th Annual Conference Companion on Genetic and Evolutionary Computation Conference, (2009), 12–14.
    [11] S. Behnia, A. Akhshani, H. Mahmodi, A. Akhavan, A novel algorithm for image encryption based on mixture of chaotic maps, Chaos Soliton Fract., 35 (2008), 408–419. https://doi.org/10.1016/j.chaos.2006.05.011. doi: 10.1016/j.chaos.2006.05.011
    [12] N. K. Pareek, V. Patidar, K. K. Sud, Image encryption using chaotic logistic map, Img. Vis. Comput., 24 (2006), 926–934. https://doi.org/10.1016/j.imavis.2006.02.021. doi: 10.1016/j.imavis.2006.02.021
    [13] X. Y. Wang, Synchronization of Chaos Systems and their Application in Secure Communication, Science Press, Beijing, 2012.
    [14] K. M. U. Maheswari, R. Kundu, H. Saxena, Pseudo random number generators algorithms and applications, Int. J. Pure Appl. Math., 118 (2018), 331–336.
    [15] L. Wang, H. Cheng, Pseudo-random number generator based on logistic chaotic system, Entropy, 21 (2019), 1–12. https://doi.org/10.3390/e21100960. doi: 10.3390/e21100960
    [16] J. Lü, X. Yu, G. Chen, Chaos synchronization of general complex dynamical networks, Phys. A, 334 (2004), 281–302. https://doi.org/10.1016/j.physa.2003.10.052. doi: 10.1016/j.physa.2003.10.052
    [17] L. Chen, K. Aihara, Chaotic simulated annealing by a neural network model with transient chaos, Neural Networks, 8 (1995), 915–930. https://doi.org/10.1016/0893-6080(95)00033-V. doi: 10.1016/0893-6080(95)00033-V
    [18] S. Kumari, R. Chugh, A novel four-step feedback procedure for rapid control of chaotic behavior of the logistic map and unstable traffic on the road, Chaos: Interdiscip. J. Nonlinear Sci., 30 (2020), 123115. doi: https://doi.org/10.1063/5.0022212. doi: 10.1063/5.0022212
    [19] S. Strogatz, Nonlinear Dynamics and Chaos: With Applications to Physics, Biology, Chemistry and Engineering, AIP Publishing, 2001.
    [20] R. L. Devaney, An Introduction to Chaotic Dynamical Systems, 2nd edition, Westview Press, USA, 2003.
    [21] M. Henon, A two-dimensional mapping with a strange attractor, Math Phys., 50 (1976), 69–77. https://doi.org/10.1007/BF01608556. doi: 10.1007/BF01608556
    [22] K. Kaneko, Theory and Application of Coupled Map Lattices, John Wiley and Sons, New York, 1993.
    [23] S. J. Cang, Z. H. Wang, Z. Q. Chen, H. Y. Jia, Analytical and numerical investigation of a new Lorenz-like chaotic attractor with compound structures, Nonlinear Dyn., 75 (2014), 745–760. https://doi.org/10.1007/s11071-013-1101-7. doi: 10.1007/s11071-013-1101-7
    [24] D. Lambi\'c, A new discrete-space chaotic map based on the multiplication of integer numbers and its application in S-box design, Nonlinear Dyn., 100 (2020), 699–711. https://doi.org/10.1007/s11071-020-05503-y. doi: 10.1007/s11071-020-05503-y
    [25] S. Kumari, R. Chugh, R. Miculescu, On the complex and chaotic dynamics of standard logistic sine square map, An. Univ. "Ovidius" Constanta - Ser.: Mat., 29 (2021), 201–227. https://doi.org/10.2478/auom-2021-0041. doi: 10.2478/auom-2021-0041
    [26] L. Liu, S. Miao, A new simple one-dimensional chaotic map and its application for image encryption, Multimedia Tools Appl., 77 (2018), 21445–21462. https://doi.org/10.1007/s11042-017-5594-9. doi: 10.1007/s11042-017-5594-9
    [27] X. Zhang, Y. Cao, A novel chaotic map and an improved chaos-based image encryption scheme, Sci. World J., 2014 (2014). https://doi.org/10.1155/2014/713541. doi: 10.1155/2014/713541
    [28] P. Verhulst, The law of population growth, Nouvelles Memories de lÀcadè mie Royale des Sciences et Belles-Lettres de Bruxelles, 18 (1845), 14–54. https://doi.org/10.3406/minf.1845.1813. doi: 10.3406/minf.1845.1813
    [29] R. A. Holmgren, A First Course in Discrete Dynamical Systems, Springer-Verlag, New York, 1996. https://doi.org/10.1007/978-1-4419-8732-7.
    [30] Ashish, J. Cao, R. Chugh, Chaotic behavior of logistic map in superior orbit and an improved chaos-based traffic control model, Nonlinear Dyn., 94 (2018), 959–975. https://doi.org/10.1007/s11071-018-4403-y. doi: 10.1007/s11071-018-4403-y
    [31] R. Chugh, M. Rani, Ashish, On the convergence of logistic map in noor orbit, Int. J. Comput. Appl., 43 (2012), 0975–8887. https://doi.org/10.5120/6200-8739. doi: 10.5120/6200-8739
    [32] S. Kumari, R. Chugh, A new experiment with the convergence and stability of logistic map via SP orbit, Int. J. Appl. Eng. Res., 14 (2019), 797–801. https://www.ripublication.com.
    [33] S. Kumari, R. Chugh, A. Nandal, Bifurcation analysis of logistic map using four step feedback procedure, Int. J. Eng. Adv. Technol., 9 (2019), 704–707. https://doi.org/10.35940/ijeat.F9166.109119. doi: 10.35940/ijeat.F9166.109119
    [34] M. Rani, R. Agarwal, A new experimental approach to study the stability of logistic map, Chaos Solitons Fract., 41 (2009), 2062–2066. https://doi.org/10.1016/j.chaos.2008.08.022. doi: 10.1016/j.chaos.2008.08.022
    [35] M. Rani, S. Goel, An experimental approach to study the logistic map in I-superior orbit, Chaos Complexity Lett., 5 (2011), 95–102.
    [36] M. T. Rosenstein, J. J. Collins, C. J. D. Luca, A practical method for calculating largest Lyapunov exponents from small data sets, Phys. D, 65 (1993), 117–134. https://doi.org/10.1016/0167-2789(93)90009-P. doi: 10.1016/0167-2789(93)90009-P
    [37] A. Wolf, J. B. Swift, H. L. Swinney, J. A. Vastano, Determining Lyapunov exponents from a time series, Phys. D, 16 (1985), 285–317. https://doi.org/10.1016/0167-2789(85)90011-9. doi: 10.1016/0167-2789(85)90011-9
    [38] A. Balestrino, A. Caiti, E. Crisostomi, Generalized entropy of curves for the analysis and classification of dynamical systems, Entropy, 11 (2009), 249–270. https://doi.org/10.3390/e11020249. doi: 10.3390/e11020249
    [39] A. Balestrino, A. Caiti, E. Crisostomi, A classification of nonlinear systems: An entropy based approach, Chem. Eng. Trans., 11 (2007), 119–124.
    [40] A. Balestrino, A. Caiti, E. Crisostomi, Entropy of curves for nonlinear systems classification, IFAC Proc. Vol., 40 (2007), 72–77. https://doi.org/10.3182/20070822-3-ZA-2920.00012. doi: 10.3182/20070822-3-ZA-2920.00012
    [41] C. E. Shannon, A mathematical theory of communication, Bell Syst. Tech. J., 27 (1948), 379–423. https://doi.org/10.1002/j.1538-7305.1948.tb01338.x. doi: 10.1002/j.1538-7305.1948.tb01338.x
    [42] S. Behnia, A. Akhshani, S. Ahadpour, H. Mahmodi, A. Akhavan, A fast chaotic encryption scheme based on piecewise nonlinear chaotic maps, Phys. Lett. A, 366 (2007), 391–396. https://doi.org/10.1016/j.physleta.2007.01.081. doi: 10.1016/j.physleta.2007.01.081
    [43] E. Fatemi-Behbahani, K. Ansari-Asl, E. Farshidi, A new approach to analysis and design of chaos-based random number generators using algorithmic converter, Circuits Syst. Signal Process., 35 (2016), 3830–3846. https://doi.org/10.1007/s00034-016-0248-0. doi: 10.1007/s00034-016-0248-0
    [44] L. Liu, S. Miao, H. Hu, Y. Deng, Pseudo-random bit generator based on non-stationary logistic maps, IET Inf. Secur., 10 (2016), 87–94. https://doi.org/10.1049/iet-ifs.2014.0192. doi: 10.1049/iet-ifs.2014.0192
    [45] A. Kanso, N. Smaoui, Logistic chaotic maps for binary numbers generations, Chaos, Solitons Fractals, 40 (2009), 2557–2568. https://doi.org/10.1016/j.chaos.2007.10.049. doi: 10.1016/j.chaos.2007.10.049
    [46] C. Liu, L. Ding, Q. Ding, Research about the characteristics of chaotic systems based on multi-scale entropy, Entropy, 21 (2019), 663. https://doi.org/10.3390/e21070663. doi: 10.3390/e21070663
    [47] L. Wang, H. Cheng, Pseudo-random number generator based on logistic chaotic system, Entropy, 21 (2019), 960. https://doi.org/10.3390/e21100960. doi: 10.3390/e21100960
    [48] J. Fridrich, Image encryption based on chaotic maps, in Proceedings of IEEE International Conference on Systems, Man and Cybernetics(ICSMC'97), 2 (1997), 1105–1110.
    [49] S. Li, Q. Li, W. Li, X. Mou, Y. Cai, Statistical properties of digital piecewise linear chaotic maps and their roles in cryptography and pseudo-random coding, in IMA International Conference on Cryptography and Coding, 2260 (2001), 205–221. https: //doi.org/10.1007/3-540-45325-3_19.
    [50] L. M. Pecora, T. L. Carroll, Synchronization in chaotic systems, Phys. Rev. Lett., 64 (1990), 821–824. https://doi.org/10.1103/PhysRevLett.64.821. doi: 10.1103/PhysRevLett.64.821
  • This article has been cited by:

    1. Serena Dipierro, Luca Lombardini, Partial differential equations from theory to applications: Dedicated to Alberto Farina, on the occasion of his 50th birthday, 2023, 5, 2640-3501, 1, 10.3934/mine.2023050
  • Reader Comments
  • © 2022 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(3531) PDF downloads(165) Cited by(15)

Figures and Tables

Figures(12)  /  Tables(2)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog