精读笔记

Problem Setting

论文标题：Obstacles and Singularities of Riemannian Distance Functions（arXiv preprint / 2026-06-24）。问题是在 R^n 中移除一个紧障碍物 O 后，研究从固定目标 k0 出发的 Riemannian 距离函数 d 在 X=closure(R^n\O) 上的奇异集。关键矛盾是：距离函数由 eikonal 方程控制，局部看像标准 Hamilton-Jacobi value function；但障碍物引入状态约束，使最短轨道可能贴边、切触、绕行、甚至以相同初速度在之后分叉。无约束 cut locus 的“多条 geodesic ⇔ 不可微”直觉在这里需要重写。论文真正解决的是：障碍物是否必然制造奇异性，以及这些奇异性是否像 semiconcave HJ 解那样传播。

Motivation

已有路线不够的点在于，它们要么关注 manifolds with boundary 中 geodesics 的存在/正则性，要么处理欧氏障碍距离的局部性质；但没有给出 constrained distance 的全局 singular set 结构。作者的核心观察是：障碍物不是一个局部扰动，而是改变了 wavefront 的拓扑演化。一个高距离等值集若完全光滑，它沿最短轨道向目标推进时应该形成一个无交叉的光滑波前；但有限障碍物会迫使这个波前在某个时刻接触边界，而接触点同时满足“轨道切边界”和“距离梯度法向于等值面/边界”的条件，两者与 eikonal 非退化性冲突。缺口不是 semiconcavity 本身，而是如何把边界约束、reachable gradients、minimizer 几何和全局拓扑拼成一个不可避免的奇异性机制。

Core Idea

核心思想可以概括为：把奇异性从“geodesic 是否唯一”转为“reachable gradient 是否唯一”，再用 geodesic flow 推动正则等值集，证明全正则波前不可能绕过障碍约束。这个转化很关键，因为在有障碍时，多条 minimizer 不一定对应多个初速度；它们可能先共用一段轨道，甚至在边界后才分叉。因此真正稳定的对象不是 minimizer 数量，而是 D*d 的多值性。

与 prior 的本质区别在于，这篇不是只把欧氏结果替换成 Riemannian metric，而是把 constrained cut-locus 问题放进 semiconcave value-function 框架中处理。障碍物带来的 inductive bias 是“波前不能无代价穿过/切过边界”：如果高能级波前没有奇异点，它就会像一个 C^1 hypersurface 一样沿 Hamiltonian characteristics 推进；但推进到障碍边界时，最短轨道的切触条件和距离梯度的 eikonal 条件不可同时满足。这是 Theorem 1.1 的核心机制。

Method

1. 边界 semiconcavity 扩张：距离函数在内部有线性 semiconcavity，但靠近障碍边界只有 1/2-fractional semiconcavity。作者使用 semiconcave extension 保留边界 reachable gradients，使边界点也能用 D*d 定义奇异性。这解决的是边界 nonsmooth analysis 不稳定的问题。

2. minimizer-gradient 对应：Proposition 3.2 证明任意最短轨道的速度满足 γdot(t) ∈ -A^{-1}D*d(γ(t))，并且轨道内部点正则；反过来，每个 reachable gradient 都对应某条 minimizer 的初速度。这一步把分析对象从曲线族转成超微分几何，是后续传播和全局拓扑论证的支点。

3. 高等值集必奇异：假设某个高能级连通分支 K0 无奇异点，则由 nonsmooth implicit function theorem 得到它是 C^1 hypersurface。沿唯一 minimizer flow 推进 K0，利用 invariance of domain 得到无交叉波前。若所有轨道能到目标不碰障碍，会和靠近目标的拓扑覆盖矛盾；若先碰障碍，则波前与边界相切，导致 Dd 同时法向于边界而 A^{-1}Dd 切向于边界，和 <A^{-1}Dd,Dd>=1 矛盾。

4. 奇异性局部传播：对内点奇异 x，利用 d_O 的 reachable gradients 与 d 的 reachable gradients 不相交，推出 D*d 不能填满 Riemannian unit co-sphere；因此 ∂D^+d(x)\D*d(x) 非空，可调用 Albano-Cannarsa 的 semiconcave singularity propagation 定理得到非平凡 Lipschitz singular arc。

Key Insight / Why It Works

最核心贡献不是 Riemannian technical generalization，而是识别出“障碍物强制波前奇异”的拓扑-微分机制。Theorem 1.1 的证明很干净：如果高能级波前全正则，那么 characteristics 从波前向目标方向形成一个无交叉流；这种流要么把一个闭 hypersurface 收缩到目标点，要么在中途碰到障碍。前者与障碍附近存在到目标的最短路径相冲突，后者与边界切触条件和 eikonal 方程冲突。因此奇异点不是 generic accident，而是 constrained geometry 的必要产物。

