Eppur si muove: Shape of topological defects -- and consequent motion -- in active nematics

Resumen Técnico

1. Planteamiento del Problema

Los defectos topológicos en sistemas con orden de cristal líquido son fundamentales para determinar sus propiedades a gran escala. En sistemas activos (donde el consumo de energía genera fuerzas internas), estos defectos exhiben comportamientos imposibles en equilibrio, como la autopropulsión de los defectos de carga $+1/2$ en nemáticos ordenados.

El problema central abordado en este trabajo es la inconsistencia en las predicciones teóricas existentes sobre la velocidad de autopropulsión de estos defectos. La literatura previa ha utilizado frecuentemente la aproximación de que la forma del defecto permanece inalterada por la actividad, asumiendo que conserva la misma forma que en cristales líquidos pasivos con una sola constante de Frank (donde las deformaciones de splay y bend tienen el mismo costo elástico).

• La contradicción: Bajo esta suposición, ciertos modelos teóricos predicen que los defectos $+1/2$ no deberían autopropulsarse en ciertos regímenes de actividad, o incluso que podrían autopropulsarse en sistemas pasivos en equilibrio (lo cual viola la conservación de la energía).
• La hipótesis: Los autores proponen que la actividad induce modificaciones en la forma del defecto (tanto en el núcleo como en el campo lejano) y que ignorar estas deformaciones conduce a conclusiones erróneas, incluso para actividades infinitesimales.

2. Metodología

Los autores emplean una combinación de teoría analítica perturbativa y simulaciones numéricas directas:

• Modelo Teórico: Utilizan un modelo genérico de nemáticos activos secos (sin acoplamiento hidrodinámico explícito con un fluido, pero con términos de actividad no lineales). La dinámica del tensor de orden nemático $Q_{ij}$ se describe mediante una ecuación de movimiento que incluye términos elásticos (con dos constantes de Frank, $K_1$ y $K_2$) y términos activos (coeficientes $L_1, L_2, K_1, K_2$).
• Desarrollo Analítico:
  • Realizan una expansión perturbativa tanto en el núcleo del defecto ($r \to 0$) como en el campo lejano ($r \to \infty$).
  • Para defectos de carga $-1/2$, derivan expresiones analíticas cerradas para la forma del defecto.
  • Para defectos de carga $+1/2$, asumen un estado estacionario en movimiento y derivan una relación directa entre la modificación del campo angular en el núcleo y la velocidad de autopropulsión.
  • Utilizan el formalismo de Halperin-Mazenko para calcular velocidades, pero corrigiendo la forma del defecto para que sea consistente con la energía libre del sistema (en lugar de imponer la forma de un sistema de una sola constante de Frank).
• Simulaciones Numéricas:
  • Implementan un solver de elementos finitos en MATLAB para resolver las ecuaciones de movimiento en un dominio circular.
  • Simulan tanto sistemas pasivos (límites de una y dos constantes de Frank) como sistemas activos con diferentes coeficientes de no linealidad.
  • Rastrean la posición del defecto en el tiempo para medir su velocidad y extraen los perfiles de orden ($S$) y ángulo ($\theta$) para verificar la forma.

3. Contribuciones Clave

1. Primera descripción analítica de la forma de defectos en nemáticos activos: Demuestran que la actividad modifica intrínsecamente la estructura espacial de los defectos, rompiendo la simetría asumida en modelos anteriores.
2. Refutación de la aproximación de forma inalterada: Muestran que asumir la forma de un sistema pasivo de una sola constante de Frank en sistemas activos (o incluso pasivos de dos constantes) lleva a predicciones físicamente incorrectas, como velocidades de defectos no nulas en equilibrio.
3. Relación forma-velocidad para defectos $+1/2$: Establecen una conexión analítica explícita entre la deformación del campo angular en el núcleo del defecto y su velocidad de autopropulsión.
4. Validación universal de la autopropulsión: Demuestran que los defectos $+1/2$ en nemáticos activos siempre se autopropulsan, independientemente de la naturaleza específica de la actividad (si los coeficientes $L_2, K_1, K_2, L_1$ son cero o no), corrigiendo predicciones previas que sugerían lo contrario.

4. Resultados Principales

• Defectos $-1/2$:

  • En el núcleo, la forma del campo de orden escalar $S$ y el campo angular $\theta$ permanece idéntica a la de un sistema pasivo de una sola constante de Frank (la corrección angular en el núcleo es cero).
  • Sin embargo, en el campo lejano, la actividad induce una deformación angular significativa descrita por funciones $A(\phi)$ y $B(\phi)$ que dependen de los coeficientes activos y elásticos. Esta deformación es claramente distinta de la de los sistemas pasivos.
  • Los defectos $-1/2$ no se autopropulsan, pero su forma alterada es crucial para calcular interacciones defecto-defecto.

• Defectos $+1/2$:

  • Velocidad: La velocidad de autopropulsión $v_d$ es no nula en todos los casos activos. La ecuación perturbativa correcta (Eq. 5 del artículo) predice:
    $$v_d \log\left(\frac{3.29}{v_d}\right) \simeq 0.07(K_1 + K_2) - 0.78L_2 + 0.64L_1$$
    Esto contrasta con la predicción errónea anterior (Eq. 4) que sugería $v_d = 0$ si $L_2 = 0$.
  • Forma del Núcleo: La modificación del campo angular en el núcleo, $b(\phi)$, está directamente ligada a la velocidad:
    $$b(\phi) = \frac{1}{4}(v_d + a_0 L_2) \sin \phi$$
    Esto demuestra que la forma del defecto y su velocidad son variables acopladas que deben resolverse simultáneamente.
  • Simulaciones: Las simulaciones confirman que la velocidad es independiente del tamaño del sistema (a diferencia de los sistemas pasivos donde el "drift" disminuye con el tamaño) y que la forma analítica predicha coincide perfectamente con los datos numéricos.

• Sistemas Pasivos de Dos Constantes de Frank:

  • El trabajo también aclara un error histórico: si se calcula la velocidad de un defecto en un sistema pasivo de dos constantes de Frank asumiendo la forma incorrecta (de una sola constante), se obtiene una velocidad no nula. Al usar la forma correcta, la velocidad se anula, respetando la termodinámica de equilibrio.

5. Significado e Impacto

• Corrección de Paradigma: El artículo establece que la forma del defecto no es un parámetro fijo en sistemas activos, sino una variable dinámica que responde a la actividad. Ignorar esto invalida teorías previas sobre la dinámica de defectos.
• Implicaciones Biológicas y Físicas: Dado que los defectos topológicos en tejidos biológicos (como colonias bacterianas o epitelios) y materiales activos artificiales impulsan procesos morfológicos, flujos caóticos (turbulencia activa) y la muerte celular, entender su forma real es crucial para modelar estos fenómenos con precisión.
• Interacciones Defecto-Defecto: La forma alterada del defecto modifica las condiciones de contorno que impone a otros defectos. Esto sugiere que las interacciones entre defectos en sistemas activos pueden ser no recíprocas y dependientes de la actividad, un tema para futuras investigaciones.
• Extensibilidad: La metodología desarrollada es aplicable a sistemas "húmedos" (con flujo de fluido), a cristales líquidos polares y a otros estados de la materia activa (hexáticos, esmécticos), proporcionando una base teórica sólida para estudiar la nucleación y dinámica de defectos en contextos más complejos.

En conclusión, el trabajo demuestra que "Eppur si muove" (y sin embargo se mueve) no solo se refiere al movimiento de los defectos, sino a la necesidad de reconsiderar su forma fundamental para entender correctamente la física de los sistemas activos.