Topological delocalisation of confined 3D active nematics

Resumen Técnico: Deslocalización Topológica de Nemáticos Activos Confinados en 3D

Planteamiento del Problema
Los nemáticos activos son sistemas fuera del equilibrio caracterizados por la creación y aniquilación espontánea de defectos topológicos móviles, lo que impulsa fenómenos como la turbulencia activa. Si bien la dinámica de las líneas de disclinación en nemáticos activos 3D en volumen (bulk) se comprende bien, exhibiendo percolación y entrelazamiento similares a los de los fundidos de polímeros, el comportamiento de estos sistemas bajo confinamiento sigue siendo esquivo. Específicamente, la interacción entre las restricciones geométricas (curvatura) y la actividad interna en entornos 3D confinados no ha sido plenamente caracterizada. Estudios previos se han centrado mayoritariamente en películas 2D sobre superficies cerradas o sistemas en volumen, dejando un vacío en la comprensión de cómo las líneas de defectos conectan puntos de frontera e impulsan la dinámica en el volumen dentro de geometrías 3D cerradas con curvatura no uniforme.

Metodología
Los autores emplean simulaciones numéricas utilizando un enfoque híbrido de Lattice Boltzmann para resolver las ecuaciones de continuidad de la nematohidrodinámica activa.

• Configuración del Sistema: Las simulaciones modelan nemáticos activos 3D confinados dentro de cilindros cerrados (y esferocilindros) con anclaje degenerado plano, asegurando que el campo del director se alinee tangencialmente a la superficie.
• Modelo de Física: El sistema se describe mediante una densidad de energía libre de campo medio que involucra un campo de fase $\phi$ para el dominio cilíndrico y un parámetro de orden tensorial $Q$ para el cristal líquido. La dinámica está gobernada por la ecuación de Cahn-Hilliard para el campo de fase y la ecuación de Beris-Edwards para el tensor $Q$, acopladas con las ecuaciones de Navier-Stokes incompresibles para la velocidad del fluido.
• Parámetros: El estudio varía el radio del cilindro ($R$) y la fuerza de la actividad ($\zeta$). Se define un número de actividad adimensional, $A = R\sqrt{\zeta/K}$, que es la relación entre el tamaño de confinamiento y la escala de longitud activa $\ell_\zeta = \sqrt{K/\zeta}$.
• Análisis: Las líneas de defectos se identifican mediante el tensor de densidad de disclinación. Las métricas clave incluyen la fracción de ocupación de defectos (para distinguir estados localizados frente a deslocalizados), longitudes de contorno de defectos, patrones de giro (caracterizados por el ángulo de torsión $\beta$) y leyes de escala para el número y la longitud de los defectos.

Contribuciones Clave y Resultados

1. Competencia entre Elasticidad y Actividad:
   El estudio revela una competencia fundamental entre la elasticidad pasiva y la actividad.

   • Estado Localizado: A baja actividad (o en equilibrio), la elasticidad pasiva predomina. Los defectos se anclan en regiones de alta curvatura (los bordes/tapones del cilindro) para minimizar las deformaciones elásticas. En equilibrio, dos líneas de disclinación cortas conectan los extremos $+1/2$ cerca de los bordes.
   • Estado Deslocalizado: A medida que la actividad aumenta, las tensiones activas superan la barrera de localización elástica. Los defectos ganan motilidad, conduciendo a un régimen de deslocalización topológica donde los defectos se mueven caóticamente por todo el volumen del cilindro, asemejándose a la turbulencia activa.

2. Firmas de la Transición de Fase:
   La transición de estados localizados a deslocalizados está marcada por un número de actividad crítico $A_c$.

   • Fluctuaciones: La transición viene acompañada de fluctuaciones pronunciadas en la fracción de ocupación de defectos.
   • Mecanismos: El mecanismo de deslocalización depende del radio del cilindro. Para radios pequeños ($R \le 20$), la deslocalización ocurre mediante la nucleación de nuevos pares de defectos a lo largo de la superficie del cilindro (impulsada por inestabilidad). Para radios mayores, los defectos escapan de los tapones y atraviesan el cilindro sin necesidad de generar nuevos defectos inmediatamente, llevando eventualmente a un estado totalmente deslocalizado.

3. Morfología de Defectos y Giro:

   • Patrones de Giro: En el estado localizado, los defectos son disclinaciones puras de $+1/2$. En el estado deslocalizado, la actividad introduce una amplia gama de patrones de giro, incluyendo torsión ($\cos \beta \approx 0$) y estados intermedios.
   • Segregación de Carga: Incluso en el estado deslocalizado, persiste un efecto de "segregación de carga" cerca de los tapones. El carácter de giro local ($\cos \beta$) decae exponencialmente desde los tapones hacia el centro del cilindro, con una amplitud que escala como $A^{-2}$, pero con una longitud de decaimiento característica que permanece independiente de la actividad.

4. Leyes de Escala en Confinamiento vs. Volumen:
   El artículo establece comportamientos de escala distintos para sistemas confinados en comparación con el volumen (bulk):

   • Número de Defectos: El número promedio de líneas de defectos escala como $\langle N_d \rangle \sim \ell_\zeta^{-2}$, similar a los sistemas 2D, debido a que las líneas están restringidas a la superficie del cilindro. Esto contrasta con los bucles de defectos en el volumen (o confinamiento débil), que escalan como $\ell_\zeta^{-3}$.
   • Longitud de Defectos: A diferencia de los sistemas en volumen donde la longitud de los defectos escala linealmente con la escala de longitud activa ($L_d \sim \ell_\zeta$), el confinamiento introduce una longitud característica independiente de la actividad. La longitud promedio de los defectos escala con el radio del cilindro como $L_d \sim R^{5/3}$.
   • Analogía con Polímeros: La distribución de probabilidad de las longitudes de los defectos en confinamiento decae exponencialmente (a diferencia de la distribución de ley de potencia en volumen). El exponente de escala $R^{5/3}$ corresponde a un paseo aleatorio de autoevitación (exponente de Flory $\nu = 3/5$), lo que sugiere que las líneas de defectos confinadas se comportan de forma análoga a "polímeros vivos" dentro de un bloque (blob) confinado.

Significancia
Los autores afirman que este trabajo establece el confinamiento de la curvatura inhomogénea como un mecanismo versátil para controlar la dinámica topológica activa. Al ajustar el número de actividad y la geometría de la frontera, se puede controlar los sitios de anclaje de defectos, la formación de vórtices y la transición a la turbulencia deslocalizada. La identificación de una transición de fase reminiscente de los aislantes/superconductores topológicos (donde las líneas de flujo permanecen ancladas o escapan) y la demostración de leyes de escala tipo polímero para las líneas de defectos proporcionan un nuevo marco teórico para comprender la materia activa 3D. Los resultados sugieren que las restricciones geométricas pueden aprovecharse para estabilizar configuraciones de defectos específicas o inducir flujos espontáneos, con relevancia potencial para sistemas biológicos como la morfogénesis epitelial y los tejidos cancerosos donde existen nemáticos activos 3D.