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

技術要約：スティフェル多様体における最短測地ループ、断面曲率、および注入半径

問題提起
スティフェル多様体（$St(n, p)$と表記）は、最適化、微分方程式の数値解法、統計学において現れる基本的な行列多様体である。この多様体上のリーマン幾何学における中心的な問題は、**注入半径（injectivity radius）**を決定することであり、これは2点間を結ぶ測地線の最小性を保証するための十分条件となる。注入半径 $inj(M)$は、長さ$r$以下のすべての測地線が、その端点間の唯一の最短曲線となる最大の半径$r$ として定義される。

完備なリーマン多様体において、注入半径は、最短の非自明な測地ループの長さの半分（$\ell/2$）と、共役半径という2つの競合する要因によって決定される。後者は、断面曲率のグローバルな上界 $K$ を用いて $\pi/\sqrt{K}$ によって下から抑えられる。注入半径は特定の計量（例：ユークリッド計量や標準的計量）については既知であるが、Hüperら（2021年）によって導入された、パラメータ $\beta > 0$ でパラメータ化された一般的な1パラメータ族の計量については、依然として未解決の問題である。具体的には、すべての $\beta$ に対して最短測地ループの長さと断面曲率の鋭い上界が完全には確立されておらず、それが全パラメータ空間にわたる注入半径の正確な計算を妨げている。

手法
著者らは、微分幾何学的解析と行列解析的手法を組み合わせることで、この問題に取り組んでいる。手法は主に以下の3つの構成要素に分かれている。

1. 測地ループの長さ（$\ell_\beta$）の解析：

   • $\beta \in (0, 2]$ の場合： 著者らは、スティフェル多様体上の測地線が、周囲のユークリッド空間 $\mathbb{R}^{n \times p}$ における曲線として捉えたとき、一定のフレネ・曲率を持つという観察を利用している。これらの曲線のユークリッド的な長さを導出し、それを $\beta$-計量へと翻訳することで、ループの長さの下限を確立している。
   • $\beta > 2$ の場合： フレネ曲率によるアプローチでは不十分である。代わりに、著者らは行列解析的な証明を用いている。彼らは行列指数の逆行列の構造と直交行列の不変群を分析している。初期速度行列を特定のスペクトル特性に対応するサブブロックに分解することで、問題を特定のブロック対角構造を持つ測地ループの解析へと還元し、長さが $2\pi$ 未満にはなり得ないことを証明している。

2. 断面曲率の境界（$K_\beta$）：
   著者らは、$\beta$-計量の下でのスティフェル多様体の断面曲率の明示的な公式を導出している。彼らは、接ベクトルを定義する行列ブロック（$A$ および $B$）のノルムを用いた、曲率式の項別推定を行っている。この解析により、$\beta \in (2/3, 1/\sqrt{2})$ における既存文献のギャップが特定され、著者らはこの区間において断面曲率が1によって上から抑えられることを証明することで、このギャップを埋めている。

3. クリンゲンベルグの定理による統合：
   計算されたループの長さ $\ell_\beta$ と曲率の境界 $K_\beta$ を用いて、著者らはクリンゲンベルグの定理を適用している。この定理は、$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$ において、$St_\beta(n, p)$ （$2 \le p \le n-1$ の場合）の最短の非自明な測地ループの長さ $\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$ を提供している。その境界は以下の通りである：
  $$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$ について最短の測地ループの正確な長さを確立し、断面曲率の境界のギャップを埋めることで、著者らはクリンゲンベルグの定理を適用するために必要な幾何学的データを提供している。

$\beta \in (1/3, 2/3)$ の区間については、正確な注入半径が依然として予想段階であることを認めつつも、控えめな記述にとどまっており、先行研究の数値実験によって支持されている。しかし、本論文はヤコビ場の明示的な表現を提供することで、標準的計量（$\beta=1/2$）に関する予想を強化し、特定された共役点が確かに最初の共役点であることを確認している。これらの結果は、多様体ベースの最適化およびデータフィッティングへの応用において、測地線の最小性を保証するための厳密な基礎を提供するものである。