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

Mappe di trasporto ottimale instabili

Informazioni Fondamentali

  • 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

Riassunto

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 ρ\rho in Rd\mathbb{R}^d. In primo luogo, si dimostra che l'illimitatezza della densità può causare instabilità: viene costruita una densità sorgente sulla sfera unitaria in Rd\mathbb{R}^d 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.

Contesto di Ricerca e Motivazione

Descrizione del Problema

Una questione centrale nella teoria del trasporto ottimale è lo studio della stabilità della mappa di trasporto ottimale TμT_\mu (dalla misura sorgente fissa ρ\rho alla misura bersaglio μ\mu) rispetto alle perturbazioni della misura bersaglio. Specificamente, si studia se esistono costanti C,α>0C, \alpha > 0 e p1p \geq 1 tali che per qualsiasi misura di probabilità μ,ν\mu, \nu su un insieme compatto YRdY \subset \mathbb{R}^d, valga: TμTνL2(ρ)CWp(μ,ν)α\|\|T_\mu - T_\nu\|\|_{L^2(\rho)} \leq C W_p(\mu, \nu)^\alpha

Importanza

L'importanza di questo problema si manifesta in molteplici aspetti:

  1. Validità del quadro del trasporto ottimale linearizzato: la mappa μTμ\mu \mapsto T_\mu fornisce un'immersione dallo spazio di Wasserstein (P(Y),W2)(P(Y), W_2) allo spazio di Hilbert L2(ρ,Rd)L^2(\rho, \mathbb{R}^d)
  2. Analisi numerica: quando la misura bersaglio μ\mu è disponibile solo attraverso un'approssimazione μ^\hat{\mu}, è necessario garantire che il Tμ^T_{\hat{\mu}} calcolato sia vicino al vero TμT_\mu
  3. Stima statistica: nel trasporto ottimale statistico, è necessario stimare la mappa di trasporto ottimale da campioni

Limitazioni della Ricerca Esistente

Sebbene siano state stabilite disuguaglianze di stabilità quantificata sotto varie ipotesi sulla misura sorgente ρ\rho, 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à.

Contributi Fondamentali

  1. 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
  2. Identificazione di due meccanismi di instabilità:
    • Esplosione superpolinomiale della densità
    • Avvicinamento a configurazioni in cui l'unicità del piano ottimale viene meno
  3. Fornitura di costruzioni esplicite di controesempi: attraverso costruzioni matematiche precise si dimostra il fenomeno di instabilità
  4. Rivelazione dei limiti di regolarità della mappa μTμ\mu \mapsto T_\mu: si dimostra che in alcuni casi questa mappa non è 12\frac{1}{2}-Hölder continua

Dettagli Metodologici

Definizione del Compito

Dati ρ,μP2(Rd)\rho, \mu \in P_2(\mathbb{R}^d), dove ρ\rho è assolutamente continua, il teorema di Brenier garantisce l'esistenza di un'unica mappa di trasporto ottimale TμL2(ρ)T_\mu \in L^2(\rho) da ρ\rho a μ\mu. Si studia la stabilità di TμT_\mu rispetto alle perturbazioni di μ\mu.

Metodo di Costruzione del Teorema 1.1

Costruzione della Misura Sorgente

Sia A=(1,0,,0)A = (1, 0, \ldots, 0), A=(1,0,,0)A' = (-1, 0, \ldots, 0), E={A,A}E = \{A, A'\}. Si definisce la funzione: f(r)=rdmin(1,(logr)2)f(r) = r^{-d} \min(1, (\log r)^{-2})

Si costruisce la densità: ρ(x)=c0f(dist(x,E))\rho(x) = c_0 f(\text{dist}(x, E)) dove c0>0c_0 > 0 è tale che ρ\rho sia una densità di probabilità sulla sfera unitaria BRd(0,1)B_{\mathbb{R}^d}(0,1).

Sequenza di Misure Bersaglio

Si considerano Bθ=(Rsinθ,Rcosθ,0,,0)B_\theta = (R\sin\theta, R\cos\theta, 0, \ldots, 0) e Bθ=(Rsinθ,Rcosθ,0,,0)B'_\theta = (-R\sin\theta, -R\cos\theta, 0, \ldots, 0), e si definisce: μθ=12(δBθ+δBθ)\mu_\theta = \frac{1}{2}(\delta_{B_\theta} + \delta_{B'_\theta})

Stime Chiave

Attraverso l'analisi geometrica si dimostra: 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)|}

Mentre Wp(μ0,μθ)Rθ/2W_p(\mu_0, \mu_\theta) \sim R\theta/2, pertanto la disuguaglianza di stabilità fallisce per qualsiasi α>0\alpha > 0.

Metodo di Costruzione del Teorema 1.2

Costruzione dell'Insieme di Supporto

Si costruisce l'unione di infinite coppie di parallelepipedi rettangolari: X=i=1+SiX = \bigcup_{i=1}^{+\infty} S_i dove 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).

Scelta dei Parametri

Si scelgono sequenze che soddisfano:

  • ri=c02ir_i = c_0 2^{-i} (decadimento superpolinomiale)
  • i=wi=c0c1i2\ell_i = w_i = c_0 c_1 i^{-2}
  • Distanze tra celle molto maggiori delle distanze all'interno delle celle

Dimostrazione dell'Instabilità

