Effect of Weak Non-Conservative Dynamics on Pattern Formation in Scalar Active Matter

Resumen Técnico: Efecto de la Dinámica No Conservativa Débil en la Formación de Patrones en Materia Activa Escalar

Planteamiento del Problema
Los sistemas de materia activa, que abarcan desde colonias bacterianas hasta coloides activos sintéticos, son inherentemente impulsados lejos del equilibrio termodinámico mediante el consumo continuo de energía. Mientras que los marcos teóricos estándar como el Modelo Activo B (AMB) describen la separación de fases en materia activa escalar bajo una estricta conservación de la masa, muchos sistemas biológicos y sintéticos exhiben violaciones débiles de la conservación de la masa debido a procesos de nacimiento-muerte de partículas, crecimiento, degradación o intercambio químico con reservorios. El estudio específico de cómo estas dinámicas no conservativas débiles modifican la evolución no lineal y las estructuras de estado estacionario de la separación de fases activa permanece poco explorado. Este trabajo aborda la brecha en la comprensión de cómo estas reacciones no conservativas modifican la evolución no lineal y las estructuras de estado estacionario de la separación de fases activa.

Metodología
El autor construye un modelo continuo mínimo extendiendo el Modelo Activo B con un término de reacción no conservativo débil. El sistema se describe mediante un campo de parámetro de orden escalar coarse-grained, $\phi(\mathbf{r}, t)$, que representa la densidad de partículas. La dinámica está gobernada por la ecuación de movimiento:
$$\frac{\partial \phi}{\partial t} = -\Gamma_\phi \nabla \cdot \mathbf{J} + \gamma \phi \left(1 - \frac{\phi}{\phi_c}\right)\left(1 + \frac{\phi}{\phi_c}\right)$$
donde $\mathbf{J} = -\nabla \mu + \sqrt{2\Gamma_\phi \phi}\eta$ es la corriente, y el potencial químico $\mu$ incluye una contribución estándar de energía de Landau y un término activo $\mu_{active} = \zeta |\nabla \phi|^2$ que rompe la simetría de inversión temporal. El parámetro $\zeta$ controla la fuerza de la actividad, mientras que $\gamma$ (con $\gamma \leq 0$) controla la fuerza de la reacción no conservativa.

El estudio emplea dos enfoques principales:

1. Análisis de Estabilidad Lineal: Para determinar la inestabilidad inicial y la relación de dispersión, analizando cómo el término de reacción y la actividad afectan la tasa de crecimiento de las perturbaciones alrededor de un estado estacionario homogéneo.
2. Simulaciones Numéricas: Simulaciones extensas en una red de $512 \times 512$ con condiciones de contorno periódicas para observar las morfologías de tiempo tardío. El sistema se inicializa con una densidad uniforme y evoluciona hacia estados estacionarios.
3. Caracterización Estructural: Análisis cuantitativo utilizando funciones de correlación de dos puntos de tiempo igual $C(r)$, factores de estructura estática $S(k)$, estadísticas de cúmulos (fracción del cúmulo más grande $L_c$, número de cúmulos aislados $N_c$, circularidad $C$), distribuciones de tamaño de gota y el parámetro de orden hexático de orientación de enlaces $\psi_6$.

Contribuciones Clave y Resultados

• Detención del Coarsening y Microseparación de Fases: El hallazgo principal es que incluso las reacciones no conservativas débiles ($\gamma < 0$) detienen fundamentalmente el proceso de coarsening (engrosamiento) típico de los sistemas activos conservados. En lugar de evolucionar hacia una separación de fases macroscópica, el sistema se estabiliza en estados de microseparación de fases fuera del equilibrio con una escala de longitud intrínseca bien definida.
• Reversión de las Transiciones Morfológicas: En el límite conservado ($\gamma = 0$), el aumento de la magnitud de la actividad negativa ($\zeta$) impulsa una transición de gotas a estructuras tipo gusano y, finalmente, a patrones laberínticos interconectados (descomposición espinodal). En contraste, la presencia de una reacción no conservativa débil revierte esta secuencia. A medida que la actividad aumenta, el sistema transiciona de patrones laberínticos a estructuras tipo gusano y, eventualmente, a gotas aisladas y polidispersas.
• Roles Distintos de los Parámetros:
  • Actividad ($\zeta$): Controla predominantemente la morfología global de los dominios y la longitud de onda característica de los patrones. El aumento de $|\zeta|$ reduce la escala de longitud característica ($\lambda^*$) e impulsa la transición de laberintos a gotas aisladas.
  • Coeficiente de Reacción ($\gamma$): Promueve principalmente la microseparación de fases y mejora el orden hexagonal local. El aumento de la magnitud de $\gamma$ suprime la coalescencia de dominios, aumenta el número de cúmulos aislados y reduce el tamaño promedio de las gotas, llevando a una distribución de tamaño más estrecha.
• Perspectivas de la Estabilidad Lineal: El análisis revela que el término de reacción entra en la relación de dispersión como una constante independiente de $q$, desplazando la curva de tasa de crecimiento verticalmente sin alterar la posición del número de onda más inestable $q^*$. La contribución activa ($\zeta |\nabla \phi|^2$) no aparece en el orden lineal, lo que significa que ninguno de los dos términos modifica la selección de la longitud de onda inicial. Las diferencias cualitativas en la morfología surgen de la evolución no lineal que sigue a la inestabilidad inicial.
• Emergencia de Orden Hexático: En la fase de gotas, el sistema exhibe un empaquetamiento hexagonal local. El término de reacción no conservativo potencia este orden hexático local ($\psi_6 \to 1$) al regular la distribución de tamaño de las gotas y facilitar la homogeneización de tamaño, lo que permite un empaquetamiento ordenado. Esto es distinto del maduramiento de Ostwald clásico, ya que la naturaleza no conservada del parámetro de orden impide el crecimiento ilimitado.
• Comportamiento de Escalamiento: Las funciones de correlación escaladas $C(r/\lambda^*)$ colapsan en una curva maestra universal para distintos valores de $\zeta$ y $\gamma$, indicando que, si bien los parámetros reescalan la longitud característica, la organización espacial normalizada subyacente permanece robusta.

Significancia
Este artículo demuestra que las violaciones débiles de la conservación de la masa no son meramente perturbaciones, sino que alteran fundamentalmente la dinámica de coarsening no lineal de la separación de fases activa. Al introducir un marco mínimo que acopla el transporte activo con reacciones no conservativas débiles, el estudio proporciona un mecanismo para estabilizar estados de microseparación de fases que, de otro modo, sufrirían coarsening hacia fases masivas. Los resultados sugieren que la dinámica no conservativa sirve como un mecanismo simple pero efectivo para controlar la autoorganización fuera del equilibrio. Este marco ofrece una base para comprender la formación de patrones en sistemas biológicos (por ejemplo, colonias bacterianas con ciclos de nacimiento-muerte) y para diseñar materiales activos sintéticos con microestructuras ajustables, tales como gotas aisladas con un ordenamiento específico, sin requerir un forzamiento externo complejo.