精读笔记
Problem Setting
论文实际解决的是：在高维 product state space 上，给定初末分布 μ、ν 和一步代价/转移约束，寻找一个连接它们的低成本 Markov process。真正困难点是每个时间步都需要表示相邻状态 coupling π^s∈P(Ω×Ω)，而 Ω 本身已经高维，π^s 的显式表示比单个分布还更不可承受。
以前路线大致卡在两端：Eulerian dynamic OT 在低维网格上 convex 但维度爆炸；neural ODE / flow matching / sample-based 方法能跑高维，但训练非凸，通常没有对原问题的 certified lower bound，也很难区分 sampling、optimization、ODE discretization 和 approximation error。Schrödinger bridge 则引入参考过程和 entropy regularization，问题结构不同。
这篇的关键矛盾是：Markov process 优化需要跨时间的一致 coupling，但高维可计算性只允许保留局部统计。论文的策略不是近似采样完整路径，而是把问题降到“局部 coupling statistics 的 convex feasibility/optimization”。
Motivation
已有路线缺的不是又一个高维 transport solver，而是一个能同时保留 convexity、动态一致性和可扩展局部表示的框架。作者看到静态 OT 中 cluster marginal / moment relaxation 已经能把高维 coupling 压缩成局部边缘或低阶矩；自然的问题是，这套思想能不能沿时间方向扩展到 Markov process。
核心观察是：Markov kernel 不是必须直接优化。只要用相邻 joint law π^s 表示一步转移，Markov evolution 就变成相邻 coupling 的左右边缘匹配。这样动态问题变成多阶段 coupling problem，而每个 coupling 都可以套用静态 local relaxation。
关键缺口是 reconstruction：local relaxation 给的是局部统计和 lower bound，不是完整 dynamics。论文因此把 BB 特例中的 dual velocity recovery 和一般 Markov setting 下的 parametric kernel fitting 放进框架里，试图把 relaxed statistics 转回可执行过程。
Core Idea
最核心的思想是把“优化 Markov kernel”改写成“优化一串相邻时间 coupling”，再对每个 coupling 做局部化。全局路径 law 被拆成 π^0,...,π^{T-1}，时间方向只通过右边缘等于下一步左边缘来连接；空间方向只保留 cluster / cluster-pair marginals 或 moments。这样原来指数维的动态优化被组织成一组低维局部变量上的 convex program。
这个做法本质上引入了 locality inductive bias：相信最优过程的低阶局部统计足以描述关心的结构，或者至少足以给出有用 lower bound。它和 neural flow 类方法的区别不是参数化更巧，而是目标仍是 convex relaxation；和传统 grid-based BB 的区别是不用 full space-time grid；和静态 OT relaxation 的区别是多了时间一致性链条，因此不是独立求一堆 static transport。
理论直觉上它有效，是因为任何真实 Markov process 都会投影成这些局部变量，所以 relaxed feasible set 是原 feasible set 的外近似；代价若局部分解，目标也能在这些变量上精确或近似表达。它牺牲的是全局 representability，换来的是 convex lower bound 和可计算的低阶中间统计。
Method
1. Sequential coupling：用 π^s(dx,dy)=ρ^s(dx)K^s(x,dy) 替代 K^s 和 ρ^s 的联合优化。它解决的是 kernel/distribution 交替依赖不好处理的问题。核心变化是 Markov dynamics 被转成相邻 coupling 的边缘匹配，所有约束都可以写在 measure/coupling 层面。
2. Marginal relaxation：离散或网格状态下，只保留每个时间步 coupling 的 cluster marginals π_A^s 和 cluster-pair marginals π_AB^s。它解决 full coupling 表示爆炸的问题。必要约束包括局部一致性、endpoint 投影、跨时间质量守恒、局部 kernel constraint、局部非负性，以及可选的全局 PSD positivity。核心变化是把全局可实现性换成一组必要条件。
3. Moment relaxation：连续状态下，用有限基函数下的 moment matrices M_A^s、M_AB^s 替代局部边缘。它解决连续空间局部边缘仍需离散化的问题。moment consistency 保证重复 monomial/product function 的条目一致，PSD 保证在 retained function span 上的二阶正性。代价若能投影到 basis span，目标就是线性的。
4. BB 特例：当无 kernel 约束且一步代价为 ||x-y||^p/(Δt)^{p-1}，sequential-coupling problem 在任意给定时间网格上与 Benamou-Brenier geodesic 的 grid marginals 等价，最优值无时间离散误差。对 p=2，dual potentials 的梯度给 velocity。
5. General dynamics recovery：对非 BB 的 constrained Markov process，relaxation 只给局部统计，因此作者拟合一个参数化 transition family，例如 Glauber kernel，使其一步预测的局部 marginals 匹配 recovered statistics。这一步解决可执行性，但不是严格的 convex recovery。
Key Insight / Why It Works
最重要的 insight 是：动态 Markov optimization 的高维困难可以拆成两个相对可控的结构：时间上只需要相邻 coupling 的一致性，空间上只需要局部 interaction 的统计表示。只要代价和约束是局部的，完整 coupling 的很多信息对目标并不直接可见；优化这些可见的局部统计就能给 lower bound，并在结构良好的问题上接近真实过程。
真正的核心贡献是 sequential-coupling + cluster relaxation 的组合。单看 local marginal / moment relaxation 并不新，来自静态 OT 和 sparse moment-SOS；单看 BB dual velocity recovery 也基于经典 duality。新增的信息在于：把这些局部 convex relaxation 放到多阶段 Markov coupling 链上，并说明 BB 和 constrained finite-state process 都可作为实例。
最可能有效的部分是 better inductive bias，而不是通用 scaling magic。它假设 distributional dynamics 的关键信息是低阶、局部、稀疏的；Gaussian sparse precision、Ginzburg-Landau lattice、Ising local update 都正好符合这个假设。因此实验成功主要说明该 inductive bias 在这些结构化系统上合理。
辅助成分包括 global PSD positivity、SOS dual velocity、max-entropy reconstruction、Glauber fitting。它们提升可用性，但不是同一层级的核心。尤其 kernel fitting 更像把 relaxed statistics 投影到手选 dynamics family；如果 family 不对或局部统计不足，过程可能局部匹配但全局偏离。
这里没有数据覆盖、retrieval、memory reuse 或 test-time compute 的成分；这是典型 convex relaxation + sparse structure exploitation。所谓 scalability 主要来自 sparse representation，而不是算法复杂度突破。文中未充分说明 relaxation tightness 的一般条件，因此“能恢复 dynamics”这个 claim 在 BB 特例较强，在一般 Markov process 中要弱得多。
Relation To Prior Work
最接近的是三条线。第一是 Benamou-Brenier dynamic OT：本文把它作为 sequential-coupling 的无约束 kinetic-cost 特例，并证明 grid-time marginals 和最优值可恢复。区别在于本文随后做空间 relaxation，而传统 BB 依赖 full spatial discretization。
第二是高维 neural OT / flow matching / TrajectoryNet 类方法：这些方法通过参数化 velocity 或 flow map 绕开网格，但通常是非凸经验优化。本文的本质差异是保留 convex program 和 lower bound，但代价是只恢复低阶统计，完整 dynamics 需要额外 reconstruction。
第三是静态 OT 的 cluster marginal / moment relaxation 和 sparse SOS。本文看似新的一部分其实是这些已有思想的动态化：对每个相邻 coupling 复制局部 relaxation，再用时间一致性约束串起来。实质创新在于把静态 local OT relaxation 变成 Markov process optimization 的统一框架，并给出 BB dual velocity 与 constrained kernel fitting 两种 recovery 路径。
和 Schrödinger bridge 的差异也明确：本文不需要 reference process，也不使用 entropic regularization；它优化一般 additive one-step cost，并允许 hard transition constraints。但这也意味着没有 SB 那种通过 reference dynamics 自然得到 path measure 的机制。
Dataset / Evaluation
实验选择是结构化、机制验证型，而不是大规模真实世界 benchmark。Gaussian-to-Gaussian 用闭式 geodesic 验证 moment 和 velocity recovery，这是最干净的 sanity check；Gaussian-to-Ginzburg-Landau 验证非高斯连续目标下的低维 marginal reconstruction；Ising single-spin process 验证 constrained finite-state Markov setting 下的 local marginal relaxation + Glauber fitting。
这些实验支持的 claim 是有限但清楚的：在局部结构强、代价局部分解、低阶统计有信息量的问题中，relaxation 能恢复中间局部统计，并可通过 dual 或 fitting 生成可执行近似 dynamics。它们没有充分验证的是：在弱局部性、长程依赖、复杂非局部转移约束下 relaxation 是否仍紧。
benchmark 没有明显 leakage 问题，因为不是学习型 benchmark；但 evaluation 存在归因限制。Gaussian 实验中真实过程本来由低阶矩完全决定，moment relaxation 天然占优。Ising 实验中 fitted Glauber family 与生成/目标结构高度匹配，因此 reproducing local marginals 不等于证明一般 kernel recovery。Ginzburg-Landau 中 maximum-entropy reconstruction 和 velocity rollout 的质量受 moment truncation 与数值后处理影响，文中没有系统拆分误差来源。
Limitation
第一，核心前提是 locality preservation。端点有局部结构并不自动意味着最优 interpolation 的中间 law 也局部可表示。论文在结论中承认这点仍未解决；这不是边缘问题，而是 relaxation tightness 的基础。
第二，relaxation 只给 lower bound 和 low-order statistics。除 BB 特例外，从局部统计到完整 Markov kernel 没有一般唯一性或最优性保证。Glauber fitting 是合理工程步骤，但它把问题转移成“选择正确 parametric family 并局部拟合”。如果真实 dynamics 不是该 family，拟合结果可能只是局部投影。
第三，SDP scaling 仍然是硬上限。cluster size、basis degree、cluster graph radius、time horizon 都会放大 PSD block 和 affine constraints。Chordal decomposition 缓解了数值规模，但方法并没有变成真正适合超高维连续系统的通用 solver。
第四，dual velocity recovery 的理论条件偏强。需要 quadratic BB、足够正则的 dual potential、connected support、绝对连续性等；moment relaxation 的 SOS dual potential 只是近似，且可能非唯一。velocity rollout 的稳定性和边界行为没有一般保证。
第五，增益归因不完全清楚。实验改善可能主要来自问题低阶结构强、cluster graph 覆盖了关键相关性、Gaussian/Ising family 与方法 bias 高度匹配，而不是 relaxation 框架在一般高维 Markov optimization 上都有效。
Takeaway
- 1. 这篇真正推动的是把 high-dimensional dynamic distribution optimization 从 full path / full grid 表示，转向 sequential local statistics optimization。
- 这个视角可迁移到任何“时间局部 + 空间稀疏”的分布演化问题。
- 2. 最值得保留的 insight 是：先用 convex relaxation 得到 certified lower bound 和局部中间统计，再用问题特定的 recovery map 转成可执行 dynamics。
- 不要把 recovery 和 relaxation 混为一谈。
一句话总结
这篇论文把静态高维 OT 的局部 convex relaxation 动态化为 sequential-coupling Markov process optimization，核心贡献是用 locality 换取 convex lower bound 和低阶中间统计，而不是提供一个通用完整 dynamics recovery solver。