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

Résumé Technique : Boucles Géodésiques Courtes, Courbure Sectionnelle et Rayon d'Injectivité de la Variété de Stiefel

Énoncé du Problème
La variété de Stiefel, notée $St(n, p)$, est une variété matricielle fondamentale apparaissant en optimisation, dans les méthodes numériques pour les équations différentielles et en statistiques. Un problème central de la géométrie riemannienne sur cette variété est la détermination du rayon d'injectivité, qui fournit une condition suffisante pour la minimalité des géodésiques reliant deux points. Le rayon d'injectivité, $inj(M)$, est défini comme le plus grand rayon $r$ tel que toute géodésique de longueur $\le r$ soit l'unique courbe la plus courte entre ses extrémités.

Pour une variété riemannienne complète, le rayon d'injectivité est déterminé par deux facteurs concurrents : la moitié de la longueur de la boucle géodésique non triviale non nulle ($\ell/2$) et le rayon conjugué. Ce dernier est borné inférieurement par $\pi/\sqrt{K}$, où $K$ est une borne supérieure globale sur la courbure sectionnelle. Bien que le rayon d'injectivité soit connu pour des métriques spécifiques (par exemple, la métrique euclidienne et la métrique canonique), il reste un problème ouvert pour la famille de métriques à un paramètre générale introduite par Hüper et al. (2021), paramétrée par $\beta > 0$. Plus précisément, la longueur exacte des boucles géodésiques les plus courtes et la borne supérieure étroite sur la courbure sectionnelle n'ont pas été pleinement établies pour tout $\beta$, empêchant le calcul exact du rayon d'injectivité à travers tout l'espace de paramètres.

Méthodologie
Les auteurs abordent ce problème en combinant l'analyse de la géométrie différentielle avec des techniques d'analyse matricielle. La méthodologie est divisée en trois composantes principales :

1. Analyse des longueurs de boucles géodésiques ($\ell_\beta$) :

   • Cas $\beta \in (0, 2]$ : Les auteurs utilisent l'observation que les géodésiques sur la variété de Stiefel, lorsqu'elles sont vues comme des courbes dans l'espace euclidien ambiant $\mathbb{R}^{n \times p}$, possèdent des courbures de Frenet constantes. En dérivant des bornes sur la longueur euclidienne de ces courbes et en les traduisant vers la métrique $\beta$, ils établissent des bornes inférieures sur la longueur de la boucle.
   • Cas $\beta > 2$ : L'approche par la courbure de Frenet est insuffisante. Au lieu de cela, les auteurs emploient une preuve d'analyse matricielle. Ils analysent la structure de l'inverse de l'exponentielle de matrice et les groupes d'invariance des matrices orthogonales. En décomposant les matrices de vitesse initiale en sous-blocs correspondant à des propriétés spectrales spécifiques, ils réduisent le problème à l'analyse de boucles géodésiques avec des structures diagonales par blocs spécifiques, prouvant que la longueur ne peut être inférieure à $2\pi$.

2. Bornes de Courbure Sectionnelle ($K_\beta$) :
   Les auteurs dérivent une formule explicite pour la courbure sectionnelle de la variété de Stiefel sous la métrique $\beta$. Ils effectuent des estimations terme à terme sur l'expression de la courbure en utilisant les normes des blocs de matrices ($A$ et $B$) définissant les vecteurs tangents. Cette analyse identifie une lacune dans la littérature existante pour $\beta \in (2/3, 1/\sqrt{2})$, que les auteurs comblent en prouvant que la courbure sectionnelle est bornée supérieurement par 1 dans cet intervalle.

3. Synthèse via le Théorème de Klingenberg :
   En utilisant les longueurs de boucles géodésiques calculées $\ell_\beta$ et les bornes de courbure $K_\beta$, les auteurs appliquent le théorème de Klingenberg, qui stipule que $inj(St_\beta(n, p)) = \min\{ \ell_\beta/2, \text{conj}_{I_{n \times p}}(St_\beta(n, p)) \}$. Ils combinent cela avec la borne inférieure sur le rayon conjugué ($\pi/\sqrt{K_\beta}$) pour déterminer le rayon d'injectivité exact ou des bornes serrées pour diverses plages de $\beta$.

Contributions Clés et Résultats

• Longueur Exacte des Boucles Géodésiques les plus Courtes : Le papier prouve que pour tout $\beta > 0$, la longueur $\ell_\beta$ de la boucle géodésique non triviale la plus courte sur $St_\beta(n, p)$ (pour $2 \le p \le n-1$) est donnée par :
  $$\ell_\beta = \min\{ \sqrt{2\beta}, 1 \} 2\pi$$
  Ce résultat unifie et étend les découvertes précédentes pour la métrique canonique ($\beta=1/2$) et la métrique euclidienne ($\beta=1$), où la longueur était connue pour être $2\pi$. Pour $\beta > 2$, la longueur est montrée comme étant exactement $2\pi$.

• Bornes de Courbure Sectionnelle Étroites : Les auteurs fournissent une borne supérieure étroite $K_\beta$ sur la courbure sectionnelle pour la plage $\beta \in (2/3, 1/\sqrt{2})$, comblant un vide dans la littérature existante. La borne est :
  $$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}$$
  (Note : le papier note que pour $\beta > 1$, la borne de courbure est 1, ce qui est cohérent avec le cas euclidien).

• Détermination du Rayon d'Injectivité :

  • Pour $\beta \in (0, 1/3] \cup [2/3, 1]$, les auteurs déterminent la valeur exacte du rayon d'injectivité :
    • Pour $\beta \in (0, 1/3]$ : $inj = \sqrt{2\beta}\pi$.
    • Pour $\beta \in [2/3, 1]$ : $inj = \pi$.
  • Pour $\beta \in (1/3, 2/3)$, la valeur exacte n'est pas déterminée en raison de l'absence d'une expression explicite pour le rayon conjugué. Cependant, le papier fournit des bornes supérieures et inférieures serrées, montrant que la valeur diffère de moins de 3 % de la valeur conjecturée dans [33].
  • Pour $\beta > 1$, le rayon d'injectivité est borné supérieurement par $\pi$.

• Cas Particuliers :

  • Pour $p=1$ (sphère unité), $inj = \pi$.
  • Pour $p=n$ (groupe orthogonal $O(n)$), le rayon conjugué est explicitement calculé comme étant $\sqrt{8\beta}\pi$ pour $n=3$ et $\sqrt{4\beta}\pi$ pour $n>3$, conduisant à un rayon d'injectivité de $\sqrt{2\beta}\pi$.

Signification
Le papier prétend résoudre le problème de la détermination du rayon d'injectivité pour la variété de Stiefel sous une large gamme de métriques $\beta$, un paramètre crucial pour les algorithmes d'optimisation sur les variétés. En établissant la longueur exacte des boucles géodésiques les plus courtes pour tout $\beta$ et en fermant la lacune des bornes de courbure sectionnelle, les auteurs fournissent les données géométriques nécessaires pour appliquer le théorème de Klingenberg.

Le travail est modeste concernant l'intervalle $\beta \in (1/3, 2/3)$, reconnaissant que le rayon d'injectivité exact reste une conjecture dans cette plage, soutenue par des expériences numériques dans la littérature antérieure ; cependant, le papier renforce la conjecture pour la métrique canonique ($\beta=1/2$) en fournissant des expressions explicites pour les champs de Jacobi, confirmant que les points conjugués identifiés sont bien les premiers points conjugués. Les résultats fournissent une base rigoureuse pour assurer la minimalité des géodésiques dans les applications d'optimisation sur les variétés et d'ajustement de données.