2025-11-22T16:07:16.564646

Unstable optimal transport maps

Letrouit
The stability of optimal transport maps with respect to perturbations of the marginals is a question of interest for several reasons, ranging from the justification of the linearized optimal transport framework to numerical analysis and statistics. Under various assumptions on the source measure, it is known that optimal transport maps are stable with respect to variations of the target measure. In this note, we focus on the mechanisms that can, on the contrary, lead to instability. We identify two of them, which we illustrate through examples of absolutely continuous source measures $ρ$ in $\mathbb{R}^d$ for which optimal transport maps are less stable, or even very unstable. We first show that instability may arise from the unboundedness of the density: we exhibit a source density on the unit ball of $\mathbb{R}^d$ which blows up superpolynomially at two points of the boundary and for which optimal transport maps are highly unstable. Then we prove that even for uniform densities on bounded open sets, optimal transport maps can be rather unstable close enough to configurations where uniqueness of optimal plans is lost.
academic

Instabile optimale Transportabbildungen

Grundinformationen

  • Papier-ID: 2510.13265
  • Titel: Unstable optimal transport maps
  • Autor: Cyril Letrouit (Université Paris-Saclay, CNRS)
  • Klassifizierung: math.OC (Mathematische Optimierung und Steuerung)
  • Veröffentlichungsdatum: 15. Oktober 2025
  • Papier-Link: https://arxiv.org/abs/2510.13265

Zusammenfassung

Die Stabilität optimaler Transportabbildungen bezüglich Störungen von Randverteilungen ist ein wichtiges Problem mit Auswirkungen auf die Rechtfertigung des linearisierten Optimalransport-Rahmens, die numerische Analyse und die Statistik. Unter verschiedenen Annahmen über die Quellmaße ist bekannt, dass optimale Transportabbildungen bezüglich Änderungen des Zielmaßes stabil sind. Dieses Papier konzentriert sich auf den gegenteiligen Mechanismus – die Faktoren, die zu Instabilität führen. Die Autoren identifizieren zwei Mechanismen und veranschaulichen die Instabilität und sogar hochgradige Instabilität optimaler Transportabbildungen durch Beispiele absolut stetiger Quellmaße ρ\rho in Rd\mathbb{R}^d. Zunächst wird bewiesen, dass Unbegrenztheit der Dichte zu Instabilität führen kann: Es wird eine Quelldichte auf der Einheitskugel in Rd\mathbb{R}^d konstruiert, die an zwei Grenzpunkten überpolynomiell explodiert, wobei die entsprechende optimale Transportabbildung hochgradig instabil ist. Dann wird bewiesen, dass selbst für gleichmäßige Dichten auf beschränkten offenen Mengen die optimale Transportabbildung erheblich instabil sein kann, wenn man sich einer Konfiguration nähert, in der die Eindeutigkeit des optimalen Plans verloren geht.

Forschungshintergrund und Motivation

Problembeschreibung

Ein Kernproblem der optimalen Transporttheorie ist die Untersuchung der Stabilität der optimalen Transportabbildung TμT_\mu (vom festen Quellmaß ρ\rho zum Zielmaß μ\mu) bezüglich Störungen des Zielmaßes. Konkret wird untersucht, ob Konstanten C,α>0C, \alpha > 0 und p1p \geq 1 existieren, so dass für beliebige Wahrscheinlichkeitsmaße μ,ν\mu, \nu auf einer kompakten Menge YRdY \subset \mathbb{R}^d gilt: TμTνL2(ρ)CWp(μ,ν)α\|\|T_\mu - T_\nu\|\|_{L^2(\rho)} \leq C W_p(\mu, \nu)^\alpha

Bedeutung

Die Bedeutung dieses Problems zeigt sich in mehreren Aspekten:

  1. Rechtfertigung des linearisierten Optimalransport-Rahmens: Die Abbildung μTμ\mu \mapsto T_\mu bietet eine Einbettung vom Wasserstein-Raum (P(Y),W2)(P(Y), W_2) in den Hilbert-Raum L2(ρ,Rd)L^2(\rho, \mathbb{R}^d)
  2. Numerische Analyse: Wenn das Zielmaß μ\mu nur durch eine Approximation μ^\hat{\mu} verfügbar ist, muss sichergestellt werden, dass das berechnete Tμ^T_{\hat{\mu}} dem wahren TμT_\mu nahe kommt
  3. Statistische Schätzung: In der statistischen optimalen Transporttheorie müssen optimale Transportabbildungen aus Stichproben geschätzt werden

