[논문 리뷰] Entropy dissipation semi-discretization schemes for Fokker-Planck equations

Shui-Nee Chow, Luca Dieci외 2명|arXiv (Cornell University)|2016. 08. 08.

Geometric Analysis and Curvature Flows참고 문헌 22인용 수 7

한 줄 요약

이 논문은 유한한 그래프 위에서 2-워샤르슈타인 거리계에서의 기울기 흐름 구조를 활용하여 비선형 Fokker-Planck 방정식에 대한 새로운 반디스크리터화 방법을 제안한다. 방법은 공간을 그래프로 디스크리터라이즈하고, 이산 2-워샤르슈타인 거리계를 정의하며, 이산 자유 에너지의 기울기 흐름으로서 형식화함으로써, 양의 해와 총 질량을 유지하면서도 이산 길버트 측도로의 지수 수렴을 보장한다.

ABSTRACT

We propose a new semi-discretization scheme to approximate nonlinear Fokker-Planck equations, by exploiting the gradient flow structures with respect to the 2-Wasserstein metric. We discretize the underlying state by a finite graph and define a discrete 2-Wasserstein metric. Based on such metric, we introduce a dynamical system, which is a gradient flow of the discrete free energy. We prove that the new scheme maintains dissipativity of the free energy and converges to a discrete Gibbs measure at exponential (dissipation) rate. We exhibit these properties on several numerical examples.

연구 동기 및 목표

• 2-워샤르슈타인 거리계에 대한 기울기 흐름 구조를 유지하는 비선형 Fokker-Planck 방정식의 반디스크리터화 방법을 개발한다.
• 유한한 그래프 위의 이산 확률 공간으로 최적 운반 이론을 확장하여, 이산 2-워샤르슈타인 거리계의 일관된 정의를 가능하게 한다.
• 해법이 자유 에너지 소산을 유지하고 이산 길버트 측도로 지수 수렴함을 보장한다.
• 그래프 기반 디스크리터라이제이션을 통해 다양한 경계 조건(예: 흐름이 없는 조건, 주기적 조건)과 일반적인 공간 영역을 다룰 수 있다.
• 양의 해를 유지하고 엔트로피의 볼록성을 계승하는 새로운 로그라플라시안 디스크리터라이제이션을 제안한다.

제안 방법

• 정점 $V$와 간선 $E$로 구성된 유한한 그래프 $G = (V, E)$를 사용하여 공간 영역을 디스크리터라이즈한다. 여기서 정점은 공간 상의 점을, 간선은 국소적 연결성을 나타낸다.
• 정점 $V$에 지지된 이산 확률 측도 공간에 대해 이산 2-워샤르슈타인 거리계를 정의하여 기울기 흐름 형식을 가능하게 한다.
• 이산 자유 에너지 $\mathcal{F}(\rho)$에 대해 이산 거리계에 대한 기울기 흐름으로서 반디스크리터라이제이션 방법을 구성한다.
• 핵심 해법으로서 $\frac{d\rho_i}{dt} = \frac{1}{\Delta x^2} \left[ \sum_{j \in N(i)} \rho_j (F_j(\rho) - F_i(\rho))_+ - \sum_{j \in N(i)} \rho_i (F_i(\rho) - F_j(\rho))_+ \right]$를 제안한다. 여기서 $F_i(\rho) = \frac{\partial}{\partial \rho_i} \mathcal{F}(\rho)$이다.
• 확산 항을 디스크리터라이즈하기 위해 로그라플라시안 항 $\frac{1}{\Delta x^2} \sum_{j \in N(i)} (\log \rho_j - \log \rho_i) g_{ij}(\rho)$를 도입하여, 양의 해와 질량 보존을 보장한다.
• 자유 에너지의 국소 최소화 조건이 엄격하게 만족될 경우, 이산 길버트 측도로의 지수 수렴을 이론적으로 확립한다.

실험 결과

연구 질문

• RQ1 이산적 환경에서 2-워샤르슈타인 거리계에 대한 기울기 흐름 구조를 유지하는 Fokker-Planck 방정식의 반디스크리터라이제이션 방법을 구성할 수 있는가?
• RQ2 유한한 그래프 위에서 기울기 흐름 역학을 지원할 수 있도록 일관된 이산 2-워샤르슈타인 거리계를 정의할 수 있는가?
• RQ3 제안된 해법이 자유 에너지 소산을 유지하고 이산 길버트 측도로 수렴함을 보장하는가?
• RQ4 해법의 수렴 속도는 무엇이며, 적절한 조건 하에서 지수 수렴임을 증명할 수 있는가?
• RQ5 확산 항의 로그라플라시안 디스크리터라이제이션은 수치적 안정성을 향상시키고, 양의 해와 질량 보존과 같은 핵심 물리적 성질을 어떻게 유지하는가?

주요 결과

• 제안된 해법은 자유 에너지가 지수적으로 감소하며, $\mathcal{F}(\rho(t)) - \mathcal{F}(\rho^\infty) \leq e^{-Ct}(\mathcal{F}(\rho^0) - \mathcal{F}(\rho^\infty))$를 만족한다. 여기서 $C > 0$은 평형 상태에서 자유 에너지의 헤시안에 의해 결정되는 상수이다.
• 해법은 이산 길버트 측도 $\rho_i^\infty = \frac{1}{K} \exp\left(-\frac{v_i + \sum_j w_{ij} \rho_j^\infty}{\beta}\right)$로 수렴하며, $K = \sum_i \exp\left(-\frac{v_i + \sum_j w_{ij} \rho_j^\infty}{\beta}\right)$이다.
• 로그라플라시안 디스크리터라이제이션은 해가 $\mathcal{P}_o(G)$에 머물도록 보장하여, 자동으로 양의 해와 총 질량 보존을 유지한다.
• 수치적 예제는 잠재력 $V(x)$가 엄격히 볼록이 아니어도 길버트 측도로의 지수 수렴을 확인하며, 엔트로피의 볼록성이 수행하는 역할을 강조한다.
• 스토케스틱한 반데르폴 및 덴징 진동자에 대해서는, 이 해법이 한계 순환과 일치하는 이중 모드 정적 측도를 성공적으로 포착하여 느린-빠른 운동 역학을 반영한다.
• 점근적 수렴 속도는 길버트 측도에서 자유 에너지의 헤시안의 가장 작은 고유값인 $\lambda_{\mathcal{H}}(\rho^\infty)$에 의해 완전히 결정된다.