Super-Arrhenius temperature dependent viscosity due to liquid-liquid phase separation in the super-cooled Kob-Andersen model

Dit onderzoek toont aan dat super-Arrhenius temperatuurafhankelijke viscositeit in het supergekoelde Kob-Andersen-model kan worden verklaard door vloeibaar-vloeibaar fase-scheiding, waarbij een nieuw ordeparameter (WCN) en een Markov-netwerkmodel worden gebruikt om de relatie tussen fasegrenzen en glasovergang te kwantificeren.

De kern

Waarom wordt vloeistof zo stroopachtig? Een verhaal over vloeibare tweeling en glas

Stel je voor dat je een pot honing hebt. Als je hem warm houdt, stroomt hij makkelijk. Maar als je hem in de koelkast zet, wordt hij steeds stroperiger, totdat hij eruitziet als steen. Dit proces heet de glasovergang. Wetenschappers weten al lang dat dit gebeurt, maar ze hebben het moeilijk om precies te begrijpen waarom het zo plotseling gebeurt.

In dit artikel kijken twee onderzoekers (Jayme en Xueyu) naar een heel speciaal soort "vloeistof" in een computermodel. Ze ontdekken iets verrassends: voordat deze vloeistof tot glas verandert, splitst hij zich eigenlijk op in twee verschillende soorten vloeistoffen die door elkaar lopen.

Hier is hoe ze dit ontdekten en wat het betekent:

1. De "Zwaartepunt-Telefoon" (De WCN-methode)

Om te zien wat er in die vloeistof gebeurt, hebben de onderzoekers een slimme truc bedacht. Normaal gesproken kijken wetenschappers naar de dichtheid (hoe dicht de deeltjes op elkaar zitten) om vloeistoffen te onderscheiden. Maar in dit geval zien de twee vloeistoffen er qua dichtheid bijna hetzelfde uit. Het is alsof je twee identieke kledingstukken hebt, maar één is van zijde en de andere van wol. Je kunt ze niet zien, maar ze voelen heel anders.

Ze gebruikten een maatstaf die ze "Gewogen Coördinatiegetal" (WCN) noemen.

• De Analogie: Stel je voor dat elke deeltje in de vloeistof een telefoon heeft. Deze telefoon telt niet alleen hoeveel buren het heeft, maar ook hoe die buren eruitzien en hoe ze gerangschikt zijn.
• Door naar deze "telefonische" informatie te kijken, konden ze de deeltjes in twee groepen indelen: Groep A en Groep B. Zelfs als ze dicht op elkaar zitten, bleek dat Groep A en Groep B heel verschillende structuren hebben.

2. De Vloeibare Tweeling (Vloeistof-Vloeistof Scheiding)

Toen ze de temperatuur verlaagden, zagen ze iets fascinerends. De vloeistof splitste zich niet in een gas en een vloeistof (zoals water dat verdampt), maar in twee verschillende vloeibare staten.

• De Analogie: Denk aan een grote feestzaal waar iedereen dance. Plotseling beginnen twee groepen mensen zich te vormen: de "Slow-Dancers" (die rustig bewegen) en de "Fast-Dancers" (die wild dansen). Ze zitten door elkaar, maar ze bewegen niet met elkaar mee.
• In het computermodel zagen ze dat deze twee groepen (de twee vloeistoffen) een grens vormen. De ene groep is iets dichter dan de andere. Ze bewezen dat dit een echte scheiding is, net zoals olie en water, maar dan tussen twee vloeistoffen die er bijna hetzelfde uitzien.

3. De Grens die "Rijpt" (De sleutel tot glas)

Waarom wordt het glas dan zo hard? De onderzoekers ontdekten dat de grens tussen deze twee groepen (de "Fast" en "Slow" dancers) niet stil staat.

• De Analogie: Stel je voor dat de "Slow-Dancers" een grote, zachte bol vormen in het midden van de zaal, en de "Fast-Dancers" eromheen zitten. Na verloop van tijd proberen de bollen groter te worden om energie te besparen (net zoals een waterdruppel rond wordt). Dit proces heet vergroting of coarsening.
• Maar omdat de temperatuur zo laag is, is het heel moeilijk voor de deeltjes om van de ene groep naar de andere te springen om die bollen groter te maken. Het is alsof de Fast-Dancers vastzitten in stroop. Ze willen wel bewegen, maar de "muur" tussen de groepen is te zwaar om te doorbreken.

