Sobolev estimates for the Keller-Segel system and applications to the JKO scheme

Resumen Técnico: Estimaciones de Sobolev para el sistema de Keller-Segel y aplicaciones al esquema JKO

Planteamiento del Problema
El artículo aborda las propiedades analíticas del sistema de Keller-Segel (KS) con difusión lineal en un dominio acotado, convexo y suave $\Omega \subset \mathbb{R}^d$. El sistema está dado por:
$$\partial_t\rho - \Delta\rho + \text{div}(\rho\nabla u[\rho]) = 0, \quad -\Delta u[\rho] = \rho$$
sujeto a condiciones de contorno de no flujo y de Dirichlet. Los desafíos principales en este contexto son:

1. Explosión en tiempo finito (blow-up): Las soluciones pueden explotar si la masa inicial excede un umbral crítico (por ejemplo, $8\pi$ en 2D). El artículo se centra en un "régimen subcrítico de datos pequeños" donde la norma $L^{d/2}$ inicial es suficientemente pequeña para evitar la explosión.
2. Existencia y Regularidad: Aunque existen soluciones débiles, obtener cotas $L^\infty$ y estimaciones de Sobolev de orden superior (específicamente $L^\infty_t W^{1,p}_x$ para todo $1 \le p < \infty$) en dimensiones $d \ge 3$ no es trivial debido a la falta de control sobre la parte negativa de la energía libre y la singularidad del núcleo de interacción.
3. Convergencia del esquema JKO: El esquema de Jordan-Kinderlehrer-Otto (JKO), una discretización temporal variacional de los flujos de gradiente en el espacio de Wasserstein, se ha establecido para el sistema KS. Sin embargo, extender los resultados de convergencia fuerte (específicamente en $L^2_t H^2_x$) conocidos para la ecuación de Fokker-Planck requiere superar dificultades específicas relacionadas con el potencial dependiente del tiempo y no suave $u[\rho]$.

Metodología
El autor emplea una combinación de análisis funcional, teoría de transporte óptimo y técnicas de flujo de gradiente discreto.

1. Desigualdad Funcional (tipo Brezis-Gallouët-Wainger):
   Una herramienta central es una nueva desigualdad funcional (Teorema 1.1) inspirada en la desigualdad de Brezis-Gallouët-Wainger. Esta proporciona una cota $L^\infty$ para la Hessiana $D^2 u$ en términos de las normas $L^\infty$ y $C^{0,\alpha}$ del Laplaciano $\Delta u$:
   $$\|D^2 u\|_{L^\infty} \le C \| \Delta u \|_{L^\infty} \left( 1 + \log \left( \frac{\|\Delta u\|_{C^{0,\alpha}}}{\|\Delta u\|_{L^\infty}} \right) \right)$$
   Este control logarítmico es crucial para cerrar las estimaciones en presencia de la no linealidad.

2. Disipación del Funcional $F_p$:
   El artículo analiza la disipación del funcional $F_p[\rho] = \int_\Omega \rho |\nabla(\log \rho - u[\rho])|^p dx$. Al calcular la derivada temporal (en el caso continuo) o la diferencia discreta (en el esquema JKO), el autor deriva una desigualdad diferencial de la forma:
   $$\frac{d}{dt} F_p \lesssim F_p (\log(F_p + 1) + 1)$$
   El crecimiento logarítmico en el lado derecho, derivado de la desigualdad funcional y las estimaciones de Calderón-Zygmund, permite obtener cotas globales en el tiempo mediante la lemma de Gronwall, evitando la explosión de las normas de Sobolev.

3. El Esquema JKO Penalizado:
   Para manejar la existencia en altas dimensiones ($d \ge 3$) donde la energía libre estándar no está acotada inferiormente, el autor introduce un funcional penalizado $\tilde{J}$. Este incluye un término artificial que controla la norma $L^{d/2}$. El artículo demuestra que, para datos iniciales suficientemente pequeños, este término de penalización permanece inactivo y el esquema coincide con el esquema JKO estándar.

