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

Technische Samenvatting: Positiviteitsbehoud in Hogere-orde Benaderingen van de Dean–Kawaki-vergelijking

Probleemstelling
De Dean–Kawaki-vergelijking beschrijft de evolutie van de empirische maat van $N$ onafhankelijke Brownse deeltjes. Hoewel deze exact is op microscopisch niveau als een martingaaloplossing voor een atomaire willekeurige maat, is de vergelijking te singular om geïnterpreteerd te worden als een dichtheids-waardige Stochastische Partiële Differentiaalvergelijking (SPDE) vanwege de niet-Lipschische wortelverdeling van de ruiscoëfficiënt en de schaal-supercriticaliteit. Om het deeltjessysteem met een continuümmodel te benaderen, is regularisatie vereist.

Een cruciale uitdaging in deze benadering is positiviteitsbehoud. Aangezien de empirische dichtheid niet-negatief moet zijn, zou een ideale numerieke of analytische benadering deze eigenschap moeten behouden. Bestaande hogere-orde methoden kampen echter met een afweging:

1. Eindige-verschil-benaderingen (bijv. [CF23, CFIR26]) bereiken een hoge-orde zwakke convergentie in hoge-densiteitsregimes, maar falen in het behoud van positiviteit vanwege de niet-lokale aard van de divergentie-operator, waardoor oplossingen negatieve waarden kunnen aannemen.
2. Positiviteitsbehoudende benaderingen (bijv. [DKP24, DJP26]) trachten de ruis te trunceren of de singulariteit te regulariseren om niet-negativiteit te garanderen, maar lijden onder zwakkere convergentiesnelheden of vereisen sterke ruis-truncatie die de kwaliteit van de benadering beperkt.

Dit artikel onderzoekt de specifieke impact van het schenden van positiviteitsbehoud op de zwakke benaderingsfout, in het bijzonder in lage-densiteitsregimes waar de initiële dichtheid lokaal kan verdwijnen. Er wordt gevraagd of niet-positiviteitsbehoudende benaderingen superpolynomiale convergentiesnelheden kunnen bereiken en hoe de fout zich gedraagt wanneer de initiële dichtheid niet af bounded is van nul.

Methodologie
De auteurs analyseren een spectrale regularisatie van de Dean–Kawaki-vergelijking:
$$du_t = \frac{1}{2}\Delta u_t dt + \frac{1}{\sqrt{N}} \nabla \cdot \left( K_\varepsilon * (f(u_t^+) dW_t) \right)$$
waarbij $K_\varepsilon$ een mollifier is met een compact ondersteund Fourier-transformatie, en $u_t^+ = \max(u_t, 0)$, waarbij $f$ de wortelruis nabootst (bijv. $f(x)=\sqrt{x}$). In tegen tegenstelling tot eerdere positiviteitsbehoudende schema's, wordt het Wiener-proces $W$ niet getrunceerd, en werkt de spectrale regularisatie op de gehele ruiscoëfficiënt. Deze keuze staat toe dat de oplossing $u_t$ negatieve waarden aanneemt, wat het behoud van positiviteit doorbreekt, maar biedt het voordeel dat de regularisatie naar de testfunctie via convolutie wordt verplaatst, waarbij de regulariteit van de testfunctie wordt benut.

De analyse steunt op:

1. Dualiteit met de Hamilton–Jacobi–Bellman (HJB)-vergelijking: Het uitbreiden van de dualiteitsresultaten van [KLvR19] naar complex-waardige testfuncties om de momentgenererende functie van de SPDE-oplossing te relateren aan de oplossing van een oneindig-dimensionale Kolmogorov backward-vergelijking.
2. Laplace-transformatie-analyse: Het meten van de zwakke fout via het verschil in Laplace-transformaties, $|E[e^{\langle u_t, \phi \rangle}] - E[e^{\langle \mu_t, \phi \rangle}]|$, wat de momenten van alle orden controleert via de Cauchy-integraalformule.
3. Decompositie van de fout: De totale fout wordt gedecomposeerd in de initiële benaderingsfout, de Fourier-afkapfout en een foutterm die specifiek voortkomt uit het negatieve deel van de oplossing ($u_t^-$).
4. Constructie van een ondergrens: Een specifieke eendimensionale voorbeeld wordt geconstrueerd waarbij de initiële dichtheid verdwijnt op een sub-interval om een ondergrens op de benaderingsfout af te leiden, wat de beperkingen van niet-positiviteitsbehoudende schema's in lage-densiteitsregimes aantoont.

Belangrijkste Bijdragen en Resultaten

