[論文レビュー] Large deviations and fluctuations of real eigenvalues of elliptic random matrices

Sung‐Soo Byun, Leslie Molag+1 more|arXiv (Cornell University)|May 4, 2023

Advanced Algebra and Geometry被引用数 5

ひとこと要約

本稿は、非エルミート性パラメータ $\tau \in [0,1)$ を持つ $N\times N$ 実楕円 Ginibre 行列における実固有値の極限挙動と揺らぎを分析する。実固有値の個数に関する中心極限定理を確立し、強非エルミート性および弱非エルミート性の両 regime において、正確な漸近的表現を導出。この表現はリーマンゼータ関数と半円重みに関する積分表現を含む普遍的スケーリング行動を明らかにする。

ABSTRACT

We study real eigenvalues of $N imes N$ real elliptic Ginibre matrices indexed by a non-Hermiticity parameter $0\leq τ<1$, in both the strong and weak non-Hermiticity regime. Here $N$ is assumed to be an even number. In both regimes, we prove a central limit theorem for the number of real eigenvalues. We also find the asymptotic behaviour of the probability $p_{N,k}^{(τ)}$ that exactly $k$ eigenvalues are real. In the strong non-Hermiticity regime, where $τ$ is fixed, we find \begin{align*} \lim_{N o\infty} \frac{1}{\sqrt{N}} \log p_{N,k_N}^{(τ)} = -\sqrt\frac{1+τ}{1-τ} \frac{ζ(3/2)}{\sqrt{2π}} \end{align*} for any sequence $(k_N)_N$ of even numbers such that $k_N = o(\frac{\sqrt N}{\log N})$ as $N o\infty$, where $ζ$ is the Riemann zeta function. In the weak non-Hermiticity regime, where $τ=1-\frac{α^2}{N}$, we obtain \begin{align*} \lim_{N o\infty} \frac{1}{N} \log p_{N,k_N}^{(τ)} \leq \frac{2}π \int_0^1 \log\left(1-e^{-α^2 s^2} ight) \sqrt{1-s^2} \, ds \end{align*} for any sequence $(k_N)_N$ of even numbers such that $k_N=o(\frac{N}{\log N})$ as $n o\infty$. This inequality is expected to be an equality.

研究の動機と目的

• 非エルミートな確率的行列集合、特にパrameter $\tau$ を持つ楕円 Ginibre エンsemble における実固有値の統計的挙動を理解すること。
• 強非エルミート性および弱非エルミート性の両 regime において、$N \to \infty$ の極限で、ちょうど $k$ 個の固有値が実数である確率 $p_{N,k}^{(\tau)}$ の漸近的表現を導出すること。
• 極限 $N \to \infty$ における実固有値個数の中心極限定理を確立し、平均回りのガウス的揺らぎを確認すること。

提案手法

• 著者らは、実楕円 Ginibre エンsemble から導かれる固有値の同時確率密度関数（JPDF）を用いる。この JPDF は Pfaffian 構造を有し、重み関数 $\omega^{(\tau)}(z) = e^{-\operatorname{Re}(z^2)/(1+\tau)} \operatorname{erfc}\left(\sqrt{\frac{2}{1-\tau^2}} \operatorname{Im} z\right)$ を含む。
• 強非エルミート性 regime（$\tau$ が固定）において、$p_{N,k_N}^{(\tau)}$ の指数的減衰率を大偏差技法を用いて分析し、リーマンゼータ関数 $\zeta(3/2)$ を含む極限を得る。
• 弱非エルミート性 regime（$\tau = 1 - \alpha^2/N$）においては、$N \to \infty$ のときの大偏差率 $\frac{1}{N}\log p_{N,k_N}^{(\tau)}$ の上界を導出し、$[0,1]$ 上の積分として表す。積分には重み $\sqrt{1-s^2}$ が含まれる。
• 特殊関数の漸近展開（特に補完誤差関数およびベッセル型積分）を用い、$\Gamma$-関数および誤差関数を含む和の評価に必要な不等式を導出する。
• 行列式の係数抽出として $p_{N,k}^{(\tau)}$ を表現するため、$\beta_{2j-1,2m}[1]$ および $\alpha_{2j-1,2m}[1]$ を含む行列式の確定的公式を用いる。
• 中心極限定理の証明にはモーメント母関数法と鞍点近似を用い、双曲線余弦およびガウス重み関数を含む積分の収束に関する補助的補題を導入する。

実験結果

リサーチクエスチョン

• RQ1 $N \to \infty$ のとき、$N \times N$ 実楕円 Ginibre エンsemble において、ちょうど $k$ 個の固有値が実数である確率 $p_{N,k}^{(\tau)}$ の漸近的挙動は何か？
• RQ2 強非エルミート性 regime（$\tau$ が固定）において、実固有値個数の揺らぎはどのように振る舞うか？
• RQ3 弱非エルミート性 regime（$\tau = 1 - \alpha^2/N$）において、$p_{N,k_N}^{(\tau)}$ の大偏差率は何か？
• RQ4 弱 regime における大偏差率の導出された上界は、実際にタイト（等号が成り立つ）であるか？
• RQ5 実固有値個数は両 regime において中心極限定理を満たすか？

主な発見

• 強非エルミート性 regime（$\tau$ が固定）では、大偏差率が $\displaystyle\lim_{N\to\infty}\frac{1}{\sqrt{N}}\log p_{N,k_N}^{(\tau)} = -\sqrt{\frac{1+\tau}{1-\tau}}\frac{\zeta(3/2)}{\sqrt{2\pi}}$ を満たす。ただし $k_N = o(\sqrt{N}/\log N)$ のとき。
• 弱非エルミート性 regime（$\tau = 1 - \alpha^2/N$）では、大偏差率は $\displaystyle\frac{2}{\pi}\int_0^1 \log\left(1 - e^{-\alpha^2 s^2}\right)\sqrt{1-s^2}\,ds$ で上から抑えられ、この不等式は実際に等号であると予想される。
• 実固有値個数は両 regime において中心極限定理を満たし、平均回りのガウス的揺らぎが確認される。
• 弱 regime では確率 $p_{N,k}^{(\tau)}$ は $N$ に指数的に減少し、強 regime では $\sqrt{N}$ に指数的に減少する。これは異なる普遍性クラスを示している。
• $\tau = 0$ の場合、正確な値は $p_{2,0}^{(0)} = 1 - \frac{\sqrt{2}}{2}$ および $p_{4,0}^{(0)} = \frac{18 - 11\sqrt{2}}{16}$ であり、既知の結果と整合する。
• 確率 $p_{N,k}^{(\tau)}$ の行列式公式は、係数 $\alpha_{2j-1,2m}[1]$ および $\beta_{2j-1,2m}[1]$ を通じて表現され、$\beta_{2j-1,2m}[1]$ は誤差関数および単項式を含む複素積分を含む。