Exact mean-field phase diagram for self-avoiding active particles in a lattice

Résumé Technique : Diagramme de Phase à Champ Moyen Exact pour des Particules Actives Auto-Évitantes sur un Réseau

Énoncé du Problème
L'article étudie la séparation de phases induite par la motilité (MIPS) dans un gaz de réseau de particules auto-propulsées soumises à une exclusion de cœur dur. Contrairement aux études précédentes qui reposent largement sur des simulations stochastiques ou des équations hydrodynamiques de type Langevin, ce travail répond au besoin d'une approche probabiliste complémentaire basée sur des équations de maître. Le système est composé de $N$ particules actives sur un réseau de Bravais générique (linéaire, carré, hexagonal, cubique simple, cubique centré et cubique à faces centrées). Chaque particule possède un directeur interne qui biaise le saut persistant le long des directions de coordination du réseau, tout en subissant une diffusion rotationnelle, parallèlement à une diffusion translationnelle de type thermique. Le défi central consiste à déterminer analytiquement la stabilité de la phase homogène (gaz) et à cartographier la surface spinodale — la limite où cette phase devient localement instable — à travers différentes géométries de réseau sans recourir à la simulation numérique pour l'analyse de la stabilité elle-même.

Méthodologie
Les auteurs adoptent un formalisme d'équation de maître dans une approximation de champ moyen générale. L'état du système est décrit par des probabilités locales $p_{r,s}(t)$, représentant la probabilité qu'un site $r$ soit occupé par une particule possédant un directeur interne pointant dans la direction $s$.

1. Linéarisation : L'équation de maître est linéarisée autour de l'état stationnaire homogène ($p^H_{r,s} = \phi/z$, où $\phi$ est la fraction de remplissage).
2. Analyse de Stabilité : La stabilité de l'état homogène est déterminée par les valeurs propres de la matrice Jacobienne $J_H$. En raison de l'invariance par translation (en supposant des conditions aux limites périodiques), les auteurs appliquent le théorème de Bloch, réduisant le problème à un problème de valeurs propres de type liaison étroite (tight-binding) de dimension $z$ pour chaque vecteur d'onde $\mathbf{k}$ dans la première zone de Brillouin.
3. Expansion de Perturbation : Pour trouver la surface spinodale, les auteurs effectuent une expansion de perturbation des valeurs propres au voisinage de $\mathbf{k} = 0$. Cette expansion permet d'identifier la dépendance d'ordre principal des valeurs propres vis-à-vis du vecteur d'onde, identifiant la condition selon laquelle la partie réelle d'une valeur propre devient positive (indiquant une instabilité).
4. Factorisation Géométrique : L'analyse révèle que l'influence de la géométrie spécifique du réseau entre exclusivement par un seul coefficient, $A$, qui est évalué exactement pour chacun des six réseaux de Bravais considérés.

Contributions Clés et Résultats

• Surface Spinodale en Forme Fermée : Le résultat principal est la dérivation de la surface spinodale sous forme analytique fermée pour six réseaux de Bravais distincts. La condition de stabilité est donnée par l'annulation du terme entre crochets dans l'Équation (19) :
  $$\lambda_0 \approx \left[ -\left(\frac{w_t z + w_a}{2d}\right) + \frac{w_a^2(1-\phi)(1-2\phi)}{w_r z} A \right] k^2$$
  où $w_a$, $w_t$, et $w_r$ sont respectivement les taux de mouvement actif, de diffusion translationnelle et de diffusion rotationnelle. Le coefficient $A$ prend des valeurs spécifiques pour chaque réseau (par exemple, $-1$ pour les réseaux linéaire et carré, $-3$ pour l'hexagonal, $-2$ pour le cubique à faces centrées).
• Rôle de la Géométrie du Réseau : L'étude démontre que si la forme fonctionnelle de la surface spinodale reste cohérente à travers les différents réseaux, la géométrie dicte la localisation précise de la limite d'instabilité uniquement via le coefficient $A$.
• Effet de la Diffusion Translationnelle : L'analyse confirme que la diffusion translationnelle ($w_t$) agit pour stabiliser l'état homogène. L'augmentation de $w_t$ « lisse » l'interface entre les phases denses et diluées et peut finalement supprimer totalement la MIPS en traversant la surface spinodale pour revenir dans la région stable.
• Courants de Probabilité Rotationnelle : Les auteurs calculent explicitement les courants de probabilité rotationnelle associés aux états stationnaires inhomogènes. Ces courants purement rotationnels non nuls constituent une signature distinctive de la rupture du détail du bilan (detailed balance) inhérente aux systèmes actifs, une caractéristique qui n'est pas directement accessible via les descriptions classiques de Langevin ou hydrodynamiques.
• Métastabilité et Temps d'Échappement : L'article caractérise la région métastable (en dehors de la spinodale mais à l'intérieur de la binodale) et analyse le temps moyen d'échappement de l'état homogène instable, montrant que le temps d'échappement diminue à mesure que le système s'enfonce plus profondément dans la région spinodale.

Signification et Revendications
L'article affirme fournir un diagramme de phase à champ moyen « exact », offrant une perspective probabiliste rigoureuse sur la MIPS qui complète les approches dynamiques existantes. En évitant les simulations stochastiques pour l'analyse de la stabilité, ce travail rend transparentes les contributions du mouvement actif et du bruit de type thermique à la séparation de phases. Une contribution théorique significative est la visualisation explicite des courants de probabilité rotationnelle dans l'espace de configuration, qui découlent directement de la rupture du détail du bilan. Les auteurs positionnent leur travail comme une étape fondamentale clarifiant comment la géométrie du réseau influence la MIPS à travers un paramètre unique et exactement calculable, fournissant ainsi un point de référence pour les études futures sur la matière active sur des structures discrètes. Le travail reste modeste dans son ampleur, se concentrant sur le niveau de champ moyen et les interactions de plus proche voisin, tout en notant que la construction pourrait théoriquement être étendue à d'autres réseaux ou à des sauts à plus longue portée.