1. Uniform-in-Tijd Zwakke Foutbounds:
   Het artikel stelt een uniform-in-tijd zwakke foutbound vast, gemeten via de Laplace-transformatie. De foutbound bestaat uit drie componenten:

   • De initiële benaderingsfout.
   • De Fourier-afkapfout (afhankelijk van de spectrale resolutie $\varepsilon$).
   • De fout van het negatieve deel: Een term proportioneel aan de waarschijnlijkheid en de grootte van het aannemen van negatieve waarden door de oplossing.

2. Hoge-densiteitsregime ($\rho_{\min} > 0$):
   Wanneer de initiële dichtheid begrensd is van onderaf door een positieve waarde ($\rho_0 \ge \rho_{\min} > 0$), is het negatieve deel van de oplossing exponentieel onderdrukt. Bijgevolg bereikt de spectrale benadering een superpolynomiale convergentiesnelheid voor gladde testfuncties ($\phi \in C^\infty$). Specifiek, met een optimale keuze van $\varepsilon \sim N^{-1/d}(\log N)^{3/d}$, vervalt de fout als $O(N^{-1-\alpha/d}(\log N)^{3\alpha/d})$ voor elke $\alpha > 0$. Dit resultaat houdt stand, ondanks het gebrek aan strikt positiviteitsbehoud, mits de dichtheid voldoende hoog is ten opzichte van de ruissterkte.

3. Lage-densiteitsregime ($\rho_{\min} = 0$):
   Wanneer de initiële dichtheid mag verdwijnen, verdwijnt de exponentiële onderdrukking van het negatieve deel. De bovengrens verslechtert naar een polynomiale snelheid: $O(N^{-(3\alpha+d)/(2\alpha+d)})$.
   Cruciaal is dat de auteurs een ondergrens leveren voor een eendimensionaal voorbeeld met $\rho_{\min}=0$. Zij bewijzen dat de zwakke fout in dit regime niet superpolynomial vervalt. Het negatieve deel van de oplossing wordt relatief groot, domineert de fluctuaties en beperkt de kwaliteit van de benadering. Dit suggereert dat benaderingen die geen positiviteit bewaren, in lage-densiteitsgebieden inherent beperkt zijn in hun nauwkeurigheid.

4. Vergelijking met Bestaande Methoden:

   • In dimensies $d=1, 2$ is de bovengrens voor het niet-positiviteitsbehoudende schema slechter dan die van de positiviteitsbehoudende SPDE-benadering in [DKP24].
   • De resultaten suggereren dat hybride schema's (die hoge-orde spectrale benaderingen combineren met deeltjesgebaseerde of positiviteitsbehoudende methoden in lage-densiteitsregio's) de meest efficiënte strategie kunnen zijn.

5. Numerieke Validatie:
   Numerieke experimenten bevestigen de theoretische bevindingen. Simulaties tonen aan dat:

   • In hoge-densiteitsregimes de spectrale Galerkin-benadering superpolynomiale convergentie vertoont voor gladde observabelen.
   • In lage-densiteitsregimes de waarschijnlijkheid dat de oplossing negatief wordt aanzienlijk toeneemt, wat correleert met de theoretische ondergrenzen op de fout.
   • Het verlies van positiviteit zeldzaam is wanneer $\rho_{\min}$ groot is, maar frequent wanneer $\rho_{\min} \approx 0$.

Betekenis en Claims
Het artikel claimt de eerste analyse te bieden van niet-positiviteitsbehoudende benaderingen voor de Dean–Kawaki-vergelijking die lage initiële dichtheden toestaat. De primaire betekenis ligt in het kwantificeren van de prijs van het schenden van positiviteitsbehoud:

• Het demonstreert dat positiviteitsbehoud essentieel is voor het bereiken van superpolynomiale convergentiesnelheden wanneer de initiële dichtheid verdwijnt.
• Het verduidelijkt dat hoewel spectrale regularisatie een hoge-orde nauwkeurigheid kan bereiken in hoge-densiteitsregimes zonder strikte positiviteit, de "negatieve excursies" van de oplossing een dominante foutbron worden in lage-densiteitsregimes, wat superpolynomiale convergentie verhindert.
• Het werk onderscheidt zich van eindige-dimensionale benaderingen (zoals [CF23]) door gebruik te maken van de oneindig-dimensionale Markoviaanse structuur en de dualiteit met de HJB-vergelijking, wat uniforme-in-tijd foutbounds en de behandeling van complex-waardige testfuncties mogelijk maakt.

De auteurs blijven bescheiden over de universaliteit van deze resultaten, waarbij zij opmerken dat de ondergrens is afgeleid van een specifiek eendimensionaal voorbeeld en dat de interactie tussen deeltjes of externe potentialen (die dynamisch lage-densiteitsregio's kunnen creëren) verdere studie vereist, hoewel de huidige analyse een fundamenteel begrip biedt van de benaderingslimieten in dergelijke scenario's.