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

Resumo Técnico: Preservação de Positividade em Aproximações de Ordem Superior da Equação de Dean–Kawasaki

Enunciado do Problema
A equação de Dean–Kawasaki descreve a evolução da medida empírica de $N$ partículas brownianas independentes. Embora seja exata ao nível microscópico como uma solução de martingala para uma medida aleatória atômica, ela é excessivamente singular para ser interpretada como uma EDP estocástica (SPDE) de densidade, devido ao coeficiente de ruído de raiz quadrada não Lipschitziano e à supercriticidade de escala. Para aproximar o sistema de partículas com um modelo contínuo, a regularização é necessária.

Um desafio crítico nesta aproximação é a preservação da positividade. Uma vez que a densidade empírica deve ser não negativa, uma aproximação numérica ou analítica ideal deve preservar essa propriedade. No entanto, métodos de ordem elevada existentes enfrentam um compromisso (trade-off):

1. Abordagens de diferenças finitas (ex: [CF23, CFIR26]) alcançam convergência fraca de ordem elevada em regimes de alta densidade, mas falham em preservar a positividade devido à natureza não local do operador divergência, permitindo que as soluções assumam valores negativos.
2. Abordagens de preservação de positividade (ex: [DKP24, DJP26]) truncam o ruído ou regularizam a singularidade para garantir a não-negatividade, mas sofrem de taxas de convergência mais fracas ou requerem um truncamento de ruído forte que limita a qualidade da aproximação.

Este artigo investiga o impacto específico da violação da preservação de positividade no erro de aproximação fraca, particularmente em regimes de baixa densidade, onde a densidade inicial pode se anular localmente. Questiona-se se aproximações que não preservam a positividade podem alcançar taxas de convergência superpolinomiais e como o erro se comporta quando a densidade inicial não é limitada inferiormente por um valor positivo.

Metodologia
Os autores analisam uma regularização espectral da equação 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)$$
onde $K_\varepsilon$ é um suavizador (mollifier) com transformada de Fourier de suporte compacto, e $u_t^+ = \max(u_t, 0)$, onde $f$ mimetiza o ruído de raiz quadrada (ex: $f(x)=\sqrt{x}$). Diferente de esquemas anteriores de preservação de positividade, o processo de Wiener $W$ não é truncado, e a regularização espectral atua em todo o coeficiente de ruído. Esta escolha permite que a solução $u_t$ assuma valores negativos, quebrando a preservação da positividade, mas oferece a vantagem de deslocar a regularização para a função de teste via convolução, explorando a regularidade da função de teste.

A análise baseia-se em:

1. Dualidade com a equação Hamilton–Jacobi–Bellman (HJB): Extensão do resultado de dualidade de [KLvR19] para funções de teste de valores complexos, para relacionar a função geradora de momentos da solução da SPDE à solução de uma equação de Kolmogorov reversa infinita-dimensional.
2. Análise de Transformada de Laplace: Medição do erro fraco através da diferença nas transformadas de Laplace, $|E[e^{\langle u_t, \phi \rangle}] - E[e^{\langle \mu_t, \phi \rangle}]|$, que controla os momentos de todas as ordens via fórmula de Cauchy.
3. Decomposição do Erro: O erro total é decomposto em erro de aproximação inicial, erro de corte de Fourier e um termo de erro decorrente especificamente da parte negativa da solução ($u_t^-$).
4. Construção de Limite Inferior: Um exemplo específico unidimensional é construído onde a densidade inicial se anula em um subintervalo para derivar um limite inferior no erro de aproximação, demonstrando as limitações de esquemas que não preservam a positividade em regimes de baixa densidade.

Principais Contribuições e Resultados

1. Limites de Erro Fraco Uniformes no Tempo:
   O artigo estabelece um limite de erro fraco uniforme no tempo medido através da transformada de Laplace. O limite de erro consiste em três componentes:

   • Erro de aproximação inicial.
   • Erro de corte de Fourier (dependente da resolução espectral $\varepsilon$).
   • Erro da parte negativa: Um termo proporcional à probabilidade e à magnitude de a solução assumir valores negativos.

