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

Resumen Técnico: Solvabilidad Global de un Modelo de Quimiotaxis con Fuente Gompertz

1. Planteamiento del Problema

El artículo aborda el sistema clásico de Keller–Segel en dos dimensiones espaciales, el cual modela el movimiento direccional de células (quimiotaxis) en respuesta a un gradiente químico auto-producido. El modelo base, conocido por exhibir fenómenos de agregación celular y "blow-up" (explosión en tiempo finito) en dimensiones $n \ge 2$, se modifica aquí incorporando una fuente de crecimiento poblacional de tipo Gompertz.

El sistema estudiado es:
$$\begin{cases} u_t = \Delta u - \chi \nabla \cdot (u \nabla v) + f(u), & \text{en } \Omega \times (0, T_{max}) \\ \tau v_t = \Delta v - v + u, & \text{en } \Omega \times (0, T_{max}) \\ \frac{\partial u}{\partial \nu} = \frac{\partial v}{\partial \nu} = 0, & \text{en } \partial\Omega \times (0, T_{max}) \\ u(x,0) = u_0(x), \quad \tau v(x,0) = \tau v_0(x), & \text{en } \Omega \end{cases}$$
Donde:

• $u(x,t)$ es la densidad celular y $v(x,t)$ la concentración de la señal química.
• $\chi > 0$ es la sensibilidad quimiotáctica.
• $\tau \in \{0, 1\}$ distingue entre el caso parabólico-elíptico ($\tau=0$) y el totalmente parabólico ($\tau=1$).
• La fuente de crecimiento es Gompertz: $f(u) = \alpha u \ln\left(\frac{K}{u}\right)$, con $\alpha, K > 0$.

Motivación:
A diferencia del término logístico clásico ($au - bu^2$), que proporciona un amortiguamiento fuerte y garantiza la acotación global en 2D, la función Gompertz presenta un efecto de amortiguamiento más débil para grandes valores de $u$. El objetivo es determinar si este amortiguamiento más débil es suficiente para prevenir el "blow-up" y asegurar la existencia global y la acotación uniforme de las soluciones.

2. Metodología

Los autores emplean un enfoque analítico basado en la teoría de ecuaciones en derivadas parciales (EDP) no lineales, utilizando las siguientes herramientas clave:

1. Existencia Local y Criterio de Extensión:

   • Se establece la existencia local de soluciones clásicas mediante procedimientos de punto fijo.
   • Se define un criterio de extensión: la solución puede extenderse globalmente ($T_{max} = \infty$) si y solo si la norma $L^\infty$ de $u$ permanece acotada. Si $T_{max} < \infty$, entonces $\limsup_{t \to T_{max}} \|u(\cdot, t)\|_{L^\infty} = \infty$.

2. Estimaciones A Priori (Lema 3.1):

   • Se demuestra que la masa total $\int_\Omega u$ está acotada uniformemente en el tiempo por una constante $M$ que depende de la masa inicial y la capacidad de carga $K$.
   • Se utiliza la desigualdad de Gagliardo-Nirenberg para relacionar la norma $L^2$ de $u$ con la integral de $\frac{|\nabla u|^2}{u}$ (entropía), lo cual es crucial para controlar los términos de difusión.

3. Funcional de Lyapunov y Desigualdades de Absorción:

   • Se define un funcional de energía combinado: $F(t) = \int_\Omega u \ln(u) + \tau b \int_\Omega |\nabla v|^2$.
   • Mediante diferenciación temporal, integración por partes y el uso de desigualdades de Young, se deriva una desigualdad diferencial para $F(t)$.
   • El objetivo es lograr una desigualdad de absorción de la forma $F'(t) \le \delta_1 - \delta_2 F(t)$, lo que implica que $F(t)$ permanece acotada.

4. Tratamiento de Casos ($\tau = 0$ y $\tau = 1$):

   • Caso $\tau=1$ (Parabólico-Parabólico): El acotamiento de $F(t)$ controla directamente tanto la entropía de $u$ como el gradiente de $v$.
   • Caso $\tau=0$ (Parabólico-Elíptico): Dado que $v$ no tiene dinámica temporal propia, se requiere una estimación adicional. Se multiplica la ecuación de $u$ por $u$ para obtener una estimación en $L^2$, utilizando nuevamente desigualdades de Gagliardo-Nirenberg y la acotación de la entropía para cerrar el sistema de estimaciones.

3. Contribuciones Clave y Resultados Principales

Teorema Principal (Teorema 1.1):
Bajo supuestos de regularidad adecuados en el dominio $\Omega \subset \mathbb{R}^2$ y datos iniciales, si se cumplen las siguientes condiciones:

1. Capacidad de carga suficiente: $K > e^{-2/\alpha}$.
2. Restricción en la sensibilidad quimiotáctica: $\chi M \le \frac{1}{2 C_{GN}^4}$, donde $M = \max\{\int_\Omega u_0, K|\Omega|\}$ y $C_{GN}$ es una constante de Gagliardo-Nirenberg.

Entonces, existe un par único de funciones $(u, v)$ que resuelve el sistema globalmente en el tiempo ($T_{max} = \infty$) y que es uniformemente acotado en $L^\infty(\Omega \times (0, \infty))$.

Comparación con la Literatura Existente:

• El artículo sitúa la fuente Gompertz ($f_G$) en una jerarquía de amortiguamiento: $f_G < f_S \lesssim f_L$.
  • $f_L$: Logística clásica (amortiguamiento fuerte).
  • $f_S$: Fuentes sub-logísticas estudiadas en Xiang (2018) (amortiguamiento intermedio).
  • $f_G$: Gompertz (amortiguamiento muy débil para grandes $u$).
• Innovación: Este trabajo demuestra que, a pesar de que la función Gompertz es más débil que las consideradas en trabajos previos (como Xiang, 2018) y no pertenece a la misma clase funcional, sigue siendo suficiente para garantizar la acotación global en 2D, siempre que se cumplan las condiciones de parámetros mencionadas.

4. Significado e Impacto

• Relevancia Biológica: El modelo Gompertz es ampliamente utilizado en oncología para describir el crecimiento tumoral, ya que captura mejor la desaceleración asimétrica del crecimiento a medida que el tumor se acerca a su capacidad de carga, en comparación con el modelo logístico simétrico.
• Avance Matemático: El resultado demuestra que los mecanismos de amortiguamiento "sub-críticos" (más débiles que el logístico cuadrático) pueden estabilizar sistemas de quimiotaxis en 2D. Esto amplía el entendimiento de las condiciones mínimas necesarias para prevenir la agregación celular descontrolada (blow-up).
• Robustez del Modelo: Se establece que la estructura del sistema de Keller–Segel es robusta frente a fuentes de crecimiento no logísticas, siempre que la sensibilidad quimiotáctica $\chi$ no sea excesiva en relación con la masa inicial y la capacidad de carga del entorno.

En conclusión, el artículo proporciona una justificación teórica rigurosa para el uso de fuentes Gompertz en modelos de quimiotaxis, garantizando que las soluciones biológicamente relevantes (densidades celulares finitas) existen globalmente bajo condiciones controlables.