2025-11-19T17:46:14.285599

The smallest eigenvalue of $Î²$-Laguerre and $Î²$-Jacobi ensembles and multivariate orthogonal polynomials

Jeong

We study the smallest eigenvalue statistics of the $Î²$-Laguerre and $Î²$-Jacobi ensembles. Using Kaneko's integral formula, we show that the smallest eigenvalue marginal density and distribution functions of the two ensembles for any $Î²>0$ can be represented in terms of multivariate Laguerre and Jacobi polynomials evaluated at a multiple of the identity, provided that the exponent of $x$ in the Laguerre and Jacobi weights is an integer. These representations are readily computable in explicit form using existing symbolic algorithms for multivariate orthogonal polynomials. From these expressions, we derive new differentiation formulas for the multivariate Laguerre and Jacobi polynomials. Furthermore, we derive explicit solutions to the Painleve V and VI differential equations associated with the smallest eigenvalue of the LUE and JUE. We provide numerical experiments and examples.

academic

β-Laguerre および β-Jacobi アンサンブルの最小固有値と多変量直交多項式

基本情報

論文ID: 2510.04422
タイトル: The smallest eigenvalue of β-Laguerre and β-Jacobi ensembles and multivariate orthogonal polynomials
著者: Sungwoo Jeong (コーネル大学)
分類: math-ph math.MP
発表日: 2025年10月9日 (arXiv v2)
論文リンク: https://arxiv.org/abs/2510.04422

要旨

本論文は、β-Laguerre および β-Jacobi アンサンブルの最小固有値の統計的性質を研究している。Kaneko 積分公式を利用して、著者は任意の β > 0 に対して、Laguerre および Jacobi 重み関数における x の指数が整数である場合、両アンサンブルの最小固有値の周辺密度および分布関数が、単位行列のスカラー倍で評価された多変量 Laguerre および Jacobi 多項式で表現できることを証明した。これらの表現は、既存の多変量直交多項式の記号計算アルゴリズムを用いて明示的な形で計算可能である。これらの表現に基づいて、著者は多変量 Laguerre および Jacobi 多項式の新しい微分公式を導出し、LUE および JUE の最小固有値に関連する Painlevé V および VI 微分方程式の明示的な解を得た。

研究背景と動機

中心的な問題

本論文が解決しようとする中心的な問題は、β-Laguerre アンサンブル (LβE) および β-Jacobi アンサンブル (JβE) の有限 n における最小固有値統計の精密な表現である。これらのランダム行列アンサンブルは、統計学では Wishart 行列および MANOVA 行列としても知られている。

重要性

理論的意義: これらのアンサンブルはランダム行列理論における最も重要なモデルの一つであり、広範な理論的価値を持つ
統計的応用: 多変量統計において、最小固有値の分布は仮説検定および統計的推論に重要な意義を持つ
物理的応用: 量子カオス、核物理などの分野で重要な応用を持つ

既存手法の限界

Fredholm 行列式表現: 精密であるが、計算複雑度が高い
超幾何関数表現: 一般的なパラメータに対して計算が困難
Painlevé 方程式: 微分方程式の数値解法が必要であり、明示的な表現が不足している

研究動機

著者の動機は、Kaneko 積分公式を利用して最小固有値統計と多変量直交多項式の間に直接的な関連を確立し、より効率的で計算可能な明示的表現を得ることである。

中核的な貢献

新しい精密表現: LβE および JβE の最小固有値の周辺密度関数 (PDF) および分布関数 (CDF) の新しい精密表現を得た。これらは多変量 Laguerre および Jacobi 多項式で表現される
微分公式: 多変量 Laguerre および Jacobi 多項式の新しい微分公式を導出した
Painlevé 方程式解: LUE および JUE の最小固有値に関連する Painlevé V および VI 方程式の明示的な有理解を発見した
計算効率: 既存の記号計算アルゴリズムを用いてこれらの関数を効率的に計算する方法を提供した

方法論の詳細

タスク定義

研究タスクは、β-Laguerre および β-Jacobi アンサンブルの最小固有値 λ_min に対して、その確率密度関数 f(x) および累積分布関数 F(x) の精密な表現式を見つけることである。

中核的な方法論フレームワーク

Kaneko 積分公式

著者の方法は、Kaneko が 1993 年に提唱した一般化 Selberg 積分公式に基づいている：

