In this paper, a novel one dimensional chaotic map K(x)=μx(1−x)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
[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(1−x)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 u∈C2,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 T∈R, 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 p≥pS=n+2n−2 (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)(n−1)2 (see Remark 1.2 below). Indeed, such triviality when n(n+2)(n−1)2≤p<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 n≥5 with nonnegative Ricci tensor. Let u∈C2,1(M×R) be a bounded eternal solution of Eq (1.1) with ut≥0 and 1<p<pS. Then, u≡0.
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 n≥3 with nonnegative Ricci tensor. Let u∈C2,1(M×R) be an eternal solution of Eq (1.1) with
1<p<n(n+2)(n−1)2. |
Then, u≡0.
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)(n−1)2, |
there exists a constant C=C(n,p) such that
∫B(x,ρ)×(t−ρ2,t+ρ2)u2p≤Cρ2−4pp−1μ(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)u2p≤Cρ2−4pp−1μ(B(x,ρ))≤Cρn+2−4pp−1. |
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 u∈C2,1(Ω×(0,T)), Ω⊂Rn, be a positive bounded solution of Eq (1.1) with 0<u≤D and let ρ>0 such that ¯B(¯x,2ρ)×[¯t−4ρ2,¯t+4ρ2]⊆Ω×(0,T). Then, for n≥1 and p>1, there exists a constant C=C(n,p)>0, such that
|∇u|u≤C(1ρ+√pDp−1)(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ρ+√pDp−1)u1−δ |
in Q.
To keep track of the dependencies of the constants, we define
K=C(1ρ+√pDp−1) | (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 w∈C2(U), any nonnegative η∈C2c(U), any real numbers d,m∈R such that d≠m+2, the following inequality holds:
2(n−m)d−(n−1)(m2+d2)4n∫Uηwm−2|∇w|4−n−1n∫Uηwm(Δw)2−2(n−1)m+(n+2)d2n∫Uηwm−1|∇w|2Δw≤m+d2∫Uwm−1|∇w|2⟨∇w,∇η⟩+∫UwmΔw⟨∇w,∇η⟩+12∫Uwm|∇w|2Δη. |
Remark 2.4. The proof is based on the Bochner formula
12Δ|∇f|2=|∇2f|2+⟨∇f,∇Δf⟩≥1n(Δ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 n≥3 and u∈C2,1(Ω×(0,T)) be a solution of Eq (1.1), with 1<p<pS and such that 0<u≤D. We assume that for some ρ>0 we have ¯B(¯x,2ρ)×[¯t−4ρ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)>−2n−2 and C=C(n,p,r,δ,D,ρ0)>0 such that
∫Qζrumu2t+∫Qζrum−2|∇u|4+∫Qζrum+p−1|∇u|2 | (2.3) |
≤C∫Qζr−4um+2(1−δ)(|∇ζ|2+|ζt|+(Δζ)2+|∇ζ|4+ζ2t)+C∫Qζrum+1−2δ|u−t| | (2.4) |
where u−t=ut∨0. 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 d≠m+2, we get
2(n−m)d−(n−1)(m2+d2)4n∫Uζrum−2|∇u|4−n−1n∫Uζrum(Δu)2−2(n−1)m+(n+2)d2n∫Uζrum−1|∇u|2Δu≤m+d2∫Uum−1|∇u|2⟨∇u,∇(ζr)⟩+∫UumΔu⟨∇u,∇(ζr)⟩+12∫Uum|∇u|2Δ(ζr). |
Hence, substituting Δu=ut−up, we obtain
2(n−m)d−(n−1)(m2+d2)4n∫Uζrum−2|∇u|4−n−1n∫UζrumutΔu+n−1n∫Uζrum+pΔu−2(n−1)m+(n+2)d2n∫Uζrum−1|∇u|2ut+2(n−1)m+(n+2)d2n∫Uζrum+p−1|∇u|2≤m+d2∫Uum−1|∇u|2⟨∇u,∇(ζr)⟩+∫Uumut⟨∇u,∇(ζr)⟩−∫Uum+p⟨∇u,∇(ζr)⟩+12∫Uum|∇u|2Δ(ζr). | (2.5) |
Integrating by parts in the two integrals in the second line above, we get
−n−1n∫UζrumutΔu+n−1n∫Uζrum+pΔu=n−1n∫Uζrum⟨∇u,∇ut⟩+m(n−1)n∫Uζrum−1|∇u|2ut+n−1n∫Uumut⟨∇u,∇(ζr)⟩−n−1n∫Uum+p⟨∇u,∇(ζr)⟩−(m+p)(n−1)n∫Uζrum+p−1|∇u|2. |
Thus, substituting and setting
X=∫Uζrum⟨∇u,∇ut⟩,Y=∫Uζrum−1|∇u|2ut, |
Z=∫Uumut⟨∇u,∇(ζr)⟩,V=∫Uum+p⟨∇u,∇(ζr)⟩, |
W=∫Uum−1|∇u|2⟨∇u,∇(ζr)⟩,R=∫Uum|∇u|2Δ(ζr), |
the inequality (2.5) becomes
2(n−m)d−(n−1)(m2+d2)4n∫Uζrum−2|∇u|4+n−1nX+m(n−1)nY+n−1nZ−n−1nV−(m+p)(n−1)n∫Uζrum+p−1|∇u|2−2(n−1)m+(n+2)d2nY+2(n−1)m+(n+2)d2n∫Uζrum+p−1|∇u|2≤m+d2W+Z−V+12R, |
Hence, rearranging and simplifying, we conclude
α∫Uζrum−2|∇u|4+β∫Uζrum+p−1|∇u|2≤−n−1nX+(n+2)d2nY+1nZ−1nV+m+d2W+12R, | (2.6) |
where
α=2(n−m)d−(n−1)(m2+d2)4nandβ=(n+2)d−2(n−1)p2n. |
We now choose
m=−dn−1 |
obtaining
α=[2(n−1)−(n−2)d]d4(n−1)andβ=(n+2)d−2(n−1)p2n. |
It is then easy to see that both constants α and β are positive if d satisfies
2(n−1)pn+2<d<2(n−1)n−2 | (2.7) |
which is a meaningful condition, since p<pS=n+2n−2. Thus, we set d to be equal to the average mean of the two values, that is
d=(n−1)pn+2+n−1n−2 |
and consequently
m=m(n,p)=−pn+2−1n−2, | (2.8) |
which implies d≠m+2. Indeed, d=m+2 if and only if p=(n+2)(n−4)n(n−2) but this value is always smaller than 1, for n∈N. Moreover, from the right inequality (2.7) we have
m=−dn−1>−2n−2. |
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+p−1|∇u|2+C(ε)∫Uζr−2um+p+1|∇ζ|2, | (2.9) |
W≤ε∫Uζrum−2|∇u|4+C(ε)∫Uζr−4um+2|∇ζ|4, | (2.10) |
R=r(r−1)∫Uζr−2um|∇u|2|∇ζ|2+r∫Uζr−1um|∇u|2Δζ≤ε∫Uζrum−2|∇u|4+C(ε)∫Uζr−4um+2|∇ζ|4+C(ε)∫Uζr−2um+2(Δζ)2≤ε∫Uζrum−2|∇u|4+C(ε)∫Uζr−4um+2|∇ζ|4+C(ε)∫Uζr−4um+2(Δζ)2=ε∫Uζrum−2|∇u|4+C(ε)∫Uζr−4um+2(|∇ζ|4+(Δζ)2) | (2.11) |
since ζ≤1 everywhere. Moreover,
X=∫Uζrum⟨∇u,∇ut⟩=dfdt−P−m2Y | (2.12) |
with
f=12∫Uζrum|∇u|2andP=12∫Uum|∇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=dgdt−S−Z−X−mY=d(g−f)dt+P−Z−S−m2Y, | (2.13) |
where
g=1m+p+1∫Uζrum+p+1andS=1m+p+1∫Uum+p+1(ζr)t, |
which implies
Z=d(g−f)dt+P−S−m2Y−∫Uζrumu2t. | (2.14) |
We bound P as follows,
P=r2∫Uζr−1ζtum|∇u|2≤ε∫Uζrum−2|∇u|4+C(ε,r)∫Uζr−2um+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ζrum−2|∇u|4+β∫Uζrum+p−1|∇u|2≤−n−1n(dfdt−P−m2Y)+(n+2)d2nY+1n(d(g−f)dt+P−S−m2Y−∫Uζrumu2t)+εn∫Uζrum+p−1|∇u|2+(m+d+1)ε2∫Uζrum−2|∇u|4+C∫Uζr−4um+2((Δζ)2+|∇ζ|4)+C∫Uζr−2um+p+1|∇ζ|2≤−dfdt+1ndgdt+((n−1)m2n+(n+2)d2n−m2n)Y−1n∫Uζrumu2t+εn∫Uζrum+p−1|∇u|2+(m+d+3)ε2∫Uζrum−2|∇u|4+C∫Uζr−4um+2((Δζ)2+|∇ζ|4+ζ2t)+C∫Uζr−2um+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
(n−1)m2n+(n+2)d2n−m2n=nd2(n−1)>0 |
as m=−dn−1.
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
α2∫Uζrum−2|∇u|4+β2∫Uζrum+p−1|∇u|2≤−dfdt+1ndgdt+nd2(n−1)Y−1n∫Uζrumu2t.+C1∫Uζr−4um+2((Δζ)2+|∇ζ|4+ζ2t)+C1∫Uζr−2um+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
α2∫Qζrum−2|∇u|4+β2∫Qζrum+p−1|∇u|2≤C1∫Qζr−4um+2((Δζ)2+|∇ζ|4+ζ2t)+C1∫Qζr−2um+p+1(|∇ζ|2+|ζt|)+nd2(n−1)∫t2t1Y−1n∫Qζ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ζrum−1|∇u|2ut≤∫Uζrum−1|∇u|2u+t≤K2∫Uζrum+1−2δu+t, |
where K is defined in formula (2.1) and we set ut=u+t+u−t, u+t=ut∧0≥0 and u−t=ut∨0≤0. It follows
∫t2t1Y≤K2∫U∫t2t1ζrum+1−2δu+t=K2∫U∫t2t1ζrum+1−2δut−K2∫U∫t2t1ζrum+1−2δu−t=K2∫U∫t2t1(ζrum+2(1−δ))t−K2∫U∫t2t1um+2(1−δ)(ζr)t−K2∫U∫t2t1ζrum+1−2δu−t≤rK2∫Qζr−1um+2(1−δ)|ζt|−K2∫Qζrum+1−2δu−t | (2.17) |
since ζ(x,ti)=0, for i=1,2 and every x∈U.
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
1n∫Qζrumu2t+α2∫Qζrum−2|∇u|4+β2∫Qζrum+p−1|∇u|2≤C2(∫Qζr−2um+p+1(|∇ζ|2+|ζt|)+∫Qζr−4um+2((Δζ)2+|∇ζ|4+ζ2t))+C2∫Qζr−1um+2(1−δ)|ζt|−C2∫Qζrum+1−2δu−t. |
Since u≤D, 0≤ζ≤1 and m+p+1>m+2>m+2(1−δ)>0, for δ sufficiently small, as m>−2n−2 and n≥3, we have um+p+1≤Cum+2(1−δ), then there exists a positive constant C3=C3(n,p,r,δ,D,ρ0) such that
∫Qζrumu2t+∫Qζrum−2|∇u|4+∫Qζrum+p−1|∇u|2≤C3∫Qζr−4um+2(1−δ)(|∇ζ|2+|ζt|+(Δζ)2+|∇ζ|4+ζ2t)−C3∫Qζrum+1−2δu−t, |
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(n−1)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)p−1+δ there exists a constant C=C(n,p,r,δ,D,ρ0)>0 such that
∫Qζru2p+m≤C∫Qζr−2(2+m)p−1+δ(|∇ζ|2+|ζt|+(Δζ)2+|∇ζ|4+ζ2t)2p+m2p−2(1−δ)+C∫Qζrum+1−2δ|u−t|. | (2.18) |
Proof. Multiplying Eq (1.1) by ζrum+p and integrating on U, we obtain
∫Uζru2p+m=∫Uζrum+put−∫Uζrum+pΔu=dgdt−S+V+(m+p)∫Uζrum+p−1|∇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+m≤C∫Qζr−4um+2(1−δ)(|∇ζ|2+|ζt|+(Δζ)2+|∇ζ|4+ζ2t)+C∫Qζrum+1−2δ|u−t| |
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 ut≥0 everywhere.
Proposition 2.8. Let Ω⊆Rn, n≥3, and u∈C2,1(Ω×(0,T)) be a positive solution of Eq (1.1) with p<pS=n+2n−2, such that 0<u≤D and ut≥0. Let 0<¯t<T, ¯x∈Ω, m>−2n−2 as in (2.8), δ∈(0,1/100), r>2(2p+m)p−1+δ 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ρ)×[¯t−4ρ2,¯t+4ρ2]⊆Ω×(0,T), letting Q(ρ)=B(¯x,ρ)×(¯t−ρ2,¯t+ρ2), there holds
∫Q(ρ/2)u2p+m≤Cρn+2−2p+mp−1+δ. |
Proof. Since u−t=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 u∈C2,1(Ω×(0,T)) be a positive bounded solution of Eq (1.1) with p<pS=n+2n−2, such that 0<u≤D and ut≥0, where Ω is an open subset of a Riemannian manifold (M,g) with nonnegative Ricci tensor of dimension n≥3. Let 0<¯t<T, ¯x∈Ω, m>−2n−2 as in (2.8), δ∈(0,1/100), r>2(2p+m)p−1+δ 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ρ)×[¯t−4ρ2,¯t+4ρ2]⊆Ω×(0,T), letting Q(ρ)=B(¯x,ρ)×(¯t−ρ2,¯t+ρ2), there holds
∫Q(ρ/2)u2p+m≤Cρ2−2p+mp−1+δμ(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, u−t 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+1−2δ>−2n−2+1−2δ>0, |
when n>4, hence, using Young's inequality, we get
C∫Qζrum+1−2δ|u−t|≤12∫Qζru2p+m+¯C∫Qζr|u−t|2p+m2p−1+2δ |
where C is the constant in formula (2.18). Then, we obtain
12∫Qζru2p+m≤C∫Qζr−2(2+m)p−1+δ(|∇ζ|2+|ζt|+(Δζ)2+|∇ζ|4+ζ2t)2p+m2p−2(1−δ)+¯C∫Q|u−t|2p+m2p−1+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ρ)×[¯t−4ρ2,¯t+4ρ2]⊆Ω×(0,T), there holds
∫Q|u−t|2p+m2p−1+2δ≤˜Cρ2−2p+mp−1+δμ(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 u∈C2,1(Ω×(0,T)) be a positive bounded solution of Eq (1.1), with p<pS=n+2n−2. Following their notation, we set
d=up−1hence, 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 u∈C2,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+22‖u‖2,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)ρ1−n+22q‖d‖1/2q,3ρ, |
hence, in the same hypotheses of this theorem, we can conclude
|u(x,t)|≤Cρ1−n+22q−n+22‖d‖1/2q,3ρ‖u‖2,3ρ≤Cρ|˜Q(3ρ)|−q+12q‖d‖1/2q,3ρ‖u‖2,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+12q‖d‖1/2q,3ρ‖u‖2,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=up−1 in this estimate we get the following corollary.
Corollary 3.2. Let (M,g) be a complete Riemannian manifold with nonnegative Ricci tensor. Let u∈C2,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+12q‖u‖p−12q(p−1),3ρ‖u‖2,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 u∈C2,1(M×R) be a bounded eternal solution of equation (1.1) such that ut≥0, with 1<p<pS=n+2n−2 and n≥3. 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)(n−1)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)(n−1)2≤p<pS. |
We choose m=−dn−1>−2n−2, as in (2.8) and we observe that
2p+m≥2n(n+2)(n−1)2−2n−2=2n3−n2−2n−1(n−2)(n−1)2=2+2(3n2−7n+1)(n−2)(n−1)2>2 |
when n≥3.
Now let γ>0 such that
q=(1+γ)2p+mp−1>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+(1−p−12(1+γ)(p−1+δ)−1p−1+δ)=1+2γ2(1+γ)−1p−1<0, |
hence, the quantity
σ=1−p−12(1+γ)(p−1+δ)−1p−1+δ | (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ρ)×[¯t−9ρ2,¯t+9ρ2]), we obtain
‖u‖2,3ρ=(∫˜Q(3ρ)u2)1/2≤(ρ2μ(B(¯x,3ρ)))2p+m−22(2p+m)(∫˜Q(3ρ)u2p+m)12p+m≤C(ρ2μ(B(¯x,3ρ)))1/2ρ−1p−1+δ | (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
‖u‖p−1q(p−1),3ρ=(∫˜Q(3ρ)uq(p−1))1/q=(∫˜Q(3ρ)u(1+γ)(2p+m))1/q≤C(∫˜Q(3ρ)u2p+m)1/q≤C(ρ2μ(B(¯x,3ρ)))1/qρ−p−1(1+γ)(p−1+δ). |
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+12q‖u‖p−12q(p−1),3ρ‖u‖2,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+12q‖u‖p−12q(p−1),3ρ‖u‖2,3ρ≤Cρ−1qμ(B(¯x,3ρ))−q+12q(ρ2μ(B(¯x,3ρ)))12qρ−p−12(1+γ)(p−1+δ)(ρ2μ(B(¯x,3ρ)))12ρ−1p−1+δ=Cρ1−p−12(1+γ)(p−1+δ)−1p−1+δ=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 ¯t∈R arbitrary, u=0 everywhere.
Finally, one can check that, if n≥5 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+m≥2 and p<n+2n−2, for any n≥3 we can choose a γ>0 such that
γ>n+24(n+2n−2−1)−1=4n−2. |
On the other hand, to fulfill condition (3.2) it is sufficient to choose γ that satisfies
11+2γ>p−2, |
which, for p<n+2n−2 and n≥6, 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
![]() |
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 |