Nonequilibrium nucleation theory for nonconserved fields: from active matter to population dynamics

Technische Samenvatting: Niet-evenwicht Nucleatie Theorie voor Niet-geconserveerde Velden

Probleemstelling
Klassieke Nucleatie Theorie (CNT) beschrijft succesvol de vorming van een stabiele fase uit een metastabiele fase in evenwichtssystemen door de vrije-energiecompetitie tussen bulkwinst en grensvlakkosten te kwantificeren. In evenwicht wordt de nucleatiebarrière bepaald door het zadelpunt van de vrije energie, en de optimale trajectorie (instanton) is de tijdsomgekeerde van het deterministische relaxatiepad. Echter, voor systemen zonder gedetailleerd evenwicht (niet-evenwicht), falen deze vereenvoudigingen: argumenten op basis van vrije energie kunnen de waarschijnlijkheid van zeldzame gebeurtenissen niet schatten, en de instanton-trajectorie is niet vooraf bekend.

Hoewel een niet-evenwichts nucleatie Theorie (NNT) is ontwikkeld voor systemen met een geconserveerde ordeparameter (bijv. Active Model B+), is een rigoureuze uitbreiding naar systemen met een niet-geconserveerde ordeparameter tot nu toe uitgebleven. Dergelijke niet-geconserveerde dynamica beheersen veel reactie-diffusiesystemen, waaronder populatiedynamica, biomoleculaire condensaten en actieve materie (bijv. Active Model A). De centrale uitdaging is dat in niet-evenwichtssystemen het grensvlakprofiel van een door ruis gedreven groeiende druppel afwijkt van het profiel dat wordt waargenomen tijdens deterministische relaxatie, wat potentieel de nucleatiebarrière op manieren kan veranderen die niet door eerdere fenomenologische toepassingen van CNT worden gevangen.

Methodologie
De auteurs construeren een rigoureuze NNT voor veldentheorieën met een enkele, scalaire, niet-geconserveerde ordeparameter $\phi(x,t)$ die wordt beheerst door de stochastische vergelijking:
$$\partial_t \phi = -\mu[\phi] + \sqrt{2T}D(\phi)\xi$$
waarbij $\mu$ de deterministische drift is (niet noodzakelijkerwijs een functionele afgeleide van een vrije energie) en $\xi$ Gaussische witte ruis is.

De kern van de methodologische innovatie is de systematische expansie van de stochastische dynamica van de druppelradius $R(t)$ naar een 1D Langevin-vergelijking:
$$\dot{R} = M(R)[- \partial_R U(R)] + \sqrt{2T M(R)}\zeta$$
waarbij $U(R)$ de quasipotentiaal is, $M(R)$ de emergente mobiliteit is, en $\zeta$ ruis is.

Belangrijke stappen in de afleiding zijn:

1. Definitie van de Reactiecoördinaat: De auteurs definiëren de druppelradius $R(t)$ niet arbitrair, maar als de unieke oplossing van een integraalvoorwaarde die een hulpfunctie $\psi'(r-R)$ betreft. Deze specifieke keuze is ontworpen om afwijkingen ($\epsilon$) in het grensvlakprofiel uit de Langevin-vergelijking voor $R$ te projecteren.
2. Omgang met Grensvlakafwijkingen: In tegenstelling tot evenwichts-CNT of geconserveerde NNT, tonen de auteurs aan dat het grensvlakprofiel van een groeiende druppel (instanton-trajectorie) verschilt van het relaxatieprofiel. Ze laten zien dat deze afwijkingen, aangeduid als $\epsilon_A$, klein zijn ($O(v_0, 1/R, \sqrt{T})$), maar cruciaal de energiebarrière beïnvloeden als ze niet correct worden behandeld.
3. Selectie van $\psi'$: Door $\psi'$ te kiezen zodanig dat deze in de kern ligt van de geadjungeerde geliniariseerde operator ($L^\dagger$), elimineren de auteurs de afhankelijkheid van de onbekende veldafwijkingen $\epsilon$ in de driftterm, wat een analytische afleiding van $M(R)$ en $U(R)$ mogelijk maakt.
4. Validatie: De analytische resultaten worden gevalideerd door de Freidlin–Wentzell actie van de oorspronkelijke veldentheorie numeriek te minimaliseren met behulp van een geometrische minimum actie methode (GMAM) gecombineerd met een Ritz-methode.