Jacobi の場合： $\int_{[0,1]^n} \prod_{j=1}^{\nu} \prod_{i=1}^{n} (x_i - y_j) w_{n,\gamma_1,\gamma_2}^{J,\beta}(x_1,\ldots,x_n) dx_1 \cdots dx_n = \frac{1}{Z_{n,\gamma_1+\nu,\gamma_2}^{J,\beta}} P_{\tilde{\beta}}^{n^{(\nu)},\tilde{\gamma_1},\tilde{\gamma_2}}(y_1,\ldots,y_\nu)$

Laguerre の場合： $\int_{\mathbb{R}_+^n} \prod_{j=1}^{\nu} \prod_{i=1}^{n} (x_i - y_j) w_{n,\gamma}^{L,\beta}(x_1,\ldots,x_n) dx_1 \cdots dx_n = D_{n,\gamma}^{L,\beta} \times L_{\tilde{\beta}}^{n^{(\nu)},\tilde{\gamma}}(\frac{2y_1}{\beta},\ldots,\frac{2y_\nu}{\beta})$

双対パラメータ

双対パラメータを導入する：

$\tilde{\beta} := \frac{4}{\beta}$
$\tilde{\gamma} := \frac{2}{\beta}(\gamma + 1) - 1$

矩形分割

矩形分割記号 $n^{(m)} := [n,\ldots,n] \vdash nm$ （m 個の n の分割）を使用する

技術的な革新点

変数置換技法: 巧妙な変数置換により最小固有値問題を Kaneko 積分に変換
単位行列のスカラー倍での多変量多項式の評価: 複雑な多変量関数を単変量関数に簡約
矩形分割の特殊性: 矩形分割インデックスを持つ多変量多項式の特殊構造を利用

実験設定

数値検証方法

著者は以下の方法を用いて数値検証を行った：

ランダム行列の生成

LβE: Dumitriu-Edelman の双対角行列モデルを使用
JβE: Killip-Nenciu の三対角行列モデルを使用

多項式計算

Dumitriu-Edelman-Shuman アルゴリズムの記号計算を使用。Gram-Schmidt 直交化プロセスに基づく

サンプルサイズ

各実験で 10^6 個の最小固有値サンプルを生成

テストパラメータ

複数のパラメータ組み合わせをテスト：

有理数 β 値（例：β = 5/2, 4/3, 1/2）
無理数 β 値（例：β = e）
異なる行列次元 n
異なるパラメータ γ, γ₁, γ₂

実験結果

主要な結果

Laguerre アンサンブル

定理 3.3: β > 0 かつ γ > -1 が整数である場合、n×n LβE の最小固有値の密度および分布は：

$F_{n,\gamma}^{L,\beta}(x) = 1 - \frac{1}{L_{\tilde{\beta}}^{n^{(\gamma)}, \frac{2}{\beta}-1}(0_\gamma)} \times e^{-\frac{nx}{2}} L_{\tilde{\beta}}^{n^{(\gamma)}, \frac{2}{\beta}-1}(-\frac{2x}{\beta}I_\gamma)$

$f_{n,\gamma}^{L,\beta}(x) = \frac{nZ_{n,\gamma}^{L,\beta}}{Z_{n-1,\beta+\gamma}^{L,\beta} \cdot L_{\tilde{\beta}}^{(n-1)^{(\gamma)}, \frac{2}{\beta}+1}(0)} \times x^\gamma e^{-\frac{nx}{2}} L_{\tilde{\beta}}^{(n-1)^{(\gamma)}, \frac{2}{\beta}+1}(-\frac{2x}{\beta}I_\gamma)$

Jacobi アンサンブル

定理 4.3: β > 0、γ₁, γ₂ > -1 かつ γ₁ が整数である場合：

$F_{n,\gamma_1,\gamma_2}^{J,\beta}(x) = 1 - (1-x)^{n(1+\gamma_1+\gamma_2+\beta(n-1)/2)} P_{\tilde{\beta}}^{n^{(\gamma_1)}, \frac{2}{\beta}-1, \tilde{\gamma_2}}(-\frac{x}{1-x}I_{\gamma_1})$

数値検証結果

すべてのテストケースで理論公式とモンテカルロシミュレーション結果が完全に一致し、理論結果の正確性を検証した。

微分公式

系 3.3.1: 多変量 Laguerre 多項式は以下の微分公式を満たす： $\frac{d}{dx}L_\beta^{n^{(\gamma)}, \frac{\beta}{2}-1}(xI_\gamma) = R_n^\beta x^\gamma L_\beta^{(n-1)^{(\gamma)}, \frac{\beta}{2}+1}(xI_\gamma) - \frac{n}{\beta}L_\beta^{n^{(\gamma)}, \frac{\beta}{2}-1}(xI_\gamma)$