Einschränkungen bisheriger Forschung

Obwohl unter verschiedenen Annahmen über das Quellmaß ρ\rho quantifizierte Stabilitätsungleichungen etabliert wurden, wurde bisher nie bewiesen, dass optimale Transportabbildungen instabil sein können. Bisherige Ergebnisse konzentrierten sich hauptsächlich auf positive Stabilitätsergebnisse und es fehlte ein tieferes Verständnis der Instabilitätsmechanismen.

Kernbeiträge

  1. Erste Konstruktion hochgradig instabiler optimaler Transportabbildungen: Es wird bewiesen, dass absolut stetige Quellmaße existieren, für die keine Form von quantifizierter Stabilitätsungleichung gilt
  2. Identifikation zweier Instabilitätsmechanismen:
    • Überpolynomielle Explosion der Dichte
    • Annäherung an Konfigurationen, in denen der optimale Plan nicht eindeutig ist
  3. Bereitstellung konkreter Gegenbeispielkonstruktionen: Durch präzise mathematische Konstruktionen wird das Instabilitätsphänomen demonstriert
  4. Offenlegung von Regularitätsbeschränkungen der Abbildung μTμ\mu \mapsto T_\mu: Es wird bewiesen, dass diese Abbildung in bestimmten Fällen nicht 12\frac{1}{2}-Hölder-stetig ist

Methodische Erläuterung

Aufgabendefinition

Gegeben ρ,μP2(Rd)\rho, \mu \in P_2(\mathbb{R}^d), wobei ρ\rho absolut stetig ist, garantiert das Brenier-Theorem die Existenz einer eindeutigen optimalen Transportabbildung TμL2(ρ)T_\mu \in L^2(\rho) von ρ\rho zu μ\mu. Die Stabilität von TμT_\mu bezüglich Störungen von μ\mu wird untersucht.

Konstruktionsmethode von Theorem 1.1

Konstruktion des Quellmaßes

Setze A=(1,0,,0)A = (1, 0, \ldots, 0), A=(1,0,,0)A' = (-1, 0, \ldots, 0), E={A,A}E = \{A, A'\}. Definiere die Funktion: f(r)=rdmin(1,(logr)2)f(r) = r^{-d} \min(1, (\log r)^{-2})

Konstruiere die Dichte: ρ(x)=c0f(dist(x,E))\rho(x) = c_0 f(\text{dist}(x, E)) wobei c0>0c_0 > 0 so gewählt ist, dass ρ\rho eine Wahrscheinlichkeitsdichte auf der Einheitskugel BRd(0,1)B_{\mathbb{R}^d}(0,1) ist.

Zielmaß-Sequenz

Betrachte Bθ=(Rsinθ,Rcosθ,0,,0)B_\theta = (R\sin\theta, R\cos\theta, 0, \ldots, 0) und Bθ=(Rsinθ,Rcosθ,0,,0)B'_\theta = (-R\sin\theta, -R\cos\theta, 0, \ldots, 0), definiere: μθ=12(δBθ+δBθ)\mu_\theta = \frac{1}{2}(\delta_{B_\theta} + \delta_{B'_\theta})

Schlüsselabschätzung

Durch geometrische Analyse wird bewiesen: TμθTμ0L2(ρ)2c0c1σd1R2log(θ/4)\|\|T_{\mu_\theta} - T_{\mu_0}\|\|^2_{L^2(\rho)} \geq \frac{c_0 c_1 \sigma_{d-1} R^2}{|\log(\theta/4)|}

Da Wp(μ0,μθ)Rθ/2W_p(\mu_0, \mu_\theta) \sim R\theta/2, scheitert die Stabilitätsungleichung für jedes α>0\alpha > 0.

