Quasi-two-dimensional dispersions of Brownian particles with competitive interactions: Dynamical clustering, non-Gaussianity and hydrodynamic correlations

Résumé technique : Dispersions quasi-bidimensionnelles de particules browniennes avec interactions compétitives

Énoncé du problème
L'article traite du comportement dynamique des dispersions quasi-bidimensionnelles (Q2D) de particules browniennes soumises à des interactions compétitives de courte portée attractives (SA) et de longue portée répulsives (LR), appelées systèmes SALR. Bien que les propriétés structurelles de ces systèmes — spécifiquement la formation de clusters de taille finie et de microphases — soient bien documentées dans les contextes 3D et Q2D, une compréhension systématique de leur dynamique reste incomplète. Les questions clés restées ouvertes sont :

1. Comment les interactions hydrodynamiques (IH) influencent-elles la stabilité des clusters et le transport collectif dans les géométries confinées Q2D ?
2. Comment l'interrelation entre le regroupement (clustering) et la diffusion collective anormale caractéristique des systèmes Q2D (où $D(q)$ et $H(q)$ divergent lorsque $q \to 0$) se manifeste-t-elle ?
3. Quels sont les mécanismes régissant la dynamique non-gaussienne et l'hétérogénéité dynamique dans ces phases de clusters ?
4. Comment les effets inertiels à court terme et le retardement hydrodynamique affectent-ils l'établissement des corrélations collectives dans les phases de clusters par rapport aux références de sphères dures ?

Méthodologie
Les auteurs utilisent une approche comparative employant deux techniques de simulation distinctes pour isoler les effets de l'hydrodynamique :

1. Dynamique de Langevin (LD) : Utilisée pour modéliser des systèmes où les interactions hydrodynamiques sont absentes. Cela sert de référence pour la dynamique surcritique (overdamped).
2. Collision Multipartiçule (MPC) : Utilisée pour résoudre explicitement l'hydrodynamique fluctuante, capturant à la fois les régimes inertiels et surcritiques. Le solvant est modélisé par des particules ponctuelles subissant une dynamique de rotation stochastique (SRD) au sein de cellules de collision, assurant la conservation de la masse, de la quantité de mouvement et de l'énergie.

Le système consiste en $N_p$ particules sphériques identiques confinées dans une monocouche planaire au sein d'un solvant Newtonien 3D. Le potentiel d'interaction combine une attraction de type Lennard-Jones généralisée abrupte avec une répulsion électrostatique (Yukawa) écrantée. L'étude fait varier systématiquement la force d'attraction ($\epsilon^*$) et la fraction de surface ($\phi_{2D}$) pour cartographier le comportement de phase.

Outils analytiques clés :

• Fonctions de corrélation de liaison ($B(t)$) et analyse de la durée de vie des clusters ($\tau_b$) pour caractériser l'arrêt dynamique.
• Paramètre d'ordre hexagonal local ($q_6$) et sa fonction de corrélation temporelle ($g_6(t)$) pour évaluer la relaxation structurelle.
• Paramètres non-gaussiens ($\alpha_2$) et la fonction auto-van Hove pour quantifier l'hétérogénéité dynamique.
• Fonctions hydrodynamiques : À la fois la fonction statique $H(q)$ et la fonction dépendante du temps $H(q,t)$, décomposées en parties auto ($H_s$) et distincte ($H_d$), pour sonder la diffusion collective et les corrélations de quantité de mouvement.
• Fonctions de corrélation courant-courant longitudinales ($J(q,t)$) pour étudier la propagation des modes sonores.

Contributions clés et résultats

• Comportement de phase et thermalisation : L'étude confirme quatre phases généralisées : Fluide Dispersé (DF), Clusters d'Équilibre (EC), Percolation Aléatoire (RP) et Percolation de Clusters (CP). La thermalisation est vérifiée, montrant que les états stationnaires sont atteints indépendamment des conditions initiales, sauf à très basses températures où des dérives énergétiques lentes se produisent.

