Ordering in Confined Two-Dimensional Nematic Systems: Mesoscopic Simulations Based on Different Mean-Field Potentials

Titolo: Ordinamento in Sistemi Nematici Bidimensionali Confinati: Simulazioni Mesoscopiche Basate su Potenziali di Campo Medio Diversi

Autori: Humberto Híjara e Apala Majumdar

1. Il Problema

I cristalli liquidi nematici (NLC) sono materiali complessi che combinano fluidità e ordine orientazionale. La ricerca contemporanea mira a comprendere come le proprietà molecolari e la micromeccanica possano essere "coarse-grained" (mediata su scala maggiore) per produrre teorie macroscopiche analiticamente trattabili ed efficienti dal punto di vista computazionale.

Il problema centrale affrontato nel lavoro è la validità e i limiti delle teorie di campo continuo, in particolare la teoria di Landau-de Gennes (LdG), quando applicata a sistemi confinati su piccola scala o in prossimità di transizioni di fase. Mentre la teoria LdG è robusta per sistemi macroscopici ben ordinati, la sua accuratezza in domini nanometrici o vicino a difetti complessi rimane parzialmente aperta.
Gli autori si chiedono:

1. Le previsioni di modelli basati su particelle (N-MPCD) dipendono dalla scelta specifica del potenziale di campo medio sottostante?
2. Esiste una corrispondenza rigorosa tra i parametri dei modelli mesoscopici (N-MPCD) e i parametri continui (LdG)?
3. I modelli mesoscopici rivelano nuovi fenomeni non catturati dalle teorie continue?

2. Metodologia

Gli autori utilizzano la dinamica di collisione multi-particella per nematici (N-MPCD), un algoritmo stocastico basato su particelle che simula posizioni, velocità e orientazioni di "pacchetti" nematici, incorporando l'idrodinamica e l'accoppiamento bidirezionale tra ordine nematico e flusso del fluido.

Potenziali di Campo Medio:
Lo studio confronta tre diversi potenziali di campo medio per descrivere le interazioni tra le molecole:

1. Maier-Saupe (MS): Il modello classico per nematici uniaxiali.
2. Marrucci-Greco (MG): Una modifica del MS che include gradienti spaziali del parametro d'ordine, adatta per sistemi frustrati o con inhomogeneità spaziali.
3. Ilg-Karlin-Öttinger (IKÖ): Un potenziale più recente che garantisce la consistenza entropica e include una dipendenza non lineare della forza di interazione dal parametro d'ordine.

Setup delle Simulazioni:

• Geometria: Domini quadrati 2D con condizioni al contorno tangenziali (l'orientazione deve essere tangente ai bordi).
• Parametri: Variazione sistematica della forza di interazione $U$ (correlata alla temperatura inversa) e della dimensione del sistema $R$.
• Analisi Teorica: Derivazione di equazioni macroscopiche a partire dalla descrizione di Fokker-Planck (FPE) utilizzando chiusure tensoriali per i momenti di ordine superiore. Questo permette di stimare la lunghezza di coerenza nematica ($\xi_N$) e definire un parametro adimensionale $\lambda = l/\xi_N$ (dove $l$ è la dimensione del sistema), che funge da ponte tra la simulazione mesoscopica e la teoria continua LdG.

3. Contributi Chiave

• Derivazione Analitica dei Parametri Critici: Gli autori derivano analiticamente la forza di interazione critica ($U_c$) necessaria per l'inizio dell'ordinamento nematico per ciascuno dei tre potenziali, utilizzando chiusure di momento basate sulla distribuzione di probabilità delle orientazioni.
  • $U_c = 4$ per i potenziali MS e MG.
  • $U_c = 4\sqrt{2} \approx 5.66$ per il potenziale IKÖ.
• Corrispondenza Mesoscopico-Continuo: Viene stabilito un metodo rigoroso per mappare i parametri N-MPCD ($R, U$) ai parametri della teoria LdG continua, definendo la lunghezza di coerenza $\xi_N$ in funzione di $U$ e delle costanti elastiche.
• Validazione delle Chiusure: Viene proposta e validata una relazione di chiusura scalare per il parametro d'ordine del quarto ordine ($S_4$) in funzione del secondo ordine ($S$), essenziale per chiudere le equazioni di evoluzione derivate da FPE.

4. Risultati

• Convergenza verso la Teoria LdG: Per sistemi grandi ($R$ elevato) e ben ordinati ($U \gg U_c$, ovvero $\lambda \gg 1$), le previsioni N-MPCD sono universalmente consistenti tra i tre potenziali e coincidono con le previsioni della teoria ridotta di Landau-de Gennes (rLdG).
  • Le configurazioni di equilibrio stabili sono le soluzioni "Diagonale" e "Ruotata".
  • Le configurazioni metastabili e i percorsi di rilassamento mostrano la stessa topologia di difetti (saddle points di indice superiore) previsti dalla teoria continua.
• Differenze in Regimi di Piccola Scala:
  • Per domini piccoli ($R$ vicino a 8a) e forze di interazione vicine al critico ($U \approx U_c$), emergono differenze significative.
  • Il potenziale MG tende a favorire soluzioni di tipo "diagonale" anche in regimi dove il MS predice la soluzione "WORS" (Well Order Reconstruction Solution, con difetti lungo le diagonali).
  • Il potenziale IKÖ, a causa della sua non-linearità, produce difetti più estesi e regioni centrali disordinate più ampie rispetto a MS e MG. In particolare, per $U \approx U_c$, il IKÖ mantiene una struttura simile alla WORS su domini più grandi rispetto agli altri modelli.
• Universalità dei Percorsi di Rilassamento: I percorsi di rilassamento da configurazioni disordinate a stati di equilibrio rivelano strutture metastabili transitorie (come coppie di difetti $\pm 1/2$ o difetti interi) che sono universali e indipendenti dal potenziale di campo medio scelto, purché il sistema sia sufficientemente grande.

5. Significato e Implicazioni

• Robustezza della Teoria Continua: Il lavoro conferma che la teoria fenomenologica di Landau-de Gennes è estremamente robusta per descrivere sistemi nematici ben ordinati e macroscopici, indipendentemente dai dettagli microscopici o mesoscopici sottostanti.
• Definizione dei Limiti di Validità: Lo studio quantifica i limiti di validità delle teorie continue. Quando $R < R_c$ o $U < U_c$ (sistemi piccoli o debolmente ordinati), le previsioni mesoscopiche possono divergere significativamente dalle teorie continue, suggerendo che gli effetti di scala finita e le fluttuazioni stocastiche diventano dominanti.
• Prospettive Multiscala: I risultati forniscono una base fondamentale per lo sviluppo di metodi numerici multiscala. Ad esempio, si potrebbe utilizzare N-MPCD (con potenziali specifici) vicino a difetti, interfacce o temperature critiche, accoppiandolo a solver deterministici LdG nelle regioni ben ordinate, utilizzando le mappe derivate ($\lambda$, $U_c$) per definire i criteri di accoppiamento.
• Implicazioni per Applicazioni: Comprendere come i dettagli microscopici influenzino la struttura dei nuclei dei difetti in domini nanometrici è cruciale per nuove applicazioni in nanotecnologia, metamateriali e dispositivi di visualizzazione avanzati.

In sintesi, il paper dimostra che mentre la fisica macroscopica dei cristalli liquidi è ben catturata da modelli continui universali, la scelta del potenziale di campo medio diventa critica e determinante quando si esplorano regimi di confinamento forte o prossimità a transizioni di fase, offrendo uno strumento quantitativo per collegare la scala molecolare a quella macroscopica.