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

Technische Zusammenfassung: Selbstkonsistente Feldtheorie für semiflexible Nematika

Problemstellung
Während die lineare Antwort von niedermolekularen thermotropen Flüssigkristallen gut verstanden ist, stellen das Verhalten von lyotropen Flüssigkristallen – bestehend aus Polymeren, mizellaren Aggregaten und Biofilamenten – erhebliche theoretische und experimentelle Herausforderungen dar. Ein zentrales Problem ist der Ursprung der beobachteten großen elastischen Anisotropie in lyotropen Systemen sowie die Kopplung zwischen Kettendichte und nematischer Ordnung. Bestehende mikroskopische Theorien, die weitgehend auf dem Onsager-Modell basieren, haben Schwierigkeiten, gleichzeitig Elastizität, topologische Defekte und die Dichte–Ordnungs-Kopplung in räumlich variierenden semiflexiblen Nematika zu adressieren. Darüber hinaus führen makroskopische Kontinuumstheorien (Oseen–Frank und Landau–de Gennes) typischerweise phänomenologische Koeffizienten ein, ohne eine explizite Verbindung zu den zugrunde liegenden molekularen Statistiken herzustellen, wodurch die mikroskopischen Ursprünge der elastischen Anisotropie und der Defektkernstrukturen ungeklärt bleiben.

Methodik
Die Autoren verwenden einen Rahmen der Selbstkonsistenten Feldtheorie (SCFT), der semiflexible Nematika als wurmähnliche Ketten (wormlike chains) modelliert. Das System wird durch einen Hamiltonian definiert, der Folgendes umfasst:

1. Biegeenergie: Gesteuert durch die Persistenzlänge ($l_p$) und die Konturlänge ($L$).
2. Maier–Saupe (MS)-Wechselwirkungen: Anisotrope orientierende Wechselwirkungen ($u_2$), die eine nematische Ordnung vorantreiben.
3. Isolisches ausgeschlossenes Volumen: Eine lösungsmittelvermittelte Wechselwirkung ($u_0$), die lyotrope Bedingungen repräsentiert (wobei Lösungsmittel-Effekte implizit über Flory–Huggins-Parameter behandelt werden).

Die Studie erfolgt in drei Stufen:

• Analyse des Gleichgewichtszustands: Die Autoren leiten ein effektives Freie-Energie-Funktionalional ab, das von der Kettensegmentdichte ($\rho$) und dem uniaxialen Ordnungsparameter ($S$) abhängt. Dies ermöglicht die Bestimmung der Isotrop–Nematisch-Koexistenz (Binodale) und der Stabilitätsgrenzen (Spinodale).
• Lineare Antwortanalyse: Schwache externe Felder werden eingeführt, um die Reaktion des Systems auf Störungen in Dichte, uniaxialer Ordnung und Direktororientierung zu untersuchen. Dies liefert die inverse Suszeptibilitätsmatrix, aus der die Frank-Elastizitätskonstanten ($K_1, K_2, K_3$) und Korrelationslängen extrahiert werden können.
• Simulation der Defektstruktur: Die vollständigen SCFT-Gleichungen werden numerisch für +1/2-Disklinationen gelöst, um Direktorprofile und Kerngrößen zu analysieren, wobei die Grenzfälle starrer Ketten mit semiflexiblen Regimen verglichen werden.

Wesentliche Beiträge und Ergebnisse

• Dichte–Nematische Kopplung und Stabilität:
  Die abgeleitete effektive Freie Energie offenbart eine starke Kopplung zwischen der Dichte ($\rho$) und dem Grad der Ordnung ($S$). Die Autoren zeigen, dass eine starke Dichte–Nematische Kopplung den nematischen Zustand destabilisieren kann. Die Lage der Spinodal-Instabilität hängt entscheidend vom Verhältnis von ausgeschlossenem Volumen zur nematischen Wechselwirkungsstärke ($u_0/u_2$) ab. Im Gegensatz dazu koppeln Direktorverzerrungen nur bei höheren Ordnungen an Dichte- und Ordnungsvariationen und bleiben im Bereich der linearen Antwort effektiv entkoppelt.

