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

Dit artikel toont aan dat actieve run-and-tumble deeltjes in een dubbelputpotentiaal meerdere stabiele stationaire toestanden kunnen vertonen, waaronder niet-symmetrische en ongekoppelde dichtheidsprofielen, wat een fundamenteel verschil is met het unieke evenwicht van passieve Brownse deeltjes.

De kern

Stel je voor dat je een grote groep mensen in een lange, rechte gang hebt. Deze mensen zijn geen gewone wandelaars; ze zijn als bacteriën (zoals E. coli). Ze hebben een heel specifiek gedrag: ze rennen een stukje in één richting, stoppen dan plotseling, draaien zich om en rennen in de andere richting. In de wetenschap noemen we dit "Run-and-Tumble" (rennen en draaien).

In dit onderzoek kijken we naar wat er gebeurt als deze rennende mensen niet alleen zijn, maar met elkaar interageren via een speciaal soort "kracht".

Het Magische Krachtveld

Stel je een onzichtbaar krachtenveld voor dat tussen deze mensen werkt:

1. Als ze te dicht bij elkaar komen: Ze duwen elkaar weg (repulsie). Net als mensen die hun persoonlijke ruimte nodig hebben.
2. Als ze ver uit elkaar zijn: Ze trekken elkaar aan (attractie). Net als magneten die elkaar willen vinden.

In de natuurkunde noemen we dit een "dubbelput-potentiaal". Het is alsof er een valkuil is: als je te dichtbij bent, word je weggestoten; als je te ver weg bent, word je teruggetrokken.

Het Grote Experiment: Wat gebeurt er als ze stilstaan?

Normaal gesproken, als je een groep mensen laat staan zonder dat ze actief rennen (alleen maar een beetje trillen door warmte), vinden ze op een gegeven moment één perfecte, stabiele vorm. Ze verdelen zich gelijkmatig. Er is maar één juiste oplossing.

Maar in dit onderzoek laten we de mensen actief rennen. En hier wordt het gek:

1. De "Gespleten" Gang (De Overgang)
Afhankelijk van hoe sterk ze elkaar aantrekken of afstoten, kan er iets vreemds gebeuren.

• Soms blijven ze allemaal in één grote groep bij elkaar in het midden van de gang.
• Maar als de krachten net iets anders zijn, splijt de groep spontaan in tweeën. De ene helft verzamelt zich links, de andere helft rechts, en het midden blijft leeg. Het is alsof de groep plotseling twee aparte "stamgroepen" vormt die elkaar vermijden.

2. Het Verrassende Nieuws: Er is niet altijd één oplossing!
Dit is het belangrijkste ontdekking van dit papier. Bij gewone, passieve mensen (die niet rennen) is er altijd maar één stabiele situatie. Maar bij deze actieve renners kan het gebeuren dat er twee verschillende, stabiele situaties mogelijk zijn voor exact dezelfde omstandigheden.

• Analogie: Stel je een bal op een heuvel voor. Normaal rollt hij altijd naar dezelfde dal. Maar hier is het alsof de heuvel zo vreemd is gevormd dat de bal soms links kan blijven liggen, en soms rechts, afhankelijk van waar je hem precies hebt neergezet. Als je de groep start met meer mensen links, blijven ze links. Start je met meer rechts, dan blijven ze rechts. Beide situaties zijn stabiel. Dit noemen we bistabiliteit.

3. De "Scheve" Groep
Bij gewone systemen is alles altijd symmetrisch (links en rechts zijn gelijk). Maar bij deze actieve renners kan het zijn dat de groep scheef wordt.

• Stel je voor dat er 60% van de mensen links zit en 40% rechts. Dit is een stabiele toestand! Ze rennen en draaien, maar ze komen nooit in evenwicht met elkaar. De groep blijft "scheef" staan. Dit is iets wat je bij passieve deeltjes nooit ziet.

Waarom is dit belangrijk?

De onderzoekers (Léo Touzo en Pierre Le Doussal) hebben dit wiskundig uitgewerkt en gecontroleerd met computersimulaties. Ze ontdekten dat:

