Global stability in a modified Leslie-Gower type predation model assuming mutual interference among generalist predators

1.   Introduction

We are interested in positive classical solutions of semilinear heat equations $u_t = \Delta u + u^p$, in ${\mathbb R}^n$ or in a complete Riemannian manifold $(M, g)$ without boundary, where $p > 1$.

In their celebrated paper ^[7], Li and Yau showed how a Harnack inequality for the classical heat equation on a manifold with nonnegative Ricci tensor can be derived from a differential inequality for the logarithm of a solution. Subsequently, in the case of a manifold with nonnegative sectional curvatures and parallel Ricci tensor, Hamilton in ^[5] proved that the Harnack estimate of Li and Yau can actually be obtained as the trace of a full matrix inequality, under some more restrictive geometric assumptions. In some cases these inequalities can be useful in proving triviality of eternal solutions, see ^[1], moreover, "geometric" versions of them appear naturally and play a key role in the analysis of mean curvature flow and of Ricci flow (which are described by much more complicated systems of parabolic PDEs), see ^[4,6].

Our aim is to extend to the semilinear setting the matrix Harnack estimate of Li & Yau–type for the heat equation developed by Hamilton in ^[5].

We set some definitions and notations. In all of the paper, the Riemannian manifolds $(M, g)$ will be smooth, complete, connected and without boundary. We will denote with $\nabla$ the Levi–Civita connection of $(M, g)$ and $\Delta$ the associated Laplace–Beltrami operator and we assume that $(M, g)$ has nonnegative sectional curvatures and parallel Ricci tensor, that is, $\nabla {{\mathrm {Ric}}} = 0$. Finally, all the solutions we will consider are classical ($C^2$ in space and $C^1$ in time, at least).

Remark 1.1. The hypothesis that $M$ has parallel Ricci tensor and nonnegative sectional curvatures is satisfied on a torus or a sphere or a complex projective space, or a product of such, or a quotient of a product by a finite group of isometries.

Definition 1.2. A quintuple of real numbers $(a, b, c, d, \theta)$ is admissible, if the following inequalities

are satisfied.

The above system turns out to have actually solutions.

Proposition 1.3. There exists a nonempty cone $\mathcal{C}$ of admissible quintuples of parameters.

We can then state our main result.

Theorem 1.4. Suppose $(M, g)$ is a complete, $n$–dimensional Riemannian manifold without boundary, with nonnegative sectional curvatures and parallel Ricci tensor. Let $f = \log u$, where $u$ is a positive classical solution of

in $M \times (0, T)$. Then, for any $(a, b, c, d, \theta) \in \mathcal{C}$ there exists a constant $\varepsilon > 0$ such that

in the sense of bilinear forms, in $M\times(0, T)$, for all $1 < p < 1+G(a, b, c, d, \theta)$, where

given

and $G_2(a, b, c, d, \theta)$ the positive solution of

Taking the trace with the metric $g$ in inequality (1.3), we get a scalar Li & Yau–type inequality (which is actually weaker than the analogous one proved in ^[1])

and substituting $u = e^f$,

in $M\times(0, T)$.

Remark 1.5. As we mentioned, inequalities like (1.3) and (1.6) are relevant for ancient (and eternal) solutions $u$, that is, solutions defined in $M\times(-\infty, T)$ for some $T\in{\mathbb R}\cup\{+\infty\}$, since by a standard argument (see [1,Section 3], for instance) choosing suitable intervals for their application, they imply

in the sense of bilinear forms and

in $M\times(-\infty, T)$.

It is possible to have an explicit bound for the function $G$ defined in Theorem 1.4 only in terms of the dimension $n$ of the manifold, as shown in the next proposition, hence giving a lower bound for the range of exponents $p$ for which the results hold.

Proposition 1.6. There holds

where

and

In particular, for any $c > 0$

and

In the next sections we show these results, while in the last one we derive, along the lines of ^[7] (and ^[5]), some consequent Harnack–type local estimates for the solutions.

2.   Proof of Theorem 1.4