Konstruktionsmethode von Theorem 1.2

Konstruktion der Trägermenge

Konstruiere eine unendliche Vereinigung von Paaren rechteckiger Parallelepipede: X=i=1+SiX = \bigcup_{i=1}^{+\infty} S_i wobei Si=T+(Ai+,i,ri)T(Ai,i,ri)S_i = T^+(A_i^+, \ell_i, r_i) \cup T^-(A_i^-, \ell_i, r_i).

Parameterauswahl

Wähle Sequenzen, die erfüllen:

  • ri=c02ir_i = c_0 2^{-i} (überpolynomieller Zerfall)
  • i=wi=c0c1i2\ell_i = w_i = c_0 c_1 i^{-2}
  • Abstände zwischen Zellen sind viel größer als Abstände innerhalb von Zellen

Beweis der Instabilität

Durch Analyse spezifischer Störungen der optimalen Transportabbildung wird bewiesen: TμTνiL2(ρ)2Wp(μ,νi)2αCwi32αpri12αp2α\frac{\|\|T_\mu - T_{\nu_i}\|\|^2_{L^2(\rho)}}{W_p(\mu, \nu_i)^{2\alpha}} \geq C w_i^{3-\frac{2\alpha}{p}} r_i^{1-\frac{2\alpha}{p}-2\alpha}

Wenn α>p2(p+1)\alpha > \frac{p}{2(p+1)}, strebt die rechte Seite gegen Unendlich.

Experimentelle Einrichtung

Theoretische Verifikation

Dieses Papier ist eine rein theoretische Arbeit, die Ergebnisse durch strenge mathematische Beweise verifiziert:

  1. Geometrische Analyse: Detaillierte Analyse der geometrischen Struktur optimaler Transportabbildungen
  2. Maßabschätzungen: Präzise Berechnung von Maßen verschiedener Regionen
  3. Distanzberechnungen: Genaue Abschätzung von Wasserstein-Distanzen und L2L^2-Distanzen

Verifikation der Konstruktion

Die Korrektheit der Konstruktion wird durch folgende Schritte verifiziert:

  1. Beweis, dass die konstruierte Dichte tatsächlich eine Wahrscheinlichkeitsdichte ist
  2. Verifikation der expliziten Ausdrücke der optimalen Transportabbildung
  3. Bestätigung, dass die geometrische Konfiguration die erforderlichen Trennungsbedingungen erfüllt

Experimentelle Ergebnisse

Hauptergebnisse

Ergebnisse von Theorem 1.1

Für die konstruierte Quelldichte ρ\rho (die an zwei Grenzpunkten der Einheitskugel überpolynomiell explodiert), wird bewiesen, dass für jeden Ball Y=BRd(0,R)Y = B_{\mathbb{R}^d}(0,R), beliebige C,α>0C, \alpha > 0 und p1p \geq 1 die Ungleichung: μ,νP(Y),TμTνL2(ρ)CWp(μ,ν)α\forall \mu, \nu \in P(Y), \|\|T_\mu - T_\nu\|\|_{L^2(\rho)} \leq C W_p(\mu, \nu)^\alpha nicht erfüllt ist.

Ergebnisse von Theorem 1.2

Für das gleichmäßige Wahrscheinlichkeitsmaß ρ\rho auf beschränkten offenen Mengen wird bewiesen, dass für beliebige C>0C > 0, p1p \geq 1 und α>p2(p+1)\alpha > \frac{p}{2(p+1)} (insbesondere α=12\alpha = \frac{1}{2}) die Stabilitätsungleichung scheitert.

Technische Erkenntnisse

  1. Kritikalität der Dichteexplosion: Polynomielle Explosion rd+δr^{-d+\delta} bewahrt noch Stabilität, aber überpolynomielle Explosion führt zu vollständiger Instabilität
  2. Hölder-Regularitätsbeschränkungen: Die Abbildung μTμ\mu \mapsto T_\mu ist in bestimmten Fällen nicht 12\frac{1}{2}-Hölder-stetig, nicht einmal 13\frac{1}{3}-Hölder-stetig
  3. Vorboten der Nicht-Eindeutigkeit: Stabilitätsverschlechterung tritt auf, bevor die Eindeutigkeit des optimalen Plans verloren geht

