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

Resumo Técnico: Teoria de campo autoconsistente para nemáticos semiflexíveis

Definição do Problema
Embora a resposta linear de cristais líquidos termotrópicos de pequenas moléculas seja bem compreendida, o comportamento de cristais líquidos liotrópicos — que compreendem polímeros, agregados micelares e biofilamentos — apresenta desafios teóricos e experimentais significativos. Uma questão central é a origem da grande anisotropia elástica observada em sistemas liotrópicos e o acoplamento entre a densidade das cadeias e a ordem nemática. As teorias microscópicas existentes, amplamente baseadas no modelo de Onsager, têm dificuldade em abordar simultaneamente a elasticidade, defeitos topológicos e o acoplamento densidade–ordem em nemáticos semiflexíveis espacialmente variáveis. Além disso, as teorias de contínuo macroscópicas (Oseen–Frank e Landau–de Gennes) tipicamente introduzem coeficientes fenomenológicos sem uma conexão explícita com a estatística molecular subjacente, deixando inexploradas as origens microscópicas da anisotropia elástica e das estruturas de núcleo de defeitos.

Metodologia
Os autores empregam um arcabouço de Teoria de Campo Autoconsistente (SCFT) modelando nemáticos semiflexíveis como cadeias de tipo wormlike chain. O sistema é definido por um Hamiltoniano que incorpora:

1. Energia de flexão: Governada pelo comprimento de persistência ($l_p$) e pelo comprimento de contorno ($L$).
2. Interações de Maier–Saupe (MS): Interações de orientação anisotrópica ($u_2$) que impulsionam o ordenamento nemático.
3. Volume excluído isotrópico: Uma interação mediada pelo solvente ($u_0$) representando condições liotrópicas (tratando implicitamente os efeitos do solvente via parâmetros de Flory–Huggins).

O estudo procede em três estágios:

• Análise do Estado Uniforme: Os autores derivam um funcional de energia livre efetivo dependente da densidade do segmento da cadeia ($\rho$) e do parâmetro de ordem uniaxial ($S$). Isso permite a determinação das fronteiras de coexistência isotrópico-nemática (binodal) e dos limites de estabilidade (espinhodal).
• Análise de Resposta Linear: Campos externos fracos são introduzidos para sondar a resposta do sistema a perturbações na densidade, na ordem uniaxial e na orientação do diretor. Isso resulta na matriz de suscetibilidade inversa, a partir da qual as constantes elásticas de Frank ($K_1, K_2, K_3$) e os comprimentos de correlação são extraídos.
• Simulação de Estrutura de Defeitos: As equações completas de SCFT são resolvidas numericamente para disclinações +1/2 para analisar perfis de diretor e tamanhos de núcleo, comparando limites de cadeias rígidas com regimes semiflexíveis.

Principais Contribuições e Resultados

• Acoplamento Densidade–Nemático e Estabilidade:
  A energia livre efetiva derivada revela um forte acoplamento entre a densidade ($\rho$) e o grau de ordem ($S$). Os autores demonstram que o forte acoplamento densidade–nemático pode desestabilizar o estado nemático. A localização da instabilidade espinhodal depende criticamente da razão entre o volume excluído e a força de interação nemática ($u_0/u_2$). Em contraste, distorções do diretor acoplam-se à densidade e às variações de ordem apenas em ordens superiores, permanecendo efetivamente desacopladas no regime de resposta linear.

• Anisotropia Elástica e Flexibilidade:
  O estudo quantifica como a flexibilidade da cadeia ($L/l_p$) altera as constantes elásticas de Frank.

  • Cadeias Rígidas: $K_3$ (flexão/bend) aumenta rapidamente com a força de alinhamento, enquanto $K_1$ (splay) e $K_2$ (torção/twist) permanecem comparáveis.
  • Cadeias Semiflexíveis: O aumento da flexibilidade suprime fortemente a elasticidade de torção ($K_2$) e de flexão ($K_3$), enquanto a elasticidade de splay ($K_1$) é afetada comparativamente de forma fraca.
  • Crossover: Essa supressão diferencial leva a um crossover de uma elasticidade dominada por flexão (bend) para uma dominada por splay conforme a flexibilidade aumenta, evidenciada por uma reversão na razão $K_1/K_3$ (de $<1$ para $>1$).

• Comprimentos de Correlação e Anisotropia:
  A análise de resposta linear revela que, embora os coeficientes de gradiente diagonal da suscetibilidade inversa sugiram comprimentos de correlação mais curtos paralelamente ao diretor ($K_{\parallel} < K_{\perp}$), o acoplamento entre densidade e ordem renormaliza os coeficientes efetivos. Isso resulta em uma anisotropia onde o comprimento de correlação paralelo ao diretor é, de fato, maior ($\xi_{\parallel} > \xi_{\perp}$). Este achado alinha-se com observações experimentais de interfaces isotrópico-nemáticas e estruturas de núcleo de defeitos.

• Tamanhos de Núcleo de Defeito:
  Simulações numéricas de disclinações +1/2 mostram que o tamanho do núcleo é anisotrópico, com o raio sendo maior ao longo da direção onde o gradiente de densidade é paralelo ao diretor. O estudo estabelece que o tamanho característico do núcleo de defeito é proporcional ao comprimento de correlação de equilíbrio ($\xi_Q$), consistente com a escala de Landau–de Gennes. Os autores observam que, embora os cromônicos liotrópicos frequentemente exibam núcleos anormalmente grandes (potencialmente devido à cisão de agregados), seu modelo atual, que exclui a polidispersidade, captura a relação de escala entre o tamanho do núcleo e o comprimento de correlação impulsionada pelas flutuações de densidade.

• Comportamento Espinhodal:
  À medida que o sistema se aproxima do espinhodal nemático (controlado pela variação de $u_0$), os comprimentos de correlação divergem seguindo a escala $\xi \sim (u_0 - u_{0,i})^{-1/2}$. No entanto, em configurações de defeitos, essa divergência é prevenida pela formação local de um domínio isotrópico (tactoide) no núcleo, que suprime a ordem nemática antes que o espinhodal seja alcançado.

Significância
O artigo fornece um arcabouço microscópico unificado que vincula a flexibilidade molecular, o acoplamento densidade–ordem nemática e a anisotropia elástica em nemáticos semiflexíveis. Ao derivar as constantes elásticas de Frank e os comprimentos de correlação diretamente da estatística das cadeias dentro do formalismo SCFT, o trabalho elucida as origens microscópicas da grande anisotropia elástica observada em sistemas liotrópicos. Os resultados validam o uso de teorias de contínuo (como Landau–de Gennes) para descrever estruturas de defeitos, ao mesmo tempo em que destacam a necessidade de incorporar flutuações de densidade e a flexibilidade das cadeias para prever com precisão as respostas elásticas e os limites de estabilidade. Os autores concluem que seu arcabouço oferece uma base para extensões futuras a cadeias auto-montadas, polidispersas e de agregação reversível, embora não proponham novos protocolos experimentais ou aplicações específicas além do escopo da atual análise teórica.