An Additive-Noise Approximation to Keller-Segel-Dean-Kawasaki Dynamics: Small-Noise Results

1. Il Problema e il Contesto

Il paper si occupa dell'analisi matematica delle fluttuazioni di sistemi di particelle interagenti che subiscono chemiotassi, modellati attraverso l'equazione di Keller–Segel–Dean–Kawasaki (KSDK).
L'equazione KSDK (1.1) descrive l'evoluzione della densità di particelle $\rho$ soggetta a:

• Diffusione e chemiotassi (termine non lineare $\nabla \cdot (\rho \nabla \Phi_\rho)$).
• Rumore stocastico moltiplicativo di tipo Dean-Kawasaki ($\nabla \cdot (\sqrt{\rho}\xi)$), che deriva dalla teoria dell'idrodinamica fluttuante (FHD) come limite di sistemi di particelle con $N$ agenti quando $N \to \infty$.

Sfide principali:

1. Instabilità: È stato dimostrato (Konarovskyi, Lehmann, von Renesse) che le equazioni stocastiche guidate dal rumore di Dean-Kawasaki sono mal poste se il rumore è moltiplicativo e l'interazione è singolare (come nel caso del potenziale di Poisson $G$). Le soluzioni esistono solo se i dati iniziali sono misure empiriche e il rumore è "allineato" perfettamente, una condizione irrealistica per l'approssimazione continua.
2. Singularità: Il potenziale di interazione $G$ (funzione di Green dell'equazione di Poisson) è singolare nell'origine, rendendo l'analisi delle interazioni non locali estremamente difficile, specialmente in dimensioni superiori a 1.
3. Supercriticità: L'equazione è supercritica nel senso della teoria delle SPDE singolari, richiedendo tecniche avanzate di rinormalizzazione.

2. Metodologia e Approccio Proposto

Gli autori propongono un'approssimazione a rumore additivo per superare le difficoltà di ben-postezza del modello originale. Invece del rumore moltiplicativo $\nabla \cdot (\sqrt{\rho}\xi)$, considerano l'equazione (1.4):
$$(\partial_t - \Delta)\rho = -\chi\nabla \cdot (\rho\nabla \Phi_\rho) - \varepsilon^{1/2}\nabla \cdot (\sqrt{\rho_{det}}\xi_\delta)$$
Dove:

• $\rho_{det}$ è la soluzione dell'equazione di Keller-Segel deterministica.
• $\xi_\delta$ è un rumore spaziale regolarizzato (mollificato) con lunghezza di correlazione $\delta$.
• Il rumore è additivo rispetto alla soluzione stocastica, ma la sua ampiezza dipende dalla soluzione deterministica di riferimento.

Strumenti Matematici:

• Decomposizione Paracontrollata (Paracontrolled Calculus): Utilizzata per gestire la bassa regolarità delle soluzioni nello spazio delle distribuzioni (spazi di Hölder-Besov $C^{\alpha+1}$ con $\alpha+1 < -1$). Questo permette di dare senso ai prodotti di distribuzioni singolari.
• Teoria delle Grandi Deviazioni (LDP) e Teorema del Limite Centrale (CLT): Dimostrati utilizzando approcci variazionali (formula di Boué-Dupuis) e tecniche di Freidlin-Wentzell generalizzate.
• Scalabilità Relativa: L'analisi dipende criticamente dalla relazione tra l'intensità del rumore $\varepsilon$ (che tende a 0) e la lunghezza di correlazione $\delta(\varepsilon)$. Gli autori distinguono due regimi:
  1. Regime "Rough" (Grezzo): Scalabilità $\varepsilon \log(\delta^{-1}) \to 0$. Permette risultati in spazi di distribuzione irregolari.
  2. Regime "Regular" (Regolare): Scalabilità più restrittiva $\varepsilon^{1/2}\delta^{-\gamma-2} \to 0$. Permette risultati in spazi di Sobolev più regolari e il controllo di valori negativi.

3. Risultati Chiave

A. Ben-Postezza e Leggi dei Grandi Numeri (LLN)

• Teorema 4.4 (LLN nel regime Rough): Sotto la scalabilità $\varepsilon \log(\delta^{-1}) \to 0$, la soluzione stocastica $\rho^{(\varepsilon)}_\delta$ converge in probabilità alla soluzione deterministica $\rho_{det}$ nello spazio delle distribuzioni $C^{\alpha+1}$.
• Teorema 3.3 (LLN nel regime Regular): Sotto la scalabilità più forte $\varepsilon^{1/2}\delta^{-\gamma-2} \to 0$, la convergenza avviene in spazi di Sobolev $H^\gamma$ con $\gamma \in (-1/2, 0]$.

B. Teorema del Limite Centrale (CLT)

• Teorema 4.6: Sotto la scalabilità $\varepsilon^{1/2} \log(\delta^{-1}) \to 0$, le fluttuazioni di primo ordine $\varepsilon^{-1/2}(\rho^{(\varepsilon)}_\delta - \rho_{det})$ convergono a un processo di Ornstein-Uhlenbeck generalizzato $v$, soluzione di un'equazione lineare stocastica guidata dal rumore $\xi$. Questo conferma che l'approssimazione a rumore additivo cattura correttamente le fluttuazioni gaussiane attorno al limite di campo medio.

C. Principi di Grandi Deviazioni (LDP)

• Teoremi 3.5 e 4.8: Vengono stabiliti principi di grandi deviazioni per la sequenza $\rho^{(\varepsilon)}_\delta$ con velocità $\varepsilon$.
• La funzione di tasso (rate function) $I(\rho)$ è data dall'infimo dell'energia $L^2$ del controllo di Cameron-Martin $h$ tale che $\rho$ soddisfi l'equazione "scheletrica" (skeleton equation):
  $$(\partial_t - \Delta)\rho = -\chi\nabla \cdot (\rho\nabla \Phi_\rho) - \nabla \cdot (\sqrt{\rho_{det}}h)$$
• Significato: Questo risultato collega l'approssimazione proposta al sistema di particelle originale, mostrando che l'errore di approssimazione risiede nella struttura dell'equazione scheletrica (additiva vs non lineare).

D. Controllo di Eventi Indesiderati

• Valori Negativi (Teorema 3.6): Poiché il rumore è additivo, la soluzione $\rho^{(\varepsilon)}_\delta$ può diventare negativa (fisicamente impossibile per una densità). Gli autori dimostrano che, sotto la scalabilità restrittiva $\varepsilon^{1/2}\delta^{-2} \to 0$, la probabilità che la parte negativa della soluzione superi una soglia $\lambda$ è esponenzialmente piccola in $\varepsilon^{-1}$.
• Esplosioni (Blow-up) (Teorema 4.12): Anche se l'equazione deterministica può esplodere in tempo finito (per $\chi$ grande), la probabilità che la soluzione stocastica esploda prima del tempo di esplosione deterministico $T^*$ è esponenzialmente piccola.

4. Contributi Originali e Significato

1. Validazione dell'Approssimazione: Il lavoro fornisce la prima giustificazione rigorosa dell'uso di un'approssimazione a rumore additivo per sistemi di particelle con interazioni singolari (oltre la condizione Ladyzhenskaya-Prodi-Serrin) e non incompressibili. Dimostra che questo approccio riproduce correttamente sia il limite di campo medio (LLN) che le fluttuazioni gaussiane (CLT).
2. Gestione delle Singolarità: L'uso della teoria delle SPDE singolari (paracontrollata) permette di trattare l'interazione di Poisson in 2D, un caso noto per essere critico e difficile da analizzare con metodi classici.
3. Distinzione tra Regimi: Il paper chiarisce la distinzione fondamentale tra il regime "rough" (dove si lavora con distribuzioni molto irregolari e si ottengono CLT e LDP) e il regime "regular" (dove si ottengono risultati più forti sulla regolarità e il controllo dei valori negativi, ma a costo di una scalabilità più restrittiva).
4. Implicazioni per la Fisica: I risultati supportano l'uso di modelli di idrodinamica fluttuante semplificati per simulare sistemi di particelle complesse, offrendo una base teorica per studiare fenomeni come metastabilità, transizioni di fase fuori equilibrio e correlazioni a lungo raggio.

In sintesi, Martini e Mayorcas colmano un gap significativo tra la teoria delle particelle interagenti e le equazioni stocastiche continue, fornendo un modello approssimato ben posto e analiticamente trattabile che conserva le proprietà statistiche fondamentali del sistema originale.