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

Technische Samenvatting: Zelfconsistente veldentheorie van semiflexibele nematische systemen

Probleemstelling
Hoewel de lineaire respons van kleine-molecuul thermotrope kristallen goed begrepen is, vormen het gedrag van lyotrope kristallen—bestaande uit polymeren, micellaire aggregaten en biofilamenten—aanzienlijke theoretische en experimentele uitdagingen. Een centraal probleem is de oorsprong van de grote elastische anisotropie die wordt waargenomen in lyotrope systemen en de koppeling tussen ketendichtheid en nematische orde. Bestaande microscopische theorieën, grotendeels gebaseerd op het Onsager-model, worstelen om tegelijkertijd elasticiteit, topologische defecten en dichtheid–orde-koppeling te adresseren in ruimtelijk variërende semiflexibele nematische systemen. Bovendien introduceren macroscopische continuümtheorieën (Oseen–Frank en Landau–de Gennes) doorgaans fenomenologische coëfficiënten zonder expliciete verbinding met onderliggende moleculaire statistieken, waardoor de microscopische oorsprong van elastische anisotropie en defectkernstructuren onverkend blijft.

Methodologie
De auteurs maken gebruik van een Self-Consistent Field Theory (SCFT) raamwerk waarbij semiflexibele nematische systemen worden gemodelleerd als wormachtige ketens. Het systeem wordt gedefinieerd door een Hamiltoniaan die de volgende elementen bevat:

1. Buigenergie: Gereguleerd door de persistentielengte ($l_p$) en de contourlengte ($L$).
2. Maier–Saupe (MS) interacties: Anisotrope oriëntatie-interacties ($u_2$) die drijven op de vorming van nematische orde.
3. Isotrope uitgesloten volume-interactie: Een door oplosmiddel gemedieerde interactie ($u_0$) die lyotrope condities representeert (waarbij effecten van het oplosmiddel impliciet worden behandeld via Flory–Huggins parameters).

De studie verloopt in drie fasen:

• Analyse van de uniforme toestand: De auteurs leiden een effectieve vrije energiefunctionaal af die afhankelijk is van de ketensegmentdichtheid ($\rho$) en de uniaxiale ordeparameter ($S$). Dit maakt de bepaling van de isotroop-nematische coëxistentie (binodaal) en de stabiliteitsgrenzen (spinodaal) mogelijk.
• Lineaire respons analyse: Zwakke externe velden worden geïntroduceerd om de respons van het systeem op perturbaties in dichtheid, uniaxiale orde en director-oriëntatie te onderzoeken. Dit levert de inverse susceptibiliteitsmatrix op, waaruit de Frank elasticiteitsconstanten ($K_1, K_2, K_3$) en correlatielengten worden geëxtraheerd.
• Simulatie van defectstructuren: De volledige SCFT-vergelijkingen worden numeriek opgelost voor +1/2 disclinaties om directorprofielen en kernafmetingen te analyseren, waarbij de regime van rigide ketens wordt vergeleken met semiflexibele regimes.

Belangrijkste bijdragen en resultaten

• Dichtheid–Nematische Koppeling en Stabiliteit:
  De afgeleide effectieve vrije energie onthult een sterke koppeling tussen dichtheid ($\rho$) en de mate van orde ($S$). De auteurs demonstreren dat sterke dichtheid–nematische koppeling de nematische toestand kan destabiliseren. De locatie van de spinodale instabiliteit hangt kritisch af van de ratio tussen het uitgesloten volume en de nematische interactiekracht ($u_0/u_2$). In tegenstelling hiertoe koppelen director-distorties aan dichtheids- en ordevariaties pas bij hogere orden en blijven ze in het lineaire responsregime effectief ontkoppeld.

• Elastische Anisotropie en Flexibiliteit:
  De studie kwantificeert hoe de flexibiliteit van de keten ($L/l_p$) de Frank elasticiteitsconstanten verandert.

  • Rigide Ketens: $K_3$ (buiging) neemt snel toe met de uitlijningssterkte, terwijl $K_1$ (splay) en $K_2$ (twist) vergelijkbaar blijven.
  • Semiflexibele Ketens: Toenemende flexibiliteit onderdrukt de twist- ($K_2$) en buigelasticiteit ($K_3$) sterk, terwijl de splay-elasticiteit ($K_1$) relatief zwak wordt beïnvloed.
  • Crossover: Deze differentiële onderdrukking leidt tot een crossover van buig-gedomineerde naar splay-gedomineerde elasticiteit naarmate de flexibiliteit toeneemt, wat wordt aangetoond door een omkering in de ratio $K_1/K_3$ (van $<1$ naar $>1$).

• Correlatielengtes en Anisotropie:
  De lineaire responsanalyse laat zien dat hoewel de diagonale gradiëntcoëfficiënten van de inverse susceptibiliteit suggereren dat er kortere correlatielengtes parallel aan de director zijn ($K_{\parallel} < K_{\perp}$), de koppeling tussen dichtheid en orde de effectieve coëfficiënten hernormaliseert. Dit resulteert in een anisotropie waarbij de correlatielengte parallel aan de director daadwerkelijk groter is ($\xi_{\parallel} > \xi_{\perp}$). Deze bevinding komt overeen met experimentele observaties van isotroop-nematische interfaces en defectkernstructuren.

• Defectkernafmetingen:
  Numerieke simulaties van +1/2 disclinaties tonen aan dat de kernafmeting anisotroop is, waarbij de straal groter is in de richting waar de dichtheidsgradiënt parallel is aan de director. De studie stelt vast dat de karakteristieke defectkernafmeting proportioneel is aan de evenwichtscorrelatielengte ($\xi_Q$), wat consistent is met de Landau–de Gennes schaling. De auteurs merken op dat hoewel lyotrope chromonica's vaak anomalistisch grote kernen vertonen (mogelijk door aggregaat-scissie), hun huidige model, dat polydispersiteit uitsluit, de schaleringsrelatie tussen kernafmeting en correlatielengte bepaald door dichtheidsfluctuaties vastlegt.

• Spinodaal Gedrag:
  Wanneer het systeem de nematische spinodaal nadert (gecontroleerd door het variëren van $u_0$), divergeren de correlatielengtes volgens de schaling $\xi \sim (u_0 - u_{0,i})^{-1/2}$. Echter, in defectconfiguraties wordt deze divergentie voortijdig gestopt door de lokale vorming van een isotrope domein (tactoïde) in de kern, die de nematische orde onderdrukt voordat de spinodaal wordt bereikt.

Significantie
Dit artikel biedt een verenigd microscopisch raamwerk dat de link legt tussen moleculaire flexibiliteit, dichtheid–nematische orde-koppeling en elastische anisotropie in semiflexibele nematische systemen. Door de Frank elasticiteitsconstanten en correlatielengten direct af te leiden uit ketenstatistieken binnen het SCFT-formalisme, verheldert dit werk de microscopische oorsprong van de grote elastische anisotropie die in lyotrope systemen wordt waargenomen. De resultaten valideren het gebruik van continuümtheorieën (zoals Landau–de Gennes) voor het beschrijven van defectstructuren, terwijl ze ook benadrukken dat het integreren van dichtheidsfluctuaties en ketenflexibiliteit noodzakelijk is om de elastische respons en stabiliteitsgrenzen accuraat te voorspellen. De auteurs concluderen dat hun raamwerk een basis biedt voor toekomstige uitbreidingen naar zelfassemblerende, polydisperse ketens en reversibele aggregatie, hoewel zij geen specifieke nieuwe experimentele protocollen of toepassingen voorstellen buiten de reikwijdte van de huidige theoretische analyse.