Similarity Solutions for the Flux limited Keller Segel System with Time Varying Chemical Decay Rate

Resumo Técnico: Soluções de Similaridade para o Sistema Keller–Segel Limitado por Fluxo com Taxa de Decaimento Químico Variável no Tempo

Enunciado do Problema
O modelo clássico de Keller–Segel (KS), embora fundamental para a quimiotaxia, sofre de duas deficiências críticas: o potencial de explosão em tempo finito das densidades celulares sob condições de massa supercrítica e a suposição de velocidade quimiotática ilimitada ($u = \chi \nabla v$), o que viola restrições biomecânicas. Os modelos de Keller–Segel limitados por fluxo (FLKS) abordam essas questões impondo uma velocidade quimiotática máxima, garantindo propagação finita e dinâmicas fisicamente significativas. No entanto, as análises existentes de FLKS geralmente assumem taxas de decaimento químico constantes. Na realidade biológica, a degradação do sinal é frequentemente mediada por enzimas (por exemplo, fosfodiesterases em Dictyostelium discoideum) cujos níveis de expressão variam ao longo dos ciclos de desenvolvimento, levando a taxas de decaimento dependentes do tempo $\kappa(t)$. Este artigo investiga um sistema FLKS unidimensional onde a taxa de decaimento químico é uma função arbitrária do tempo, visando determinar quais padrões temporais de decaimento admitem reduções analíticas de similaridade.

Metodologia
Os autores empregam uma análise de simetria de Lie sistemática combinada com a teoria de transformações de equivalência.

1. Formulação do Sistema Estendido: A função de decaimento arbitrária $\kappa(t)$ é tratada como uma variável dependente adicional sujeita a restrições ($\kappa_x = \kappa_u = \kappa_v = 0$), estendendo o espaço de fases para $(x, t, u, v, \kappa)$.
2. Classificação de Grupos: Os autores derivam as equações determinantes para os geradores infinitesimais de transformações pontuais que mapeiam a classe de equações FLKS em si mesma. Isso envolve resolver um sistema de equações diferenciais parciais para os coeficientes do gerador e a função classificadora $\kappa(t)$.
3. Construção do Sistema Ótimo: Para cada caso identificado de $\kappa(t)$, os autores constroem um sistema ótimo de subálgebras unidimensionais utilizando o método da representação adjunta. Isso garante uma classificação não redundante de todas as possíveis reduções de similaridade.
4. Reduções de Simetria: As simetrias admitidas são utilizadas para reduzir o sistema original de equações diferenciais parciais (EDPs) para equações diferenciais ordinárias (EDOs) ou equações integro-diferenciais.
5. Soluções Explícitas: Funções específicas de limitação de fluxo (formas tangente hiperbólica e Weber-Fechner) são utilizadas para derivar soluções exatas ou semi-analíticas (quadraturas) para os sistemas reduzidos.

Principais Contribuições e Resultados

1. Classificação Completa de Grupos:
   O estudo identifica que, para uma função de decaimento arbitrária $\kappa(t)$, o sistema admite apenas a simetria de núcleo de translação espacial ($X_1 = \partial_x$). No entanto, três formas funcionais distintas de $\kappa(t)$ estendem a álgebra de simetria:

   • Decaimento Constante ($\kappa(t) = \kappa_0$): Admite translação temporal ($X_2 = \partial_t$) e translação espacial.
   • Decaimento Inverso-Tempo (Lei de Potência) ($\kappa(t) = \mu/t$): Admite uma simetria de escala ($X_3 = t\partial_t + \frac{1}{2}x\partial_x + \dots$), permitindo reduções auto-similares.
   • Decaimento Exponencial ($\kappa(t) = \kappa_0 e^{\lambda t}$): Admite uma simetria combinada de translação temporal e reescalonamento químico ($X_4 = \partial_t - \lambda v \partial_v$).

2. Derivação de Reduções de Similaridade:

   • Homogeneidade Espacial: Para todos os casos, soluções espacialmente homogêneas ($u=U(t), v=V(t)$) são derivadas, reduzindo o sistema a EDOs acopladas solucionáveis via fatores integrantes.
   • Ondas Viajantes: Para decaimento constante e exponencial, soluções de ondas viajantes ($u=U(y), v=V(y)$ com $y=t-\alpha x$) são construídas. Para o caso de decaimento constante com um limitador de fluxo tangente hiperbólico, soluções explícitas em forma fechada para a densidade celular e concentração de quimioatraente são obtidas.
   • Soluções Auto-Similares: Para o caso de decaimento em lei de potência, o sistema reduz-se a um sistema acoplado de EDOs envolvendo um limitador tangente hiperbólico ao quadrado limitado. Os autores derivam uma representação exata de quadratura para a densidade celular, embora o sistema resultante para o gradiente envolva termos não locais (integro-diferenciais).

3. Soluções Analíticas Explícitas:
   O artigo fornece fórmulas explícitas para:

   • Decaimento Constante: Perfis de ondas viajantes onde a densidade celular satura devido à limitação de fluxo, equilibrando produção e degradação constante.
   • Decaimento Exponencial: Soluções envolvendo a função integral exponencial ($Ei$), demonstrando como a regulação positiva ($\lambda > 0$) pode extinguir campos de quimioatraentes, enquanto a regulação negativa ($\lambda < 0$) leva a sinais de longa duração.
   • Decaimento em Lei de Potência: Uma representação formal de perfis auto-similares, destacando a transição entre ondas viajantes transitórias e comportamento difusivo de longo prazo.

Significado e Alegações
Os autores afirmam que seu trabalho fornece uma fundação matemática rigorosa para entender quais padrões temporais de decaimento em sinalização quimiotática admitem reduções de similaridade. O estudo demonstra que:

• O decaimento dependente do tempo não é meramente uma regularização matemática, mas codifica estratégias de controle biologicamente significativas (por exemplo, regulação enzimática).
• Leis de decaimento específicas (constante, lei de potência, exponencial) correspondem a regimes dinâmicos distintos: frentes viajantes limitadas, redistribuição auto-similar e supressão ou amplificação rápida de sinais, respectivamente.
• A abordagem baseada em simetria permite progresso analítico em modelos limitados por fluxo com mecanismos de degradação realistas e variáveis no tempo, oferecendo soluções de referência para validar esquemas numéricos em dimensões superiores.

Limitações e Direções Futuras
O artigo nota modestamente que a redução auto-similar para o caso de decaimento em lei de potência leva a um sistema integro-diferencial que não admite uma redução puramente local de EDO sob funções de limitação de fluxo genéricas. Consequentemente, as fórmulas de solução explícita para este caso são interpretadas como representações formais que exigem suposições adicionais de fechamento ou integração numérica. Os autores identificam uma análise completa deste caso específico como uma direção aberta para investigação futura.