Verwandte Arbeiten

Positive Stabilitätsergebnisse

  • Gigli (2011): Theorie der 12\frac{1}{2}-Hölder-Stetigkeit der Abbildung μTμ\mu \mapsto T_\mu
  • Letrouit & Mérigot (2024): Stabilitätsungleichungen für beschränkte Dichten auf John-Domänen
  • Weitere Arbeiten: Stabilitätsergebnisse für log-konkave Dichten, polynomiell explodierende Dichten, etc.

Instabilitätsforschung

Bisher gab es nur Beispiele für Instabilität von Kantorovich-Potentialfunktionen; dieses Papier konstruiert erstmals instabile Transportabbildungen.

Schlussfolgerungen und Diskussion

Hauptschlussfolgerungen

  1. Optimale Transportabbildungen können hochgradig instabil sein, wobei jede quantifizierte Stabilitätsungleichung scheitern kann
  2. Überpolynomielle Explosion der Dichte ist eine hinreichende Bedingung für Instabilität
  3. Selbst für gleichmäßige Dichten tritt Instabilität auf, wenn man sich Nicht-Eindeutigkeitskonfigurationen nähert

Einschränkungen

  1. Theorem 1.1 erfordert unbegrenzte Dichte, was in praktischen Anwendungen möglicherweise unrealistisch ist
  2. Die Konstruktion in Theorem 1.2 erfordert unendlich viele zusammenhängende Komponenten mit komplexer geometrischer Struktur
  3. Es konnte nicht bestimmt werden, ob beschränkte Dichten die starke Instabilität von Theorem 1.1 erreichen können

Zukünftige Richtungen

  1. Untersuchung der stärksten Instabilität bei beschränkter Dichte
  2. Entwicklung statistischer Schätzmethoden für instabile Quellmaße
  3. Erforschung der tieferen Verbindung zwischen Instabilität und geometrischer Struktur

Tiefgreifende Bewertung

Stärken

  1. Theoretischer Durchbruch: Erstmaliger Beweis der Instabilität optimaler Transportabbildungen, schließt eine theoretische Lücke
  2. Raffinierte Konstruktion: Durch geschickte geometrische Konstruktion und Parameterwahl werden theoretische Ziele erreicht
  3. Strenge Analyse: Mathematische Beweise sind vollständig und technische Details sind angemessen behandelt
  4. Tiefe Einsichten: Offenlegung der inneren Verbindung zwischen Stabilität, Eindeutigkeit und Dichteregularität

Schwächen

  1. Praktische Einschränkungen: Die konstruierten Beispiele treten in praktischen Anwendungen möglicherweise selten auf
  2. Offene Probleme: Bestimmte natürliche Fragen (wie starke Instabilität bei beschränkter Dichte) bleiben ungelöst
  3. Fehlende numerische Verifikation: Als theoretische Arbeit fehlen numerische Experimente zur Verifikation

Einfluss

  1. Theoretischer Wert: Bietet wichtige Gegenbeispiele und Einsichten für die Optimalransporttheorie
  2. Anwendungswarnung: Warnt Praktiker vor Instabilitätsrisiken in bestimmten Konfigurationen
  3. Forschungsanregung: Bietet neue Richtungen und Werkzeuge für nachfolgende Forschung

Anwendbare Szenarien

  1. Tiefgreifende Forschung der Optimalransporttheorie
  2. Theoretische Analyse des linearisierten Optimalransport-Rahmens
  3. Stabilitätsanalyse in der statistischen optimalen Transporttheorie
  4. Stabilitätsbewertung numerischer Optimalransport-Algorithmen

Literaturverzeichnis

Das Papier zitiert wichtige Literatur aus dem Optimalransportbereich, einschließlich grundlegender Arbeiten von Brenier, Regularitätsergebnisse von Gigli und aktuelle Forschung zu quantifizierten Stabilitätsungleichungen, die eine solide theoretische Grundlage für die Forschung bieten.