Measure preserving maps with bounded total variation

1. 研究背景与问题定义

核心问题：
考虑定义在开集 $\Omega \subset \mathbb{R}^d$ 上的分段仿射 Lipschitz 映射 $\phi: \Omega \to \mathbb{R}$。假设对于几乎所有的 $t > 0$，映射 $x \mapsto x + t\nabla\phi(x)$ 是单射（injective）。
在 Liu 和 Pego 的先前工作 [LP25] 中，作者提出了一个猜想：

猜想 1.1： 满足上述条件的任何映射 $\phi$ 必须在 $\Omega$ 内是局部凸的（locally convex）。

研究动机：
在最优传输理论中，凸势函数 $\phi$ 的梯度 $\nabla\phi$ 生成的映射 $T_t(x) = x + t\nabla\phi(x)$ 天然具有测度保持性质（即保持 Lebesgue 测度）。本文旨在从逆问题的角度证明：如果 $T_t$ 保持测度且满足一定的正则性条件，那么生成它的势函数 $\phi$ 必须是凸的。

主要挑战：
证明 $\nabla\phi$ 不在 $BV$（有界变差）空间中的情况目前几乎完全是一个开放问题。本文专注于$\nabla\phi \in BV_{loc}(\Omega)$ 的情形。

2. 主要结果

定理 1.2 (及定理 2.3)：
设 $\phi \in W^{1,1}(\Omega)$ 满足测度保持条件（MPC），即存在满测度集 $F \subset \Omega$，使得 $T_t(x) = x + t\nabla\phi(x)$ 在 $F$ 上对几乎所有 $t \ge 0$ 保持测度。
若假设 $\nabla\phi \in BV_{loc}(\Omega, \mathbb{R}^d)$，则 $\phi$ 在 $\Omega$ 中是局部凸的。

关键推论：
该结果证明了在 $\nabla\phi$ 具有有界变差的假设下，Liu-Pego 猜想成立。

3. 方法论与证明思路

证明过程分为两个主要步骤，结合了测度论、$BV$ 函数理论和最优传输工具。

步骤一：证明次调和性 (Subharmonicity)

• 命题 2.4 / 命题 3.3： 无论是否假设 $BV$条件，只要$T_t$ 是“扩张的”（expanding，即推前测度小于等于 Lebesgue 测度），则 $\phi$ 是次调和函数（subharmonic）。
• 数学表述： 对于所有非负测试函数 $\psi \in C^1_c(\Omega)$，满足：
  $$-\int_\Omega \nabla\phi(x) \cdot \nabla\psi(x) \, dx \ge 0 \quad \Rightarrow \quad \Delta\phi \ge 0 \quad \text{in } \mathcal{D}'(\Omega)$$
• 证明方法：
  1. 定义测度 $\rho_t = (T_t)_\sharp \mathcal{L}^d$。
  2. 利用扩张性质 $\rho_t \le \mathcal{L}^d$，构造辅助函数 $g_\psi(t) = \int \psi d\rho_t$。
  3. 通过计算 $g_\psi(t)$ 在 $t=0$ 处的导数，利用控制收敛定理和 Lipschitz 性质，导出 $\int \nabla\phi \cdot \nabla\psi \le 0$。
  4. 此外，第 4 节还提供了基于最优传输的另一种证明：利用 Brenier 的极分解定理（Polar Factorization），将 $T_t$ 分解为凸梯度映射 $\nabla f_t$ 和保测映射 $\Phi_t$ 的复合，通过分析 $t \to 0$ 时的极限行为得出 $\Delta\phi \ge 0$。

步骤二：利用 $BV$ 性质证明凸性

这是证明的核心，依赖于 $\nabla\phi \in BV_{loc}$ 的假设。

• 导数分解： 根据 $BV$函数理论，$\nabla\phi$的分布导数$D(\nabla\phi)$可分解为绝对连续部分$D_a(\nabla\phi)$和奇异部分$D_s(\nabla\phi)$。
  $$D(\nabla\phi) = D_a(\nabla\phi) + D_s(\nabla\phi)$$
