精读笔记

Problem Setting

论文实际解决的是 regularized OT 在 ε→0 时相对于 unregularized quadratic OT 的 sharp bias asymptotics。这里的困难不在于 ROT/EOT 的 Γ-convergence 或 vanishing-regularization limit，而在于首个非平凡 correction term 的尺度、常数和适用条件。

关键矛盾是：局部上，最优 coupling 应该集中在 Brenier map 图附近，并在横向方向呈现一个可显式求解的 profile；但全局上，任何 admissible coupling 必须严格满足两个边缘约束。前者是局部分析问题，后者是全局可行性问题。已有 sharp 结果卡在后者：EOT 可用 Gaussian/Markov-type correction，QOT 可用 Barenblatt 加复杂边界修补，但这些修补高度依赖具体 regularizer 和几何条件。

因此本文真正要解决的是：能否用一个统一的机制，把局部最优 profile 转化为真实 coupling，同时不破坏 sharp 常数。

Motivation

已有路线的问题是把两件本应分离的事情绑在一起：一边计算局部最优 profile，一边强行让该 profile 本身满足边缘约束。这样会导致上界构造变得非常脆弱，尤其在非紧支撑、密度退化、边界不规则或 p≠2 时很难推广。

作者的核心观察是：sharp 常数其实只由 Brenier 图附近的局部 profile 决定；剩余边缘质量只要足够小、补得足够便宜，就不应影响 leading order。换句话说，global marginal correction 不需要继承局部最优结构，它只需要在正确尺度下 negligible。

这个观察把问题从“构造一个全局近似最优 coupling”改写成“构造一个局部最优 subcoupling + 一个低成本 remainder coupling”。这是本文的关键缺口填补。

Core Idea

核心思想是 profile/completion separation。先在源-源坐标中看 coupling π∈Π(μ,μ)，再通过 (Id,∇φ) 推到 Π(μ,ν)。Brenier 图被拉回为对角线，运输代价相对 OT 的增量正好由 Bregman divergence 控制。由于 A2 给出 Hessian 一致椭圆，Bregman divergence 在对角线附近二阶近似为局部二次型。

局部 profile 来自对偶的一维凸共轭结构：entropy 给 Gaussian，L^p 给 compactly supported Barenblatt。真正新的是，作者不要求这些 profile 拼起来直接成为 full coupling，而是只构造 dominated subcoupling，其边缘小于 μ；剩下的 μ-缺口用 quantization/block coupling 补齐。这样 sharp 常数完全由 leading subcoupling 贡献，remainder 只作为可行性修补。

与 prior 的本质区别不在 Gaussian 或 Barenblatt 本身，这些局部形状在 EOT/QOT 中已有；区别在于边缘约束的组织方式从 problem-specific correction 变成 regularizer-aware quantization completion。

Method

下界机制：使用 regularized OT 的弱对偶。以 Kantorovich potential 为基准，只在源端加入 ε 尺度或 ε^{2/D} 尺度的 perturbation。对偶 penalty 在 Brenier 图附近局部化，Bregman divergence 用 Hessian 二阶展开，积分极限化为 Gaussian/Barenblatt profile 的显式积分。这一步解决的是 sharp 常数的不可低于性。

局部 profile 机制：EOT 的 Ξ′ 是指数，因此归一化条件产生 Gaussian kernel；L^p 的 Ξ′ 是 positive-part polynomial，因此产生 Barenblatt kernel。尺度由“局部运输扩散成本”和“正则化密度爆炸成本”平衡得到：L^p 中半径约 ε^{1/(d(p-1)+2)}，偏差约 ε^{2/(d(p-1)+2)}。

上界机制：用简单函数从下逼近 f_μ，在每个 cell 内冻结 Hessian 和密度，放入局部 profile，得到 symmetric subcoupling。因为 profile 被截在 cell×cell 内，边缘自动 dominated by μ。剩余边缘质量 λ_rem 用 quantization coupling 补齐。

regime split：p≤2 与 EOT 的 remainder 可通过 concavity/Jensen 控制，估计只依赖总剩余质量和 cell 数，因此支持非紧支撑和弱密度条件。p>2 时 r^{p-1} convex，极小 cell 会放大 L^p penalty，因此必须构造 comparable-mass cells；这解释了额外的 Lipschitz boundary 和密度上下界假设。

Key Insight / Why It Works

最重要的 insight 是：sharp asymptotic constant 是局部对象，不是全局 coupling correction 的对象。只要 leading mass 在 Brenier 图附近以正确 profile 分布，常数就确定了；marginal feasibility 可以由一个 asymptotically negligible remainder 单独处理。这一点把问题从“全局精细构造”降维成“局部 profile + 粗粒度 quantization”。

方法成立的根本原因有三层。第一，quadratic cost 与 Brenier map 让相对 OT 的 excess cost 精确等于 Bregman divergence 积分；这给了局部二次化的入口。第二，convex duality 把 profile 选择变成点态凸共轭计算，因此 Gaussian/Barenblatt 不是猜出来的，而是由 regularizer 的 Legendre geometry 强制产生。第三，quantization 的误差尺度恰好能和 regularization 尺度平衡，并且可被压到 leading term 以下。

