[論文レビュー] The eigenvectors of Gaussian matrices with an external source

Romain Allez, Joël Bun+1 more|arXiv (Cornell University)|Dec 22, 2014

Random Matrices and Applications参考文献 15被引用数 14

ひとこと要約

本稿は、Dyson Brownian motionを介して確定的外部源によって摂動される大規模なガウス行列の固有ベクトルの動的挙動を研究する。固有ベクトルの重なりに関する自己同調的時間発展方程式を導出し、孤立した固有ベクトルに対しては、適切に中心化・スケーリングした結果、$ \sqrt{N} $ で収束するガウス分布に近づくことを証明した。その分散は時間 $ t $ に依存する明示的な式で与えられる。

ABSTRACT

We consider a diffusive matrix process $(X_t)_{t\ge 0}$ defined as $X_t:=A+H_t$ where $A$ is a given deterministic Hermitian matrix and $(H_t)_{t\ge 0}$ is a Hermitian Brownian motion. The matrix $A$ is the "external source" that one would like to estimate from the noisy observation $X_t$ at some time $t>0$. We investigate the relationship between the non-perturbed eigenvectors of the matrix $A$ and the perturbed eigenstates at some time $t$ for the three relevant scaling relations between the time $t$ and the dimension $N$ of the matrix $X_t$. We determine the asymptotic (mean-squared) projections of any given non-perturbed eigenvector $|ψ_j^0 angle$, associated to an eigenvalue $a_j$ of $A$ which may lie inside the bulk of the spectrum or be isolated (spike) from the other eigenvalues, on the orthonormal basis of the perturbed eigenvectors $|ψ_i^t angle,i eq j$. We derive a Burgers type evolution equation for the local resolvent $(z-X_t)_{ii}^{-1}$, describing the evolution of the local density of a given initial state $|ψ_j ^0 angle$. We are able to solve this equation explicitly in the large $N$ limit, for any initial matrix $A$. In the case of one isolated eigenvector $|ψ_j^0 angle$, we prove a central limit Theorem for the overlap $\langle ψ_j^0|ψ_j^t angle$. When properly centered and rescaled by a factor $\sqrt{N}$, this overlap converges in law towards a centered Gaussian distribution with an explicit variance depending on $t$. Our method is based on analyzing the eigenvector flow under the Dyson Brownian motion.

研究の動機と目的

• 確定的行列 $ A $ の固有ベクトルが、ヘルミート・ブラウン運動による摂動を受けるときの時間発展を理解すること。
• 大 $ N $ における非摂動状態と摂動状態の固有ベクトルの重なりの漸近的挙動を分析すること。
• 局所的リゾルベントおよび固有ベクトル射影に関する閉形式の時間発展方程式を導出すること。
• 孤立した固有ベクトルとその摂動された対応物との重なりに関する中心極限定理を確立すること。

提案手法

• 行列過程を $ X_t = A + H_t $ とモデル化し、$ H_t $ は i.i.d. ガウス成分を $ 1/\sqrt{N} $ でスケーリングしたヘルミート・ブラウン運動である。
• 固有値および固有ベクトルの確率微分方程式を用いて、Dyson Brownian motion下での固有ベクトルの流れを分析する。
• 固定された非摂動固有ベクトル $ |\psi_j^0\rangle $ が摂動固有ベクトル基底 $ |\psi_i^t\rangle $ に射影される量に関する自己同調的時間発展方程式 (2.5) を導出する。
• 局所的リゾルベント $ (z - X_t)^{-1}_{ii} $ を用いて、初期状態の局所的密度の時間発展を記述するBurgers型方程式を導出する。
• Stieltjes変換の技法を適用し、一般の初期行列 $ A $ に対して大 $ N $ における方程式を解く。
• 再帰的積分表現およびChaslesの関係を用いて、重なり分布のモーメントが $ g^{(2n)}(t) = (2n-1)!! \cdot g^{(2)}(t)^n $ を満たすことを証明し、孤立固有ベクトルにおけるCLTに至る。

実験結果

リサーチクエスチョン

• RQ1 確定的行列 $ A $ の固有ベクトルは、時間に依存するヘルミート・ブラウン運動の加法によってどのように変化するか？
• RQ2 大 $ N $ における非摂動固有ベクトルが摂動固有ベクトル基底に射影されるときの漸近的挙動は何か？
• RQ3 孤立固有ベクトルとその摂動された対応物との重なりは時間経過とともにどのように振る舞い、既知の分布に収束するか？
• RQ4 局所的リゾルベントの時間発展を記述する閉形式の時間発展方程式を導出できるか？
• RQ5 孤立固有ベクトルの重なり $ \langle \psi_j^0 | \psi_j^t \rangle $ の分布は何か？また、大 $ N $ におけるその分散は何か？

主な発見

• 孤立固有ベクトルの重なり $ \langle \psi_j^0 | \psi_j^t \rangle $ は、適切に中心化・スケーリングした後、$ \sqrt{N} $ で確率的にガウス分布に収束する。
• 重なりの極限分散は時間 $ t $ に依存し、初期固有値がスペクトルのバルク領域内にある位置に依存する。
• 局所的リゾルベント $ (z - X_t)^{-1}_{ii} $ に対してBurgers型の時間発展方程式を導出し、初期状態の局所的密度の時間発展を記述する。
• 一般の初期行列 $ A $ に対して、非摂動固有ベクトルの（平均二乗）射影は、固有ベクトルの流れのダイナミクスから導かれた自己同調的方程式 (2.5) を用いて計算可能である。
• 重なり分布のモーメントは $ g^{(2n)}(t) = (2n-1)!! \cdot g^{(2)}(t)^n $ を満たし、これはガウス極限を示唆し、完全なモーメント構造を提供する。
• 本手法はバルク固有値およびスパイク固有値の両方を適切に扱い、初期固有値がバルク内にあるか孤立しているかに応じて異なる挙動を示す。