Self-consistent field theory of semiflexible nematics: Density-nematic coupling, anisotropic elasticity, and defect core sizes

Resumen Técnico: Teoría de campo autoconsistente de nemáticos semiflexibles

Planteamiento del problema
Si bien la respuesta lineal de los cristales líquidos termotrópicos de moléculas pequeñas está bien comprendida, el comportamiento de los cristales líquidos liotrópicos —que comprenden polímeros, agregados micelares y biofilamentos— presenta desafíos teóricos y experimentales significativos. Un problema central es el origen de la gran anisotropía elástica observada en sistemas liotrópicos y el acoplamiento entre la densidad de la cadena y el orden nemático. Las teorías microscópicas existentes, basadas mayoritariamente en el modelo de Onsager, tienen dificultades para abordar simultáneamente la elasticidad, los defectos topológicos y el acoplamiento densidad-orden en nemáticos semiflexibles con variaciones espaciales. Además, las teorías de continuo macroscópicas (Oseen–Frank y Landau–de Gennes) introducen típicamente coeficientes fenomenológicos sin una conexión explícita con la estadística molecular subyacente, dejando inexplorados los orígenes microscópicos de la anisotropía elástica y las estructuras del núcleo de los defectos.

Metodología
Los autores emplean un marco de Teoría de Campo Autoconsistente (SCFT, por sus siglas en inglés) que modela los nemáticos semiflexibles como cadenas de tipo gusano (wormlike chains). El sistema se define mediante un Hamiltoniano que incorpora:

1. Energía de flexión: Gobernada por la longitud de persistencia ($l_p$) y la longitud de contorno ($L$).
2. Interacciones de Maier–Saupe (MS): Interacciones de orientación anisotrópica ($u_2$) que impulsan el orden nemático.
3. Volumen excluido isotrópico: Una interacción mediada por el solvente ($u_0$) que representa las condiciones liotrópicas (tratando implícitamente los efectos del solvente mediante parámetros de Flory–Huggins).

El estudio procede en tres etapas:

• Análisis del estado uniforme: Los autores derivan un funcional de energía libre efectivo dependiente de la densidad de los segmentos de la cadena ($\rho$) y del parámetro de orden uniaxial ($S$). Esto permite la determinación de los límites de coexistencia isótropo-nemático (binodal) y de estabilidad (espínodal).
• Análisis de respuesta lineal: Se introducen campos externos débiles para sondear la respuesta del sistema a perturbaciones en la densidad, el orden uniaxial y la orientación del director. Esto genera la matriz de susceptibilidad inversa, a partir de la cual se extraen las constantes elásticas de Frank ($K_1, K_2, K_3$) y las longitudes de correlación.
• Simulación de la estructura de defectos: Se resuelven numéricamente las ecuaciones completas de SCFT para disclinaciones de +1/2 para analizar los perfiles del director y los tamaños de núcleo, comparando los límites de cadenas rígidas con los regímenes semiflexibles.

Contribuciones clave y resultados

• Acoplamiento Densidad–Nemático y Estabilidad:
  La energía libre efectiva derivada revela un fuerte acoplamiento entre la densidad ($\rho$) y el grado de orden ($S$). Los autores demuestran que un fuerte acoplamiento densidad-nemático puede desestabilizar el estado nemático. La ubicación de la inestabilidad espínodal depende críticamente de la relación entre el volumen excluido y la fuerza de interacción nemática ($u_0/u_2$). En contraste, las distorsiones del director se acoplan a las variaciones de densidad y orden solo en órdenes superiores, permaneciendo efectivamente desacopladas en el régimen de respuesta lineal.

• Anisotropía Elástica y Flexibilidad:
  El estudio cuantifica cómo la flexibilidad de la cadena ($L/l_p$) altera las constantes elásticas de Frank.

  • Cadenas rígidas: $K_3$ (flexión o bend) aumenta rápidamente con la fuerza de alineación, mientras que $K_1$ (divergencia o splay) y $K_2$ (torsión o twist) permanecen comparables.
  • Cadenas semiflexibles: El aumento de la flexibilidad suprime fuertemente la elasticidad de torsión ($K_2$) y de flexión ($K_3$), mientras que la elasticidad de divergencia ($K_1$) se ve afectada comparativamente de forma débil.
  • Cruce (Crossover): Esta supresión diferencial conduce a un cruce de una elasticidad dominada por la flexión a una dominada por la divergencia a medida que aumenta la flexibilidad, evidenciado por una reversión en la relación $K_1/K_3$ (de $<1$ a $>1$).

• Longitudes de Correlación y Anisotropía:
  El análisis de respuesta lineal revela que, si bien los coeficientes de gradiente diagonales de la susceptibilidad inversa sugieren longitudes de correlación más cortas paralelas al director ($K_{\parallel} < K_{\perp}$), el acoplamiento entre la densidad y el orden renormaliza los coeficientes efectivos. Esto resulta en una anisotropía donde la longitud de correlación paralela al director es en realidad mayor ($\xi_{\parallel} > \xi_{\perp}$). Este hallazgo concuerda con las observaciones experimentales de interfaces isótropo-nemáticas y estructuras de núcleo de defectos.

• Tamaños del Núcleo de Defectos:
  Las simulaciones numéricas de disclinaciones de +1/2 muestran que el tamaño del núcleo es anisotrópico, siendo el radio mayor en la dirección donde el gradiente de densidad es paralelo al director. El estudio establece que el tamaño característico del núcleo del defecto es proporcional a la longitud de correlación de equilibrio ($\xi_Q$), consistente con la escala de Landau–de Gennes. Los autores señalan que, aunque los cromónicos liotrópicos a menudo exhiben núcleos anomalamente grandes (potencialmente debido a la escisión de agregados), su modelo actual, que excluye la polidispersidad, captura la relación de escala entre el tamaño del núcleo y la longitud de correlación impulsada por las fluctuaciones de densidad.

• Comportamiento Espínodal:
  A medida que el sistema se acerca al espínodal nemático (controlado mediante la variación de $u_0$), las longitudes de correlación divergen siguiendo la escala $\xi \sim (u_0 - u_{0,i})^{-1/2}$. Sin embargo, en configuraciones de defectos, esta divergencia es prevenida por la formación local de un dominio isótropo (tactoide) en el núcleo, lo que suprime el orden nemático antes de alcanzar el espínodal.

Significancia
El artículo proporciona un marco microscópico unificado que vincula la flexibilidad molecular, el acoplamiento densidad-orden nemático y la anisotropía elástica en nemáticos semiflexibles. Al derivar las constantes elásticas de Frank y las longitudes de correlación directamente de la estadística de las cadenas dentro del formalismo SCFT, el trabajo elucida los orígenes microscópicos de la gran anisotropía elástica observada en sistemas liotrópicos. Los resultados validan el uso de teorías de continuo (como Landau–de Gennes) para describir estructuras de defectos, al tiempo que resaltan la necesidad de incorporar las fluctuaciones de densidad y la flexibilidad de la cadena para predecir con precisión las respuestas elásticas y los límites de estabilidad. Los autores concluyen que su marco ofrece una base para futuras extensiones a cadenas autoensambladas, polidispersas y de agregación reversible, aunque no proponen protocolos experimentales o aplicaciones específicas más allá del alcance del presente análisis teórico.