We follow the line of Hamilton in ^[5] and, by simplicity and clarity, we show the proof when $(M, g)$ is ${\mathbb R}^n$ with its canonical metric. In the general case of a complete $n$–dimensional Riemannian manifold, the extra curvature terms which appear in the computations because of the operations of interchanging covariant derivatives "have the right sign" in the final inequality, since we assumed nonnegative sectional curvatures and parallel Ricci tensor (see the end of this section). Furthermore, thanks to standard localization arguments (as explicitly shown in ^[1], see also ^[4]), in applying the maximum principle – on which the proof is based – we can argue as if we were in a compact case (all the maximum/minimum points there exist).

Let $u:M\times[0, T)\to{\mathbb R}$ be a positive solution of $u_t = \Delta u + u^p$. Setting $f = \log u$, we have

Moreover, by Eq (1.2), we also have the following relations:

where $f_i = \nabla_i f$, $f_{ij} = \nabla_j \nabla_i f$, $f_{ijk} = \nabla_k \nabla_j\nabla_i f$ and so on. Let $(a, b, c, d, \theta)\in \mathcal{C}$ and define the symmetric bilinear form

Observe that

and, by a simple computation,

where

Define ${{\mathop f\limits^\circ} _{ij}}: = f_{ij}-(\Delta f /n) g_{ij}$ to be the tracefree part of the Hessian of $f$. Since

we have

In order to estimate the last three terms, we use the following algebraic inequality for traceless symmetric bilinear forms

where ${{\mathop {\rho}\limits^\circ} } = \max\{|\lambda_i|: \ \lambda_i \;{\rm{ eigenvalue \;of }} \;{{\mathop f\limits^\circ} _{i, j}}\}$ denotes the spectral radius of ${{\mathop f\limits^\circ} _{ij}}$. By Young's inequality, for every $K > 0$, one has

We set

where $K > 0$, by the first two inequalities of system (1.1) since $(a, b, c, d, \theta)\in \mathcal{C}$. Using inequality (2.5) to estimate the last three terms of Eq (2.4) we achieve

where the last inequality follows from the third inequality in system (1.1) and since $(a, b, c, d, \theta)\in \mathcal{C}$. The computation above yields

Now we claim that there exists $\varepsilon > 0$ such that if $F_{ij}\le 0$ then

where $F^2_{ij} = F_{ik} \, F_{kj}$. Thus, let us suppose that $F_{ij}\le 0$. Combining the estimate (2.6) with the definition of $Q_{ij}$ given in (2.3) we get

Recalling that for $(a, b, c, d, \theta)\in \mathcal{C}$ there holds $(\theta-b) > 0$, we only have to show that

Setting $x = p-1 \ge 0$, we can recast the previous inequalities as

Being $4cd(\theta-b) > 0$, $(d-a)\theta^2 > 0$ and $(d-c)\theta^2 > 0$, Eq (1.5) admits two solutions, of which only one is positive. Hence, defining $G_1(b, d, \theta)$ and $G_2(a, b, c, d, \theta)$ as in the statement of the Theorem and $G(a, b, c, d, \theta)$ as in (1.4), the fact that $0 \le x \le G(a, b, c, d, \theta)$ implies that $x$ satisfies both the inequalities in (2.8). Thus, Eq (2.7) is satisfied with $\varepsilon = \frac{1}{2(\theta-b)}$.

The application of maximum principle (see for example Theorem C.1.3 in ^[8] or Lemma 8.2 in ^[3]) concludes the proof.

We point out that in the case of a general Riemannian manifold extra curvature terms appear in the computations. For instance, the second and the fifth line in formulas (2.1) would become respectively

and

where $R_{ij}$ and $R_{ikjl}$ denote the components of the Ricci and the curvature tensor, respectively (we follow the convention in ^[2]). At this point it is straightforward to verify that the extra curvature terms which would appear in Eq (2.2) either disappear or have the right sign under our geometric assumptions, as in [5,Theorem 4.3].

3.   Proof of Proposition 1.3

