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

Resumo Técnico: Laços Geodésicos Curtos, Curvatura Seccional e Raio de Injetividade da Variedade de Stiefel

Enunciado do Problema
A variedade de Stiefel, denotada por $St(n, p)$, é uma variedade matricial fundamental que aparece em otimização, métodos numéricos para equações diferenciais e estatística. Um problema central de geometria Riemanniana nesta variedade é determinar o raio de injetividade, que fornece uma condição suficiente para a minimalidade de geodésicas conectando dois pontos. O raio de injetividade, $inj(M)$, é definido como o maior raio $r$ tal que toda geodésica de comprimento $\le r$ é a curva única mais curta entre seus pontos finais.

Para uma variedade Riemanniana completa, o raio de injetividade é determinado por dois fatores concorrentes: metade do comprimento do laço geodésico não trivial mais curto ($\ell/2$) e o raio conjugado. Este último é limitado inferiormente por $\pi/\sqrt{K}$, onde $K$ é um limite superior global na curvatura seccional. Embora o raio de injetividade seja conhecido para métricas específicas (por exemplo, a métrica euclidiana e a métrica canônica), ele permanece um problema aberto para a família geral de métricas de um parâmetro introduzida por Hüper et al. (2021), parametrizada por $\beta > 0$. Especificamente, o comprimento exato dos laços geodésicos mais curtos e o limite superior agudo da curvatura seccional não foram totalmente estabelecidos para todo $\beta$, impedindo o cálculo exato do raio de injetividade em todo o espaço de parâmetros.

Metodologia
Os autores abordam este problema combinando análise de geometria diferencial com técnicas de análise matricial. A metodologia é dividida em três componentes principais:

1. Análise de Comprimentos de Laços Geodésicos ($\ell_\beta$):

   • Caso $\beta \in (0, 2]$: Os autores utilizam a observação de que as geodésicas na variedade de Stiefel, quando vistas como curvas no espaço euclidiano ambiente $\mathbb{R}^{n \times p}$, possuem curvaturas de Frenet constantes. Ao derivar limites para o comprimento euclidiano dessas curvas e traduzi-los para a métrica $\beta$, eles estabelecem limites inferiores para o comprimento do laço.
   • Caso $\beta > 2$: A abordagem de curvatura de Frenet é insuficiente. Em vez disso, os autores empregam uma prova de análise matricial. Eles analisam a estrutura da inversa da exponencial de matriz e os grupos de invariância de matrizes ortogonais. Ao decompor as matrizes de velocidade inicial em subblocos correspondentes a propriedades espectrais específicas, eles reduzem o problema à análise de laços geodésicos com estruturas bloco-diagonais específicas, provando que o comprimento não pode ser menor que $2\pi$.

2. Limites de Curvatura Seccional ($K_\beta$):
   Os autores derivam uma fórmula explícita para a curvatura seccional da variedade de Stiefel sob a métrica $\beta$. Eles realizam estimativas termo a termo na expressão da curvatura usando as normas dos blocos de matrizes ($A$ e $B$) que definem os vetores tangentes. Esta análise identifica uma lacuna na literatura existente para $\beta \in (2/3, 1/\sqrt{2})$, a qual os autores fecham ao provar que a curvatura seccional é limitada superiormente por 1 neste intervalo.

3. Síntese via Teorema de Klingenberg:
   Utilizando os comprimentos de laço calculados $\ell_\beta$ e os limites de curvatura $K_\beta$, os autores aplicam o teorema de Klingenber, que afirma que $inj(St_\beta(n, p)) = \min\{ \ell_\beta/2, \text{conj}_{I_{n \times p}}(St_\beta(n, p)) \}$. Eles combinam isso com o limite inferior do raio conjugado ($\pi/\sqrt{K_\beta}$) para determinar o raio de injetividade exato ou limites estreitos para vários intervalos de $\beta$.

Principais Contribuições e Resultados

• Comprimento Exato dos Laços Geodésicos Mais Curtos: O artigo prova que para todo $\beta > 0$, o comprimento $\ell_\beta$ do laço geodésico não trivial mais curto em $St_\beta(n, p)$ (para $2 \le p \le n-1$) é dado por:
  $$\ell_\beta = \min\{ \sqrt{2\beta}, 1 \} 2\pi$$
  Este resultado unifica e estende achados anteriores para a métrica canônica ($\beta=1/2$) e a métrica euclidiana ($\beta=1$), onde o comprimento era conhecido como $2\pi$. Para $\beta > 2$, o comprimento é mostrado como sendo exatamente $2\pi$.

• Limites Agudos de Curvatura Seccional: Os autores fornecem um limite superior agudo $K_\beta$ sobre a curvatura seccional para o intervalo $\beta \in (2/3, 1/\sqrt{2})$, preenchendo uma lacuna na literatura existente. O limite é:
  $$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: O artigo observa que para $\beta > 1$, o limite de curvatura é 1, consistente com o caso euclidiano).

• Determinação do Raio de Injetividade:

  • Para $\beta \in (0, 1/3] \cup [2/3, 1]$, os autores determinam o valor exato do raio de injetividade:
    • 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)$, o valor exato não é determinado devido à falta de uma expressão explícita para o raio conjugado. No entanto, o artigo fornece limites superiores e inferiores estreitos, mostrando que o valor difere em no máximo 3% do valor conjecturado em [33].
  • Para $\beta > 1$, o raio de injetividade é limitado superiormente por $\pi$.

• Casos Especiais:

  • Para $p=1$ (esfera unitária), $inj = \pi$.
  • Para $p=n$ (grupo ortogonal $O(n)$), o raio conjugado é explicitamente calculado como $\sqrt{8\beta}\pi$ para $n=3$ e $\sqrt{4\beta}\pi$ para $n>3$, levando a um raio de injetividade de $\sqrt{2\beta}\pi$.

Significância
O artigo afirma resolver o problema de determinar o raio de injetividade para a variedade de Stiefel sob uma ampla gama de métricas $\beta$, um parâmetro crucial para algoritmos de otimização em variedades. Ao estabelecer o comprimento exato dos laços geodésicos mais curtos para todo $\beta$ e fechar a lacuna nos limites de curvatura seccional, os autores fornecem os dados geométricos necessários para aplicar o teorema de Klingenberg.

O trabalho é modesto em relação ao intervalo $\beta \in (1/3, 2/3)$, reconhecendo que o raio de injetividade exato permanece uma conjectura nesse intervalo, apoiada por experimentos numéricos em literatura prévia. No entanto, o artigo fortalece a conjectura para a métrica canônica ($\beta=1/2$) ao fornecer expressões explícitas para campos de Jacobi, confirmando que os pontos conjugados identificados são de fato os primeiros pontos conjugados. Os resultados fornecem uma base rigorosa para garantir a minimalidade de geodésicas em aplicações de otimização em variedades e ajuste de dados.