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.
- 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
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 ρ in Rd. Zunächst wird bewiesen, dass Unbegrenztheit der Dichte zu Instabilität führen kann: Es wird eine Quelldichte auf der Einheitskugel in Rd 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.
Ein Kernproblem der optimalen Transporttheorie ist die Untersuchung der Stabilität der optimalen Transportabbildung Tμ (vom festen Quellmaß ρ zum Zielmaß μ) bezüglich Störungen des Zielmaßes. Konkret wird untersucht, ob Konstanten C,α>0 und p≥1 existieren, so dass für beliebige Wahrscheinlichkeitsmaße μ,ν auf einer kompakten Menge Y⊂Rd gilt:
∥∥Tμ−Tν∥∥L2(ρ)≤CWp(μ,ν)α
Die Bedeutung dieses Problems zeigt sich in mehreren Aspekten:
- Rechtfertigung des linearisierten Optimalransport-Rahmens: Die Abbildung μ↦Tμ bietet eine Einbettung vom Wasserstein-Raum (P(Y),W2) in den Hilbert-Raum L2(ρ,Rd)
- Numerische Analyse: Wenn das Zielmaß μ nur durch eine Approximation μ^ verfügbar ist, muss sichergestellt werden, dass das berechnete Tμ^ dem wahren Tμ nahe kommt
- Statistische Schätzung: In der statistischen optimalen Transporttheorie müssen optimale Transportabbildungen aus Stichproben geschätzt werden
Obwohl unter verschiedenen Annahmen über das Quellmaß ρ 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.
- 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
- Identifikation zweier Instabilitätsmechanismen:
- Überpolynomielle Explosion der Dichte
- Annäherung an Konfigurationen, in denen der optimale Plan nicht eindeutig ist
- Bereitstellung konkreter Gegenbeispielkonstruktionen: Durch präzise mathematische Konstruktionen wird das Instabilitätsphänomen demonstriert
- Offenlegung von Regularitätsbeschränkungen der Abbildung μ↦Tμ: Es wird bewiesen, dass diese Abbildung in bestimmten Fällen nicht 21-Hölder-stetig ist
Gegeben ρ,μ∈P2(Rd), wobei ρ absolut stetig ist, garantiert das Brenier-Theorem die Existenz einer eindeutigen optimalen Transportabbildung Tμ∈L2(ρ) von ρ zu μ. Die Stabilität von Tμ bezüglich Störungen von μ wird untersucht.
Setze A=(1,0,…,0), A′=(−1,0,…,0), E={A,A′}. Definiere die Funktion:
f(r)=r−dmin(1,(logr)−2)
Konstruiere die Dichte:
ρ(x)=c0f(dist(x,E))
wobei c0>0 so gewählt ist, dass ρ eine Wahrscheinlichkeitsdichte auf der Einheitskugel BRd(0,1) ist.
Betrachte Bθ=(Rsinθ,Rcosθ,0,…,0) und Bθ′=(−Rsinθ,−Rcosθ,0,…,0), definiere:
μθ=21(δBθ+δBθ′)
Durch geometrische Analyse wird bewiesen:
∥∥Tμθ−Tμ0∥∥L2(ρ)2≥∣log(θ/4)∣c0c1σd−1R2
Da Wp(μ0,μθ)∼Rθ/2, scheitert die Stabilitätsungleichung für jedes α>0.
Konstruiere eine unendliche Vereinigung von Paaren rechteckiger Parallelepipede:
X=⋃i=1+∞Si
wobei Si=T+(Ai+,ℓi,ri)∪T−(Ai−,ℓi,ri).
Wähle Sequenzen, die erfüllen:
- ri=c02−i (überpolynomieller Zerfall)
- ℓi=wi=c0c1i−2
- Abstände zwischen Zellen sind viel größer als Abstände innerhalb von Zellen
Durch Analyse spezifischer Störungen der optimalen Transportabbildung wird bewiesen:
Wp(μ,νi)2α∥∥Tμ−Tνi∥∥L2(ρ)2≥Cwi3−p2αri1−p2α−2α
Wenn α>2(p+1)p, strebt die rechte Seite gegen Unendlich.
Dieses Papier ist eine rein theoretische Arbeit, die Ergebnisse durch strenge mathematische Beweise verifiziert:
- Geometrische Analyse: Detaillierte Analyse der geometrischen Struktur optimaler Transportabbildungen
- Maßabschätzungen: Präzise Berechnung von Maßen verschiedener Regionen
- Distanzberechnungen: Genaue Abschätzung von Wasserstein-Distanzen und L2-Distanzen
Die Korrektheit der Konstruktion wird durch folgende Schritte verifiziert:
- Beweis, dass die konstruierte Dichte tatsächlich eine Wahrscheinlichkeitsdichte ist
- Verifikation der expliziten Ausdrücke der optimalen Transportabbildung
- Bestätigung, dass die geometrische Konfiguration die erforderlichen Trennungsbedingungen erfüllt
Für die konstruierte Quelldichte ρ (die an zwei Grenzpunkten der Einheitskugel überpolynomiell explodiert), wird bewiesen, dass für jeden Ball Y=BRd(0,R), beliebige C,α>0 und p≥1 die Ungleichung:
∀μ,ν∈P(Y),∥∥Tμ−Tν∥∥L2(ρ)≤CWp(μ,ν)α
nicht erfüllt ist.
Für das gleichmäßige Wahrscheinlichkeitsmaß ρ auf beschränkten offenen Mengen wird bewiesen, dass für beliebige C>0, p≥1 und α>2(p+1)p (insbesondere α=21) die Stabilitätsungleichung scheitert.
- Kritikalität der Dichteexplosion: Polynomielle Explosion r−d+δ bewahrt noch Stabilität, aber überpolynomielle Explosion führt zu vollständiger Instabilität
- Hölder-Regularitätsbeschränkungen: Die Abbildung μ↦Tμ ist in bestimmten Fällen nicht 21-Hölder-stetig, nicht einmal 31-Hölder-stetig
- Vorboten der Nicht-Eindeutigkeit: Stabilitätsverschlechterung tritt auf, bevor die Eindeutigkeit des optimalen Plans verloren geht
- Gigli (2011): Theorie der 21-Hölder-Stetigkeit der Abbildung μ↦Tμ
- 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.
Bisher gab es nur Beispiele für Instabilität von Kantorovich-Potentialfunktionen; dieses Papier konstruiert erstmals instabile Transportabbildungen.
- Optimale Transportabbildungen können hochgradig instabil sein, wobei jede quantifizierte Stabilitätsungleichung scheitern kann
- Überpolynomielle Explosion der Dichte ist eine hinreichende Bedingung für Instabilität
- Selbst für gleichmäßige Dichten tritt Instabilität auf, wenn man sich Nicht-Eindeutigkeitskonfigurationen nähert
- Theorem 1.1 erfordert unbegrenzte Dichte, was in praktischen Anwendungen möglicherweise unrealistisch ist
- Die Konstruktion in Theorem 1.2 erfordert unendlich viele zusammenhängende Komponenten mit komplexer geometrischer Struktur
- Es konnte nicht bestimmt werden, ob beschränkte Dichten die starke Instabilität von Theorem 1.1 erreichen können
- Untersuchung der stärksten Instabilität bei beschränkter Dichte
- Entwicklung statistischer Schätzmethoden für instabile Quellmaße
- Erforschung der tieferen Verbindung zwischen Instabilität und geometrischer Struktur
- Theoretischer Durchbruch: Erstmaliger Beweis der Instabilität optimaler Transportabbildungen, schließt eine theoretische Lücke
- Raffinierte Konstruktion: Durch geschickte geometrische Konstruktion und Parameterwahl werden theoretische Ziele erreicht
- Strenge Analyse: Mathematische Beweise sind vollständig und technische Details sind angemessen behandelt
- Tiefe Einsichten: Offenlegung der inneren Verbindung zwischen Stabilität, Eindeutigkeit und Dichteregularität
- Praktische Einschränkungen: Die konstruierten Beispiele treten in praktischen Anwendungen möglicherweise selten auf
- Offene Probleme: Bestimmte natürliche Fragen (wie starke Instabilität bei beschränkter Dichte) bleiben ungelöst
- Fehlende numerische Verifikation: Als theoretische Arbeit fehlen numerische Experimente zur Verifikation
- Theoretischer Wert: Bietet wichtige Gegenbeispiele und Einsichten für die Optimalransporttheorie
- Anwendungswarnung: Warnt Praktiker vor Instabilitätsrisiken in bestimmten Konfigurationen
- Forschungsanregung: Bietet neue Richtungen und Werkzeuge für nachfolgende Forschung
- Tiefgreifende Forschung der Optimalransporttheorie
- Theoretische Analyse des linearisierten Optimalransport-Rahmens
- Stabilitätsanalyse in der statistischen optimalen Transporttheorie
- Stabilitätsbewertung numerischer Optimalransport-Algorithmen
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.