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

기술 요약: Dean–Kawaki 방정식의 고차 근사법에서의 양의 값 보존성 (Positivity Preservation)

문제 정의
Dean–Kawaki 방정식은 $N$개의 독립적인 브라운 운동 입자의 경험적 측도(empirical measure)의 진화를 기술한다. 이 방정식은 미시적 수준에서는 원자적 무작위 측도(atomic random measure)를 갖는 마팅게일 해로서 정확하지만, 비립시츠(non-Lipschitz) 제곱근 노이즈 계수와 스케일링 초임계성(supercriticality)으로 인해 밀도 값을 가진 확률 편미분 방정식(SPDE)으로 해석하기에는 너무 특이하다(singular). 이를 연속체 모델로 근사하기 위해서는 정규화가 필요하다.

이러한 근사 과정에서 핵심적인 과제는 **양의 값 보존(positivity preservation)**이다. 경험적 밀도는 반드시 0 이상이어야 하므로, 이상적인 수치적 또는 해석적 근사는 이 성질을 보존해야 한다. 그러나 기존의 고차 방법들은 다음과 같은 트레이드오프에 직면해 있다:

1. 유한 차분 접근법 (예: [CF23, CFIR26])은 고밀도 영역에서 고차 약수렴(weak convergence)을 달성하지만, 발산 연산자(divergence operator)의 비국소적 성질로 인해 해가 음수 값을 가질 수 있어 양의 값 보존에 실패한다.
2. 양의 값 보존 접근법 (예: [DKP24, DJP26])은 양의 값을 보장하기 위해 노이즈를 절단하거나 특이성을 정규화하지만, 수렴 속도가 낮거나 근사 품질을 제한하는 강한 노이즈 절단이 필요하다.

본 논문은 특히 초기 밀도가 국소적으로 사라질 수 있는 저밀도 영역에서, 양의 값 보존 위반이 약 근사 오차(weak approximation error)에 미치는 구체적인 영향을 조사한다. 즉, 양의 값 보존을 하지 않는 근사법이 초다항식(superpolynomial) 수렴 속도를 달erm할 수 있는지, 그리고 초기 밀도가 0에서 떨어져 있지 않을 때 오차가 어떻게 거동하는지를 질문한다.

방법론
저자들은 다음과 같은 **스펙트럼 정규화(spectral regularization)**된 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)$$
여기서 $K_\varepsilon$는 푸리에 변환이 컴팩트 지지(compactly supported)를 갖는 몰리파이어(mollifier)이며, $u_t^+ = \max(u_t, 0)$이고, $f$는 제곱근 노이즈(예: $f(x)=\sqrt{x}$)를 모사한다. 이전의 양의 값 보존 스킴들과 달리, 위너 프로세스(Wiener process) $W$는 절단되지 않으며, 스펙트럼 정규화는 전체 노이즈 계수에 작용한다. 이러한 선택은 해 $u_t$가 음수 값을 가질 수 있게 하여 양의 값 보존을 깨뜨리지만, 정규화를 테스트 함수에 대한 컨볼루션(convolution)으로 이동시켜 테스트 함수의 정규성을 활용할 수 있다는 장점을 제공한다.

분석은 다음을 기반으로 한다:

1. Hamilton–Jacobi–Bellman (HJB) 방정식과의 쌍대성(Duality): [KLvR19]의 쌍대성 결과를 복소수 값 테스트 함수로 확장하여, SPDE 해의 적률 생성 함수(moment generating function)를 무한 차원 콜모고로프 역방향 방정식(Kolmogorov backward equation)의 해와 연결한다.
2. 라플라스 변환 분석: 라플라스 변환의 차이인 $|E[e^{\langle u_t, \phi \rangle}] - E[e^{\langle \mu_t, \phi \rangle}]|$를 통해 약 오차를 측정하며, 이는 코시의 적분 공식을 통해 모든 차수의 적률을 제어한다.
3. 오차 분해: 전체 오차를 초기 근사 오차, 푸리에 컷오프(Fourier cut-off) 오차, 그리고 음의 부분(negative part)($u_t^-$)에서 기인하는 오차 항으로 분해한다.
4. 하한(Lower Bound) 구축: 초기 밀도가 특정 구간에서 0이 되는 특정 1차원 예시를 구축하여 근사 오차의 하한을 도출함으로써, 양의 값 보존을 하지 않는 스킴의 한계를 입증한다.

주요 기여 및 결과

