Resumo Técnico: Estabilidade de Mapas de Transporte Ótimo e Segunda Variação da Distância 2-Monge–Kantorovich

1. Formulação do Problema

O artigo aborda a estabilidade quantitativa de mapas e potenciais de transporte ótimo no contexto do custo quadrático ($p=2$). Especificamente, investiga como o mapa de transporte ótimo $T = \nabla \phi$ e o potencial de Brenier associado $\phi$ mudam quando a densidade fonte $f$ e a densidade alvo $g$ são perturbadas.

Enquanto a estabilidade qualitativa (convergência em medida ou convergência uniforme sob fortes hipóteses) está bem estabelecida via argumentos de compacidade, os autores focam na derivação de estimativas quantitativas. Essas estimativas limitam a distância entre mapas/potenciais ótimos correspondentes a diferentes densidades em termos da distância entre as próprias densidades (medida em $L^p$, $L^2$ ou distância 2-Monge–Kantorovich). O trabalho visa especificamente cenários onde as densidades são não degeneradas (limitadas inferiormente de zero) e possuem regularidade de Hölder, em domínios uniformemente convexos.

2. Metodologia

A estratégia técnica central envolve linearizar a equação de Monge–Ampère ao longo de interpolações suaves dos dados. A metodologia procede da seguinte forma:

1. Interpolação: Dadas duas pares de densidades $(f_0, g_0)$ e $(f_1, g_1)$, os autores constroem uma curva suave de densidades $\{(f_t, g_t)\}_{t \in [0,1]}$.
   • Para o Teorema 1.1, utilizam geodésicas 2-Monge–Kantorovich (formulação de Benamou–Brenier).
   • Para os Teoremas 1.2 e 1.3, utilizam interpolação linear em espaços $L^p$.
2. Diferenciação de Potenciais: Eles definem uma curva de potenciais de Brenier $\phi_t$ correspondente a $(f_t, g_t)$. A diferença $\phi_1 - \phi_0$ é expressa como a integral de sua derivada temporal: $\phi_1 - \phi_0 = \int_0^1 \partial_t \phi_t \, dt$.
3. Equação Linearizada: A derivada $\xi_t = \partial_t \phi_t$ é caracterizada como a solução única de uma equação de Monge–Ampère linearizada derivada na forma de divergência:
   $$-\text{div}\left[ f_t (D^2 \phi_t)^{-1} \nabla \xi_t \right] = \left( \frac{\partial_t g_t(\nabla \phi_t)}{g_t(\nabla \phi_t)} - \frac{\partial_t f_t}{f_t} \right) f_t \quad \text{em } \Omega,$$
   sujeita a condições de contorno de Neumann.
4. Regularidade e Elipticidade: Um passo crucial é estabelecer que a matriz de coeficientes $(D^2 \phi_t)^{-1}$ permanece uniformemente elíptica e que a solução $\xi_t$ mantém regularidade suficiente (e.g., $C^{1,\alpha}$) ao longo do caminho de interpolação. Isso depende de:
   • Lema 1.7: Um argumento de compacidade provando que, se as densidades são $C^{0,\beta}$ e não degeneradas, o potencial de Brenier é $C^{2,\alpha}$ até a fronteira.
   • Teorema da Função Implícita: Usado para provar a suavidade do mapa $t \mapsto \phi_t$ em espaços de Banach apropriados (Teorema 1.5).
5. Estimação: Ao estimar $\xi_t$ usando a teoria de regularidade elíptica (estimativas de Schauder, limites $L^\infty$) e integrando sobre $t$, os autores derivam as limites finais de estabilidade.

3. Contribuições e Resultados Chave

A. Estimativas de Estabilidade Quantitativa

O artigo estabelece três teoremas principais de estabilidade sob a hipótese de densidades não degeneradas e contínuas de Hölder em domínios uniformemente convexos com fronteiras suaves:

• Teorema 1.1 (Estabilidade Lipschitz $L^2$):
  Para densidades contínuas de Hölder, os autores provam a continuidade Lipschitz dos mapas de transporte ótimo na norma $L^2$ em relação à distância 2-Monge–Kantorovich ($d_2$) entre densidades:
  $$\|\nabla \phi_1 - \nabla \phi_0\|_{L^2} \leq C [d_2(f_0, f_1) + d_2(g_0, g_1)].$$
  Isso melhora resultados anteriores (e.g., Berman [20]) que frequentemente envolviam expoentes dependentes da dimensão ou exigiam regularidade mais forte. Os autores observam que esse comportamento Lipschitz é uma característica especial de medidas regulares e não degeneradas.