First of all, observe that if $(a, b, c, d, \theta)$ is admissible and $\lambda > 0$ is a constant, then $(\lambda a, \lambda b, \lambda c, \lambda d, \lambda \theta)$ is still admissible. Thus, if we show that an admissible quintuple of parameters exists, then we have a cone of admissible parameters. Let us consider a quintuple of the form $(a, b, c, a, \theta)$ where $a-c > 0$ and $\theta-b > 0$. We want to find $a, b, c, \theta$ such that the third inequality of system (1.1) is satisfied. To do this, let us set $a-c = \delta$ for some $\delta > 0$. Then, the third inequality of system (1.1) becomes

Now let us assume that $\theta-b = c$ and $b = kc$ for some $k > 0$ to achieve

Being $c > 0$, we can set $z = \frac{\delta}{c}$ and recast the previous inequality as

Let us set

We want to find some suitable values for $z$ and $k$ such that $H(z, k)\ge 0$. First of all, let us determine some necessary conditions on $k$ in such a way that $H(\cdot, k)$ admits positive roots.

Proposition 3.1. There always exists a negative number $z_- < 0$ such that $H(z_-, k) = 0$. If there exist two positive numbers $z_1, z_2 > 0$ such that $H(z_i, k) = 0$, then $k > \sqrt{3n}$.

Proof. Let us observe that $H(0, k) = -(n-1) < 0$ and $\lim_{z \to -\infty}H(z, k) = +\infty$, hence we obviously have $z_-$. Moreover, let us observe that, by Descartes' rule of the signs, the maximum number of positive solutions of the equation $H(z, k) = 0$ for fixed $k$ is $0$ if $k^2-3n\le 0$ and $2$ if $k^2-3n > 0$, thus, if two positive solutions $z_1, z_2$ exist, then $k > \sqrt{3n}$.

We can use the previous necessary condition to achieve a necessary and sufficient condition on the existence of two positive roots.

Proposition 3.2. Fix $n \ge 2$. Then $H(z, k)$ admits two (possibly equal) positive roots if and only if $k \ge k(n)$, where

Proof. As we have shown before, a necessary condition to have two positive roots is that $k > \sqrt{3n}$. On the other hand, let us observe that

hence, by Descartes' rule of the signs we know that $H(\cdot, z)$ admits at most one negative root. Moreover, since $H(0, k) = -(n-1) < 0$ and $\lim_{z \to -\infty}H(z, k) = +\infty$, then we know that $H(z, k)$ admits exactly one negative root for any $k > \sqrt{3n}$. Thus, if $H(z, k)$ admits three real roots, two of them have to be positive.

The discriminant of $H(\cdot, k)$ is given by

We have that $\Delta_H(k)\ge 0$ if and only if $P_1(k): = k^3-\frac{27}{4}kn-\frac{27}{4}n\ge 0$. $P_1(k)$ is a depressed cubic polynomial with $p_1(n) = q_1(n) = -\frac{27}{4}n$. By Descartes' rule of the signs and the fact that $P_1(0) = -\frac{27}{4}n < 0$ we know that $P_1(k)$ always admits a unique positive root. The discriminant of $P_1(k)$ is given by

thus $P_1(k)$ admits three real roots. Since we are under the casus irreducibilis, we have to provide trigonometric solutions to recognize what is the real solution we are interested in. To do this, we will use Viéte's procedure. Consider the equation $P_1(k) = 0$ and set $k = u\cos(\theta)$ to achieve

Multiplying everything by $\frac{4}{u^3}$ we get

Now set $\frac{27n}{u^2} = 3$, that is to say $u = 3\sqrt{n}$, to achieve

Recalling that $4\cos^3(\theta)-3\cos(\theta) = \cos(3\theta)$, we get

and then

where $j = 0, 1, 2$. Finally, we get the three real roots of $P_1(k)$ as

However, we know that $P_1(k)$ admits only one positive root. Being

we have that $k_0 > 0$ is the solution we are searching for. Thus, let us relabel

to conclude that $P_1(k)\ge 0$ if and only if $k \ge k(n)$.

