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

Resumo Técnico: Efeito de Dinâmicas Não Conservativas Fracas na Formação de Padrões em Matéria Ativa Escalar

Problema
Sistemas de matéria ativa, que variam de colônias bacterianas a coloides ativos sintéticos, são inerentemente impulsionados longe do equilíbrio termodinâmico por meio do consumo contínuo de energia. Embora os arcabouços teóricos padrão, como o Modelo Ativo B (AMB), descrevam a separação de fases em matéria ativa escalar sob conservação estrita de massa, muitos sistemas biológicos e sintéticos exibem violações fracas da conservação de massa devido a processos de nascimento-morte de partículas, crescimento, degradação ou troca química com reservatórios. A influência específica dessas dinâmicas não conservativas fracas sobre a morfologia de tempo tardio, a seleção de padrões e o comportamento de coalescência (coarsening) em sistemas ativos escalares permanece subexplorada. Este trabalho aborda a lacuna na compreensão de como reações não conservativas modificam a evolução não linear e as estruturas de estado estacionário da separação de fases ativa.

Metodologia
O autor constrói um modelo contínuo mínimo estendendo o Modelo Ativo B com um termo de reação não conservativa fraca. O sistema é descrito por um campo de parâmetro de ordem escalar coarse-grained, $\phi(\mathbf{r}, t)$, representando a densidade de partículas. A dinâmica é governada pela equação de movimento:
$$\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)$$
onde $\mathbf{J} = -\nabla \mu + \sqrt{2\Gamma_\phi \phi}\eta$ é a corrente, e o potencial químico $\mu$ inclui uma contribuição padrão de energia de Landau e um termo ativo $\mu_{active} = \zeta |\nabla \phi|^2$ que quebra a simetria de reversão temporal. O parâmetro $\zeta$ controla a força da atividade, enquanto $\gamma$ (com $\gamma \leq 0$) controla a força da reação não conservativa.

O estudo emprega duas abordagens primárias:

1. Análise de Estabilidade Linear: Para determinar a instabilidade inicial e a relação de dispersão, analisando como o termo de reação e a atividade afetam a taxa de crescimento de perturbações em torno de um estado estacionário homogêneo.
2. Simulações Numéricas: Simulações extensivas em uma rede $512 \times 512$ com condições de contorno periódicas para observar morfologias de tempo tardio. O sistema é inicializado com uma densidade uniforme e evoluído para estados estacionários.
3. Caracterização Estrutural: Análise quantitativa usando funções de correlação de dois pontos de tempo igual $C(r)$, fatores de estrutura estáticos $S(k)$, estatísticas de cluster (fração do maior cluster $L_c$, número de clusters isolados $N_c$, circularidade $C$), distribuições de tamanho de gota e o parâmetro de ordem hexática de orientação de ligação $\psi_6$.

Principais Contribuições e Resultados

• Interrupção da Coalescência e Microfase de Separação: O principal achado é que mesmo reações não conservativas fracas ($\gamma < 0$) interrompem fundamentalmente o processo de coalescência típico de sistemas ativos conservados. Em vez de evoluírem para uma separação de fase em massa (bulk), o sistema se estabiliza em estados de microfase separada fora do equilíbrio com uma escala de comprimento intrínseca bem definida.
• Reversão de Transições Morfológicas: No limite conservado ($\gamma = 0$), o aumento da magnitude da atividade negativa ($\zeta$) conduz a uma transição de gotas para estruturas do tipo verme (worm-like) e, finalmente, para padrões labirínticos interconectados (decomposição espinodal). Em contraste, a presença de uma reação não conservativa fraca inverte essa sequência. À medida que a atividade aumenta, o sistema transita de padrões labirínticos para estruturas do tipo verme e, eventualmente, para gotas isoladas e polidispersas.
• Papéis Distintos dos Parâmetros:
  • Atividade ($\zeta$): Controla predominantemente a morfologia global dos domínios e o comprimento de onda característico dos padrões. O aumento de $|\zeta|$ reduz a escala de comprimento característica ($\lambda^*$) e impulsiona a transição de labirintos para gotas isoladas.
  • Coeficiente de Reação ($\gamma$): Promove primariamente a microfase de separação e aumenta o ordenamento hexagonal local. O aumento da magnitude de $\gamma$ suprime a coalescência de domínios, aumenta o número de clusters isolados e reduz o tamanho médio das gotas, levando a uma distribuição de tamanho mais estreita.
• Insights de Estabilidade Linear: A análise revela que o termo de reação entra na relação de dispersão como uma constante independente de $q$, deslocando a curva de taxa de crescimento verticalmente sem alterar a posição do número de onda mais instável $q^*$. A contribuição ativa ($\zeta |\nabla \phi|^2$) não aparece na ordem linear, o que significa que nenhum dos termos modifica a seleção do comprimento de onda inicial. As diferenças qualitativas na morfologia surgem da evolução não linear após a instabilidade inicial.
• Emergência de Ordem Hexática: Na fase de gotas, o sistema exibe empacotamento hexagonal local. O termo de reação não conservativa aumenta esse ordenamento hexático local ($\psi_6 \to 1$) ao regular a distribuição de tamanho das gotas e facilitar a homogeneização de tamanho, o que permite um empacotamento ordenado. Isso é distinto do amadurecimento de Ostwald clássico, pois a natureza não conservada do parâmetro de ordem impede o crescimento ilimitado.
• Comportamento de Escala: Funções de correlação escalonadas $C(r/\lambda^*)$ colapsam em uma curva mestra universal para variações de $\zeta$ e $\gamma$, indicando que, embora os parâmetros reescalonem o comprimento característico, a organização espacial normalizada subjacente permanece robusta.

Significância
O artigo demonstra que violações fracas da conservação de massa não são meramente perturbações, mas alteram fundamentalmente a dinâmica de coalescência não linear da separação de fases ativa. Ao introduzir um arcabouço minimalista que acopla o transporte ativo com reações não conservativas fracas, o estudo fornece um mecanismo para estabilizar estados de microfase separada que, de outra forma, coalesceriam em fases de massa. Os resultados sugerem que a dinâmica não conservativa serve como um mecanismo simples, porém eficaz, para controlar a auto-organização fora do equilíbrio. Este arcabouço oferece uma base para compreender a formação de padrões em sistemas biológicos (ex: colônias bacterianas com ciclos de nascimento-morte) e para projetar materiais ativos sintéticos com microestruturas ajustáveis, como gotas isoladas com ordenamento específico, sem a necessidade de forças externas complexas.