[論文レビュー] Convergence Rate for Spectral Distribution of Addition of Random Matrices

Zhigang Bao, László Erdős+1 more|arXiv (Cornell University)|Jun 9, 2016

Random Matrices and Applications参考文献 19被引用数 3

ひとこと要約

本稿は、$ A $ と $ B $ が確定的エルミート行列で、$ U $ がハール分布を従うユニタリ行列であるとき、和 $ H = A + UBU^* $ のスペクトル分布に対する最適収束速度 $ N^{-1+\gamma} $ を確立する。著者らは、特別に構築された補助量 $ Z_i $ を用いて、新たなフラクチュエーション平均化機構を導入し、バルクスペクトルにおいてこの収束速度を達成した。これは対数因子を除いて通常の固有値間隔と一致する。

ABSTRACT

Let $A$ and $B$ be two $N$ by $N$ deterministic Hermitian matrices and let $U$ be an $N$ by $N$ Haar distributed unitary matrix. It is well known that the spectral distribution of the sum $H=A+UBU^*$ converges weakly to the free additive convolution of the spectral distributions of $A$ and $B$, as $N$ tends to infinity. We establish the optimal convergence rate ${\frac{1}{N}}$ in the bulk of the spectrum.

研究の動機と目的

• バルクスペクトルにおいて、$ H = A + UBU^* $ のスペクトル分布が自由加法畳み込み $ \mu_A \boxplus \mu_B $ に収束する最適速度を確立すること。
• ハールユニタリ行列に要素ごとの独立性が欠如していることによる、標準的なフラクチュエーション平均化手法の複雑さを克服すること。
• 高次モーメントの計算を回避するため、部分積分に基づく再帰的モーメント推定を用いる、フラクチュエーション平均化の新手法を開発すること。
• グリーン関数の正規化トレース $ \frac{1}{N}\mathrm{Tr}\,G $ が、高々 $ N^{-1+\gamma} $ のオーダーで変動することを示し、これが最適収束を示唆すること。
• 経験的スペクトル測度 $ \mu_H $ が $ N^{-1} $ のレートで集中することを証明すること。これは固有値間隔のスケールと一致する。

提案手法

• グリーン関数の対角成分 $ G_{ii} $ と $ (UBU^*G)_{ii} $ の線形結合である補助量 $ Z_i $ を導入し、$ N^{-1/2} $ スケールのフラクチュエーションの代数的キャンセレーションを活用すること。
• 残留する $ N^{-1/2} $ スケールのフラクチュエーションを吸収する補正項 $ \Upsilon $ を定義し、それが実際に $ N^{-1/2} $ ではなく $ N^{-1+\gamma} $ のオーダーであることを示すこと。
• ガウス変数に基づく部分積分に基づく再帰的モーメント推定技術を用い、複雑な高次モーメント展開を回避すること。
• ハールユニタリ行列の部分的確率的分解を活用し、グリーン関数成分内の依存構造を制御すること。
• ヘルファ・シュトローム関数計算を適用し、グリーン関数トレース推定からスペクトル測度の収束速度を導出すること。
• スケール $ N^{-1+\gamma} $ における局所法を基盤として、バルク領域 $ z \in \mathcal{S}_\mathcal{I}(a,b) $ で作業すること。

実験結果

リサーチクエスチョン

• RQ1 バルクスペクトルにおいて、$ H = A + UBU^* $ の経験的スペクトル分布が自由加法畳み込み $ \mu_A \boxplus \mu_B $ に収束する最適速度は何か？
• RQ2 ハールユニタリ行列に要素ごとの独立性がない状況で、フラクチュエーション平均化はどのように達成できるか？
• RQ3 収束速度を $ N^{-2/3} $ を超えて向上させることは可能か？ もしそうなら、どのような新技術が必要か？
• RQ4 経験的スペクトル測度 $ \mu_H $ がその期待値の周りでどの程度集中するか？
• RQ5 グリーン関数トレース $ \frac{1}{N}\mathrm{Tr}\,G $ の振るまいは、$ \mu_H $ の収束にどのように関係するか？

主な発見

• バルク領域において、$ H = A + UBU^* $ のスペクトル分布が $ \mu_A \boxplus \mu_B $ に収束する速度は、任意の $ \gamma > 0 $ に対して $ N^{-1+\gamma} $ であり、$ N^\gamma $ 要素を除いて最適である。
• グリーン関数の正規化トレース $ \frac{1}{N}\mathrm{Tr}\,G $ は、高々 $ N^{-1+\gamma} $ のオーダーで変動し、最適スケールと一致する。
• 新規補助量 $ Z_i $ を用いたフラクチュエーション平均化機構が、代数的構造を通じて $ N^{-1/2} $ スケールのフラクチュエーションをキャンセルする形で、成功裏に実装された。
• 残留フラクチュエーションを吸収するために導入された補正項 $ \Upsilon $ は、直観的なサイズとは異なり、$ N^{-1/2} $ ではなく $ N^{-1+\gamma} $ のオーダーであることが示された。
• $ \mu_H $ がその期待値の周りに集中するレートは $ N^{-1} $ であり、これはバルク領域で最適である。
• 結果は $ \mathcal{S}_\mathcal{I}(a,b) $ 内のすべての $ z $ に対して一様に成り立ち、収束速度は典型的な固有値間隔 $ N^{-1} $ と一致するため、鋭いものである。