Abstract
We study the existence and uniqueness of the solution of a non-linear coupled system constituted of a degenerate diffusion-growth-fragmentation equation and a differential equation, resulting from the modeling of bacterial growth in a chemostat. This system is derived, in a large population approximation, from a stochastic individual-based model where each individual is characterized by a non-negative trait whose dynamics is described by a diffusion process. Two uniqueness results are highlighted. They differ in their hypotheses related to the influence of the resource on individual trait dynamics, the main difficulty being the non-linearity due to this dependence and the degeneracy of the diffusion coefficient. Further we show by probabilistic arguments that the semi-group of the stochastic trait dynamics admits a density. We deduce that the diffusion-growth-fragmentation equation admits a function solution with a certain Besov regularity.
Similar content being viewed by others
References
Aldous, D.: Stopping times and tightness. Ann. Probab. 6, 335–340 (1978)
Bansaye, V., Méléard, S.: Stochastic Models for Structured Populations: Scaling Limits and Long Time Behavior. Springer (2015)
Bansaye, V., Tran, V.C.: Branching feller diffusion for cell division with parasite infection. ALEA Lat. Am. J. Probab. Math. Stat. 8, 95–127 (2011)
Bertoin, J.: The asymptotic behavior of fragmentation processes. J. Eur. Math. Soc. (JEMS) 5(4), 395–416 (2003)
Bertoin, J.: On a Feynman-Kac approach to growth-fragmentation semigroups and their asymptotic behaviors. J. Funct. Anal. 277(11), 108270 (2019)
Bertoin, J., Watson, A.R.: Probabilistic aspects of critical growth-fragmentation equations. Adv. Appl. Prob. 48, 37–61 (2016)
Bertoin, J., Watson, A.R.: A probabilistic approach to spectral analysis of growth-fragmentation equations. J. Funct. Anal. 274(8), 2163–2204 (2018)
Bouleau, N., Hirsch, F.: Dirichlet Forms and Analysis on Wiener Space. de Gruyter Studies in Mathematics, vol. 14 (1991)
Campillo, F., Fritsch, C.: Weak convergence of a mass-structured individual-based model. Appl. Math. Optim. 72, 37–73 (2015)
Champagnat, N., Jabin, P.-E., Méléard, S.: Adaptation in a stochastic multi-resources chemostat model. Journal de Mathématiques Pures et Appliquées 101(6), 755–788 (2014)
Champagnat, N., Méléard, S.: Invasion and adaptative evolution for individual-based spatially structured populations. J. Math. Biol. 55, 147–188 (2007)
Collet, P., Martínez, S., Méléard, S.: Stochastic models for a chemostat and long-time behavior. Adv. Appl. Prob. 45, 822–836 (2013)
Dȩbiec, T., Doumic, M., Gwiazda, P., Wiedemann, E.: Relative entropy method for measure solutions of the growth-fragmentation equation. SIAM J. Math. Anal. 50(6), 5811–5824 (2018)
Debussche, A., Romito, M.: Existence of densities for the 3D Navier–Stokes equations driven by Gaussian noise. Probab. Theory Relat. Fields 158, 575–596 (2014)
Doumic, M., Escobedo, M.: Time asymptotics for a critical case in fragmentation and growth-fragmentation equations. Kinet. Relat. Models 9(2), 251–297 (2016)
Fontbona, J., Méléard, S.: Non local Lotka-Volterra system with cross-diffusion in an heterogeneous medium. J. Math. Biol. 70(4), 829–854 (2015)
Fournier, N., Méléard, S.: A microscopic probabilistic description of a locally regulated population and macroscopic approximation. Ann. Appl. Probab. 14(4), 1880–1919 (2004)
Friedman, A.: Differential Equations and Applications, vol. 1. Academic Press, Inc. (LONDON) LTD. (1975)
Harmand, J., Lobry, C., Rapaport, A., Sari, T.: Le Chémostat: Théorie mathématique de la culture de micro-organismes, vol. 1 (2017)
Ikeda, N., Watanabe, S.: Stochastic Differential Equations and Diffusion Processes. North-Holland Mathematical Library (1988)
Jauffret, M.D., Gabriel, P.: Eigen elements of a general aggregation-fragmentation model. Math. Models Methods Appl. Sci. 20(5), 757–783 (2010)
Jourdain, B., Méléard, S., Woyczynski, W.A.: Lévy flights in evolutionary ecology. J. Math. Biol. 65, 677–707 (2012)
Kunita, H.: Stochastic differential equations and stochastic flows of diffeomorphisms. In: ’École d’été de probabilités de Saint-Flour, XII-1982’. Lectures Notes in Mathematics, vol. 1097, pp. 143–303. Springer, Berlin (1984)
Lions, P.-L., Musiela, M.: Convexity of solutions of parabolic equations. C. R. Acad. Sci. Paris 342, 915–921 (2006)
Méléard, S., Roelly, S.: Sur les convergences étroite ou vague de processus à valeurs mesures. C. R. Acad. Sci. Paris Sér. I Math. 317, 785–788 (1993)
Monod, J.: La technique de culture continue, théorie et applications. Ann. Inst. Pasteur 79(4), 390–410 (1950)
Novick, A., Szilard, L.: Description of the chemostat. Science 112(2920), 715–716 (1950)
Revuz, D., Yor, M.: Continuous Martingales and Brownian Motion. Springer (1998)
Roelly-Coppoletta, S.: A criterion of convergence of measure-valued processes: application to measure branching processes. Stoch. Stoch. Rep. 17, 43–65 (1986)
Romito, M.: A simple method for the existence of a density for stochastic evolutions with rough coefficients. Electron. J. Probab. 23, 1–43 (2018)
Sato, K., Ueno, T.: Multi-dimensional diffusion and the Markov process on the boundary. J. Math. Kyoto Univ. 4(3), 529–605 (1965)
Tran, V.C.: Large population limit and time behaviour of a stochastic particle model describing an age-structured population. ESAIM Probab. Stat. EDP Sci. 12, 345–386 (2008)
Triebel, H.: Theory of Function Spaces. Monographs in Mathematics, vol. 78. Birkhäuser Verlag, Basel (1983)
Triebel, H.: Theory of Function Spaces II. Monographs in Mathematics, vol. 84. Birkhäuser Verlag, Basel (1992)
Yamada, T., Watanabe, S.: On the uniqueness of solutions of stochastic differential equations. J. Math. Kyoto Univ. (JMKYAZ) 11(1), 553–563 (1971)
Acknowledgements
I would like to thank Sylvie Méléard and Carl Graham for their continual guidance during this work. I also thank Sylvie Méléard for fruitful exchanges and discussions. I will end by thanking an anonymous reviewer for his relevant remarks that allowed me to perform the results and the quality of this article. This work has been supported by the Chair “Modélisation Mathématique et Biodiversité” of Veolia-École Polytechnique-Muséum national d’Histoire naturelle-Fondation X.
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.
Appendices
Appendix
Proof of Theorem 2.1
Let us recall that Assumptions 2.1 and 2.2 hold.
1.1 Step 1: Control of tails
Since the individual trait is unbounded, we need to get a control on the measure valued process tails. We then have the following lemma that is adapted from [22].
Lemma A.1
Let \(T>0\), if a subsequence of \(\left\{ \left( \nu ^K_t\right) _{t\ge 0},K>1 \right\} \) converges in law in the Skorohod space \(\mathbb {D}([0,T],(\mathcal {M}_F(\mathbb {R}_+),v))\) then its limit \((\nu _t)_{t\in [0,T]}\) satisfies
In addition, there exists a decreasing sequence of increasing functions \((f_n)_{n\ge 1}\subset \mathcal {C}_b^2(\mathbb {R}_+,\mathbb {R})\) such that \(f_n=0\) on [0, n/2] and \(f_n = 1\) on \([n,\infty )\), which satisfies
Proof of Lemma A.1
In order to show (A.1), we approximate the functions \(x\mapsto 1\) by an increasing sequence of continuous and compact supported functions in order to use the convergence assumption for the vague topology. We then introduce the uniformly bounded sequence of continuous functions defined for \(n\in \mathbb {N}^{\star }\) by
that converges increasingly pointwise towards the constant function \(x\mapsto 1\). It follows from the Fatou lemma that holds thanks to the convergence of the sequence of measure valued process for the vague topology, that we have for all \(n\in \mathbb {N}\)
As \(n\rightarrow \infty \) in this inequality, we obtain by monotony
In order to show (A.2), we approximate the indicator of the outside of the compacts by a non-increasing and uniformly bounded sequence of regular functions. This allows us to control the tails of the distributions \(\{\mathcal {Z}^K,K>1\}\). Let us then introduce an increasing function \(f\in \mathcal {C}^2(\mathbb {R}_+,\mathbb {R})\) such that \(f\equiv 0\) on [0, 1/2] and \(f\equiv 1\) on \([1,\infty )\). We set
which belongs to \(\mathcal {C}^2_b(\mathbb {R}_+,\mathbb {R})\) and whose derivatives have compact supports. It then follows from (2.12) that for all \(s\le t\le T\)
where
The sequence \((A^{K,f_n}_T)_{K>1}\) is uniformly integrable according to K and converges in law towards the random variable \(C\int _0^T \left\langle \nu _u,(1+x)(f'_n + \left| f''_n\right| )\right\rangle du\) as \(K\rightarrow \infty \) because the function \(x\mapsto (1+x)(f'_n + \left| f''_n\right| )(x)\) is continuous with compact support. We have
that induces thanks to Gronwall’s lemma that
As \(K\rightarrow \infty \), it follows from (2.15) and the Doob inequality that the right term with the martingale converges towards 0. In addition the remaining terms in the expectation are uniformly integrable and converge in law, then
\(\square \)
1.2 Step 2: Tightness
Given the above Lemma A.1, the tightness of the sequence of laws \(\Big \{\mathcal {Z}^K = \mathcal {L}(\nu ^K,R^K),K>1 \Big \}\) follows. It is adapted from [2], Theorem 7.4.
Proposition A.1
The sequence of laws \(\left\{ \mathcal {Z}^K = \mathcal {L}(\nu ^K,R^K),K>1 \right\} \) is tight in the space of probability measures \(\mathcal {P}(\mathbb {D}([0,T],(\mathcal {M}_F(\mathbb {R}_+)\times \mathbb {R},w\otimes |.|)))\) for all \(T>0\).
Proof
Let us first endow the measure space \(\mathcal {M}_F(\mathbb {R}_+)\) with its vague topology and prove the tightness property for the sequence \(\left\{ \mathcal {Z}^K=\mathcal {L}\left( \nu ^K, R^K\right) ,K>1 \right\} \). By noticing that \(\mathcal {C}^2_c(\mathbb {R}_+,\mathbb {R})\) is a dense subspace of \(\mathcal {C}_0(\mathbb {R}_+,\mathbb {R})\) for the topology of uniform convergence on compact sets, it suffices according to [29] to show that for each test function \(f\in \mathcal {C}_c^2(\mathbb {R}_+,\mathbb {R})\) the sequence of processes \(\left\{ \left( \left\langle \nu ^K_t,f\right\rangle , R^K_t\right) _{t\ge 0},K>1 \right\} \) is tight in the Skorohod space \(\mathbb {D}([0,T],\mathbb {R}^2)\). First, remember that we have the decomposition
(see Lemma 2.1) and thanks to (2.15),
It then directly follows from Doob’s inequality that the sequence of martingales \(\left\{ M^{K,f},K>1 \right\} \) converges in \(\mathbb {L}^2\) toward 0 locally uniformly in time. This sequence is then tight in the Skorohod space \(\mathbb {D}([0,T],\mathbb {R})\), and all we need to prove now is the tightness of the other sequences \(\{V^{K,f},K>1\}\) and \(\{R^K,K>1\}\). For that purpose, we use the Aldous criterion (see [1]), that consists in showing the uniform controls
and showing that the Aldous condition holds for \((R^K_t,V^{K,f}_t)_{t\ge 0}\). The assertions (a) and (b) are immediate. Indeed, \(R^K_t \le r\vee R_0 = \bar{R}\) for any \(t\ge 0\) and thanks to the boundedness of the test function,
that implies (a) thanks to (2.10) and (2.15). Let us now consider \(\varepsilon >0\) and a stopping time \(\tau \) that satisfies \(\tau +\epsilon \le T\), then
By a similar computation, we obtain
that implies the Aldous condition thanks to Lemma 2.1. The sequence of law \((\mathcal {Z}^K)_K\) is then tight and by the Prokhorov’s theorem one can extract from each sub-sequence, a sub-sub-sequence that converges in \(\mathcal {P}(\mathbb {D}([0,T],(\mathcal {M}_F(\mathbb {R}_+)\times \mathbb {R},v\otimes |.|)))\). Let us denote by \(\mathcal {Z}\) a limit value of this sequence, and for simplicity of the notation, \((\mathcal {Z}^K)_K\) a sub-sequence that converges towards \(\mathcal {Z}\). For each \(K>1\) we have
We then deduce from the continuity of the mapping \(\mu \mapsto \sup _{0\le t\le T}\left| \left\langle \mu _{t},f\right\rangle - \left\langle \mu _{t-},f\right\rangle \right| \) for \(f\in \mathcal {C}_c(\mathcal {E}\times \mathbb {R}_+,\mathbb {R})\), that the limiting law \(\mathcal {Z}\) only charges the set \(\mathcal {C}([0,T],(\mathcal {M}_F(\mathbb {R}_+)\times \mathbb {R},v\otimes |.|))\). We are aiming to prove the same result for the weak topology on \(\mathcal {M}_F(\mathbb {R}_+)\) and this is where the sequence of functions introduced in Lemma A.1 is usefull.
Let us denote by \((\nu ,R)\) a process of law \(\mathcal {Z}\), then according to the above argument its trajectories are in the space \(\mathcal {C}([0,T],(\mathcal {M}_F(\mathbb {R}_+)\times \mathbb {R},v\otimes |.|))\). As in [22], it follows from the Fatou lemma that for \(n,l\in \mathbb {N}\),
As \(l\rightarrow \infty \) it follows from Beppo-Levi’s theorem applied to the left term, that
and then by the Lemma A.1,
One can therefore extract a sub-sequence \(\left( \sup _{0\le t\le T}\langle \nu _t,f_{n_k}\rangle \right) _k\) that converges almost surely towards 0, and then the limit law \(\mathcal {Z}\) only charges the set \(\mathcal {C}([0,T],(\mathcal {M}_F(\mathbb {R}_+)\times \mathbb {R},w\otimes |.|))\). We deduce that the convergence of \(\{\mathcal {Z}^K,K>1\}\) towards \(\mathcal {Z}\) for the weak topology on \(\mathcal {M}_F(\mathbb {R}_+)\) holds if the sequence of total masses \(\left\{ \left\langle \nu ^K,1\right\rangle ,K>1\right\} \) converges in law towards \(\left\langle \nu ,1\right\rangle \) in the Skorohod space \(\mathbb {D}([0,T],\mathbb {R})\) (see [25]). So let us consider a Lipschitz continuous function \(\Phi : \mathbb {D}([0,T],\mathbb {R})\rightarrow \mathbb {R}\), and \(n\in \mathbb {N}\).
The function \(1-f_n\) being continuous with a compact support and the sequence \(\{\nu ^K,K>1\}\) converging in law toward the limit process \(\nu \) in \(\mathbb {D}([0,T],(\mathcal {M}_F(\mathbb {R}_+)\times \mathbb {R},v\otimes |.|))\), we have
It follows that for all \(n\in \mathbb {N}\),
Thanks to Lemma A.1, that implies that as \(n\rightarrow \infty \) in the right term,
The sub-sequence \(\{(\langle \nu ^K_t,1\rangle )_{t\in [0,T]},K>1 \}\) then converges in law towards the process \((\langle \nu _t,1\rangle )_{t\in [0,T]}\), and then the subsequence \((\mathcal {Z}^K)_K\) converges towards \(\mathcal {Z}\) in the space of probability distributions \(\mathcal {P}(\mathbb {D}([0,T],(\mathcal {M}_F(\mathbb {R}_+)\times \mathbb {R},w\otimes |.|)))\). That is enough to conclude that our initial sequence of laws is tight. \(\square \)
1.3 Step 3: Identification of a limiting process
We are here interested in the characterization of a process \((\nu _t,R_t)_{t\in [0,T]}\) whose law is a limiting value of the sequence \(\{\mathcal {Z}^K,K>1\}\). Given the operator defined in (2.9), we have the following result
Proposition A.2
Let us denote by \(\mathcal {Z}\) a limiting value of the sequence of distributions \(\left\{ \mathcal {Z}^K,K>1\right\} \). Then a process \((\nu ,R)\) of law \(\mathcal {Z}\) satisfies (2.16).
Proof
Let \(f\in \mathcal {C}_c^2(\mathbb {R}_+,\mathbb {R})\) be a test function and \(t\in [0,T]\) a fixed time. We introduce the following real-valued functions defined for \((\mu ,r)\) in the subset of \(\mathcal {C}([0,T],(\mathcal {M}_F(\mathbb {R}_+)\times \mathbb {R},w\otimes |.|))\) constitued of processes that satisfy \(|r| \le \bar{R}\)
We are aiming to show that
For this purpose, we first note that the function \(\Phi _{t}\) defined as above is continuous and
It then follows by the Fatou lemma that
Futhermore, the terms in the integrals in the expression of the function \(\Phi ^f_t\) are continuous and bounded, then this function is continuous. In addition, it follows from the decomposition (2.12) that
and then thanks to the Fatou lemma and Doob’s inequality,
We now deduce from (2.15) that
\(\square \)
Some properties of a specific SDE
Let us introduce two continuous functions \(h(x,s),q(x,s)\in \mathcal {C}(\mathbb {R}_+\times [0,T],\mathbb {R})\) that satisfy the following assumption.
Assumption B.1
-
(A.1) Boundary condition: for all \(s\le [0, T]\),
$$\begin{aligned} q(0,s) = 0 \text { and }h(0,s) \ge 0. \end{aligned}$$ -
(A.2) Lipschitz condition: for all \(x,y\ge 0\) and \(s\in [0,T]\),
$$\begin{aligned} |h(x,s) - h(y,s)| + |q(x,s) - q(y,s)| \le C|x-y|. \end{aligned}$$ -
(A.3) The function q(x, s) is non negative.
Then we have the following results
Lemma B.1
There is weak existence and uniqueness, and strong existence for the stochastic differential equation
The solution \((Z^x_{s,u})_{u\in [s,T]}\) that starts at time \(s<T\) at the position \(x\ge 0\), has a zero local time at 0, is almost surely non negative at any time and satifies the moment estimate
In addition, for all \(x_1,x_2\ge 0\) and \(s_1,s_2\le T\), we have for all \(u\le T\)
where \((Z^{x_1}_{s_1,u})_{u\in [s_1,T]}, (Z^{x_2}_{s_2,u})_{u\in [s_2,T]}\) are assumed to be built with the same Brownian motion, and satisfies \(Z^{x_j}_{s_j,u} = x_j\) for all \(u\le s_j\) with \(j=1,2\).
Lemma B.2
Let \(p(x)\in \mathcal {C}^1_b(\mathbb {R}_+,\mathbb {R})\) and \(a(x,s)\in \mathcal {C}_b(\mathbb {R}_+\times [0,T],\mathbb {R})\) such that
then for a given \(t\le T\), the function defined by
satisfy for all \(s_1,s_2\le t\) and \(x_1,x_2\ge 0\) the following inequality
1.1 Proof of Lemma B.1
We split this proof into the several steps.
Existence and uniqueness According to [20], Theorem 1.1 it suffices to show that we have weak existence and pathwise uniqueness to conclude that there is weak existence and uniqueness, and strong existence. In addition it follows from the previous paragraph that this solution will be well defined at any time. The stochastic differential equation (B.1) is non homogeneous with continuous coefficients, then the coupled process \((Z_u,u)_{u\ge 0}\) satisfies an homogeneous stochastic differential equation with continuous coefficients. We then have weak existence for each initial condition \(x\ge 0\) and initial time \(s\in [0,T]\) thanks to [20], Theorem 2.3, p173. Furthermore, for \(x,y\ge 0\) and \(s\le u\le T\) we have by hypothesis
thanks to the inequality \(|\sqrt{a}-\sqrt{b}|\le \sqrt{|a-b|}\,,\, \forall a,b\ge 0\). We deduce from ( [35]) that there is pathwise uniqueness.
Moment estimate and positiveness Notice that the solution that we constructed in the previous paragraph is well defined until it becomes negative or explodes. In order to get well defined terms for all times, we introduce the new equation
for which there is weak existence and pathwise uniqueness by an argument similar to the one in the previous paragraph, but only until explosion time. It follows from the continuity of the function h(x, s) and hypothesis (A.2) of Assumption B.1 that
For given \(x\ge 0\) and \(s\le t\), we consider the sequence of stopping times defined by
that converges increasingly. Then for \(p\ge 2\) and \(v\in [s,T]\) we have
that implies that
It follows from the Doob inequality that
and then
Hence, we obtain for all \(v\in [s,T]\)
and thanks to the Gronwall lemma,
If we assume that there exists \(T_{0}<T\) such that \(\mathbb {P}\left\{ \lim _N\lambda _N< T_0 \right\} > 0\), then for all \(v\in [T_0,T]\), we will have
That is impossible because the left side term tends towards infinity as \(N\rightarrow \infty \). It follows that \(\lambda _N\xrightarrow []{a.s} T\) and by the monotone convergence theorem,
It then suffices to take \(v = T\). The case \(0< p< 2\) is obviously solved thanks to the Hölder inequality
The process \((\tilde{Z}^x_{s,u})_{u\in [s,T]}\) is then well defined at any time and satisfies a moment estimate similar to (B.2). We will now show that this process is almost surely non negative, and then by uniqueness this will imply that it is a modification of the process \((Z^x_{s,u})_{u\in [s,T]}\). Indeed, the following decomposition holds
and it follows from the Tanaka formula that
For \(\varepsilon >0\),
since \((x,u)\mapsto q(x,u)/x\) is bounded. Then the local time \(L^0_{t}(\tilde{Z}^x_{s,\cdot })\) is zero almost surely (see [28], Ch. IX, Lemma 3.3). Equation (B.8) becomes
where the second term in the right side of the equality is a true martingale thanks to the moment estimate that we previously stated. It then follows from the hypothesis (A.1) and (A.2) in Assumption B.1 that
and hence by the Gronwall lemma, \(\mathbb {E}\left[ (\tilde{Z}^x_{s,t})_- \right] = 0\) for any \(t\in [s,T]\), thus \((\tilde{Z}^x_{s,t})_- = 0\) p.s. We deduce that the stochastic process \((\tilde{Z}^x_{s,t})_{t\in [s,T]}\) is almost surely non negative at any time, and it is possible to remove the absolute values in (B.6) to conclude by a uniqueness argument that
The stochastic process \((Z^x_{s,t})_{t\in [s,T]}\) is then almost surely non negative at any time, satisfies the moment estimate (B.2) and has a zero local time at 0.
Dependence on initial conditions Thanks to the strong existence, it is possible to construct two solutions \((Z^{x_1}_{s_1,v})_{v\in [s_1,T]}\) and \((Z^{x_2}_{s_2,v})_{v\in [s_2,T]}\) of (B.1) on the same filtered probability space and with the same Brownian motion, for all \(x_1,x_2\ge 0\) and \(0\le s_1\le s_2\le T\). Those solutions are well extended on [0, T] if we set
Hence, it follows from the Tanaka formula that for all \(s_2\le v\le T\),
where \(l^0_v(Z^{x_1}_{s_1,\cdot }-Z^{x_2}_{s_2,\cdot })\) is the local time of the process \(Z^{x_1}_{s_1,\cdot }-Z^{x_2}_{s_2,\cdot }\) at 0 up to time v. Furthermore, thanks to the inequality \(|\sqrt{a}-\sqrt{b}|\le \sqrt{|a-b|}\,,\forall a,b\ge 0\) again, we have for \(\varepsilon >0\),
We deduce thanks to [28], Ch. IX, lemma 3.3 that the local time \(l^0_{v}(Z^{x_1}_{s_1,\cdot }-Z^{x_2}_{s_2,\cdot })\) is zero almost surely for all \(v\in [s_2,T]\), and it follows from (B.9) that
which implies by the Gronwall lemma that
We are now interested in what happens for \(s_1\le v\le s_2\). Indeed,
and thanks to the Doob inequality,
Hence, we have for all \(v\in [s_1,s_2]\),
which implies thanks to (B.2) that
Notice that this inequality gives (B.3) for \(v\in [s_1,s_2]\). Taking in particular \(v=s_2\) and considering (B.10), we deduce that (B.3) also holds for \(v\in [s_2,T]\). We finally obtain without any order on \(s_1,s_2\in [0,T], x_1,x_2\ge 0\) that for all \(v\in [s_1\wedge s_2, T]\)
that is easily extended for \(v\in [0,T]\).
1.2 Proof of Lemma B.2
Let us first introduce the decomposition
thanks to the triangular inequality. Furthermore, by considering that \(s_1\le s_2\), we have
Let us set for all \(s_1\le v\le t\)
then
that implies that,
and thanks to the Gronwall lemma,
We show by a similar argument that
Hence, one can write for all \(x_1\ge 0, s_1\le s_2 \le t\)
and thanks to Lemma B.1, we have for all \(s_1,s_2\le t\)
In addition, for \(x_1,x_2\ge 0\) and \(s_2\le t\),
Let us now set for all \(s_2\le v\le t\)
then in the same way as for \(V^{x_1}_{s_1,s_2}(v)\) we obtain
Hence for all \(x_1,x_2\ge 0, s_2\le t\) we have
and thanks to Lemma B.1,
Finally, the inequality (B.5) follows from (B.12), (B.13) and (B.14).
Difference operator and Besov spaces \(B^s_{1,\infty }(\mathbb {R})\)
For \(\alpha \in \, ]0,1[\), we denote by \(\mathscr {C}^{\alpha }_b(\mathbb {R})\) the Hölder-Zygmund space that is the set of real valued functions on \(\mathbb {R}\) such that
We also introduce the difference operator defined for all \(f:\mathbb {R}\rightarrow \mathbb {R}\), for given \(m\ge 1\) integer and \(h\in \mathbb {R}\) as
that is, by a recurrence argument,
with the following properties.
Lemma C.1
Let \(m\ge 1\) be an integer, \(\alpha \in \, ]0,1[\), \(h\in \mathbb {R}\) and \(f: \mathbb {R}\rightarrow \mathbb {R}\), then:
-
1.
If \(f\in \mathscr {C}^{\alpha }_b(\mathbb {R})\),
$$\begin{aligned} \Vert \Delta ^m_hf\Vert _{\infty } \le C_m |h|^{\alpha }\Vert f\Vert _{\mathscr {C}^{\alpha }_b(\mathbb {R})}. \end{aligned}$$(C.2) -
2.
If \(f\in \mathscr {C}^{\alpha }_b(\mathbb {R})\),
$$\begin{aligned} |\Delta ^m_hf(x) - \Delta ^m_hf(y)| \le C_m \Vert f\Vert _{\mathscr {C}^{\alpha }_b(\mathbb {R})} |x-y|^{\alpha } , \forall x,y\in \mathbb {R}. \end{aligned}$$(C.3) -
3.
If \(f\in \mathcal {C}^m(\mathbb {R},\mathbb {R})\),
$$\begin{aligned} \Vert \Delta ^m_hf\Vert _{\mathbb {L}^1} \le C_m |h|^m\Vert \partial _x^mf\Vert _{\mathbb {L}^1} \end{aligned}$$(C.4) -
4.
If \(f\in \mathcal {C}_b(\mathbb {R},\mathbb {R})\) and \(g: \mathbb {R}\rightarrow \mathbb {R}\) is an integrable function, then for all \(a\in \mathbb {R}\),
$$\begin{aligned} \int _{\mathbb {R}}\Delta ^m_hf(x + a) g(x)dx = \int _{\mathbb {R}}f(x + a)\Delta ^m_{-h}g(x)dx \end{aligned}$$(C.5)
Proof
For \(0<\alpha <1\), \(h\in \mathbb {R}\) and an integer \(m\ge 1\),
-
1.
If \(f\in \mathscr {C}^{\alpha }_b(\mathbb {R})\), we proceed by a recurrence argument on \(m\ge 1\). Indeed,
$$\begin{aligned} \Vert \Delta ^1_hf\Vert _{\infty } = \sup _{x\in \mathbb {R}}\left| f(x+h)-f(x) \right| \le |h|^{\alpha }\Vert f\Vert _{\mathscr {C}^{\alpha }_b(\mathbb {R})} \end{aligned}$$and if we assume that the property holds at range m, we get
$$\begin{aligned} \Vert \Delta ^{m+1}_hf\Vert _{\infty } = \sup _{x\in \mathbb {R}}\left| \Delta ^m_hf(x+h)-\Delta ^m_hf(x) \right| \le 2C_m|h|^{\alpha }\Vert f\Vert _{\mathscr {C}^{\alpha }_b(\mathbb {R})} \end{aligned}$$that corresponds to the property with \(C_{m+1} = 2C_m\). We deduce that it holds for each integer \(m\ge 1\).
-
2.
If \(f\in \mathscr {C}^{\alpha }_b(\mathbb {R})\), then thanks to (C.1) we get for all \(x,y\in \mathbb {R}\),
$$\begin{aligned} |\Delta ^m_hf(x) - \Delta ^m_hf(y)|\le & {} \sum _{j=0}^m \begin{pmatrix} m \\ j \end{pmatrix}|f(x+jh) - f(y+jh)| \\\le & {} 2^m\Vert f\Vert _{\mathscr {C}^{\alpha }_b(\mathbb {R})}|x-y|^{\alpha }. \end{aligned}$$ -
3.
We proceed by a recurrence argument to show firstly that
$$\begin{aligned} \mathbf {(P): }~~\forall m\ge 1, ~\Delta ^m_hf(x) = h^m\int _{0}^mH_m(t)\partial _x^mf(x+th)dt, \text { if }f\in \mathcal {C}^m(\mathbb {R},\mathbb {R}) \end{aligned}$$where \(H_m\) is bounded and does not depend on f. Indeed for \(f\in \mathcal {C}^1(\mathbb {R},\mathbb {R})\), then
$$\begin{aligned} \Delta ^{1}_hf(x) = f(x+h)-f(x) = \int _x^{x+h}f'(t)dt = h\int _0^1f'(x+th)dt \end{aligned}$$and hence \(H_1(t) = 1\). We now assume that the property holds at range m, then if \(f\in \mathcal {C}^{m+1}(\mathbb {R},\mathbb {R})\),
$$\begin{aligned} \displaystyle \Delta ^{m+1}_hf(x)= & {} \Delta _h^{m}f(x+h) - \Delta ^m_hf(x)\\ \displaystyle= & {} h^m\int _0^mH_m(t)\left[ \partial _x^mf(x+(1+t)h) - \partial _x^mf(x+th) \right] dt\\ \displaystyle= & {} h^m\int _0^mH_m(t)\int _{x+th}^{x+(1+t)h} \partial _x^{m+1}f(u)dudt\\ \displaystyle= & {} h^{m+1}\int _0^mH_m(t)\int _{t}^{1+t} \partial _x^{m+1}f(x+uh)dudt\\ \displaystyle= & {} h^{m+1}\int _0^mH_m(t)\int _{0}^{m+1} \partial _x^{m+1}f(x+uh)1_{\{u\in [t,1+t]\}}dudt\\ \displaystyle= & {} h^{m+1}\int _0^{m+1}\underbrace{\left( \int _0^mH_m(t)1_{\{u\in [t,1+t]\}}dt\right) }_{=\, H_{m+1}(u)} \partial _x^{m+1}f(x+uh)du \end{aligned}$$where the function \(H_{m+1}\) is bounded and does not depend on f. We deduce that the property \(\mathbf {(P)}\) holds, and then for all \(f\in \mathcal {C}^m(\mathbb {R},\mathbb {R})\),
$$\begin{aligned} \Vert \Delta ^m_hf\Vert _{\mathbb {L}^1} \le |h|^m\Vert H_m\Vert _{\infty }\int _{\mathbb {R}}\int _0^m|\partial _x^mf(x+th)|dtdx = m\Vert H_m\Vert _{\infty }|h|^m\Vert \partial _x^mf\Vert _{\mathbb {L}^1}. \end{aligned}$$ -
4.
Let \(f\in \mathcal {C}_b(\mathbb {R},\mathbb {R})\), and \(g:\mathbb {R}\rightarrow \mathbb {R}\) an integrable function and \(a\in \mathbb {R}\), then by changing variables in the integral
$$\begin{aligned} \int _{\mathbb {R}}\Delta ^m_hf(x + a) g(x)dx= & {} \sum _{j=0}^m (-1)^{m-j} \begin{pmatrix} m \\ j \end{pmatrix}\int _{\mathbb {R}}f(x+a+jh)g(x)dx \\= & {} \sum _{j=0}^m (-1)^{m-j} \begin{pmatrix} m \\ j \end{pmatrix}\int _{\mathbb {R}}f(x+a)g(x - jh)dx \\= & {} \int _{\mathbb {R}}f(x+a)\Delta ^m_{-h}g(x)dx \end{aligned}$$
\(\square \)
The difference operator allows to define the Besov spaces \(\mathcal {B}^s_{1,\infty }(\mathbb {R})\) as the set of function \(f:\mathbb {R}\rightarrow \mathbb {R}\) that satisfy for \(m>s\)
(see Triebel [33] Theorem 2.5.12, or [34] Theorem 2.6.1). The following result follows
Lemma C.2
Let \((E,\mathcal {F}_E,\mu )\) be a measured space where \(\mu \) is a positive measure. Let \((x,y)\mapsto f(x,y)\) be an integrable function on \(\mathbb {R}\times E\), and
Then we have for all \(s>0\) and \(m>s\) an integer,
Proof
It follows from Fubini’s theorem that
Similarly, for any \(s>0\), \(m>s\) an integer and \(h\in [-1,1]\),
The result directly follows by taking the supremum in h and summing the above inequalities. \(\square \)
Furthermore, the above definition of Besov spaces allows a sufficient condition for the existence of a density for random variables given by the following result
Lemma C.3
([30], Lemma A.1) Let X be a real valued random variable. If there are an integer \(m\ge 1\), a real number \(\theta >0\), a real \(\alpha >0\), with \(\alpha<\theta <m\), and a constant \(K>0\) such that for every \(\phi \in \mathscr {C}^{\alpha }_b(\mathbb {R})\) and \(h\in \mathbb {R}\) with \(|h|\le 1\),
then X has a density \(f_X\) with respect to the Lebesgue measure on \(\mathbb {R}\). Moreover \(f_X\in \mathcal {B}^{\theta -\alpha }_{1,\infty }(\mathbb {R})\) and
The consequence that follows gives a sufficient condition that will is used in the paper
Lemma C.4
Let X be a real valued random variable. For a given non negative real valued function \(\sigma (x)\) such that \(\sigma (X)\in \mathbb {L}^1\) with \(\mathbb {E}\big [ \sigma (X) \big ] \ne 0\), if there exists an integer \(m\ge 1\), a real number \(\theta >0\), a real \(\alpha >0\) with \(\alpha<\theta <m\), and a constant \(K>0\) such that for every \(\phi \in \mathscr {C}^{\alpha }_b(\mathbb {R})\) and \(h\in \mathbb {R}\) with \(|h|\le 1\),
then X admits a density \(f_X\) on \(\{ x: \sigma (x)\ne 0 \}\). In addition, if we denote by \(f_X\) its density on this set, then the function \(x\mapsto \sigma (x)f_X(x)\) is in the Besov space \(\mathcal {B}^{\theta -\alpha }_{1,\infty }(\mathbb {R})\) and satisfies the bound
Proof
Let us denote by \(\mu (dx)\) the law of X and set
that is a non negative finite measure by assumption. Then we have
where Y is a real valued random variable of law \(\nu (dy)/\Vert \nu \Vert _{\mathbb {L}^1}\). This implies with (C.7) that
and it follows from Lemma C.3 that the probability measure \(\nu (dy)/\Vert \nu \Vert _{\mathbb {L}^1}\) admits a density \(g_X\in \mathcal {B}^{\theta -\alpha }_{1,\infty }(\mathbb {R})\) with respect to the Lebesgue measure on \(\mathbb {R}\) that satisfies the bound
Then \(\mu (dx)\) admits a density \(f_X\) on the set \(\{x: \sigma (x)\ne 0\}\) that is defined by
and the bound (C.8) follows. \(\square \)
\(\mathbb {L}^1\) estimates of derivatives of the Gaussian density
We are interested in the density of a centered Gaussian random variable with variance \(\sigma ^2\), that is
We have the following result
Lemma D.1
For each integer \(m\ge 1\),
where the constant \(C_m\) does not depend on \(\sigma \).
Proof
Because of its exponential factor, let us set
Then \(P^{\sigma }_{1}(x) = - x/\sigma \) and
that implies that each \(P^{\sigma }_m\) is a polynomial, thus
with
We easily verify by a recurrence argument on \(m\ge 1\) that
and then
\(\square \)
Rights and permissions
Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
About this article
Cite this article
Tchouanti, J. Well posedness and stochastic derivation of a diffusion-growth-fragmentation equation in a chemostat. Stoch PDE: Anal Comp 12, 466–524 (2024). https://doi.org/10.1007/s40072-023-00288-8
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s40072-023-00288-8
Keywords
- Diffusion-growth-fragmentation equation coupled with resource
- Stochastic Feller-type diffusion
- Mild formulation
- Large population approximation
- Existence of density