技术摘要：最优传输映射的稳定性与 2-Monge–Kantorovich 距离的二阶变分

1. 问题陈述

本文探讨了二次成本（$p=2$）背景下最优传输映射和势函数的定量稳定性。具体而言，它研究了当源密度 $f$ 和目标密度 $g$ 受到扰动时，最优传输映射 $T = \nabla \phi$ 及其相关的 Brenier 势函数 $\phi$ 如何变化。

虽然定性稳定性（在强假设下的测度收敛或一致收敛）已通过紧性论证得到充分确立，但作者致力于推导定量估计。这些估计将对应于不同密度的最优映射/势函数之间的距离，以密度本身之间的距离（在 $L^p$、$L^2$ 或 2-Monge–Kantorovich 距离中衡量）为界进行界定。该工作专门针对密度非退化（有界远离零）且具有Hölder 正则性的一致凸域情形。

2. 方法论

核心技术策略涉及沿数据的平滑插值线性化 Monge–Ampère 方程。方法论步骤如下：

1. 插值：给定两对密度 $(f_0, g_0)$ 和 $(f_1, g_1)$，作者构造了一族平滑密度曲线 $\{(f_t, g_t)\}_{t \in [0,1]}$。
   • 对于定理 1.1，他们利用2-Monge–Kantorovich 测地线（Benamou–Brenier 公式）。
   • 对于定理 1.2 和 1.3，他们使用 $L^p$ 空间中的线性插值。
2. 势函数微分：他们定义了对应于 $(f_t, g_t)$ 的 Brenier 势函数曲线 $\phi_t$。差值 $\phi_1 - \phi_0$ 表示为其时间导数的积分：$\phi_1 - \phi_0 = \int_0^1 \partial_t \phi_t \, dt$。
3. 线性化方程：导数 $\xi_t = \partial_t \phi_t$ 被表征为散度形式导出的线性化 Monge–Ampère 方程的唯一解：
   $$-\text{div}\left[ f_t (D^2 \phi_t)^{-1} \nabla \xi_t \right] = \left( \frac{\partial_t g_t(\nabla \phi_t)}{g_t(\nabla \phi_t)} - \frac{\partial_t f_t}{f_t} \right) f_t \quad \text{在 } \Omega \text{ 中},$$
   并满足 Neumann 边界条件。
4. 正则性与椭圆性：关键步骤是确立系数矩阵 $(D^2 \phi_t)^{-1}$ 保持一致椭圆性，且解 $\xi_t$ 沿插值路径保持足够的正则性（例如 $C^{1,\alpha}$）。这依赖于：
   • 引理 1.7：一个紧性论证，证明了如果密度是 $C^{0,\beta}$ 且非退化的，则 Brenier 势函数在边界处也是 $C^{2,\alpha}$ 的。
   • 隐函数定理：用于证明映射 $t \mapsto \phi_t$ 在适当的 Banach 空间中的光滑性（定理 1.5）。
5. 估计：通过使用椭圆正则性理论（Schauder 估计、$L^\infty$ 界）估计 $\xi_t$ 并对 $t$ 进行积分，作者推导出了最终的稳定性界。

3. 主要贡献与结果

A. 定量稳定性估计

在假设非退化、Hölder 连续密度定义在具有光滑边界的均匀凸域上的前提下，本文确立了三个主要的稳定性定理：

• 定理 1.1（$L^2$ Lipschitz 稳定性）：
  对于 Hölder 连续密度，作者证明了最优传输映射在 $L^2$ 范数下关于密度间的 2-Monge–Kantorovich 距离（$d_2$）具有Lipschitz 连续性：
  $$\|\nabla \phi_1 - \nabla \phi_0\|_{L^2} \leq C [d_2(f_0, f_1) + d_2(g_0, g_1)].$$
  这改进了以往的结果（例如 Berman [20]），以往结果通常涉及依赖于维数的指数或需要更强的正则性。作者指出，这种 Lipschitz 行为是正则非退化测度的一个特殊特征。

