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

Résumé technique : Théorie du champ auto-cohérent pour les nématiques semiflexibles

Énoncé du problème
Bien que la réponse linéaire des cristaux liquides thermotropes à petites molécules soit bien comprise, le comportement des cristaux liquides lyotropes — comprenant des polymères, des agrégats micellaires et des biofilaments — présente des défis théoriques et expérimentaux importants. Un problème central est l'origine de la grande anisotropie élastique observée dans les systèmes lyotropes et le couplage entre la densité des chaînes et l'ordre nématique. Les théories microscopiques existantes, largement basées sur le modèle d'Onsager, peinent à traiter simultanément l'élasticité, les défauts topologiques et le couplage densité-ordre dans les nématiques semiflexibles spatialement variables. De plus, les théories de continuum macroscopiques (Oseen–Frank et Landau–de Gennes) introduisent généralement des coefficients phénoménologiques sans connexion explicite avec les statistiques moléculaires sous-jacentes, laissant inexplorées les origines microscopiques de l'anisotropie élastique et des structures de cœur de défauts.

Méthodologie
Les auteurs utilisent un cadre de théorie du champ auto-cohérent (SCFT) modélisant les nématiques semiflexibles comme des chaînes de type « wormlike chain » (chaînes de type filament). Le système est défini par un hamiltonien incorporant :

1. L'énergie de flexion : régie par la longueur de persistance ($l_p$) et la longueur de contour ($L$).
2. Les interactions de Maier–Saupe (MS) : interactions d'orientation anisotropes ($u_2$) pilotant l'ordre nématique.
3. Le volume exclu isotrope : une interaction médiée par le solvant ($u_0$) représentant les conditions lyotropes (traitant implicitement les effets de solvant via les paramètres de Flory–Huggins).

L'étude se déroule en trois étapes :

• Analyse de l'état uniforme : Les auteurs dérivent une fonctionnelle d'énergie libre effective dépendant de la densité de segment de chaîne ($\rho$) et du paramètre d'ordre uniaxial ($S$). Cela permet de déterminer les frontières de coexistence isotrope-nématique (binodale) et les limites de stabilité (spinodale).
• Analyse de la réponse linéaire : Des champs externes faibles sont introduits pour sonder la réponse du système aux perturbations de densité, d'ordre uniaxial et d'orientation du directeur. Cela produit la matrice de susceptibilité inverse, à partir de laquelle les constantes élastiques de Frank ($K_1, K_2, K_3$) et les longueurs de corrélation sont extraites.
• Simulation de la structure des défauts : Les équations complètes de la SCFT sont résolues numériquement pour des disclinaisons +1/2 afin d'analyser les profils de directeur et les tailles de cœur, en comparant les limites de chaînes rigides avec les régimes semiflexibles.

Contributions clés et résultats

• Couplage Densité–Nématique et Stabilité :
  La fonctionnelle d'énergie libre effective dérivée révèle un fort couplage entre la densité ($\rho$) et le degré d'ordre ($S$). Les auteurs démontrent qu'un fort couplage densité-nématique peut déstabiliser l'état nématique. L'emplacement de l'instabilité spinodale dépend de manière critique du rapport entre le volume exclu et la force d'interaction nématique ($u_0/u_2$). En revanche, les distorsions du directeur se couplent à la densité et aux variations d'ordre uniquement à des ordres supérieurs, restant effectivement découplées dans le régime de réponse linéaire.

• Anisotropie Élastique et Flexibilité :
  L'étude quantifie comment la flexibilité des chaînes ($L/l_p$) modifie les constantes élastiques de Frank.

  • Chaînes rigides : $K_3$ (flexion/bend) augmente rapidement avec la force d'alignement, tandis que $K_1$ (évasement/splay) et $K_2$ (torsion/twist) restent comparables.
  • Chaînes semiflexibles : L'augmentation de la flexibilité supprime fortement l'élasticité de torsion ($K_2$) et de flexion ($K_3$), tandis que l'élasticité d'évasement ($K_1$) est affectée de manière relativement faible.
  • Crossover : Cette suppression différentielle conduit à un passage (crossover) d'une élasticité dominée par la flexion à une élasticité dominée par l'évasement à mesure que la flexibilité augmente, comme en témoigne l'inversion du rapport $K_1/K_3$ (de $<1$ à $>1$).

• Longueurs de Corrélation et Anisotropie :
  L'analyse de la réponse linéaire révèle que, bien que les coefficients de gradient diagonaux de la susceptibilité inverse suggèrent des longueurs de corrélation plus courtes parallèlement au directeur ($K_{\parallel} < K_{\perp}$), le couplage entre la densité et l'ordre renormalise les coefficients effectifs. Cela résulte en une anisotropie où la longueur de corrélation parallèle au directeur est en réalité plus grande ($\xi_{\parallel} > \xi_{\perp}$). Cette conclusion concorde avec les observations expérimentales des interfaces isotrope-nématique et des structures de cœur de défauts.

• Tailles de Cœur de Défaut :
  Les simulations numériques de disclinaisons +1/2 montrent que la taille du cœur est anisotrope, le rayon étant plus grand là où le gradient de densité est parallèle au directeur. L'étude établit que la taille caractéristique du cœur de défaut est proportionnelle à la longueur de corrélation d'équilibre ($\xi_Q$), ce qui est cohérent avec l'échelle de Landau–de Gennes. Les auteurs notent que si les chromoniques lyotropes présentent souvent des cœurs anormalement larges (potentiellement dus à la scission d'agrégats), leur modèle actuel, qui exclut la polydispersité, capture la relation d'échelle entre la taille du cœur et la longueur de corrélation pilotée par les fluctuations de densité.

• Comportement Spinodal :
  À mesure que le système approche du spinodal nématique (contrôlé par la variation de $u_0$), les longueurs de corrélation divergent selon la loi d'échelle $\xi \sim (u_0 - u_{0,i})^{-1/2}$. Cependant, dans les configurations de défauts, cette divergence est préemptée par la formation locale d'un domaine isotrope (tactoïde) au cœur, ce qui supprime l'ordre nématique avant que le spinodal ne soit atteint.

Signification
Ce travail fournit un cadre microscopique unifié reliant la flexibilité moléculaire, le couplage densité-ordre nématique et l'anisotropie élastique dans les nématiques semiflexibles. En dérivant les constantes élastiques de Frank et les longueurs de corrélation directement à partir des statistiques de chaînes au sein du formalisme SCFT, l'étude élucide les origines microscopiques de la grande anisotropie élastique observée dans les systèmes lyotropes. Les résultats valident l'utilisation des théories de continuum (comme Landau–de Gennes) pour décrire les structures de défauts, tout en soulignant la nécessité d'incorporer les fluctuations de densité et la flexibilité des chaînes pour prédire avec précision les réponses élastiques et les limites de stabilité. Les auteurs concluent que leur cadre offre une base pour de futures extensions aux chaînes auto-assemblées, polydisperses et à l'agrégation réversible, bien qu'ils ne proposent pas de nouveaux protocoles expérimentaux ou d'applications spécifiques au-delà de la portée de la présente analyse théorique.