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

Cette étude utilise le nombre de coordination pondéré et un modèle de réseau de Markov pour démontrer que la séparation de phase liquide-liquide dans le modèle de Kob-Andersen sur-refroidi, caractérisée par un coarsening dépendant de la température, constitue un mécanisme plausible expliquant la transition vitreuse et la viscosité de type super-Arrhenius.

L'essentiel

Le Secret du Verre : Quand deux liquides se disputent la place

Imaginez que vous avez un verre d'eau très chaud. Si vous le laissez refroidir doucement, il gèle et devient de la glace (un solide cristallin, ordonné). Mais si vous refroidissez certains liquides très vite, ils ne gèlent pas. Ils deviennent de plus en plus visqueux, comme du miel, puis comme du goudron, jusqu'à devenir un verre : solide mais sans structure ordonnée.

La question que les scientifiques se posent depuis des décennies est : pourquoi ces liquides deviennent-ils si "lents" et si durs juste avant de devenir du verre ?

Dans cet article, Jayme Brickley et Xueyu Song proposent une réponse fascinante : ce n'est pas un simple ralentissement, c'est une séparation.

1. La découverte : Deux liquides dans un seul verre

Habituellement, on pense qu'un liquide est un liquide. Mais les chercheurs ont découvert que, dans certaines conditions de froid extrême (sans être encore solide), le liquide Kob-Andersen (un modèle informatique de liquide) se comporte comme s'il contenait deux liquides différents qui ne veulent pas se mélanger.

C'est un peu comme si vous versiez de l'huile et de l'eau dans un même verre. Au début, c'est un mélange trouble, mais rapidement, l'huile forme des gouttes et l'eau forme le reste. Ici, les chercheurs ont vu que le liquide se séparait en deux "états" distincts, comme deux équipes qui se disputent le terrain.

Pour voir cela, ils ont utilisé un outil spécial appelé le Nombre de Coordination Pondéré (WCN).

• L'analogie : Imaginez que chaque molécule a une "carte d'identité" basée sur ses voisins immédiats. La plupart des méthodes de mesure sont comme une photo floue qui ne voit que la foule globale. Le WCN, lui, est une loupe puissante qui permet de distinguer les "gros" voisins des "petits" voisins. Grâce à cette loupe, ils ont pu voir que certaines molécules formaient des groupes compacts (l'équipe A) et d'autres des groupes plus lâches (l'équipe B).

2. La frontière qui bouge : Le jeu de la coalescence

Une fois que ces deux liquides se sont séparés, ils forment des îles (des domaines) qui flottent dans le liquide. Entre ces îles, il y a une frontière.

• L'analogie : Imaginez des bulles de savon dans une baignoire. Au début, il y a beaucoup de petites bulles. Avec le temps, les petites bulles fusionnent pour en former de plus grandes. C'est ce qu'on appelle la coalescence.
• Dans ce liquide, les chercheurs ont vu que les "îles" de l'équipe A et de l'équipe B essayaient de fusionner pour devenir plus grandes et plus stables. Mais plus il fait froid, plus ce processus devient lent.

3. Pourquoi le verre devient-il si dur ? (La viscosité)

C'est ici que la magie opère. La viscosité, c'est la résistance d'un liquide à s'écouler.