4. Estimaciones Discretas y Desigualdad de Cinco Gradientes:
   Para probar la convergencia, el autor adapta la "desigualdad de cinco gradientes" (originalmente para estimaciones BV en transporte óptimo) al entorno discreto de JKO. Esto permite controlar la disipación del funcional de información de Fisher $F_2$ a nivel discreto. Los pasos clave incluyen:

   • Establecer una cota $L^\infty$ uniforme para las soluciones discretas utilizando un método de iteración de Alikakos adaptado al esquema JKO.
   • Probar una cota inferior estricta (positividad) para las densidades discretas, lo cual es esencial para controlar los potenciales de Kantorovich y sus hessianas.
   • Mostrar que los potenciales de Kantorovich convergen a cero con una tasa específica, cerrando la brecha entre la disipación discreta y los límites continuos.

Contribuciones Clave y Resultados

• Teorema 1.1 (Desigualdad Funcional): Demuestra una nueva desigualdad que acota $\|D^2 u\|_{L^\infty}$ usando la interpolación logarítmica entre las normas $L^\infty$ y $C^{0,\alpha}$ de $\Delta u$. Esta es una refinación de las inclusiones de Sobolev clásicas adecuada para la estructura específica del sistema KS.
• Teorema 1.5 (Estimaciones de Sobolev): Establece que, para datos iniciales pequeños (norma $L^{d/2}$), la solución débil al sistema KS existe globalmente y está acotada en $L^\infty((0,T) \times \Omega) \cap L^\infty(0, T; W^{1,p}(\Omega))$ para todo $1 \le p < \infty$. La prueba se basa en la disipación de $F_p$ y la desigualdad diferencial logarítmica derivada.
• Teorema 1.7 (Convergencia del Esquema JKO):
  • Existencia: Demuestra la existencia de la secuencia del esquema JKO en altas dimensiones utilizando el funcional penalizado.
  • Regularidad: Muestra que la curva discreta $\rho_\tau$ está uniformemente acotada en $L^\infty((0,T) \times \Omega)$ y $L^\infty(0, T; W^{1,p}(\Omega))$.
  • Convergencia Fuerte: Demuestra que la curva discreta $\rho_\tau$ converge fuertemente en $L^2(0, T; H^2(\Omega))$ a la solución débil del sistema KS a medida que el paso de tiempo $\tau \to 0$.

Significado y Reivindicaciones
El artículo afirma extender los resultados recientes de Santambrogio y Toshpulatov (quienes demostraron la convergencia $L^2_t H^2_x$ para la ecuación de Fokker-Planck) al contexto de agregación-difusión del sistema de Keller-Segel. El autor destaca tres dificultades específicas superadas:

1. Existencia en Altas Dimensiones: Superar la falta de control de la energía libre en $d \ge 3$ mediante el truco de la penalización.
2. Cotas Inferiores: Establecer una cota inferior uniforme para la densidad discreta $\rho_\tau$ sin asumir cotas a priori sobre $D^2 u$, lo cual no es trivial ya que el potencial $u$ depende de $\rho$.
3. Regularidad de los Potenciales: Probar estimaciones $C^{0,\alpha}$ uniformes para las curvas JKO para controlar la convergencia de los potenciales de Kantorovich, un paso que anteriormente dependía de estimaciones $L^\infty_t W^{1,p}_x$ que solo se conocían para $p < 2$ en el contexto de KS.

El autor señala que, mientras que la cota $L^\infty$ para la PDE continua se consideraba un problema abierto en parte de la literatura (por ejemplo, [23]), este trabajo proporciona una prueba constructiva vía el esquema JKO. Los resultados se presentan como una justificación rigurosa de la convergencia del esquema JKO en normas de Sobolev altas para el sistema de Keller-Segel, apoyándose en la interacción entre la desigualdad funcional y la estructura de flujo de gradiente.