精读笔记

Problem Setting

【Projected McKean--Vlasov Dynamics for Entropic Weak Optimal Transport】（arXiv preprint / 2026）

这篇论文实际处理的是：给定 µ、ν，最小化 J(π)=∫ c(x,π_x)µ(dx)+εH(π|µ⊗ν)，其中 c 对条件分布 π_x 非线性。目标不是提出一个新的 weak OT 静态定理，而是为这个熵正则问题构造一条连续时间、保两个边际、并收敛到 minimizer 的动力学。

关键困难有两个。第一，weak cost 对 conditional kernel x↦π_x 敏感，而 ordinary W_2 收敛不控制 disintegration；所以 classical Wasserstein gradient flow 的状态空间错了。第二，边际约束 Π(µ,ν) 要求 X 边际和 Y 边际都固定，但 weak force 本身依赖 L(Y|X)，投影修正又必须通过 L(X|Y) 表达；这使漂移成为 forward/backward conditioning 交织的 McKean--Vlasov 对象。

以前 projected Langevin for entropic OT 可以依赖 joint-density regularity 和 classical cost 线性结构；这里 X 不扩散，只有 Y 扩散，导致无法直接利用全空间椭圆正则来让 conditioning 连续。核心矛盾是：优化目标需要保留 π_x 的非线性信息，但动力学又必须持续保持 Y 边际 ν 不变。

Motivation

已有路线缺的是一个和 weak cost 匹配的几何。weak OT 的基本变量不是普通 joint law，而是 joint law 加上其 conditional information structure。若只在 P_2(X×R^d) 上看，两个 joint laws 可以 Wasserstein 很近但 conditional kernels 差异很大，因此 J 不连续，梯度流更谈不上稳定。

作者的观察是：adapted Wasserstein / nested distance 正好把 conditional law 纳入拓扑。于是 projected Langevin 的思想可以迁移，但必须改写成 adapted geometry 下的投影梯度流。这个方向自然来自两个缺口的交叉：一边是 weak OT 静态理论已有 entropy minimizer 和 dual representation；另一边是 entropic OT 已有 projected Langevin dynamics。缺的是把二者连接起来的动态机制。

Core Idea

论文的核心思想可以压缩成一句：把 weak OT 的条件核非线性看成 adapted Wasserstein 空间中的一阶结构，然后把梯度投影到固定边际耦合集合的切空间上。

未投影的 y-force 是 \hat c(x,y,π_x)，它知道当前 x 对应的 conditional distribution；这保留了 weak transport 的本质。为了不破坏第二边际 ν，必须减去在同一 y 位置上的平均 force，即 \tilde c_π(y)=E[\hat c(X,Y,π_X)|Y=y]。这个 conditional centering 是全篇最重要的结构：它不是普通 mean-field correction，而是对固定 ν 约束的局部拉格朗日投影。

和 prior 的本质区别在于：classical entropic OT 的 force 只依赖 (x,y)，投影项是线性对象；这里 force 先依赖 L(Y|X)，再按 L(X|Y) 做投影。因此信息流变成双向条件化：forward kernel 决定 weak force，backward kernel 决定 marginal-preserving projection。这是新的 inductive bias，也是理论复杂度的来源。

Method

1. Adapted gradient：解决 ordinary Wasserstein gradient 不适配 conditional cost 的问题。形式上，测试函数 φ 的 adapted gradient 为 (E[∂_xφ|X=x], ∂_yφ)。第一阶段运动只能依赖 X，第二阶段运动可依赖 (X,Y)。这一步把 weak OT 的信息结构编码进切空间。

2. Projection onto Π(µ,ν)：解决梯度下降会破坏边际的问题。作者采用 P^π(u_1,u_2)=(0,u_2-E[u_2|Y])。它冻结 X，并从 Y 速度中去掉只依赖 Y 的分量，从而保持 Y 边际。机制上这是 conditional centering，不是简单归一化。

3. Projected McKean--Vlasov SDE：由投影梯度得到 dY_t=-(\hat c(X,Y_t,(π_t)_X)-\tilde c_{π_t}(Y_t)+ε∇V(Y_t))dt+√(2ε)dB_t，X frozen。熵项给出 Langevin drift/noise，weak cost 给出 McKean--Vlasov dependence，projection 给出保边际修正。

4. Well-posedness proof strategy：用 MV / FP / conditional FP 三种表述切换。L(Y|X) 的非线性通过 adapted compactness 控制；L(X|Y) 的条件期望通过 parabolic regularity 获得 Hölder / Sobolev 稳定性。mimicking theorem 用来从 Fokker--Planck 解回到 weak SDE 解。

5. Long-time convergence：能量耗散为 \tilde I(π)=E|\hat c-\tilde c+ε∇log(dπ/dµ⊗ν)|^2。证明路线是紧性 + lower semicontinuity + zero-dissipation characterization。零耗散推出 Euler--Lagrange 条件，再由 convexity 和唯一性识别 optimizer。

Key Insight / Why It Works

最核心贡献是 conditional projection term \tilde c_π(y)。它解决的不是数值稳定性，而是几何一致性：在每个 y 上减去 conditional mean force，使得净速度对 Y 边际没有一阶影响。与此同时，未平均的 \hat c(x,y,π_x) 仍保留每个 x-fiber 的弱输运驱动。这个“先按 X 条件生成力、再按 Y 条件投影”的结构，是文章真正新增的信息。

为什么这套东西能成立：entropy 提供扩散和 Fisher information 控制；固定 ν 让 ∫_X ∂_y logρ(z,y)π^y(dz)=-∇V(y) 这样的抵消成立；adapted topology 让 π_x 的极限行为可控；convexity 让零耗散可以转化为全局最优性，而不是局部 stationarity。

