A uniform rate of convergence for the entropic potentials in the quadratic Euclidean setting

この論文「A UNIFORM RATE OF CONVERGENCE FOR THE ENTROPIC POTENTIALS IN THE QUADRATIC EUCLIDEAN SETTING（二次ユークリッド空間におけるエントロピーポテンシャルの一様収束率）」は、最適輸送問題の正則化手法であるエントロピー正則化（エントロピー輸送）において、正則化パラメータ $\varepsilon$ が 0 に近づくときのポテンシャル関数及其勾配の収束速度を、コンパクト集合上の一様収束（uniform convergence）の観点から定量化したものです。

以下に、問題設定、手法、主要な貢献、結果、および意義について詳細にまとめます。

1. 問題設定 (Problem)

• 背景: 2 つの確率測度 $\mu, \nu$ 間の二次コスト（ユークリッド距離の二乗）に基づく最適輸送問題（Brenier 写像 $\nabla \phi_0$ の計算）は、非線形偏微分方程式（Monge-Ampère 方程式）を解く必要があり、数値的に困難です。
• エントロピー正則化: この問題を解くための実用的な手法として、目的関数にエントロピー項を加えた「エントロピー正則化最適輸送問題」が提案されています。この問題は Sinkhorn アルゴリズムを用いて効率的に計算可能です。
• 既存の知見: $\varepsilon \to 0$ のとき、エントロピー輸送の解は最適輸送の解に収束することは知られていますが、その収束速度に関する定量的な評価は限定的でした。
  • 値関数（コスト）の収束：$\varepsilon \log \varepsilon$ などの展開が知られています。
  • 輸送プラン（結合）の収束：$W_2$ 距離で $O(\sqrt{\varepsilon})$ など。
  • ポテンシャルの勾配の収束：$L^2(\mu)$ ノルムでの評価は $O(\varepsilon)$ や $O(\varepsilon^2)$ などが研究されていますが、コンパクト集合上の一様収束（$L^\infty$ ノルム）の速度については、連続設定（連続測度）において明確な評価がなされていませんでした。

本研究の目的: 二次コストのもとで、絶対連続な測度 $\mu, \nu$ に対して、エントロピーポテンシャル $\phi_\varepsilon$ およびその勾配 $\nabla \phi_\varepsilon$ が、Brenier ポテンシャル $\phi_0$ およびその勾配 $\nabla \phi_0$ に、コンパクト集合上でどの程度の速度で一様に収束するかを評価することです。

2. 手法と仮定 (Methodology & Assumptions)

著者は、ガウス測度の場合の明示的な計算を「アペタイザー」として提示し、一般的なケースに対して以下の仮定を置いた上で解析を行っています。

主要な仮定:

1. (A1) 測度の構造と凸性:
   • $\mu(dx) = e^{-V(x)}dx$, $\nu(dy) = e^{-W(y)}dy$ と表せる。
   • $V$ は滑らかで、そのヘッセ行列が上方有界（$\nabla^2 V \preceq \alpha I_d$）。
   • $W$ は滑らかで、そのヘッセ行列が下方有界（$\nabla^2 W \succeq \beta I_d$）。
   • この仮定により、Brenier 写像 $\nabla \phi_0$ とエントロピー写像 $\nabla \phi_\varepsilon$ がリプシッツ連続であり、そのリプシッツ定数が制御可能になります。
2. (A2) ポアンカレ不等式: 測度 $\mu$ がポアンカレ不等式を満たす（これにより、関数の $L^2$ ノルムをその勾配の $L^2$ ノルムで制御できます）。
3. (A3) 微分エントロピーの有限性: $\mu$ の微分エントロピーが有限である。

解析手法:

• ガリャルド・ニレンベルグ不等式 (Gagliardo-Nirenberg inequality) の応用:
  関数の $L^\infty$ ノルム（最大値）を、高階微分の $L^\infty$ ノルムと関数自体の $L^2$ ノルムで評価する不等式を利用します。
  $$\|\nabla u\|_\infty \le C \|\nabla^2 u\|_\infty^\theta \|u\|_2^{1-\theta}$$
  ここで、$\theta$ は次元 $d$ に依存するパラメータです。