Now we can conclude the proof of Proposition 1.3. Indeed, let us consider $k > k(n)$, in such a way that $\Delta_H(k) > 0$. For such fixed $k$, $H(\cdot, k)$ admits two positive roots $0 < z_1 < z_2$. Consider any $z \in [z_1, z_2]$. Setting, without loss of generality, $c = 1$, we then know that the quintuple $(z+1, k, 1, z+1, k+1)$ is admissible.

4.   Proof of Proposition 1.6

To prove Proposition 1.6, we want to exhibit an admissible quintuple $(a, b, c, d, \theta)$ such that $G(a, b, c, d, \theta)$ can be explicitly calculated. To do this, let us consider again a quintuple of the form $(z+1, k, 1, z+1, k+1)$ where $k \ge k(n)$. Moreover, let us consider $z(k, n)$ to be a positive local maximum point of $H(z, k)$, for fixed $k \ge k(n)$, such that $H(z(k, n), k) \ge 0$. If $k > k(n)$, then this maximum always exists, by a simple application of Rolle's theorem on the interval $[z_1, z_2]$, together with the fact that both $z_1, z_2$ are simple roots. If $k = k(n)$, then $z_1 = z_2$ and it coincides with such local maximum of the polynomial $H(z, k(n))$. Let us evaluate it explicitly.

Proposition 4.1. For fixed $n \ge 2$ and $k > \sqrt{3n}$, $H(z, k)$ admits two (eventually equal) critical points if and only if these critical points are positive and $k \ge k_0(n)$ where

If $k > k_0(n)$, the local maximum point is given by

In particular, $k_0(n) \le k(n)$ and, if $k \ge k(n)$, there holds $H(z(k, n), k)\ge 0$.

Proof. Let us first observe that

Since $k > \sqrt{3n}$, Descartes' rule of the signs tells us that, if the solutions exist, they must be positive. Hence we have only to show that the discriminant is nonnegative. Let us determine the discriminant of the polynomial (4.1):

Being $k > 0$, $\Delta_z(k)\ge 0$ if and only if

Let us first show that this polynomial admits a unique positive root. Observe that $P_2(0) = -6n < 0$ and $\lim_{k \to +\infty}P_2(k) = +\infty$, thus $P_2(k)$ admits a positive root $k_0(n)$. Moreover, Descartes' rule of the signs tells us that $P_2(k)$ admits at most one positive root, hence $k_0(n)$ is the unique positive root.

Now let us determine $k_0(n)$. First of all, let us observe that, since $P_2(k)$ is a depressed cubic polynomial, its discriminant is given by

since $n \ge 2$. Thus we know $P_2(k)$ admits three different roots. Since we are under the casus irreducibilis, we have to provide trigonometric solutions to recognize what is the real solution we are interested in. Arguing again by Viéte's procedure and selecting the unique positive root, we get

Hence, as $k > k_0(n)$, we can find two solutions to equation $H'(z, k) = 0$. Being $H(\cdot, k)$ a polynomial with $H(0, k) = -(n-1)$ and $\lim_{z \to -\infty}H(z, k) = +\infty$, we know that $H(z, k)$ is decreasing as $z \le 0$. Thus, the first critical point has to be a local minimum and the second critical point a local maximum. Therefore, writing explicitly the second solution of $H'(z, k) = 0$, we get

Finally, observe that, by Rolle's theorem, if $k > k(n)$, then the interval $[z_1, z_2]$ has to admit one of the two critical points. Moreover, since $H(z, k) > 0$ as $z \in (z_1, z_2)$ and $H(z_1, k) = H(z_2, k) = 0$, then such critical point is the local maximum and $H(z(k, n), k)$ is positive. If $k = k(n)$, then $z_1 = z_2$ is a double root and $H(z_1, k(n)) = H(z_2, k(n)) = 0$. Being $H(z, k(n))\le 0$ for any $z \ge 0$ (since the other simple root is negative), we have that $z_1 = z_2 = z(k(n), n)$ and $H(z(k(n), n), k(n)) = 0$. This also obviously implies $k_0(n) < k(n)$.

