2025-11-22T16:07:16.564646

Unstable optimal transport maps

Letrouit
The stability of optimal transport maps with respect to perturbations of the marginals is a question of interest for several reasons, ranging from the justification of the linearized optimal transport framework to numerical analysis and statistics. Under various assumptions on the source measure, it is known that optimal transport maps are stable with respect to variations of the target measure. In this note, we focus on the mechanisms that can, on the contrary, lead to instability. We identify two of them, which we illustrate through examples of absolutely continuous source measures $ρ$ in $\mathbb{R}^d$ for which optimal transport maps are less stable, or even very unstable. We first show that instability may arise from the unboundedness of the density: we exhibit a source density on the unit ball of $\mathbb{R}^d$ which blows up superpolynomially at two points of the boundary and for which optimal transport maps are highly unstable. Then we prove that even for uniform densities on bounded open sets, optimal transport maps can be rather unstable close enough to configurations where uniqueness of optimal plans is lost.
academic

Unstable optimal transport maps

基本信息

  • 论文ID: 2510.13265
  • 标题: Unstable optimal transport maps
  • 作者: Cyril Letrouit (Université Paris-Saclay, CNRS)
  • 分类: math.OC (Mathematical Optimization and Control)
  • 发表时间: 2025年10月15日
  • 论文链接: https://arxiv.org/abs/2510.13265

摘要

最优传输映射相对于边际分布扰动的稳定性是一个重要问题，涉及线性化最优传输框架的合理性、数值分析和统计学等多个方面。在源测度的各种假设下，已知最优传输映射相对于目标测度的变化是稳定的。本文关注相反的机制——导致不稳定性的因素。作者识别了两种机制，并通过Rd\mathbb{R}^d中绝对连续源测度ρ\rho的例子说明了最优传输映射的不稳定性甚至高度不稳定性。首先证明了密度的无界性可能导致不稳定：构造了Rd\mathbb{R}^d单位球上的源密度，该密度在边界的两点处超多项式爆炸，对应的最优传输映射高度不稳定。然后证明即使对于有界开集上的均匀密度，当接近最优计划唯一性丧失的配置时，最优传输映射也可能相当不稳定。

研究背景与动机

问题描述

最优传输理论中的一个核心问题是研究最优传输映射TμT_\mu（从固定源测度ρ\rho到目标测度μ\mu）相对于目标测度扰动的稳定性。具体来说，研究是否存在常数C,α>0C, \alpha > 0和p≥1p \geq 1，使得对于紧集Y⊂RdY \subset \mathbb{R}^d上的任意概率测度μ,ν\mu, \nu，有： ∥∥Tμ−Tν∥∥L2(ρ)≤CWp(μ,ν)α\|\|T_\mu - T_\nu\|\|_{L^2(\rho)} \leq C W_p(\mu, \nu)^\alpha

重要性

这个问题的重要性体现在多个方面：

  1. 线性化最优传输框架的合理性：映射μ↦Tμ\mu \mapsto T_\mu提供了从Wasserstein空间(P(Y),W2)(P(Y), W_2)到Hilbert空间L2(ρ,Rd)L^2(\rho, \mathbb{R}^d)的嵌入
  2. 数值分析：当目标测度μ\mu只能通过近似μ^\hat{\mu}获得时，需要保证计算得到的Tμ^T_{\hat{\mu}}接近真实的TμT_\mu
  3. 统计估计：在统计最优传输中，需要从样本估计最优传输映射

现有研究局限性

虽然在源测度ρ\rho的各种假设下已经建立了量化稳定性不等式，但此前从未证明最优传输映射可能是不稳定的。现有结果主要关注稳定性的正面结果，缺乏对不稳定性机制的深入理解。

核心贡献

  1. 首次构造高度不稳定的最优传输映射：证明存在绝对连续源测度，使得任何形式的量化稳定性不等式都不成立
  2. 识别两种不稳定性机制：
    • 密度的超多项式爆炸
    • 接近最优计划非唯一的配置
  3. 提供具体的反例构造：通过精确的数学构造展示不稳定性现象
  4. 揭示映射μ↦Tμ\mu \mapsto T_\mu的正则性限制：证明在某些情况下该映射不是12\frac{1}{2}-Hölder连续的

