[论文解读] Information Geometry of q-Gaussian Densities and Behaviors of Solutions to Related Diffusion Equations

Atsumi Ohara, Tatsuaki Wada|ArXiv.org|Oct 3, 2008

Statistical Mechanics and Entropy参考文献 34被引用 7

一句话总结

该论文运用信息几何方法分析多孔介质方程（PME）及其关联的非线性福克-普朗克方程（NFPE）的解，通过将解嵌入q-高斯密度的流形中进行研究。结果表明，q-高斯流形是不变且吸引的，解轨迹通过保持矩的m-投影被投影到该流形上，且沿测地线演化；自相似（巴伦布拉特-帕特尔）解在几何上具有特殊性，表现为双重自平行曲线，且收敛到该流形的速率是$1/\sqrt{1+t}$，快于收敛到自相似解的$1/t^\beta$速率。

ABSTRACT

This paper presents new geometric aspects of the behaviors of solutions to the porous medium equation (PME) and its associated equation. First we discuss the Legendre structure with information geometry on the manifold of generalized exponential densities. Next by considering such a structure in particular on the q-Gaussian densities, we derive several physically and geometrically interesting properties of the solutions. They include, for example, characterization of the moment-conserving projection of a solution, evaluation of evolutional velocities of the second moments and the convergence rate to the manifold in terms of the geodesic curves, divergence and so on.

研究动机与目标

• 通过信息几何分析多孔介质方程（PME）及其关联的非线性福克-普朗克方程（NFPE）解的长期行为。
• 确立q-高斯族$\mathcal{M}$作为两个方程的不变且吸引的流形。
• 利用m-投影、测地线和Bregman散度等几何工具表征解的演化过程。
• 揭示自相似（巴伦布拉特-帕特尔）解在q-高斯流形内的几何意义。
• 通过信息论散度量化解收敛到流形$\mathcal{M}$的速率。

提出的方法

• 为q-高斯密度的流形配备与广义系综统计相容的对偶信息几何结构。
• 使用m-投影以保持一阶与二阶矩，实现对演化解的矩保持近似。
• 通过m-联络与e-联络定义$\mathcal{M}$上的测地线曲线，证明自相似解为双重自平行（m-与e-测地线）。
• 采用Bregman散度量化解收敛到流形$\mathcal{M}$的速率。
• 应用信息几何中的毕达哥拉斯关系，将解与投影之间、投影与自相似轨迹之间的散度进行分解。
• 利用Csiszár-Kullback不等式推导出解收敛到$\mathcal{M}$的速率$1/\sqrt{1+t}$。

实验结果

研究问题

• RQ1 PME与NFPE的解相对于q-高斯族$\mathcal{M}$的行为如何？
• RQ2 q-高斯流形$\mathcal{M}$是否对PME与NFPE的动力学具有不变性与吸引性？
• RQ3 能否通过Bregman散度等几何工具表征二阶矩的演化及解收敛到$\mathcal{M}$的过程？
• RQ4 自相似（巴伦布拉特-帕特尔）解在$\mathcal{M}$内的几何性质是什么？
• RQ5 解收敛到$\mathcal{M}$的速率与收敛到自相似解的速率相比如何？

主要发现

• PME与NFPE解的q-高斯密度流形$\mathcal{M}$是不变且吸引的。
• 解在$\mathcal{M}$上的m-投影保持其一阶与二阶矩，且对应于$\mathcal{M}$上的一条测地线曲线。
• 解收敛到$\mathcal{M}$的速率为$1/\sqrt{1+t}$，快于$1<t\leq2$时收敛到自相似解的$1/t^\beta$速率。
• 自相似解在$\mathcal{M}$上同时是m-测地线与e-测地线，因此为双重自平行，具有几何上的特殊性。
• 存在Bregman散度的毕达哥拉斯分解：$\mathcal{D}[u\|u^{\rm BP}] = \mathcal{D}[u\|\hat{u}] + \mathcal{D}[\hat{u}\|u^*] + \mathcal{D}[u^*\|u^{\rm BP}]$，其中$\hat{u}$为m-投影，$u^*$为$\hat{u}$在自相似轨迹上的m-投影。
• 二阶矩的演化速度与从$\mathcal{M}$出发的Bregman散度成正比，表明其具有独特的慢扩散行为。