4. De Viscositeit (De Stroop-meting)

De onderzoekers berekenden hoe stroperig (viskeus) de vloeistof werd naarmate het kouder werd.

• Normaal gesproken wordt iets stroperiger naarmate het kouder wordt, maar dit gaat vaak lineair (zoals een rechte lijn).
• In dit model zagen ze iets extreems: de stroperigheid nam exponentieel toe. Het werd niet alleen harder, het werd ontzettend snel onbeweeglijk.
• De conclusie: Dit extreme "vastlopen" (de glasovergang) wordt veroorzaakt door de moeite die het kost om de grenzen tussen die twee vloeibare groepen te laten groeien. De vloeistof probeert te "rijpen", maar de temperatuur is te laag om de grenzen snel genoeg te laten bewegen. Het resultaat? De hele vloeistof bevriest in plaats van te stromen.

Samenvatting in één zin

Deze studie toont aan dat glas niet zomaar "bevroren water" is, maar eerder een vloeistof die in paniek raakt omdat hij probeert zich te splitsen in twee verschillende soorten vloeistoffen, maar te koud is om die splitsing te voltooien, waardoor hij vastloopt in een glasachtige staat.

Waarom is dit belangrijk?
Als we begrijpen dat glasovergang eigenlijk een onvoltooide scheiding van vloeistoffen is, kunnen we misschien nieuwe materialen maken die niet zo snel breken, of beter begrijpen hoe water in extreme kou zich gedraagt (wat belangrijk is voor het begrijpen van het weer op andere planeten of in de ruimte).

Technische samenvatting

Titel: Super-Arrhenius temperatuurafhankelijke viscositeit veroorzaakt door vloeibaar-vloeibaar faseafscheiding in het afgekoelde Kob-Andersen-model

1. Het Probleem

Het fundamentele mechanisme achter de glasovergang (glass transition) blijft een van de grootste uitdagingen in de statistische mechanica en materiaalkunde. Een kenmerkend aspect van glasvormende vloeistoffen is de dramatische toename van de viscositeit bij afkoeling in het onderkoelde gebied, vaak beschreven als "super-Arrhenius" gedrag (een exponentiële groei die sterker is dan de lineaire Arrhenius-relatie).
Hoewel er veel onderzoek is gedaan naar dynamische heterogeniteit en vloeibaar-vloeibaar faseovergangen (LLT), is de relatie tussen thermodynamisch gedreven vloeibaar-vloeibaar faseafscheiding en de dynamische vertraging die leidt tot glasvorming nog niet volledig opgehelderd. Bestaande methoden om viscositeit te berekenen (zoals de Green-Kubo-methode) falen vaak bij zeer lage temperaturen omdat de relaxatietijden te lang worden voor directe simulaties.

2. Methodologie

De auteurs gebruiken een combinatie van moleculaire dynamica (MD) simulaties en geavanceerde data-analyse om dit probleem aan te pakken:

• Systeem: Er wordt gebruikgemaakt van het Kob-Andersen (KA) binair Lennard-Jones-model (80% type A, 20% type B deeltjes). Dit model staat bekend om zijn vermogen om kristallisatie te voorkomen en is ideaal voor het bestuderen van glasvorming.
• Ordeparameter (WCN): In plaats van traditionele parameters zoals dichtheid, gebruiken de auteurs de Gewogen Coördinatiegetallen (Weighted Coordination Numbers - WCN).
  • WCN's worden afgeleid van de radiale verdelingsfunctie $g(r)$ door Gaussische verdelingen op de kenmerkende pieken te plaatsen.
  • Dit creëert een vector voor elk deeltje dat zijn lokale omgeving beschrijft.
  • Principal Component Analysis (PCA) wordt toegepast om de dimensionaliteit te reduceren, gevolgd door K-means clustering om de deeltjes in twee thermodynamisch verschillende vloeibare toestanden te classificeren.