Next, we want to evaluate $G$ on the quintuple $(z(k, n)+1, k, 1, z(k, n)+1, k+1)$ as $k \ge k(n)$. To do this, we first need to exploit a simple property of $z(k, n)$.

Proposition 4.2. For any $k \ge k(n)$ there holds $z(k, n)\ge 2$.

Proof. By Eq (4.3), we have that $z(k, n)\ge 2$ if and only if

Being $k \ge k(n)$, we already know that the quantity under the square root is nonnegative. On the other hand, if $k\ge 3\sqrt{n}$, then inequality (4.4) holds true. Let us consider $k(n)\le k < 3\sqrt{n}$. Then inequality (4.4) is equivalent to

that is to say

Since $k(n) > 0$, the previous inequality is verified as

To conclude the proof, we have to show that $k_1(n)\le k(n)$. Being $k_0(n)\le k(n)$, this is obvious if $k_1(n)\le k_0(n)$, thus let us suppose $k_1(n) > k_0(n)$.

Let us consider the function $g(k): = H(z(k, n), k)$ for $k \ge k_0(n)$. Being $z(k, n)$ a local maximum of $H$, there holds $\frac{\partial \, H}{\partial \, z}(z(k, n), k) = 0$. Thus, we have

In particular $g'(k)\ge 0$ if and only if $kz(k, n)\ge 1$. The latter holds if and only if

where $\Delta_z(k)$ is defined in Eq (4.2). Being $k \ge k_0(n)$, we have that $\Delta_z(k)\ge 0$ and then inequality (4.5) is verified. This implies, in particular, that $g(k)$ is increasing as $k \ge k_0(n)$. Moreover, by definition of $k(n)$, there holds $g(k(n))\ge 0$. On the other hand, being $z(k_1(n), n) = 2$ and $k_1(n) > \frac{1}{2}$, we get

Hence, we have $k_1(n) < k(n)$, concluding the proof.

With this property in mind, we can evaluate

To do this, let us first observe that

Concerning $G_2(z(k, n)+1, k, 1, z(k, n)+1, k+1)$, it is the unique solution of

that is to say

Hence, we have

However, being $z(k, n)\ge 2$ by the previous proposition, there holds

then

Now we want to optimize on $k \ge k(n)$. To do this, let us show the following Proposition.

Proposition 4.3. The function $k \mapsto z(k, n)$ is increasing on $[k(n), +\infty)$.

Proof. Just observe that

being $\Delta_z(k)\ge 0$ by the fact that $k \ge k(n)\ge k_0(n)$.

The previous proposition implies that the function $k \mapsto G(z(k, n)+1, k, 1, z(k, n)+1, k+1)$ is decreasing as $k \ge k(n)$, thus it achieve its maximum value as $k = k(n)$. Setting

we conclude the proof of Proposition 1.6.

5.   Matrix Harnack inequalities

We begin the section with the following technical lemma.

Lemma 5.1. There exists a cone $\mathcal{C}'$ of an admissible quintuples $(a, b, c, d, \theta)$ such that $a = d$ and $b-a+c < 0$.

Proof. Let us consider $k \geq k(n)$ with $k(n)$ defined in Proposition 1.6 and let $(z(k, n)+1, k, 1, z(k, n)+1, k+1) \in \mathcal{C}$ with $z(k, n)$ defined in (4.3). The condition $b-a+c < 0$ becomes $k < z(k, n)$, i.e.,

This is obviously true if

With this choice of $k$ we conclude the proof.

Remark 5.2. It is easy to see that if $(a, b, c, d, \theta) \in \mathcal{C}$ then also $(a, b, c, a, \theta) \in \mathcal{C}$.

We will now derive some Harnack–type inequalities as a consequence of Theorem 1.4.