方法详解

任务定义

给定ρ,μ∈P2(Rd)\rho, \mu \in P_2(\mathbb{R}^d)，其中ρ\rho绝对连续，Brenier定理保证存在唯一的最优传输映射Tμ∈L2(ρ)T_\mu \in L^2(\rho)从ρ\rho到μ\mu。研究TμT_\mu相对于μ\mu扰动的稳定性。

定理1.1的构造方法

源测度构造

设A=(1,0,…,0)A = (1, 0, \ldots, 0)，A′=(−1,0,…,0)A' = (-1, 0, \ldots, 0)，E={A,A′}E = \{A, A'\}。定义函数： f(r)=r−dmin⁡(1,(log⁡r)−2)f(r) = r^{-d} \min(1, (\log r)^{-2})

构造密度： ρ(x)=c0f(dist(x,E))\rho(x) = c_0 f(\text{dist}(x, E)) 其中c0>0c_0 > 0使得ρ\rho是单位球BRd(0,1)B_{\mathbb{R}^d}(0,1)上的概率密度。

目标测度序列

考虑Bθ=(Rsin⁡θ,Rcos⁡θ,0,…,0)B_\theta = (R\sin\theta, R\cos\theta, 0, \ldots, 0)和Bθ′=(−Rsin⁡θ,−Rcos⁡θ,0,…,0)B'_\theta = (-R\sin\theta, -R\cos\theta, 0, \ldots, 0)，定义： μθ=12(δBθ+δBθ′)\mu_\theta = \frac{1}{2}(\delta_{B_\theta} + \delta_{B'_\theta})

关键估计

通过几何分析证明： ∥∥Tμθ−Tμ0∥∥L2(ρ)2≥c0c1σd−1R2∣log⁡(θ/4)∣\|\|T_{\mu_\theta} - T_{\mu_0}\|\|^2_{L^2(\rho)} \geq \frac{c_0 c_1 \sigma_{d-1} R^2}{|\log(\theta/4)|}

而Wp(μ0,μθ)∼Rθ/2W_p(\mu_0, \mu_\theta) \sim R\theta/2，因此稳定性不等式对任何α>0\alpha > 0都失效。

定理1.2的构造方法

支撑集构造

构造无穷个矩形平行六面体对的并集： X=⋃i=1+∞SiX = \bigcup_{i=1}^{+\infty} S_i 其中Si=T+(Ai+,ℓi,ri)∪T−(Ai−,ℓi,ri)S_i = T^+(A_i^+, \ell_i, r_i) \cup T^-(A_i^-, \ell_i, r_i)。

参数选择

选择序列满足：

  • ri=c02−ir_i = c_0 2^{-i}（超多项式衰减）
  • ℓi=wi=c0c1i−2\ell_i = w_i = c_0 c_1 i^{-2}
  • 细胞间距离远大于细胞内距离

不稳定性证明

通过分析特定扰动下的最优传输映射变化，证明： ∥∥Tμ−Tνi∥∥L2(ρ)2Wp(μ,νi)2α≥Cwi3−2αpri1−2αp−2α\frac{\|\|T_\mu - T_{\nu_i}\|\|^2_{L^2(\rho)}}{W_p(\mu, \nu_i)^{2\alpha}} \geq C w_i^{3-\frac{2\alpha}{p}} r_i^{1-\frac{2\alpha}{p}-2\alpha}

当α>p2(p+1)\alpha > \frac{p}{2(p+1)}时，右边趋于无穷。

实验设置

理论验证

本文为纯理论工作，通过严格的数学证明验证结果：

  1. 几何分析：详细分析最优传输映射的几何结构
  2. 测度估计：精确计算各区域的测度
  3. 距离计算：准确估计Wasserstein距离和L2L^2距离

