Effect of Weak Non-Conservative Dynamics on Pattern Formation in Scalar Active Matter

Résumé Technique : Effet d'une dynamique non conservative faible sur la formation de motifs dans la matière active scalaire

Énoncé du Problème
Les systèmes de matière active, allant des colonies bactériennes aux colloïdes actifs synthétiques, sont intrinsèquement pilotés loin de l'équilibre thermodynamique par une consommation continue d'énergie. Alors que les cadres théoriques standards tels que le Modèle Actif B (AMB) décrivent la séparation de phases dans la matière active scalaire sous une conservation stricte de la masse, de nombreux systèmes biologiques et synthétiques présentent de faibles violations de la conservation de la masse dues à des processus de naissance-mort, de croissance, de dégradation ou d'échange chimique avec des réservoirs. L'influence spécifique de ces dynamiques non conservatives faibles sur la morphologie à temps long, la sélection de motifs et le comportement de mûrissement des systèmes actifs scalaires reste peu explorée. Ce travail comble la lacune de compréhension de la manière dont ces réactions non conservatives modifient l'évolution non linéaire et les structures en régime stationnaire de la séparation de phases active.

Méthodologie
L'auteur construit un modèle continu minimal en étendant le Modèle Actif B avec un terme de réaction non conservatif faible. Le système est décrit par un champ d'ordre scalaire macroscopique, $\phi(\mathbf{r}, t)$, représentant la densité de particules. La dynamique est régie par l'équation de mouvement :
$$\frac{\partial \phi}{\partial t} = -\Gamma_\phi \nabla \cdot \mathbf{J} + \gamma \phi \left(1 - \frac{\phi}{\phi_c}\right)\left(1 + \frac{\phi}{\phi_c}\right)$$
où $\mathbf{J} = -\nabla \mu + \sqrt{2\Gamma_\phi \phi}\eta$ est le courant, et le potentiel chimique $\mu$ inclut une contribution standard de l'énergie de Landau ainsi qu'un terme actif $\mu_{active} = \zeta |\nabla \phi|^2$ qui brise la symétrie de renversement du temps. Le paramètre $\zeta$ contrôle l'intensité de l'activité, tandis que $\gamma$ (avec $\gamma \leq 0$) contrôle l'intensité de la réaction non conservative.

L'étude emploie deux approches principales :

1. Analyse de Stabilité Linéaire : Pour déterminer l'instabilité initiale et la relation de dispersion $\sigma(q)$, en analysant comment le terme de réaction et l'activité affectent le taux de croissance des perturbations autour d'un état stationnaire homogène.
2. Simulations Numériques : Des simulations étendues sur un réseau de $512 \times 512$ avec des conditions limites périodiques pour observer les morphologies à temps long. Le système est initialisé avec une densité uniforme et évolue vers des états stationnaires.
3. Caractérisation Structurelle : Analyse quantitative utilisant des fonctions de corrélation à deux points à temps égal $C(r)$, des facteurs de structure statiques $S(k)$, des statistiques de clusters (fraction du plus grand cluster $L_c$, nombre de clusters isolés $N_c$, circularité $C$), des distributions de taille de gouttelettes, et le paramètre d'ordre hexatique de type bond-orientational $\psi_6$.

Contributions Clés et Résultats

• Arrêt du Mûrissement et Séparation de Microphases : La principale conclusion est que même de faibles réactions non conservatives ($\gamma < 0$) arrêtent fondamentalement le processus de mûrissement typique des systèmes actifs conservés. Au lieu d'évoluer vers une séparation de phases globale, le système se stabilise dans des états de microphase séparée hors équilibre avec une échelle de longueur intrinsèque bien définie.
• Renversement des Transitions Morphologiques : Dans la limite conservée ($\gamma = 0$), l'augmentation de la magnitude de l'activité négative ($\zeta$) conduit à une transition des gouttelettes vers des structures de type "vers" (worm-like) et enfin vers des motifs labyrinthiques interconnectés (décomposition spinodale). En revanche, la présence d'une réaction non conservative faible inverse cette séquence. À mesure que l'activité augmente, le système passe de motifs labyrinthiques à des structures de type "vers" et finalement à des gouttelettes isolées et polydispersées.
• Rôles Distincts des Paramètres :
  • Activité ($\zeta$) : Contrôle principalement la morphologie globale des domaines et la longueur d'onde caractéristique des motifs. L'augmentation de $|\zeta|$ réduit l'échelle de longueur caractéristique ($\lambda^*$) et conduit la transition des labyrinthes vers les gouttelettes isolées.
  • Coefficient de Réaction ($\gamma$) : Favorise principalement la séparation de microphases et renforce l'ordre hexagonal local. L'augmentation de la magnitude de $\gamma$ supprime la coalescence des domaines, augmente le nombre de clusters isolés et réduit la taille moyenne des gouttelettes, conduisant à une distribution de taille plus étroite.
• Aperçus de la Stabilité Linéaire : L'analyse révèle que le terme de réaction entre dans la relation de dispersion comme une constante indépendante de $q$, décalant la courbe de taux de croissance verticalement sans modifier la position du nombre d'onde le plus instable $q^*$. La contribution active ($\zeta |\nabla \phi|^2$) n'apparaît pas à l'ordre linéaire, ce qui signifie que ni le terme de réaction ni l'activité ne modifient la sélection de la longueur d'onde initiale. Les différences qualitatives de morphologie proviennent de l'évolution non linéaire suivant l'instabilité initiale.
• Émergence de l'Ordre Hexatique : Dans la phase de gouttelettes, le système présente un empilement hexagonal local. Le terme de réaction non conservatif renforce cet ordre hexatique local ($\psi_6 \to 1$) en régulant la distribution de taille des gouttelettes et en facilitant l'homogénéisation de la taille, ce qui permet un empilement ordonné. Ceci est distinct du mûrissement d'Ostwald classique, car la nature non conservée du paramètre d'ordre empêche une croissance illimitée.
• Comportement d'Échelle : Les fonctions de corrélation mises à l'échelle $C(r/\lambda^*)$ se superposent sur une courbe maîtresse universelle pour diverses valeurs de $\zeta$ et $\gamma$, indiquant que si les paramètres rescalent la longueur caractéristique, l'organisation spatiale normalisée sous-jacente reste robuste.

Signification
Cet article démontre que les faibles violations de la conservation de la masse ne sont pas de simples perturbations mais altèrent fondamentalement la dynamique de mûrissement non linéaire de la séparation de phases active. En introduisant un cadre minimal couplant le transport actif avec de faibles réactions non conservatives, l'étude fournit un mécanisme pour stabiliser des états de microphase séparée qui, autrement, mûriraient vers des phases globales. Les résultats suggèrent que les dynamiques non conservatives servent de mécanisme simple mais efficace pour contrôler l'auto-organisation hors équilibre. Ce cadre offre une base pour comprendre la formation de motifs dans les systèmes biologiques (par exemple, les colonies bactériennes avec des cycles de naissance-mort) et pour concevoir des matériaux actifs synthétiques avec des microstructures ajustables, tels que des gouttelettes isolées avec un ordre spécifique, sans nécessiter de forçage externe complexe.