On the role of positivity preservation for high order approximations of the Dean--Kawasaki equation

Sintesi Tecnica: Preservazione della Positività nelle Approssimazioni di Ordine Elevato dell'Equazione di Dean–Kawaki

Enunciato del Problema
L'equazione di Dean–Kawaki descrive l'evoluzione della misura empirica di $N$ particelle browniane indipendenti. Sebbene sia esatta a livello microscopico come soluzione di martingala per una misura casuale atomica, è troppo singolare per essere interpretata come una SPDE (Equazione Differenziale Stocastica Parziale) con densità come valore a causa del coefficiente di rumore radice quadrata non Lipschitziano e della supercriticità di scala. Per approssimare il sistema di particelle con un modello continuo, è necessaria una regolarizzazione.

Una sfida critica in questa approssimazione è la preservazione della positività. Poiché la densità empirica deve essere non negativa, un'approssimazione numerica o analitica ideale dovrebbe preservare questa proprietà. Tuttavia, gli esistenti metodi ad alto ordine affrontano un compromesso:

1. Approcci a differenze finite (ad es., [CF23, CFIR26]) ottengono un'alta convergenza debole nei regimi ad alta densità, ma falliscono nel preservare la positività a causa della natura non locale dell'operatore divergenza, permettendo alle soluzioni di assumere valori negativi.
2. Approcci che preservano la positività (ad es., [DKP24, DJP26]) troncano il rumore o regolarizzano la singolarità per garantire la non negatività, ma soffrono di tassi di convergenza più deboli o richiedono un forte troncamento del rumore che limita la qualità dell'approssimazione.

Questo articolo investiga l'impatto specifico della violazione della preservazione della positività sull'errore di approssimazione debole, in particolare nei regimi a bassa densità dove la densità iniziale può annullarsi localmente. Si pone la domanda se le approssimazioni che non preservano la positività possano raggiungere tassi di convergenza superpolinomiali e come l'errore si comporti quando la densità iniziale non è limitata inferiormente da uno zero.

Metodologia
Gli autori analizzano una regolarizzazione spettrale dell'equazione di Dean–Kawaki:
$$du_t = \frac{1}{2}\Delta u_t dt + \frac{1}{\sqrt{N}} \nabla \cdot \left( K_\varepsilon * (f(u_t^+) dW_t) \right)$$
dove $K_\varepsilon$ è un mollificatore con trasformata di Fourier a supporto compatto, e $u_t^+ = \max(u_t, 0)$, con $f$ che imita la radice quadrata (ad es., $f(x)=\sqrt{x}$). A differenza dei precedenti schemi che preservano la positività, il processo di Wiener $W$ non viene troncato, e la regolarizzazione spettrale agisce sull'intero coefficiente di rumore. Questa scelta permette alla soluzione $u_t$ di assumere valori negativi, rompendo la preservazione della positività, ma offre il vantaggio di spostare la regolarizzazione sulla funzione di test tramite convoluzione, sfruttando la regolarità della funzione di test stessa.

L'analisi si basa su:

1. Dualità con l'equazione di Hamilton–Jacobi–Bellman (HJB): Estensione del risultato di dualità di [KLvR19] a funzioni di test a valori complessi per relazionare la funzione generatrice dei momenti della soluzione della SPDE alla soluzione di un'equazione Kolmogorov backward infinita-dimensionale.
2. Analisi della Trasformata di Laplace: Misurazione dell'errore debole tramite la differenza nelle trasformate di Laplace, $|E[e^{\langle u_t, \phi \rangle}] - E[e^{\langle \mu_t, \phi \rangle}]|$, che controlla i momenti di ogni ordine tramite la formula integrale di Cauchy.
3. Decomposizione dell'Errore: L'errore totale è scomposto in errore di approssimazione iniziale, errore di cut-off di Fourier e un termine di errore derivante specificamente dalla parte negativa della soluzione ($u_t^-$).
4. Costruzione del Limite Inferiore: Viene costruito un esempio specifico monodimensionale in cui la densità iniziale svanisce su un sotto-intervallo per derivare un limite inferiore sull'errore di approssimazione, dimostrando i limiti degli schemi che non preservano la positività.

Contributi Chiave e Risultati

