Ordering and Defect Dynamics in Passive and Active Nematopolars

Resumo Técnico: Dinâmica de Ordenação e Defeitos em Nematopolares Passivos e Ativos

Enunciado do Problema
O artigo aborda a compreensão teórica do comportamento coletivo fora do equilíbrio em sistemas "nematopolares", onde coexistem interações polares (assimétricas cabeça-cauda) e nemáticas (simétricas cabeça-cauda). Tais sistemas são observados em diversos contextos biológicos e sintéticos, incluindo cristais líquidos ferroelétricos, cristais líquidos vivos, colônias bacterianas e misturas de microtúbulos-motores. Embora abordagens teóricas anteriores tenham utilizado modelos de dois campos (acoplando um vetor polar a um tensor nemático) ou focado em turbulência ativa, há uma falta de caracterização sistemática quanto à ordenação de fase e dinâmica de coalescência (coarsening) em modelos mínimos de campo único para sistemas nematopolares secos. Especificamente, a interação entre simetrias concorrentes e o papel da autopropulsão (atividade) na dinâmica de defeitos e nas leis de crescimento de domínios permanece inexplorada.

Metodologia
Os autores propõem e investigam numericamente um modelo "seco" minimalista de campo único para sistemas nematopolares.

• Modelo: O sistema é descrito por um campo vetorial polar $\mathbf{p}$ evoluindo via dinâmica puramente relaxacional (tipo Modelo A) no limite passivo. O funcional de energia livre inclui um potencial de poço duplo para a polarização volumétrica, um termo elástico polar ($\propto |\nabla \mathbf{p}|^2$) e um termo de tipo nemático derivado de um tensor construído a partir de $\mathbf{p}$ ($\hat{P} \propto \mathbf{p}\mathbf{p}^\top - |\mathbf{p}|^2 I/2$).
• Atividade: O comportamento fora do equilíbrio é introduzido via um termo de autoadvecção, $\Lambda(\mathbf{p} \cdot \nabla)\mathbf{p}$, na equação de movimento. Este termo permite que a polarização atue simultaneamente como um parâmetro de ordem e um campo de velocidade local, consistente com as teorias de Toner-Tu para matéria ativa seca.
• Simulação: Os autores utilizam esquemas de diferença finita em redes quadradas 2D ($N=512, 1024$) com condições de contorno periódicas.
• Observáveis: O estudo rastreia a localização, o número e o tipo de defeitos topológicos (cargas inteiras e semipolares/fracionárias), funções de correlação polar e nemática ($C_p, C_n$), e comprimentos de correlação característicos ($L_p, L_n$).
• Abordagem: A investigação procede em duas etapas: (1) simulações controladas de pares de defeitos isolados e loops para analisar mecanismos de relaxação elementares, e (2) simulações de larga escala partindo de condições iniciais desordenadas para estudar a dinâmica de coalescência global.

Principais Contribuições e Resultados

1. Morfologia de Defeitos e Interações Mediadas por Cordas (Caso Passivo):

   • No regime passivo ($\Lambda=0$) com uma competição equilibrada entre acoplamentos polares e nemáticos, o sistema desenvolve espontaneamente cordas de despolarização (linhas de polarização nula) conectando defeitos de carga semipolar ($\pm 1/2$) e loops de despolarização fechados.
   • Pares de Defeitos: Defeitos semipolares de cargas opostas conectados por uma corda exibem velocidades de aniquilação não monotônicas como função da razão de força relativa polar/nemática ($\bar{k} = k_p/k_n$). Isso decorre de uma competição entre a redução das constantes elásticas efetivas e o custo energético da própria corda.
   • Defeitos de Mesma Carga: Diferente de sistemas nemáticos puros, onde cargas iguais se repelem indefinidamente, defeitos semipolares de mesma carga neste modelo atingem uma separação de equilíbrio finita, determinada pelo equilíbrio entre as forças de repulsão do tipo Coulomb e a tensão da corda conectora.
   • Dinâmica de Loops: Loops de despolarização relaxam via dois mecanismos distintos dependendo dos parâmetros do sistema: (i) encolhimento (dirigido pela curvatura, dominante em regimes dominados por nemáticos) e (ii) relaxação rotacional, onde o campo de polarização dentro do loop se reorienta para coincidir com o exterior (dominante em regimes dominados por polares). Em regimes polares extremos, os loops podem romper-se em cordas de carga inteira.

2. Dinâmica de Coalescência (Caso Passivo):

   • Partindo de estados desordenados, o sistema passa por um processo de coalescência. A escala de comprimento característica $L(t)$ cresce como $L(t) \sim (t/\ln t)^{1/2}$.
   • Esta escala é consistente com sistemas que possuem simetria rotacional contínua e defeitos pontuais (análogos ao modelo XY 2D e sistemas nemáticos O(2)/Z2), confirmando que a presença de cordas e loops não altera a lei fundamental de escala dinâmica, embora diminua a taxa absoluta de ordenação em comparação com sistemas nemáticos puros.
   • O número total de defeitos decai como $n_{tot} \sim \ln t / t$.

3. Dinâmica Ativa e Coalescência Arrestada:

   • Ao introduzir a autoadvecção ($\Lambda \neq 0$), o sistema exibe uma Quebra de Simetria de Carga Induzida pela Motilidade (MICS). A simetria entre defeitos positivos e negativos é quebrada, levando à coexistência de defeitos inteiros positivos (tipo aster) e defeitos semipolares negativos (tipo trevo).
   • Defeitos inteiros negativos e defeitos semipolares positivos são dinamicamente suprimidos e se decompõem nas espécies estáveis.
   • Coalescência Arrestada: Atividade suficientemente forte leva a um estado estacionário onde a escala de comprimento característica satura (estagna) e a densidade de defeitos permanece finita. O sistema entra em um regime de "coalescência arrestada", onde a evolução dos defeitos é fortemente suprimida, impedindo que o sistema atinja um estado ordenado uniforme.

Significância
O artigo afirma fornecer um arcabouço unificado simples para entender a ordenação e a dinâmica de defeitos em sistemas nematopolares. Ao utilizar um modelo de campo único minimalista, os autores demonstram que:

• Estruturas complexas de defeitos (cordas, loops e coexistência de cargas mistas) podem emergir da competição entre termos de energia livre polar e nemática sem a necessidade de múltiplos parâmetros de ordem.
• A atividade, implementada via um termo de autoadvecção simples, é um mecanismo fundamental capaz de remodelar a ordenação mediada por defeitos, estabilizar estados fora do equilíbrio e induzir uma quebra de simetria que não possui paralelo em sistemas passivos.
• Os resultados oferecem uma base teórica para interpretar o comportamento de diversos sistemas biológicos e sintéticos, como cristais líquidos ferroelétricos, cristais líquidos vivos e misturas de microtúbulos-motores, onde ambos os tipos de ordem emergem dos mesmos graus de liberdade.

Os autores observam explicitamente que, embora seu modelo seja minimalista, ele captura a fenomenologia essencial desses sistemas e desvincula os papéis das simetrias concorrentes e da atividade. Eles sugerem que este arcabouço serve como ponto de partida para investigar classes mais amplas de materiais ativos com ordens orientacionais concorrentes.