Technische Samenvatting: Stabiliteit van Optimale Transportafbeeldingen en Tweede Variatie van de 2-Monge–Kantorovich-Afstand

1. Probleemstelling

Het artikel behandelt de kwantitatieve stabiliteit van optimale transportafbeeldingen en potentialen in de context van de kwadratische kost ($p=2$). Specifiek onderzoekt het hoe de optimale transportafbeelding $T = \nabla \phi$ en het bijbehorende Brenier-potentieel $\phi$ veranderen wanneer de bron-dichtheid $f$ en de doeldichtheid $g$ worden verstoord.

Hoewel kwalitatieve stabiliteit (convergentie in maat of uniforme convergentie onder sterke aannames) goed is vastgesteld via compactheidsargumenten, richten de auteurs zich op het afleiden van kwantitatieve schattingen. Deze schattingen begrenzen de afstand tussen optimale afbeeldingen/potentialen die corresponderen met verschillende dichtheden, uitgedrukt in termen van de afstand tussen de dichtheden zelf (gemeten in $L^p$, $L^2$ of 2-Monge–Kantorovich-afstand). Het werk richt zich specifiek op situaties waarin dichtheden niet-degeneraat zijn (afgezet van nul) en Hölder-regulariteit bezitten, op uniform convex domeinen.

2. Methodologie

De kern van de technische strategie bestaat uit het lineariseren van de Monge–Ampère-vergelijking langs gladde interpolaties van de data. De methodologie verloopt als volgt:

1. Interpolatie: Gegeven twee paren dichtheden $(f_0, g_0)$ en $(f_1, g_1)$, construeren de auteurs een gladde kromme van dichtheden $\{(f_t, g_t)\}_{t \in [0,1]}$.
   • Voor Stelling 1.1 gebruiken ze 2-Monge–Kantorovich-geodeten (Benamou–Brenier-formulering).
   • Voor Stellingen 1.2 en 1.3 gebruiken ze lineaire interpolatie in $L^p$-ruimten.
2. Differentiatie van Potentialen: Ze definiëren een kromme van Brenier-potentialen $\phi_t$ die corresponderen met $(f_t, g_t)$. Het verschil $\phi_1 - \phi_0$ wordt uitgedrukt als het integraal van zijn tijdsafgeleide: $\phi_1 - \phi_0 = \int_0^1 \partial_t \phi_t \, dt$.
3. Gelineariseerde Vergelijking: De afgeleide $\xi_t = \partial_t \phi_t$ wordt gekarakteriseerd als de unieke oplossing van een gelineariseerde Monge–Ampère-vergelijking afgeleid in divergentievorm:
   $$-\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{in } \Omega,$$
   onderworpen aan Neumann-randvoorwaarden.
4. Regulariteit en Ellipticiteit: Een cruciale stap is het vaststellen dat de coëfficiëntenmatrix $(D^2 \phi_t)^{-1}$ uniform elliptisch blijft en dat de oplossing $\xi_t$ voldoende regulariteit behoudt (bijv. $C^{1,\alpha}$) langs het interpolatiepad. Dit berust op:
   • Lemma 1.7: Een compactheidsargument dat bewijst dat als dichtheden $C^{0,\beta}$ en niet-degeneraat zijn, het Brenier-potentieel $C^{2,\alpha}$ is tot aan de rand.
   • Stelling van de Implíciete Functie: Gebruikt om de gladheid van de afbeelding $t \mapsto \phi_t$ te bewijzen in geschikte Banachruimten (Stelling 1.5).
5. Schatting: Door $\xi_t$ te schatten met behulp van elliptische regulariteitstheorie (Schauder-schattingen, $L^\infty$-grenzen) en te integreren over $t$, leiden de auteurs de uiteindelijke stabiliteitsgrenzen af.

3. Belangrijkste Bijdragen en Resultaten

A. Kwantitatieve Stabiliteitsschattingen

Het artikel stelt drie belangrijke stabiliteitsstellingen vast onder de aanname van niet-degeneraat, Hölder-continue dichtheden op uniform convex domeinen met gladde randen:

• Stelling 1.1 (Lipschitz $L^2$-stabiliteit):
  Voor Hölder-continue dichtheden bewijzen de auteurs Lipschitz-continuïteit van de optimale transportafbeeldingen in de $L^2$-norm met betrekking tot de 2-Monge–Kantorovich-afstand ($d_2$) tussen dichtheden:
  $$\|\nabla \phi_1 - \nabla \phi_0\|_{L^2} \leq C [d_2(f_0, f_1) + d_2(g_0, g_1)].$$
  Dit verbetert eerdere resultaten (bijv. Berman [20]) die vaak dimensie-afhankelijke exponenten omvatten of sterkere regulariteit vereisten. De auteurs merken op dat dit Lipschitz-gedrag een specifiek kenmerk is van regelmatige, niet-degeneraat maatstaven.

