Uniform boundedness for the two-dimensional Keller-Segel system with Gompertz growth

Titolo: Risolubilità globale e limitatezza uniforme per un modello chemiotattico di Keller-Segel con crescita di Gompertz

1. Il Problema

Il documento analizza il sistema di Keller-Segel in due dimensioni spaziali, un modello matematico fondamentale per descrivere il movimento direzionale di cellule (chemiotassi) in risposta a gradienti chimici. Il modello classico è noto per la sua tendenza a sviluppare fenomeni di aggregazione cellulare che possono portare a "blow-up" (esplosione della densità cellulare in tempo finito), specialmente in dimensioni $n \ge 2$.

Per prevenire questo comportamento non fisico e ottenere dinamiche più realistiche, il termine di trasporto è spesso modificato includendo una fonte di crescita logistica ($au - bu^2$). Tuttavia, l'obiettivo di questo studio è investigare l'effetto di una fonte di crescita più debole, la funzione di Gompertz:
$$f(u) = \alpha u \ln\left(\frac{K}{u}\right)$$
dove $\alpha > 0$ è il tasso di crescita intrinseco e $K > 0$ è la capacità portante.
Il sistema studiato è:
$$\begin{cases} u_t = \Delta u - \chi \nabla \cdot (u \nabla v) + \alpha u \ln(K/u) \\ \tau v_t = \Delta v - v + u \end{cases}$$
con condizioni al contorno di Neumann omogenee e dati iniziali non negativi. Qui $\chi > 0$ è la sensibilità chemiotattica e $\tau \in \{0, 1\}$ distingue tra il caso parabolico-ellittico ($\tau=0$) e il caso completamente parabolico ($\tau=1$).

La sfida principale risiede nel fatto che la funzione di Gompertz ha un effetto di smorzamento (mortalità) più debole rispetto alla crescita logistica classica e anche rispetto alle fonti "sotto-logistiche" studiate in letteratura recente (es. Xiang, 2018), rendendo più difficile garantire la limitatezza uniforme delle soluzioni.

2. Metodologia

Gli autori adottano un approccio analitico rigoroso basato su stime a priori e criteri di estensione delle soluzioni. La metodologia si articola nei seguenti passaggi:

• Esistenza Locale: Si stabilisce l'esistenza di una soluzione classica locale in un intervallo di tempo $[0, T_{max})$ utilizzando procedure di punto fisso standard. Viene introdotto un criterio di estensione (Lemma 2.1) che afferma che la soluzione può essere prolungata globalmente a meno che la norma $L^\infty$ di $u$ non diverga.
• Criterio di Limitatezza: Si fa riferimento a un criterio noto (Lemma 2.2) che lega la limitatezza globale di $u$ alla limitatezza a priori di due quantità: l'entropia $\int_\Omega u \ln u$ e l'energia del segnale chimico $\int_\Omega |\nabla v|^2$.
• Stime A Priori:
  • Massa Totale: Si dimostra che la massa totale $\int_\Omega u$ è limitata superiormente da una costante $M$ che dipende dai dati iniziali e dalla capacità portante $K$, sfruttando la concavità della funzione di Gompertz e la disuguaglianza di Jensen.
  • Disuguaglianza di Gagliardo-Nirenberg: Si utilizza per legare la norma $L^2$ di $u$ al gradiente e alla massa, ottenendo una disuguaglianza differenziale chiave.
  • Funzionale di Lyapunov: Viene definito un funzionale energetico $F(t) = \int_\Omega u \ln u + \tau b \int_\Omega |\nabla v|^2$. Derivando $F(t)$ rispetto al tempo e utilizzando disuguaglianze di Young e le stime sulla massa, gli autori costruiscono una disuguaglianza differenziale del tipo $F'(t) \le \delta_1 - \delta_2 F(t)$.
• Gestione dei Casi $\tau=0$ e $\tau=1$:
  • Per $\tau=1$, il controllo diretto sul funzionale $F(t)$ garantisce la limitatezza di entrambi i termini.
  • Per $\tau=0$ (caso parabolico-ellittico), poiché il termine $\int |\nabla v|^2$ non appare nella derivata temporale di $F$, gli autori devono derivare una stima separata per $\|u\|_{L^2}$ moltiplicando l'equazione per $u$ e utilizzando nuovamente le disuguaglianze di Gagliardo-Nirenberg e la limitatezza dell'entropia.

3. Risultati Chiave

Il risultato principale è enunciato nel Teorema 1.1. Gli autori dimostrano che, sotto opportune condizioni sui parametri, il sistema ammette una soluzione globale e uniformemente limitata.

Le condizioni sufficienti sono:

1. Condizione sulla capacità portante: $K > e^{-2/\alpha}$. Questa condizione garantisce che il termine di crescita non diventi troppo negativo o instabile per densità elevate.
2. Condizione sulla sensibilità chemiotattica: $\chi M \le \frac{1}{2 C_{GN}^4}$, dove $M = \max \{ \int_\Omega u_0, K|\Omega| \}$ e $C_{GN}$ è la costante di Gagliardo-Nirenberg.

Conclusione del Teorema: Se queste condizioni sono soddisfatte, esiste una coppia unica di funzioni $(u, v)$ che risolve il problema per ogni $t \in [0, \infty)$, con $u$ e $v$ limitate uniformemente in $L^\infty(\Omega \times (0, \infty))$.

4. Contributi e Significato

• Estensione della letteratura: Questo lavoro estende i risultati noti sulla limitatezza delle soluzioni di Keller-Segel a una classe di fonti di crescita più debole (Gompertz). Mentre la crescita logistica quadratica ($u^2$) è un potente smorzatore, la funzione di Gompertz ($u \ln(1/u)$) ha un effetto di smorzamento asintoticamente più debole per grandi valori di $u$.
• Rilevanza Biologica: La scelta della funzione di Gompertz è motivata biologicamente, in quanto descrive più accuratamente la proliferazione delle cellule tumorali rispetto al modello logistico, catturando la decelerazione asimmetrica della crescita man mano che ci si avvicina alla capacità portante.
• Soglia Critica: Il lavoro identifica chiaramente come la stabilità del sistema dipenda da un bilancio tra la forza della chemiotassi ($\chi$), la massa iniziale e i parametri della crescita di Gompertz ($\alpha, K$). Dimostra che anche con uno smorzamento debole, è possibile prevenire il blow-up se la chemiotassi non è troppo forte rispetto alla massa totale.
• Metodologia Robusta: La dimostrazione fornisce un quadro tecnico solido per trattare termini di reazione non lineari complessi in sistemi di reazione-diffusione-chemiotassi, utilizzando combinazioni sofisticate di disuguaglianze funzionali e stime energetiche.

In sintesi, il paper conferma che l'introduzione di una sorgente di crescita di tipo Gompertz è sufficiente a garantire la globalità e la limitatezza delle soluzioni in due dimensioni, a patto che la sensibilità chemiotattica sia sufficientemente piccola rispetto alla massa totale del sistema, offrendo così un modello matematico più realistico per lo studio della dinamica tumorale.