• Teorema 1.2 (Estabilidade sem Hipótese de Continuidade):
  Removendo a hipótese de continuidade de Hölder, mas mantendo a integrabilidade $L^2$ e a não degenerescência, os autores obtêm uma estimativa de estabilidade Lipschitz para os potenciais de Brenier na norma $L^2$:
  $$\|\phi_0 - \phi_1\|_{L^2} \leq C (\|f_0 - f_1\|_{L^2} + \|g_0 - g_1\|_{L^2}).$$
  Eles também derivam uma estimativa correspondente para os gradientes com um expoente de potência $1/3$. Os autores conjecturam que o limite $L^2$ nos potenciais poderia potencialmente ser melhorado para um limite $L^\infty$ se uma desigualdade de Sobolev específica para o operador linearizado pudesse ser provada.

• Teorema 1.3 (Estabilidade Uniforme $C^{1,\alpha}$):
  Em um cenário de métrica uniforme, os autores provam que, para densidades $C^{0,\beta}$, os mapas ótimos são estáveis na norma $C^{1,\alpha}$ em relação às distâncias $L^p$ das densidades (onde $p = d/(1-\alpha)$):
  $$\|\phi_1 - \phi_0\|_{C^{1,\alpha}} \leq C (\|f_1 - f_0\|_{L^p} + \|g_1 - g_0\|_{L^p}).$$
  Diferentemente de estimativas interiores anteriores (e.g., [28, 29]), este resultado é global, aplica-se a perturbações ao redor de mapas não identidade e não sofre de expoentes que decaem com a dimensão. Os autores enfatizam que a condição de limite inferior ($g_i \geq a > 0$) é fundamental; sem ela, mesmo em 1D, não existe expoente de estabilidade de Hölder $\eta < 1$.

B. Segunda Variação da Distância

• Teorema 1.4: Os autores derivam uma fórmula explícita para a segunda variação do quadrado da distância 2-Monge–Kantorovich ao longo de uma perturbação de densidades.
  $$\frac{1}{2} \frac{d^2}{dt^2}\bigg|_{t=0} [d_2(f_t, g_t)^2] = \int_\Omega \langle (D^2 \phi_0)^{-1} \nabla \xi, \nabla \xi \rangle f \, dx,$$
  onde $\xi$ é a solução da equação linearizada. Esta fórmula quantifica a convexidade linear da distância e estende cálculos formais encontrados anteriormente em [30] (que eram restritos a perturbações ao redor da identidade) para mapas de transporte ótimo regulares gerais.

C. Ferramenta Técnica: Diferenciabilidade Temporal

• Teorema 1.5 & Corolário 1.6: Uma parte significativa do trabalho é dedicada a provar que o potencial de Brenier $\phi_t$ é diferenciável em relação ao parâmetro de interpolação $t$. Os autores caracterizam a derivada $\partial_t \phi_t$ como a solução fraca da equação de Monge–Ampère linearizada na forma de divergência. Este resultado é de interesse independente para a análise de problemas de transporte ótimo dependentes do tempo.

4. Significado e Alegações

Os autores posicionam seu trabalho como um avanço significativo na estabilidade quantitativa do transporte ótimo, uma área de pesquisa atualmente ativa.

• Melhoria nos Expoentes: O artigo alega fornecer estimativas de estabilidade Lipschitz (expoente 1) em cenários onde a literatura anterior frequentemente produzia expoentes sub-lineares (e.g., $1/2$, $1/5$, ou expoentes decrescentes dependentes da dimensão).
• Global vs. Interior: O Teorema 1.3 é destacado como uma estimativa global, contrastando com resultados anteriores que frequentemente eram restritos a estimativas interiores ou cenários geométricos específicos (e.g., superfícies de codimensão um).
• Segunda Variação Explícita: A derivação da fórmula de segunda variação é apresentada como uma ferramenta que quantifica explicitamente a convexidade da distância de Monge–Kantorovich, estendendo a compreensão teórica além do mapa identidade.
• Limitações: Os autores são modestos quanto ao escopo, afirmando explicitamente que seus resultados dependem de densidades não degeneradas (limitadas inferiormente de zero) e regularidade de Hölder. Eles reconhecem que, para medidas discretas ou densidades que se anulam, a natureza Lipschitz das estimativas se quebra, e os expoentes ótimos podem ser menores (conjecturados como $1/2$ para alvos discretos). Eles não alegam resolver a estabilidade para medidas gerais irregulares ou alvos singulares.

Em resumo, o artigo fornece um quadro rigoroso para linearizar a equação de Monge–Ampère na forma de divergência para obter limites de estabilidade quantitativa agudos para mapas e potenciais de transporte ótimo sob hipóteses de regularidade e não degenerescência, fornecendo também uma nova fórmula para a segunda variação da distância de transporte.