[論文レビュー] Distribution Functions for Edge Eigenvalues in Orthogonal and Symplectic Ensembles: Painlevé Representations

Momar Dieng|ArXiv.org|Nov 18, 2004

Random Matrices and Applications被引用数 4

ひとこと要約

本稿は、エッジスケーリング極限におけるガウス正規（GOE）およびガウスクaternion（GSE）アンサンブルの $m^{\text{th}}$ 番目の最大固有値の分布について、明示的なペアノレ型表現を導出する。ペアノレ II 方程式の解と行列値カーネルを用いた新規なフレドホルム行列式表現を用い、トレーシー＝ワイドム法則を一般化する再帰関係を確立し、GOEとGSE間の固有値の入れ違いに関する長年の予想を解決する。具体的には、$D_1(s,\lambda)$ と $D_4(s,\lambda)$ 関数が単純なパrameter代入によって関連付けられていないことを示し、その予想が誤りであることを示す。

ABSTRACT

We derive Painlevé--type expressions for the distribution of the $m^{th}$ largest eigenvalue in the Gaussian Orthogonal and Symplectic Ensembles in the edge scaling limit. The work of Johnstone and Soshnikov (see [7], [10]) implies the immediate relevance of our formulas for the $m^{th}$ largest eigenvalue of the appropriate Wishart distribution.

研究の動機と目的

• エッジスケーリング下でのガウス正規およびクaternionアンサンブルにおける $m^{\text{th}}$ 異常最大固有値の分布について、正確なペアノレ型表現を導出すること。
• GSEにおける最大固有値の分布がGOEにおける第二最大固有値に対応するという長年の予想を、フレドホルム行列式関数 $D_\beta(s,\lambda)$ の構造を分析することで解明すること。
• ジョンストーンおよびソシュニコフの結果を拡張し、RMT普遍性を用いてWishart分布における $m^{\text{th}}$ 異常最大固有値の漸近的挙動を厳密な枠組みで提示すること。
• 常微分方程式系と境界値問題ソルバを用いた数値安定な手法を開発し、実験的Wishartデータとの正確な百分位数比較を可能にすること。

提案手法

• フレドホルム行列式 $D_\beta(s,\lambda)$ を用いて、GOEおよびGSEにおける $m^{\text{th}}$ 固有値分布の再帰関係を導出。式は $F_\beta(s,m+1) - F_\beta(s,m) = \frac{(-1)^m}{m!} \frac{d^m}{d\lambda^m} D_\beta^{1/2}(s,\lambda)\big|_{\lambda=1}$ で表され、$D_\beta(s,\lambda)$ はフレドホルム行列式を意味する。
• $D_1(s,\lambda)$ と $D_4(s,\lambda)$ が、トレーシー＝ワイドムの公式に $q(x,\lambda)$ を単純に代入することで得られるものではないことを確立。これにより、過去の予想が誤りであることが修正される。
• ペアノレ II 方程式の解 $q(x)$ 及びその摂動を含む積分を計算するための常微分方程式系 $q_0, q_0', I_0, I_0', J_0$ および $q_1, q_1', I_1, I_1', J_1$ を導入。
• 初期近似には MATLAB の ODE ソルバ `ode45` を、左側への正確な延長には境界値問題ソルバ `bvp4c` を用い、$x = -8$ を超える領域では漸近展開を貼り合わせる。
• 大きな $x$ における $q(x)$ の漸近展開を用いて、ペアノレ II 解の接続問題に近接する不安定な遷移領域を処理。
• 100 × 100 および 100 × 400 の行列に対する実験的Wishart分布と比較することで、計算された $F_1(s,m)$ 百分位数の妥当性を検証。標準的信頼水準において良好な一致を示す。

実験結果

リサーチクエスチョン

• RQ1 GOEおよびGSEにおける $m^{\text{th}}$ 異常最大固有値の分布は、ペアノレ II 超越関数を用いた明示的な関数形で表現可能か？
• RQ2 GSEにおける最大固有値がGOEにおける第二最大固有値に対応するという予想は、フレドホルム行列式 $D_\beta(s,\lambda)$ のレベルで成立するか？
• RQ3 $\beta = 1,4$ の $D_\beta(s,\lambda)$ 関数を正しく構築するには、単純なパrameter代入による誤った導関数関係を避けるにはどうすればよいか？
• RQ4 ペアノレ II 解およびその摂動を含む積分を安定的かつ正確に計算するための数値的手法は何か？
• RQ5 導出された $F_1(s,m)$ 分布は、有限サイズの設定における実験的Wishart固有値百分位数とどの程度一致するか？

主な発見

• 本稿は、$D_1(s,\lambda) = D_2(s,\tilde{\lambda}) \cdot \frac{1}{2} \left(1 + \frac{1}{2} \int_s^\infty q_0(x) \, dx \right)$ を証明し、$D_1(s,\lambda)$ が $D_2(s,\lambda)$ に対して単純なパrameter代入によって関連付けられていないことを示し、長年の予想が誤りであることを否定する。
• GSEの場合、$D_4(s,\lambda) = D_2(s,\lambda) \cdot \cosh^2\left(-\frac{1}{2} \int_s^\infty q_0(x) \, dx \right)$ が成り立つ。これはGSEの予想が正当であることを確認するが、GOEの場合は成り立たない。
• 再帰関係 $F_\beta(s,m+1) - F_\beta(s,m) = \frac{(-1)^m}{m!} \frac{d^m}{d\lambda^m} D_\beta^{1/2}(s,\lambda)\big|_{\lambda=1}$ は $\beta = 1,4$ に対して有効であり、トレーシー＝ワイドム枠組みが高次のエッジ固有値に拡張されることを示す。
• `bvp4c` と漸近展開による貼り合わせを用いた数値的妥当性評価により、計算された $F_1(s,m)$ 分布は、100 × 100 および 100 × 400 行列の両方において、90番目および95番目の百分位数で0.01未満のずれで実験的Wishart百分位数とよく一致している。
• 常微分方程式系 $\frac{d}{ds} \begin{pmatrix} q_0 \\ q_0' \\ I_0 \\ I_0' \\ J_0 \end{pmatrix} = \begin{pmatrix} q_0' \\ s q_0 + 2 q_0^3 \\ I_0' \\ q_0^2 \\ -q_0 \end{pmatrix}$ は閉じており、$F_1(s,m)$ を求めるために必要な $I_0$ および $J_0$ 積分の安定的計算を可能にする。
• 初期近似に `ode45` を、境界値問題の精密化に `bvp4c` を用いることで、高い精度を達成。$x = -8$ 近傍の数値的に不安定な領域を超えて解を延長する際には、漸近展開を用いる。