• Stelling 1.2 (Stabiliteit zonder Continuïteitsaanname):
  Door de Hölder-continuïteitsaanname te verwijderen maar $L^2$-integreerbaarheid en niet-degeneratie te behouden, verkrijgen de auteurs een Lipschitz-stabiliteitsschatting voor de Brenier-potentialen in de $L^2$-norm:
  $$\|\phi_0 - \phi_1\|_{L^2} \leq C (\|f_0 - f_1\|_{L^2} + \|g_0 - g_1\|_{L^2}).$$
  Ze leiden ook een overeenkomstige schatting voor de gradienten af met een machtsexponent van $1/3$. De auteurs vermoeden dat de $L^2$-grens voor potentialen mogelijk kan worden verbeterd tot een $L^\infty$-grens als een specifieke Sobolev-ongelijkheid voor de gelineariseerde operator kan worden bewezen.

• Stelling 1.3 (Uniforme $C^{1,\alpha}$-stabiliteit):
  In een uniforme metriek-omgeving bewijzen de auteurs dat voor $C^{0,\beta}$-dichtheden de optimale afbeeldingen stabiel zijn in de $C^{1,\alpha}$-norm met betrekking tot $L^p$-afstanden van de dichtheden (waarbij $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}).$$
  In tegenstelling tot eerdere interieure schattingen (bijv. [28, 29]), is dit resultaat globaal, van toepassing op verstoringen rond niet-identiteitsafbeeldingen, en lijdt het niet aan dimensie-afnemende exponenten. De auteurs benadrukken dat de ondergrensvoorwaarde ($g_i \geq a > 0$) fundamenteel is; zonder deze bestaat er zelfs in 1D geen Hölder-stabiliteitsexponent $\eta < 1$.

B. Tweede Variatie van de Afstand

• Stelling 1.4: De auteurs leiden een expliciete formule voor de tweede variatie van de gekwadrateerde 2-Monge–Kantorovich-afstand langs een verstoring van dichtheden af.
  $$\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,$$
  waarbij $\xi$ de oplossing is van de gelineariseerde vergelijking. Deze formule kwantificeert de lineaire convexiteit van de afstand en breidt formele berekeningen uit die eerder in [30] werden gevonden (die beperkt waren tot verstoringen rond de identiteit) naar algemene regelmatige optimale transportafbeeldingen.

C. Technisch Hulpmiddel: Tijdsdifferentieerbaarheid

• Stelling 1.5 & Corollarium 1.6: Een aanzienlijk deel van het werk is gewijd aan het bewijzen dat het Brenier-potentieel $\phi_t$ differentieerbaar is met betrekking tot de interpolatieparameter $t$. De auteurs karakteriseren de afgeleide $\partial_t \phi_t$ als de zwakke oplossing van de gelineariseerde Monge–Ampère-vergelijking in divergentievorm. Dit resultaat is van onafhankelijk belang voor de analyse van tijd-afhankelijke optimale transportproblemen.

4. Betekenis en Claims

De auteurs positioneren hun werk als een significante vooruitgang in de kwantitatieve stabiliteit van optimaal transport, een momenteel actief onderzoeksgebied.

• Verbetering van Exponenten: Het artikel claimt Lipschitz- (exponent 1) stabiliteitsschattingen te bieden in situaties waar eerdere literatuur vaak sub-lineaire exponenten opleverde (bijv. $1/2$, $1/5$, of dimensie-afhankelijke afnemende exponenten).
• Globaal versus Interieur: Stelling 1.3 wordt geïntroduceerd als een globale schatting, in contrast met eerdere resultaten die vaak beperkt waren tot interieure schattingen of specifieke geometrische situaties (bijv. oppervlakken van codimensie één).
• Expliciete Tweede Variatie: De afleiding van de formule voor de tweede variatie wordt gepresenteerd als een hulpmiddel dat de convexiteit van de Monge–Kantorovich-afstand expliciet kwantificeert, en het theoretische begrip uitbreidt voorbij de identiteitsafbeelding.
• Beperkingen: De auteurs zijn bescheiden over de reikwijdte en stellen expliciet dat hun resultaten afhankelijk zijn van niet-degeneraat dichtheden (afgezet van nul) en Hölder-regulariteit. Ze erkennen dat voor discrete maatstaven of verdwijnende dichtheden de Lipschitz-natuur van de schattingen ineenstort, en dat de optimale exponenten lager kunnen zijn (vermoedelijk $1/2$ voor discrete doelen). Ze claimen niet de stabiliteit voor algemene ruwe maatstaven of singuliere doelen op te lossen.

Kort samengevat biedt het artikel een rigoureus kader voor het lineariseren van de Monge–Ampère-vergelijking in divergentievorm om scherpe, kwantitatieve stabiliteitsgrenzen te verkrijgen voor optimale transportafbeeldingen en potentialen onder regulariteits- en niet-degeneratie-aannames, terwijl het ook een nieuwe formule biedt voor de tweede variatie van de transportafstand.