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

Sintesi Tecnica: Teoria del campo auto-consistente per nematici semiflessibili

Problematica
Mentre la risposta lineare dei cristalli liquidi termotropici a piccole molecole è ben compresa, il comportamento dei cristalli liquidi liotropici — che comprendono polimeri, aggregati micellari e biofilamenti — presenta sfide teoriche e sperimentali significative. Una questione centrale è l'origine della grande anisotropia elastica osservata nei sistemi liotropici e l'accoppiamento tra la densità delle catene e l'ordine nematico. Le teorie microscopiche esistenti, ampiamente basate sul modello di Onsager, faticano a affrontare simultaneamente l'elasticità, i difetti topologici e l'accoppiamento densità-ordine in nematici semiflessibili spazialmente variabili. Inoltre, le teorie del continuo macroscopico (Oseen–Frank e Landau–de Gennes) introducono tipicamente coefficienti fenomenologici senza un esplicito collegamento con la statistica molecolare sottostante, lasciando inesplorate le origini microscopiche dell'anisotropia elastica e delle strutture del nucleo dei difetti.

Metodologia
Gli autori utilizzano un quadro di Teoria del Campo Auto-Consistente (SCFT) modellando i nematici semiflessibili come catene tipo wormlike (vermicellari). Il sistema è definito da un Hamiltoniano che incorpora:

1. Energia di flessione: governata dalla lunghezza di persistenza ($l_p$) e dalla lunghezza di contorno ($L$).
2. Interazioni Maier–Saupe (MS): interazioni orientative anisotrope ($u_2$) che guidano l'ordinamento nematico.
3. Volume escluso isotropo: un'interazione mediata dal solvente ($u_0$) che rappresenta le condizioni liotropiche (trattando implicitamente gli effetti del solvente tramite i parametri di Flory–Huggins).

Lo studio procede in tre fasi:

• Analisi dello Stato Uniforme: Gli autori derivano un funzionale di energia libera efficace dipendente dalla densità dei segmenti di catena ($\rho$) e dal parametro d'ordine unassiale ($S$). Ciò consente la determinazione dei confini di coesistenza isotropo-nematico (binodale) e dei limiti di stabilità (spinodale).
• Analisi della Risposta Lineare: Vengono introdotti deboli campi esterni per sondare la risposta del sistema a perturbazioni nella densità, nell'ordine unassiale e nell'orientamento del direttore. Ciò fornisce la matrice di suscettibilità inversa, da cui vengono estratti i costanti elastiche di Frank ($K_1, K_2, K_3$) e le lunghezze di correlazione.
• Simulazione della Struttura dei Difetti: Le equazioni SCFT complete sono risolte numericamente per disclinazioni +1/2 per analizzare i profili del direttore e le dimensioni del nucleo, confrontando i limiti delle catene rigide con i regimi semiflessibili.

Contributi Chiave e Risultati

• Accoppiamento Densità–Nematico e Stabilità:
  L'energia libera efficace derivata rivela un forte accoppiamento tra densità ($\rho$) e grado di ordine ($S$). Gli autori dimostrano che un forte accoppiamento densità-nematico può destabilizzare lo stato nematico. La posizione dell'instabilità spinodale dipende criticamente dal rapporto tra il volume escluso e la forza dell'interazione nematica ($u_0/u_2$). Al contrario, le distorsioni del direttore si accoppiano alla densità e alle variazioni di ordine solo ad ordini superiori, rimanendo efficacemente disaccoppiate nel regime di risposta lineare.

• Anisotropia Elastica e Flessibilità:
  Lo studio quantifica come la flessibilità della catena ($L/l_p$) alteri le costanti elastiche di Frank.

  • Catene Rigide: $K_3$ (bend/flessione) aumenta rapidamente con la forza di allineamento, mentre $K_1$ (splay/espansione) e $K_2$ (twist/torsione) rimangono comparabili.
  • Catene Semiflessibili: L'aumento della flessibilità sopprime fortemente l'elasticità di torsione ($K_2$) e di flessione ($K_3$), mentre l'elasticità di espansione ($K_1$) è influenzata in modo relativamente debole.
  • Crossover: Questa soppressione differenziale porta a un crossover da un'elasticità dominata dalla flessione a una dominata dall'espansione, evidenziato da un'inversione nel rapporto $K_1/K_3$ (da $<1$ a $>1$).

• Lunghezze di Correlazione e Anisotropia:
  L'analisi della risposta lineare rivela che, sebbene i coefficienti di gradiente diagonale della suscettibilità inversa suggeriscano lunghezze di correlazione più brevi parallelamente al direttore ($K_{\parallel} < K_{\perp}$), l'accoppiamento tra densità e ordine rinormalizza i coefficienti effettivi. Ciò risulta in un'anisotropia in cui la lunghezza di correlazione parallela al direttore è in realtà maggiore ($\xi_{\parallel} > \xi_{\perp}$). Questo risultato è in linea con le osservazioni sperimentali delle interfacce isotropo-nematico e delle strutture del nucleo dei difetti.

• Dimensioni del Nucleo dei Difetti:
  Le simulazioni numeriche di disclinazioni +1/2 mostrano che la dimensione del nucleo è anisotropa, con il raggio che è maggiore lungo la direzione in cui il gradiente di densità è parallelo al direttore. Lo studio stabilisce che la dimensione caratteristica del nucleo del difetto è proporzionale alla lunghezza di correlazione all'equilibrio ($\xi_Q$), coerentemente con la scalatura di Landau–de Gennes. Gli autori osservano che, sebbene i cromonici liotropici spesso presentino nuclei anomalamente grandi (potenzialmente a causa della scissione degli aggregati), il loro modello attuale, che esclude la polidispersità, cattura la relazione di scalatura tra la dimensione del nucleo e la lunghezza di correlazione guidata dalle fluttuazioni di densità.

• Comportamento Spinodale:
  Man mano che il sistema si avvicina allo spinodale nematico (controllato variando $u_0$), le lunghezze di correlazione divergono seguendo la scalatura $\xi \sim (u_0 - u_{0,i})^{-1/2}$. Tuttavia, nelle configurazioni di difetto, questa divergenza è preclusa dalla formazione locale di un dominio isotropo (tattoide) al nucleo, che sopprime l'ordine nematico prima che lo spinodale venga raggiunto.

Significato
Il lavoro fornisce un quadro microscopico unificato che lega la flessibilità molecolare, l'accoppiamento densità-ordine nematico e l'anisotropia elastica nei nematici semiflessibili. Derivando le costanti elastiche di Frank e le lunghezze di correlazione direttamente dalla statistica delle catene all'interno del formalismo SCFT, il lavoro elucida le origini microscopiche della grande anisotropia elastica osservata nei sistemi liotropici. I risultati validano l'uso di teorie del continuo (come Landau–de Gennes) per descrivere le strutture dei difetti, sottolineando però la necessità di incorporare le fluttuazioni di densità e la flessibilità delle catene per prevedere accuratamente le risposte elastiche e i limiti di stabilità. Gli autori concludono che il loro quadro offre una base per estensioni future a catene auto-assemblate, polidisperse e all'aggregazione reversibile, sebbene non propongano specifici nuovi protocolli sperimentali o applicazioni oltre l'ambito della presente analisi teorica.