2. Regime de Alta Densidade ($\rho_{\min} > 0$):
   Quando a densidade inicial é limitada inferiormente por um valor positivo ($\rho_0 \ge \rho_{\min} > 0$), a parte negativa da solução é exponencialmente suprimida. Consequentemente, a aproximação espectral alcança uma taxa de convergência superpolinomial para funções de teste suaves ($\phi \in C^\infty$). Especificamente, com uma escolha ótima de $\varepsilon \sim N^{-1/d}(\log N)^{3/d}$, o erro decai como $O(N^{-1-\alpha/d}(\log N)^{3\alpha/d})$ para qualquer $\alpha > 0$. Este resultado sustenta-se apesar da falta de preservação estrita da positividade, desde que a densidade seja suficientemente alta em relação à força do ruído.

3. Regime de Baixa Densidade ($\rho_{\min} = 0$):
   Quando a densidade inicial pode se anular, a supressão exponencial da parte negativa desaparece. O limite superior deteriora-se para uma taxa polinomial: $O(N^{-(3\alpha+d)/(2\alpha+d)})$.
   Crucialmente, os autores fornecem um limite inferior para um exemplo unidimensional com $\rho_{\min}=0$. Eles provam que o erro fraco não pode decair superpolinomialmente neste regime. A parte negativa da solução torna-se relativamente grande, dominando as flutuações e limitando a qualidade da aproximação. Isso sugere que, em regiões de baixa densidade, aproximações que não preservam a positividade são inerentemente limitadas em precisão.

4. Comparação com Métodos Existentes:

   • Em dimensões $d=1, 2$, o limite superior para o esquema de não-preservação de positividade no regime de baixa densidade é pior do que o da aproximação de SPDE de preservação de positividade em [DKP24].
   • Os resultados sugerem que esquemas híbridos (combinando aproximações espectrais de ordem elevada em regiões de alta densidade com métodos baseados em partículas ou de preservação de positividade em regiões de baixa densidade) podem ser a estratégia mais eficiente.

5. Validação Numérica:
   Experimentos numéricos confirmam os achados teóricos. As simulações mostram que:

   • Em regimes de alta densidade, a aproximação de Galerkin espectral exibe convergência superpolinomial para observáveis suaves.
   • Em regimes de baixa densidade, a probabilidade de a solução se tornar negativa aumenta significadamente, correlacionando-se com os limites inferiores teóricos do erro.
   • A perda de positividade é rara quando $\rho_{\min}$ é grande, mas frequente quando $\rho_{\min} \approx 0$.

Significância e Alegações
O artigo afirma fornecer a primeira análise de aproximações que não preservam a positividade para a equação de Dean–Kawasaki que permite baixas densidades iniciais. Sua principal significância reside na quantificação do custo de violar a preservação da positividade:

• Demonstra que a preservação da positividade é essencial para alcançar taxas de convergência superpolinomiais quando a densidade inicial se anula.
• Esclarece que, embora a regularização espectral possa alcançar alta precisão em regimes de alta densidade sem positividade estrita, as "excursões negativas" da solução tornam-se uma fonte de erro dominante em regiões de baixa densidade, impedindo a convergência superpolinomial.
• O trabalho distingue-se das abordagens de dimensão finita (como [CF23]) ao utilizar a estrutura Markoviana infinita-dimensional e a dualidade com a equação HJB, permitindo limites de erro uniformes no tempo e o tratamento de funções de teste de valores complexos.

Os autores mantêm modéstia quanto à universalidade destes resultados, observando que o limite inferior é derivado para um exemplo unidimensional específico e que a interação entre partículas ou potenciais externos (que podem criar regiões de baixa densidade dinamicamente) requer estudo adicional, embora a análise atual forneça uma compreensão fundamental dos limites de aproximação nesses cenários.