Ordering and Defect Dynamics in Passive and Active Nematopolars

Resumen Técnico: Dinámica de Ordenamiento y Defectos en Nematopolares Pasivos y Activos

Planteamiento del Problema
El artículo aborda la comprensión teórica del comportamiento colectivo fuera del equilibrio en sistemas "nematopolares", donde coexisten interacciones polares (asimetría cabeza-cola) y nemáticas (simetría cabeza-cola). Tales sistemas se observan en diversos contextos biológicos y sintéticos, incluyendo cristales líquidos ferroeléctricos, cristales líquidos vivos, colonias bacterianas y mezclas de microtúbulos y motores. Si bien los enfoques teóricos previos han utilizado modelos de dos campos (acoplando un vector polar a un tensor nemático) o se han centrado en la turbulencia activa, existe una falta de caracterización sistemática respecto al ordenamiento de fase y la dinámica de engrosamiento en modelos mínimos de un solo campo para sistemas nematopolares secos. Específicamente, la interacción entre simetrías competidoras y el papel de la autopropulsión (actividad) en la dinámica de defectos y las leyes de crecimiento de dominios permanece inexplorado.

Metodología
Los autores proponen e investigan numéricamente un modelo mínimo de un solo campo "seco" para sistemas nematopolares.

• Modelo: El sistema se describe mediante un único campo vectorial polar $\mathbf{p}$ que evoluciona a través de una dinámica puramente relajacional (tipo Modelo A) en el límite pasivo. El funcional de energía libre incluye un potencial de doble pozo para la polarización de volumen, un término elástico polar ($\propto |\nabla \mathbf{p}|^2$) y un término de tipo nemático derivado de un tensor construido a partir de $\mathbf{p}$ ($\hat{P} \propto \mathbf{p}\mathbf{p}^\top - |\mathbf{p}|^2 I/2$).
• Actividad: El comportamiento fuera del equilibrio se introduce mediante un término de autoadvección, $\Lambda(\mathbf{p} \cdot \nabla)\mathbf{p}$, en la ecuación de movimiento. Este término permite que la polarización actúe simultáneamente como un parámetro de orden y como un campo de velocidad local, de forma consistente con las teorías de Toner-Tu para materia activa seca.
• Simulación: Los autores emplean esquemas de diferencias finitas en redes cuadradas 2D ($N=512, 1024$) con condiciones de contorno periódicas.
• Observables: El estudio rastrea la ubicación, el número y el tipo de defectos topológicos (cargas enteras y de semienteras), las funciones de correlación polar y nemática ($C_p, C_n$), y las longitudes de correlación características ($L_p, L_n$).
• Enfoque: La investigación procede en dos etapas: (1) simulaciones controladas de pares de defectos aislados y bucles para analizar los mecanismos de relajación elementales, y (2) simulaciones a gran escala partiendo de condiciones iniciales desordenadas para estudiar la dinámica de engrosamiento global.

Contribuciones Clave y Resultados

1. Morfología de Defectos e Interacciones Mediadas por Cuerdas (Caso Pasivo):

   • En el régimen pasivo ($\Lambda=0$) con una competencia equilibrada entre acoplamientos polares y nemáticos, el sistema desarrolla espontáneamente cuerdas de despolarización (líneas de polarización nula) que conectan defectos de semientera ($\pm 1/2$) y bucles de despolarización cerrados.
   • Pares de Defectos: Los defectos de semientera de carga opuesta conectados por una cuerda exhiben velocidades de aniquilación no monotónicas en función de la relación de fuerza relativa polar/nemática ($\bar{k} = k_p/k_n$). Esto surge de una competencia entre la reducción de las constantes elásticas efectivas y el costo energético de la cuerda misma.
   • Defectos de la Misma Carga: A diferencia de los sistemas nemáticos puros donde las cargas iguales se repelen indefinidamente, los defectos de semientera de la misma carga en este modelo alcanzan una separación de equilibrio finita determinada por el equilibrio entre las fuerzas de tipo Coulomb repulsivas y la tensión de la cuerda conectora.
   • Dinámica de Bucles: Los bucles de despolarización se relajan mediante dos mecanismos distintos dependiendo de los parámetros del sistema: (i) encogimiento (impulsado por la curvatura, dominante en regímenes dominados por lo nemático) y (ii) relajación rotacional, donde el campo de polarización dentro del bucle se reorienta para coincidir con el exterior (dominante en regímenes dominados por lo polar). En regímenes polares extremos, los bucles pueden romperse en cuerdas de carga entera.

2. Dinámica de Engrosamiento (Caso Pasivo):

   • Partiendo de estados desordenados, el sistema experimenta un engrosamiento. La escala de longitud característica $L(t)$ crece como $L(t) \sim (t/\ln t)^{1/2}$.
   • Esta escala es consistente con sistemas que poseen simetría rotacional continua y defectos puntuales (análogos al modelo XY 2D y sistemas nemáticos O(2)/Z2), confirmando que la presencia de cuerdas y bucles no altera la ley de escala dinámica fundamental, aunque ralentiza la tasa absoluta de ordenamiento en comparación con los sistemas puramente nemáticos.
   • El número total de defectos decae como $n_{tot} \sim \ln t / t$.

3. Dinámica Activa y Engrosamiento Arrestado:

   • Al introducir la autoadvección ($\Lambda \neq 0$), el sistema exhibe una Ruptura de Simetría de Carga Inducida por Motilidad (MICS). La simetría entre defectos positivos y negativos se rompe, llevando a la coexistencia de defectos enteros positivos (tipo aster) y defectos de semientera negativos (tipo trébol).
   • Los defectos enteros negativos y los defectos de semientera positivos son dinámicamente suprimidos y se descomponen en las especies estables.
   • Engrosamiento Arrestado: Una actividad suficientemente fuerte conduce a un estado estacionario donde la escala de longitud característica se satura (meseta) y la densidad de defectos permanece finita. El sistema entra en un régimen de "engrosamiento arrestado" donde la evolución de los defectos es fuertemente suprimida, impidiendo que el sistema alcance un estado ordenado uniforme.

Significancia
El artículo afirma proporcionar un marco unificado simple para comprender el ordenamiento y la dinámica de defectos en sistemas nematopolares. Al utilizar un modelo minimalista de un solo campo, los autores demuestran que:

• Estructuras complejas de defectos (cuerdas, bucles y coexistencia de cargas mixtas) pueden emerger de la competencia entre los términos de energía libre polar y nemática sin requerir múltiples parámetros de orden.
• La actividad, implementada mediante un simple término de autoadvección, es un mecanismo fundamental capaz de remodelar el ordenamiento mediado por defectos, estabilizar estados fuera del equilibrio (engrosamiento arrestado) e inducir una ruptura de simetría que no tiene contraparte en los sistemas pasivos.
• Los resultados ofrecen una base teórica para interpretar el comportamiento de diversos sistemas biológicos y sintéticos, tales como cristales líquidos ferroeléctricos, cristales líquidos vivos y mezclas de microtúbulos y motores, donde ambos tipos de orden emergen de los mismos grados de libertad.

Los autores señalan explícitamente que, aunque su modelo es minimalista, captura la fenomenología esencial de estos sistemas y desvincula los roles de las simetrías competidoras y la actividad. Sugieren que este marco sirve como punto de partida para investigar clases más amplias de materiales activos con órdenes orientacionales en competencia.