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

Resumen Técnico: Soluciones de Similitud para el Sistema Keller–Segel Limitado por Flujo con Tasa de Decaimiento Químico Variable en el Tiempo

Enunciado del Problema
El modelo clásico de Keller–Segel (KS), aunque fundamental para la quimiotaxis, adolece de dos deficiencias críticas: la posibilidad de explosión en tiempo finito de las densidades celulares bajo condiciones de masa supercrítica y la asunción de una velocidad quimiotáctica ilimitada ($u = \chi \nabla v$), lo cual viola las restricciones biomecánicas. Los modelos de Keller–Segel limitados por flujo (FLKS) abordan estos problemas imponiendo una velocidad quimiotáctica máxima, asegurando una propagación finita y dinámicas físicamente significativas. Sin embargo, los análisis existentes de FLKS suelen asumir tasas de decaimiento químico constantes. En la realidad biológica, la degradación de señales a menudo está mediada por enzimas (por ejemplo, fosfodiesterasas en Dictyostelium discoideum) cuyos niveles de expresión varían a lo largo de los ciclos de desarrollo, lo que conduce a tasas de decaimiento dependientes del tiempo $\kappa(t)$. Este artículo investiga un sistema FLKS unidimensional donde la tasa de decaimiento químico es una función arbitraria del tiempo, con el objetivo de determinar qué patrones de decaimiento temporal admiten reducciones analíticas de similitud.

Metodología
Los autores emplean un análisis de simetría de Lie sistemático combinado con la teoría de transformaciones de equivalencia.

1. Formulación del Sistema Extendido: La función de decaimiento arbitraria $\kappa(t)$ se trata como una variable dependiente adicional sujeta a restricciones ($\kappa_x = \kappa_u = \kappa_v = 0$), extendiendo el espacio de fases a $(x, t, u, v, \kappa)$.
2. Clasificación de Grupos: Los autores derivan las ecuaciones determinantes para los generadores infinitesimales de transformaciones puntuales que mapean la clase de ecuaciones FLKS en sí misma. Esto implica resolver un sistema de ecuaciones diferenciales parciales para los coeficientes del generador y la función clasificadora $\kappa(t)$.
3. Construcción del Sistema Óptimo: Para cada caso identificado de $\kappa(t)$, los autores construyen un sistema óptimo de subálgebras unidimensionales utilizando el método de la representación adjunta. Esto asegura una clasificación no redundante de todas las posibles reducciones de similitud.
4. Reducciones por Simetría: Las simetrías admitidas se utilizan para reducir el sistema original de ecuaciones diferenciales parciales (EDP) a ecuaciones diferenciales ordinarias (EDO) o ecuaciones integro-diferenciales.
5. Soluciones Explícitas: Se utilizan funciones específicas de limitación de flujo (formas tangente hiperbólica y Weber-Fechner) para derivar soluciones exactas o semianalíticas (cuadraturas) para los sistemas reducidos.

Contribuciones y Resultados Clave

1. Clasificación Completa de Grupos:
   El estudio identifica que, para una función de decaimiento arbitraria $\kappa(t)$, el sistema admite únicamente la simetría de núcleo de traslación espacial ($X_1 = \partial_x$). Sin embargo, tres formas funcionales distinguidas de $\kappa(t)$ extienden el álgebra de simetría:

   • Decaimiento Constante ($\kappa(t) = \kappa_0$): Admite traslación temporal ($X_2 = \partial_t$) y traslación espacial.
   • Decaimiento Inverso-Temporal (Ley de Potencia) ($\kappa(t) = \mu/t$): Admite una simetría de escalamiento ($X_3 = t\partial_t + \frac{1}{2}x\partial_x + \dots$), permitiendo reducciones autosimilares.
   • Decaimiento Exponencial ($\kappa(t) = \kappa_0 e^{\lambda t}$): Admite una simetría combinada de traslación temporal y reescalado químico ($X_4 = \partial_t - \lambda v \partial_v$).

2. Derivación de Reducciones de Similitud:

   • Homogeneidad Espacial: Para todos los casos, se derivan soluciones homogéneas espacialmente ($u=U(t), v=V(t)$), reduciendo el sistema a EDO acopladas resolubles mediante factores integrantes.
   • Ondas Viajeras: Para decaimiento constante y exponencial, se construyen soluciones de ondas viajeras ($u=U(y), v=V(y)$ con $y=t-\alpha x$). Para el caso de decaimiento constante con un limitador de flujo tangente hiperbólico, se obtienen soluciones explícitas en forma cerrada para la densidad celular y la concentración de quimioatrayente.
   • Soluciones Autosimilares: Para el caso de decaimiento de ley de potencia, el sistema se reduce a un sistema acoplado de EDO que involucra un limitador tanh-cuadrado acotado. Los autores derivan una representación exacta por cuadratura para la densidad celular, aunque el sistema resultante para el gradiente involucra términos no locales (integro-diferenciales).

3. Soluciones Analíticas Explícitas:
   El artículo proporciona fórmulas explícitas para:

   • Decaimiento Constante: Perfiles de ondas viajeras donde la densidad celular se satura debido a la limitación de flujo, equilibrando la producción y la degradación constante.
   • Decaimiento Exponencial: Soluciones que involucran la función integral exponencial ($Ei$), demostrando cómo la sobre-regulación ($\lambda > 0$) puede extinguir los campos de quimioatrayente, mientras que la sub-regulación ($\lambda < 0$) conduce a señales de larga duración.
   • Decaimiento de Ley de Potencia: Una representación formal de perfiles autosimilares, destacando la transición entre ondas viajeras transitorias y comportamiento difusivo a largo plazo.

Significado y Afirmaciones
Los autores afirman que su trabajo proporciona una base matemática rigurosa para comprender qué patrones de decaimiento temporal en la señalización quimiotáctica admiten reducciones de similitud. El estudio demuestra que:

• El decaimiento dependiente del tiempo no es meramente una regularización matemática, sino que codifica estrategias de control biológicamente significativas (por ejemplo, regulación enzimática).
• Leyes de decaimiento específicas (constante, ley de potencia, exponencial) corresponden a regímenes dinámicos distintos: frentes viajeros acotados, redistribución autosimilar y supresión o amplificación rápida de señales, respectivamente.
• El enfoque basado en simetrías permite avances analíticos en modelos limitados por flujo con mecanismos de degradación realistas y variables en el tiempo, ofreciendo soluciones de referencia para validar esquemas numéricos en dimensiones superiores.

Limitaciones y Direcciones Futuras
El artículo nota modestamente que la reducción autosimilar para el caso de decaimiento de ley de potencia conduce a un sistema integro-diferencial que no admite una reducción puramente local a EDO bajo funciones de limitación de flujo genéricas. En consecuencia, las fórmulas de solución explícitas para este caso se interpretan como representaciones formales que requieren suposiciones adicionales de cierre o integración numérica. Los autores identifican un análisis completo de este caso específico como una dirección abierta para futuras investigaciones.