We discuss linear autonomous evolution equations on function spaces which have the property that a positive initial value leads to a solution which initially changes sign, but then becomes – and stays – positive again for sufficiently large times. This eventual positivity phenomenon has recently been discovered for various classes of differential equations, but so far a general theory to explain this type of behaviour exists only under additional spectral assumptions.

1 Evolution equations and positivity

To set the stage, we start with a reminder about linear evolution equations whose solutions are positive whenever the initial value is.

Linear ODEs and positivity

For a matrix ARd×dA\in\mathbb{R}^{d\times d}, the linear and autonomous initial value problem

{u˙(t)=Au(t)for t[0,),u(0)=u0,\begin{cases}\dot{u}(t)=Au(t)\quad&\textrm{for}\ t\in[0,\infty),\\ u(0)=u_{0},\end{cases}

where u0Rdu_{0}\in\mathbb{R}^{d}, is well-known to be solved by the function

u ⁣:[0,)tetAu0Rd.u\colon[0,\infty)\ni t\mapsto e^{tA}u_{0}\in\mathbb{R}^{d}.

We say that the matrix family (etA)t[0,)(e^{tA})_{t\in[0,\infty)} is positive if etAu00e^{tA}u_{0}\geq 0 for all t[0,)t\in[0,\infty) whenever u00u_{0}\geq 0; equivalently, etA0e^{tA}\geq 0 for all t[0,)t\in[0,\infty). Here, we use the notation 0\geq 0 for a vector or a matrix to say that all its entries are 0\geq 0.

Remark 1. There is some terminological inconsistency in the literature with respect to this notion: in matrix analysis and in some parts of PDE theory, it is common to use the word non-negativity; we use the notion positivity instead, which is more common in functional analysis.

To get an intuition for this concept, it is useful to recall that positivity of the matrix exponential function is easy to characterise in terms of AA.

Theorem 2.

For ARd×dA\in\mathbb{R}^{d\times d}, the family (etA)t[0,)(e^{tA})_{t\in[0,\infty)} is positive if and only if every off-diagonal entry of AA is 0\geq 0.

Proof. \Rightarrow” For indices jkj\neq k, one has

Ajk=limt0ej,etAidtek=limt01tej,etAek0,A_{jk}=\lim_{t\downarrow 0}\Bigl\langle e_{j},\frac{e^{tA}-\mathrm{id}}{t}e_{k}\Bigr\rangle =\lim_{t\downarrow 0}\frac{1}{t}\langle e_{j},e^{tA}e_{k}\rangle\geq 0,

where ej,ekRde_{j},e_{k}\in\mathbb{R}^{d} are the canonical unit vectors and ,\langle{\,\cdot\,},{\cdot\,}\rangle denotes the standard inner product on Rd\mathbb{R}^{d}.

\Leftarrow” By assumption, one has, for a sufficiently large number c0c\geq 0, the inequality A+cid0A+c\,\mathrm{id}\geq 0, and hence

etA=etcet(A+cid)0e^{tA}=e^{-tc}e^{t(A+c\,\mathrm{id})}\geq 0

for all t[0,)t\in[0,\infty), where the inequality at the end follows from the series expansion of the matrix exponential function. ∎

A typical situation where positivity of matrix exponential functions occurs is the study of Markov processes on finite state spaces.

Example 3. Assume that all off-diagonal entries of ARd×dA\in\mathbb{R}^{d\times d} are 0\geq 0 and that all rows of AA sum up to 00. Then (etA)t[0,)(e^{tA})_{t\in[0,\infty)} is positive, and the vector 1Rd\operatorname{\mathbb{1}}\in\mathbb{R}^{d} whose entries are all equal to 11 satisfies A1=0A\operatorname{\mathbb{1}}=0 and thus etA1=1e^{tA}\operatorname{\mathbb{1}}=\operatorname{\mathbb{1}} for all t[0,)t\in[0,\infty). This shows that each of the matrices etAe^{tA}, t0t\geq 0, is row stochastic, so (etA)t[0,)(e^{tA})_{t\in[0,\infty)} describes a continuous-time Markov process on the finite state space {1,,d}\{1,\dots,d\}.

Infinite-dimensional equations

In infinite dimensions, we are still interested in initial value problems of the form

{u˙(t)=Au(t)for t[0,),u(0)=u0,\begin{cases}\dot{u}(t)=Au(t)\quad&\textrm{for}\ t\in[0,\infty),\\ u(0)=u_{0},\end{cases}

but this time, u0u_{0} is an element of a Banach space EE, and A ⁣:Edom(A)EA\colon E\supseteq\operatorname{dom}(A)\to E is a linear operator which is defined on a vector subspace dom(A)\operatorname{dom}(A) of EE. The initial value problem is well-posed if and only if AA is a generator of a C0C_{0}-semigroup(etA)t[0,)(e^{tA})_{t\in[0,\infty)}. Such a C0C_{0}-semigroup is a family of bounded linear operators on EE which is a suitable infinite-dimensional substitute of the matrix exponential function and has similar properties, but it is not given by an exponential series in general. The solution uu to the initial value problem is then given, again, by the formula u(t)=etAu0u(t)=e^{tA}u_{0} for t[0,)t\in[0,\infty). The generator and the C0C_{0}-semigroup determine each other uniquely, and the relation between semigroup and generator can in general be expressed by the formula

dom(A)={vE:limt01t(etAid)v exists in E},\displaystyle\operatorname{dom}(A)=\Bigl\{v\in E:\lim_{t\downarrow 0}\frac{1}{t}(e^{tA}-\mathrm{id})v\ \textrm{exists in}\ E\Bigr\},
Av=limt01t(etAid)v.\displaystyle Av=\lim_{t\downarrow 0}\frac{1}{t}(e^{tA}-\mathrm{id})v.

The following notion will be used several times later on. For a linear operator A ⁣:Edom(A)XA\colon E\supseteq\operatorname{dom}(A)\to X on a Banach space XX, the quantity

s(A):sup{Reλ:λσ(A)}[,],s(A)\coloneq\sup\{\operatorname{Re}\lambda:\lambda\in\sigma(A)\}\in[-\infty,\infty],

where σ(A)\sigma(A) denotes the spectrum of AA, is called the spectral bound of AA. If AA generates a C0C_{0}-semigroup, then s(A)<s(A)<\infty (see [20 K.-J. Engel and R. Nagel, One-parameter semigroups for linear evolution equations. Grad. Texts in Math. 194, Springer, New York (2000) , Theorem II.1.10 (ii)]). More information about C0C_{0}-semigroup theory can be found, for instance, in the monographs [31 A. Pazy, Semigroups of linear operators and applications to partial differential equations. Appl. Math. Sci. 44, Springer, New York (1983) , 20 K.-J. Engel and R. Nagel, One-parameter semigroups for linear evolution equations. Grad. Texts in Math. 194, Springer, New York (2000) ].

Let us briefly illustrate the concept of a C0C_{0}-semigroup by two very classical examples.

Examples 4. (a) Let p(1,)p\in(1,\infty) and let the operator AA be the Laplace operator on the space Lp(Rn)L^{p}(\mathbb{R}^{n}), i.e.

dom(A)=W2,p(Rn),\displaystyle\operatorname{dom}(A)=W^{2,p}(\mathbb{R}^{n}),
Av=Δv:j=1nj2vfor vdom(A).\displaystyle Av=\Delta v\coloneq\sum_{j=1}^{n}\partial_{j}^{2}v\quad\textrm{for}\ v\in\operatorname{dom}(A).

Then AA generates a C0C_{0}-semigroup (etA)t[0,)(e^{tA})_{t\in[0,\infty)} on Lp(Rn)L^{p}(\mathbb{R}^{n}) that is given by the formula

(etAu0)(x)=1(4πt)n/2Rnexp(xy224t)u0(y) ⁣dy(e^{tA}u_{0})(x)=\frac{1}{(4\pi t)^{n/2}}\int_{\mathbb{R}^{n}}\exp\Bigl(-\frac{\lVert x-y\rVert_{2}^{2}}{4t}\Bigr)u_{0}(y)\mathop{}\!\mathrm{d}y

for u0Lp(Rn)u_{0}\in L^{p}(\mathbb{R}^{n}) and xRnx\in\mathbb{R}^{n}. The semigroup is called the heat semigroup since it describes the solutions to the heat equation

u˙(t)=Δu(t).\dot{u}(t)=\Delta u(t).

Similar observations can be made on the space L1(Rn)L^{1}(\mathbb{R}^{n}), but the domain of the Laplace operator cannot be chosen to be a Sobolev space in that case, due to the lack of elliptic regularity.

(b) Let p[1,)p\in[1,\infty) and let the operator AA be the negative first derivative on Lp(0,)L^{p}(0,\infty), given by

dom(A)={vW1,p(0,):v(0)=0},Av=v.\operatorname{dom}(A)=\{v\in W^{1,p}(0,\infty):v(0)=0\},\quad Av=-v^{\prime}.

Then AA generates the so-called right shift semigroup(etA)t[0,)(e^{tA})_{t\in[0,\infty)} on Lp(0,)L^{p}(0,\infty) given by

(etAu0)(x)={u0(xt)if tx,0if t>x,(e^{tA}u_{0})(x)=\begin{cases}u_{0}(x-t)\quad&\textrm{if}\ t\leq x,\\ 0\quad&\textrm{if}\ t>x,\end{cases}

for u0Lp(0,)u_{0}\in L^{p}(0,\infty). The mapping u ⁣:[0,)tetAu0Lp(0,)u\colon[0,\infty)\ni t\mapsto e^{tA}u_{0}\in L^{p}(0,\infty) is a so-called mild solution to the transport equation

{u˙(t,x)=xu(t,x)for t,x>0,u(0,x)=u0(x)for x>0,u(t,0)=0for t>0;\left\{\begin{aligned} \dot{u}(t,x)&=-\partial_{x}u(t,x)&&\textrm{for}\ t,x>0,\\ u(0,x)&=u_{0}(x)&&\textrm{for}\ x>0,\\ u(t,0)&=0&&\textrm{for}\ t>0;\end{aligned}\right.

see [20 K.-J. Engel and R. Nagel, One-parameter semigroups for linear evolution equations. Grad. Texts in Math. 194, Springer, New York (2000) , Definition II.6.3] for the definition of mild solutions. This example is an easy illustration of the general principle that boundary conditions of a PDE are encoded in the domain of the corresponding operator AA.

Positive C₀-semigroups

In order to discuss positiveC0C_{0}-semigroups, one needs an order structure on the underlying Banach space EE. This can be for instance a partial order induced by a general closed convex cone, or more specifically the order structure of a Banach lattice. To facilitate the exposition here, we will restrict our attention to the illustrative case of function spaces, most importantly to LpL^{p}-spaces (over σ\sigma-finite measure spaces).

For a function fLpf\in L^{p}, we write f0f\geq 0 to indicate that f(ω)0f(\omega)\geq 0 for almost all ω\omega. In accordance with the terminology used above, we call a function ffpositive if it satisfies f0f\geq 0. A C0C_{0}-semigroup (etA)t[0,)(e^{tA})_{t\in[0,\infty)} on LpL^{p} is called positive if etAu00e^{tA}u_{0}\geq 0 for all t[0,)t\in[0,\infty) whenever u00u_{0}\geq 0. Equivalently, each of the operators etAe^{tA} is positive – which we denote by etA0e^{tA}\geq 0 – in the sense that it maps positive functions to positive functions.

We have already encountered two examples of positive C0C_{0}-semigroups: as is easy to see, both semigroups in Examples 4 are positive.

2 Positivity for large times

Let us now proceed to a more surprising situation, where positive initial values lead to solutions which might change sign at first, but again become – and stay – positive for sufficiently large times. In this section, we illustrate by means of two easy examples that this kind of behaviour can indeed occur; a more systematic account is presented in the subsequent section.

A matrix example

Let us start with a simple three-dimensional example.

Example 5. Consider the orthonormal basis B\mathcal{B} of R3\mathbb{R}^{3} consisting of the three vectors

v1=13(111),v2=12(101),v3=16(121).v_{1}=\frac{1}{\sqrt{3}}\begin{pmatrix}1\\ 1\\ 1\end{pmatrix},\quad v_{2}=\frac{1}{\sqrt{2}}\begin{pmatrix}-1\\ \hphantom{-}0\\ \hphantom{-}1\end{pmatrix},\quad v_{3}=\frac{1}{\sqrt{6}}\begin{pmatrix}\hphantom{-}1\\ -2\\ \hphantom{-}1\end{pmatrix}.

Let AR3×3A\in\mathbb{R}^{3\times 3} be such that its representation matrix with respect to the basis B\mathcal{B} is given by

R:[000011011],R\coloneq\begin{bmatrix}0&\hphantom{-}0&\hphantom{-}0\\ 0&-1&-1\\ 0&\hphantom{-}1&-1\end{bmatrix},

i.e., we let A=VRV1A=VRV^{-1}, where VR3×3V\in\mathbb{R}^{3\times 3} consists of the columns v1,v2,v3v_{1},v_{2},v_{3}. A direct computation shows that AA has some strictly negative off-diagonal entries, so (etA)t[0,)(e^{tA})_{t\in[0,\infty)} is not positive according to Theorem 2. On the other hand, AA has the eigenvalue 00 (with eigenvector v1v_{1}) as well as the further eigenvalues 1±i-1\pm i, so etAe^{tA} converges to the matrix v1v1Tv_{1}\cdot v_{1}^{\mathrm{T}}, whose entries are all equal to 1/31/3, as tt\to\infty; this shows that etAe^{tA} is a positive matrix for all sufficiently large times tt.

A fourth order PDE

Let us now discuss an infinite-dimensional example where eventual positivity occurs.

Example 6. Let us consider the biharmonic heat equation with periodic boundary conditions on L2(0,1)L^{2}(0,1). It is given by

u˙(t)=Au(t)for t[0,),\dot{u}(t)=Au(t)\quad\textrm{for}\ t\in[0,\infty),

where A ⁣:L2(0,1)dom(A)L2(0,1)A\colon L^{2}(0,1)\supseteq\operatorname{dom}(A)\to L^{2}(0,1) has domain

dom(A)={vH4(0,1):v(k)(0)=v(k)(1) for k=0,1,2,3}\operatorname{dom}(A)=\{v\in H^{4}(0,1):v^{(k)}(0)=v^{(k)}(1)\ \textrm{for}\ k=0,1,2,3\}

and is given by Av=v(4)Av=-v^{(4)} for each vdom(A)v\in\operatorname{dom}(A). The C0C_{0}-semigroup (etA)t[0,)(e^{tA})_{t\in[0,\infty)} is not positive; this can for instance be seen by associating a sesquilinear form to A-A and using the so-called Beurling–Deny criterion [30 E. M. Ouhabaz, Analysis of heat equations on domains. London Math. Soc. Monogr. Ser. 31, Princeton University Press, Princeton, NJ (2005) , Corollary 2.18].

However, we can prove positivity for large times. To this end, note that the operator AA is self-adjoint, and its spectrum consists of isolated eigenvalues only since dom(A)\operatorname{dom}(A) embeds compactly into L2(0,1)L^{2}(0,1). The largest eigenvalue of AA is 00, and the constant function 1\operatorname{\mathbb{1}} spans the corresponding eigenspace. Hence we conclude, for instance from the spectral theorem for self-adjoint operators with compact resolvent, that

etAu0u0,11:01u0(x) ⁣dx1in L2(0,1)e^{tA}u_{0}\to\langle u_{0},\operatorname{\mathbb{1}}\rangle\operatorname{\mathbb{1}}\coloneq\int_{0}^{1}u_{0}(x)\mathop{}\!\mathrm{d}x\cdot\operatorname{\mathbb{1}}\quad\textrm{in}\ L^{2}(0,1)

for each u0L2(0,1)u_{0}\in L^{2}(0,1) as tt\to\infty. Since AA is self-adjoint, the operators etAe^{tA} have the property that for t>0t>0 they map L2(0,1)L^{2}(0,1) into dom(A)\operatorname{dom}(A) and thus into L(0,1)L^{\infty}(0,1). Moreover, they are even continuous from L2(0,1)L^{2}(0,1) to L(0,1)L^{\infty}(0,1) (this follows for instance from the closed graph theorem), so for u0L2(0,1)u_{0}\in L^{2}(0,1), we even have

etAu0=e1Ae(t1)Au0u0,1e1A1=u0,11e^{tA}u_{0}=e^{1\cdot A}e^{(t-1)A}u_{0}\to\langle u_{0},\operatorname{\mathbb{1}}\rangle\,e^{1\cdot A}\operatorname{\mathbb{1}}=\langle u_{0},\operatorname{\mathbb{1}}\rangle\operatorname{\mathbb{1}}

as tt\to\infty, where the convergence takes place with respect to the norm in L(0,1)L^{\infty}(0,1). This implies that if u00u_{0}\geq 0, then etAu00e^{tA}u_{0}\geq 0 for all sufficiently large times tt.

3 A systematic theory

After the previous ad hoc examples, we now present a few excerpts of a more systematic account of eventual positivity.

Eventually positive matrix semigroups

Example 5 already gives quite a straightforward idea of how to obtain a sufficient condition for a matrix exponential function to be eventually positive: if a matrix ARd×dA\in\mathbb{R}^{d\times d} has a simple real eigenvalue that dominates the real parts of all other eigenvalues and if the corresponding eigenvectors of AA and the transposed matrix ATA^{\mathrm{T}} have strictly positive entries only, then we expect etAe^{tA} to be positive – and in fact to even have strictly positive entries only – for all sufficiently large tt. A bit more surprising is the Perron–Frobenius-like fact that the converse implication also holds. This was proved by Noutsos and Tsatsomeros in [29 D. Noutsos and M. J. Tsatsomeros, Reachability and holdability of nonnegative states. SIAM J. Matrix Anal. Appl. 30, 700–712 (2008) , Theorem 3.3], who thus obtained the following theorem (in a slightly different form; see [17 D. Daners, J. Glück and J. B. Kennedy, Eventually positive semigroups of linear operators. J. Math. Anal. Appl. 433, 1561–1593 (2016) , Theorem 6.1] for the following version of the result).

Theorem 7.

For a matrix ARd×dA\in\mathbb{R}^{d\times d}, the following assertions are equivalent.

  1. There exists a time t00t_{0}\geq 0 such that all entries of etAe^{tA} are strictly positive for all t>t0t>t_{0}.

  2. The spectral bound s(A)s(A) is a geometrically simple eigenvalue of AA and strictly larger than the real part of every other eigenvalue of AA. Moreover, both AA and ATA^{\mathrm{T}} have a strictly positive eigenvector for s(A)s(A), respectively.

Here, a strictly positive vector means a vector whose entries are all strictly positive.

Individual vs. uniform behaviour

In infinite dimensions, there is a subtlety that we have not properly discussed yet. Let (etA)t[0,)(e^{tA})_{t\in[0,\infty)} be a C0C_{0}-semigroup on a function space EE. If, for every 0u0E0\leq u_{0}\in E, there exists a time t00t_{0}\geq 0 such that etAu00e^{tA}u_{0}\geq 0 for all tt0t\geq t_{0}, it is natural to call the semigroup individually eventually positive since t0t_{0} might depend on u0u_{0}. If in addition t0t_{0} can be chosen to be independent of u0u_{0}, then we call the semigroup uniformly eventually positive.

In finite dimensions, the two concepts can be easily seen to coincide (just apply the semigroup to all canonical unit vectors), but in infinite dimensions, there exist semigroups which are individually but not uniformly eventually positive [17 D. Daners, J. Glück and J. B. Kennedy, Eventually positive semigroups of linear operators. J. Math. Anal. Appl. 433, 1561–1593 (2016) , Examples 5.7 and 5.8].

Conditions for eventual positivity in infinite dimensions

The arguments given in Example 6 show individual eventual positivity of the semigroup, and the same argument can easily be generalised to a more abstract setting. There is one important issue to note, though: if the leading eigenfunction is not bounded away from 00, but might be equal to 00 on the boundary of the underlying domain (as in the case of Dirichlet boundary conditions), then it no longer suffices for the argument that e1AL2e^{1\cdot A}L^{2} be contained in LL^{\infty}; instead, one needs the condition that every vector in e1AL2e^{1\cdot A}L^{2} is dominated by a multiple of the leading eigenfunction. This property is closely related to Sobolev embedding theorems, and can be used to give a characterisation of a certain strong version of individual eventual positivity that is reminiscent of Theorem 7.

On the other hand, giving conditions for uniform rather than individual eventual positivity is more subtle. It requires a domination condition not only on the vectors in the image of e1AL2e^{1\cdot A}L^{2}, but also on the image of the dual operator. If the semigroup is self-adjoint, though, this dual condition becomes redundant and one ends up with the following sufficient condition for uniform eventual positivity.

Theorem 8.

Let (Ω,μ)(\Omega,\mu) be a σ\sigma-finite measure space, let (etA)t[0,)(e^{tA})_{t\in[0,\infty)} be a self-adjoint C0C_{0}-semigroup on L2:L2(Ω,μ)L^{2}\coloneq L^{2}(\Omega,\mu) which leaves the set of real-valued functions invariant, and let uL2u\in L^{2} be a function which is strictly positive almost everywhere. Assume that the following assumptions hold.

  1. The spectral bound s(A)s(A) is a simple eigenvalue of AA, and the corresponding eigenspace contains a function vv satisfying vcuv\geq\nobreak cu for a number c>0c>0.

  2. There exists a time t10t_{1}\geq 0 such that the modulus of every vector in et1AL2e^{t_{1}A}L^{2} is dominated by a multiple of uu.

Then (etA)t[0,)(e^{tA})_{t\in[0,\infty)} is uniformly eventually positive.

The really interesting part in the conclusion of the theorem is the word uniformly, and this is more involved than the argument presented in Example 6. Two different proofs of the theorem are known: the first one is based on an eigenvalue estimate and the theory of Hilbert–Schmidt operators [24 J. Glück, Invariant sets and long time behaviour of operator semigroups. PhD thesis, Universität Ulm (2016) , Theorem 10.2.1] (the assumptions in the reference are slightly different, but the same argument works under the assumptions presented above); the second one employs a duality argument and can thus be generalised to non-self-adjoint semigroups on more general spaces [14 D. Daners and J. Glück, A criterion for the uniform eventual positivity of operator semigroups. Integral Equations Operator Theory 90, Paper No. 46 (2018) , Theorem 3.3 and Corollary 3.5]. This reference also shows that the theorem can be adjusted to even yield a characterisation of a stronger type of eventual positivity.

Theorem 8 implies the non-trivial observation that the semigroup in Example 6 is even uniformly eventually positive.

Spectral properties

Positive semigroups are known to have surprising structural properties, in particular with regard to their spectrum. For some of these properties, it can be shown that they are shared by eventually positive semigroups, though some of the proofs are different from the classical proofs for the positive case. Here are two examples.

  • If the spectrum of the generator AA of an individually eventually positive semigroup is non-empty, then it follows that the spectral bound s(A)s(A) is itself a spectral value [17 D. Daners, J. Glück and J. B. Kennedy, Eventually positive semigroups of linear operators. J. Math. Anal. Appl. 433, 1561–1593 (2016) , Theorem 7.6].

  • For uniformly eventually positive semigroups on LpL^{p}-spaces, the spectral bound s(A)s(A) coincides with the so-called growth bound of the semigroup (see e.g. [20 K.-J. Engel and R. Nagel, One-parameter semigroups for linear evolution equations. Grad. Texts in Math. 194, Springer, New York (2000) , Definition I.5.6] for a definition); this was recently shown by Vogt [35 H. Vogt, Stability of uniformly eventually positive C0-semigroups on Lp-spaces. arXiv:2110.02310v2 (2021) , Theorem 2]. The same can be shown, even for individually eventually positive semigroups, on spaces of continuous functions [6 S. Arora and J. Glück, Stability of (eventually) positive semigroups on spaces of continuous functions. arXiv:2110.04581v1 (2021) , Theorem 4].

More results on the spectrum of eventually positive C0C_{0}-semigroups can be found in [5 S. Arora and J. Glück, Spectrum and convergence of eventually positive operator semigroups. Semigroup Forum 103, 791–811 (2021) ].

4 More examples

The biharmonic heat equation

Example 6 can be adjusted in the following way: we replace the unit interval with a ball BB in Rd\mathbb{R}^{d}, the fourth derivative with the square Δ2\Delta^{2} of the Laplace operator, and the periodic boundary conditions with so-called clamped plate boundary conditions, which require both the function and its normal derivative to vanish at the boundary. On L2(B)L^{2}(B), this yields the operator AA given by

dom(A)=H4(B)H02(B),\displaystyle\operatorname{dom}(A)=H^{4}(B)\cap H^{2}_{0}(B),
Av=Δ2v,\displaystyle Av=-\Delta^{2}v,

where H4(B)H^{4}(B) and H02(B)H^{2}_{0}(B) denote Sobolev spaces. The operator AA is self-adjoint and has negative spectral bound. It thus generates a C0C_{0}-semigroup (etA)t[0,)(e^{tA})_{t\in[0,\infty)} which describes the solutions to the so-called bi-harmonic heat equation

u˙(t)=Au(t)for t[0,).\dot{u}(t)=Au(t)\quad\textrm{for}\ t\in[0,\infty).

We have the following result.

Theorem 9.

The bi-harmonic heat semigroup (etA)t[0,)(e^{tA})_{t\in[0,\infty)} on L2(B)L^{2}(B) is uniformly eventually positive.

Rough outline of the proof. Since BB is a ball, the inverse operator (A)1(-A)^{-1} – or rather its integral kernel, the so-called Green function of AA – can be computed explicitly, and this was in fact done by Boggio over a hundred years ago [9 T. Boggio, Sulle funzioni di Green d’ordine m. Rend. Circ. Mat. Palermo 20, 97–135 (1905) ] (see also [27 H.-C. Grunau and G. Sweers, The maximum principle and positive principal eigenfunctions for polyharmonic equations. In Reaction diffusion systems (Trieste, 1995), Lecture Notes in Pure and Appl. Math. 194, Dekker, New York, 163–182 (1998) , Section 2]). The explicit formula shows that (A)1(-A)^{-1} maps positive functions to positive functions, and even strengthens their positivity in an appropriate sense. Hence, by a Krein–Rutman type result, we obtain that the leading eigenfunction of AA is strictly positive inside BB. Given the specific boundary conditions, it is not too surprising that we also get that assumptions (1) and (2) of Theorem 8 are satisfied if we choose u=d2u=d^{2}, where d ⁣:B[0,)d\colon B\to[0,\infty) describes the distance of each point in BB to the boundary B\partial B. Hence, Theorem 8 gives the desired eventual positivity. ∎

For more details, we refer to [16 D. Daners, J. Glück and J. B. Kennedy, Eventually and asymptotically positive semigroups on Banach lattices. J. Differential Equations 261, 2607–2649 (2016) , second subsection of Section 6] and [14 D. Daners and J. Glück, A criterion for the uniform eventual positivity of operator semigroups. Integral Equations Operator Theory 90, Paper No. 46 (2018) , third subsection of Section 4]. A few comments are in order.

Remark 10. (a) The argument sketched above breaks down for general domains in Rd\mathbb{R}^{d}, since the inverse (A)1(-A)^{-1} need no longer be positive in this case. This is a very well-studied topic in PDE theory; see for instance the surveys [33 G. Sweers, When is the first eigenfunction for the clamped plate equation of fixed sign? In Proceedings of the USA-Chile Workshop on Nonlinear Analysis (Viña del Mar-Valparaiso, 2000), Electron. J. Differ. Equ. Conf. 6, Southwest Texas State Univ., San Marcos, TX, 285–296 (2001) ] by Sweers and [11 A. Dall’Acqua and G. Sweers, On domains for which the clamped plate system is positivity preserving. In Partial differential equations and inverse problems, Contemp. Math. 362, Amer. Math. Soc., Providence, RI, 133–144 (2004) ] by Dall’Acqua and Sweers for more information.

(b) However, if we replace BB with a domain which is sufficiently close to a ball, we still obtain the same result. The main point here is that positivity of (A)1(-A)^{-1} or, under slightly larger perturbations, at least positivity of the leading eigenfunction of AA remains true on such domains as shown by Grunau and Sweers in [27 H.-C. Grunau and G. Sweers, The maximum principle and positive principal eigenfunctions for polyharmonic equations. In Reaction diffusion systems (Trieste, 1995), Lecture Notes in Pure and Appl. Math. 194, Dekker, New York, 163–182 (1998) , Theorem 5.2]. So Theorem 9 holds on this more general class of domains, too.

(c) Theorem 9 remains true on general LpL^{p}-spaces rather than on L2L^{2}; see for instance [14 D. Daners and J. Glück, A criterion for the uniform eventual positivity of operator semigroups. Integral Equations Operator Theory 90, Paper No. 46 (2018) , Theorem 4.4].

(d) If we replace the clamped plate boundary conditions with so-called hinged boundary conditions, which require u=Δu=0u=\Delta u=0 on the boundary, the situation becomes much easier because the operator can then be written as minus the square of the Dirichlet Laplace operator. In this case, we have eventual positivity of the semigroup on general domains; on the space of continuous functions, this example is worked out in [16 D. Daners, J. Glück and J. B. Kennedy, Eventually and asymptotically positive semigroups on Banach lattices. J. Differential Equations 261, 2607–2649 (2016) , Theorem 6.1].

Non-local boundary conditions

Let us now go back to the unit interval and consider the Laplace operator, i.e. the second spatial derivative. If we impose local boundary conditions – such as for instance Dirichlet, Neumann or mixed Dirichlet and Neumann boundary condition, the Laplace operator is well-known to generate a positive semigroup (also on general domains in arbitrary dimension); see for instance [30 E. M. Ouhabaz, Analysis of heat equations on domains. London Math. Soc. Monogr. Ser. 31, Princeton University Press, Princeton, NJ (2005) , Corollary 4.3]. However, let us consider an example of non-local boundary conditions instead. More specifically, we consider the operator AA on L2(0,1)L^{2}(0,1) given by

dom(A)={vH2(0,1):v(0)=v(1)=v(0)+v(1)},\displaystyle\operatorname{dom}(A)=\{v\in H^{2}(0,1):v^{\prime}(0)=-v^{\prime}(1)=v(0)+v(1)\},
Av=v.\displaystyle Av=v^{\prime\prime}.

This is a self-adjoint operator; the operator, and in particular its relation to the Dirichlet and the Neumann Laplace operator, is discussed in more detail in [2 K. Akhlil, Locality and domination of semigroups. Results Math. 73, Paper No. 59 (2018) , Section 3]. The spectral bound of AA is negative, and the inverse (A)1(-A)^{-1} can be computed explicitly [16 D. Daners, J. Glück and J. B. Kennedy, Eventually and asymptotically positive semigroups on Banach lattices. J. Differential Equations 261, 2607–2649 (2016) , proof of Theorem 6.11 (i)]; from this formula and the spectral theory of positive operators, we can conclude that s(A)s(A) is a simple eigenvalue and that there is a corresponding eigenfunction which is strictly positive on the closed interval [0,1][0,1]; see [16 D. Daners, J. Glück and J. B. Kennedy, Eventually and asymptotically positive semigroups on Banach lattices. J. Differential Equations 261, 2607–2649 (2016) , Theorem 6.11] for details. Moreover, we have e1AL2(0,1)dom(A)L(0,1)e^{1\cdot A}L^{2}(0,1)\subseteq\operatorname{dom}(A)\subseteq L^{\infty}(0,1), so the assumptions of Theorem 8 are satisfied for u=1u=\operatorname{\mathbb{1}}, and we obtain the following result.

Theorem 11.

The semigroup (etA)t[0,)(e^{tA})_{t\in[0,\infty)} on L2(0,1)L^{2}(0,1) generated by the Laplace operator with the non-local boundary conditions given above is uniformly eventually positive.

Compare also [4 S. Arora, R. Chill and J.-D. Djida, Domination of semigroups generated by regular forms. arXiv:2111.15489v1 (2021) , Section 4.2] for a related discussion. An example of eventual positivity for different non-local boundary conditions which lead to a non-self-adjoint realisation of the Laplace operator can be found in [14 D. Daners and J. Glück, A criterion for the uniform eventual positivity of operator semigroups. Integral Equations Operator Theory 90, Paper No. 46 (2018) , Theorem 4.3].

Further examples

Today, eventual positivity, and closely related properties as for instance asymptotic positivity, are known for various further C0C_{0}-semigroups, including the semigroup generated by the Dirichlet-to-Neumann operator on the unit circle for various parameter choices [12 D. Daners, Non-positivity of the semigroup generated by the Dirichlet-to-Neumann operator. Positivity 18, 235–256 (2014) ] (which was the initial motivation for the development of the general theory), several delay differential equations ([17 D. Daners, J. Glück and J. B. Kennedy, Eventually positive semigroups of linear operators. J. Math. Anal. Appl. 433, 1561–1593 (2016) , Section 6.5], [24 J. Glück, Invariant sets and long time behaviour of operator semigroups. PhD thesis, Universität Ulm (2016) , Section 11.6] and [14 D. Daners and J. Glück, A criterion for the uniform eventual positivity of operator semigroups. Integral Equations Operator Theory 90, Paper No. 46 (2018) , Theorem 4.6]), the semigroup generated by a bi-Laplacian with certain Wentzell boundary conditions [19 R. Denk, M. Kunze and D. Ploß, The bi-Laplacian with Wentzell boundary conditions on Lipschitz domains. Integral Equations Operator Theory 93, Paper No. 13 (2021) , Section 7], various semigroups on metric graphs ([26 F. Gregorio and D. Mugnolo, Higher-order operators on networks: hyperbolic and parabolic theory. Integral Equations Operator Theory 92, Paper No. 50 (2020) , Proposition 3.7], [25 F. Gregorio and D. Mugnolo, Bi-Laplacians on graphs and networks. J. Evol. Equ. 20, 191–232 (2020) , Section 6] and [8 S. Becker, F. Gregorio and D. Mugnolo, Schrödinger and polyharmonic operators on infinite graphs: parabolic well-posedness and p-independence of spectra. J. Math. Anal. Appl. 495, Paper No. 124748 (2021) , Proposition 5.5]) and semigroups generated by Laplacians coupled by point interactions [28 A. Hussein and D. Mugnolo, Laplacians with point interactions – expected and unexpected spectral properties. In Semigroups of operators – theory and applications, Springer Proc. Math. Stat. 325, Springer, Cham, 47–67 (2020) , Proposition 2].

5 Unbounded domains and local properties

The biharmonic heat equations on unbounded domains

A major drawback of Theorem 8 is that it can only be applied if the leading spectral value is even an eigenvalue of the operator AA. This makes it impossible to apply the theorem to various differential operators that live on unbounded domains. For instance, consider the biharmonic operator AA on L2(Rd)L^{2}(\mathbb{R}^{d}) given by

dom(A)=H4(Rd),\displaystyle\operatorname{dom}(A)=H^{4}(\mathbb{R}^{d}),
Av=Δ2v.\displaystyle Av=-\Delta^{2}v.

The spectrum of AA, which is the set (,0](-\infty,0], does not contain eigenvalues, so Theorem 8 cannot be applied. Still, the semigroup (etA)t[0,)(e^{tA})_{t\in[0,\infty)} exhibits a certain local eventual positivity property: for every compact set KRdK\subseteq\mathbb{R}^{d} and every initial value 0u0L2(Rd)L1(Rd)0\leq u_{0}\in L^{2}(\mathbb{R}^{d})\cap L^{1}(\mathbb{R}^{d}), the solution u ⁣:tetAu0u\colon t\mapsto e^{tA}u_{0} to the biharmonic heat equation u˙(t)=Au(t)\dot{u}(t)=Au(t) becomes eventually positive on KK. This was proved, under slightly different assumptions on u0u_{0} in [23 F. Gazzola and H.-C. Grunau, Eventual local positivity for a biharmonic heat equation in ℝn. Discrete Contin. Dyn. Syst. Ser. S 1, 83–87 (2008) , Theorem 1 (i)] and [22 A. Ferrero, F. Gazzola and H.-C. Grunau, Decay and eventual local positivity for biharmonic parabolic equations. , Theorem 1.1 (ii)] by explicit kernel estimates; for more general powers of Δ\Delta, a similar result was recently shown in [21 L. C. F. Ferreira and V. A. Ferreira, Jr., On the eventual local positivity for polyharmonic heat equations. Proc. Amer. Math. Soc. 147, 4329–4341 (2019) , Theorem 1.1]. Under the assumptions described above, the result was proved by Fourier transform methods in [18 D. Daners, J. Glück and J. Mui, Local uniform convergence and eventual positivity of solutions to biharmonic heat equations. arXiv:2111.02753v1 (2021) , Theorem 2.1].

If we replace the whole space Rd\mathbb{R}^{d} with an infinite cylinder – for instance of the form R×B\mathbb{R}\times B, where BRd1B\subseteq\mathbb{R}^{d-1} is a ball – and again impose clamped plate boundary conditions, the same local eventual positivity result remains true. The proof is technically more involved, though, and relies on a detailed analysis of the specific partial differential equation under consideration; see [18 D. Daners, J. Glück and J. Mui, Local uniform convergence and eventual positivity of solutions to biharmonic heat equations. arXiv:2111.02753v1 (2021) , Theorem 2.3 and Section 4].

However, despite the successful analysis of the aforementioned concrete differential equations, an abstract and general theory as outlined in Section 3 for operators with leading eigenvalue is not yet in sight for the case without eigenvalues.

Open problem. Develop a theory of locally eventually positive C0C_{0}-semigroups (etA)t[0,)(e^{tA})_{t\in[0,\infty)} which is applicable in situations where the generator AA does not have a leading eigenvalue.

Eigenvalues revisited

Getting back to operators which do have a leading eigenvalue, results such as Theorem 8 might still not be applicable in some cases due to conditions (1) and (2) which are sometimes particularly subtle at the boundary of Ω\Omega (if Ω\Omega is, say, a domain in Rd\mathbb{R}^{d} and AA a differential operator). When all functions are restricted to compact subsets of Ω\Omega, though, conditions of the type (1) and (2) might still be satisfied.

This motivates the development of a theory of locally eventually positive semigroups for generators that do have a leading eigenvalue with strictly positive eigenfunction. Such a theory was presented by Arora in [3 S. Arora, Locally eventually positive operator semigroups. J. Oper. Theory, to appear ]. An application of the theory to certain fourth order operators with unbounded coefficients on Rd\mathbb{R}^{d} (which sometimes have eigenvalues due to the growth of the coefficients) was given in [1 D. Addona, F. Gregorio, A. Rhandi and C. Tacelli, Bi-Kolmogorov  type operators and weighted Rellich’s inequalities. arXiv:2104.03811v1 (2021) , Section 3.2].

6 Related topics and results

We close the article by discussing a few related concepts.

Perturbation theory

If AA generates a positive C0C_{0}-semigroup on a function space EE, it is quite easy to see that if BB is a positive and bounded linear operator on EE and MM is a bounded and real-valued multiplication operator on EE, then the perturbed semigroup (et(A+B+M))t[0,)(e^{t(A+B+M)})_{t\in[0,\infty)} is positive, too: if M=0M=0, this follows for instance from the so-called Dyson–Phillips series representation of perturbed semigroups [20 K.-J. Engel and R. Nagel, One-parameter semigroups for linear evolution equations. Grad. Texts in Math. 194, Springer, New York (2000) , Theorem III.1.10], and if MM is non-zero, it follows from the previous case by using the formula

et(A+B+M)=etcet(A+B+M+cid)e^{t(A+B+M)}=e^{-tc}e^{t(A+B+M+c\,\mathrm{id})}

for a real number c0c\geq 0 that is sufficiently large to ensure that M+cidM+c\,\mathrm{id} is positive.

For eventual positivity, though, the situation is much more subtle. Under quite general conditions, one can show that eventual positivity of a semigroup cannot be preserved by all positive perturbations of the generator. This was proved in [15 D. Daners and J. Glück, Towards a perturbation theory for eventually positive semigroups. J. Operator Theory 79, 345–372 (2018) , Theorem 2.3]; related results in finite dimensions had earlier been obtained in [32 F. Shakeri and R. Alizadeh, Nonnegative and eventually positive matrices. Linear Algebra Appl. 519, 19–26 (2017) , Theorem 3.5 and Proposition 3.6]. On the other hand, sufficiently small positive perturbations can be shown not to destroy eventual positivity under appropriate assumptions [15 D. Daners and J. Glück, Towards a perturbation theory for eventually positive semigroups. J. Operator Theory 79, 345–372 (2018) , Section 4].

Maximum and anti-maximum principles

One abstract way to formulate that a linear operator A ⁣:Edom(A)EA\colon E\supseteq\operatorname{dom}(A)\to E on a function space EE satisfies a maximum principle is to require that (A)1(-A)^{-1} be a positive operator, i.e. maps positive functions to positive functions. If 00 is in the spectrum of EE, or more generally if the spectral bound of AA satisfies s(A)0s(A)\geq 0, it is often more natural to consider the resolvent(λidA)1(\lambda\,\mathrm{id}-A)^{-1} for real numbers λ>s(A)\lambda>s(A). If the resolvent at one such point λ0\lambda_{0} is positive, then the same is true for all λ(s(A),λ0)\lambda\in(s(A),\lambda_{0}), too, and we say that AA satisfies a maximum principle. More precisely, this is a uniform maximum principle, while we say that AA satisfies an individual maximum principle if, for each 0fE0\leq f\in E, there exists an (ff-dependent) number λ0>s(A)\lambda_{0}>s(A) such that (λidA)1f0(\lambda\,\mathrm{id}-A)^{-1}f\geq 0 for all λ(s(A),λ0)\lambda\in(s(A),\lambda_{0}).

Similarly, it is common to say that AA satisfies a uniform anti-maximum principle if s(A)s(A) is, say, an isolated spectral value and for all λ\lambda in a left neighbourhood of s(A)s(A) the resolvent (λidA)1(\lambda\,\mathrm{id}-A)^{-1} maps positive functions to negative functions. Likewise, we can define an individual anti-maximum principle (and clearly, the same concepts can be defined at isolated spectral values different from s(A)s(A), too).

Anti-maximum principles have a considerable history and have, for instance, been studied for various elliptic differential operators; see e.g. [10 P. Clément and L. A. Peletier, An anti-maximum principle for second-order elliptic operators. J. Differential Equations 34, 218–229 (1979) ] for a seminal paper on this topic. For biharmonic and polyharmonic operators the validity of (anti-)maximum principles is closely related to the boundary conditions and the geometry of the underlying domain, as explained in Remark 10.

The argument sketched after Theorem 9 can be generalised (and partially reversed) to obtain a correspondence between the following three types of properties:

  1. eventual positivity of the semigroup (etA)t[0,)(e^{tA})_{t\in[0,\infty)},

  2. spectral properties of AA and positivity of the leading eigenfunction,

  3. an individual (anti-)maximum principle for AA.

This correspondence was discussed in [16 D. Daners, J. Glück and J. B. Kennedy, Eventually and asymptotically positive semigroups on Banach lattices. J. Differential Equations 261, 2607–2649 (2016) , Sections 3–5], where the terminology eventual positivity and negativity of the resolvent was used to describe maximum and anti-maximum principles. Indeed, equivalence between the three properties (a)–(c) is true under a number of technical restrictions which have been analysed in more detail in [13 D. Daners and J. Glück, The role of domination and smoothing conditions in the theory of eventually positive semigroups. Bull. Aust. Math. Soc. 96, 286–298 (2017) ].

Uniform (anti-)maximum principles are more difficult to analyse than their individual counterparts – a phenomenon that occurs, as pointed out above, for semigroups, too, but becomes even more pronounced when studying (anti-)maximum principles. An abstract operator theoretic approach to uniform anti-maximum principles was first presented by Takáč in [34 P. Takáč, An abstract form of maximum and anti-maximum principles of Hopf’s type. J. Math. Anal. Appl. 201, 339–364 (1996) , Section 5], and recent progress on the topic was made in [7 S. Arora and J. Glück, An operator theoretic approach to uniform (anti-)maximum principles. J. Differential Equations 310, 164–197 (2022) ]. As a sample result, let us discuss the following special case of [7 S. Arora and J. Glück, An operator theoretic approach to uniform (anti-)maximum principles. J. Differential Equations 310, 164–197 (2022) , Corollary 5.4] for self-adjoint operators on L2L^{2}.

Theorem 12.

Let (Ω,μ)(\Omega,\mu) be a σ\sigma-finite measure space, and let A ⁣:L2dom(A)L2A\colon L^{2}\supseteq\operatorname{dom}(A)\to L^{2} be a real and self-adjoint operator on L2:L2(Ω,μ)L^{2}\coloneq L^{2}(\Omega,\mu). Let uL2u\in L^{2} be a function which is >0>0 almost everywhere, and assume that there exists an integer m0m\geq 0 such that every vector in dom(Am)\operatorname{dom}(A^{m}) is dominated in modulus by a multiple of uu. Assume moreover that λ0R\lambda_{0}\in\mathbb{R} is an isolated spectral value of AA and a simple eigenvalue whose eigenspace contains a function vv that satisfies vcuv\geq cu for a number c>0c>0.

If μ1>λ0\mu_{1}>\lambda_{0} is in the resolvent set of AA and (μ1idA)10(\mu_{1}\,\mathrm{id}-A)^{-1}\geq 0, then the following assertions are equivalent.

  1. One has (μidA)10(\mu\,\mathrm{id}-A)^{-1}\leq 0 for all μ\mu in a left neighbourhood of λ0\lambda_{0}.

  2. There exists a real number d>0d>0 such that

    (μ1idA)1fdf,uufor all 0fL2.(\mu_{1}\,\mathrm{id}-A)^{-1}f\leq d\,\langle f,u\rangle u\quad\textrm{for all}\ 0\leq f\in L^{2}.

The assumption that AA be a real operator means that the domain dom(A)\operatorname{dom}(A) is spanned by real-valued functions and that AA maps real-valued functions to real-valued functions. Assertion (i) of the theorem is a uniform anti-maximum principle, while assertion (ii) can be considered as an upper kernel estimate for the resolvent (in other words: as an upper Green function estimate) of AA. Simple consequences of this theorem are the classical results that the Dirichlet Laplace operator on an interval does not satisfy a uniform anti-maximum principle, while the Neumann Laplace operator on an interval does (see [7 S. Arora and J. Glück, An operator theoretic approach to uniform (anti-)maximum principles. J. Differential Equations 310, 164–197 (2022) , Proposition 6.1 (a) and (b)] for a few more details). More involved examples where the theorem (or more general versions thereof) can be applied are discussed in [7 S. Arora and J. Glück, An operator theoretic approach to uniform (anti-)maximum principles. J. Differential Equations 310, 164–197 (2022) , Section 6].

Jochen Glück is a postdoc at Universität Passau, Germany. His research focusses on order structures in functional analysis, operator semigroups and their applications to differential equations. In 2021 he was awarded the Jaroslav and Barbara Zemánek Prize for achievements in functional analysis, with special emphasis on operator theory. jochen.glueck@uni-passau.de

    References

    1. D. Addona, F. Gregorio, A. Rhandi and C. Tacelli, Bi-Kolmogorov  type operators and weighted Rellich’s inequalities. arXiv:2104.03811v1 (2021)
    2. K. Akhlil, Locality and domination of semigroups. Results Math. 73, Paper No. 59 (2018)
    3. S. Arora, Locally eventually positive operator semigroups. J. Oper. Theory, to appear
    4. S. Arora, R. Chill and J.-D. Djida, Domination of semigroups generated by regular forms. arXiv:2111.15489v1 (2021)
    5. S. Arora and J. Glück, Spectrum and convergence of eventually positive operator semigroups. Semigroup Forum 103, 791–811 (2021)
    6. S. Arora and J. Glück, Stability of (eventually) positive semigroups on spaces of continuous functions. arXiv:2110.04581v1 (2021)
    7. S. Arora and J. Glück, An operator theoretic approach to uniform (anti-)maximum principles. J. Differential Equations 310, 164–197 (2022)
    8. S. Becker, F. Gregorio and D. Mugnolo, Schrödinger and polyharmonic operators on infinite graphs: parabolic well-posedness and p-independence of spectra. J. Math. Anal. Appl. 495, Paper No. 124748 (2021)
    9. T. Boggio, Sulle funzioni di Green d’ordine m. Rend. Circ. Mat. Palermo 20, 97–135 (1905)
    10. P. Clément and L. A. Peletier, An anti-maximum principle for second-order elliptic operators. J. Differential Equations 34, 218–229 (1979)
    11. A. Dall’Acqua and G. Sweers, On domains for which the clamped plate system is positivity preserving. In Partial differential equations and inverse problems, Contemp. Math. 362, Amer. Math. Soc., Providence, RI, 133–144 (2004)
    12. D. Daners, Non-positivity of the semigroup generated by the Dirichlet-to-Neumann operator. Positivity 18, 235–256 (2014)
    13. D. Daners and J. Glück, The role of domination and smoothing conditions in the theory of eventually positive semigroups. Bull. Aust. Math. Soc. 96, 286–298 (2017)
    14. D. Daners and J. Glück, A criterion for the uniform eventual positivity of operator semigroups. Integral Equations Operator Theory 90, Paper No. 46 (2018)
    15. D. Daners and J. Glück, Towards a perturbation theory for eventually positive semigroups. J. Operator Theory 79, 345–372 (2018)
    16. D. Daners, J. Glück and J. B. Kennedy, Eventually and asymptotically positive semigroups on Banach lattices. J. Differential Equations 261, 2607–2649 (2016)
    17. D. Daners, J. Glück and J. B. Kennedy, Eventually positive semigroups of linear operators. J. Math. Anal. Appl. 433, 1561–1593 (2016)
    18. D. Daners, J. Glück and J. Mui, Local uniform convergence and eventual positivity of solutions to biharmonic heat equations. arXiv:2111.02753v1 (2021)
    19. R. Denk, M. Kunze and D. Ploß, The bi-Laplacian with Wentzell boundary conditions on Lipschitz domains. Integral Equations Operator Theory 93, Paper No. 13 (2021)
    20. K.-J. Engel and R. Nagel, One-parameter semigroups for linear evolution equations. Grad. Texts in Math. 194, Springer, New York (2000)
    21. L. C. F. Ferreira and V. A. Ferreira, Jr., On the eventual local positivity for polyharmonic heat equations. Proc. Amer. Math. Soc. 147, 4329–4341 (2019)
    22. A. Ferrero, F. Gazzola and H.-C. Grunau, Decay and eventual local positivity for biharmonic parabolic equations.
    23. F. Gazzola and H.-C. Grunau, Eventual local positivity for a biharmonic heat equation in ℝn. Discrete Contin. Dyn. Syst. Ser. S 1, 83–87 (2008)
    24. J. Glück, Invariant sets and long time behaviour of operator semigroups. PhD thesis, Universität Ulm (2016)
    25. F. Gregorio and D. Mugnolo, Bi-Laplacians on graphs and networks. J. Evol. Equ. 20, 191–232 (2020)
    26. F. Gregorio and D. Mugnolo, Higher-order operators on networks: hyperbolic and parabolic theory. Integral Equations Operator Theory 92, Paper No. 50 (2020)
    27. H.-C. Grunau and G. Sweers, The maximum principle and positive principal eigenfunctions for polyharmonic equations. In Reaction diffusion systems (Trieste, 1995), Lecture Notes in Pure and Appl. Math. 194, Dekker, New York, 163–182 (1998)
    28. A. Hussein and D. Mugnolo, Laplacians with point interactions – expected and unexpected spectral properties. In Semigroups of operators – theory and applications, Springer Proc. Math. Stat. 325, Springer, Cham, 47–67 (2020)
    29. D. Noutsos and M. J. Tsatsomeros, Reachability and holdability of nonnegative states. SIAM J. Matrix Anal. Appl. 30, 700–712 (2008)
    30. E. M. Ouhabaz, Analysis of heat equations on domains. London Math. Soc. Monogr. Ser. 31, Princeton University Press, Princeton, NJ (2005)
    31. A. Pazy, Semigroups of linear operators and applications to partial differential equations. Appl. Math. Sci. 44, Springer, New York (1983)
    32. F. Shakeri and R. Alizadeh, Nonnegative and eventually positive matrices. Linear Algebra Appl. 519, 19–26 (2017)
    33. G. Sweers, When is the first eigenfunction for the clamped plate equation of fixed sign? In Proceedings of the USA-Chile Workshop on Nonlinear Analysis (Viña del Mar-Valparaiso, 2000), Electron. J. Differ. Equ. Conf. 6, Southwest Texas State Univ., San Marcos, TX, 285–296 (2001)
    34. P. Takáč, An abstract form of maximum and anti-maximum principles of Hopf’s type. J. Math. Anal. Appl. 201, 339–364 (1996)
    35. H. Vogt, Stability of uniformly eventually positive C0-semigroups on Lp-spaces. arXiv:2110.02310v2 (2021)

    Cite this article

    Jochen Glück, Evolution equations with eventually positive solutions. Eur. Math. Soc. Mag. 123 (2022), pp. 4–11

    DOI 10.4171/MAG/65
    This open access article is published by EMS Press under a CC BY 4.0 license, with the exception of logos and branding of the European Mathematical Society and EMS Press, and where otherwise noted.