• Elastische Anisotropie und Flexibilität:
  Die Studie quantifiziert, wie die Kettenflexibilität ($L/l_p$) die Frank-Elastizitätskonstanten verändert.

  • Starre Ketten: $K_3$ (Bend/Beugung) steigt mit der Ausrichtungsstärke schnell an, während $K_1$ (Splay/Aufweitung) und $K_2$ (Twist/Drehung) vergleichbar bleiben.
  • Semiflexible Ketten: Eine zunehmende Flexibilität unterdrückt die Twist- ($K_2$) und Bend-Elastizität ($K_3$) stark, während die Splay-Elastizität ($K_1$) vergleichsweise schwach beeinflusst wird.
  • Crossover: Diese differenzielle Unterdrückung führt zu einem Übergang von einer Bend-dominierten zu einer Splay-dominierten Elastizität mit zunehmender Flexibilität, was durch eine Umkehrung des Verhältnisses $K_1/K_3$ (von $<1$ zu $>1$) belegt wird.

• Korrelationslängen und Anisotropie:
  Die lineare Antwortanalyse zeigt, dass, obwohl die diagonalen Gradientenkoeffizienten der inversen Suszeptibilität kürzere Korrelationslängen parallel zum Direktor suggerieren ($K_{\parallel} < K_{\perp}$), die Kopplung zwischen Dichte und Ordnung die effektiven Koeffizienten renormiert. Dies führt zu einer Anisotropie, bei der die Korrelationslänge parallel zum Direktor tatsächlich größer ist ($\xi_{\parallel} > \xi_{\perp}$). Dieser Befund steht im Einklang mit experimentellen Beobachtungen von Isotrop–Nematisch-Grenzflächen und Defektkernstrukturen.

• Defektkerngrößen:
  Numerische Simulationen von +1/2-Disklinationen zeigen, dass die Kerngröße anisotrop ist, wobei der Radius entlang der Richtung größer ist, in der der Dichtegradient parallel zum Direktor verläuft. Die Studie stellt fest, dass die charakteristische Defektkerngröße proportional zur äquilibrierten Korrelationslänge ($\xi_Q$) ist, was konsistent mit der Landau–de Gennes-Skalierung ist. Die Autoren merken an, dass lyotrope Chromonika oft anomal große Kerne aufweisen (potenziell aufgrund von Aggregat-Spaltung), ihr aktuelles Modell jedoch, das Polydispersität ausschließt, die Skalierungsbeziehung zwischen Kerngröße und Korrelationslänge unter dem Einfluss von Dichtefluktuationen erfasst.

• Spinodales Verhalten:
  Wenn sich das System dem nematischen Spinodale nähert (gesteuert durch Variation von $u_0$), divergieren die Korrelationslängen nach der Skalierung $\xi \sim (u_0 - u_{0,i})^{-1/2}$. In Defektkonfigurationen wird diese Divergenz jedoch durch die lokale Bildung eines isotropen Bereichs (Taktoid) im Kern verhindert, welcher die nematische Ordnung unterdrückt, bevor das Spinodale erreicht wird.

Bedeutung
Das Paper liefert einen vereinheitlichten mikroskopischen Rahmen, der Kettenflexibilität, Dichte–Nematische Ordnung-Kopplung und elastische Anisotropie in semiflexiblen Nematika verknüpft. Durch die Ableitung der Frank-Elastizitätskonstanten und Korrelationslängen direkt aus der Kettenstatistik innerhalb des SCFT-Formalismus erläutert die Arbeit die mikroskopischen Ursprünge der großen elastischen Anisotropie, die in lyotropen Systemen beobachtet wird. Die Ergebnisse validieren die Verwendung von Kontinuumstheorien (wie Landau–de Gennes) zur Beschreibung von Defektstrukturen, betonen jedoch die Notwendigkeit, Dichtefluktuationen und Kettenflexibilität einzubeziehen, um elastische Antworten und Stabilitätsgrenzen präzise vorherzusagen. Die Autoren kommen zu dem Schluss, dass ihr Rahmenwerk eine Basis für zukünftige Erweiterungen auf selbstorganisierte, polydisperse Ketten und reversible Aggregation bietet, sie jedoch keine spezifischen neuen experimentellen Protokolle oder Anwendungen über den Umfang der aktuellen theoretischen Analyse hinaus vorschlagen.