Shortest Geodesic Loops, Sectional Curvature, and Injectivity Radius of the Stiefel Manifold

技术摘要：斯蒂费尔流形（Stiefel Manifold）的最短测地回路、截面曲率与注入半径

问题陈述
斯蒂费尔流形，记作 $St(n, p)$，是出现在优化、微分方程数值方法和统计学中的一种基础矩阵流形。黎曼几何中关于该流形的一个核心问题是确定其注入半径（injectivity radius），这为连接两点的测地线是否为最短曲线提供了充分条件。注入半径 $inj(M)$定义为最大的半径$r$，使得每条长度$\le r$ 的测地线都是连接其端点的唯一最短曲线。

对于一个完备的黎曼流形，注入半径由两个相互竞争的因素决定：最短非平凡测地回路长度的一半 ($\ell/2$) 以及共轭半径（conjugate radius）。后者由 $\pi/\sqrt{K}$ 给出下界，其中 $K$ 是截面曲率（sectional curvature）的全局上界。虽然注入半径在特定度量下（例如欧几里得度量和规范度量）是已知的，但对于 Hüper 等人 (2021) 引入的由参数 $\beta > 0$ 参数化的广义一参数族度量，该问题仍是一个开放问题。具体而言，最短测地回路的精确长度以及截面曲率的锐利上界在所有 $\beta$ 取值下均未被完全建立，这阻碍了在整个参数空间内精确计算注入半径。

方法论
作者通过结合微分几何分析与矩阵分析技术来解决这一问题。该方法分为三个主要部分：

1. 测地回路长度分析 ($\ell_\beta$)：

   • 情形 $\beta \in (0, 2]$： 作者利用了这样一个观察：当斯蒂费尔流形上的测地线被视为其环境欧几里得空间 $\mathbb{R}^{n \times p}$ 中的曲线时，它们具有常数 Frenet 曲率。通过推导这些曲线在欧几里得长度上的界限并将其转化为 $\beta$-度量下的长度，作者建立了回路长度的下界。
   • 情形 $\beta > 2$： Frenet 曲率方法在此不再适用。相反，作者采用了矩阵分析证明法。他们分析了矩阵指数逆的结构以及正交矩阵的不变群。通过将初始速度矩阵分解为对应于特定谱性质的子块，他们将问题简化为分析具有特定对角块结构的测地回路，从而证明回路长度不会小于 $2\pi$。

2. 截面曲率界限 ($K_\beta$)：
   作者推导了 $\beta$-度量下斯蒂费尔流形截面曲率的显式公式。他们利用定义切向量的矩阵块（$A$ 和 $B$）的范数，对曲率表达式进行了逐项估计。该分析识别出既往文献中 $\beta \in (2/3, 1/\sqrt{2})$ 区间的空白，并通过证明在该区间内截面曲率上界为 1 来填补了这一空白。

3. 通过 Klingenberg 定理进行综合：
   利用计算出的回路长度 $\ell_\beta$ 和曲率界限 $K_\beta$，作者应用了 Klingenberg 定理，该定理指出 $inj(St_\beta(n, p)) = \min\{ \ell_\beta/2, \text{conj}_{I_{n \times p}}(St_\beta(n, p)) \}$。他们结合共轭半径的下界 ($\pi/\sqrt{K_\beta}$) 来确定各种 $\beta$ 范围内的精确注入半径或紧凑界限。

主要贡献与结果

• 最短测地回路的精确长度： 本文证明了对于所有 $\beta > 0$，对于 $2 \le p \le n-1$ 的 $St_\beta(n, p)$，其最短非平凡测地回路的长度 $\ell_\beta$ 为：
  $$\ell_\beta = \min\{ \sqrt{2\beta}, 1 \} 2\pi$$
  这一结果统一并扩展了以往关于规范度量 ($\beta=1/2$) 和欧几里得度量 ($\beta=1$) 的发现，后者的长度已知为 $2\pi$。对于 $\beta > 2$，其长度被证明精确为 $2\pi$。

• 锐利截面曲率界限： 作者为 $\beta \in (2/3, 1/\sqrt{2})$ 范围内的截面曲率提供了一个锐利上界，填补了现有文献中的空白。该界限为：
  $$K_\beta = \begin{cases} \frac{1}{4\beta} & 0 < \beta \le \frac{1}{6}(4-\sqrt{10}) \\ \frac{4-3\beta}{2} & \frac{1}{6}(4-\sqrt{10}) < \beta \le \frac{2}{3} \\ 1 & \frac{2}{3} \le \beta \le 1 \end{cases}$$
  （注：文中指出，对于 $\beta > 1$，曲率界限为 1，这与欧几里得情形一致）。

• 注入半径的确定：

  • 对于 $\beta \in (0, 1/3] \cup [2/3, 1]$，作者确定了精确值：
    • 对于 $\beta \in (0, 1/3]$：$inj = \sqrt{2\beta}\pi$。
    • 对于 $\beta \in [2/3, 1]$：$inj = \pi$。
  • 对于 $\beta \in (1/3, 2/3)$，由于缺乏共轭半径的显式表达式，精确值尚未确定。然而，论文提供了紧凑的上界和下界，表明其值与文献 [33] 中的猜想值偏差不超过 3%。
  • 对于 $\beta > 1$，注入半径的上界为 $\pi$。

• 特殊情形：

  • 对于 $p=1$（单位球面），$inj = \pi$。
  • 对于 $p=n$（正交群 $O(n)$），共轭半径被明确计算为 $n=3$ 时为 $\sqrt{8\beta}\pi$，对于 $n>3$ 为 $\sqrt{4\beta}\pi$，从而导致注入半径为 $\sqrt{2\beta}\pi$。

意义
本文声称解决了在广泛的 $\beta$-度量下确定斯蒂费尔流形注入半径的问题，这对于基于流形的优化算法至关重要。通过为所有 $\beta$ 建立最短测地回路的精确长度并填补截面曲率界限的空白，作者为应用 Klingenberg 定理提供了必要的几何数据。

尽管在 $\beta \in (1/3, 2/3)$ 区间内的研究较为保守，承认该范围内的精确注入半径仍是一个猜想，但论文通过提供 Jacobi 场的显式表达式，加强了对该猜想的支持，并证实了所识别的共轭点确实是第一共轭点。这些结果为确保基于流形的优化和数据拟合应用中的测地线最小性提供了坚实的理论基础。