精读笔记
Problem Setting
这篇论文实际处理的是 QOT 解映射的定量稳定性，而不是重新分析 QOT 的存在性或 sparsity。对象是四元组 (P,Q,c,ε) 变化时，dual direct sum h=f⊕g、optimal density ζ=(h−c)+/ε、optimal coupling π 以及 support Σ={h−c>0} 如何变化。
真正困难点有两个。第一，QOT 的自然 dual 几何在 H⊂L2(P⊗Q) 中，但 P,Q 改变后空间本身变了；fixed-space perturbation theory 不能直接用于 Wasserstein perturbation。第二，support 是一个水平集问题，L2 稳定不足以控制 Hausdorff 距离；必须得到 L∞ 稳定，并且还要排除 slack 在边界附近 flat 的情况。
以前方法卡在这里：EOT 可以利用 Schrödinger system 的 smooth implicit map 和 strong concavity-like 结构，但 QOT dual 有正部非光滑和非强凹；divergence-regularized OT 的 shadow/data-processing 路线能给 coupling 或 value 的稳定性，却看不到 potentials，更看不到 support。关键矛盾是：QOT 的 sparsity 是它相对 EOT 的主要优势，但 sparsity 的几何对象 support 又是最不连续、最需要 pointwise control 的对象。
Motivation
已有路线缺的是对 sparse structure 本身的稳定性控制。EOT 的 full support 使 support 问题退化；classical OT 的 stability 多是 weak convergence 或 map stability，机制不同；一般 divergence regularization 的稳定性结果通常停留在 primal metric，不提供 boundary-level 信息。
作者的核心观察是 QOT 的 primal 对象完全由 dual slack σ=h−c 决定：density 是 σ+ 的线性函数，support 是 σ 的正集。因此，真正该控制的是 h 的 L∞ 扰动，而不是只控制 coupling 的 Wasserstein 距离。
为什么会想到这个方向：近期 [20] 已经给出 fixed marginals 下 QOT dual 的 local PL/error bound，这相当于提供了局部可逆性的弱替代。本文要做的是把这个 fixed-space 局部几何转译到 moving marginals，并进一步升级到 L∞，最终服务于 support stability。缺口非常明确：从 L2 error bound 到 varying-marginal L∞ potential stability，再到 support Hausdorff stability。
Core Idea
核心思想是把 QOT 稳定性组织成一条 dual-to-primal-to-support 的链：先控制 h=f⊕g，再由 ζ=(h−c)+/ε 控制 density/coupling，最后由 Σ={σ>0} 控制 support。这个视角把 sparse support 的不连续性问题转化为 slack 零水平集的稳定性问题。
理论上它能成立，是因为 QOT 的 FOC 虽然不像 EOT Schrödinger system 那样显式给出 log-integral 自映射，但仍然隐含了局部单调性：对固定 x，a↦∫(a+g(y)−c(x,y))+dQ(y) 在最优 f(x) 附近有正斜率下界，只要在某个 y-neighborhood 上 slack 保持正。这个斜率下界把 residual 小转化为 f 点态小；再配合 L2 已经控制了 g 的平均误差，就能闭合成 L∞ 控制。
和 prior 的本质区别在于，它不是从 primal coupling metric 入手，也不是依赖 EOT 的 smooth exponential structure，而是用 QOT 自身的正部积分方程恢复 pointwise stability。这是对 QOT 几何的一次更细粒度使用：不是只说 sparse，而是证明 sparse support 的边界在非退化条件下可稳定追踪。
Method
1. L2 稳定性：解决 moving-marginal 下 fixed PL/error bound 不能直接用的问题。作者把 h′ 作为 reference 空间 H⊂L2(P⊗Q) 中的测试点，估计 DΦ(h′)。marginal 变化项通过 potentials/slack 的 Lipschitz 性和 W1 duality 控制，cost 与 ε 的变化直接进入 residual。核心变化是把“空间变了”的问题改写为 reference objective 上的 gradient residual 问题。
2. locality 自动闭合：error bound 本身需要 L∞ locality。作者用 potentials 的 Lipschitz 性和 reference marginals 的 ball mass lower bound，把若 L∞ 偏差达到 δ 则 L2 偏差至少有固定质量的反证结构建立起来。也就是说，局部 L2 小在 Lipschitz + mass lower bound 下能阻止 L∞ 爆点。这一步是把 local condition 从假设变成可验证 small perturbation condition。
3. L∞ 升级：这是论文的关键机制。通过 FOC 构造 Fx(a)=∫(a+g−c)+dQ，并证明在 f(x) 附近 Fx 有统一斜率下界；再利用 h−h′ 的 L2 控制、cost/measure/ε 扰动控制 Fx(f′(x))−ε。斜率下界迫使 |f−f′| 点态小，g 同理。这个步骤解决了 support 需要 pointwise slack control 的问题。
4. primal 与 support：density 稳定性直接来自正部函数 1-Lipschitz 和 ε 缩放；coupling 的 TV/W1 稳定性还要加上 product reference measure P⊗Q 本身的变化。support 稳定性则额外引入非退化条件 σ(z)≤−a dist(z,Σ) outside support；这把 L∞ slack error 转成 Hausdorff boundary displacement。二次 cost 下通过 σ 对每个坐标的凹性得到 a≈ε/diam(ΩP)。
Key Insight / Why It Works
最核心贡献是 L2→L∞ 的机制，而不是后面的 density/coupling 推论。QOT 的正部结构看起来比 EOT 的指数结构更不光滑，但它也带来一个局部线性斜率：只要在某个区域 slack 为正，正部积分对 potential shift 的响应就是该区域的质量。因此，ball mass lower bound 给出统一斜率，FOC 给出 self-bound。这是整篇论文真正有效的原因。
L2 stability 本身主要是已有 [20] error bound 的迁移和包装，但 moving-marginal residual estimate 是必要桥梁。它的技术价值在于避免要求 P′,Q′ 相对某个固定 reference 绝对连续，从而能用 W1 比较 marginals。不过这一步的概念增量小于 L∞ self-bound。
support stability 的核心不是“potential 稳定推出 support 稳定”，这句话通常是假的。真正的 insight 是：必须把 support stability 看成 level-set stability，并显式加入边界非退化。Example 6.2 很重要，因为它说明没有 linear detachment 时，即使 potentials uniform 收敛、marginal supports 不动，support 也可以不收敛。这排除了很多过度乐观的说法。
哪些部分可能只是辅助：TV/W1 coupling stability 基本是由 density bound + reference measure perturbation 机械推出；它们验证链条完整，但不是最原创部分。显式常数也更多是 theoretical bookkeeping，实际 tightness 增益来源不清。真正可迁移的是“PL/error bound + self-bound + level-set nondegeneracy”的组合框架。
Relation To Prior Work
最接近的两条线是 divergence-regularized OT stability 和 EOT potential stability。前者如 Bayraktar–Eckstein–Zhang 给 coupling 的 Hölder stability，但技术路线是 shadow/data processing，天然偏 primal，不能解析 potentials 和 support。本文本质上是 dual-geometric 路线。后者如 Carlier–Laborde、IPFP/Sinkhorn stability 使用 Schrödinger system 或 strong smooth structure；本文明确指出 QOT 无法照搬，因为 dual 非强凹且正部结构非光滑。
和 QOT sparsity 文献相比，已有工作主要描述固定 marginals、小 ε 下 support 的几何或 sparsity bounds；本文问的是这些 support 在数据扰动下是否稳定。它不是重新证明 sparse，而是给 sparse object 加上 perturbation theory。
看似新但其实是已有思想重组的部分：local error bound/PL → stability 是标准 perturbation 分析思想；density/coupling 由 potential 推出也是 QOT duality 的直接后果。实质创新在两个地方：一是 moving L2 space 下的 local Lipschitz estimate；二是从 QOT FOC 中抽取 L∞ self-bound，并进一步连接到 support Hausdorff stability。它属于 regularized OT 的 deterministic sensitivity analysis，但比以往 value/coupling stability 更接近 free-boundary/level-set stability。
Dataset / Evaluation
这篇是理论论文，没有 dataset、实验 benchmark 或真实系统 evaluation。所谓 evaluation 是定理覆盖范围和反例充分性。
核心 claim 的验证方式比较干净：先给一般 Lipschitz cost 下的 potential/coupling 稳定，再在 support 部分明确指出需要非退化条件；随后用二次 cost 证明该条件成立，并用反例证明该条件不可省。这个证据结构比单纯给 sufficient condition 更有说服力。
但从应用角度看，定理的可用性没有被数值验证。常数依赖维度、ε、ball mass、density lower bound、support diameter 等，实际是否 informative 文中未充分说明。尤其在 weak regularization ε→0 或高维稀疏数据场景，常数可能退化到不能指导实践。因此 evaluation 支撑了数学 claim，但没有支撑“实际稳定性强”或“算法层面可用”的 claim。
Limitation
第一，结果是局部的，而且 locality threshold 依赖非常强的 regularity constants。对真实 empirical/discrete marginals，Assumption 2.3 中 P 有密度且在凸紧支撑上上下有界并不自然；虽然 Q 可以更一般，但整体仍偏连续几何设定。
第二，小 ε regime 是 QOT 最关心的 sparse/OT-like regime，但常数随 ε 可能严重恶化。support stability 的斜率 a≈ε/diam(ΩP)，意味着 ε 越小，边界对 slack error 越敏感；这和实际 sparse structure 可能更脆弱是一致的，但也限制了结果强度。
第三，support 稳定性必须额外输入 marginal supports 的 Hausdorff displacement。W1 小并不控制 support 远端少量质量的移动，因此这不是缺陷而是事实；但它说明该结果不能仅作为 distributional stability 使用，需要几何 support control。
第四，一般 cost 下非退化条件是核心瓶颈。二次 cost 可以靠 concavity 解决，但其他 cost 的可验证条件文中未充分说明。换言之，support stability 的 generality 实际上由 slack boundary geometry 决定，而不是由 QOT regularization 自动保证。
第五，显式常数很可能非常保守；文中没有 tightness 分析。这里的增益不是 scaling 或 data coverage，而是更好的 inductive bias/latent geometry extraction；但实际数值上是否有可观稳定半径，仍不清楚。
Takeaway
- 1. QOT 的 sparse support 不能只靠 coupling metric 研究；正确对象是 dual slack σ=h−c 及其零水平集。
- 2. 对 QOT 来说，L∞ potential stability 是连接算法/统计扰动和几何 support 稳定性的关键中间层；L2/weak stability 不够。
- 3. support stability 本质是 level-set stability，必须有非退化/linear detachment。
- 二次 cost 的坐标凹性给了一个非常有用的可迁移模板。
一句话总结
这篇论文把 QOT 稳定性从 value/coupling 层推进到 dual potential 的 L∞ 控制和 sparse support 的 Hausdorff 稳定，是 QOT 从“有稀疏性”走向“稀疏结构可扰动分析”的一篇理论奠基工作。