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

Resumen Técnico: Preservación de la Positividad en Aproximaciones de Alto Orden de la Ecuación de Dean–Kawasaki

Planteamiento del Problema
La ecuación de Dean–Kawasaki describe la evolución de la medida empírica de $N$ partículas brownianas independientes. Aunque es exacta a nivel microscópico como una solución de Márkov para una medida aleatoria atómica, es demasiado singular para ser interpretada como una EDP estocástica (SPDE) con valores de densidad debido al coeficiente de ruido de raíz cuadrada no Lipschitz y a la supercriticidad de escala. Para aproximar el sistema de partículas con un modelo continuo, se requiere regularización.

Un desafío crítico en esta aproximación es la preservación de la positividad. Dado que la densidad empírica debe ser no negativa, una aproximación ideal, ya sea numérica o analítica, debería preservar esta propiedad. Sin embargo, los métodos de alto orden existentes enfrentan un compromiso:

1. Los enfoques de diferencias finitas (p. ej., [CF23, CFIR26]) logran una convergencia débil de alto orden en regímenes de alta densidad, pero fallan en preservar la positividad debido a la naturaleza no local del operador divergencia, lo que permite que las soluciones tomen valores negativos.
2. Los enfoques que preservan la positividad (p. ej., [DKP24, DJP26]) truncan el ruido o regularizan la singularidad para asegurar la no negatividad, pero sufren de tasas de convergencia más débiles o requieren un truncamiento de ruido fuerte que limita la calidad de la aproximación.

Este artículo investiga el impacto específico de violar la preservación de la positividad en el error de aproximación débil, particularmente en regímenes de baja densidad donde la densidad inicial puede anularse localmente. Se pregunta si las aproximaciones que no preservan la positividad pueden alcanzar tasas de convergencia superpolinomiales y cómo se comporta el error cuando la densidad inicial no está acotada lejos de cero.

Metodología
Los autores analizan una regularización espectral de la ecuación de Dean–Kawasaki:
$$du_t = \frac{1}{2}\Delta u_t dt + \frac{1}{\sqrt{N}} \nabla \cdot \left( K_\varepsilon * (f(u_t^+) dW_t) \right)$$
donde $K_\varepsilon$ es un suavizador con transformada de Fourier de soporte compacto, y $u_t^+ = \max(u_t, 0)$, donde $f$ imita el ruido de raíz cuadrada (por ejemplo, $f(x)=\sqrt{x}$). A diferencia de los esquemas anteriores que preservan la positividad, el proceso de Wiener $W$ no se trunca, y la regularización espectral actúa sobre todo el coeficiente de ruido. Esta elección permite que la solución $u_t$ tome valores negativos, rompiendo la preservación de la positividad, pero ofrece la ventaja de trasladar la regularización a la función de prueba mediante la convolución, explotando la regularidad de la función de prueba.

El análisis se basa en:

1. Dualidad con la ecuación de Hamilton–Jacobi–Bellman (HJB): Extender el resultado de dualidad de [KLvR19] a funciones de prueba de valores complejos para relacionar la función generadora de momentos de la solución de la SPDE con la solución de una ecuación de Kolmogorov hacia atrás de dimensión infinita.
2. Análisis de la Transformada de Laplace: Medir el error débil mediante la diferencia en las transformadas de Laplace, $|E[e^{\langle u_t, \phi \rangle}] - E[e^{\langle \mu_t, \phi \rangle}]|$, lo cual controla los momentos de todos los órdenes mediante la fórmula integral de Cauchy.
3. Descomposición del Error: El error total se descompone en el error de aproximación inicial, el error de corte de Fourier y un término de error que surge específicamente de la parte negativa de la solución ($u_t^-$).
4. Construcción de la Cota Inferior: Se construye un ejemplo específico unidimensional donde la densidad inicial se anula en un subintervalo para derivar una cota inferior en el error de aproximación, demostrando las limitaciones de los esquemas que no preservan la positividad.

Contribuciones Clave y Resultados