1. Errori Deboli Uniformi nel Tempo:
   Il documento stabilisce un limite superiore dell'errore debole uniforme nel tempo misurato attraverso la trasformata di Laplace. Il limite dell'errore consiste di tre componenti:

   • Errore di approssimazione iniziale.
   • Errore di cut-off di Fourier (dipendente dalla risoluzione spettrale $\varepsilon$).
   • Errore della parte negativa: Un termine proporzionale alla probabilità e all'entità del fatto che la soluzione assuma valori negativi.

2. Regime ad Alta Densità ($\rho_{\min} > 0$):
   Quando la densità iniziale è limitata inferiormente da uno zero ($\rho_0 \ge \rho_{\min} > 0$), la parte negativa della soluzione è esponenzialmente soppressa. Di conseguenza, l'approssimazione spettrale raggiunge un tasso di convergenza superpolinomiale per funzioni di test regolari ($\phi \in C^\infty$). Nello specifico, con una scelta ottimale di $\varepsilon \sim N^{-1/d}(\log N)^{3/d}$, l'errore decade come $O(N^{-1-\alpha/d}(\log N)^{3\alpha/d})$ per ogni $\alpha > 0$. Questo risultato è valido nonostante la mancanza di una stretta preservazione della positività, a condizione che la densità sia sufficientemente alta rispetto alla forza del rumore.

3. Regime a Bassa Densità ($\rho_{\min} = 0$):
   Quando la densità iniziale può annullarsi, la soppressione esponenziale della parte negativa scompare. Il limite superiore peggiora in un tasso polinomiale: $O(N^{-(3\alpha+d)/(2\alpha+d)})$.
   Fondamentalmente, gli autori forniscono un limite inferiore per un esempio monodimensionale con $\rho_{\min}=0$. Dimostrano che l'errore debole non può decadere superpolinomialmente in questo regime. La parte negativa della soluzione diventa relativamente grande, dominando le fluttuazioni e limitando la qualità dell'approssimazione. Ciò suggerisce che, nelle regioni a bassa densità, le approssimazioni che non preservano la positività sono intrinsecamente limitate in termini di accuratezza.

4. Confronto con Metodi Esistenti:

   • In dimensioni $d=1, 2$, il limite superiore per lo schema che non preserva la positività nel regime a bassa densità è peggiore di quello dell'approssimazione SPDE che preserva la positività in [DKP24].
   • I risultati suggeriscono che gli schemi ibridi (che combinano approssimazioni spettrali ad alto ordine nelle regioni ad alta densità con metodi basati su particelle o che preservano la positività nelle regioni a bassa densità) potrebbero essere la strategia più efficiente.

5. Validazione Numerica:
   Gli esperimenti numerici confermano i risultati teorici. Le simulazioni mostrano che:

   • Nei regimi ad alta densità, l'approssimazione di Galerkin spettrale esibisce una convergenza superpolinomiale per osservabili regolari.
   • Nei regimi a bassa densità, la probabilità che la soluzione diventi negativa aumenta significativamente, correlando con i limiti inferiori teorici dell'errore.
   • La perdita di positività è rara quando $\rho_{\min}$ è grande, ma frequente quando $\rho_{\min} \approx 0$.

Significato e Rivendicazioni
L'articolo sostiene di fornire la prima analisi delle approssimazioni che non preservano la positività per l'equazione di Dean–Kawaki che permetta densità iniziali basse. Il suo significato primario risiede nel quantificare il costo della violazione della preservazione della positività:

• Dimostra che la preservazione della positività è essenziale per raggiungere tassi di convergenza superpolinomiali quando la densità iniziale svanisce.
• Chiarisce che, sebbene la regolarizzazione spettrale possa raggiungere un'accuratezza elevata nei regimi ad alta densità senza una stretta positività, le "escursioni negative" della soluzione diventano una fonte di errore dominante nelle regioni a bassa densità, impedendo la convergenza superpolinomiale.
• Il lavoro si distingue dagli approcci a dimensione finita (come [CF23]) utilizzando la struttura markoviana infinita-dimensionale e la dualità con l'equazione HJB, permettendo limiti di errore uniformi nel tempo e il trattamento di funzioni di test a valori complessi.

Gli autori rimangono modesti riguardo alla universalità di questi risultati, notando che il limite inferiore è derivato da un esempio specifico monodimensionale e che l'interazione tra particelle o potenziali esterni (che possono creare regioni a bassa densità dinamicamente) richiede ulteriori studi, sebbene l'attuale analisi fornisca una comprensione fondamentale dei limiti di tali approssimazioni in tali scenari.