Non-equilibrium coupling to a diffusing density breaks Ising universality

Technische Samenvatting: Niet-evenwichtscoppeling aan een diffunderende dichtheid verbreekt Ising-universaliteit

Probleemstelling
De Ising-universaliteitsklasse wordt algemeen beschouwd als opmerkelijk robuust tegen niet-evenwichtverstoringen. In standaard renormalisatiegroep (RG)-analyses stromen termen die het principe van gedetailleerd evenwicht schenden doorgaans naar nul bij het vaste punt, waardoor de evenwichts-kritische exponenten behouden blijven. De auteurs onderzoeken echter of deze robuustheid standhoudt wanneer een ordeparameter ($\psi$) niet-reciprook gekoppeld is aan een geconserveerde, diffunderende dichtheid ($\rho$). Een dergelijke opstelling komt voor in systemen zoals bistabiele mobiele deeltjes (bijv. actieve Ising-modellen) of reactienetwerken met mobiele katalysatoren, waarbij de dichtheid onafhankelijk van de ordeparameter evolueert en daarmee het gedetailleerd evenwicht expliciet schendt. De centrale vraag is of deze specifieke niet-evenwichtscoppeling de grootschalige kritische gedragingen verandert, of dat het systeem binnen de Ising-universaliteitsklasse blijft.

Methodologie
De auteurs introduceren het Brownian Ising Model (BIM), een minimaal veldentheorie die de hydrodynamische koppeling beschrijft tussen een $Z_2$ ordeparameter en een geconserveerde scalaire dichtheid. De dynamica wordt beheerst door:

1. Een Model A-achtige vergelijking voor $\psi$ met een massaterm die afhankelijk is van de lokale dichtheid $\rho$.
2. Een puur diffusiële vergelijking voor de dichtheidsfluctuatie $\delta\rho$.
3. Niet-reciprooke koppeling: de dichtheid beïnvloedt de ordeparameter, maar de ordeparameter heeft geen terugkoppeling op de dichtheidsstroom (in tegenstelling tot Model C).

De analyse maakt gebruik van de dynamische renormalisatiegroep met behulp van het Martin-Siggia-Rose padintegraalformalisme. De studie richt zich op het gedrag nabij de bovenste kritische dimensie $d_c = 4$ met behulp van een $\epsilon = 4-d$ expansie. Belangrijke stappen zijn:

• Het afleiden van de bèta-functies voor de dimensieloze koppelingsconstanten: de niet-lineariteit $u_R$, de dichtheid-magnetisatie koppeling $f_R$, en de tijdschaalratio $w_R$.
• Het analyseren van de stroom in het $(u_R, \tilde{f}_R)$ vlak, waarbij $\tilde{f}_R$ een geschaalde koppeling is die eindig blijft in zowel de snelle-diffusielimiet ($w_R \to 0$) als de gekwelde limiet ($w_R \to \infty$).
• Het uitvoeren van één-lus en twee-lus berekeningen om de stabiliteit van vaste punten en kritische exponenten te bepalen.
• Het toepassen van een all-orders gemodificeerde Harris-criterium, afgeleid van de verschuivingssymmetrie van het veldentheorie, om de stabiliteit verder te gaan dan de perturbatietheorie te beoordelen.

Belangrijkste Bijdragen en Resultaten

1. Ontdekking van het BIM-vaste punt:
   De RG-stroom laat zien dat het standaard Ising (Model A) vaste punt instabiel is voor de dichtheidskoppeling. Het systeem stroomt in plaats daarvan naar een nieuw, IR-stabiel vast punt (het BIM-vaste punt) gekenmerkt door:

   • $X_R^* = 0$ (overeenkomend met de snelle-diffusielimiet waarbij de dichtheid veel sneller relaxeert dan de ordeparameter).
   • Niet-nul vaste waarden voor de koppelingsconstanten: $u_R^* = \epsilon/6$ en $\tilde{f}_R^* = \epsilon/8$.