1. Límites de Error Débil Uniformes en el Tiempo:
   El artículo establece un límite de error débil uniforme en el tiempo medido a través de la transformada de Laplace. El límite de error consta de tres componentes:

   • Error de aproximación inicial.
   • Error de corte de Fourier (dependiente de la resolución espectral $\varepsilon$).
   • Error de la parte negativa: Un término proporcional a la probabilidad y magnitud de que la solución tome valores negativos.

2. Régimen de Alta Densidad ($\rho_{\min} > 0$):
   Cuando la densidad inicial está acotada lejos de cero ($\rho_0 \ge \rho_{\min} > 0$), la parte negativa de la solución se suprime exponencialmente. En consecuencia, la aproximación espectral logra una tasa de convergencia superpolinomial para funciones de prueba suaves ($\phi \in C^\infty$). Específicamente, con una elección óptima de $\varepsilon \sim N^{-1/d}(\log N)^{3/d}$, el error decae como $O(N^{-1-\alpha/d}(\log N)^{3\alpha/d})$ para cualquier $\alpha > 0$. Este resultado se mantiene a pesar de la falta de preservación estricta de la positividad, siempre que la densidad sea suficientemente alta en relación con la fuerza del ruido.

3. Régimen de Baja Densidad ($\rho_{\min} = 0$):
   Cuando se permite que la densidad inicial se anule, la supresión exponencial de la parte negativa desaparece. El límite superior se deteriora a una tasa polinomial: $O(N^{-(3\alpha+d)/(2\alpha+d)})$.
   Crucialmente, los autores proporcionan una cota inferior para un ejemplo unidimensional con $\rho_{\min}=0$. Demuestran que el error débil no puede decaer superpolinomialmente en este régimen. La parte negativa de la solución se vuelve relativamente grande, dominando las fluctuaciones y limitando la calidad de la aproximación. Esto sugiere que, en regiones de baja densidad, las aproximaciones que no preservan la positividad están inherentemente limitadas en su precisión.

4. Comparación con Métodos Existentes:

   • En dimensiones $d=1, 2$, el límite superior para el esquema que no preserva la positividad en el régimen de baja densidad es peor que el de la aproximación SPDE que preserva la positividad en [DKP24].
   • Los resultados sugieren que los esquemas híbridos (que combinan aproximaciones espectrales de alto orden en regiones de alta densidad con métodos basados en partículas o que preservan la positividad en regiones de baja densidad) pueden ser la estrategia más eficiente.

5. Validación Numérica:
   Los experimentos numéricos confirman los hallazgos teóricos. Las simulaciones muestran que:

   • En regímenes de alta densidad, la aproximación de Galerkin espectral exhibe una convergencia superpolinomial para observables suaves.
   • En regímenes de baja densidad, la probabilidad de que la solución se vuelva negativa aumenta significamente, correlacionándose con los límites teóricos inferiores del error.
   • La pérdida de positividad es rara cuando $\rho_{\min}$ es grande, pero frecuente cuando $\rho_{\min} \approx 0$.

Significancia y Reivindicaciones
El artículo afirma proporcionar el primer análisis de aproximaciones que no preservan la positividad para la ecuación de Dean–Kawasaki que permite densidades iniciales bajas. Su principal significancia radica en cuantificar el costo de violar la preservación de la positividad:

• Demuestra que la preservación de la positividad es esencial para lograr tasas de convergencia superpolinomiales cuando la densidad inicial se anula.
• Aclara que, si bien la regularización espectral puede lograr una alta precisión en regímenes de alta densidad sin una positividad estricta, las "excursiones negativas" de la solución se convierten en una fuente de error dominante en regiones de baja densidad, impidiendo la convergencia superpolinomial.
• El trabajo se distingue de los enfoques de dimensión finita (como [CF23]) al utilizar la estructura de Márkov de dimensión infinita y la dualidad con la ecuación HJB, permitiendo límites de error uniformes en el tiempo y el tratamiento de funciones de prueba de valores complejos.

Los autores mantienen la modestia respecto a la universalidad de estos resultados, señalando que la cota inferior se deriva para un ejemplo unidimensional específico y que la interacción entre partículas o potenciales externos (que pueden crear regiones de baja densidad dinámicamente) requiere mayor estudio, aunque el análisis actual proporciona una comprensión fundamental de los límites de la aproximación en tales escenarios.