• ステップ:
  1. コンパクト集合 $K$ 上で、差 $\phi_\varepsilon - \phi_0$ を切断関数 $w$ を用いて局所化します。
  2. 仮定 (A1) より、$\nabla^2 \phi_\varepsilon$ と $\nabla^2 \phi_0$ の $L^\infty$ ノルムが有界であることを示し、高階微分項を制御します。
  3. 既存の結果（MS25 など）およびポアンカレ不等式 (A2) を用いて、$L^2(\mu)$ ノルムにおける差 $\|\phi_\varepsilon - \phi_0\|_{L^2(\mu)}$ が $O(\sqrt{\varepsilon})$ 程度であることを利用し、低次の項を評価します。
  4. これらを組み合わせて、一様ノルム $\|\nabla \phi_\varepsilon - \nabla \phi_0\|_{K, \infty}$ の収束率を導出します。

3. 主要な結果 (Key Results)

定理 1.2 (勾配の一様収束率):
仮定 (A1)-(A3) を満たす測度 $\mu, \nu$ に対して、任意のコンパクト集合 $K \subset \mathbb{R}^d$ に対し、ある定数 $C_{grad} > 0$ が存在し、以下の不等式が成り立ちます。
$$\|\nabla \phi_\varepsilon - \nabla \phi_0\|_{K, \infty} \le C_{grad} \, \varepsilon^{\frac{1}{d+4}}$$
ここで、$\|\cdot\|_{K, \infty}$ は集合 $K$ 上での最大値ノルムです。

定理 1.3 (ポテンシャルの一様収束率):
さらに、ポテンシャルの正規化条件（$\int \phi_\varepsilon d\mu = \int \phi_0 d\mu = 0$）を満たす場合、コンパクトかつ連結な集合 $K$ に対して、
$$\|\phi_\varepsilon - \phi_0\|_{K, \infty} \le C_{pot} \left( \varepsilon^{\frac{1}{d+4}} + \varepsilon \right)$$
が成り立ちます。

補題 1.1 (ガウス測度の場合):
$\mu, \nu$ が非退化なガウス測度である場合、収束率は $O(\varepsilon)$ であり、コンパクト集合の半径 $R$ と共分散行列のノルムに依存する明示的な定数で評価可能です。これは、一般的な結果が最適ではない可能性を示唆しています。

4. 意義と貢献 (Significance)

1. 連続設定における初の一様収束率:
   従来の研究では、離散設定や半離散設定、あるいは $L^2$ ノルムでの収束評価が中心でした。本論文は、連続設定（絶対連続測度）において、コンパクト集合上での一様収束（$L^\infty$）の速度を初めて定量化した点で画期的です。
2. 数値解析への示唆:
   エントロピー正則化は数値計算（Sinkhorn アルゴリズム）で広く利用されていますが、$\varepsilon$ を小さくしたときの誤差の挙動、特に空間的な一様誤差について理論的な保証を提供しました。これは、最適輸送マップの近似精度を評価する上で重要です。
3. 凸解析と関数不等式の統合:
   凸関数の性質（凸性による勾配の収束）と、ガリャルド・ニレンベルグ不等式、ポアンカレ不等式を巧妙に組み合わせることで、微分可能性の仮定が厳密な場合でも、測度の構造（対数凹性など）から収束性を導出する枠組みを構築しました。
4. 次元依存性の明示:
   収束率が $\varepsilon^{\frac{1}{d+4}}$ となることは、次元 $d$ が高くなるほど収束が遅くなることを示しており、高次元問題におけるエントロピー正則化の限界や、適切な $\varepsilon$ の選択に関する指針を与えます。

まとめ

この論文は、エントロピー正則化最適輸送の理論的基盤を強化し、特に「勾配の一様収束」という重要な性質について、具体的な収束速度の上限を導出しました。ガウス測度の場合には $O(\varepsilon)$ の高速な収束が得られる一方、一般の測度に対しては次元に依存する $O(\varepsilon^{\frac{1}{d+4}})$ の評価を与え、連続最適輸送の近似理論における重要な進展となっています。