Metastability and ripening of multi-component liquid mixtures

Technische Samenvatting: Metastabiliteit en Rijping van Meercomponenten Vloeistofmengsels

Probleemstelling
Fasescheiding is een fundamenteel mechanisme in de zachte materie, cruciaal voor biofysische organisatie (bijv. membraanloze organellen), het ontstaan van leven en voedselwetenschap. Hoewel de thermodynamica van meercomponentenmengsels, met name met betrekking tot spinodale instabiliteiten, uitgebreid is bestudeerd, blijven hun uit-evenwichtsdynamica onderbelicht. Specifiek de laatstadium-coarsening (rijping) en metastabiliteit van mengsels met een groot aantal componenten ($M$) zijn niet goed begrepen. De heersende intuïtie suggereert dat meercomponentensystemen de klassieke Ostwald-rijpingswetten erven die worden waargenomen in binaire mengsels (druppelstraal $R \sim t^{1/3}$ en druppelgetal $n \sim t^{-2/3}$ in 2D). Echter, numerieke studies die dit bevestigen of weerleggen in aanwezigheid van gedesordreerde interacties zijn schaars.

Methodologie
De auteurs maken gebruik van een combinatie van analytische veldentheorie en grootschalige numerieke simulaties om incompressibele meercomponenten vloeistofmengsels met willekeurige interacties te bestuderen.

• Theoretisch Kader: Het systeem wordt beschreven in het $T$-$V$-$N_i$-ensemble met behulp van een veldentheorie voor de spatio-temporele evolutie van de volumefracties van de opgeloste stoffen $\phi_i$. De vrije energiedichtheid maakt gebruik van het Flory-Huggins-model, waarbij de paarinteracties $\chi_{ij}$ worden gesampled uit een Gaussische verdeling met een gemiddelde $b$ en een standaarddeviatie die schaalt als $\sigma/\sqrt{M}$. Deze schaling zorgt ervoor dat thermodynamische grootheden onafhankelijk blijven van het aantal opgeloste stoffen wanneer $M \to \infty$.
• Analytische Benadering: De auteurs passen formalismen toe uit de zachte spin-glazen en de willekeurige matrixtheorie (random matrix theory). Ze voeren een lineaire stabiliteitsanalyse uit op de Hessiaan van de vrije energie om de spinodale curves te bepalen en analyseren de resolvent van de relaxatieoperator. Dit stelt hen in staat om de voorwaarden voor condensatie versus demixingsinstabiliteiten af te leiden.
• Numerieke Simulaties: Grootschalige simulaties van de gekoppelde veldvergelijkingen (Eq. 2) worden uitgevoerd. De auteurs initialiseren systemen met homogene concentraties en introduceren puntvormige perturbaties om druppelvorming te initiëren. Ze volgen de evolutie van de totale volumefracties en individuele concentraties van de opgeloste stoffen om onderscheidende fasen te identificeren en de oppervlakte en het aantal druppels over de tijd te meten.

Belangrijkste Bijdragen en Resultaten

1. Drie Dynamische Regimes: Het artikel schetst drie verschillende dynamische stadia van fasescheiding:

   • Vroege Fase (Spinodale Decompositie): De auteurs leiden de spinodale curve analytisch af. Ze tonen aan dat de stabiliteit van de homogene fase wordt beheerst door het eigenwaardespectrum van een willekeurige relaxatieoperator. Voor grote totale volumefracties wordt de spinodale temperatuur bepaald door de voorwaarde dat de onderste rand van het eigenwaardespectrum nul raakt, wat leidt tot kritische vertraging (critical slowing down). Het verval van perturbaties volgt een machtswet ($t^{-1}$) die karakteristiek is voor Wigner-semicirkel spectrale dichtheden.
   • Intermediaire Fase (Aanvang van Multifase Coexistentie): Na de initiële demixing vormt het systeem twee lokaal stabiele fasen (oplosmiddel-rijke en opgeloste stof-rijke). De auteurs demonstreren dat binnen deze fasen de eigenwaardes van de Hessiaan-matrix afwijken van de eenvoudige superpositie van een Wigner-semicirkel en een enkele eigenwaarde. Dit leidt tot een complexe distributie van relaxatietijdsschalen, waarbij sommige modi snel equidibreren terwijl andere traag blijven, wat het systeem potentieel in metastabiele configuraties kan vangen.
   • Laatste Fase (Rijping): De studie onderzoekt de coarsening van druppels. In tegenstelling tot de verwachting dat meercomponentenmengsels de binaire schaling erven, onthullen de simulaties significante afwijkingen.