1. 시간 균등 약 오차 경계 (Uniform-in-Time Weak Error Bounds):
   논문은 라플라스 변환을 통해 측정된 시간 균등 약 오차 경계를 확립한다. 오차 경계는 세 가지 구성 요소로 이루어진다:

   • 초기 근사 오차.
   • 푸리어 컷오프 오차 (스펙트럼 해상도 $\varepsilon$에 의존).
   • 음의 부분 오차: 해가 음수 값을 가질 확률과 그 크기에 비례하는 항.

2. 고밀도 영역 ($\rho_{\min} > 0$):
   초기 밀도가 0에서 떨어져 있을 때 ($\rho_0 \ge \rho_{\min} > 0$), 해의 음의 부분은 지수적으로 억제된다. 결과적으로 스펙트럼 근사는 매끄러운 테스트 함수($\phi \in C^\infty$)에 대해 초다항식 수렴 속도를 달성한다. 구체적으로, $\varepsilon \sim N^{-1/d}(\log N)^{3/d}$로 최적화하여 선택하면, 오차는 임의의 $\alpha > 0$에 대해 $O(N^{-1-\alpha/d}(\log N)^{3\alpha/d})$로 감소한다. 이 결과는 엄격한 양의 값 보존이 결여되어 있음에도 불구하고, 밀도가 노이즈 강도에 비해 충분히 높다면 성립한다.

3. 저밀도 영역 ($\rho_{\min} = 0$):
   초기 밀도가 0이 될 수 있는 경우, 음의 부분의 지수적 억제가 사라진다. 상한(upper bound)은 다항식 속도인 $O(N^{-(3\alpha+d)/(2\alpha+d)})$로 악화된다.
   결정적으로, 저자들은 $\rho_{\min}=0$인 1차원 예시에 대한 하한을 제시한다. 이들은 이 영역에서 약 오차가 초다항식적으로 감소할 수 없음을 증명한다. 해의 음의 부분이 상대적으로 커지며 변동을 지배하고 근사의 품질을 제한하게 된다. 이는 양의 값 보존을 하지 않는 근사법이 저밀도 영역에서는 본질적으로 한계가 있음을 시사한다.

4. 기존 방법과의 비교:

   • 차원 $d=1, 2$에서, 저밀도 영역에서의 비-양의 값 보존 스킴의 상한은 [DKP24]의 양의 값 보존 SPDE 근사보다 나쁘다.
   • 결과는 하이브리드 스킴(고밀도 영역에서는 고차 스펙트럼 근사를 사용하고, 저밀도 영역에서는 입자 기반 또는 양의 값 보존 방법을 사용하는 방식)이 가장 효율적인 전략이 될 수 있음을 시사한다.

5. 수치적 검증:
   수치 실험은 이론적 발견을 확인해 준다. 시뮬레이션 결과는 다음과 같다:

   • 고밀도 영역에서 스펙트럼 갈레르킨(spectral Galerkin) 근사는 매끄러운 관측량에 대해 초다항식 수렴을 보인다.
   • 저밀도 영역에서는 해가 음수가 될 확률이 크게 증가하며, 이는 이론적 하한과 상관관계를 갖는다.
   • $\rho_{\min}$이 클 때는 양의 값 손실이 드물지만, $\rho_{\min} \approx 0$일 때는 빈번하게 발생한다.

의의 및 주장
본 논문은 저밀도 초기 밀도를 허용하는 Dean–Kawaki 방정식에 대한 최초의 비-양의 값 보존 근사 분석을 제공한다고 주장한다. 주요 의의는 양의 값 보존 위반의 비용을 정량화한 데 있다:

• 양의 값 보존이 필수적임을 입증했다: 초기 밀도가 0이 될 때 초다항식 수렴 속도를 달성하기 위해서는 양의 값 보존이 필수적이다.
• 스펙트럼 정규화가 고밀도 영역에서는 양의 값 보존 없이도 높은 정확도를 달성할 수 있지만, **"음의 편차(negative excursions)"**가 저밀도 영역에서는 지배적인 오차 원인이 되어 초다한 수렴을 방해한다는 점을 명확히 했다.
• 이 연구는 무한 차원 마르코프 구조와 HJB 방정식과의 쌍대성을 활용함으로써, 유한 차원 접근법([CF23] 등)과 차별화되어 시간 균등 오차 경계와 복소수 값 테스트 함수의 처리를 가능하게 했다.

저자들은 이 결과의 보편성에 대해 신중한 태도를 유지하며, 하한이 특정 1차원 예시로부터 도출되었다는 점과 입자 간의 상호작용이나 외부 포텐셜(이들이 동적으로 저밀도 영역을 생성할 수 있음)에 대한 추가 연구가 필요함을 언급하였다. 다만, 현재의 분석이 그러한 시나리오에서의 근사 한계에 대한 기초적인 이해를 제공한다는 점을 강조하였다.