真正有效的工具是 semiconcavity + reachable gradients，而不是 geodesic uniqueness 本身。因为有障碍时，“多条 minimizer”这个信号太粗：Section 5 的点 P 有两条 minimizer 但初速度相同，因此 D*d 仍是 singleton，距离函数正则。这说明论文对奇异性的定义抓住了正确变量：不是路径条数，而是最优反馈方向的多值性。

传播结果的关键 insight 是把“离障碍最近方向”当作一个缺失的 covector：Lemma 6.1 说明 D*d 与 D*d_O 不相交，几何上因为到障碍的最短轨道正交撞边界，而到目标的 constrained minimizer 触边界只能切向。这制造了 D^+d 边界上不属于 D*d 的点，从而触发经典 semiconcave singularity propagation。这里不是 scaling，也不是 engineering；增益完全来自正确的 nonsmooth geometric formulation。

辅助性较强的是二维欧氏 Section 7。它给出更直观的“两侧绕障碍”刻画，但依赖 n=2、A=I、障碍连通，推广性有限。其价值更多是解释 sharpness 和直觉，而不是主定理的必要组件。

Relation To Prior Work

最接近的是三条谱系：Riemannian manifolds with boundary 中的 geodesic/minimizer 理论；Hamilton-Jacobi/state constraint value functions；semiconcave function singularity propagation。论文的新增信息在于把这三条线合在 constrained distance 的 singular set 上，而不是单独研究 minimizer 或 PDE 正则性。

相对于无约束 cut locus 理论，本质差异是 singularity 不再等价于 minimizer 数量，而等价于初始最优方向/ reachable gradient 的多值性。障碍物允许“同初速度多 minimizer”这种无约束情形不存在的退化，因此传统 cut locus 语言不足。

相对于作者前作和欧氏障碍结果，实质创新有两点：一是推广到一般 smooth Riemannian metric，并保持 minimizer-gradient 对应；二是给出每个足够高等值集连通分支必含奇异点的全局结论。看似新的 propagation 形式部分借用了 Albano-Cannarsa 的 semiconcave 传播定理，真正新增的是验证 constrained setting 下触发该定理的几何条件。

Section 5 的反例也很重要：它否定了“障碍边界必有奇异点”的自然猜想，说明 prior 中 convex Euclidean obstacle 的边界奇异结论不能简单推广到非凸情形。

Dataset / Evaluation

这是一篇纯数学论文，没有 dataset、benchmark 或实验评估。所谓 evaluation 是定理覆盖范围与 sharpness examples：Theorem 1.1 验证“障碍必然产生远场奇异性”；Theorem 1.2 验证“内点奇异性非孤立”；Section 5 验证“边界奇异性非必然”；torus 例子验证“局部传播不能自动升级为全局无界传播”。这些例子足够支持论文的结构性 claim，但不支持更强的 conjecture，例如任意内点奇异弧最终必然通向无穷远。文中对非 C^2 障碍、多个障碍、高维复杂拓扑的覆盖没有展开。

Limitation

主要限制不是技术细节，而是结论的适用边界。第一，障碍物需要 C^2 边界；semiconcavity 扩张、切触条件和 minimizer 正则性都依赖这个光滑性，非光滑障碍下机制可能出现角点反射/多法向，结论需重做。第二，Theorem 1.2 只给局部 Lipschitz singular arc，不保证全局延伸；torus 例子说明奇异性可以撞到障碍后停止。第三，二维欧氏的精细结构明显依赖“只有两侧可绕”的平面拓扑，高维或一般 metric 下不成立。

另一个隐含前提是 metric 只通过 A(x) 的 smooth uniformly positive definite 结构进入；如果考虑 Finsler、sub-Riemannian 或一般控制系统，minimizer-gradient 对应和正交/切触矛盾未必保持。文中未充分说明多障碍、多连通 X、或存在 conjugate phenomena 时 singular set 的全局拓扑如何变化。增益来源很清楚，不是 engineering/scaling，而是几何约束与 semiconcavity 理论的正确耦合；但上限也明显：它给的是存在性和局部传播，不是 singular set 的完整分类。

Takeaway

• 1. 对 constrained distance，正确的奇异性变量是 reachable gradients / optimal initial covectors，而不是 minimizer 数量；这点可迁移到更一般 state-constrained control。
• 2. 障碍物强制远场波前产生奇异性：这不是 genericity 结论，而是由无交叉 characteristic flow 与边界切触条件之间的结构性矛盾导致的。
• 3. semiconcave singularity propagation 在有状态约束时仍可用，但必须先证明存在 ∂D^+d\D*d 的 covector；这里 d_O 与 d 的 reachable gradients disjoint 是关键 trick。
• 4. 后续真正值得做的是全局奇异弧理论：什么时候局部 propagation 会延伸到无穷远，什么时候会被障碍吸收，以及高维 constrained cut locus 是否存在可计算/可分类的拓扑结构。

一句话总结

这篇论文把带障碍 Riemannian 距离函数的奇异性从 classical cut locus 问题重构为 semiconcave Hamilton-Jacobi value function 的 constrained wavefront 问题，核心贡献是证明障碍物必然制造远场奇异性且内点奇异性局部传播。