Global existence and stability in a class of chemotaxis systems with lethal interactions, nonlinear diffusion and production

Résumé Technique : Existence Globale et Stabilité dans une Classe de Systèmes de Chimiotaxie avec Interactions Létales

Formulation du Problème
Cet article étudie une classe de systèmes de chimiotaxie modélisant l'interaction entre des populations bactériennes (spécifiquement E. coli) et une substance chimique toxique auto-produite (peroxyde d'hydrogène, $H_2O_2$). Le système prend en compte des interactions létales où le produit chimique agit à la fois comme un chimiorépulsif et un agent de mortalité. L'étude se concentre sur deux formulations mathématiques distinctes au sein d'un domaine borné et lisse $\Omega \subset \mathbb{R}^n$ ($n \ge 1$) soumis à des conditions aux limites de Neumann homogènes :

1. Système Totalement Parabolique ($\tau=1$) : La densité de population $u$ et la concentration chimique $v$ évoluent toutes deux selon des équations paraboliques.
2. Système Parabolique-Elliptique ($\tau=0$) : La densité de population $u$ suit une équation parabolique, tandis que la concentration chimique $v$ est régie par une équation elliptique (approximation de l'état quasi-stationnaire).

Le système directeur (1.1) incorpore :

• Diffusion Non Linéaire : Représentée par $d_1 \nabla \cdot (D(u)\nabla u)$, où $D(s) \ge (s+1)^\alpha$.
• Répulsion Chimiotactique Non Linéaire : Modélisée par $\chi \nabla \cdot (S(u)\nabla v)$, où $0 \le S(s) \le s(s+1)^\beta$.
• Croissance Logistique et Létalité : Un terme de croissance logistique $ru(1-u^{\kappa-1})$ et un terme de mortalité $-\mu uv$.
• Dynamique Chimique : Production proportionnelle à $u^m$, taux de décomposition $b$, et une source externe $f(x,t)$.

Méthodologie
Les auteurs emploient un cadre analytique rigoureux combinant plusieurs techniques avancées en équations aux dérivées partielles (EDP) :

• Existence Locale et Régularité : Des arguments de point fixe standards et la théorie de la régularité parabolique sont utilisés pour établir l'existence locale de solutions classiques.
• Estimations A Priori : Pour prouver l'existence globale, les auteurs dérivent des bornes $L^p$ pour la solution $u$. Cela implique de tester les équations avec des puissances spécifiques de $(u+1)$ et d'utiliser l'inégalité de Gagliardo-Nirenberg pour gérer les termes non linéaires.
• Régularité de Sobolev Maximale : Cette théorie est appliquée à l'équation chimique pour contrôler le gradient de $v$ et le laplacien $\Delta v$, ce qui est crucial pour gérer les termes de autodiffusion chimiotactique.
• Fonctionnelles de Lyapunov : Pour analyser le comportement à long terme, des fonctionnelles d'énergie soigneusement construites (fonctionnelles de Lyapunov) sont utilisées. Ces fonctionnelles combinent des termes de type entropie pour la population et des normes $L^2$ pour la concentration chimique.
• Itération de Moser-Alikakos : Cette technique est employée pour élever les bornes $L^p$ en bornes $L^\infty$, assurant le caractère uniformément borné des solutions.

Contributions Clés et Résultats

1. Existence Globale et Bornage (Théorème 1.1) :
   L'article établit l'existence de solutions classiques uniques et globalement bornées pour des dimensions spatiales arbitraires $n \ge 1$. La condition principale pour l'existence globale est une contrainte sur les exposants régissant la diffusion, la chimiotaxie et la production :
   $$\beta + m < \alpha + \frac{2}{n}$$
   Cette inégalité garantit que la diffusion non linéaire ($\alpha$) est suffisamment forte pour contrer les effets déstabilisateurs de la répulsion chimiotactique non linéaire ($\beta$) et de la production chimique ($m$). Sous cette condition, les solutions satisfont :
   $$\|u(\cdot, t)\|_{L^\infty(\Omega)} + \|v(\cdot, t)\|_{W^{1,\infty}(\Omega)} \le C$$
   pour tout $t > 0$.

2. Stabilité Asymptotique des États de Coexistence (Théorème 1.2) :
   Lorsque l'apport externe $f$ est constant ($f_c$) et que les paramètres d'interaction satisfont $\mu f_c < br$ (permettant un équilibre positif), le système converge exponentiellement vers un état stationnaire de coexistence spatialement homogène $(u^*, v^*)$. Cette stabilité est maintenue si :

   • $\kappa \ge 2$ et $m=1$.
   • La sensibilité chimiotactique est suffisamment faible : $\chi^2 < \frac{4d_1 d_2 \mu}{a u^*}$.
   • La diffusion domine la sensibilité chimiotactique dans un sens spécifique : $2\beta \le \alpha$.

3. Stabilité Asymptotique des États d'Extinction (Théorème 1.3) :
   Dans le régime où l'interaction létale est dominante ($\mu f_c \ge br$), le système converge exponentiellement vers un « état de disparition » où la population bactérienne devient éteinte ($\bar{u}=0$) et la concentration chimique se stabilise à $\bar{v} = f_c/b$. Ce résultat suppose $\kappa=2$ et $m=1$.

Signification et Revendications
Les auteurs positionnent ce travail comme la première étude traitant des interactions létales dans les systèmes de chimiotaxie qui incorporent simultanément une diffusion non linéaire et des mécanismes de production non linéaires.

• Extension de la Théorie de Keller-Segel : L'article revendique une extension de l'analyse de l'exposant critique trouvée dans les modèles classiques de Keller-Segel (où le seuil critique est typiquement $2/n$). En introduisant l'inégalité $\beta + m < \alpha + 2/n$, les auteurs démontrent comment la diffusion non linéaire ($\alpha$) peut stabiliser le système même dans des dimensions plus élevées ou avec une sensibilité chimiotactique plus forte, récupérant et étendant efficacement les résultats précédents de Liu, Tao et Winkler.
• Pertinence Biologique : Les résultats mathématiques fournissent une base théorique pour comprendre comment les populations d'E. coli peuvent persister ou s'effondrer en présence de toxines auto-générées. L'analyse suggère que la diffusion non linéaire agit comme un mécanisme de stabilisation crucial, empêchant l'explosion de population en temps fini (blow-up) et garantissant que le système se stabilise soit dans une coexistence stable, soit dans un état d'extinction contrôlée, selon l'équilibre entre les taux de croissance et les interactions létales.
• Robustesse : Les résultats sont valables tant pour les contextes totalement paraboliques que parabolique-elliptiques, offrant une vue complète de la dynamique du système sous différentes hypothèses d'échelles de temps pour le signal chimique.

L'article conclut que sous les hypothèses structurelles dérivées, l'interaction complexe entre la répulsion chimiotactique, la diffusion non linéaire et les interactions létales ne conduit pas à une croissance non bornée, mais plutôt à des comportements à long terme bien définis et stables.