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

Resumen Técnico: Bucles Geodésicos Cortos, Curvatura Seccional y Radio de Inyectividad de la Variedad de Stiefel

Planteamiento del Problema
La variedad de Stiefel, denotada como $St(n, p)$, es una variedad matricial fundamental que aparece en la optimización, los métodos numéricos para ecuaciones diferenciales y la estadística. Un problema central de la geometría riemanniana en esta variedad es determinar el radio de inyectividad, el cual proporciona una condición suficiente para la minimidad de las geodésicas que conectan dos puntos. El radio de inyectividad, $inj(M)$, se define como el radio más grande $r$ tal que toda geodésica de longitud $\le r$ es la curva más corta única entre sus extremos.

Para una variedad riemanniana completa, el radio de inyectividad está determinado por dos factores contrapuestos: la mitad de la longitud del bucle geodésico no trivial más corto ($\ell/2$) y el radio conjugado. Este último está acotado inferiormente por $\pi/\sqrt{K}$, donde $K$ es un límite superior global de la curvatura seccional. Aunque el radio de inyectividad se conoce para métricas específicas (por ejemplo, la métrica euclídea y la métrica canónica), sigue siendo un problema abierto para la familia general de métricas de un parámetro introducida por Hüper et al. (2021), parametrizada por $\beta > 0$. Específicamente, la longitud exacta de los bucles geodésicos más cortos y el límite superior agudo de la curvatura seccional no se establecieron completamente para todo $\beta$, lo que impidió el cálculo exacto del radio de inyectividad a través de todo el espacio de parámetros.

Metodología
Los autores abordan este problema combinando el análisis de la geometría diferencial con técnicas de análisis matricial. La metodología se divide en tres componentes principales:

1. Análisis de las Longitudes de los Bucles Geodésicos ($\ell_\beta$):

   • Caso $\beta \in (0, 2]$: Los autores utilizan la observación de que las geodésicas en la variedad de Stiefel, cuando se ven como curvas en el espacio euclídeo ambiente $\mathbb{R}^{n \times p}$, poseen curvaturas de Frenet constantes. Al derivar límites para la longitud euclídea de estas curvas y trasladarlos a la métrica $\beta$, establecen límites inferiores para la longitud del bucle.
   • Caso $\beta > 2$: El enfoque de la curvatura de Frenet es insuficiente. En su lugar, los autores emplean una prueba de análisis matricial. Analizan la estructura de la inversa de la exponencial de la matriz y los grupos de invariancia de las matrices ortogonales. Al descomponer las matrices de velocidad inicial en subbloques correspondientes a propiedades espectrales específicas, reducen el problema al análisis de bucles geodésicos con estructuras de bloques diagonales específicas, demostrando que la longitud no puede ser menor que $2\pi$.

2. Límites de la Curvatura Seccional ($K_\beta$):
   Los autores derivan una fórmula explícita para la curvatura seccional de la variedad de Stiefel bajo la métrica $\beta$. Realizan estimaciones término a término de la expresión de la curvatura utilizando las normas de los bloques de matrices ($A$ y $B$) que definen los vectores tangentes. Este análisis identifica un vacío en la literatura existente para $\beta \in (2/3, 1/\sqrt{2})$, el cual los autores cierran demostrando que la curvatura seccional está acotada superiormente por 1 en este intervalo.

3. Síntesis mediante el Teorema de Klingenberg:
   Utilizando las longitudes de bucle calculadas $\ell_\beta$ y los límites de curvatura $K_\beta$, los autores aplican el teorema de Klingenberg, que establece que $inj(St_\beta(n, p)) = \min\{ \ell_\beta/2, \text{conj}_{I_{n \times p}}(St_\beta(n, p)) \}$. Combinan esto con el límite inferior del radio conjugado ($\pi/\sqrt{K_\beta}$) para determinar el radio de inyectividad exacto o límites ajustados para varios rangos de $\beta$.

Contribuciones Clave y Resultados

• Longitud Exacta de los Bucles Geodésicos Cortos: El artículo demuestra que para todo $\beta > 0$, la longitud $\ell_\beta$ del bucle geodésico no trivial más corto en $St_\beta(n, p)$ (para $2 \le p \le n-1$) viene dada por:
  $$\ell_\beta = \min\{ \sqrt{2\beta}, 1 \} 2\pi$$
  Este resultado unifica y extiende hallazgos previos para la métrica canónica ($\beta=1/2$) y la métrica euclídea ($\beta=1$), donde la longitud se conocía como $2\pi$. Para $\beta > 2$, se demuestra que la longitud es exactamente $2\pi$.

• Límites Agudos de la Curvatura Seccional: Los autores proporcionan un límite superior agudo $K_\beta$ sobre la curvatura seccional para el rango $\beta \in (2/3, 1/\sqrt{2})$, llenando un vacío en la literatura existente. El límite es:
  $$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}$$
  (Nota: El artículo señala que para $\beta > 1$, el límite de la curvatura es 1, consistente con el caso euclídeo).

• Determinación del Radio de Inyectividad:

  • Para $\beta \in (0, 1/3] \cup [2/3, 1]$, los autores determinan el valor exacto del radio de inyectividad:
    • Para $\beta \in (0, 1/3]$: $inj = \sqrt{2\beta}\pi$.
    • Para $\beta \in [2/3, 1]$: $inj = \pi$.
  • Para $\beta \in (1/3, 2/3)$, el valor exacto no se determina debido a la falta de una expresión explícita para el radio conjugado. Sin embargo, el artículo proporciona límites superiores e inferiores ajustados, mostrando que el valor difiere en un máximo del 3% del valor conjeturado en [33].
  • Para $\beta > 1$, el radio de inyectividad está acotado superiormente por $\pi$.

• Casos Especiales:

  • Para $p=1$ (esfera unitaria), $inj = \pi$.
  • Para $p=n$ (grupo ortogonal $O(n)$), el radio conjugado se calcula explícitamente como $\sqrt{8\beta}\pi$ para $n=3$ y $\sqrt{4\beta}\pi$ para $n>3$, lo que conduce a un radio de inyectividad de $\sqrt{2\beta}\pi$.

Significancia
El artículo afirma resolver el problema de determinar el radio de inyectividad para la variedad de Stiefel bajo una amplia gama de las métricas $\beta$, un parámetro crucial para los algoritmos de optimización en variedades. Al establecer la longitud exacta de los bucles geodésicos más cortos para todo $\beta$ y cerrar el vacío en los límites de la curvatura seccional, los autores proporcionan los datos geométricos necesarios para aplicar el teorema de Klingenberg.

El trabajo es modesto respecto al intervalo $\beta \in (1/3, 2/3)$, reconociendo que el radio de inyectividad exacto sigue siendo una conjetura en este rango, apoyada por experimentos numéricos en la literatura previa; sin embargo, el artículo fortalece la conjetura para la métrica canónica ($\beta=1/2$) al proporcionar expresiones explícitas para los campos de Jacobi, confirmando que los puntos conjugados identificados son, de hecho, los primeros puntos conjugados. Los resultados proporcionan una base rigurosa para asegurar la minimidad de las geodésicas en aplicaciones de optimización basada en variedades y ajuste de datos.