Non-uniqueness of the steady state for run-and-tumble particles with a double-well interaction potential

Cette étude démontre que, contrairement aux particules browniennes, les particules en marche aléatoire avec retournement (RTP) interagissant via un potentiel double puits peuvent présenter une non-unicité de l'état stationnaire à l'infini, caractérisée par des solutions multiples, une transition de connectivité du support et des états asymétriques.

L'essentiel

Imaginez une foule de petites créatures microscopiques, comme des bactéries, qui décident de se déplacer seules dans un long couloir. C'est le monde des particules "Run-and-Tumble" (courir et chavirer). Elles avancent tout droit pendant un moment, puis, de manière aléatoire, elles changent brusquement de direction.

Dans ce papier, les auteurs étudient ce qui se passe quand on met des milliers de ces créatures dans un couloir, en les obligeant à interagir entre elles via une force spéciale : une sorte de "potentiel à double puits".

Voici l'explication de leurs découvertes, simplifiée et imagée :

1. Le décor : Un paysage en forme de "W"

Imaginez le sol du couloir n'est pas plat, mais ressemble à la lettre W (ou deux vallées séparées par une petite colline au milieu).

• Près du centre : Les particules se repoussent (comme si la colline les empêchait de se toucher).
• Loin du centre : Elles s'attirent (comme si les vallées les appelaient).