2. Vertraagde Ostwald-rijping:

   • In binaire mengsels ($M=1$) volgt het systeem de klassieke Ostwald-schaling ($\langle A \rangle \sim t^{2/3}$, $n \sim t^{-2/3}$).
   • In meercomponentenmengsels ($M=5, 10, 15$) is de initiële rijping ernstig vertraagd. De schalingswetten wijken significant af van de $t^{1/3}$ en $t^{-2/3}$ exponenten.
   • De auteurs schrijven deze vertraging toe aan glasachtige relaxatie en het ontstaan van langdurige metastabiele toestanden. De aanwezigheid van meerdere onderscheidende fasen (tot $M+2$) creëert ruimtelijke heterogeniteiten. De cruciale aanname van Ostwald-rijping — dat druppels geïsoleerd zijn en een uniforme verre-veld concentratie ervaren — wordt in deze systemen geschonden.
   • Hoewel de binaire schaling uiteindelijk wordt teruggevonden na een lange transient (zodra een significant deel van de initiële druppels is opgelost), worden de intermediaire dynamieken gedomineerd door niet-klassiek gedrag.

3. Metastabiliteit en Initiële Condities:

   • Een kenmerk van glasachtige dynamica is de afhankelijkheid van de eindtoestand van de initiële condities. De auteurs demonstreren dat identieke controleparameters (temperatuur, interactiematrix) kunnen leiden tot verschillende langdurige metastabiele toestanden, afhankelijk van het initiële aantal druppel-zaden.
   • Starten met een hoog aantal druppels ($N=500$) leidt tot een toestand waarin onderscheidende fasen ruimtelijk gescheiden blijven.
   • Starten met minder druppels ($N=200$) leidt tot een toestand waarin fasen groeien, instabiel worden en nieuwe fasen aan hun interfaces nucleëren, wat resulteert in bevochtiging tussen de fasen.
   • Deze gevoeligheid suggereert dat het systeem gevangen is in lokale minima van de vrije energielandschap, gekenmerkt door verschillende bevochtigingshoeken.

Betekenis en Claims
Het artikel claimt de eerste uitgebreide karakterisering te bieden van de fasescheidingsdynamica in meercomponentenmengsels met gedesordreerde interacties, waarbij de kloof tussen spinodale analyse en laatstadium-coarsening wordt overbrugd.

• Afwijking van Binaire Intuïtie: Het werk daagt de heersende intuïtie uit dat meercomponentensystemen simpelweg binaire schalingswetten erven. Het stelt vast dat de complexiteit van interacties en het potentieel voor meerdere coëxisterende fasen de rijping-kinetica fundamenteel veranderen.
• Glasachtige Dynamica in Fasescheiding: Door de vertraagde rijping en de afhankelijkheid van initiële condities te koppelen aan glasachtige relaxatie en metastabiliteit, suggereert het artikel dat fasescheiding in complexe biologische mengsels (zoals het cytoplasma) mogelijk wordt beheerst door niet-evenwichtsdynamica die typisch is voor spin-glazen, in plaats van eenvoudige thermodynamische minimalisatie.
• Theoretische Instrumenten: De auteurs stellen voor dat theoretische instrumenten ontwikkeld voor spin-glazen, zoals dynamische gemiddelde-veldtheorieën en mode-koppelingstheorieën, noodzakelijk zijn om de rijke dynamische regimes van deze systemen volledig te verklaren, aangezien eenvoudige lineaire stabiliteitsanalyse onvoldoende is om het laatstadium-gedrag te vatten.

De auteurs blijven bescheiden over de uiteindelijke stationaire toestand en merken op dat, hoewel hun simulaties langdurige metastabiele toestanden laten zien, de ware evenwichtstoestand (waarschijnlijk met niet-nul bevochtigingshoeken) mogelijk pas bereikt wordt op tijdschalen die buiten de huidige numerieke capaciteiten liggen.