[论文解读] Optimal Transportation Theory with Repulsive Costs

Simone Di Marino, Augusto Gerolin+1 more|arXiv (Cornell University)|Jun 15, 2015

Mathematical Inequalities and Applications参考文献 108被引用 16

一句话总结

本综述为具有排斥性代价的多边际最优传输建立了理论与数值框架，重点研究库仑、谐波及行列式代价函数。结果表明，这些代价下的最优传输计划通过集中质量以实现最大空间分离，且在径向对称条件下会涌现出唯一确定性解，该结论通过使用正则化IPFP与离散优化的数值模拟得到验证。

ABSTRACT

This paper intents to present the state of art and recent developments of the optimal transportation theory with many marginals for a class of repulsive cost functions. We introduce some aspects of the Density Functional Theory (DFT) from a mathematical point of view, and revisit the theory of optimal transport from its perspective. Moreover, in the last three sections, we describe some recent and new theoretical and numerical results obtained for the Coulomb cost, the repulsive harmonic cost and the determinant cost.

研究动机与目标

• 为具有排斥性代价函数的多边际最优传输建立全面的理论与数值框架。
• 分析具有库仑、谐波与行列式代价的多边际传输问题最小化器的结构。
• 建立Monge型解存在且唯一的条件，尤其在径向对称条件下。
• 通过使用正则化迭代比例拟合程序（IPFP）的数值实验验证理论发现。
• 将最优传输解与密度泛函理论（DFT）及量子多体系统中的应用相联系。

提出的方法

• 通过下半连续、凸递减函数形式化排斥性代价，例如库仑代价 $ c(x_1,\dots,x_N) = \sum_{i<j} \frac{1}{|x_j - x_i|} $ 与谐波代价 $ c_w(x_1,\dots,x_N) = -\sum_{i<j} |x_j - x_i|^2 $。
• 应用对偶理论重构多边际Kantorovich问题，并推导必要最优性条件。
• 通过径向对称性将高维问题简化为一维，将 $ \gamma $-耦合转化为径向范数 $ r_i = \|x_i\| $ 上的测度。
• 采用熵正则化实现多边际问题的正则化版本，以确保唯一性与数值稳定性。
• 使用迭代比例拟合程序（IPFP）在离散网格上求解正则化问题，控制 $ \epsilon $-正则化参数。
• 将计算得到的耦合投影以可视化其支撑与密度结构，确认确定性映射与扩散平面行为。

实验结果

研究问题

• RQ1 在何种条件下，具有排斥性代价的多边际最优传输存在Monge型解？
• RQ2 库仑代价与谐波代价下的最优传输计划在结构与空间分布上如何不同？
• RQ3 径向对称性在降低多边际传输问题维度方面起何作用？
• RQ4 在对称的边际条件下，多边际传输中的行列式代价是否能导致唯一且确定性的耦合？
• RQ5 如正则化IPFP等数值方法如何准确逼近奇异排斥性代价函数的最小化器？

主要发现

• 对于库仑代价，在 $[-1,1]^{3}$ 上的数值模拟中，$ M_d = 1000 $ 且 $ \epsilon = 0.0005 $ 时，$ (x,y) $-投影显示出扩散平面支撑，与最大空间分离一致。
• 对于谐波代价，最优耦合表现出确定性结构，其投影集中于一个重新排列质量以最小化排斥力的映射上。
• 在 $ \mathcal{B}_{0.5} $ 上的均匀径向边际下，行列式代价的最优耦合被支持在确定性映射的图像上，经 $ M_d = 100 $ 与 $ \epsilon = 0.01 $ 验证。
• 对于指数径向边际 $ \lambda_i(r) = 4\pi r^2 e^{-4r} $，最优耦合保持确定性且数值稳定，模拟参数为 $ M_d = 300 $，$ \epsilon = 0.05 $。
• 对于三种不同的径向边际，耦合投影集中于单调重排映射上，经 $ M_d = 300 $，$ \epsilon = 0.07 $ 验证，表明在非相同密度下仍具鲁棒性。
• 正则化问题的 $ \Gamma $-极限收敛于原始最优传输问题，相对熵项捕捉了极限耦合的结构。