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

Resumo Técnico: Adesão célula-célula e o problema de Hele–Shaw multifásico como o limite singular de um sistema de Keller–Segel

Enunciado do Problema
Este artigo investiga o limite singular $\varepsilon \to 0$ de um sistema de equações acopladas de Patlak–Keller–Seger (PKS) que modelam a evolução de $N$ espécies distintas de células interagentes. A motivação principal é a Hipótese da Adesão Diferencial (DAH), proposta por Malcolm Steinberg, que postula que populações celulares se auto-organizam minimizando a energia de adesão interfacial, analogamente a fluidos que minimizam a tensão superficial. Os autores visam derivar rigorosamente esta hipótese a partir de um modelo contínuo não local onde as espécies interagem via forças de atração não locais (quimiotaxia) e são balanceadas por uma pressão não linear que impede a concentração.

O sistema descreve a densidade $\rho_i^\varepsilon$ da $i$-ésima espécie evoluindo de acordo com:
$$\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$$
onde $\phi_j^\varepsilon$ satisfaz uma equação elíptica $\sigma \phi_j^\varepsilon - \varepsilon^2 \Delta \phi_j^\varepsilon = \rho_j^\varepsilon$. A matriz de interação $A = (\alpha_{ij})$ codifica a força das interações, com $\alpha_{ii} > 0$ representando a autoatração. O parâmetro $\varepsilon$ representa um limite de escala correspondente a interações de curto alcance e grandes populações.

Metodologia
A análise procede em duas etapas principais:

1. $\Gamma$-Convergência do Funcional de Energia:
   O sistema PKS é visto como um fluxo de gradiente de um funcional de energia $E_\varepsilon$. Os autores primeiro estabelecem condições suficientes sobre a matriz de interação $A$ (simetria, semidefinitividade positiva e não degenerescência das colunas) para garantir que a energia seja limitada inferiormente. Eles realizam uma análise assintótica de $E_\varepsilon$ quando $\varepsilon \to 0$.

   • Decomposição: A energia é reescrita como uma soma de um termo de energia potencial com uma estrutura de múltiplos poços e um termo de energia de transição.
   • Conexão com Allen–Cahn: Ao introduzir uma transformação envolvendo a raiz quadrada da matriz de interação $A = Q^T Q$, o problema é mapeado para um funcional do tipo Allen–Cahn vetorial com um potencial de múltiplos poços $h(\psi)$.
   • Construção da Sequência de Recuperação: Um desafio técnico crítico abordado é a construção de uma sequência de recuperação que satisfaça as restrições de massa para todas as espécies simultaneamente. Os autores desenvolvem um mecanismo de correção específico usando "regiões de vácuo" para ajustar as massas sem aumentar significativamente o custo energético.

2. Convergência do Fluxo de Gradiente (Dinâmica):
   Baseando-se no resultado estático de $\Gamma$-convergência, os autores analisam as soluções dependentes do tempo. Assumindo uma condição de convergência de energia (que a energia integrada no tempo da sequência converge para a energia integrada no tempo do limite), eles provam a convergência das soluções para um problema de contorno livre multifásico. Isso envolve estabelecer estimativas uniformes sobre os fluxos e pressões, e utilizar argumentos de compacidade (lémas do tipo Lions–Aubin) para passar para o limite nas equações de continuidade.

Principais Contribuições e Resultados

1. Derivação Rigorosa da Energia Interfacial:
   O artigo prova que $E_\varepsilon$ $\Gamma$-converge para uma energia interfacial limite $E_0$ da forma:
   $$E_0(\vec{\rho}) = \sum_{i,j=0}^N \sigma_{ij} \mathcal{H}^{d-1}(\Sigma_{ij})$$
   onde $\Sigma_{ij}$ são as interfaces entre as espécies $i$ e $j$ (e o vácuo). Crucialmente, os coeficientes de tensão superficial $\sigma_{ij}$ não são arbitrários, mas são determinados por um problema geodésico no plano de fase:
   $$\sigma_{ij} = d_A(\theta_i \vec{e}_i, \theta_j \vec{e}_j)$$
   onde $d_A$ é uma distância definida pelo comprimento do caminho mais curto no espaço de fase ponderado pelo potencial $h$. Isto fornece um elo rigoroso entre os coeficientes de interação microscópicos $\alpha_{ij}$ e as tensões superficiais macroscópicas.

2. Classificação de Regimes de Adesão (Validação da DAH):
   Os autores analisam as condições sob as quais a desigualdade triangular para tensões superficiais torna-se uma igualdade ou é estrita. Para o caso de duas espécies ($N=2$), eles provam:

   • Ordenação Completa (Complete Sorting): Ocorre quando as interações cruzadas são nulas ($\alpha_{12}=0$), levando a $\sigma_{12} = \sigma_{01} + \sigma_{02}$.
   • Engolfamento Parcial (Partial Engulfment): Ocorre para pequenas interações cruzadas positivas, caracterizado por desigualdades triangulares estritas ($\sigma_{ij} < \sigma_{ik} + \sigma_{kj}$), resultando em ângulos de contato não tangenciais em junções triplas.
   • Engolfamento Total (Full Engulfment): Ocorre quando as interações cruzadas são suficientemente fortes em relação às autointerações, levando a $\sigma_{01} = \sigma_{02} + \sigma_{21}$, onde uma espécie envolve completamente a outra.
     Estes resultados validam matematicamente os regimes qualitativos previstos pela DAH de Steinberg.

3. Convergência para o Fluxo de Hele–Shaw Multifásico:
   O segundo resultado principal estabelece que a evolução temporal dos agregados celulares segregados converge para um problema de contorno livre de Hele–Shaw multifásico com tensão superficial. O sistema limite é caracterizado por:

   • Incompressibilidade dentro de cada fase ($\text{div}(v_i) = 0$).
   • Uma condição de salto de pressão através das interfaces: $[p]_{\Sigma_{ij}} = 2\sigma_{ij} \kappa_{ij}$, onde $\kappa_{ij}$ é a curvatura média.
   • Condições cinemáticas relacionando a velocidade da interface ao gradiente de pressão.
   • Condições de ângulo de contato de Herring em junções triplas e contato ortogonal com o limite do domínio.

Significância e Alegações
O artigo alega fornecer uma base matemática rigorosa para a Hipótese da Adesão Diferencial, preenchendo a lacuna entre modelos de quimiotaxia não locais e a minimização de energia interfacial local. Diferente de trabalhos anteriores que dependiam de aproximações de Cahn–Hilliard para derivar modelos de adesão locais, este trabalho retém a estrutura genuinamente não local do sistema PKS e deriva a energia interfacial diretamente no limite de interface nítida.

Os autores enfatizam que sua derivação dos coeficientes de tensão superficial via um problema geodésico no plano de fase é uma contribuição inovadora, oferecendo uma fórmula precisa que relaciona as forças de interação microscópicas à tensão superficial macroscópica. Além disso, eles identificam rigorosamente os regimes de parâmetros específicos (baseados na matriz de interação $A$) que levam a distintos fenômenos biológicos, como ordenação celular e engolfamento, indo além de descrições fenomenológicas para provas analíticas. A convergência para o problema de Hele–Shaw multifásico com tensão superficial e condições de ângulo de contato fornece uma descrição contínua da dinâmica impulsionada pela adesão que é consistente com modelos estabelecidos de mecânica de fluidos.