• Dit gedrag specifiek is voor actieve materie (dingen die zelf energie verbruiken om te bewegen).
• Het laat zien dat de natuur veel meer "keuzes" kan maken dan we dachten, als de deeltjes actief zijn.
• Het is alsof de groep een soort "super-molecuul" vormt dat kan kiezen tussen verschillende vormen, afhankelijk van hoe het begon.

Samenvatting in één zin

Wanneer een groep actieve rennende deeltjes interacteert met een kracht die ze zowel aantrekt als afstoot, kunnen ze spontaan in twee groepen splitsen, en kunnen ze kiezen tussen verschillende stabiele vormen (soms zelfs scheve vormen), wat betekent dat er geen één "perfecte" eindtoestand is, maar meerdere mogelijkheden die bestaan naast elkaar.

Het is een beetje alsof je een klasje kinderen laat rennen in een gang: afhankelijk van hoe je begint, kunnen ze eindigen als één grote kluwen, of als twee aparte groepen die elkaar vermijden, en beide situaties kunnen even lang blijven bestaan!

Technische samenvatting

Titel: Non-univociteit van de stationaire toestand voor run-and-tumble deeltjes met een dubbelput-interactiepotentiaal

Auteurs: Léo Touzo en Pierre Le Doussal
Context: Laboratoire de Physique de l'École Normale Supérieure en University of Chicago.

---

1. Probleemstelling en Motivatie

Het artikel onderzoekt het gedrag van een systeem van $N$ interactieve Run-and-Tumble Particles (RTPs) in één dimensie. RTPs zijn een model voor actieve materie (zoals E. coli-bacteriën) die zich voortbewegen met een constante snelheid $v_0$ en willekeurig van richting veranderen ("tumble") met een snelheid $\gamma$.

• Interactie: De deeltjes ondergaan een paarsgewijze interactiepotentiaal $W(r) = -k_0 r^2/2 + g r^4/4$. Dit is een dubbelputpotentiaal: afstotend op korte afstanden (voor $k_0 > 0$) en aantrekkend op grote afstanden.
• Doel: Het analytisch bepalen van de stationaire dichtheid $\rho_s(x)$ in de limiet van een groot aantal deeltjes ($N \to \infty$).
• Kernvraag: In tegenstelling tot passieve Brownse deeltjes (die een unieke Gibbs-evenwichtstoestand hebben), vertonen actieve systemen vaak nieuwe fasen. De auteurs willen onderzoeken of er meerdere stationaire toestanden (bistabiliteit) kunnen bestaan en of de symmetrie van de dichtheid kan worden verbroken, zelfs in één dimensie.

2. Methodologie

De auteurs gebruiken een combinatie van analytische theorie en numerieke simulaties:

• Zelfconsistent Vergelijking: Voor $N \to \infty$ en zonder thermische ruis ($T=0$) wordt gebruik gemaakt van de Dean-Kawasaki-methode (uitgebreid naar RTPs). Dit leidt tot een niet-lineaire, zelfconsistent vergelijking voor de stationaire dichtheid $\rho_s(x)$:
  $$\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)$$
  Hierin is $\tilde{F}(x)$ een effectief krachtveld dat afhangt van de dichtheid zelf via de interactiekracht.
• Renormalisatie: Door de specifieke vorm van de potentiaal ($W \propto -k_0 r^2 + g r^4$) kan het krachtveld worden uitgedrukt in termen van een gerenormaliseerde parameter $k = k_0 - 3m_2$, waarbij $m_2$ het tweede moment (variantie) van de dichtheid is. Dit vereenvoudigt het probleem aanzienlijk.
• Analytische Oplossing: De vergelijking wordt opgelost voor verschillende waarden van $k$, wat leidt tot verschillende ondersteuningsgebieden (support) voor de dichtheid.
• Numerieke Validatie: Directe simulaties van de bewegingsvergelijkingen voor $N=100$ deeltjes worden uitgevoerd om de analytische voorspellingen te testen.

3. Belangrijkste Resultaten

A. Overgang van Geconnecteerd naar Gescheiden Support

