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

Sintesi Tecnica: Teoria della Nucleazione Fuori Equilibrio per Campi Non Conservati

Definizione del Problema
La Teoria della Nucleazione Classica (CNT) descrive con successo la formazione di una fase stabile da una metastabile in sistemi in equilibrio, quantificando la competizione di energia libera tra il guadagno di volume e il costo interfacciale. In equilibrio, la barriera di nucleazione è determinata dal punto di sella dell'energia libera e la traiettoria ottimale (instanton) è la inversione temporale della traiettoria di rilassamento deterministica. Tuttavia, per sistemi privi di dettaglio di bilancio (fuori equilibrio), queste semplificazioni falliscono: gli argomenti basati sull'energia libera non possono stimare le probabilità di eventi rari, e la traiettoria dell'instanton non è nota a priori.

Sebbene una teoria della nucleazione fuori equilibrio (NNT) sia stata sviluppata per sistemi con un parametro d'ordine conservato (ad esempio, Modello Attivo B+), un'estensione rigorosa a sistemi con un parametro d'ordine non conservato è rimasta elusiva. Tali dinamiche non conservate governano molti sistemi di reazione-diffusione, inclusi la dinamica delle popolazioni, i condensati biomolecolari e la materia attiva (ad esempio, Modello Attivo A). La sfida centrale è che, nei sistemi fuori equilibrio, il profilo interfacciale di una goccia in crescita guidata dal rumore devia dal profilo osservato durante il rilassamento deterministico, alterando potenzialmente la barriera di nucleazione in modi non catturati dalle precedenti applicazioni fenomenologiche della CNT.

Metodologia
Gli autori costruiscono una NNT rigorosa per teorie di campo con un singolo parametro d'ordine scalare $\phi(x,t)$ governato dall'equazione stocastica:
$$\partial_t \phi = -\mu[\phi] + \sqrt{2T}D(\phi)\xi$$
dove $\mu$ è il drift deterministico (non necessariamente una derivata funzionale di un'energia libera) e $\xi$ è il rumore bianco gaussiano.

L'innovazione metodologica principale consiste nella sistematica espansione della dinamica stocastica del raggio della goccia $R(t)$ in un'equazione di Langevin 1D:
$$\dot{R} = M(R)[- \partial_R U(R)] + \sqrt{2T M(R)}\zeta$$
dove $U(R)$ è il quasi-potenziale, $M(R)$ è la mobilità emergente e $\zeta$ è il rumore.

I passaggi chiave della derivazione includono:

1. Definizione della Coordinata di Reazione: Gli autori definiscono il raggio della goccia $R(t)$ non arbitrariamente, ma come l'unica soluzione a una condizione integrale che coinvolge una funzione ausiliaria $\psi'(r-R)$. Questa specifica scelta è progettata per proiettare fuori le deviazioni ($\epsilon$) nel profilo interfacciale dall'equazione di Langevin per $R$.
2. Gestione delle Deviazioni Interfacciali: A differenza della CNT in equilibrio o della NNT conservata, gli autori dimostrano che il profilo interfacciale di una goccia in crescita (traiettoria dell'instanton) differisce dal profilo di rilassamento. Essi mostrano che queste deviazioni, denominate $\epsilon_A$, sono piccole ($O(v_0, 1/R, \sqrt{T})$) ma influenzano criticamente la barriera energetica se non gestite correttamente.
3. Selezione di $\psi'$: Scegliendo $\psi'$ in modo che appartenga al nucleo dell'operatore lineare aggiunto ($L^\dagger$), gli autori eliminano la dipendenza dalle incognite deviazioni del campo $\epsilon$ dal termine di drift, permettendo una derivazione analitica di $M(R)$ e $U(R)$.
4. Validazione: I risultati analitici sono validati minimizzando numericamente l'azione di Freidlin–Wentzell della teoria di campo originale utilizzando un metodo di azione minima geometrica (GMAM) combinato con un metodo di Ritz.

Contributi Chiave e Risultati

1. Formule Analitiche per la NNT: L'articolo deriva espressioni esplicite per il raggio critico $R_c$, la mobilità $M(R)$ e la barriera del quasi-potenziale $U(R)$ per sistemi fuori equilibrio non conservati.
   • Il raggio critico è dato da $R_c = \frac{(d-1)\sigma}{\int \psi' \mu_0[\phi_0]}$.
   • La barriera del quasi-potenziale $U(R_c)$ dipende dalla struttura del rumore $D(\phi)$ e dalla specifica scelta di $\psi'$, il che significa che non può essere inferita esclusivamente dalla dinamica priva di rumore o dalla tensione interfacciale $\sigma$ da sola.
2. Fallimento dell'Ansatz di Rilassamento a Tempo Invertito (TRR): Gli autori dimostrano che, per sistemi generici fuori equilibrio, assumere che l'instanton segua una traiettoria di rilassamento invertita nel tempo (TRR) porta a risultati errati. Nello specifico, l'ansatz TRR sovrastima la barriera di nucleazione $U(R_c)$ di un fattore $O(1)$ (ad esempio, di oltre un fattore 2 nel Modello Attivo A), portando a previsioni dei tempi di nucleazione infinitamente troppo lunghi nel limite di basso rumore.
3. Applicazione a Modelli Specifici:
   • Dinamica delle Popolazioni: Applicata a equazioni di tipo Fisher con forti effetti Allee. I risultati mostrano che le barriere di nucleazione dipendono dalla forma specifica del rumore moltiplicativo $D(\phi)$, distinguendosi dalle previsioni di equilibrio.
   • Modello Attivo A (AMA): Applicata a una teoria di campo minima per l'ordinamento di fase fuori equilibrio. La teoria predice correttamente la dipendenza della barriera dal parametro di attività $\lambda$ e dal campo esterno $h$. La minimizzazione numerica dell'azione conferma le previsioni analitiche sia per $R_c$ che per $U(R_c)$, mentre l'ansatz TRR fallisce significativamente.
4. Distorsione del Profilo Interfacciale: Il documento quantifica la deviazione tra il profilo dell'instanton e il profilo di rilassamento, mostrando che scala come $\epsilon_A \sim O(v_0)$, in linea con gli assunti analitici.

Significato
Il lavoro sostiene di fornire la prima estensione ben fondata della teoria della nucleazione alle teorie di campo non conservate fuori equilibrio. La sua importanza primaria risiede nel correggere il malinteso secondo cui la barriera di nucleazione nei sistemi fuori equilibrio possa essere approssimata invertendo temporalmente il percorso di rilassamento. Gli autori mostrano che le deviazioni nel profilo interfacciale sono un effetto di primo ordine nel calcolo della barriera per la dinamica non conservata.

Definendo accuratamente la coordinata di reazione per proiettare fuori queste deviazioni, gli autori recuperano un'equazione di Langevin chiusa che fornisce il corretto quasi-potenziale. Questo framework permette la previsione analitica dei tassi di nucleazione in sistemi che spaziano dalla materia attiva alla dinamica ecologica delle popolazioni, dove il dettaglio di bilancio è violato. Il lavoro colma il divario tra le applicazioni fenomenologiche della CNT nella materia attiva e la meccanica statistica rigorosa, fornendo un benchmark per gli algoritmi numerici di eventi rari e suggerendo che approcci simili potrebbero essere applicabili ai condensati biomolecolari e ai sistemi a più parametri d'ordine.