Cartes de transport optimal instables

Informations de base

• ID de l'article: 2510.13265
• Titre: Unstable optimal transport maps
• Auteur: Cyril Letrouit (Université Paris-Saclay, CNRS)
• Classification: math.OC (Optimisation mathématique et contrôle)
• Date de publication: 15 octobre 2025
• Lien de l'article: https://arxiv.org/abs/2510.13265

Résumé

La stabilité des cartes de transport optimal par rapport aux perturbations des distributions marginales est une question importante qui touche à la validité du cadre du transport optimal linéarisé, à l'analyse numérique et aux statistiques. Sous diverses hypothèses sur la mesure source, on sait que la variation des cartes de transport optimal par rapport à la mesure cible est stable. Cet article se concentre sur le mécanisme inverse — les facteurs qui conduisent à l'instabilité. Les auteurs identifient deux mécanismes et illustrent l'instabilité, voire l'instabilité extrême, des cartes de transport optimal par des exemples de mesures sources absolument continues $\rho$ dans $\mathbb{R}^d$. Premièrement, il est démontré que l'absence de borne de la densité peut conduire à l'instabilité : une densité source est construite sur la boule unité de $\mathbb{R}^d$ qui explose de manière surpolynomiale en deux points de la frontière, et la carte de transport optimal correspondante est extrêmement instable. Ensuite, il est démontré que même pour une densité uniforme sur un ensemble ouvert borné, la carte de transport optimal peut être considérablement instable lorsqu'on s'approche d'une configuration où l'unicité du plan optimal est perdue.

Contexte et motivation de la recherche

Description du problème

Une question centrale de la théorie du transport optimal est d'étudier la stabilité de la carte de transport optimal $T_\mu$ (de la mesure source fixe $\rho$ vers la mesure cible $\mu$) par rapport aux perturbations de la mesure cible. Concrètement, il s'agit de savoir s'il existe des constantes $C, \alpha > 0$ et $p \geq 1$ telles que pour toutes mesures de probabilité $\mu, \nu$ sur un ensemble compact $Y \subset \mathbb{R}^d$, on ait : $\|\|T_\mu - T_\nu\|\|_{L^2(\rho)} \leq C W_p(\mu, \nu)^\alpha$

Importance

L'importance de cette question se manifeste sous plusieurs aspects :

1. Validité du cadre du transport optimal linéarisé : l'application $\mu \mapsto T_\mu$ fournit un plongement de l'espace de Wasserstein $(P(Y), W_2)$ dans l'espace de Hilbert $L^2(\rho, \mathbb{R}^d)$
2. Analyse numérique : lorsque la mesure cible $\mu$ n'est accessible que par une approximation $\hat{\mu}$, il est nécessaire de garantir que le $T_{\hat{\mu}}$ calculé est proche du vrai $T_\mu$
3. Estimation statistique : dans le transport optimal statistique, il est nécessaire d'estimer la carte de transport optimal à partir d'échantillons

Limitations de la recherche existante

Bien que des inégalités de stabilité quantifiées aient été établies sous diverses hypothèses sur la mesure source $\rho$, on n'a jamais démontré auparavant que les cartes de transport optimal pouvaient être instables. Les résultats existants se concentrent principalement sur les résultats positifs de stabilité, manquant d'une compréhension approfondie des mécanismes d'instabilité.

Contributions principales

1. Première construction de cartes de transport optimal extrêmement instables : démonstration de l'existence de mesures sources absolument continues pour lesquelles aucune inégalité de stabilité quantifiée ne peut être établie
2. Identification de deux mécanismes d'instabilité :
   • Explosion surpolynomiale de la densité
   • Approche d'une configuration où l'unicité du plan optimal est perdue
3. Fourniture de constructions de contre-exemples explicites : démonstration du phénomène d'instabilité par des constructions mathématiques précises
4. Révélation des limitations de régularité de l'application $\mu \mapsto T_\mu$ : démonstration que dans certains cas, cette application n'est pas $\frac{1}{2}$-Hölder continue