Normalement, si ces particules étaient des billes classiques (chauffées par l'agitation thermique), elles se répartiraient de manière uniforme et prévisible. Mais ici, elles sont "actives" : elles ont leur propre énergie et leur propre volonté de bouger.

2. La grande surprise : La foule peut se diviser en deux camps

Les auteurs ont découvert que, selon la force de l'attraction et la vitesse à laquelle les particules changent de direction, le système peut se comporter de deux manières très différentes :

• Le groupe uni : Toutes les particules restent ensemble dans une seule grande masse, occupant tout le couloir disponible.
• La séparation (Bifurcation) : Soudainement, la foule se scinde en deux groupes distincts qui s'éloignent l'un de l'autre, laissant un espace vide au milieu. C'est comme si une foule de gens, au lieu de se tenir ensemble, décidait soudainement de former deux rangées séparées par une grande distance, sans aucune raison extérieure.

3. Le mystère de l'histoire (La non-unicité)

C'est ici que ça devient vraiment fascinant. Dans le monde classique (les billes), l'état final est toujours le même, peu importe comment vous commencez. Mais ici, l'histoire compte.

• Si vous commencez avec les particules regroupées au centre, elles peuvent rester ensemble.
• Si vous commencez avec deux petits groupes déjà séparés, ils peuvent rester séparés.

Même si les conditions extérieures (la température, la force d'attraction) sont exactement les mêmes, le système peut choisir l'un ou l'autre état. C'est comme si vous aviez deux états stables possibles pour la même situation, et que le choix dépendait de votre "mémoire" initiale. C'est ce qu'on appelle la bistabilité.

4. La rupture de symétrie : L'inégalité des groupes

Dans le cas classique, si vous avez deux groupes, ils sont toujours de la même taille (50/50). Mais avec ces particules actives, les auteurs ont découvert quelque chose de nouveau : les deux groupes peuvent avoir des tailles différentes !

Imaginez une foule qui se divise en deux, mais où un groupe a 60% des gens et l'autre seulement 40%. Ce déséquilibre peut être stable et durer indéfiniment. C'est une rupture de la symétrie parfaite, rendue possible par le "bruit" actif de ces particules.

5. Pourquoi est-ce important ?

Ce papier montre que le monde des particules "actives" (qui consomment de l'énergie pour bouger, comme les bactéries ou les robots microscopiques) est beaucoup plus riche et imprévisible que celui des objets passifs.

• Analogie finale : Pensez à une réunion de famille.
  • Cas classique (billes) : Tout le monde se rassemble autour de la table à manger, peu importe qui est arrivé en premier.
  • Cas actif (ce papier) : Selon l'humeur initiale et la façon dont les gens sont arrivés, la famille peut décider de rester autour de la table, OU se diviser en deux groupes qui discutent dans des pièces différentes, OU former un grand groupe d'un côté et un petit de l'autre. Et une fois cette décision prise, elle peut rester stable pendant très longtemps.

En résumé : Ce travail révèle que les systèmes actifs peuvent avoir plusieurs "destins" stables différents pour les mêmes conditions, et qu'ils peuvent briser la symétrie naturelle, créant des structures complexes et déséquilibrées qui n'existent pas dans la physique classique. C'est une nouvelle façon de comprendre comment la matière vivante ou active s'organise.

Résumé technique

1. Problématique et Contexte

L'article étudie le comportement d'un système de $N$ particules actives de type « marche et retournement » (Run-and-Tumble Particles, RTP) en une dimension, interagissant via un potentiel pair à double puits. Ce modèle vise à comprendre les phénomènes hors équilibre dans les systèmes actifs, en particulier la formation d'états liés et les transitions de phase qui n'ont pas d'équivalent dans les systèmes browniens passifs.

Le potentiel d'interaction choisi est de la forme :
$$W(r) = -\frac{k_0 r^2}{2} + \frac{g r^4}{4}$$
Ce potentiel est répulsif à courte distance (si $k_0 > 0$) et attractif à longue distance. L'objectif principal est de déterminer la densité stationnaire $\rho_s(x)$ dans la limite thermodynamique ($N \to +\infty$) et d'analyser l'unicité de cet état stationnaire.

Contrairement aux particules browniennes qui convergent vers un unique état d'équilibre de Gibbs, les systèmes actifs peuvent présenter des comportements complexes tels que la séparation de phase induite par la motilité (MIPS) et, comme le montre cet article, une non-unicité de l'état stationnaire (bistabilité) et une brisure de symétrie de parité.

2. Méthodologie

Les auteurs utilisent une approche analytique basée sur l'équation auto-cohérente pour la densité stationnaire, dérivée précédemment dans la limite $N \to \infty$ pour les RTPs en l'absence de bruit thermique ($T=0$).

• Équation auto-cohérente : La densité totale $\rho_s(x)$ satisfait une équation intégrale où le champ de force effectif $\tilde{F}(x)$ dépend lui-même de la densité :
  $$\rho_s(x) = \frac{K}{v_0^2 - \tilde{F}(x)^2} \exp\left( 2\gamma \int^x dz \frac{\tilde{F}(z)}{v_0^2 - \tilde{F}(z)^2} \right)$$
  avec $\tilde{F}(x) = -\int dy W'(x-y)\rho_s(y)$.
• Paramètre renormalisé : Pour le potentiel polynomial donné, le champ de force s'exprime en fonction d'un paramètre d'interaction « renormalisé » $k = k_0 - 3m_2$, où $m_2$ est le second moment de la densité. Cela permet de résoudre analytiquement l'équation pour un $k$ donné, puis de déterminer $k_0$ a posteriori.
• Analyse des supports : L'étude se concentre sur la structure du support de $\rho_s(x)$ (intervalle unique ou union d'intervalles disjoints) en fonction de $k$ et du taux de retournement $\gamma$.
• Validation numérique : Les prédictions analytiques sont validées par des simulations numériques directes des équations du mouvement pour $N=100$ particules, montrant une excellente concordance même pour des systèmes de taille modeste.

3. Résultats Principaux

A. Transition de Support (Connecté vs Déconnecté)

Il existe une valeur critique $k_c = 3/2^{2/3}$ du paramètre renormalisé qui marque une transition dans la topologie du support de la densité :

• Régime connecté ($k < k_c$) : Le support est un intervalle unique $[-y_1, y_1]$. La densité est continue et symétrique.
• Régime déconnecté ($k > k_c$) : Le support se scinde en deux intervalles disjoints $[-y_1, y_3] \cup [-y_3, y_1]$. Les particules forment spontanément deux groupes spatialement séparés.
• Comportement critique : À la transition ($k=k_c$), la densité au centre $x=0$ s'annule de manière essentielle (singularité essentielle), et le support reste déconnecté mais avec une densité nulle au centre.

B. Bistabilité et Non-Unicité de l'État Stationnaire

C'est la découverte majeure de l'article. La relation entre le paramètre physique $k_0$ et le paramètre renormalisé $k$ n'est pas toujours bijective.

• Pour des taux de retournement $\gamma$ faibles ($\gamma < \gamma_c \approx 0.179$), la fonction $k_0(k)$ devient non monotone près de $k_c$.
• Conséquence : Pour une même valeur de $k_0$ (dans un certain intervalle), il existe deux solutions stables distinctes pour la densité stationnaire : l'une avec un support connecté et l'autre avec un support déconnecté.
• L'état atteint dépend des conditions initiales et de la réalisation du bruit, ce qui contraste fortement avec l'unicité de l'équilibre thermodynamique pour les particules browniennes.

C. États Stationnaires Asymétriques (Brisure de Parité)

Dans le régime à support déconnecté ($k > k_c$), les auteurs montrent qu'il est possible d'obtenir des états stationnaires non symétriques ($\rho_s(-x) \neq \rho_s(x)$).

• Ces états sont caractérisés par un troisième moment non nul $m_3$, correspondant à une répartition inégale des particules entre les deux groupes (par exemple, plus de particules à gauche qu'à droite).
• Le paramètre d'asymétrie peut prendre une continuité de valeurs dans un intervalle déterminé, dépendant des conditions initiales.
• Cela démontre une brisure spontanée de la symétrie de parité dans l'état stationnaire, un phénomène absent dans les systèmes passifs à température finie.

D. Comportement aux Bords

L'article analyse en détail le comportement de la densité aux bords du support ($\pm y_1$ et $\pm y_3$). Selon la valeur de $\gamma$ par rapport à des courbes critiques $\gamma_1(k)$ et $\gamma_3(k)$, la densité peut soit diverger (singularité de puissance), soit s'annuler algébriquement.

4. Signification et Implications

• Nature de l'activité : La non-unicité de l'état stationnaire et la brisure de symétrie sont des signatures directes du bruit actif (bruit télégraphique borné) et de la persistance du mouvement. Ces phénomènes disparaissent dans la limite brownienne ($\gamma, v_0 \to \infty$ avec $T_{eff}$ fixe).
• Robustesse : Bien que le modèle utilise un potentiel polynomial qui diverge à l'infini (peu réaliste physiquement), les auteurs montrent via des simulations (Annexe C) que la phénoménologie (transition connecté/déconnecté, bistabilité) persiste pour des potentiels plus réalistes décroissant à l'infini.
• Nouveaux états de la matière : Le système forme des « super-molécules actives » qui peuvent se scinder en deux groupes. La possibilité d'avoir plusieurs états stables pour les mêmes paramètres de contrôle ouvre la voie à des mémoires hystérétiques dans les systèmes actifs.
• Limites de taille finie : Les simulations montrent que ces effets (bistabilité, asymétrie) sont observables dès $N=100$, et que la brisure de symétrie apparaît même pour $N=5$, suggérant que ces phénomènes ne sont pas de simples artefacts de la limite thermodynamique.

En conclusion, cet article fournit une solution analytique exacte pour un modèle d'interactions actives complexes, révélant une richesse phénoménologique inédite (bistabilité, états asymétriques continus) qui défie les intuitions issues de la physique statistique d'équilibre.