The study of harvesting mechanisms in predator-prey systems with an Allee effect on prey has always garnered significant attention. In this paper, the dynamics and harvesting strategies of a predator-prey system are investigated, where the prey is subject to the Allee effect. The positivity and boundedness of solutions, the existence and stability of equilibria are further studied. The existence of a Hopf bifurcation at the interior equilibrium point of the system is investigated and verified by numerical simulations. Furthermore, we investigate the maximum sustainable yield (MSY), maximum sustainable total yield (MSTY) and the optimal economic profit of the proposed system. We find that MSY can be attained through predator harvesting, while MSTY is observed when harvesting efforts are uniform across species. In these situations, the biological system maintains stability. Using the method of control parametrization, the optimal economic profit and harvesting strategy are obtained. The results show that the harvesting efforts can affect the stability of the system, resulting in several interesting biological phenomena. This research provides a theoretical basis for biological resource management.
Citation: Chengchong Lu, Xinxin Liu, Zhicheng Li. The dynamics and harvesting strategies of a predator-prey system with Allee effect on prey[J]. AIMS Mathematics, 2023, 8(12): 28897-28925. doi: 10.3934/math.20231481
[1] | Xiangqi Zheng . On the extinction of continuous-state branching processes in random environments. AIMS Mathematics, 2021, 6(1): 156-167. doi: 10.3934/math.2021011 |
[2] | Tawatchai Petaratip, Piyapong Niamsup . Stability analysis of an unemployment model with time delay. AIMS Mathematics, 2021, 6(7): 7421-7440. doi: 10.3934/math.2021434 |
[3] | Kexin Zhang, Caihui Yu, Hongbin Wang, Xianghong Li . Multi-scale dynamics of predator-prey systems with Holling-IV functional response. AIMS Mathematics, 2024, 9(2): 3559-3575. doi: 10.3934/math.2024174 |
[4] | Huda Abdul Satar, Raid Kamel Naji, Mainul Haque . A food chain model with Allee effect: Analysis on the behaviors of equilibria. AIMS Mathematics, 2025, 10(5): 12598-12618. doi: 10.3934/math.2025568 |
[5] | Yuanfu Shao . Dynamics and optimal harvesting of a stochastic predator-prey system with regime switching, S-type distributed time delays and Lévy jumps. AIMS Mathematics, 2022, 7(3): 4068-4093. doi: 10.3934/math.2022225 |
[6] | Shufan Wang, Zhihui Ma, Xiaohua Li, Ting Qi . A generalized delay-induced SIRS epidemic model with relapse. AIMS Mathematics, 2022, 7(4): 6600-6618. doi: 10.3934/math.2022368 |
[7] | Bingxian Wang, Mei Xu . Asymptotic behavior of some differential inequalities with mixed delays on time scales and their applications. AIMS Mathematics, 2024, 9(6): 16453-16467. doi: 10.3934/math.2024797 |
[8] | Huizhen Qu, Jianwen Zhou . S-asymptotically ω-periodic dynamics in a fractional-order dual inertial neural networks with time-varying lags. AIMS Mathematics, 2022, 7(2): 2782-2809. doi: 10.3934/math.2022154 |
[9] | Linxia Hu, Yonghong Shen, Xiumei Jia . Global behavior of a discrete population model. AIMS Mathematics, 2024, 9(5): 12128-12143. doi: 10.3934/math.2024592 |
[10] | Adel M. Al-Mahdi, Maher Noor, Mohammed M. Al-Gharabli, Baowei Feng, Abdelaziz Soufyane . Stability analysis for a Rao-Nakra sandwich beam equation with time-varying weights and frictional dampings. AIMS Mathematics, 2024, 9(5): 12570-12587. doi: 10.3934/math.2024615 |
The study of harvesting mechanisms in predator-prey systems with an Allee effect on prey has always garnered significant attention. In this paper, the dynamics and harvesting strategies of a predator-prey system are investigated, where the prey is subject to the Allee effect. The positivity and boundedness of solutions, the existence and stability of equilibria are further studied. The existence of a Hopf bifurcation at the interior equilibrium point of the system is investigated and verified by numerical simulations. Furthermore, we investigate the maximum sustainable yield (MSY), maximum sustainable total yield (MSTY) and the optimal economic profit of the proposed system. We find that MSY can be attained through predator harvesting, while MSTY is observed when harvesting efforts are uniform across species. In these situations, the biological system maintains stability. Using the method of control parametrization, the optimal economic profit and harvesting strategy are obtained. The results show that the harvesting efforts can affect the stability of the system, resulting in several interesting biological phenomena. This research provides a theoretical basis for biological resource management.
A basic but important problem in number theory is the triple correlation sums problem, which concerns the non-trivial bounds for
∑h≤H∑n≤Xa(n)b(l1n+l2h)c(h)or∑h≤H∑n≤Xa(n)b(n+l1h)c(n+l2h). |
Here, a(n),b(n) and c(n) are three arithmetic functions, H,X≥2 and l1,l2∈Z. These type of sums play the vital roles of their own in many topics, such as the moments of L-functions (or zeta-functions), subconvexity, the Gauss circle problem and the Quantum Unique Ergodicity (QUE) conjecture, etc (see for instance [1,5,7,8,9,11,12,13,15,19] and the references therein). In the case of all the arithmetic functions being the divisor functions, Browning [4] found that
∑h≤H,n≤Xa(n)b(n+h)c(n+2h)=118Υ(h)∏p(1−1p)2(1+2p)HXlog3X+o(HXlog3X) |
for certain function Υ(h), provided that H≥X3/4+ε. It is remarkable that Blomer [2] used the spectral decomposition for partially smoothed triple correlation sums to prove that
∑h≥1∑n≤XW(hH)τ(n)a(n+h)τ(n+2h)=HˆW(1)∑n≤Xa(n)∑d≥1S(2n,0;d)d2×(logn+2γ−2logd)2+O((H2√X+HX14+√XH+X√H)‖a‖2). |
Here, the gamma constant γ≈0.57721, W is a bump function supported on [1/2,5/2], satisfying that xjW(j)(x)≪1 for any j∈N+; while, ˆW is the Mellin transform of W and
‖a‖2=√∑n≤X|a2(n)| |
is the ℓ2-norm. Notice that, the parameter H is reduced to H≥X1/3+ε. Let k,k′∈2N+. Let f1∈B∗k(1) and f2∈B∗k′(1) be two Hecke newforms on GL2, with λf1 and λf2 being their n-th Hecke eigenvalues, respectively (see §2 for definitions). Subsequently, Lin [20] claimed that
∑h≥1∑n≤XW(hH)λf1(n)a(n+h)λf2(n+2h)≪Xε(√XH+X√H)‖a‖2, |
which is non-trivial, provided that H≥X2/3+ε. Recently, Singh [28], however, were able to attain that, for any f1,f2,f3∈B∗k(1) (or B∗λ(1)) and some constant η>0,
∑h≥1∑n≥1W1(hH)W2(nX)λf1(n)λf2(n+h)λf3(n+2h)≪X1−η+εH, |
where W1,W2 are any bump functions compactly supported on [1/2,5/2] with bounded derivatives. Here, for any N∈N+, B∗λ(N) denotes the collection of the primitive newforms of Laplacian eigenvalue λ on Γ0(N) (see §2 for backgrounds). Until now, the best result is due to Lü-Xi [21,22] who achieved that
∑h≥1∑n≤XW(hH)a(n)b(n+h)λf1(n+2h)≪XεΔ1(X,H)‖a‖2‖b‖2, |
which allows one to take H≥X2/5+ε. Here, the definition of Δ1(X,H) can be referred to [22,Theorem 3.1]. More recently, Hulse et al. [14] successfully attained
∑h≥1∑n≥1λg1(n)λg2(h)λg3(2n−h)exp(−hH−nX)≪Xκ−1+ϑ+12+εHκ−12−ϑ+12+ε, | (1.1) |
where ϑ<7/64 denotes the currently best record for the Generalized Ramanujan Conjecture. Here, λg1(n) (resp. λg1(n) and λg1(n)) denote the n-th non-normalized coefficients of holomorphic cusp forms g1 (resp. g2 and g3), each of weight κ≥2 and level M≥2.
In the present paper, we aim to consider the level aspect for the triple correlation sums. It is noticeable that, just lately, Munshi [26] obtained that, whenever X1/3+ε≤p≤X, one has the inequality
∑n≥1λf(n)λf(n+pm)≪p14X34+ε | (1.2) |
for any f∈B∗k(p) and fixed integer m such that |m|≤X/p. In this paper, we would like to go further obtaining the following quantitive estimate.
Theorem 1.1. Let X,H≥2, and p be a prime such that p≤X. Let U,V be two smooth weight functions supported [1/2,5/2] with bounded derivatives. Then, for any sequence a={an∈C} and any newforms g∈B∗k(1) (or B∗λ(1)) and f∈B∗k(p) (or B∗λ(p)), we have
∑m≥1∑n≥1anλg(m)λf(m+pn)U(mX)V(nH)≪Xεmax(√XHp,X)‖a‖2, | (1.3) |
where the implied constant depends only on the weight k (or the spectral parameter λ) and ε.
Remark 1.2. Our main result (1.3) is non-trivial for any given parameters X and H; particularly, for any automorphic cusp form π of any rank N, N≥2, with λπ(n) being its n-th normoailzed Fourier coefficient, we find
∑m≥1∑n≥1λπ(n)λg(m)λf(m+pn)U(mX)V(nH)≪Xεmax(H√Xp,X√H) |
by the Rankin-Selberg's bound, which says that ∑n≤X|λπ(n)|2≪π,εX1+ε (see for instance [6,Remark 12.1.8]).
Remark 1.3. The merits that comes from (1.3) is that the implied constant does not depend on the level parameter anymore. One may verify that our result (1.3), however, is exhibited to be a strengthened upper-bound whenever H≥√X/p, compared with an application of Munshi's estimate (1.2). Indeed, Munshi's estimate implies the upper-bound ≪p1/4X3/4+ε√H‖a‖2 for the triple sum above. Moreover, one might save roughly a magnitude of √H in the interesting case of pH≍X; in the average sense, the main result improves upon the estimate due to Munshi. One, on the other hand, wanders whether or not the non-trivial bounds for the scenarios where the cusp forms f,g being of higher rank could be achieved; we shall plan to investigates this problem in the future work.
Notations. Throughout the paper, ε always denotes an arbitrarily small positive constant. n∼X means that X<n≤2X. For any integers m,n, (m,n) means the great common divisor of m,n. Finally, μ(n) denotes the Möbius function of n.
For any k∈2N+ and N∈N+, let us denote by Sk(N) the vector space of the normalized holomorphic cusp forms on Γ0(N) of weight k and trivial nebentypus. Whenever f∈Sk(N), one has
f(z)=∑n≥1λf(n)nk−12e(nz) |
for Im(z)>0. We also denote by Sλ(N) the vector space of the normalized Maaß forms on Γ0(N) of weight 0, Laplacian eigenvalue λ=1/4+r2 (r∈R) and trivial nebentypus. For any f∈Sλ(N), there exists the following Fourier expansion
f(z)=2√|y|∑n≠0λf(n)Kir(2π|ny|)e(nx), |
where z=x+iy. The set of the primitive forms B∗k(N) (or B∗λ(N)) consists of common eigenfunctions of all the Hecke operators Tn for any n≥1. Regarding the individual bounds for λf(n), we have
λf(n)≪(nN)ε, | (2.1) |
whenever f∈B∗k(N) (or B∗λ(N)).
We will need the Voronoi summation formula in the analysis; see [18,Theorem A.4].
Lemma 2.1. Let k,N and the form f be as before. For any a,q∈N+ such that (a,q)=1, we set N2:=N/(N,q). Let h be a bump function of bounded derivatives. Then, there exists a constant ϱ of modulus one and a newform f⋆∈B∗k(N) (or B∗λ(N)) such that
∑n≥1λf(n)e(anq)h(nX)=2πϱq√N2∑n≥1λf⋆(n)e(−¯aN2nq)H♭(nXq2N2;h)+2πϱq√N2∑n≥1λf⋆(n)e(¯aN2nq)H♮(nXq2N2;h), |
where
H♭(x;h)=∫∞0h(ξ)Jf(4π√xξ)dξ,andH♮(x;h)=∫∞0h(ξ)Kf(4π√xξ)dξ. |
Here, if f is holomorphic
Jf(x)=2πikJk−1(x),Kf(x)=0; |
while, if f is a Maaß form
Jf(x)=−πsin(πir)(J2ir(x)−J−2ir(x)),Kf(x)=4cosh(πr)K2ir(x). |
For any s∈R, one may write
Jk−1(s)=s−12(F+k(s)e(x)+F−k(s)e(−s)) | (2.2) |
for some smooth functions F± satisfying that
sjF±k(j)(s)≪k,js(1+s)32 |
for any j≥0; the resource might be referred to [29,Section 6.5] if s<1 and [29,Section 3.4] if s≥1. One thus sees that, for H♭, the n-variable is essentially truncated at n≪q2N2/X1−ε, by repeated integration by parts. Furthermore, notice that, by Appendix of [3],
K2ir(s)≪r,ε{s−12exp(−s),s>1+π|r|,(1+|r|)ε,0<s≤1+π|r|; | (2.3) |
one will find the n-variable enjoys the analogous truncation range for H♮ with that for H♭.
Now, let us recall the following Wilton-type bound; the resource, however, may be referred to [17], together with [10,27].
Lemma 2.2. Let X≥2 and W be a smooth function, compactly supported on [1/2,5/2] such that xjW(j)(x)≪1 for any j∈N+. For any newform f∈B∗k(N) (or B∗λ(N)) and α∈R, we thus have
∑n≥1λf(n)e(nα)W(nX)≪√XN13+ε, | (2.4) |
where the implied ≪-constant depends merely on k (or λ) and ε.
As a variant of the circle method, the δ-symbol method plays a focal role in number theory. We will now briefly recall a version of the circle method; see for instance [16,Chapter 20].
Lemma 2.3. Fix X,Q≥1. For any n≤X, one has
δ(n)=1Q∑q≤Q1q∑amodq(a,q)=1e(anq)∫Rg(q,τ)e(nτqQ)dτ, |
where
g(q,τ)=1+h(q,τ)withh(q,τ)=O(1qQ(|τ|+qQ))A,τj∂j∂τjg(q,τ)≪logQmin(Qq,1|τ|), |
and g(q,τ)≪|τ|−A for any sufficiently large A. In particular, the effective range ofthe τ-integral is [−Xε,Xε].
We will have a need of the following lemmas which will be applied in §3.
Lemma 2.4. Let Q≥2. Let F(x,y) be a smooth bump functionsupported on [1/2,5/2]×[1/2,5/2], which satisfies that
XiYj∂i∂xi∂j∂yjF(xX,yY)≪i,j1 |
for any integers i,j∈N+ and any X,Y≥1. Then, for any c∈Z, sequence a={an∈C} and newform f∈B∗k(1) (or B∗λ(1)), there holds that
∑q≥1∑n≥1anS(n,c;q)F(nX,qQ)≪Xε(√XQ+1Q2>XQ2)‖a‖2, | (2.5) |
where the symbol 1P equals 1 if the assertion P is true, and 0 otherwise.
Proof. First, via the Cauchy-Schwarz inequality, we might evaluate the double sum as
≤(∑q1,q2≥1∑n≥1S(n,c;q1)¯S(n,c;q1)F(nX,q1Q)F(nX,q2Q))12‖a‖2. |
It thus follows from the Weil bound that the non-generic terms q1=q2 shall contribute a upper-bound ≪X1+εQ2 to the parentheses above, which gives the term √XQ‖a‖2 on the RHS of (2.5); while, for the generic terms q1≠q2, if one writes q1=ˆq1δ, q2=ˆq2δ with δ=(q1,q2) satisfying that (δ,ˆq1)=1, Poisson summation formula with the modulus ˆq1ˆq2δ thus might produce the following bound for the triple sum that
≪∑δ≤Q∑ˆq1,ˆq2≤Q/δsup0<|l|≪ˆq1ˆq2δ/X1−ε|∑αmodˆq1ˆq2δS(α,c;ˆq1δ)¯S(α,c;ˆq2δ)e(αlˆq1ˆq2δ)|. |
Notice, here, the inner-most sum vanishes, if l=0, and it is necessary that Q2>X as well. At this point, on applying Chinese remainder theorem, the sum over α turns out to be
ˆq1ˆq2δe(−c¯δl⋅ˆq2ˆq1)∑∑∗smodδe(¯aˆq1sδ−ˆq1c⋅¯(ˆq2s+l)ˆq2δ) |
with ¯δδ≡1modˆq1, ¯ll≡1modˆq1, ¯ˆq1ˆq1≡1modδ and ¯ˆq2s+l(ˆq2s+l)≡1modˆq2δ; trivially evaluating everything thus exactly leads to the term Q2‖a‖2 in (2.5).
Lemma 2.5. Let the parameters X,Q,c, the form f and the sequence a={an∈C} be as in Lemma 2.4. Let W(x,y,z) be a smooth bump functionsupported on [1/2,5/2]×[1/2,5/2]×[1/2,5/2], with the partial derivatives satisfying
XiYjZk∂i∂xi∂j∂yj∂k∂zkW(xX,yY,zZ)≪i,j,k1 |
for ever integers i,j,k∈N+ and any X,Y,Z≥1. There thus holds that
∑q≥1∑n≥1an∑m≥1λf(m)√mS(m−np,c;q)W(mX,nH,qQ)≪Xε(√HQ+Q2)‖a‖2. | (2.6) |
Proof. To show the lemma, the initial procedure is to invoke the Cauchy-Schwarz inequality; we are thus led to evaluating
∑q1,q2≥1∑n≥1∑m1,m2≥1λf(m1)¯λf(m2)√m1m2S(m1−np,c;q1)¯S(m2−np,c;q2)W(m1X,nH,q1Q)W(m2X,nH,q2Q). | (2.7) |
(1) First, let us begin with considering the generic terms q1=q2=q, say. In this moment, Poisson is applicable, which yields an alternative form for (2.7):
H∑l1∈Z∑q∼Q1q∑m1,m2≥1λf(m1)¯λf(m2)√m1m2Y0(l1,m1,m2,p,c;q)I0(l1,m1,m2), |
where
Y0(l,m1,m2,p,c;q)=∑αmodqS(m1−αp,c;q)¯S(m2−αp,c;q)e(−αlq), |
and
I0(l,m1,m2)=∫RW(m1X,ξ,qQ)W(m2X,ξ,qQ)e(lHξq)dξ. |
Notice that the exponential sum modulo q asymptotically equals
qe(−m1¯plq)∑∑∗αmodqe((m1−m2)α+n⋅¯α−¯pl−n¯αq), |
where we have employed the relation involving Ramanujan sum that
S(n,0;q)=∑ab=qμ(a)∑βmodqe(βnb). | (2.8) |
Upon combining with Lemma 2.2, one thus sees that the zero-frequency shall contribute a bound by ≪HQXε for any ε>0. While, on the other hand, if Q>H, one may find the non-zero frequencies will be indispensable to contribute a magnitude to (2.7). It can be demonstrated that, in this situation, the contribution, however, is estimated as ≪Q2Xε; this gives totally a quantity by
≪HQXε+1Q>HQ2Xε. | (2.9) |
(2) Now, we are left with the non-generic case where q1≠q2. One writes q1=q′1h, q2=q′2h, with (q1,q2)=h and (q′1,q′2)=1. Notice that h is co-prime with one of factors q′1,q′2; without loss of generality, one assumes that (h,q′1)=1. The expression in (2.7) thus becomes
∑h≪Q∑q1,q2≤Q/h∑n≥1∑m1,m2≥1λf(m1)¯λf(m2)√m1m2S(m1−np,c;q′1h)ׯS(m2−np,c;q′2h)W(m1X,nH,q′1hQ)W(m2X,nH,q′2hQ). | (2.10) |
Upon exploiting the Poisson twice, we thus arrive at
H∑l2∈Z∑h≪Q∑q′1,q′2≤Q/h1q′1q′2h∑m1,m2≥1λf(m1)¯λf(m2)√m1m2Y†(l2,m1,m2,p,c;q′1,q′2,h)I†(l2,m1,m2,q′1,q′2h), |
where the exponential sum Y† and the resulting integral I† are defined as
Y†(l,m1,m2,p,c;q′1,q′2,h)=∑αmodq′1q′2hS(m1−αp,c;q′1h)¯S(m2−αp,c;q′2h)e(−αlq′1q′2h), |
and
I†(l,m1,m2,q′1,q′2,h)=∫RW(m1X,ξ,q′1hQ)W(m2X,ξ,q′2hQ)e(lHξq′1q′2h)dξ. |
Here, one finds that the zero-frequency l2=0 does exist anymore. Indeed, upon recalling that (h,q′1)=1, one writes α=q′1¯q′1x+q′2h¯q′2y, with xmodq′2h and ymodq′1 such that (x,q′2h) and (y,q′1)=1; applying Chinese remainder theorem and (2.8), the sum over α thus essentially turns out to be
q′1q′2he(m1¯q′2pl+nq′2p¯hlq′1+m2¯q′1plq′2h)∑∑∗αmodhe((m1−m2)¯q′1q′2α−nq′1⋅¯q′2α−¯pl+nq′2¯q′1¯αq′2h). |
Via Lemma 2.2, it thus follows that the display (2.10) is dominated by ≪εXεQ4, upon opening the Kloosterman sum above. This, together with (2.9), shows the desired estimates in the parentheses of (2.6).
In this part, let us focus on the proof of Theorem 1.1. We shall first manage to separate the variables n,m by applying Lemma 2.3; in this paper, we shall employ a vital trick, that is, the 'conductor lowering mechanism' (see [23,24] or the survey [25]). One may see that actually there holds the following
δ(n)=1pQ∑q≤Q1q∑amodqp(a,q)=1e(anqp)∫Rg(q,τ)e(nτqQp)dτ; | (3.1) |
while, the parameter Q shall be taken as Q=√X/p. Now, for three smooth functions U,V,R, supported [1/2,5/2] with bounded derivatives, we shall detect the shift l=m+pn via (3.1), which yields an alternative form for the triple sum in (1.3) as follows
S(X,p,H)=X√HpQ∫R∑q≤Qg(q,τ)q∑γmodpq(γ,q)=1∑l≥1al√le(−γlq)Vτ(lH)∑m≥1λg(m)√m×e(mγpq)Uτ(mX)∑n≥1λf(n)√ne(−nγpq)Rτ(nX)dτ, | (3.2) |
where
Uτ(m)=U(m)e(mXτpqQ),Rτ(n)=R(n)e(−nXτpqQ),Vτ(u)=V(u)e(−uHτqQ). |
We shall proceed to distinguish whether (γ,p)=1 or not in the analysis, so that we are led to three parts, i.e., the non-degenerate term SNon−de., the degenerate term SDeg. and the error term SErr., which are respectively given by
SNon−de.(X,p,H)=X√HpQ∫R∑q≤Q(q,p)=1g(q,τ)q∑∑∗γmodpq∑l≥1al√le(−lγpq)Vτ(lH)×∑m≥1λg(m)√me(mγpq)Uτ(mX)∑n≥1λf(n)√ne(−nγq)Rτ(nX)dτ, | (3.3) |
SDeg.(X,p,H)=X√HpQ∫R∑q≤Q(q,p)=1g(q,τ)q∑∑∗γmodq∑l≥1al√ue(−lpγq)Vτ(lH)×∑m≥1λπ(1,m)√me(mγq)Uτ(mX)∑n≥1λf(n)√ne(−nγq)Rτ(nX)dτ, | (3.4) |
and
SErr.(X,p,H)=X√HpQ∫R∑q≤Qp|qg(q,τ)q∑∑∗γmodpq∑l≥1al√le(−lγpq)Vτ(lH)×∑m≥1λg(m)√me(mγpq)Uτ(mX)∑n≥1λf(n)√ne(−nγq)Rτ(nX)dτ. | (3.5) |
One might see that, here, it suffices to consider SNon−de.; the same argument works for SErr. which serves as a noisy term and for which we save more. We shall now now begin with SNon−de.; the analysis of the term SDeg. will be postponed to the end of this paper.
In this part, let us concentrate on the analysis of SNon−de.. For any ι,ν,υ,ρ∈R, write
Wτ(ι,ν,υ,ρ)=Uτ(ι)Rτ(ν)Vτ(υ)ηQ(ρQ), |
where ηQ is a smooth bump function supported on [1Q/2,5Q/2], satisfying that η(j)Q(x)≪1 for any j∈N+. One might find that the quantity we are focusing on is the following
X√HpQsupτ≪XεsupQ≤Q∑q≥1(q,p)=1g(q,τ)q∑∑∗γmodpq∑u≥1al√le(−lγq)∑m≥1λg(m)√m×e(mγpq)∑n≥1λf(n)√ne(−nγpq)Wτ(mX,nX,lH,qQ). | (3.6) |
We intend to invoking the Voronoi formula, Lemma 2.1; we thus arrive at
X√HpQsupτ≪XεsupQ≤Q∑q≥1(q,p)=1g(q,τ)q∑l≥1al√l∑m≥1λg(m)√m∑n≥1λf(n)√n×S(n−m,pl;pq){^W♭τ(mXpQ2,nXpQ2,lH,qQ)+^W♮τ(mXpQ2,nXpQ2,lH,qQ)}, | (3.7) |
where, for any ⋆,∗∈{♭,♮}, each integral ^W⋆,∗τ is defined as
^W⋆,∗τ(ι,ν,υ,ρ)=ηQ(ρQ)Vτ(υ)H⋆(Q2ιq2;Rτ)H∗(Q2νq2;Uτ). |
It is remarkable that, here, from (2.2) and (2.3), we have the identical crude estimate that ^W⋆,∗τ≪Xε for any ε>0; in this sense, one sees that it suffices to deal simply with ^W♭,♭τ, upon noticing that the argument of the other terms (i.e., ^W♭,♮τ, ^W♮,♭τ and ^W♮,♮τ) can follow similarly with it. One, however, on the other hand, sees that the inner-most sum modulo pq can be converted into
pS(¯p(n−m),l;q) |
with n≡mmodp. Now, if one writes n=m+pk with k≪Xε, we find that (3.6) is no more than
≪X√HQsupτ≪XεsupQ≤Q∑m≪pXε∑k≪Xελf(m)λf(pk+m)√m(pk+m)×∑q≥1(q,p)=1g(q,τ)q∑l≥1al√lS(l,k;q)^W♭,♭τ(mXpQ2,nXpQ2,lH,qQ). |
At this point, an application of Lemma 2.3 shows that the RHS of the expression above is bounded by
≪X1+εQsupτ≪XεsupQ≤Q(√H+1Q2>HQ)‖a‖2≪Xεmax(√XHp,X)‖a‖2, | (3.8) |
upon recalling the value of Q.
Now, let us have a look at the multiple-sum SDeg.. One might verify that SDeg. is of the form
X√HpQsupτ≪XεsupQ≤Q∑q≥1(q,p)=1g(q,τ)q∑∑∗γmodq∑l≥1al√ue(−lpγq)∑m≥1λg(m)√m×e(mγq)∑n≥1λf(n)√ne(−nγq)Wτ(nX,mX,uH,qQ) | (3.9) |
with Wτ being as before. We will now proceed by appealing to the Voronoi formula, Lemma 2.1, again to transform the sums over n into the dualized form, so that we infer that the expression above should be controlled by
X√HpQsupτ≪XεsupQ≤Qsupn≪pQ2X−1+ε|∑q≥1(q,p)=1g(q,τ)q∑l≥1al√l∑m≥1λg(m)√m=×S(m−lp,n;q){~W♭τ(mX,nXpQ2,lH,qQ)+~W♮τ(mX,nXpQ2,lH,qQ)}|, |
where, for ⋆∈{♭,♮}, each integral transform ~W⋆τ is given by
~W⋆τ(ι,ν,υ,ρ)=ηQ(ρQ)Vτ(υ)Uτ(ι)H⋆(Q2νq2;Rτ). |
Via Lemma 2.5, one thus deduces
SDeg.(X,p,H)≪X1+εpQsupτ≪XεsupQ≤Q(√H+Q)‖a‖2≪Xε(√XHp+Xp)‖a‖2. |
This leads to the estimates we would like to prove in Theorem 1.1, upon combining with (3.8).
In this paper, we investigate the triple correlations sums of Fourier coefficient of newforms on GL2, with the levels aspects being explicitly determined; our method is flexible enough to deal with the Maaß new forms. It is also remarkable that more recently, the authors are able to establish a sharp bound in the scenario where one of the froms f,g in (1.3) is a Maaß cuspidal form on GL3 (not necessarily self-dual) with the trivial level.
This research was funded by National Natural Science Foundation of China (grant numbers 413618022 and 61402335).
The authors declare that they have no conflicts of interest.
[1] |
R. J. Han, L. N. Guin, B. X. Dai, Consequences of refuge and diffusion in a spatiotemporal predator-prey model, Nonlinear Anal. Real, 60 (2021), 103311. https://doi.org/10.1016/j.nonrwa.2021.103311 doi: 10.1016/j.nonrwa.2021.103311
![]() |
[2] |
L. N. Guin, P. K. Mandal, Effect of prey refuge on spatiotemporal dynamics of the reaction-diffusion system, Comput. Math. Appl., 68 (2014), 1325–1340. https://doi.org/10.1016/j.camwa.2014.08.025 doi: 10.1016/j.camwa.2014.08.025
![]() |
[3] |
X. X. Liu, Q. D. Huang, Dynamics of a predator-prey system with inducible defense and disease in the prey, Nonlinear Anal. Real, 71 (2023), 103802. https://doi.org/10.1016/j.nonrwa.2022.103802 doi: 10.1016/j.nonrwa.2022.103802
![]() |
[4] |
D. M. Johnson, A. M. Liebhold, P. C. Tobin, O. N. Bjørnstad, Allee effects and pulsed invasion by the gypsy moth, Nature, 444 (2006), 361–363. https://doi.org/10.1038/nature05242 doi: 10.1038/nature05242
![]() |
[5] |
L. Hammann, B. Silliman, B. Blasius, Optimal planting distance in a simple model of habitat restoration with an Allee effect, Front. Mar. Sci., 7 (2021), 610412. https://doi.org/10.3389/fmars.2020.610412 doi: 10.3389/fmars.2020.610412
![]() |
[6] |
R. A. Morris, A. Hernández-Flores, F. Aranceta-Garza, A. Cuevas-Jiménez, K. Kitson-Walters, Allee effect considerations to support the spatial management of a sedentary marine species: The queen conch (Aliger gigas) on the Pedro Bank, Jamaica, Mar. Freshwater Res., 74 (2023), 375–386. https://doi.org/10.1071/MF22140 doi: 10.1071/MF22140
![]() |
[7] |
D. W. Morris, Measuring the Allee effect: positive density dependence in small mammals, Ecology, 83 (2002), 14–20. https://doi.org/10.1890/0012-9658(2002)083[0014:MTAEPD]2.0.CO;2 doi: 10.1890/0012-9658(2002)083[0014:MTAEPD]2.0.CO;2
![]() |
[8] |
G. Mandal, S. Das, L. N. Guin, S. Chakravarty, Dynamic response of a system of interactive species influenced by fear and Allee consequences, Eur. Phys. J. Plus, 138 (2023), 661. https://doi.org/10.1140/epjp/s13360-023-04246-0 doi: 10.1140/epjp/s13360-023-04246-0
![]() |
[9] |
J. Gascoigne, L. Berec, S. Gregory, F. Courchamp, Dangerously few liaisons: A review of mate-finding Allee effects, Popul. Ecol., 51 (2009), 355–372. https://doi.org/10.1007/s10144-009-0146-4 doi: 10.1007/s10144-009-0146-4
![]() |
[10] |
D. Sen, S. Ghorai, S. Sharma, M. Banerjee, Allee effect in prey's growth reduces the dynamical complexity in prey-predator model with generalist predator, Appl. Math. Model., 91 (2021), 768–790. https://doi.org/10.1016/j.apm.2020.09.046 doi: 10.1016/j.apm.2020.09.046
![]() |
[11] |
J. P. Tripathi, P. S. Mandal, A. Poonia, V. P. Bajiya, A widespread interaction between generalist and specialist enemies: The role of intraguild predation and Allee effect, Appl. Math. Model., 89 (2021), 105–135. https://doi.org/10.1016/j.apm.2020.06.074 doi: 10.1016/j.apm.2020.06.074
![]() |
[12] |
D. W. Goodsman, M. A. Lewis, The minimum founding population in dispersing organisms subject to strong Allee effects, Methods Ecol. Evol., 7 (2016), 1100–1109. https://doi.org/10.1111/2041-210X.12573 doi: 10.1111/2041-210X.12573
![]() |
[13] |
C. Zhang, X. Y. Li, Dynamics of a discrete Leslie-Gower model with harvesting and Holling-Ⅱ functional response, Mathematics, 11 (2023), 3303. https://doi.org/10.3390/math11153303 doi: 10.3390/math11153303
![]() |
[14] |
Z. W. Xiao, X. D. Xie, Y. L. Xue, Stability and bifurcation in a Holling type Ⅱ predator-prey model with Allee effect and time delay, Adv. Differ. Equ., 2018 (2018), 288. https://doi.org/10.1186/s13662-018-1742-4 doi: 10.1186/s13662-018-1742-4
![]() |
[15] |
W. X. Wang, Y. B. Zhang, C. Z Liu, Analysis of a discrete-time predator prey system with Allee effect, Ecol. Complex., 8 (2011), 81–85. https://doi.org/10.1016/j.ecocom.2010.04.005 doi: 10.1016/j.ecocom.2010.04.005
![]() |
[16] |
X. Q. Wang, Y. L. Cai, H. H. Ma, Dynamics of a diffusive predator-prey model with Allee effect on predator, Discrete Dyn. Nat. Soc., 2013 (2013), 984960. https://doi.org/10.1155/2013/984960 doi: 10.1155/2013/984960
![]() |
[17] |
R. D. Parshad, E. Quansah, K. Black, R. K. Upadhyay, S. K. Tiwari, N. Kumari, Long time dynamics of a three-species food chain model with Allee effect in the top predator, Comput. Math. Appl., 71 (2016), 503–528. https://doi.org/10.1016/j.camwa.2015.12.015 doi: 10.1016/j.camwa.2015.12.015
![]() |
[18] |
E. González-Olivares, H. Meneses-Alcay, B. González-Yañez, J. Mena-Lorca, A. Rojas-Palma, R. Ramos-Jiliberto, Multiple stability and uniqueness of the limit cycle in a Gause-type predator-prey model considering the Allee effect on prey, Nonlinear Anal. Real, 12 (2011), 2931–2942. https://doi.org/10.1016/j.nonrwa.2011.04.003 doi: 10.1016/j.nonrwa.2011.04.003
![]() |
[19] |
H. Y. Liu, H. G. Yu, C. J. Dai, Z. L. Ma, Q. Wang, M. Zhao, Dynamical analysis of an aquatic amensalism model with non-selective harvesting and Allee effect, Math. Biosci. Eng., 18 (2021), 8857–8882. https://doi.org/10.3934/mbe.2021437 doi: 10.3934/mbe.2021437
![]() |
[20] |
J. F. Wang, J. P. Shi, J. J. Wei, Dynamics and pattern formation in a diffusive predator-prey system with strong Allee effect in prey, J. Differ. Equations, 251 (2021), 1276–1304. https://doi.org/10.1016/j.jde.2011.03.004 doi: 10.1016/j.jde.2011.03.004
![]() |
[21] |
M. Y. Hamada, T. El-Azab, H. El-Metwally, Allee effect in a Ricker type predator-prey model, J. Math. Comput. Sci-JM., 29 (2023), 239–251. https://doi.org/10.22436/jmcs.029.03.03 doi: 10.22436/jmcs.029.03.03
![]() |
[22] |
S. Rana, S. Bhattacharya, S. Samanta, Spatiotemporal dynamics of Leslie-Gower predator-prey model with Allee effect on both populations, Math. Comput. Simulat., 200 (2022), 32–49. https://doi.org/10.1016/j.matcom.2022.04.011 doi: 10.1016/j.matcom.2022.04.011
![]() |
[23] |
J. Zu, M. Mimura, The impact of Allee effect on a predator-prey system with Holling type Ⅱ functional response, Appl. Math. Comput., 217 (2010), 3542–3556. https://doi.org/10.1016/j.amc.2010.09.029 doi: 10.1016/j.amc.2010.09.029
![]() |
[24] |
J. A. Dunne, R. J. Williams, N. D. Martinez, Network structure and biodiversity loss in food webs: Robustness increases with connectance, Ecol. Lett., 5 (2002), 558–567. https://doi.org/10.1046/j.1461-0248.2002.00354.x doi: 10.1046/j.1461-0248.2002.00354.x
![]() |
[25] |
B. Ghosh, T. K. Kar, Sustainable use of prey species in a prey-predator system: Jointly determined ecological thresholds and economic trade-offs, Ecol. Model., 272 (2014), 49–58. https://doi.org/10.1016/j.ecolmodel.2013.09.013 doi: 10.1016/j.ecolmodel.2013.09.013
![]() |
[26] |
X. X. Liu, Q. D. Huang, Analysis of optimal harvesting of a predator-prey model with Holling type Ⅳ functional response, Ecol. Complex., 42 (2020), 100816. https://doi.org/10.1016/j.ecocom.2020.100816 doi: 10.1016/j.ecocom.2020.100816
![]() |
[27] |
X. X. Liu, Q. D. Huang, Dynamics and optimal control of a Monod-Haldane predator-prey system with mixed harvesting, Int. J. Bifurcat. Chaos, 30 (2020), 2050243. https://doi.org/10.1142/S0218127420502430 doi: 10.1142/S0218127420502430
![]() |
[28] |
T. Das, R. N. Mukherjee, K. S. Chaudhuri, Harvesting of a prey-predator fishery in the presence of toxicity, Appl. Math. Model., 33 (2009), 2282–2292. https://doi.org/10.1016/j.apm.2008.06.008 doi: 10.1016/j.apm.2008.06.008
![]() |
[29] |
S. G. Mortoja, P. Panja, S. K. Mondal, Stability analysis of plankton-fish dynamics with cannibalism effect and proportionate harvesting on fish, Mathematics, 11 (2023), 3011. https://doi.org/10.3390/math11133011 doi: 10.3390/math11133011
![]() |
[30] |
C. J. Wei, L. S. Chen, Periodic solution and heteroclinic bifurcation in a predator-prey system with Allee effect and impulsive harvesting, Nonlinear Dyn., 76 (2014), 1109–1117. https://doi.org/10.1007/s11071-013-1194-z doi: 10.1007/s11071-013-1194-z
![]() |
[31] |
Z. C. Shang, Y. H. Qiao, Multiple bifurcations in a predator-prey system of modified Holling and Leslie type with double Allee effect and nonlinear harvesting, Math. Comput. Simulat., 205 (2023), 745–764. https://doi.org/10.1016/j.matcom.2022.10.028 doi: 10.1016/j.matcom.2022.10.028
![]() |
[32] |
D. Y. Wu, H. Y. Zhao, Y. Yuan, Complex dynamics of a diffusive predator-prey model with strong Allee effect and threshold harvesting, J. Math. Anal. Appl., 469 (2019), 982–1014. https://doi.org/10.1016/j.jmaa.2018.09.047 doi: 10.1016/j.jmaa.2018.09.047
![]() |
[33] |
S. S. Yu, J. Liu, X. J. Lin, Multiple positive periodic solutions of a Gause-type predator-prey model with Allee effect and functional responses, Aims Math., 5 (2020), 6135–6148. https://doi.org/10.3934/math.2020394 doi: 10.3934/math.2020394
![]() |
[34] |
R. P. Gupta, D. K. Yadav, Role of Allee effect and harvesting of a food-web system in the presence of scavengers, J. Biol. Syst., 30 (2022), 149–181. https://doi.org/10.1142/S021833902250005X doi: 10.1142/S021833902250005X
![]() |
[35] |
B. Barman, B. Ghosh, Dynamics of a spatially coupled model with delayed prey dispersal, Int. J. Model. Simul., 42 (2021), 400–414. https://doi.org/10.1080/02286203.2021.1926048 doi: 10.1080/02286203.2021.1926048
![]() |
[36] | C. Roberts, Ordinary differential equations: Applications, models, and computing, New York: Chapman and Hall, 2011. https://doi.org/10.1201/b15111 |
[37] | Y. A. Kuznetsov, Elements of applied bifurcation theory, New York: Springer, 2004. https://doi.org/10.1007/978-1-4757-3978-7 |
[38] | L. Perko, Differential equations and dynamical systems, New York: Springer, 1991. https://doi.org/10.1007/978-1-4684-0392-3 |
[39] |
B. Ghosh, T. K. Kar, T. Legović, Sustainability of exploited ecologically interdependent species, Popul. Ecol., 56 (2014), 527–537. http://doi.org/10.1007/s10144-014-0436-3 doi: 10.1007/s10144-014-0436-3
![]() |
[40] |
B. Ghosh, D. Pal, T. Legović, T. K. Kar, Harvesting induced stability and instability in a tri-trophic food chain, Math. Biosci., 304 (2018), 89–99. http://doi.org/10.1016/j.mbs.2018.08.003 doi: 10.1016/j.mbs.2018.08.003
![]() |
[41] |
B. Ghosh, T. K. Kar, T. Legović, Relationship between exploitation, oscillation, MSY and extinction, Math. Biosci., 256 (2014), 1–9. http://doi.org/10.1016/j.mbs.2014.07.005 doi: 10.1016/j.mbs.2014.07.005
![]() |
[42] | L. S. Pontryagin, Mathematical theory of optimal processes, London: Routledge, 1987. https://doi.org/10.1201/9780203749319 |
[43] | C. W. Clark, Mathematical bioeconomics: The optimal management of renewable resources, 1990. |
[44] |
A. Rojas-Palma, E. González-Olivares, Optimal harvesting in a predator-prey model with Allee effect and sigmoid functional response, Appl. Math. Model., 36 (2012), 1864–1874. https://doi.org/10.1016/j.apm.2011.07.081 doi: 10.1016/j.apm.2011.07.081
![]() |