• Autodiffusion et sous-diffusion : À mesure que la force d'attraction augmente, l'autodiffusion est significativement supprimée.

  • Dans les phases fluides (DF, EC, RP), les particules présentent des diffusions browniennes distinctes à court et long terme, ralenties par un regroupement transitoire plutôt que par un encagement (caging).
  • Dans la phase de Percolation de Clusters (CP), une relaxation en deux étapes émerge, correspondant à deux régimes de sous-diffusion distincts, coïncidant avec la formation d'un gel.

• Durée de vie des clusters et critères de transition :

  • On constate que la durée de vie des clusters ($\tau_b$) est plus fortement influencée par la force d'attraction que par la concentration de particules.
  • Deux critères dynamiques sont identifiés pour la transition des phases non-regroupées vers les phases de clusters :
    1. Un seuil dans la durée de vie moyenne des clusters : $\tau_b / \tau_D \approx 10$ sépare les phases DF et EC.
    2. Le temps de relaxation de l'ordre hexagonal local ($\tau_6$) : le croisement $\tau_6 \approx \tau_D$ marque la transition DF–EC.
  • Dans la phase EC, la durée de vie des clusters s'aligne sur la théorie du franchissement de barrière de Kramers, tandis que dans le régime DF, elles suivent la théorie de dissociation de Smoluchowski.

• Dynamique non-gaussienne : Les systèmes SALR Q2D à clusters présentent une dynamique non-gaussienne prononcée. La fonction auto-van Hove dans la phase de clusters d'équilibre affiche une forme approximativement exponentielle, cohérente avec un mécanisme de « diffusivité-diffusante ». Cela indique une hétérogénéité dynamique significative où les particules dans les clusters ont une mobilité restreinte tandis que celles dans la phase diluée restent mobiles.

• Interactions Hydrodynamiques (IH) et Diffusion Collective :

  • Les simulations MPC révèlent que la diffusion collective à grande échelle anormalement augmentée, caractéristique des systèmes Q2D interagissant hydrodynamiquement, est qualitativement préservée dans les dispersions Q2D-SALR.
  • Cependant, la présence de clustering supprime le pic d'ordre à moyen rayon de la fonction hydrodynamique par rapport aux homologues 3D.
  • Crucialement, l'analyse des corrélations hydrodynamiques dépendantes du temps ($H(q,t)$) et des modes sonores indique que le clustering conduit à un début plus précoce des interactions hydrodynamiques que dans les systèmes de référence de sphères dures Q2D. Cela implique que les IH deviennent pertinentes sur des échelles de temps inertielles dans les systèmes SALR à clusters.

Signification et affirmations
L'article affirme fournir une analyse dynamique complète qui comble le fossé entre la séparation de microphases structurelles et les propriétés de transport dans les systèmes SALR confinés. En combinant la LD et la MPC, les auteurs démontrent que :

1. La force d'attraction est le paramètre de contrôle dominant pour le regroupement dynamique, surpassant les effets de la concentration de particules.
2. Les interactions hydrodynamiques jouent un rôle critique dans la dynamique collective des dispersions Q2D-SALR, préservant la diffusion à longue portée anormale tout en modifiant les caractéristiques spectrales de la fonction hydrodynamique.
3. Le clustering accélère la pertinence de l'hydrodynamique, faisant en sorte que les IH se manifestent sur des échelles de temps inertielles, une conclusion qui remet en question l'hypothèse d'IH instantanées dans les descriptions théoriques de ces systèmes.
4. L'étude propose des critères dynamiques (basés sur la durée de vie des clusters et la relaxation de l'ordre hexagonal) pour distinguer les phases dispersées des phases de clusters, complétant ainsi les critères structurels existants.

Le travail souligne la nécessité de prendre en compte à la fois les interactions compétitives et le couplage hydrodynamique pour comprendre pleinement les mécanismes microscopiques du regroupement dynamique et du transport collectif dans les dispersions colloïdales confinées.