Abstract
We present a generation theorem for positive semigroups on an \(L^1\) space. It provides sufficient conditions for the existence of positive and integrable solutions of initial-boundary value problems. An application to a two-phase cell cycle model is given.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
1 Introduction
We study well-posedness of linear evolution equations on \(L^1\) of the form
where \(\varPsi _0,\varPsi \) are positive and possibly unbounded linear operators on \(L^1\), the linear operator A is such that Eq. (1) with \(\varPsi =0\) generates a positive semigroup on \(L^1\), i.e., a \(C_0\)-semigroup of positive operators on \(L^1\). We present sufficient conditions for the operators \(A, \varPsi _0,\) and \(\varPsi \) under which there is a unique positive semigroup on \(L^1\) providing solutions of the initial-boundary value problem (1). For a general theory of positive semigroups and their applications we refer the reader to [4, 7, 11, 14, 34]. An overview of different approaches used in studying initial-boundary value problems is presented in [13].
Our result is an extension of Greiner’s [19] by considering unbounded \(\varPsi \) and positive semigroups. Unbounded perturbations of the boundary conditions of a generator were studied recently in [1, 2] by using extrapolated spaces and various admissibility conditions. In the proof of our perturbation theorem we apply a result about positive perturbations of resolvent positive operators [3] with non-dense domain in AL-spaces in the form given in [37, Theorem 1.4]. It is an extension of the well known perturbation result due to Desch [15] and by Voigt [41]. For positive perturbations of positive semigroups in the case when the space is not an AL-space we refer to [5, 10]. We also present a result about stationary solutions of (1). We illustrate our general results with an age-size-dependent cell cycle model generalizing the discrete time model of [22, 25, 38]. This model can be described as a piecewise deterministic Markov process (see Sect. 5 and [34]). Our approach can also be used in transport equations [8, 23].
2 General results
Let \((E,\mathcal {E},m)\) and \((E_{\partial },\mathcal {E}_{\partial },m_{\partial })\) be two \(\sigma \)-finite measure spaces. Denote by \(L^1=L^1(E,\mathcal {E},m)\) and \(L^1_{\partial }=L^1(E_{\partial },\mathcal {E}_{\partial },m_{\partial })\) the corresponding spaces of integrable functions. Let \(\mathcal {D}\) be a linear subspace of \( L^1\). We assume that \(A:\mathcal {D}\rightarrow L^1\) and \(\varPsi _0,\varPsi :\mathcal {D}\rightarrow L^1_{\partial }\) are linear operators satisfying the following conditions:
-
(i)
for each \(\lambda >0\), the operator \(\varPsi _0:\mathcal {D}\rightarrow L^1_{\partial }\) restricted to the nullspace \(\mathcal {N}(\lambda I-A)=\{f\in \mathcal {D}:\lambda f-Af=0\}\) of the operator \((\lambda I-A,\mathcal {D})\) has a positive right inverse, i.e., there exists a positive operator \(\varPsi (\lambda ):L^1_{\partial } \rightarrow \mathcal {N}(\lambda I-A)\) such that \(\varPsi _0\varPsi (\lambda )f_{\partial }=f_{\partial }\) for \(f_{\partial }\in L^1_{\partial }\);
-
(ii)
the operator \(\varPsi :\mathcal {D} \rightarrow L^1_{\partial }\) is positive and there exists \(\omega \in {\mathbb {R}}\) such that the operator \(I_{\partial }-\varPsi \varPsi (\lambda ):L^1_{\partial }\rightarrow L^1_{\partial }\) is invertible with positive inverse for all \(\lambda >\omega \), where \(I_{\partial }\) is the identity operator on \(L^1_{\partial }\);
-
(iii)
the operator \(A_0\subseteq A\) with \(\mathcal {D}(A_0)=\{f\in \mathcal {D}:\varPsi _0f=0\}\) is the generator of a positive semigroup on \(L^1\);
-
(iv)
for each nonnegative \(f\in \mathcal {D}\)
$$\begin{aligned} \int _E Af(x)\,m(dx)-\int _{E_\partial }\varPsi _0f(x)\,m_{\partial }(dx)\le 0. \end{aligned}$$(2)
Theorem 1
Assume conditions (i)–(iv). Then the operator \((A_{\varPsi },\mathcal {D}(A_\varPsi ))\) defined by
is the generator of a positive semigroup on \(L^1\). Moreover, the resolvent operator of \(A_{\varPsi }\) at \(\lambda >\omega \) is given by
Proof
The space \(\mathcal {X}=L^1\times L^1_{\partial }\) is an AL-space with norm
We define operators \(\mathcal {A},\mathcal {B}:\mathcal {D}(\mathcal {A}) \rightarrow L^1\times L^1_{\partial }\) with \(\mathcal {D}(\mathcal {A})=\mathcal {D} \times \{0\}\) by (see e.g. [34])
We have \(\mathcal {D}(A_0)\times \{0\}\subset \mathcal {D}(\mathcal {A})\subset L^1\times \{0\}\) and \(\mathcal {D}(A_0)\) is dense in \(L^1\). Hence, \(\overline{\mathcal {D}(\mathcal {A})}=L^1\times \{0\}\). For every \(\lambda >0\) the resolvent of the operator \(\mathcal {A}\) at \(\lambda >0\) is given by
Thus \((\mathcal {A},\mathcal {D}(\mathcal {A}))\) is resolvent positive, i.e., its resolvent operator \(R(\lambda ,\mathcal {A})\) is positive for all sufficiently large \(\lambda >0\). We now show that\(\Vert \lambda R(\lambda , \mathcal {A})\Vert \le 1\) for all \(\lambda >0\). Since the operator \(\lambda R(\lambda , \mathcal {A})\) is positive, it is enough to show that
The operator \(R(\lambda ,A_0)\) is positive, \(R(\lambda ,A_0)f\in \mathcal {D}(A_0)\subseteq \mathcal {D}\) and \( \varPsi _0 R(\lambda ,A_0)f=0\) for \(f\in L^1\). From this and (2) we see that
for all nonnegative \(f\in L^1\). We have \(AR(\lambda ,A_0)f=\lambda R(\lambda ,A_0)f-f\) for all \(f\in L^1\), by (iii). Thus, we get
By assumption (i), \(A\varPsi (\lambda )f_{\partial }=\lambda \varPsi (\lambda )f_{\partial }\) and \(\varPsi _0 \varPsi (\lambda )f_{\partial }=f_{\partial }\) for \(f_{\partial }\in L^1_{\partial }\). This together with condition (2) implies that
for all nonnegative \(f_{\partial }\in L_{\partial }^1\), completing the proof of (6).
Let \(\mathcal {I}\) be the identity operator on \(\mathcal {X}=L^1\times L^1_{\partial }\). We have \(\mathcal {B}R(\lambda ,\mathcal {A})(f,f_{\partial })=(0, \varPsi R(\lambda ,A_0)f+\varPsi \varPsi (\lambda )f_{\partial })\) for any \((f,f_{\partial })\). Thus, \(\mathcal {I}-\mathcal {B}R(\lambda ,\mathcal {A})\) is invertible if and only if \( I_{\partial }-\varPsi \varPsi (\lambda )\) is invertible. In that case
Combining this with (ii) we conclude that \(\mathcal {I}-\mathcal {B}R(\lambda ,\mathcal {A})\) is invertible with positive inverse \((\mathcal {I}-\mathcal {B}R(\lambda ,\mathcal {A}))^{-1}\) for all \(\lambda >\omega \). Hence, the spectral radius of the positive operator \(\mathcal {B}R(\lambda ,\mathcal {A})\) is strictly smaller than 1 for some \(\lambda >\omega \). It follows from [37, Theorem 1.4] that the part of \((\mathcal {A}+\mathcal {B},\mathcal {D}(\mathcal {A}))\) in \(\mathcal {X}_0=\overline{\mathcal {D}(\mathcal {A})}\) denoted by \(((\mathcal {A}+\mathcal {B})_{|},\mathcal {D}((\mathcal {A}+\mathcal {B})_{|}))\) generates a positive semigroup on \(\mathcal {X}_0\). We have \( \mathcal {D}((\mathcal {A}+\mathcal {B})_{|})=\mathcal {D}(A_{\varPsi })\times \{0\}\) and \( (\mathcal {A}+\mathcal {B})_{|}(f,0)=(A_{\varPsi }f,0)\), \(f\in \mathcal {D}(A_{\varPsi }).\) Consequently, the operator \((A_{\varPsi },\mathcal {D}(A_{\varPsi }))\) is densely defined and generates a positive semigroup on \(L^1\). Finally, the operator \((\mathcal {A}+\mathcal {B},\mathcal {D}(\mathcal {A}))\) is resolvent positive with resolvent given by \( R(\lambda ,\mathcal {A}+\mathcal {B})=R(\lambda ,\mathcal {A})(\mathcal {I}-\mathcal {B}R(\lambda ,\mathcal {A}))^{-1} \) for \(\lambda >\omega \). Hence, the formula for \(R(\lambda ,A_{\varPsi })\) is also valid. \(\square \)
Remark 1
Condition (iv) ensures that the operator \((A_0,\mathcal {D}(A_0))\) satisfies
for all nonnegative \(f\in \mathcal {D}(A_0)\). If, additionally,
-
(v)
\((A_0,\mathcal {D}(A_0))\) is densely defined and resolvent positive,
then \((A_0,\mathcal {D}(A_0))\) is the generator of a substochastic semigroup on \(L^1 \), i.e., a positive semigroup of contractions on \(L^1\). This is a consequence of the Hille–Yosida theorem, see e.g. [34, Theorem 4.4]. Thus it is enough to assume condition (v) instead of (iii). Observe also that (iii) and (iv) imply that \((0,\infty )\subseteq \rho (A_0)\).
Remark 2
Note that if \((A_{\varPsi },\mathcal {D}(A_{\varPsi }))\) is the generator of a positive semigroup and
then \((A_{\varPsi },\mathcal {D}(A_{\varPsi }))\) generates a stochastic semigroup, i.e., a positive semigroup of operators preserving the \(L^1\) norm of nonnegative elements (see e.g. [7, Section 6.2] and [34, Corollary 4.1]).
Remark 3
If we assume that
-
(a)
\((A,\mathcal {D})\) is closed,
-
(b)
\(\varPsi _0\) is onto and continuous with respect to the graph norm \(\Vert f\Vert _A=\Vert f\Vert +\Vert Af\Vert \),
then \(\varPsi (\lambda )\) exists for each \(\lambda >0\) and is bounded, by [19, Lemma 1.2]. If \(\varPsi _0\) is positive, then \(\varPsi (\lambda )\) is positive. Thus condition (i) can be replaced by conditions (a) and (b).
Remark 4
Greiner [19, Theorem 2.1] establishes that \((A_{\varPsi },\mathcal {D}(A_{\varPsi }))\) is the generator of a \(C_0\)-semigroup for any bounded \(\varPsi \) provided that conditions (a) and (b) hold true, \((A_0,\mathcal {D}(A_0))\) is the generator of a \(C_0\)-semigroup, and that there exist constants \(\gamma >0\) and \(\lambda _0\) such that
This is condition (2.1) of Greiner [19, Theorem 2.1]. Some extensions of this result are provided in [20, 29] for unbounded \(\varPsi \), as well as in [1, 2].
Remark 5
Recall that a positive operator on an AL-space defined everywhere is automatically bounded. Thus our assumption (i) implies that \(\varPsi (\lambda )\) is bounded for each \(\lambda >0\). Moreover, its norm is determined through its values on the positive cone. From assumptions (i) and (iv) it follows that \(\lambda \Vert \varPsi (\lambda )\Vert \le 1\) for each \(\lambda >0\), as was shown in the proof of Theorem 1. Thus, for \(f=\varPsi (\lambda )f_{\partial }\), we get (9) with \(\gamma =1\). Now suppose, as in [19], that \(\varPsi \) is bounded. Then \(\Vert \varPsi \varPsi (\lambda )\Vert \le \Vert \varPsi \Vert /\lambda \) for all \(\lambda >0\). Hence, the operator \(I_{\partial }-\varPsi \varPsi (\lambda )\) is invertible for \(\lambda > \Vert \varPsi \Vert \). Since \(I-\varPsi (\lambda )\varPsi \) is also invertible, we have \((I-\varPsi (\lambda )\varPsi )^{-1}=I+\varPsi (\lambda )(I_{\partial }-\varPsi \varPsi (\lambda ))^{-1}\varPsi \) and, by (4),
Consequently, if \(\varPsi \) is bounded and positive, then we get the same result as in [19].
We now look at a simple example where Theorem 1 can be easily applied and it should be compared with [1, Corollary 25].
Example 1
Consider the space \(L^1=L^1[0,1]\) and the first derivative operator \(A= \frac{d}{dx}\) with domain \(\mathcal {D}=W^{1,1}[0,1]\). Let \(E_{\partial }\) be the one point set \(\{1\}\) and \(m_{\partial }\) be the point measure \(\delta _1\) at 1, so that the boundary space is \(L^1_{\partial }=\{f_{\partial }:\{1\}\rightarrow {\mathbb {R}}:f_{\partial }(1)\in {\mathbb {R}}\} \) and can be identified with \({\mathbb {R}}\), by writing \(f_{\partial }=f_{\partial }(1)\). Let the boundary operators \(\varPsi _0\) and \(\varPsi \) be defined by
where \(\mu \) is a finite Borel measure. Note that for each \(\lambda >0\) and \(f\in \mathcal {N}(\lambda I-A)\) we have \(f'=\lambda f\). Thus \(f'\) is a continuous function. Consequently, for each \(f_{\partial }\in L^1_{\partial }\) and \(\lambda >0\), the solution \(f=\varPsi (\lambda )f_\partial \) of equation \(f'=\lambda f\) satisfying \(\varPsi _0(\lambda )f=f_{\partial }\) is of the form
Hence condition (i) holds true. We have
and the restriction \(A_0\) of the operator A to
is the generator of a positive semigroup. Thus conditions (iii) and (iv) hold true. If there exists \(\lambda >0\) such that
then condition (ii) holds true and the operator \(A_\varPsi \subseteq \frac{d}{dx}\) with domain
is the generator of a positive semigroup, by Theorem 1. Now suppose that \(\mu \) is a probability measure, so that \(\mu ([0,1])=1\). Then
for all \(\lambda >0\). Thus if (10) does not hold for any \(\lambda >0\) then \(e^{\lambda (x-1)}=1\) for all \(\lambda >0\) and \(\mu \)-almost every \(x\in [0,1]\) implying that \(\mu \{x\in [0,1]: x=1\}=1\). Consequently, if \(\mu \) is a probability measure such that \(\mu \ne \delta _1\) then \((A_\varPsi ,\mathcal {D}(A_\varPsi ))\) is the generator of a positive semigroup.
It should be noted that in [34, Theorem 4.6] the assumption that the domain \(\mathcal {D}(A_{\varPsi })\) of the operator \(A_{\varPsi }\) is dense is missing. Making use of Theorem 1, we get the following result.
Theorem 2
Assume conditions (i)–(iv). If B is a bounded positive operator such that
then \((A_{\varPsi }+B,\mathcal {D}(A_{\varPsi }))\) is the generator of a substochastic semigroup.
We conclude this section with a result concerning the existence of steady states of the positive semigroup from Theorem 1. Note that given any \(\lambda ,\mu \in \rho (A_0)\) we have \( \varPsi (\lambda )=\varPsi (\mu )+(\mu -\lambda )R(\lambda ,A_0)\varPsi (\mu ) \), see [19, Lemma 1.3]. Thus \(\varPsi (\lambda )\ge \varPsi (\mu )\) for \(\lambda \le \mu \). Consequently, for each nonnegative \(f_{\partial }\in L^1_{\partial }\) the pointwise limit
exists and \(\varPsi (0)f_{\partial }\) is nonnegative.
Theorem 3
Assume conditions (i)–(iv). Let \(\varPsi (0)\) be as in (11). If a nonnegative \(f_{\partial }\in L^1_{\partial }\) satisfies \(\varPsi (0)f_{\partial }\in L^1\) and \(f_{\partial }=\varPsi \varPsi (0)f_{\partial }\), then \(\varPsi (0)f_{\partial }\in \mathcal {D}(A_\varPsi )\) and \(A_\varPsi \varPsi (0)f_{\partial }=0\). Conversely, if \(A_\varPsi f=0\) for a nonnegative \(f\in \mathcal {D}(A_\varPsi )\) then \(f_{\partial }=\varPsi f\) satisfies \(\varPsi \varPsi (\lambda )f_{\partial }\le f_{\partial }\) for all \(\lambda >\max \{0,\omega \}\), where \(\omega \) is as in (ii).
Proof
It follows from condition (i) that \(\varPsi (\lambda )f_{\partial }\in \mathcal {D}\), \(\varPsi _0\varPsi (\lambda )f_{\partial }=f_{\partial }\), and \(A\varPsi (\lambda )f_{\partial }=\lambda f_{\partial }\) for all \(\lambda >0\). We have \(\varPsi (\lambda )f_{\partial }\rightarrow \varPsi (0)f_{\partial }\) in \(L^1\), as \(\lambda \rightarrow 0\). Thus \(A\varPsi (\lambda )f_{\partial }\rightarrow 0\) in \(L^1\), as \(\lambda \rightarrow 0\). Recall from the proof of Theorem 1 that the operator \((\mathcal {A}+\mathcal {B})(f,0)=(Af,\varPsi f-\varPsi _0 f)\), \(f\in \mathcal {D}\), is a closed operator in the space \(L^1\times L^1_{\partial }\). The operators \(\varPsi \) and \(\varPsi _0\) are positive and we have \(\varPsi \varPsi (\lambda )f_{\partial }\rightarrow \varPsi \varPsi (0)f_{\partial }=f_{\partial }=\varPsi _0\varPsi (0)f_{\partial }\). Thus, \((\mathcal {A}+\mathcal {B})(\varPsi (\lambda )f_{\partial },0)\rightarrow (0,0)\) as \(\lambda \rightarrow 0\). This implies that \(\varPsi (0)f_{\partial }\in \mathcal {D}(A_{\varPsi })\) and \(A_{\varPsi }\varPsi (0)f_{\partial }=0\).
For the converse, suppose that \(f\in \mathcal {D}(A_\varPsi )\) and \(A_{\varPsi }f=0\). We have \(R(\lambda ,A_\varPsi )(\lambda f-A_\varPsi f)=0\). Thus \(\lambda R(\lambda ,A_\varPsi )f=f\) and \(\varPsi f=\varPsi R(\lambda ,A_\varPsi )(\lambda f)=\varPsi R(\lambda ,A_0)(\lambda f)+\varPsi \varPsi (\lambda ) (I_{\partial }-\varPsi \varPsi (\lambda ))^{-1}\varPsi R(\lambda ,A_0)(\lambda f)\), by (4). Since
we conclude that \(\varPsi f=(I_{\partial }-\varPsi \varPsi (\lambda ))^{-1}\varPsi R(\lambda ,A_0)(\lambda f)\). This implies that \((I_{\partial }-\varPsi \varPsi (\lambda ))\varPsi f=\varPsi R(\lambda ,A_0)(\lambda f)\ge 0\) for \(\lambda >\max \{0,\omega \}\) and completes the proof. \(\square \)
3 A model of a two phase cell cycle in a single cell line
The cell cycle is the period from cell birth to its division into daughter cells. It contains four major phases: \(G_1\) phase (cell growth before DNA replicates), S phase (DNA synthesis and replication), \(G_2\) phase (post DNA replication growth period), and M (mitotic) phase (period of cell division). The Smith–Martin model [36] divides the cell cycle into two phases: A and B. The A phase corresponds to all or part of \(G_1\) phase of the cell cycle and has a variable duration, while the B phase covers the rest of the cell cycle. The cell enters the phase A after birth and waits for some random time \(T_A\) until a critical event occurs that is necessary for cell division. Then the cell enters the phase B which lasts for a finite fixed time \(T_B\). At the end of the B-phase the cell splits into two daughter cells. We assume that individual states of the cell are characterized by age \(a\ge 0\) in each phase and by size \(x> 0\), which can be volume, mass, DNA content or any quantity conserved trough division. We assume that individual cells of size x increase their size over time in the same way, with growth rate g(x) so that \(dx/dt=g(x)\), and all cells age over time with unitary velocity so that \(da/dt=1\). We assume that the probability that a cell is still being in the phase A at age a is equal to \(H(a)\), so the rate of exit from the phase A at age a is \(\rho (a)\) given by
where \(h\) is a probability density function defined on \([0,\infty )\), describing the distribution of the time \(T_A\), the duration of the phase A. We make the following assumptions:
-
(I)
The function \(h\) in (12) is a probability density function so that \(h:[0,\infty )\rightarrow [0,\infty )\) is Borel measurable and the function \(H\) in (12) satisfies: \(H(0)=1\), \(H(\infty )=0\).
-
(II)
The growth rate function \(g:(0,\infty )\rightarrow (0,\infty )\) is globally Lipschitz continuous and \(g(x)>0\) for \(x>0\).
The Smith and Martin hypothesis [36] states that \(h\) is exponentially distributed with parameter \(p>0\), so that \(\rho (a)=p\) for all \(a>0\). However, this does not agree with experimental data, see e.g. [18, 43] for recent results. The generation time of a cell, i.e. the time from birth to division, can be written as \(T=T_{A}+T_B\). Thus the distribution of the generation time has a probability density of the form
Cell generation times can have lognormal or bimodal distribution (see [35]), exponentially modified Gaussian [17], or tempered stable distributions [30].
To describe the growth of cells we define
where \({\bar{x}}> 0\) or \({\bar{x}}=0\), if the integral is finite. The value \(\mathfrak {Q}(x)\) has a simple biological interpretation. If \({\bar{x}}\) is the size of a cell, then \(\mathfrak {Q}(x)\) is the time it takes the cell to reach the size x. It follows from assumption (II) that the function \(\mathfrak {Q}\) is strictly increasing and continuous. We denote by \(\mathfrak {Q}^{-1}\) the inverse of \(\mathfrak {Q}\). Define
for \(t\ge 0\) and \(x_0>0\). Then \(\pi _tx_0\) satisfies the initial value problem
If \(\mathfrak {Q}(0)=-\infty \) then \(\mathfrak {Q}^{-1}\) is defined on \({\mathbb {R}}\). Hence, formula (14) extends to all \(t\in {\mathbb {R}}\) and \(x_0>0\). We also set \(\pi _t0=0\) for \(t>0\) in this case. If \(\mathfrak {Q}(0)=0\) then \(\mathfrak {Q}^{-1}\) is defined only on \((0,\infty )\) and we set \(\pi _t0=\mathfrak {Q}^{-1}(t)\) for \(t>0\). We can extend formula (14) to all negative t satisfying \(\mathfrak {Q}(x_0)+t>0\); otherwise we set \(\pi _{t}x_0=0\). Note that at time \(t=T\), the generation time, a “mother cell” of size \(\pi _{T}x_0\) divides into two daughter cells of equal size \(\frac{1}{2}\pi _{T}x_0\).
In the probabilistic model of [22, 25, 38, 39] a sequence of consecutive descendants of a single cell was studied. Let f be the probability density function of the size distribution at birth at time \(t_0\) of mother cells and let \(t_1>t_0\) be a random time of birth of daughter cells. Then the probability density function of the size distribution of daughter cells is given by [25, 38]
where
The iterates \(P^2f, P^3f,\ldots \) denote densities of the size distribution of consecutive descendants born at random times \(t_2, t_3,\ldots \). The operator P defined by (15) is a positive contraction on \(L^1(0,\infty )\), the space of Borel measurable functions defined on \((0,\infty )\) and integrable with respect to the Lebesgue measure. Here we extend the probabilistic model to a continuous time situation by examining what happens at all times t and not only at \(t_0,t_1,t_2,\ldots \).
We denote by \(p_1(t,a,x)\) and \(p_2(t,a,x)\) the densities of the age and size distribution of cell in the A-phase and in the B-phase at time t, age a, and size x, respectively. Neglecting cell deaths the equations can be written as
with boundary and initial conditions
In this model, cells in the A-phase enter the B-phase at rate \(\rho \). This is taken into account by the boundary condition (18). All cells stay in the B-phase until they reach the age \(T_B\). Then they divide their size into half (17). The model is complemented with initial conditions (19). The model we propose is different as compared to mass/maturity structured models [16, 21, 31, 40] where a cell leaves the phase A with intensity being dependent on maturity, not age. In the case of \(T_B=0\) there is only one phase present; a maturity structured model being a continuous time extension of [24] is studied in [27], while age and volume/maturity structured population models of growth and division were studied extensively since the seminal work of [12, 26, 33]. We refer the reader to [28] for historical remarks concerning modeling of age structured populations and to [35, 42] for recent reviews.
We look for positive solutions of (16)–(19) in the space \(L^1=L^1(E,\mathcal {E},m)\) with \(E=E_1\times \{1\}\cup E_2\times \{2\}\), where
and
m is the product of the two-dimensional Lebesgue measure and the counting measure on \(\{1,2\}\), and \(\mathcal {E}\) is the \(\sigma \)-algebra of all Borel subsets of E. We identify \(L^1=L^1(E,\mathcal {E},m)\) with the product of the spaces \(L^1(E_1)\) and \(L^1(E_2)\) of functions defined on the sets \(E_1\) and \(E_2\), respectively, and being integrable with respect to the two-dimensional Lebesgue measure. We say that the operator P has a steady state in \(L^1(0,\infty )\) if there exists a probability density function f such that \(Pf=f\). Similarly, a semigroup \(\{S(t)\}_{t\ge 0}\) has a steady state in \(L^1\) if there exists a nonnegative \(f\in L^1\) such that \(S(t)f=f\) for all \(t>0\) and \(\Vert f\Vert _1=1\) where \(\Vert \cdot \Vert _1\) is the norm in \(L^1\).
Theorem 4
Assume conditions (I) and (II). There exists a unique positive semigroup \(\{S(t)\}_{t\ge 0}\) on \(L^1\) which provides solutions of (16)–(19) and \(\{S(t)\}_{t\ge 0}\) is stochastic. If \(H\in L^1(0,\infty )\) then the semigroup \(\{S(t)\}_{t\ge 0}\) has a steady state in \(L^1\) if and only if the operator P in (15) has a steady state in \(L^1(0,\infty )\).
We give the proof of Theorem 4 in the next section. Theorem 4 combined with [9] implies the following sufficient conditions for the existence of steady states of (16)–(19).
Corollary 1
Assume conditions (I) and (II). Suppose that \(H\in L^1(0,\infty )\) and that \(|\mathfrak {Q}(0)|<\infty \). If
then (16)–(19) has a steady state and it is unique if, additionally, \(h(a)>0\) for all sufficiently large a. Conversely, if there is \(x_0\ge 0\) such that \(H(Q(\lambda (x_0)))>0\) and \(\mathbb {E}(T_A)>\sup _{x\ge x_0} (\mathfrak {Q}(\lambda (x))-\mathfrak {Q}(x)) \), then (16)–(19) has no steady states.
If the cell growth is exponential so that we have \(g(x)=kx\) for all \(x>0\), where k is a positive constant, then it is known [22, 38, 39] that the operator P has no steady state. We now consider a linear cell growth and assume that \(g(x)=k\) for all \(x>0\). We see that \(\mathfrak {Q}(x)=x/k\), the operator P is of the form (see [39] or the last section)
and condition (20) holds if and only if \(\mathbb {E}(T_A)<\infty \). Combining Corollary 1 with Theorem 4 implies the following.
Corollary 2
Assume that \(g(x)=k\) for \(x>0\) and that \(h(a)>0\) for all sufficiently large \(a>0\). If \(\mathbb {E}(T_A)<\infty \) then the semigroup \(\{S(t)\}_{t\ge 0}\) has a unique steady state.
4 Proof of Theorem 4
We will show that Theorem 4 can be deduced from Theorems 1 and 3. To this end, we introduce some notation. Let us define
where \(\pi _t\) is given by (14). Then \(t\mapsto \pi (t,a_0,x_0)\) solves the system of equations \(a'(t)=1\) and \(x'(t)=g(x(t))\) with initial condition \(a(0)=a_0\) and \(x(0)=x_0\). Recall that \(E_1\) is an open set. For any \(x_0,a_0\in E_1\) we define
and the incoming part of the boundary \(\partial E_1\)
Observe that \(t_{-}(a_0,x_0)=a_0\) for all \((a_0,x_0)\in E_1\) and that \(\varGamma _1^{-}=\{0\}\times (0,\infty )\). We consider on \(\varGamma _1^{-}\) the Borel measure \(m_1^{-}\) being the product of the point measure \(\delta _{0}\) at 0 and the Lebesgue measure on \((0,\infty )\). We define the operator \(\text {T}_{\max }\) on \(L^1(E_1)\) by [6]
with domain
where the differentiation is understood in the sense of distributions. Then it follows from [6] that for \(f\in \mathcal {D}(\text {T}_{\max })\) the following limit
exists for almost every \(z\in \varGamma _1^{-}\) with respect to the measure \(m_1^{-}\) on \(\varGamma _1^{-}\). According to [6, Theorem 4.4] the operator \(T_0=\text {T}_{\max }\) with domain
is the generator of a substochastic semigroup on \(L^1(E_1)\) given by
By [6, Proposition 5.1], the operator \((\text {T},\mathcal {D}(\text {T}))\) defined by
is the generator of a substochastic semigroup on \(L^1(E_1)\) of the from
Note that we can identify the space \(L^1(E_2)\) with the subspace
of \(L^1(E_1)\) and we have \(\text {T}_{\max }(\mathcal {D}(\text {T}_{\max })\cap L^1(E_2))\subseteq L^1(E_2)\). We set
and
We also define the exit time from the set \(E_2\) by
and the outgoing part of the boundary \(\partial E_2\)
We have \(t_+(a_0,x_0)=T_B-a_0\) and \(\varGamma _2^{+}=\{(T_B,x): x>\pi _{T_{B}}0\}\). We define the Borel measure \(m_2^{-}\) on \(\varGamma _2^{-}\) as the measure \(m_1^{-}\) and the \(m_2^{+}\) on \(\varGamma _2^{+}\) as the product of the point measure at \(T_B\) and the one dimensional Lebesgue measure. Since \(U_0(t)(L^1(E_2))\subseteq L^1(E_2)\), the part of the operator \((\text {T}_{0},\mathcal {D}(\text {T}_0))\) in \(L^1(E_2)\) is the generator of a substochastic semigroup \(\{U_2(t)\}_{t\ge 0}\) in \(L^1(E_2)\). Moreover, the following pointwise limits exist
for almost every \(z\in \varGamma _2^{\pm }\) with respect to the Borel measure \(m_2^{\pm }\) on \(\varGamma _2^{\pm }\).
Let \(E_{\partial }= \varGamma _1^{-}\times \{1\}\cup \varGamma _2^{-}\times \{2\}\), \(\mathcal {E}_{\partial }\) be the \(\sigma \)-algebra of Borel subsets of \(E_{\partial }\) and \(m_{\partial }\) be the product of the Lebesgue measure on the line \(\{0\}\times (0,\infty )\) and the counting measure on \(\{1,2\}\). To simplify the notation we identify \(L^1_{\partial }=L^1(E_{\partial },\mathcal {E}_{\partial },m_{\partial })\) with the product space \(L^1(0,\infty )\times L^1(0,\infty )\). We define operators \(A_1\) and \(A_2\) by
We set
and we define the operator A on \(\mathcal {D}\) by setting \(Af=(A_1f_1,A_2f_2)\) for \(f=(f_1,f_2)\in \mathcal {D}\). We take operators \(\varPsi _0,\varPsi :\mathcal {D}\rightarrow L^1_{\partial }\) of the form
and
for \(f=(f_1,f_2)\in \mathcal {D}\). We show that the operator \((A_{\varPsi },\mathcal {D}(A_{\varPsi }))\) is the generator of a positive semigroup on \(L^1\), where \(A_\varPsi f=Af\) for \(f\in \mathcal {D}(A_{\varPsi })=\{f\in \mathcal {D}: \varPsi _0f=\varPsi f\}\). To this end, we check that assumptions (i)–(iv) of Theorem 1 from Sect. 2 are satisfied.
We first show that conditions (iii) and (iv) hold. The operator A restricted to \(\mathcal {D}(A_0)=\{(f_1,f_2)\in \mathcal {D}_1\times \mathcal {D}_2: \text {B}^{-}f_1=0, \text {B}^{-}f_2=0\}\) is the generator of the semigroup \(\{S_0(t)\}_{t\ge 0}\) given by
since \(\{U_1(t)\}_{t\ge 0}\) and \(\{U_2(t)\}_{t\ge 0}\) are semigroups on the spaces \(L^1(E_1)\) and \(L^1(E_2)\) with the corresponding generators. The semigroup \(\{S_0(t)\}_{t\ge 0}\) is substochastic. For all nonnegative \(f=(f_1,f_2)\in \mathcal {D}\) we have
By [6, Proposition 4.6], this reduces to
implying that condition (iv) holds.
For \(f=(f_1,f_2)\in \mathcal {D}\) we can rewrite the equation \(\lambda f-Af=0\) as
Hence, we see that the right inverse of \(\varPsi _0\) when restricted to the nullspace of \(\lambda I-A\) is given by
for \((a,x)\in E_1\) and \(f_{\partial }=(f_{\partial ,1},f_{\partial ,2})\in L^1_{\partial }\). Moreover, if \((f_1,f_2)=\varPsi (\lambda )f_{\partial }\) then
Thus \(f_1\in \mathcal {D}_1\). Similarly, \(f_2\in \mathcal {D}_2\). Hence, condition (i) holds.
To check condition (ii) take \(\lambda >0\) and \(f_{\partial }\in L^1_{\partial }\). For \((f_1,f_2)=\varPsi (\lambda )f_{\partial }\) we have
This implies that
Hence,
where
and, by (12),
For \(f_\partial \in L^1_{\partial }\) we obtain
showing that \(\Vert \varPsi \varPsi (\lambda )\Vert <1\) for all \(\lambda >0\). Consequently, it follows from Theorem 1 that the operator \((A_{\varPsi },\mathcal {D}(A_{\varPsi }))\) is the generator of a positive semigroup \(\{S(t)\}_{t\ge 0}\) on \(L^1\). The semigroup \(\{S(t)\}_{t\ge 0}\) is stochastic, since (8) holds by (25).
Next assume that \(H\in L^1(0,\infty )\). By Theorem 3, it remains to look for fixed points of the operator \(\varPsi \varPsi (0)\). Here \(\varPsi (0)\) defined as in (11) is, by (26), of the form
for \((a,x)\in E_1\). Observe that \(\varPsi (0)f_{\partial }\in L^1\) for \(f_{\partial }\in L^1_{\partial }\), since \(e^{-\int _0^a\rho (r)dr}=H(a)\), by (12), and
We have \(\pi _{-T_B}(2x)=\mathfrak {Q}^{-1}(\mathfrak {Q}(2x)-T_B)=\lambda (x)\) if \(2x>\pi _{T_B}0\) and
Hence
and
If \(f_\partial =\varPsi \varPsi (0)f_{\partial }\) then \(f_{\partial ,2}(x)=(\varPsi \varPsi (0)f_{\partial })_2(x)\) and \(f_{\partial ,1}\) satisfies
By changing the variables \(r=\pi _{-a}(\lambda (x))\), we arrive at the equation
Equivalently, \( f_{\partial ,1} =Pf_{\partial ,1} \) where P is as in (15). Consequently, equation \(\varPsi \varPsi (0)f_{\partial }=f_{\partial }\) has a solution in \(L^1_{\partial }\) if and only if the equation \(Pf_{\partial ,1}=f_{\partial ,1}\) has a solution in \(L^1(0,\infty )\). Observe also that the operator \(\varPsi \varPsi (0)\) preserves the \(L^1_{\partial }\) norm on nonnegative elements. Hence, if \(f_{\partial }\in L^1_{\partial }\) is such that \(\varPsi \varPsi (0)f_{\partial }\le f_{\partial }\) then \(\varPsi \varPsi (0)f_{\partial }= f_{\partial }\). Thus the assertion follows from Theorem 3.
5 Final remarks
Our model can be described as a piecewise deterministic Markov process \(\{X(t)\}_{t\ge 0}\). We considered three variables (a, x, i), where \(i=1\) if a cell is in the phase A, \(i=2\) if it is in the phase B, the variable x describes the cell size, and a describes the time which elapsed since the cell entered the ith phase. Let \(t_0=0\). If we observe consecutive descendants of a given cell and the nth generation time is denoted by \(t_n\), then \(t_{n+1}=s_{n}+T_{B}\) where \(s_n\) is the time when the cell from the nth generation enters the phase B, \(n\ge 0\). A newborn cell at time \(t_n\) is with age \(a(t_n)=0\) and with initial size equal to \(x(t_n^-)/2\), where \(x(t_n^-)\) is the size of its mother cell. The cell ages with velocity 1 and its size grows according to the equations \(x'(t)=g(x(t))\) for \(t\in (t_n,s_n)\). If the cell enters the phase B then its age is reset to 0 and its size still grows according to \(x'(t)=g(x(t))\) for \(t\in (s_n,s_n+T_B)\). We have
and at the end of the second phase the cell divides into two cells, so that we have
Thus the process \(X(t)=(a(t),x(t),i(t))\) satisfies the following system of ordinary differential equations
between consecutive times \(t_0,s_0,t_1,s_1,\ldots \), called jump times. At jump times the process is given by (30) and (31). If the distribution of X(0) has a density f then X(t) has a density S(t)f, i.e.,
for any Borel set \(B_i\subset E_i\), where \(\{S(t)\}_{t\ge 0}\) is the stochastic semigroup from Theorem 4.
If \(f_{\partial ,1}\) is the density of the size distribution at time \(t_0=0\) and \(f_{\partial ,2}\) is the density of the distribution of size at time \(s_1\), then the distribution of size at time \(t_{1}\) is given by
and
where
is the density of the size x(a) of the cell at time a, if x(0) has a density \(f_{\partial ,1}\). Thus the density of the mass \(x(t_{1})\) is given by
for Lebesgue almost every \(x\in (0,\infty )\), where P is as in (15). Now, if the operator P has a steady state \(f_{\partial ,1}\in L^1(0,\infty )\) so that \(f_{\partial ,1}\) satisfies (29) and if \(f_{\partial ,2}\) is as in (32), then \(f^*=(f_1^*,f_2^*)\) given by
is the steady state for the semigroup \(\{S(t)\}_{t\ge 0}\) existing by Theorem 4. Moreover, it is unique if P has a unique steady state.
Remark 6
It should be noted that in the two-phase cell cycle model in [31] the rate of exit from the phase A depends on x, not on a, and that there is no such equivalence between the existence of steady states as presented in Theorem 4. Our results remain true if we assume as in [31] that division into unequal parts takes place. Methods as in [31, 34] can also be used in our model to study asymptotic behaviour of the semigroup \(\{S(t)\}_{t\ge 0}\). For a different approach to study positivity and asymptotic behaviour of solutions of population equations in \(L^1\) we refer to [32].
We conclude this section with an extension of the age-size dependent model from [12] to a model with two phases. Let \(p_i(t,a,x)\) be the function representing the distribution of cells over all individual states a and x at time t in the phase A for \(i=1\) or B for \(i=2\), i.e., \(\int _{a_1}^{a_2}\int _{x_1}^{x_2}p_i(t,a,x)dadx\) is the number of cells with age between \(a_1\) and \(a_2\) and size between \(x_1\) and \(x_2\) at time t in the given phase. Then \(p_1\) and \(p_2\) satisfy Eqs. (16), (18), (19) while the boundary condition (17) takes the form
since a mother cell at the moment of division \(T_B\) has size 2x and gives birth to two daughters of size x entering the phase A at age 0.
Theorem 5
Assume conditions (I) and (II). Then there exists a unique positive semigroup on \(L^1\) which provides solutions of (16), (34), (18), (19).
This follows from Theorem 1 in the same way as Theorem 4, where now to check condition (ii) we note that
for all \(f_\partial \in L^1_{\partial }\) and \(\lambda >0\), implying that \(\Vert \varPsi \varPsi (\lambda )\Vert <1\) for all \(\lambda >\omega \) with \(\omega =\log 2/T_B\).
References
Adler, M., Bombieri, M., Engel, K.J.: On perturbations of generators of \(C_0\)-semigroups. Abstr. Appl. Anal. Art. ID 213020, 13 (2014)
Adler, M., Bombieri, M., Engel, K.J.: Perturbation of analytic semigroups and applications to partial differential equations. J. Evol. Equ. 17(4), 1183–1208 (2017)
Arendt, W.: Resolvent positive operators. Proc. Lond. Math. Soc. (3) 54(2), 321–349 (1987)
Arendt, W., Grabosch, A., Greiner, G., Groh, U., Lotz, H.P., Moustakas, U., Nagel, R., Neubrander, F., Schlotterbeck, U.: One-Parameter Semigroups of Positive Operators. Lecture Notes in Mathematics, vol. 1184. Springer, Berlin (1986)
Arendt, W., Rhandi, A.: Perturbation of positive semigroups. Arch. Math. (Basel) 56(2), 107–119 (1991)
Arlotti, L., Banasiak, J., Lods, B.: On transport equations driven by a non-divergence-free force field. Math. Methods Appl. Sci. 30(17), 2155–2177 (2007)
Banasiak, J., Arlotti, L.: Perturbations of Positive Semigroups with Applications. Springer Monographs in Mathematics. Springer, London (2006)
Banasiak, J., Falkiewicz, A.: Some transport and diffusion processes on networks and their graph realizability. Appl. Math. Lett. 45, 25–30 (2015)
Baron, K., Lasota, A.: Asymptotic properties of Markov operators defined by Volterra type integrals. Ann. Polon. Math. 58(2), 161–175 (1993)
Bátkai, A., Jacob, B., Voigt, J., Wintermayr, J.: Perturbations of positive semigroups on AM-spaces. Semigroup Forum 96(2), 333–347 (2018)
Bátkai, A., Kramar Fijavž, M., Rhandi, A.: Positive operator semigroups: from finite to infinite dimensions. In: Operator Theory: Advances and Applications, vol. 257. Birkhäuser/Springer, Cham (2017)
Bell, G.I., Anderson, E.C.: Cell growth and division: I. A mathematical model with applications to cell volume distributions in mammalian suspension cultures. Biophys. J. 7(4), 329–351 (1967)
Bobrowski, A.: Boundary conditions in evolutionary equations in biology. In: Evolutionary Equations with Applications in Natural Sciences. Lecture Notes in Mathematics, vol. 2126, pp. 47–92. Springer, Cham (2015)
Bobrowski, A.: Convergence of one-parameter operator semigroups. In: Models of Mathematical Biology and Elsewhere. New Mathematical Monographs, vol. 30. Cambridge University Press, Cambridge (2016)
Desch, W.: Perturbations of positive semigroups in AL-spaces (1988, unpublished)
Diekmann, O., Heijmans, H.J.A.M., Thieme, H.R.: On the stability of the cell size distribution. J. Math. Biol. 19(2), 227–248 (1984)
Golubev, A.: Exponentially modified Gaussian (EMG) relevance to distributions related to cell proliferation and differentiation. J. Theor. Biol. 262(2), 257–266 (2010)
Golubev, A.: Applications and implications of the exponentially modified gamma distribution as a model for time variabilities related to cell proliferation and gene expression. J. Theor. Biol. 393, 203–217 (2016)
Greiner, G.: Perturbing the boundary conditions of a generator. Houst. J. Math. 13(2), 213–229 (1987)
Hadd, S., Manzo, R., Rhandi, A.: Unbounded perturbations of the generator domain. Discrete Contin. Dyn. Syst. 35(2), 703–723 (2015)
Hannsgen, K.B., Tyson, J.J.: Stability of the steady-state size distribution in a model of cell growth and division. J. Math. Biol. 22(3), 293–301 (1985)
Hannsgen, K.B., Tyson, J.J., Watson, L.T.: Steady-state size distributions in probabilistic models of the cell division cycle. SIAM J. Appl. Math. 45(4), 523–540 (1985)
Kramar, M., Sikolya, E.: Spectral properties and asymptotic periodicity of flows in networks. Math. Z 249(1), 139–162 (2005)
Lasota, A., Mackey, M.C.: Globally asymptotic properties of proliferating cell populations. J. Math. Biol. 19(1), 43–62 (1984)
Lasota, A., Mackey, M.C., Tyrcha, J.: The statistical dynamics of recurrent biological events. J. Math. Biol. 30(8), 775–800 (1992)
Lebowitz, J.L., Rubinow, S.I.: A theory for the age and generation time distribution of a microbial population. J. Math. Biol. 1(1), 17–36 (1974/1975)
Mackey, M.C., Tyran-Kamińska, M.: Dynamics and density evolution in piecewise deterministic growth processes. Ann. Polon. Math. 94(2), 111–129 (2008)
Metz, J.A., Diekmann, D.: The Dynamics of Physiologically Structured Populations. Lecture Notes in Biomathematics, vol. 68. Springer, Berlin (1986)
Nickel, G.: A new look at boundary perturbations of generators. Electron. J. Differ. Equ. No. 95, 14 (2004)
Palmer, K.J., Ridout, M.S., Morgan, B.J.T.: Modelling cell generation times by using the tempered stable distribution. J. R. Stat. Soc. Ser. C 57(4), 379–397 (2008)
Pichór, K., Rudnicki, R.: Applications of stochastic semigroups to cell cycle models (2018, preprint). arXiv:1806.00091
Rhandi, A.: Positivity and stability for a population equation with diffusion on \(L^1\). Positivity 2(2), 101–113 (1998)
Rubinow, S.: A maturity-time representation for cell populations. Biophys. J. 8(10), 1055–1073 (1968)
Rudnicki, R., Tyran-Kamińska, M.: Piecewise Deterministic Processes in Biological Models. Springer Briefs in Applied Sciences and Technology. Mathematical Methods. Springer, Cham (2017)
Sherer, E., Tocce, E., Hannemann, R., Rundell, A., Ramkrishna, D.: Identification of age-structured models: cell cycle phase transitions. Biotechnol. Bioeng. 99(4), 960–974 (2008)
Smith, J., Martin, L.: Do cells cycle? Proc. Natl. Acad. Sci. USA 70(4), 1263–1267 (1973)
Thieme, H.R.: Positive perturbations of dual and integrated semigroups. Adv. Math. Sci. Appl. 6(2), 445–507 (1996)
Tyrcha, J.: Age-dependent cell cycle models. J. Theor. Biol. 213(1), 89–101 (2001)
Tyson, J.J., Hannsgen, K.B.: Global asymptotic stability of the size distribution in probabilistic models of the cell cycle. J. Math. Biol. 22(1), 61–68 (1985)
Tyson, J.J., Hannsgen, K.B.: Cell growth and division: a deterministic/probabilistic model of the cell cycle. J. Math. Biol. 23(2), 231–246 (1986)
Voigt, J.: On resolvent positive operators and positive \(C_0\)-semigroups on \(AL\)-spaces. Semigroup Forum 38(2), 263–266 (1989)
Webb, G.F.: Population models structured by age, size, and spatial position. In: Structured Population Models in Biology and Epidemiology. Lecture Notes in Mathematics, vol. 1936, pp. 1–49. Springer, Berlin (2008)
Yates, C.A., Ford, M.J., Mort, R.L.: A multi-stage representation of cell proliferation as a Markov process. Bull. Math. Biol. 79(12), 2905–2928 (2017)
Author information
Authors and Affiliations
Corresponding author
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
This work was partially supported by the Polish NCN Grant 2017/27/B/ST1/00100 and by the Grant 346300 for IMPAN from the Simons Foundation and the matching 2015–2019 Polish MNiSW fund.
Rights and permissions
OpenAccess This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
About this article
Cite this article
Gwiżdż, P., Tyran-Kamińska, M. Positive semigroups and perturbations of boundary conditions. Positivity 23, 921–939 (2019). https://doi.org/10.1007/s11117-019-00644-w
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11117-019-00644-w