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.
- ID Articolo: 2510.13265
- Titolo: Unstable optimal transport maps
- Autore: Cyril Letrouit (Université Paris-Saclay, CNRS)
- Classificazione: math.OC (Ottimizzazione Matematica e Controllo)
- Data di Pubblicazione: 15 ottobre 2025
- Link Articolo: https://arxiv.org/abs/2510.13265
La stabilità delle mappe di trasporto ottimale rispetto alle perturbazioni delle distribuzioni marginali è una questione importante che coinvolge la validità del quadro del trasporto ottimale linearizzato, l'analisi numerica e la statistica. Sotto varie ipotesi sulla misura sorgente, è noto che la variazione delle mappe di trasporto ottimale rispetto ai cambiamenti della misura bersaglio è stabile. Questo articolo si concentra sul meccanismo opposto — i fattori che conducono all'instabilità. Gli autori identificano due meccanismi e illustrano l'instabilità e persino l'instabilità estrema delle mappe di trasporto ottimale attraverso esempi di misure sorgente assolutamente continue ρ in Rd. In primo luogo, si dimostra che l'illimitatezza della densità può causare instabilità: viene costruita una densità sorgente sulla sfera unitaria in Rd che esplode in modo superpolinomiale in due punti del bordo, e la corrispondente mappa di trasporto ottimale è estremamente instabile. Successivamente, si dimostra che anche per densità uniformi su insiemi aperti limitati, la mappa di trasporto ottimale può essere piuttosto instabile quando ci si avvicina a configurazioni in cui l'unicità del piano ottimale viene meno.
Una questione centrale nella teoria del trasporto ottimale è lo studio della stabilità della mappa di trasporto ottimale Tμ (dalla misura sorgente fissa ρ alla misura bersaglio μ) rispetto alle perturbazioni della misura bersaglio. Specificamente, si studia se esistono costanti C,α>0 e p≥1 tali che per qualsiasi misura di probabilità μ,ν su un insieme compatto Y⊂Rd, valga:
∥∥Tμ−Tν∥∥L2(ρ)≤CWp(μ,ν)α
L'importanza di questo problema si manifesta in molteplici aspetti:
- Validità del quadro del trasporto ottimale linearizzato: la mappa μ↦Tμ fornisce un'immersione dallo spazio di Wasserstein (P(Y),W2) allo spazio di Hilbert L2(ρ,Rd)
- Analisi numerica: quando la misura bersaglio μ è disponibile solo attraverso un'approssimazione μ^, è necessario garantire che il Tμ^ calcolato sia vicino al vero Tμ
- Stima statistica: nel trasporto ottimale statistico, è necessario stimare la mappa di trasporto ottimale da campioni
Sebbene siano state stabilite disuguaglianze di stabilità quantificata sotto varie ipotesi sulla misura sorgente ρ, non è mai stato provato che le mappe di trasporto ottimale possono essere instabili. I risultati esistenti si concentrano principalmente su risultati positivi di stabilità, mancando di una comprensione profonda dei meccanismi di instabilità.
- Prima costruzione di mappe di trasporto ottimale estremamente instabili: si dimostra l'esistenza di misure sorgente assolutamente continue per cui nessuna forma di disuguaglianza di stabilità quantificata è valida
- Identificazione di due meccanismi di instabilità:
- Esplosione superpolinomiale della densità
- Avvicinamento a configurazioni in cui l'unicità del piano ottimale viene meno
- Fornitura di costruzioni esplicite di controesempi: attraverso costruzioni matematiche precise si dimostra il fenomeno di instabilità
- Rivelazione dei limiti di regolarità della mappa μ↦Tμ: si dimostra che in alcuni casi questa mappa non è 21-Hölder continua
Dati ρ,μ∈P2(Rd), dove ρ è assolutamente continua, il teorema di Brenier garantisce l'esistenza di un'unica mappa di trasporto ottimale Tμ∈L2(ρ) da ρ a μ. Si studia la stabilità di Tμ rispetto alle perturbazioni di μ.
Sia A=(1,0,…,0), A′=(−1,0,…,0), E={A,A′}. Si definisce la funzione:
f(r)=r−dmin(1,(logr)−2)
Si costruisce la densità:
ρ(x)=c0f(dist(x,E))
dove c0>0 è tale che ρ sia una densità di probabilità sulla sfera unitaria BRd(0,1).
Si considerano Bθ=(Rsinθ,Rcosθ,0,…,0) e Bθ′=(−Rsinθ,−Rcosθ,0,…,0), e si definisce:
μθ=21(δBθ+δBθ′)
Attraverso l'analisi geometrica si dimostra:
∥∥Tμθ−Tμ0∥∥L2(ρ)2≥∣log(θ/4)∣c0c1σd−1R2
Mentre Wp(μ0,μθ)∼Rθ/2, pertanto la disuguaglianza di stabilità fallisce per qualsiasi α>0.
Si costruisce l'unione di infinite coppie di parallelepipedi rettangolari:
X=⋃i=1+∞Si
dove Si=T+(Ai+,ℓi,ri)∪T−(Ai−,ℓi,ri).
Si scelgono sequenze che soddisfano:
- ri=c02−i (decadimento superpolinomiale)
- ℓi=wi=c0c1i−2
- Distanze tra celle molto maggiori delle distanze all'interno delle celle
Attraverso l'analisi di perturbazioni specifiche della mappa di trasporto ottimale, si dimostra:
Wp(μ,νi)2α∥∥Tμ−Tνi∥∥L2(ρ)2≥Cwi3−p2αri1−p2α−2α
Quando α>2(p+1)p, il membro destro tende all'infinito.
Questo articolo è un lavoro puramente teorico, che verifica i risultati attraverso dimostrazioni matematiche rigorose:
- Analisi geometrica: analisi dettagliata della struttura geometrica della mappa di trasporto ottimale
- Stime di misura: calcolo esatto della misura di varie regioni
- Calcolo delle distanze: stima accurata della distanza di Wasserstein e della distanza L2
Attraverso i seguenti passaggi si verifica la correttezza della costruzione:
- Dimostrazione che la densità costruita è effettivamente una densità di probabilità
- Verifica dell'espressione esplicita della mappa di trasporto ottimale
- Conferma che la configurazione geometrica soddisfa le condizioni di separazione richieste
Per la densità sorgente costruita ρ (che esplode in modo superpolinomiale in due punti del bordo della sfera unitaria), si dimostra che per qualsiasi sfera Y=BRd(0,R), qualsiasi C,α>0 e p≥1, la disuguaglianza:
∀μ,ν∈P(Y),∥∥Tμ−Tν∥∥L2(ρ)≤CWp(μ,ν)α
non è valida.
Per la misura di probabilità uniforme ρ su un insieme aperto limitato, si dimostra che per qualsiasi C>0, p≥1 e α>2(p+1)p (in particolare, α=21), la disuguaglianza di stabilità fallisce.
- Criticità dell'esplosione della densità: l'esplosione polinomiale r−d+δ mantiene ancora la stabilità, ma l'esplosione superpolinomiale conduce a instabilità completa
- Limitazioni della regolarità Hölder: la mappa μ↦Tμ in alcuni casi non è 21-Hölder continua, e nemmeno 31-Hölder continua
- Presagi della non-unicità: la perdita di unicità del piano ottimale è preceduta da un deterioramento della stabilità
- Gigli (2011): teoria del limite superiore della 21-continuità Hölder della mappa μ↦Tμ
- Letrouit & Mérigot (2024): disuguaglianze di stabilità per densità limitate su domini di John
- Altri lavori: risultati di stabilità per densità log-concave, densità con esplosione polinomiale, ecc.
Precedentemente esistevano solo esempi di instabilità del potenziale di Kantorovich; questo articolo costruisce per la prima volta mappe di trasporto instabili.
- Le mappe di trasporto ottimale possono essere estremamente instabili, e qualsiasi disuguaglianza di stabilità quantificata può fallire
- L'esplosione superpolinomiale della densità è una condizione sufficiente per l'instabilità
- Anche per densità uniformi, l'avvicinamento a configurazioni non-uniche produce instabilità
- Il Teorema 1.1 richiede densità illimitate, che potrebbero non essere realistiche nelle applicazioni pratiche
- La costruzione del Teorema 1.2 richiede infinite componenti connesse, con una struttura geometrica complessa
- Non è stato determinato se densità limitate possono raggiungere la forte instabilità del Teorema 1.1
- Ricerca sulla massima instabilità nel caso di densità limitate
- Sviluppo di metodi di stima statistica per misure sorgente instabili
- Esplorazione dei legami più profondi tra instabilità e struttura geometrica
- Rottura teorica: prima dimostrazione dell'instabilità delle mappe di trasporto ottimale, colmando un vuoto teorico
- Costruzione ingegnosa: realizzazione dell'obiettivo teorico attraverso costruzioni geometriche e scelte di parametri astute
- Analisi rigorosa: dimostrazione matematica completa, gestione appropriata dei dettagli tecnici
- Intuizioni profonde: rivelazione dei legami intrinseci tra stabilità, unicità e regolarità della densità
- Limitazioni di praticità: gli esempi costruiti potrebbero non essere comuni nelle applicazioni pratiche
- Problemi aperti: alcuni problemi naturali (come la forte instabilità per densità limitate) rimangono irrisolti
- Assenza di verifica numerica: come lavoro teorico, mancano esperimenti numerici di verifica
- Valore teorico: fornisce importanti controesempi e intuizioni alla teoria del trasporto ottimale
- Avvertimento applicativo: avverte i praticanti dei rischi di instabilità in certe configurazioni
- Ispirazione per la ricerca: fornisce nuove direzioni e strumenti per ricerche successive
- Ricerca approfondita della teoria del trasporto ottimale
- Analisi teorica del quadro del trasporto ottimale linearizzato
- Analisi di stabilità nel trasporto ottimale statistico
- Valutazione della stabilità degli algoritmi numerici di trasporto ottimale
L'articolo cita importanti letteratura nel campo del trasporto ottimale, inclusi i lavori fondamentali di Brenier, i risultati di regolarità di Gigli, e ricerche recenti sulla stabilità quantificata, fornendo una base teorica solida per la ricerca.