2. Nieuwe Signaturen van de Universaliteitsklasse:
   Het BIM-vaste punt definieert een nieuwe universaliteitsklasse die op twee fundamentele manieren verschilt van Ising:

   • Schending van de Vlootatie-Dissipatie Stelling (FDT): De anomalieën voor de correlatiefunctie ($\eta$) en de responsfunctie ($\eta'$) splitsen uit elkaar.
     • $\eta = -\epsilon/8 + O(\epsilon^2)$ (negatief, wat wijst op lang reikende correlaties die worden versterkt door de ruis).
     • $\eta' = O(\epsilon^2)$ (positief).
     • De ongelijkheid $\eta \neq \eta'$ is een definitieve signatuur van niet-evenwichtskritikaliteit.
   • Gemodificeerde Kritische Exponenten: De correlatielengte-exponent $\nu$ wijkt zelfs bij de één-lus orde af van de Ising-waarde:
     • $\nu = 1/2 + \epsilon/8 + O(\epsilon^2)$.
     • In $d=3$ levert de twee-lus berekening $\nu \approx 0.65$ op.

3. Mechanisme van Instabiliteit:
   De instabiliteit van het Ising-vaste punt komt voort uit het feit dat de snelle-diffusielimiet ($w_R \to 0$) de dichtheidskoppeling transformeert in een lang reikende multiplicatieve ruis met ruimtelijke correlaties die vervallen als $r^{-(d-2)}$. In tegenstelling tot kort reikende ruis, is deze lang reikende ruis RG-relevant nabij $d_c=4$ en drijft het het systeem weg van het Ising-vaste punt.

4. Stabiliteitsanalyse:

   • Bevestiging via twee-lus: Hoewel de één-lus analyse de stabiliteit van het BIM-vaste punt ambigu laat (door een verdwijnende lineaire eigenwaarde in de $X_R$-richting), bevestigt een twee-lus berekening de stabiliteit met een positieve correctie-naar-schaling exponent $\omega$.
   • All-orders criterium: Met behulp van een gegeneraliseerd Harris-criterium, $\nu > 2/(d + z - 2)$, bewijzen de auteurs dat noch het standaard Ising-, noch de gedilueerde Ising-vaste punten stabiel zijn in $d=3$. Het BIM-vaste punt is de enige levensvatbare attractor voor $d \leq 4$.
   • Dimensionale afhankelijkheid: De analyse voorspelt dat Ising-universaliteit wordt hersteld in $d=2$, wat consistent is met bestaande numerieke resultaten.

5. Eindige-grootte Effecten:
   De correctie-naar-schaling exponent $\omega$ blijkt anomal klein te zijn in $d=3$ ($\omega \approx 0.020$). Dit impliceert dat eindige-grootte correcties zeer traag vervallen ($L^{-\omega}$), wat betekent dat numerieke simulaties uitzonderlijk grote systeemgroottes of expliciete correctieanalyses vereisen om asymptotische exponenten betrouwbaar te extraheren.

Betekenis
Dit artikel stelt vast dat de Ising-universaliteitsklasse niet robuust is tegen niet-reciprooke koppeling met een geconserveerde diffunderende dichtheid. Deze koppeling creëert een nieuwe universaliteitsklasse (BIM) die wordt gekenmerkt door:

• Een echte niet-evenwichtsnatuur, gemanifesteerd door de splitsing van correlatie- en respons-exponenten ($\eta \neq \eta'$).
• Een negatieve anomalie $\eta$, wat reflecteert op versterkte ruimtelijke correlaties door lang reikende ruis.
• Een afwijkende set kritische exponenten ($\nu \approx 0.65$ in $d=3$).

De auteurs betogen dat dit mechanisme waarschijnlijk relevant is voor een brede klasse van levende en gedreven systemen, waaronder celpopulaties, neurale assemblages en actieve materie, waarbij bistabiele ordeparameters interageren met fluctuerende omgevingen zonder terugkoppeling. De splitsing van $\eta$ en $\eta'$ wordt geïdentificeerd als de primaire meetbare signatuur voor het detecteren van dergelijke niet-evenwichtskritikaliteit in simulaties en experimenten.