Proposition 5.3. Let $(M, g)$ be an $n$–dimensional complete Riemannian manifold with nonnegative sectional curvatures and parallel Ricci tensor. Consider $(a, b, c, a, \theta) \in \mathcal{C'}$ and $1 < p < 1 + G(a.b, c, a, \theta)$. Let $u:M\times[0, T)\to{\mathbb R}$ be a classical positive solution of the equation $u_t = \Delta u + u^p$, then there exists $\varepsilon = \varepsilon(n, p, a, b, c, \theta)$ such that, given any $0 < t_1 < t_2\leq T$ and $x_1, x_2 \in M$, the following inequality holds

where

with $\Gamma(x_1, x_2)$ given by all the paths in $M$ parametrized by $[0, 1]$ joining $x_2$ to $x_1$, $f = \log u$ and $\rho = \max \{ |\lambda_i |: \, \lambda_i eigenvalue of f_{ij} \}$.

Proof. By Theorem 1.4 we know that there exists $\varepsilon$ such that

Recalling that $f$ satisfies

we get

Since $\rho$ is the spectral radius of $f_{ij}$ we get $f_{ij} \leq \rho g_{ij}$. On the other hand the matrix $f_i f_j$ has rank one and the only nonzero eigenvalue is $|\nabla f|^2$, then $f_i f_j \leq |\nabla f|^2 g_{ij}$. Plugging these inequalities into (5.1) we obtain

Now let us consider any $\gamma \in \Gamma(x_1, x_2)$ as well as $\eta:[0, 1]\to M\times[t_1, t_2]$ defined as $\eta(s) = (\gamma(s), (1-s)t_2+st_1)$. Evaluating (5.2) in $\eta(s)$ and applying it to $\dot{\gamma}(s)$ we see that

Noticing that $\eta(0) = (x_2, t_2)$ and $\eta(1) = (x_1, t_1)$, we have

Being $\gamma \in \Gamma(x_1, x_2)$ arbitrary, we conclude the proof.

With the same strategy we can also prove the following variant.

Proposition 5.4. Let $(M, g)$ be an $n$–dimensional complete Riemannian manifold with nonnegative sectional curvatures and parallel Ricci tensor. Consider $(a, b, c, a, \theta) \in \mathcal{C'}$ and $1 < p < 1 + G(a.b, c, a, \theta)$. Let $u:M\times[0, T)\to{\mathbb R}$ a classical positive solution of the equation $u_t = \Delta u + u^p$, then there exists $\varepsilon = \varepsilon(n, p, a, b, c, \theta)$ such that, given any $0 < t_1 < t_2\leq T$ and $x_1, x_2 \in M$, the following inequality holds

where

with $\Gamma(x_1, x_2)$ given by all the paths in $M$ parametrized by $[0, 1]$ joining $x_2$ to $x_1$, $f = \log u$, ${{\mathop f\limits^\circ} _{ij}} = f_{ij}- \frac{\Delta f}{n} g_{ij}$ and ${{\mathop {\rho}\limits^\circ} } = \max \{ |\lambda_i |: \, \lambda_i \;eigenvalue \;of\; {{\mathop f\limits^\circ} _{ij}} \}$.

Finally, we also have another one coming from the scalar (trace) version of the Li & Yau inequality (1.6). We underline that it holds up to an exponent $p$ lower than the analogous one for the Harnack inequality obtained in ^[1], moreover, this latter holds under the only hypothesis of nonnegative Ricci tensor.

Proposition 5.5. Let $(M, g)$ be an $n$–dimensional complete Riemannian manifold with nonnegative sectional curvatures and parallel Ricci tensor. Consider $(a, b, c, a, \theta) \in \mathcal{C}$ and $1 < p < 1 + G(a.b, c, a, \theta)$. Let $u:M\times[0, T)\to{\mathbb R}$ a classical positive solution of the equation $u_t = \Delta u + u^p$, then there exists $\varepsilon = \varepsilon(n, p, a, b, c, \theta)$ such that, given any $0 < t_1 < t_2\leq T$ and $x_1, x_2 \in M$, the following inequality holds

where

with $\Gamma(x_1, x_2)$ given by all the paths in $M$ parametrized by $[0, 1]$ joining $x_2$ to $x_1$.

Proof. Let $f = \log u$, by Theorem 1.4 we know that there exists $\varepsilon > 0$ such that inequality (1.6) holds, then the proof proceeds as in Proposition 5.3.

Conflict of interest

The authors declare no conflict of interest.