[論文レビュー] Matrix models for circular ensembles

Rowan Killip, Irina Nenciu|ArXiv.org|Oct 2, 2004

Data Management and Algorithms被引用数 29

ひとこと要約

本稿では、βに依存する確率変数から構成される新しい三重対角構造を用いて、任意の逆温度βにおける円形エンsembles Gibbs分布に従う固有値を持つスパースなランダム行列モデルを構築する。主な貢献は、直交多項式を用いたユニタリ変換によって、ジャコビ・エンsembleの三重対角モデルの未解決問題を完全に解明したことである。

ABSTRACT

We describe an ensemble of (sparse) random matrices whose eigenvalues follow the Gibbs distribution for n particles of the Coulomb gas on the unit circle at inverse temperature beta. Our approach combines elements from the theory of orthogonal polynomials on the unit circle with ideas from recent work of Dumitriu and Edelman. In particular, we resolve a question left open by them: find a tri-diagonal model for the Jacobi ensemble.

研究の動機と目的

• 任意の逆温度β > 0における単位円上のクーロンガスの行列モデルを構築すること。
• ジャコビ・エンsembleの三重対角行列モデルを求める未解決問題を解くこと。
• 直交多項式理論を用いて、ランダム行列の固有値分布と円形エンsemblesのGibbs測度を統一すること。
• 円形エンsemblesを計算的に効率的かつスパースな表現（約4n個の非ゼロ要素）で表現すること。

提案手法

• 単位円板内に分布するβに依存する確率変数αk ∼ Θβ(n−k−1)+1を定義し、ほとんど確実に|αk|² < 1を満たす。
• ρk = √(1 - |αk|²)として、2×2ブロック行列Ξk = [[conj(αk), ρk], [ρk, -αk]]を構築し、スパースなブロック対角行列LおよびMを形成する。
• 行列積LM + MLが、Mを対角化するユニタリ変換Sによって、2つのジャコビ行列の直和にユニタリ同値であることを示す。
• スペクトル定理を用いて、得られたジャコビ行列のスペクトル測度が重み|Δ|βの円形エンsembles Gibbs測度と一致することを示す。
• 単位円上での直交多項式理論を活用し、スペクトル測度の再帰係数が、パラメータαkと関係することを示す。
• LM + MLの固有値分布がGibbs測度(1.5)と一致することを、ユニタリ同値性と既知のスペクトル測度を用いて確認する。

実験結果

リサーチクエスチョン

• RQ1 任意のβ > 0に対して、固有値が円形エンsembles Gibbs分布に従うスパースかつ三重対角なランダム行列モデルを構築可能か？
• RQ2 このようなモデルが、ジャコビ・エンsembleの三重対角表現を求める未解決問題を解消できるか？
• RQ3 スペクトル測度の再帰係数は、行列構築における基本的確率パラメータαkとどのように関係するか？
• RQ4 単位円上での直交多項式は、行列モデルのスペクトル測度の構築および解析において果たす役割は何か？
• RQ5 ユニタリ変換およびスペクトル理論を用いて、行列積LM + MLの固有値分布が|Δ|β重みのGibbs測度と一致することを示せるか？

主な発見

• 行列モデルLM + MLは、任意のβ > 0に対して正確にn個の固有値を持ち、重み|Δ|βのGibbs測度(1.5)に従う。
• 本研究は、ジャコビ・エンsembleの最初の明示的三重対角行列モデルを提供し、DumitriuとEdelmanが提起した未解決問題を解決した。
• 得られたジャコビ行列のスペクトル測度は、一方の枝では(1/2)(1 + x)dν(x)、他方の枝では(1/2)(1 - x)dν(x)に等しく、逆温度βにおける円形エンsemblesと一致する。
• ジャコビ行列の再帰係数は、確率的パラメータαkを用いて明示的に与えられ、b_{k+1}およびa_{k+1}はα2k、α2k+1、およびα2k-1の関数として表現される。
• 行列モデルはスパースであり、非ゼロ要素が約4n個であるため、固有値統計の効率的な数値計算が可能である。
• ユニタリ変換SはMを対角化し、LM + MLを2つのジャコビ行列の直和に変換する。これにより、円形エンsemblesとのスペクトル同値性が証明される。