Cell-cell adhesion and multiphase Hele-Shaw problem as the singular limit of a Keller-Segel system

Resumen Técnico: La adhesión célula-célula y el problema de Hele–Shaw multifásico como el límite singular de un sistema de Keller–Segel

Planteamiento del Problema
Este artículo investiga el límite singular $\varepsilon \to 0$ de un sistema de ecuaciones acopladas de Patlak–Keller–Seger (PKS) que modela la evolución de $N$ especies de células distintas que interactúan entre sí. La motivación principal es la Hipótesis de la Adhesión Diferencial (DAH), propuesta por Malcolm Steinberg, la cual postula que las poblaciones celulares se autoorganizan minimizando la energía de adhesión interfacial, de forma análoga a los fluidos que minimizan la tensión superficial. Los autores pretenden derivar rigurosamente esta hipótesis a partir de un modelo continuo no local donde las especies interactúan mediante fuerzas de atracción no locales (quimiotaxis) y son equilibradas por una presión no lineal que evita la concentración.

El sistema describe la densidad $\rho_i^\varepsilon$ de la $i$-ésima especie evolucionando según:
$$\varepsilon \partial_t \rho_i^\varepsilon - \text{div} \left( \rho_i^\varepsilon \nabla f_m'(\sum \rho_j^\varepsilon) \right) + \text{div} \left( \rho_i^\varepsilon \sum_{j=1}^N \alpha_{ij} \nabla \phi_j^\varepsilon \right) = 0$$
donde $\phi_j^\varepsilon$ satisface una ecuación elíptica $\sigma \phi_j^\varepsilon - \varepsilon^2 \Delta \phi_j^\varepsilon = \rho_j^\varepsilon$. La matriz de interacción $A = (\alpha_{ij})$ codifica la fuerza de las interacciones, con $\alpha_{ii} > 0$ representando la autoatracción. El parámetro $\varepsilon$ representa un límite de escala correspondiente a interacciones de corto alcance y poblaciones grandes.

Metodología
El análisis procede en dos etapas principales:

1. $\Gamma$-Convergencia del Funcional de Energía:
   El sistema PKS se visualiza como un flujo de gradiente de un funcional de energía $E_\varepsilon$. Los autores establecen primero condiciones suficientes sobre la matriz de interacción $A$ (simetría, semidefinición positiva y no degeneración de las columnas) para asegurar que la energía esté acotada inferiormente. Realizan un análisis asintótico de $E_\varepsilon$ cuando $\varepsilon \to 0$.

   • Descomposición: La energía se reescribe como una suma de un término de energía potencial con una estructura de múltiples pozos y un término de energía de transición.
   • Conexión con Allen-Cahn: Mediante la introducción de una transformación que involucra la raíz cuadrada de la matriz de interacción $A = Q^T Q$, el problema se mapea a un funcional de tipo Allen-Cahn vectorial con un potencial de múltiples pozos $h(\psi)$.
   • Construcción de la Secuencia de Recuperación: Un desafío técnico crítico abordado es la construcción de una secuencia de recuperación que satisfaga las restricciones de masa para todas las especies simultáneamente. Los autores desarrollan un mecanismo de corrección específico utilizando "regiones de vacío" para ajustar las masas sin aumentar significativamente el costo energético.

2. Convergencia del Flujo de Gradiente (Dinámica):
   Basándose en el resultado estático de $\Gamma$-convergencia, los autores analizan las soluciones dependientes del tiempo. Asumiendo una condición de convergencia de energía (que la energía integrada en el tiempo de la secuencia converge a la energía límite integrada en el tiempo), demuestran la convergencia de las soluciones a un problema de frontera libre multifásico. Esto implica establecer estimaciones uniformes sobre los flujos y las presiones, y utilizar argumentos de compacidad (lemmas de tipo Lions-Aubin) para pasar al límite en las ecuaciones de continuidad.

Contribuciones Clave y Resultados

1. Derivación Rigurosa de la Energía Interfacial:
   El artículo demuestra que $E_\varepsilon$ realiza una $\Gamma$-convergencia hacia una energía interfacial límite $E_0$ de la forma:
   $$E_0(\vec{\rho}) = \sum_{i,j=0}^N \sigma_{ij} \mathcal{H}^{d-1}(\Sigma_{ij})$$
   donde $\Sigma_{ij}$ son las interfaces entre las especies $i$ y $j$ (y el vacío). Crucialmente, los coeficientes de tensión superficial $\sigma_{ij}$ no son arbitrarios, sino que están determinados por un problema geodésico en el plano de fase:
   $$\sigma_{ij} = d_A(\theta_i \vec{e}_i, \theta_j \vec{e}_j)$$
   donde $d_A$ es una distancia definida por la longitud de la trayectoria más corta en el espacio de fase ponderada por el potencial $h$. Esto proporciona un vínculo riguroso entre los coeficientes de interacción microscópicos $\alpha_{ij}$ y las tensiones superficiales macroscópicas.

2. Clasificación de Regímenes de Adhesión (Validación de la DAH):
   Los autores analizan las condiciones bajo las cuales la desigualdad triangular para las tensiones superficiales se convierte en una igualdad o es estricta. Para el caso de dos especies ($N=2$), demuestran:

   • Ordenamiento Completo (Complete Sorting): Ocurre cuando las interacciones cruzadas son cero ($\alpha_{12}=0$), lo que conduce a $\sigma_{12} = \sigma_{01} + \sigma_{02}$.
   • Engullimiento Parcial (Partial Engulfment): Ocurre para interacciones cruzadas positivas pequeñas, caracterizado por desigualdades triangulares estrictas ($\sigma_{ij} < \sigma_{ik} + \sigma_{kj}$), lo que resulta en ángulos de contacto no tangenciales en las uniones triples.
   • Engullimiento Total (Full Engulfment): Ocurre cuando las interacciones cruzadas son suficientemente fuertes en relación con las autointeracciones, lo que lleva a $\sigma_{01} = \sigma_{02} + \sigma_{21}$, donde una especie rodea completamente a la otra.
     Estos resultados validan matemáticamente los regímenes cualitativos predichos por la DAH de Steinberg.

3. Convergencia al Flujo de Hele–Shaw Multifásico:
   El segundo resultado principal establece que la evolución temporal de los agregados celulares segregados converge a un problema de frontera libre de Hele–Shaw multifásico con tensión superficial. El sistema límite se caracteriza por:

   • Incompresibilidad dentro de cada fase ($\text{div}(v_i) = 0$).
   • Una condición de salto de presión a través de las interfaces: $[p]_{\Sigma_{ij}} = 2\sigma_{ij} \kappa_{ij}$, donde $\kappa_{ij}$ es la curvatura media.
   • Condiciones cinemáticas que relacionan la velocidad de la interfaz con el gradiente de presión.
   • Condiciones de ángulo de contacto de Herring en las uniones triples y contacto ortogonal con el límite del dominio.

Significado y Reivindicaciones
El artículo afirma proporcionar una base matemática rigurosa para la Hipótesis de la Adhesión Diferencial, cerrando la brecha entre los modelos de quimiotaxis no locales y la minimización de la energía interfacial local. A diferencia de trabajos previos que dependían de aproximaciones de Cahn-Hilliard para derivar modelos de adhesión locales, este trabajo retiene la estructura genuinamente no local del sistema PKS y deriva la energía interfacial directamente en el límite de interfaz aguda.

Los autores enfatizan que su derivación de los coeficientes de tensión superficial mediante un problema geodésico en el plano de fase es una contribución novedosa, ofreciendo una fórmula precisa que relaciona la tensión superficial macroscópica con las fuerzas de interacción microscópicas. Además, identifican rigurosamente los regímenes de parámetros específicos (basados en la matriz de interacción $A$) que conducen a distintos fenómenos biológicos como el ordenamiento celular y el engullimiento, yendo más allá de las descripciones fenomenológicas hacia las pruebas analíticas. La convergencia al problema de Hele–Shaw multifásico con tensión superficial y condiciones de ángulo de contacto proporciona una descripción continua de la dinámica impulsada por la adhesión que es consistente con los modelos de mecánica de fluidos establecidos.