Kernbijdragen en Resultaten

1. Analytische Formules voor NNT: Het artikel leidt expliciete expressies af voor de kritische radius $R_c$, de mobiliteit $M(R)$, en de quasipotentiaalbarrière $U(R)$ voor niet-geconserveerde niet-evenwichtssystemen.
   • De kritische radius wordt gegeven door $R_c = \frac{(d-1)\sigma}{\int \psi' \mu_0[\phi_0]}$.
   • De quasipotentiaalbarrière $U(R_c)$ hangt af van de ruisstructuur $D(\phi)$ en de specifieke keuze van $\psi'$, wat betekent dat deze niet uitsluitend kan worden afgeleid uit de ruisloze dynamica of de grensvlakspanning $\sigma$ alleen.
2. Falen van de Time-Reversed Relaxation (TRR) Ansatz: De auteurs demonstreren dat voor generieke niet-evenwichtssystemen, het aannemen dat de instanton een tijdsomgekeerde relaxatie-trajectorie volgt (TRR), leidt tot onjuiste resultaten. Specifiek leidt de TRR-ansatz tot een overschatting van de nucleatiebarrière $U(R_c)$ met een $O(1)$-factor (bijv. door meer dan een factor 2 in Active Model A), wat leidt tot voorspellingen van nucleatietijden die oneindig lang zijn in het lage-ruislimiet.
3. Toepassing op Specifieke Modellen:
   • Populatiedynamica: Toegepast op Fisher-type vergelijkingen met sterke Allee-effecten. De resultaten laten zien dat nucleatiebarrières afhangen van de specifieke vorm van de multiplicatieve ruis $D(\phi)$, wat hen onderscheidt van evenwichtsvoorspellingen.
   • Active Model A (AMA): Toegepast op een minimale veldentheorie voor niet-evenwicht fase-ordening. De theorie voorspelt correct de afhankelijkheid van de barrière van de activiteitsparameter $\lambda$ en het externe veld $h$. Numerieke minimalisatie van de actie bevestigt de analytische voorspellingen voor zowel $R_c$ als $U(R_c)$, terwijl de TRR-ansatz significant faalt.
4. Distortie van het Grensvlakprofiel: Het artikel kwantificeert de afwijking tussen de instanton-profiel en het relaxatieprofiel, waarbij wordt getoond dat dit schaalt als $\epsilon_A \sim O(v_0)$, wat consistent is met de analytische aannames.

Betekenis
Het artikel claimt de eerste goed gefundeerde uitbreiding van nucleatietheorie naar niet-geconserveerde niet-evenwichtsveldtheorieën te bieden. De primaire betekenis ligt in het corrigeren van de misconceptie dat de nucleatiebarrière in niet-evenwichtssystemen benaderd kan worden door de relaxatiepad te tijdsomkeren. De auteurs tonen aan dat afwijkingen in het grensvlakprofiel een leidende orde-effect zijn in de berekening van de barrière voor niet-geconserveerde dynamica.

Door de reactiecoördinaat zorgvuldig te definiëren om deze afwijkingen te projecteren, herstellen de auteurs een gesloten Langevin-vergelijking die de correcte quasipotentiaal oplevert. Dit raamwerk maakt de analytische voorspelling van nucleatiesnelheden mogelijk in systemen variërend van actieve materie tot ecologische populatiedynamica, waar gedetailleerd evenwicht gebroken is. Het werk overbrugt de kloof tussen fenomenologische toepassingen van CNT in actieve materie en rigoureuze statistische mechanica, en biedt een benchmark voor numerieke zeldzame-gebeurtenis algoritmen en suggereert dat vergelijkbare benaderingen van toepassing kunnen zijn op biomoleculaire condensaten en systemen met meerdere ordeparameters.