本文最核心贡献是 upper bound 的 quantization completion，而不是下界或局部 profile 计算。下界基本沿经典 dual ansatz；Gaussian/Barenblatt 也不是全新对象。真正新增的是证明这些局部对象可以在弱条件下拼成 admissible coupling，同时不污染 leading constant。

这不是 scaling，也不是 data coverage，也不是 representation alignment。它是一个证明架构上的 decomposition：把 exact marginal constraint 从局部最优 profile 中剥离出来。可以迁移的 insight 是：在带硬约束的变分渐近问题中，leading-order optimizer 不一定需要自身满足全部硬约束；可以先构造 dominated leading object，再用低阶 completion 恢复可行性。

Relation To Prior Work

最接近的是 Pal 2024 的 EOT asymptotics、Garriz-Molina/González-Sanz/Mordant 2024 的 QOT asymptotics，以及 Eckstein-Nutz 2024 的 quantization-rate bounds。本文处在 Schrödinger/EOT asymptotics、Tsallis/L^p sparse regularization、regularized OT sample/algorithmic theory 的交叉位置。

与 Pal 的差异：Pal 的 EOT 上界依赖 Gaussian perturbation 和 Markov construction，适合 entropy 的指数结构；本文把 Gaussian 只作为 local leading block，然后用 quantization 补边缘，因此支持非紧支撑和更弱密度条件。

与 QOT asymptotics 的差异：此前 p=2 的 Barenblatt profile 已出现，但上界需要 transport correction、boundary frame、stained-glass couplings 等专门结构。本文把这些 problem-specific patching 替换为统一的 dominated subcoupling + completion。

与 Eckstein-Nutz 的关系更微妙：quantization 用于 regularized OT convergence rate 已有，但那些结果是 rate-level、非 sharp constant；本文把 quantization 嵌入 sharp asymptotic proof，并让它只承担 feasibility completion，而不是决定主常数。

看似新的 Gaussian-to-Barenblatt interpolation，本质上来自 Tsallis/KL convex conjugate 的连续过渡；实质创新在于这个过渡被放进同一套上下界证明，而不是只作为形式类比。

Dataset / Evaluation

这是一篇纯数学分析论文，没有 dataset、实验、benchmark 或真实系统评估。评价标准是定理覆盖范围、假设强度、上下界是否匹配、常数是否显式、是否统一已有特殊情形。

从理论 evidence 看，论文的主 claim 得到了直接支持：对 EOT 给出二阶 sharp asymptotics；对所有 p>1 的 L^p-ROT 给出首阶 sharp asymptotics；并解释 p→1 时 Tsallis regularization 与 EOT 常数的一致性。

但 evaluation 的局限也很清楚：它验证的是 asymptotic ε→0 的解析极限，不直接说明有限 ε 下的数值误差、算法稳定性、sample complexity 常数或实际 sparse coupling 质量。若把这篇论文用于 computational OT 的实践论证，需要额外数值或非渐近结果；本文本身不提供这部分证据。

Limitation

A2 是主要结构性限制：要求 Brenier potential 全局 C^2 且 Hessian 上下界一致。这相当于强凸强光滑的全局 Monge map，排除了许多实际分布、非光滑支撑、退化密度、低维流形支撑和多解/非唯一结构。论文说假设 mild，更多是相对 prior 的边界/密度条件而言；从 OT regularity 角度看并不弱。

p>2 的处理暴露了方法上限。由于 L^p penalty 对小 mass cell 极其敏感，remainder completion 必须控制 cell mass comparability。因此 finite-constant 情形仍需 bounded Lipschitz support 与密度上下界。也就是说，quantization 并没有完全消除全局几何条件，而是把它们集中到了 p>2 的 remainder 控制里。

EOT 的 K_EOT 可为 -∞，L^p 的 K 可为 ∞。论文允许这些 extended limits，这是理论上完整的，但在应用解释上意义有限：当常数发散时，sharp expansion 更多是在说明 chosen normalization 不够捕捉有限偏差，而不是给出可用近似。

p→1 的一致性只是 fixed ε 下常数层面的展开。文中明确没有从 p-ROT 严格推出 EOT asymptotics，因为双极限需要额外 uniformity。这里的桥接是很有启发性的，但不是完整的 unified limiting theorem。

此外，所有结果围绕 quadratic cost 和 Brenier map 的 Bregman structure。若 cost 非二次、map 正则性差或 unregularized OT optimizer 不唯一，核心机制是否仍可保 sharp 常数，文中未充分说明。

Takeaway

• 第一，regularized OT 的 small-ε bias 本质上由 Brenier 图附近的局部 profile 决定；global marginal constraint 不是主常数来源，而是 feasibility noise，只要能被低阶处理。
• 第二，entropy 与 L^p regularization 的差异可以统一理解为 convex conjugate 诱导的局部 profile 差异：Gaussian 对应指数尾，Barenblatt 对应 compact support polynomial profile。
• p 控制的不只是 sparsity，也控制 asymptotic scale。
• 第三，quantization completion 是可迁移的证明模式：在其他正则化、其他 hard-constrained variational problems 中，可以尝试先构造 dominated local optimizer，再用 block/quantization 补齐约束。

一句话总结

这篇论文把 regularized OT 的 sharp small-regularization 渐近从 EOT/QOT 特例推进到统一的 L^p/Tsallis 图景，真正贡献是用“局部 profile + quantization completion”重构了满足边缘约束的上界证明。