Attraverso l'analisi di perturbazioni specifiche della mappa di trasporto ottimale, si dimostra: 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}

Quando α>p2(p+1)\alpha > \frac{p}{2(p+1)}, il membro destro tende all'infinito.

Configurazione Sperimentale

Verifica Teorica

Questo articolo è un lavoro puramente teorico, che verifica i risultati attraverso dimostrazioni matematiche rigorose:

  1. Analisi geometrica: analisi dettagliata della struttura geometrica della mappa di trasporto ottimale
  2. Stime di misura: calcolo esatto della misura di varie regioni
  3. Calcolo delle distanze: stima accurata della distanza di Wasserstein e della distanza L2L^2

Verifica della Costruzione

Attraverso i seguenti passaggi si verifica la correttezza della costruzione:

  1. Dimostrazione che la densità costruita è effettivamente una densità di probabilità
  2. Verifica dell'espressione esplicita della mappa di trasporto ottimale
  3. Conferma che la configurazione geometrica soddisfa le condizioni di separazione richieste

Risultati Sperimentali

Risultati Principali

Risultati del Teorema 1.1

Per la densità sorgente costruita ρ\rho (che esplode in modo superpolinomiale in due punti del bordo della sfera unitaria), si dimostra che per qualsiasi sfera Y=BRd(0,R)Y = B_{\mathbb{R}^d}(0,R), qualsiasi C,α>0C, \alpha > 0 e p1p \geq 1, la disuguaglianza: μ,ν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 non è valida.

Risultati del Teorema 1.2

Per la misura di probabilità uniforme ρ\rho su un insieme aperto limitato, si dimostra che per qualsiasi C>0C > 0, p1p \geq 1 e α>p2(p+1)\alpha > \frac{p}{2(p+1)} (in particolare, α=12\alpha = \frac{1}{2}), la disuguaglianza di stabilità fallisce.

Scoperte Tecniche

  1. Criticità dell'esplosione della densità: l'esplosione polinomiale rd+δr^{-d+\delta} mantiene ancora la stabilità, ma l'esplosione superpolinomiale conduce a instabilità completa
  2. Limitazioni della regolarità Hölder: la mappa μTμ\mu \mapsto T_\mu in alcuni casi non è 12\frac{1}{2}-Hölder continua, e nemmeno 13\frac{1}{3}-Hölder continua
  3. Presagi della non-unicità: la perdita di unicità del piano ottimale è preceduta da un deterioramento della stabilità

Lavori Correlati

Risultati Positivi di Stabilità

  • Gigli (2011): teoria del limite superiore della 12\frac{1}{2}-continuità Hölder della mappa μTμ\mu \mapsto T_\mu
  • 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.

Ricerca sull'Instabilità

Precedentemente esistevano solo esempi di instabilità del potenziale di Kantorovich; questo articolo costruisce per la prima volta mappe di trasporto instabili.

Conclusioni e Discussione

Conclusioni Principali

  1. Le mappe di trasporto ottimale possono essere estremamente instabili, e qualsiasi disuguaglianza di stabilità quantificata può fallire
  2. L'esplosione superpolinomiale della densità è una condizione sufficiente per l'instabilità
  3. Anche per densità uniformi, l'avvicinamento a configurazioni non-uniche produce instabilità

Limitazioni

  1. Il Teorema 1.1 richiede densità illimitate, che potrebbero non essere realistiche nelle applicazioni pratiche
  2. La costruzione del Teorema 1.2 richiede infinite componenti connesse, con una struttura geometrica complessa
  3. Non è stato determinato se densità limitate possono raggiungere la forte instabilità del Teorema 1.1

Direzioni Future

  1. Ricerca sulla massima instabilità nel caso di densità limitate
  2. Sviluppo di metodi di stima statistica per misure sorgente instabili
  3. Esplorazione dei legami più profondi tra instabilità e struttura geometrica

Valutazione Approfondita

Punti di Forza

  1. Rottura teorica: prima dimostrazione dell'instabilità delle mappe di trasporto ottimale, colmando un vuoto teorico
  2. Costruzione ingegnosa: realizzazione dell'obiettivo teorico attraverso costruzioni geometriche e scelte di parametri astute
  3. Analisi rigorosa: dimostrazione matematica completa, gestione appropriata dei dettagli tecnici
  4. Intuizioni profonde: rivelazione dei legami intrinseci tra stabilità, unicità e regolarità della densità

Carenze

  1. Limitazioni di praticità: gli esempi costruiti potrebbero non essere comuni nelle applicazioni pratiche
  2. Problemi aperti: alcuni problemi naturali (come la forte instabilità per densità limitate) rimangono irrisolti
  3. Assenza di verifica numerica: come lavoro teorico, mancano esperimenti numerici di verifica

Impatto

  1. Valore teorico: fornisce importanti controesempi e intuizioni alla teoria del trasporto ottimale
  2. Avvertimento applicativo: avverte i praticanti dei rischi di instabilità in certe configurazioni
  3. Ispirazione per la ricerca: fornisce nuove direzioni e strumenti per ricerche successive

Scenari Applicabili

  1. Ricerca approfondita della teoria del trasporto ottimale
  2. Analisi teorica del quadro del trasporto ottimale linearizzato
  3. Analisi di stabilità nel trasporto ottimale statistico
  4. Valutazione della stabilità degli algoritmi numerici di trasporto ottimale

Bibliografia

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.