• Simulatieprotocol: Systemen worden voorbereid bij hoge temperatuur ($T^*=1.2$) en vervolgens instantaan afgekoeld (quenched) naar lagere temperaturen ($T^* = 0.3$ tot $0.58$).
• Viscositeitsberekening: Voor het diep onderkoelde gebied waar directe simulatie onmogelijk is, wordt een Markov Netwerk Model (MNM) gebruikt. Dit model benut de informatie over de vloeibaar-vloeibaar faseafscheiding om de viscositeit te schatten op basis van de overdrachtssnelheden tussen de twee geclusterde domeinen.

3. Belangrijkste Bijdragen en Resultaten

A. Constructie van het Fase-diagram

• De auteurs hebben succesvol de gas-vloeistof binodale lijn gereconstrueerd en deze methode uitgebreid naar de vloeistof-vloeistof binodale lijn.
• Ze hebben het kritieke punt van de vloeistof-vloeistof overgang bepaald: $T^*_c \approx 0.6$ en $\rho^*_c \approx 1.1$.
• De hefboomregel (lever rule) werd verifieerd voor verschillende bulk-dichtheden, wat sterke bewijzen levert voor lokaal thermodynamisch evenwicht tussen de twee vloeibare fasen.

B. Druk- en Dichtheidsprofielen

• Door druk- en dichtheidsprofielen te analyseren over het interface tussen de twee vloeibare fasen, tonen de auteurs mechanisch evenwicht aan.
• Ze onderscheiden tussen convexe (positieve kromming) en concave (negatieve kromming) interfaces.
• De concave interfaces vertonen een lichte schending van het mechanisch evenwicht (een toename in dichtheid en druk net voor de spanningsoppervlakte), wat fungeert als drijvende kracht voor de vergrotingskinetiek (coarsening kinetics) van de faseafscheiding.
• Oppervlaktespanning ($\sigma_s$) werd berekend via integratie van drukprofielen en via de Laplace-vergelijking, waarbij beide methoden vergelijkbare resultaten opleverden.

C. Viscositeit en Super-Arrhenius Gedrag

• De kern van het artikel is het modelleren van de viscositeit ($\eta$) met het Markov Netwerk Model. De auteurs stellen dat de overgang tussen de twee vloeibare domeinen wordt geactiveerd door de interfaciële vrije energie (gerelateerd aan de oppervlaktespanning).
• De formule voor viscositeit wordt afgeleid als:
  $$\eta(T) \propto \frac{1}{T} e^{\phi / k_B T}$$
  waarbij $\phi$ de som is van de lokale oppervlaktespanning vermenigvuldigd met het oppervlak (de totale interfaciële vrije energie).
• Resultaat: De berekende viscositeit toont een toename van bijna 16 orde van grootte in het temperatuurbereik van $T^*=0.55$ tot $0.36$.
• Dit gedrag is duidelijk super-Arrhenius, wat overeenkomt met literatuurresultaten voor het KA-model, maar nu verklaard via een specifiek thermodynamisch mechanisme.

4. Betekenis en Conclusie

De studie biedt een nieuw en overtuigend perspectief op de glasovergang:

1. Mechanisme: De dramatische vertraging van de dynamiek (glasvorming) wordt niet alleen toegeschreven aan entropische barrières of lokale structuurveranderingen, maar direct gekoppeld aan de kinetiek van de vloeibaar-vloeibaar faseafscheiding. De "vergroting" (coarsening) van de domeinen wordt vertraagd door de hoge oppervlaktespanning en de complexe morfologie van het interface.
2. Validatie: Het feit dat het Markov Netwerk Model, gebaseerd op slechts twee toestanden (de twee vloeibare fasen), de complexe super-Arrhenius viscositeit nauwkeurig kan reproduceren, suggereert dat de vloeibaar-vloeibaar faseafscheiding een fundamenteel mechanisme is voor glasvorming in dit systeem.
3. Toekomstperspectief: De auteurs wijzen erop dat dit mechanisme mogelijk ook van toepassing is op andere systemen, zoals het ST2-watermodel, wat zou kunnen leiden tot een universeel begrip van glasovergangen.

Samenvattend bewijst dit werk dat vloeibaar-vloeibaar faseafscheiding niet alleen een thermodynamisch fenomeen is, maar de drijvende kracht kan zijn achter de extreme viscositeitsstijging die kenmerkend is voor de overgang van vloeistof naar glas.