Er treedt een fase-overgang op bij een kritieke waarde van de gerenormaliseerde parameter $k_c = 3/2^{2/3}$:

1. Geconnecteerd Support ($k < k_c$): De deeltjes vormen één enkele cluster met een continue dichtheid op het interval $[-y_1, y_1]$.
2. Gescheiden Support ($k > k_c$): De binding breekt spontaan op. De deeltjes splitsen zich in twee groepen met een disjunct support (twee gescheiden intervallen: $[-y_1, y_3] \cup [-y_3, y_1]$). Er is een "gat" in het midden waar geen deeltjes zijn.
   • Dit gedrag is uniek voor actieve deeltjes; passieve Brownse deeltjes vertonen dit niet bij $T>0$.

B. Non-Univociteit en Bistabiliteit

Een van de meest opmerkelijke bevindingen is dat de relatie tussen de fysieke parameter $k_0$ en de gerenormaliseerde parameter $k$ niet-uniek kan zijn:

• Voor lage waarden van de tumble-snelheid $\gamma$ (onder een kritieke drempel $\gamma_c \approx 0.179$) is de functie $k(k_0)$ niet-monotoon.
• Dit betekent dat voor één specifieke waarde van $k_0$ er twee verschillende stabiele stationaire toestanden kunnen bestaan (bijvoorbeeld één met geconnecteerd support en één met gescheiden support).
• Het systeem kan in welke van deze toestanden terechtkomt, hangt af van de beginvoorwaarden en de realisatie van het ruisproces. Dit contrasteert sterk met het unieke evenwicht van thermische systemen.

C. Breuk van Pariteitssymmetrie

In het regime met gescheiden support ($k > k_c$) is de stationaire toestand niet noodzakelijk symmetrisch ($\rho_s(x) \neq \rho_s(-x)$):

• De deeltjes kunnen ongelijk verdeeld zijn over de twee clusters (bijv. 60% links, 40% rechts).
• Dit wordt gekarakteriseerd door een derde moment $m_3$ dat een continu bereik van waarden kan aannemen.
• De asymmetrie wordt bepaald door de begincondities en de geschiedenis van het systeem. Numerieke simulaties tonen aan dat dit al optreedt voor kleine $N$ (bijv. $N=5$).

D. Randgedrag en Singulariteiten

De dichtheid vertoont specifieke gedragingen aan de randen van het support:

• Afhankelijk van de verhouding tussen $\gamma$ en de randpositie, kan de dichtheid divergeren (een piek vormen) of verdwijnen (naar nul gaan) aan de randen.
• Er zijn specifieke curves in het $(k, \gamma)$-diagram die deze overgangen markeren.

4. Significatie en Conclusie

• Fundamenteel Inzicht: Het artikel demonstreert dat actieve systemen, zelfs in één dimensie, een veel rijkere fenomenologie vertonen dan passieve systemen. De beperking van het telegrafische ruisproces (bounded noise) is cruciaal voor het ontstaan van meervoudige stationaire toestanden en symmetriebreking.
• Bistabiliteit: De ontdekking van bistabiliteit in een eenvoudig 1D-model met korte- en langetermijninteracties is een belangrijk concept voor het begrijpen van niet-evenwichtsfasenovergangen.
• Relevantie voor Realistische Systemen: Hoewel het model een potentiaal gebruikt die divergeert op oneindig (ongeacht realistisch), tonen aanvullende simulaties (Appendix C) aan dat de fenomenologie (gescheiden support, bistabiliteit) ook geldt voor realistischere interacties die op oneindig verdwijnen.
• Toekomstperspectief: De resultaten suggereren dat "actieve supermoleculen" (gebonden toestanden die spontaan splitsen) en complexe dynamische basins of aantrekkingsgebieden mogelijk zijn in actieve materie.

Samenvattend: Dit werk biedt een exacte analytische beschrijving van een actief systeem dat meervoudige stationaire toestanden, bistabiliteit en spontane symmetriebreking vertoont, wat een scherp contrast vormt met de uniciteit van het thermodynamisch evenwicht in passieve systemen. De voorspellingen worden bevestigd door simulaties voor systemen met slechts 100 deeltjes.