最可能的核心贡献：不是粒子算法，也不是 entropy regularization 本身，而是把 weak transport force 与 marginal projection 组合成一个 well-posed projected McKean--Vlasov dynamics，并证明其 adapted-topology 收敛。

可能只是辅助的部分：DRO sensitivity derivation 更像合理的几何动机，不是完全严密的 Riemannian calculus；数值实验主要是展示动力学不会明显跑偏。粒子离散与 grid projection 属于 engineering prototype，不能证明方法在高维或复杂 weak cost 下有效。

这不是 scaling、retrieval、data coverage 类型的工作；它本质上是 better inductive bias / latent structure alignment：把条件核结构作为状态空间几何的一部分，而不是把 weak cost 硬塞进 ordinary Wasserstein dynamics。

Relation To Prior Work

最接近的是 Conforti--Lacker--Pal 的 projected Langevin dynamics for entropic OT。本文继承了“Langevin diffusion + conditional projection 保持边际 + 长时收敛到 entropic optimizer”的谱系，但把 linear cost 推到 weak cost，因此动力学从普通 projected diffusion 变成 coupled McKean--Vlasov SDE。

和 classical entropic OT / Sinkhorn 的区别：这里没有给出离散 scaling iteration；最优密度虽有 Gibbs form，但 δ_mc(x,π_x) 依赖未知 kernel，固定点是非线性的。本文选择连续时间扩散路线，而不是推广 Sinkhorn。

和 weak OT 静态理论的关系：静态 existence/duality/optimizer representation 主要来自 fundamental theorem of weak OT；本文新增的是 dynamic / gradient-flow 层面的结构。

和 adapted transport / nested distance 的关系：adapted Wasserstein 在这里不是背景装饰，而是必要拓扑。真正不同点是把 adapted topology 的信息结构推进到梯度和投影公式中。

需要注意：d=1 的 projection 公式和 [26] 的正交投影非常接近；d≥2 中文中承认严格正交等式不成立，但仍使用同一公式。这部分更像从一维和保边际需求外推出来的结构，而不是完全闭合的高维几何定理。

Dataset / Evaluation

这是一篇概率/分析论文，evaluation 主要是数值 illustration，不是 benchmark-driven。数值任务覆盖两个典型 sanity check：一维 quadratic OT，以及带 martingale penalty 的 MOT-like 问题。它们验证了动力学在简单场景中能恢复已知结构：monotone rearrangement 和 two-graph / left-monotone coupling。

这些实验支持“动力学方向与理论预期一致”，但不支持更强 claim，例如高维可扩展性、相较 Sinkhorn/LP/MOT solver 的效率优势、复杂 weak cost 的鲁棒性、粒子近似误差或 propagation of chaos。grid projection、finite penalty λ_M、很小 ε 下的行为都带有明显 numerical engineering 成分。

没有真实世界数据，也没有多任务 benchmark。对本文而言这不是致命问题，因为核心 claim 是理论 well-posedness 和 convergence；但若把它解读成实用算法论文，证据明显不足。

Limitation

1. 高维投影的理论缺口：文中在 Remark 3.4 明确说 d≥2 时 (TΠ)^⊥=L^2(µ)⊕L^2(ν) 不再成立，但后续仍采用一维投影公式。这是最实质的 limitation。它可能仍是正确的保边际动力学，但是否真是 adapted Wasserstein steepest descent，文中未充分说明。

2. 假设较强：需要 convexity、linear functional derivative、sub-Gaussian bound、bounded divergence、ν=e^{-V}dy 且 ∇V 有足够 integrability。很多 weak/MOT 成本或奇异边际不一定满足。

3. ε>0 是结构性依赖，不只是数值平滑。熵项提供扩散、密度、Fisher information 和紧性；ε→0 后是否保留有效动力学，文中未解决。

4. well-posedness 是 weak solution / uniqueness in law，未给 strong well-posedness 或粒子系统严格收敛。算法层面的 finite-N bias、grid bias、time discretization error 都没有系统理论。

5. scalability 上限不清楚。每步需要估计 L(Y|X) 和 L(X|Y)，在高维连续空间中这本身就是困难的 conditional density / kernel estimation 问题。论文把静态 weak OT 的非线性优化转移成动态条件化估计问题。

6. MOT 示例使用 penalty λ_M 而非硬 martingale constraint；λ_M→∞ 与 ε→0 的双极限没有证明。数值上看到 two-graph 不等于解决 entropic MOT 或 WMOT 的一般情形。

Takeaway

• 1. weak OT 的动态优化必须尊重 conditional kernel 拓扑；ordinary Wasserstein geometry 不是轻微不合适，而是结构性错误。
• 2. “forward weak force + backward conditional projection”是可以迁移的机制：任何目标依赖 L(Y|X)、但约束固定 Y 边际的问题，都可能需要类似的条件中心化。
• 3. entropy 在这里不仅是 regularizer，也是证明机器：没有 diffusion/Fisher information，conditioning 的稳定性很难闭合。
• 4. 未来真正值得做的是严密高维 adapted tangent projection、ε→0 极限、粒子近似误差，以及不依赖网格条件化的 scalable solver。

一句话总结

这篇论文把 projected Langevin 的保边际梯度流思想推进到熵正则 weak OT，真正贡献是识别并分析了“依赖 L(Y|X) 的弱输运力 + 按 L(X|Y) 条件中心化投影”的 McKean--Vlasov 动力学，但其高维几何投影和算法可扩展性仍是主要未闭合环节。