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

不安定な最適輸送写像

基本情報

論文ID: 2510.13265
タイトル: Unstable optimal transport maps
著者: Cyril Letrouit (Université Paris-Saclay, CNRS)
分類: math.OC (数学的最適化と制御)
発表日: 2025年10月15日
論文リンク: https://arxiv.org/abs/2510.13265

要約

最適輸送写像の周辺分布摂動に対する安定性は、線形化最適輸送フレームワークの妥当性、数値解析、統計学など複数の分野に関わる重要な問題である。ソース測度の様々な仮定の下で、最適輸送写像のターゲット測度の変化に対する安定性は既知である。本論文は逆のメカニズム、すなわち不安定性をもたらす要因に焦点を当てる。著者は2つのメカニズムを特定し、 $\mathbb{R}^d$ における絶対連続ソース測度 $\rho$ の例を通じて、最適輸送写像の不安定性、さらには高度な不安定性を示す。まず、密度の無界性が不安定性をもたらす可能性があることを証明する。 $\mathbb{R}^d$ の単位球上のソース密度を構成し、この密度は境界の2点で超多項式的に爆発し、対応する最適輸送写像は高度に不安定である。次に、有界開集合上の均一密度に対してさえ、最適計画の一意性が失われる配置に近づくとき、最適輸送写像は相当な不安定性を示す可能性があることを証明する。

研究背景と動機

問題の記述

最適輸送理論の中心的な問題の1つは、固定ソース測度 $\rho$ からターゲット測度 $\mu$ への最適輸送写像 $T_\mu$ のターゲット測度摂動に対する安定性を研究することである。具体的には、定数 $C, \alpha > 0$ と $p \geq 1$ が存在し、コンパクト集合 $Y \subset \mathbb{R}^d$ 上の任意の確率測度 $\mu, \nu$ に対して以下が成立するかどうかを研究する： $\|\|T_\mu - T_\nu\|\|_{L^2(\rho)} \leq C W_p(\mu, \nu)^\alpha$

重要性

この問題の重要性は複数の側面に現れる：

線形化最適輸送フレームワークの妥当性：写像 $\mu \mapsto T_\mu$ はWasserstein空間 $(P(Y), W_2)$ からHilbert空間 $L^2(\rho, \mathbb{R}^d)$ への埋め込みを提供する
数値解析：ターゲット測度 $\mu$ が近似 $\hat{\mu}$ を通じてのみ得られる場合、計算された $T_{\hat{\mu}}$ が真の $T_\mu$ に近いことを保証する必要がある
統計推定：統計的最適輸送では、サンプルから最適輸送写像を推定する必要がある

既存研究の限界

ソース測度 $\rho$ の様々な仮定の下で定量的安定性不等式が確立されているが、最適輸送写像が不安定である可能性は以前証明されていない。既存の結果は主に安定性の肯定的な結果に焦点を当てており、不安定性メカニズムの深い理解が欠けている。

核心的貢献

高度に不安定な最適輸送写像の初めての構成：定量的安定性不等式のいかなる形式も成立しない絶対連続ソース測度が存在することを証明
2つの不安定性メカニズムの特定：
- 密度の超多項式的爆発
- 最適計画の非一意性に近づく配置
具体的な反例構成の提供：正確な数学的構成を通じて不安定性現象を示す
写像 $\mu \mapsto T_\mu$ の正則性限界の解明：特定の場合にこの写像が $\frac{1}{2}$ -Hölder連続でないことを証明

方法論の詳細

タスク定義

$\rho, \mu \in P_2(\mathbb{R}^d)$ が与えられ、 $\rho$ は絶対連続である。Brenierの定理は $\rho$ から $\mu$ への一意の最適輸送写像 $T_\mu \in L^2(\rho)$ の存在を保証する。 $\mu$ の摂動に対する $T_\mu$ の安定性を研究する。

定理1.1の構成方法

ソース測度の構成

$A = (1, 0, \ldots, 0)$ 、 $A' = (-1, 0, \ldots, 0)$ 、 $E = \{A, A'\}$ とする。関数を定義する： $f(r) = r^{-d} \min(1, (\log r)^{-2})$

密度を構成する： $\rho(x) = c_0 f(\text{dist}(x, E))$ ここで $c_0 > 0$ は $\rho$ が単位球 $B_{\mathbb{R}^d}(0,1)$ 上の確率密度となるようにする。

ターゲット測度の列

$B_\theta = (R\sin\theta, R\cos\theta, 0, \ldots, 0)$ と $B'_\theta = (-R\sin\theta, -R\cos\theta, 0, \ldots, 0)$ を考え、以下を定義する： $\mu_\theta = \frac{1}{2}(\delta_{B_\theta} + \delta_{B'_\theta})$

主要な推定

幾何学的分析を通じて以下を証明する： $\|\|T_{\mu_\theta} - T_{\mu_0}\|\|^2_{L^2(\rho)} \geq \frac{c_0 c_1 \sigma_{d-1} R^2}{|\log(\theta/4)|}$

$W_p(\mu_0, \mu_\theta) \sim R\theta/2$ であるため、安定性不等式はすべての $\alpha > 0$ に対して失効する。

定理1.2の構成方法

サポート集合の構成

無限個の矩形平行六面体対の和集合を構成する： $X = \bigcup_{i=1}^{+\infty} S_i$ ここで $S_i = T^+(A_i^+, \ell_i, r_i) \cup T^-(A_i^-, \ell_i, r_i)$ である。

パラメータ選択

列は以下を満たすように選択される：

$r_i = c_0 2^{-i}$ （超多項式的減衰）
$\ell_i = w_i = c_0 c_1 i^{-2}$
セル間距離はセル内距離よりはるかに大きい

不安定性の証明

特定の摂動下での最適輸送写像の変化を分析することで、以下を証明する： $\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}$

$\alpha > \frac{p}{2(p+1)}$ のとき、右辺は無限大に趨く。

実験設定

理論的検証

本論文は純粋な理論的研究であり、厳密な数学的証明を通じて結果を検証する：

幾何学的分析：最適輸送写像の幾何学的構造を詳細に分析
測度推定：各領域の測度を正確に計算
距離計算：Wasserstein距離と $L^2$ 距離を正確に推定

構成の検証

以下のステップを通じて構成の正確性を検証する：

構成された密度が実際に確率密度であることを証明
最適輸送写像の明示的表現を検証
幾何学的配置が必要な分離条件を満たすことを確認

密度爆発の臨界性：多項式爆発 $r^{-d+\delta}$ は依然として安定性を保つが、超多項式的爆発は完全な不安定性をもたらす
Hölder正則性の限界：写像 $\mu \mapsto T_\mu$ は特定の場合に $\frac{1}{2}$ -Hölder連続ではなく、さらに $\frac{1}{3}$ -Hölder連続でさえない
非一意性の前兆：最適計画の一意性が失われる前に安定性の悪化が生じる