• Le mécanisme : Pour que le liquide s'écoule, les molécules doivent changer de place. Mais dans ce modèle, pour qu'une molécule change de place, elle doit traverser la frontière entre les deux "liquides".
• L'analogie : Imaginez que vous essayez de traverser une rivière. Si la rivière est calme, c'est facile. Mais ici, la rivière est remplie de remous et de courants opposés (la tension de surface entre les deux liquides). Plus il fait froid, plus ces remous deviennent violents et difficiles à traverser.
• Les chercheurs ont calculé que la difficulté à traverser cette frontière (l'énergie nécessaire) augmente énormément quand la température baisse. C'est ce qui explique pourquoi la viscosité explose (elle devient des milliards de fois plus grande) juste avant que le liquide ne se fige en verre.

4. La conclusion : Une nouvelle vision du verre

Avant, on pensait que le verre se formait parce que les molécules s'arrêtaient simplement de bouger. Cette étude suggère une autre histoire :
Le verre se forme parce que le liquide est pris dans un jeu de séparation entre deux états. Les molécules sont bloquées en train de essayer de fusionner leurs "îles" respectives, mais le froid les empêche de bouger assez vite.

C'est comme si le liquide était coincé dans une danse où les partenaires essaient de se rapprocher, mais que la musique (la température) ralentit tellement qu'ils finissent par se figer sur place, formant une structure rigide.

En résumé :
Cette recherche montre que la transformation d'un liquide en verre n'est pas un simple ralentissement, mais le résultat d'une lutte interne entre deux formes de liquide qui tentent de se séparer. En mesurant la "tension" entre ces deux formes, les scientifiques ont pu prédire avec précision pourquoi le liquide devient aussi dur et lent, offrant une nouvelle clé pour comprendre la nature mystérieuse du verre.

Résumé technique

Titre : Viscosité dépendante de la température de type Super-Arrhenius due à la séparation de phase liquide-liquide dans le modèle Kob-Andersen surfondue

1. Problématique

La transition vitreuse reste l'un des phénomènes les plus débattus en physique de la matière condensée. Un de ses traits caractéristiques est l'augmentation exponentielle (super-Arrhenius) de la viscosité des liquides surfondus lorsque la température diminue, bien avant d'atteindre la température de transition vitreuse ($T_g$).

• Le défi : Comprendre le mécanisme physique sous-jacent à ce ralentissement dynamique drastique.
• Le contexte : Bien que l'hétérogénéité dynamique soit observée dans les systèmes vitreux, le lien entre cette hétérogénéité et la thermodynamique d'une séparation de phase liquide-liquide (LLT) potentielle n'est pas clairement établi.
• L'objectif : Investiguer si la cinétique de grossissement (coarsening) d'une séparation de phase liquide-liquide peut expliquer la dynamique de la transition vitreuse dans le modèle binaire de Kob-Andersen (KA), un système de référence pour les liquides vitrifiables.

2. Méthodologie

Les auteurs ont employé une approche combinant la dynamique moléculaire (DM) et des modèles statistiques avancés :

• Système simulé : Le modèle binaire de Kob-Andersen (80% de particules A, 20% de particules B) utilisant un potentiel de Lennard-Jones. Les simulations ont été réalisées sous ensembles NVT et NPT avec des conditions aux limites périodiques.
• Paramètre d'ordre : Le Nombre de Coordination Pondéré (WCN) :
  • Au lieu d'utiliser la densité (insuffisante pour distinguer deux phases liquides amorphes similaires), les auteurs utilisent le WCN.
  • Le WCN est un vecteur par particule basé sur les fonctions de distribution radiale $g(r)$, pondéré par des distributions gaussiennes normalisées autour des pics caractéristiques.
  • Une analyse en composantes principales (PCA) est appliquée aux vecteurs WCN pour réduire la dimensionnalité, suivie d'un clustering par la méthode des $k$-means ($k=2$) pour identifier deux états thermodynamiques distincts (liquide 1 et liquide 2).
• Construction du Diagramme de Phase :
  • Reconstruction de la ligne binodale gaz-liquide pour valider la méthode.
  • Extension à la ligne binodale liquide-liquide en quenchant (refroidissement instantané) le système de $T^*=1.2$ vers des températures plus basses ($T^* \in [0.3, 0.585]$).
  • Vérification de la règle des leviers (lever rule) sur les densités et compositions des phases coexistantes pour confirmer l'équilibre thermodynamique local.
• Calcul de la Tension Superficielle et de la Viscosité :
  • Analyse des profils de densité et de pression à travers l'interface pour calculer la tension superficielle ($\sigma_s$) via l'intégration de la différence entre les pressions normale et tangentielle.
  • Utilisation d'un Modèle de Réseau de Markov (MNM) : Le système est modélisé comme deux bassins (domaines) séparés par une barrière d'activation. La viscosité est calculée en fonction de l'énergie d'activation, qui correspond à l'énergie libre interfaciale totale ($\phi$) nécessaire pour le transfert de particules entre les domaines.

3. Contributions Clés

• Développement d'un classificateur robuste : L'utilisation du WCN couplé à la PCA permet de distinguer avec précision deux phases liquides amorphes ayant des densités très proches, là où les paramètres d'ordre classiques échouent.
• Cartographie de la séparation de phase LL : Première reconstruction complète de la ligne binodale liquide-liquide et estimation du point critique ($T_c^* \approx 0.6$, $\rho_c^* \approx 1.1$) pour le modèle KA, démontrant l'existence d'une transition de phase du premier ordre entre deux états liquides.
• Lien mécanique entre interface et viscosité : Démonstration que la viscosité super-Arrhenius peut être expliquée par la cinétique de grossissement de l'interface liquide-liquide. L'énergie d'activation pour le mouvement des particules est directement liée à la tension superficielle de l'interface.
• Validation thermodynamique : Preuve de l'équilibre mécanique et thermodynamique local à l'interface (via les profils de pression et la règle des leviers), même dans le régime de surfusion profonde.

4. Résultats

• Diagramme de Phase : Une ligne binodale liquide-liquide a été établie avec succès. La règle des leviers est respectée pour différentes densités globales, confirmant la coexistence de deux phases liquides distinctes.
• Profils d'Interface :
  • Les profils de densité et de pression montrent un équilibre mécanique pour les interfaces convexes.
  • Pour les interfaces concaves, une légère rupture de l'équilibre mécanique est observée (augmentation locale de la densité et de la pression avant la surface de tension), fournissant une force motrice pour la cinétique de grossissement.
  • La tension superficielle est positive pour les courbures convexes et négative pour les concaves, agissant comme une interférence constructive ou destructive.
• Viscosité et Comportement Super-Arrhenius :
  • En utilisant le modèle de réseau de Markov avec l'énergie d'activation dérivée de la tension interfaciale, les auteurs ont calculé la viscosité.
  • La viscosité augmente d'environ 16 ordres de grandeur lorsque la température passe de $T^*=0.55$ à $T^*=0.36$.
  • Cette augmentation suit un comportement super-Arrhenius, en accord avec les résultats de la littérature (méthode Green-Kubo là où elle est applicable et données expérimentales).
  • Le modèle montre que le ralentissement dynamique est gouverné par la barrière énergétique de la séparation de phase liquide-liquide.

5. Signification et Implications

• Nouveau paradigme pour la transition vitreuse : L'article propose que la transition vitreuse n'est pas simplement un gel cinétique, mais est intrinsèquement liée à la thermodynamique d'une séparation de phase liquide-liquide sous-jacente. Le ralentissement dynamique est une conséquence directe de la difficulté à faire évoluer l'interface entre ces deux phases.
• Universalité potentielle : Les auteurs suggèrent que ce mécanisme pourrait s'appliquer à d'autres systèmes, comme l'eau (modèle ST2), où une séparation de phase liquide-liquide est suspectée. Si confirmé, cela offrirait une image unifiée et simple pour comprendre les transitions vitreuses.
• Méthodologie pour les systèmes complexes : L'approche utilisant le WCN et le modèle de réseau de Markov offre un outil puissant pour étudier la dynamique lente dans des systèmes où les méthodes traditionnelles (comme Green-Kubo) deviennent numériquement impossibles en raison des temps de relaxation extrêmement longs.

En conclusion, cette étude établit un lien quantitatif solide entre la thermodynamique de la séparation de phase liquide-liquide et la dynamique de viscosité super-Arrhenius, suggérant que la cinétique de grossissement des domaines liquides est le moteur principal du ralentissement dynamique menant à la transition vitreuse.