构造验证

通过以下步骤验证构造的正确性：

  1. 证明构造的密度确实是概率密度
  2. 验证最优传输映射的显式表达式
  3. 确认几何配置满足所需的分离条件

实验结果

主要结果

定理1.1结果

对于构造的源密度ρ\rho（在单位球边界两点处超多项式爆炸），证明了对任何球Y=BRd(0,R)Y = B_{\mathbb{R}^d}(0,R)、任何C,α>0C, \alpha > 0和p≥1p \geq 1，不等式： ∀μ,ν∈P(Y),∥∥Tμ−Tν∥∥L2(ρ)≤CWp(μ,ν)α\forall \mu, \nu \in P(Y), \|\|T_\mu - T_\nu\|\|_{L^2(\rho)} \leq C W_p(\mu, \nu)^\alpha 都不成立。

定理1.2结果

对于有界开集上的均匀概率测度ρ\rho，证明了对任何C>0C > 0、p≥1p \geq 1和α>p2(p+1)\alpha > \frac{p}{2(p+1)}（特别地，α=12\alpha = \frac{1}{2}），稳定性不等式失效。

技术发现

  1. 密度爆炸的临界性：多项式爆炸r−d+δr^{-d+\delta}仍保持稳定性，但超多项式爆炸导致完全不稳定
  2. Hölder正则性限制：映射μ↦Tμ\mu \mapsto T_\mu在某些情况下不是12\frac{1}{2}-Hölder连续，甚至不是13\frac{1}{3}-Hölder连续
  3. 非唯一性的预兆：最优计划唯一性丧失前会出现稳定性恶化

相关工作

稳定性正面结果

  • Gigli (2011)：映射μ↦Tμ\mu \mapsto T_\mu的12\frac{1}{2}-Hölder连续性理论上界
  • Letrouit & Mérigot (2024)：John域上有界密度的稳定性不等式
  • 其他工作：对数凹密度、多项式爆炸密度等的稳定性结果

不稳定性研究

此前仅有Kantorovich势函数的不稳定性例子，本文首次构造了不稳定的传输映射。

结论与讨论

主要结论

  1. 最优传输映射可能高度不稳定，任何量化稳定性不等式都可能失效
  2. 密度的超多项式爆炸是不稳定性的充分条件
  3. 即使对于均匀密度，接近非唯一配置时也会出现不稳定性

局限性

  1. 定理1.1需要无界密度，实际应用中可能不现实
  2. 定理1.2的构造需要无穷多个连通分量，几何结构复杂
  3. 未能确定有界密度是否能达到定理1.1的强不稳定性结论

未来方向

  1. 研究有界密度情况下的最强不稳定性
  2. 开发针对不稳定源测度的统计估计方法
  3. 探索不稳定性与几何结构的深层联系

深度评价

优点

  1. 理论突破性：首次证明最优传输映射的不稳定性，填补了理论空白
  2. 构造精妙：通过巧妙的几何构造和参数选择实现了理论目标
  3. 分析严谨：数学证明完整，技术细节处理得当
  4. 洞察深刻：揭示了稳定性与唯一性、密度正则性的内在联系

不足

  1. 实用性限制：构造的例子在实际应用中可能不常见
  2. 开放问题：某些自然问题（如有界密度的强不稳定性）仍未解决
  3. 数值验证缺失：作为理论工作，缺乏数值实验验证

影响力

  1. 理论价值：为最优传输理论提供了重要的反例和洞察
  2. 应用警示：提醒实践者注意某些配置下的不稳定性风险
  3. 研究启发：为后续研究提供了新的方向和工具

适用场景

  1. 最优传输理论的深入研究
  2. 线性化最优传输框架的理论分析
  3. 统计最优传输中的稳定性分析
  4. 数值最优传输算法的稳定性评估

参考文献

论文引用了最优传输领域的重要文献，包括Brenier的奠基性工作、Gigli的正则性结果、以及近期的量化稳定性研究等，为研究提供了坚实的理论基础。