Détails méthodologiques

Définition de la tâche

Étant donné $\rho, \mu \in P_2(\mathbb{R}^d)$, où $\rho$ est absolument continue, le théorème de Brenier garantit l'existence d'une unique carte de transport optimal $T_\mu \in L^2(\rho)$ de $\rho$ vers $\mu$. On étudie la stabilité de $T_\mu$ par rapport aux perturbations de $\mu$.

Méthode de construction du Théorème 1.1

Construction de la mesure source

Soit $A = (1, 0, \ldots, 0)$, $A' = (-1, 0, \ldots, 0)$, $E = \{A, A'\}$. Définissons la fonction : $f(r) = r^{-d} \min(1, (\log r)^{-2})$

Construisons la densité : $\rho(x) = c_0 f(\text{dist}(x, E))$ où $c_0 > 0$ est choisi de sorte que $\rho$ soit une densité de probabilité sur la boule unité $B_{\mathbb{R}^d}(0,1)$.

Séquence de mesures cibles

Considérons $B_\theta = (R\sin\theta, R\cos\theta, 0, \ldots, 0)$ et $B'_\theta = (-R\sin\theta, -R\cos\theta, 0, \ldots, 0)$, définissons : $\mu_\theta = \frac{1}{2}(\delta_{B_\theta} + \delta_{B'_\theta})$

Estimations clés

Par analyse géométrique, on démontre : $\|\|T_{\mu_\theta} - T_{\mu_0}\|\|^2_{L^2(\rho)} \geq \frac{c_0 c_1 \sigma_{d-1} R^2}{|\log(\theta/4)|}$

Tandis que $W_p(\mu_0, \mu_\theta) \sim R\theta/2$, par conséquent l'inégalité de stabilité échoue pour tout $\alpha > 0$.

Méthode de construction du Théorème 1.2

Construction de l'ensemble support

Construction d'une union infinie de paires de parallélépipèdes rectangulaires : $X = \bigcup_{i=1}^{+\infty} S_i$ où $S_i = T^+(A_i^+, \ell_i, r_i) \cup T^-(A_i^-, \ell_i, r_i)$.

Choix des paramètres

Sélection de séquences satisfaisant :

• $r_i = c_0 2^{-i}$ (décroissance surpolynomiale)
• $\ell_i = w_i = c_0 c_1 i^{-2}$
• Distance inter-cellules bien supérieure à la distance intra-cellules

Démonstration de l'instabilité

Par analyse des variations de la carte de transport optimal sous perturbations spécifiques, on démontre : $\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}$

Lorsque $\alpha > \frac{p}{2(p+1)}$, le membre droit tend vers l'infini.

Configuration expérimentale

Vérification théorique

Cet article est un travail purement théorique, dont les résultats sont vérifiés par démonstration mathématique rigoureuse :

1. Analyse géométrique : analyse détaillée de la structure géométrique des cartes de transport optimal
2. Estimation de mesures : calcul précis des mesures dans diverses régions
3. Calcul de distances : estimation précise de la distance de Wasserstein et de la distance $L^2$

Vérification de la construction

Vérification de la correction de la construction par les étapes suivantes :

1. Démonstration que la densité construite est effectivement une densité de probabilité
2. Vérification des expressions explicites des cartes de transport optimal
3. Confirmation que la configuration géométrique satisfait les conditions de séparation requises

Résultats expérimentaux

Résultats principaux

Résultats du Théorème 1.1

Pour la densité source construite $\rho$ (explosion surpolynomiale en deux points de la frontière de la boule unité), il est démontré que pour toute boule $Y = B_{\mathbb{R}^d}(0,R)$, tous $C, \alpha > 0$ et $p \geq 1$, l'inégalité : $\forall \mu, \nu \in P(Y), \|\|T_\mu - T_\nu\|\|_{L^2(\rho)} \leq C W_p(\mu, \nu)^\alpha$ ne peut pas être satisfaite.

Résultats du Théorème 1.2