• 定理 1.2（无连续性假设的稳定性）：
  在去除 Hölder 连续性假设但保留 $L^2$ 可积性和非退化性的情况下，作者获得了 Brenier 势函数在 $L^2$ 范数下的Lipschitz 稳定性估计：
  $$\|\phi_0 - \phi_1\|_{L^2} \leq C (\|f_0 - f_1\|_{L^2} + \|g_0 - g_1\|_{L^2}).$$
  他们还推导出了梯度的相应估计，其指数为 $1/3$ 次幂。作者推测，如果能为线性化算子证明特定的 Sobolev 不等式，则势函数的 $L^2$ 界有可能被改进为 $L^\infty$ 界。

• 定理 1.3（一致 $C^{1,\alpha}$ 稳定性）：
  在一致度量设定下，作者证明了对于 $C^{0,\beta}$ 密度，最优映射在 $C^{1,\alpha}$ 范数下关于密度的 $L^p$ 距离（其中 $p = d/(1-\alpha)$）是稳定的：
  $$\|\phi_1 - \phi_0\|_{C^{1,\alpha}} \leq C (\|f_1 - f_0\|_{L^p} + \|g_1 - g_0\|_{L^p}).$$
  与以往的内点估计（例如 [28, 29]）不同，该结果是全局的，适用于恒等映射以外的扰动，且不受维度衰减指数的影响。作者强调，下界条件（$g_i \geq a > 0$）是根本性的；如果没有它，即使在 1D 中，也不存在 Hölder 稳定性指数 $\eta < 1$。

B. 距离的二阶变分

• 定理 1.4：作者推导出了沿密度扰动方向，平方 2-Monge–Kantorovich 距离的二阶变分的显式公式。
  $$\frac{1}{2} \frac{d^2}{dt^2}\bigg|_{t=0} [d_2(f_t, g_t)^2] = \int_\Omega \langle (D^2 \phi_0)^{-1} \nabla \xi, \nabla \xi \rangle f \, dx,$$
  其中 $\xi$ 是线性化方程的解。该公式量化了距离的线性凸性，并将之前发现于 [30] 中的形式计算（仅限于恒等映射附近的扰动）推广到一般的正则最优传输映射。

C. 技术工具：时间可微性

• 定理 1.5 与推论 1.6：工作的很大一部分致力于证明 Brenier 势函数 $\phi_t$ 关于插值参数 $t$ 是可微的。作者将导数 $\partial_t \phi_t$ 表征为散度形式线性化 Monge–Ampère 方程的弱解。这一结果对于时变最优传输问题的分析具有独立的兴趣。

4. 意义与主张

作者将其工作定位为当前活跃研究领域中最优传输定量稳定性的重要进展。

• 指数改进：本文声称在以往文献通常给出次线性指数（例如 $1/2$、$1/5$ 或依赖于维数的衰减指数）的设定中，提供了Lipschitz（指数为 1）稳定性估计。
• 全局与内点：定理 1.3 被强调为全局估计，与以往通常局限于内点估计或特定几何设定（例如余维数为 1 的曲面）的结果形成对比。
• 显式二阶变分：二阶变分公式的推导被呈现为一种工具，它显式地量化了 Monge–Kantorovich 距离的凸性，将理论理解从恒等映射扩展开来。
• 局限性：作者对其范围持谦逊态度，明确指出其结果依赖于非退化密度（有界远离零）和Hölder 正则性。他们承认，对于离散测度或消失密度，估计的 Lipschitz 性质会失效，且最优指数可能更低（推测离散目标的指数为 $1/2$）。他们并未声称解决了一般粗糙测度或奇异目标的稳定性问题。

总之，本文提供了一个严格的框架，用于在散度形式下线性化 Monge–Ampère 方程，从而在正则性和非退化性假设下获得最优传输映射和势函数的尖锐、定量稳定性界，同时也提供了传输距离二阶变分的新公式。