• 奇异部分分析 ($D_s$)：
  • 利用 Alberti 秩 1 定理，奇异部分的极矩阵形式为 $a(x) \otimes b(x)$。
  • 由于 Hessian 矩阵必须对称，推导出 $D_s(\nabla\phi) = \lambda(x) \vec{n}(x) \otimes \vec{n}(x)$。
  • 结合步骤一得出的次调和性（$\Delta\phi \ge 0$），推导出系数 $\lambda(x) \ge 0$。
  • 结论：$D_s(\nabla\phi)$ 是一个正定的矩阵值测度。
• 绝对连续部分分析 ($D_a$)：
  • 利用 Lusin-Lipschitz 性质，构造 Lipschitz 逼近函数 $f_\Lambda$。
  • 在集合 $G_\Lambda$ 上，$D_a(\nabla\phi)$ 等同于 $Df_\Lambda$。
  • 利用测度保持条件 $T_t$ 的雅可比行列式性质：$|\det(\text{Id} + t D_a(\nabla\phi))| = 1$ 几乎处处成立。
  • 设 $\lambda_i(x)$ 为 $D_a(\nabla\phi)$ 的特征值，则 $\prod |1 + t\lambda_i(x)| = 1$ 对所有 $t \ge 0$ 成立。
  • 这迫使所有特征值 $\lambda_i(x) = 0$。
  • 结论：$D_a(\nabla\phi) = 0$。
• 最终结论：
  由于 $D_a(\nabla\phi) = 0$ 且 $D_s(\nabla\phi) \ge 0$（正定），整个分布 Hessian $D(\nabla\phi)$ 是正定的。在分布意义下，正定 Hessian 等价于函数的局部凸性。

4. 关键贡献与创新点

1. 解决了特定情形下的猜想： 首次证明了在 $\nabla\phi \in BV_{loc}$ 的假设下，Liu-Pego 猜想成立。这填补了从连续梯度（$C^0$）到一般 $BV$ 梯度之间的理论空白。
2. **去除了 $BV$假设的次调和性证明：** 证明了即使没有$BV$假设，满足测度保持条件的势函数也是次调和的（$\Delta\phi \ge 0$）。这是一个独立的强结果。
3. 多方法融合：
   • 结合了几何测度论（$BV$ 函数的奇异分解、Alberti 秩 1 定理）。
   • 结合了最优传输理论（Brenier 极分解、Wasserstein 距离、凸势函数的性质）。
   • 第 4 节提供的基于最优传输的次调和性证明，为处理非 $BV$ 情况提供了新的视角和潜在工具。
4. 揭示了绝对连续部分的消失： 通过测度保持条件对雅可比行列式的约束，严格证明了在 $BV$ 框架下，梯度的绝对连续部分必须为零，从而将问题简化为对奇异部分的分析。

5. 局限性与未来方向

• **$BV$假设的必要性：** 目前证明依赖于$\nabla\phi \in BV_{loc}$。如果$\nabla\phi$不属于$BV$（例如具有更奇异的性质），证明中的$D_a$和$D_s$ 分解及相应的分析工具不再直接适用。
• 开放问题： 论文明确指出，证明 $\nabla\phi \notin BV$ 时的猜想仍然是一个几乎完全开放的问题。
• 未来工作： 需要开发新的分析工具来处理非 $BV$ 梯度的测度保持映射，或者寻找反例。

6. 总结与意义

这篇论文在最优传输和流体动力学的交叉领域做出了重要贡献。它通过严谨的 $BV$ 分析，确认了“刚性破缺势流”（Rigidly breaking potential flows）在特定正则性条件下必然对应于凸势函数。这不仅验证了物理直觉（凸势对应稳定的流体运动），也为理解测度保持映射的几何结构提供了深刻的数学依据。其证明技巧，特别是将测度保持条件转化为 Hessian 矩阵特征值的约束，具有广泛的推广价值。