Pour une mesure de probabilité uniforme $\rho$ sur un ensemble ouvert borné, il est démontré que pour tous $C > 0$, $p \geq 1$ et $\alpha > \frac{p}{2(p+1)}$ (en particulier, $\alpha = \frac{1}{2}$), l'inégalité de stabilité échoue.

Découvertes techniques

1. Caractère critique de l'explosion de densité : l'explosion polynomiale $r^{-d+\delta}$ maintient la stabilité, mais l'explosion surpolynomiale entraîne une instabilité complète
2. Limitations de régularité Hölder : l'application $\mu \mapsto T_\mu$ n'est pas $\frac{1}{2}$-Hölder continue dans certains cas, voire pas $\frac{1}{3}$-Hölder continue
3. Présage de non-unicité : la dégradation de la stabilité précède la perte d'unicité du plan optimal

Travaux connexes

Résultats positifs de stabilité

• Gigli (2011) : théorie de la continuité $\frac{1}{2}$-Hölder de l'application $\mu \mapsto T_\mu$
• Letrouit & Mérigot (2024) : inégalités de stabilité pour densités bornées sur domaines de John
• Autres travaux : résultats de stabilité pour densités log-concaves, densités avec explosion polynomiale, etc.

Recherche sur l'instabilité

Seuls des exemples d'instabilité du potentiel de Kantorovich existaient auparavant ; cet article est le premier à construire des cartes de transport instables.

Conclusion et discussion

Conclusions principales

1. Les cartes de transport optimal peuvent être extrêmement instables, et toute inégalité de stabilité quantifiée peut échouer
2. L'explosion surpolynomiale de la densité est une condition suffisante pour l'instabilité
3. Même pour une densité uniforme, l'instabilité apparaît lorsqu'on s'approche d'une configuration de non-unicité

Limitations

1. Le Théorème 1.1 nécessite une densité non bornée, ce qui peut ne pas être réaliste dans les applications pratiques
2. La construction du Théorème 1.2 nécessite une infinité de composantes connexes, avec une structure géométrique complexe
3. Impossibilité de déterminer si une densité bornée peut atteindre la forte instabilité du Théorème 1.1

Directions futures

1. Étude de l'instabilité maximale pour les densités bornées
2. Développement de méthodes d'estimation statistique adaptées aux mesures sources instables
3. Exploration des liens profonds entre instabilité et structure géométrique

Évaluation approfondie

Avantages

1. Percée théorique : première démonstration de l'instabilité des cartes de transport optimal, comblant une lacune théorique
2. Construction ingénieuse : réalisation de l'objectif théorique par des constructions géométriques astucieuses et des choix de paramètres
3. Analyse rigoureuse : démonstrations mathématiques complètes avec traitement approprié des détails techniques
4. Intuitions profondes : révélation des connexions intrinsèques entre stabilité, unicité et régularité de densité

Insuffisances

1. Limitations pratiques : les exemples construits peuvent ne pas être courants dans les applications réelles
2. Problèmes ouverts : certaines questions naturelles (comme l'instabilité forte pour densités bornées) restent non résolues
3. Absence de vérification numérique : en tant que travail théorique, il manque des expériences numériques de vérification

Impact

1. Valeur théorique : fournit des contre-exemples importants et des intuitions pour la théorie du transport optimal
2. Avertissement pratique : alerte les praticiens sur les risques d'instabilité dans certaines configurations
3. Inspiration pour la recherche : fournit de nouvelles directions et outils pour les recherches ultérieures

Domaines d'application

1. Recherche approfondie en théorie du transport optimal
2. Analyse théorique du cadre du transport optimal linéarisé
3. Analyse de stabilité dans le transport optimal statistique
4. Évaluation de la stabilité des algorithmes numériques de transport optimal

Références

L'article cite les travaux importants du domaine du transport optimal, incluant les travaux fondateurs de Brenier, les résultats de régularité de Gigli, ainsi que les recherches récentes sur la stabilité quantifiée, fournissant une base théorique solide pour la recherche.