Level statistics in the fractal phase of generalized Rosenzweig--Porter models

巨大で混沌としたダンスフロアを想像してみてください。そこでは、何千もの目に見えないダンサー（量子粒子）が、自分の居場所を見つけようとしています。あるシナリオでは、彼らは完璧に同期した混沌の中で、あらゆる場所に瞬時に飛び回ります。また別のシナリオでは、隅の方で動けなくなり、その場に凍りついてしまいます。しかし、その中間には、完全に凍りついているわけでも、完全に自由なわけでもない、奇妙な中間領域が存在します。彼らは「フラクタル」なのです。彼らは広がってはいますが、海岸線や雪の結晶のように、ギザギザとした自己相似的なパターンを描いて広がり、ダンスフロアのごく一部しか占有しません。

この論文は、これらダンサーたちの「音楽」をマッピングすることについてのものです。具体的には、著者たちはこれらのダンサーが到達できるエネルギー準位（音符）の間隔に耳を傾けています。彼らは知りたいのです。もし特定の音域をズームアップした場合、それらはどのように集まるのでしょうか？磁石のように互いに反発し合うのでしょうか、それとも群衆の中の他人同士のように、互いに無関心なのでしょうか？

普遍的な「フラクタル」のリズム

研究者たちは、ローゼンツヴァイク・ポーター（RP）モデルと呼ばれる有名な数学的モデルからスタートしました。このモデルを、混沌としたダンスフロアのレシピだと考えてください。これは2つの材料を混ぜ合わせます：ランダムな初期位置のリスト（対角行列）と、ダンサーが新しい場所にジャンプするためのランダムなルール（非対角行列）です。

レシピの中のたった一つのつまみ（$\gamma$ と呼ばれるパラメータ）を調整することで、ダンスフロアを完全に混沌とした状態から、完全に凍りついた状態へと変化させることができます。中間ゾーン（$1 < \gamma < 2$）では、そのフラクタル相が得られます。

大きな疑問は、**「流れている音楽の種類が重要なのか？」**ということでした。標準的なレシピでは、「ジャンプのルール」は非常に特定の滑らかなガウス分布（完璧なベルカーブのようなもの）から描かれます。しかし、もし音楽を、トゲのあるものや、ヘビーテイルを持つもの、あるいは異なる統計的ルールに従うものに変えたらどうなるでしょうか？フラクタルのダンサーたちは、依然として同じリズムで踊り続けるのでしょうか？

主な知見：
著者たちは、「イエス」であると述べています。高度な数学的ツール（レプリカ法と自由確率論）と大規模なコンピュータ・シミュレーションを用いて、このフラクタル相における「レベル統計（音符のリズム）」は普遍的であることを発見しました。

ジャンプのルールが滑らかなガウス曲線から描かれようと、「ウィグナー」行列（すべてのジャンプが独立しているもの）であろうと、あるいは「直交不変」行列（ジャンプに特定の回転対称性があるもの）であろうと、ダンサーたちは皆、全く同じスケーリング則に従って踊り終えるのです。

彼らは、**レベル圧縮率（level compressibility）**の「マスターカーブ」を発見しました。これは、群衆がいかに「ぷにぷに」しているか（弾力性があるか）の尺度だと想像してください。

• 極めて小さなスケール（極小のエネルギー窓）では、群衆は混沌とした流体（ランダム行列理論）のように振る舞います。
• 極めて大きなスケールでは、群衆は凍りついた独立したガス（ポアソン統計）のように振る舞います。
• 中間では、群衆は**サウレス・エネルギー（Thouless energy, $E_T$）**と呼ばれる特別なエネルギー・スケールに対する窓のサイズのみに依存する、特定の滑らかな遷移曲線に従います。

彼らは、この遷移曲線がテストされたすべてのRPモデルのバリエーションに対して同一であることを証明しました。ノイズがガウス型であろうと、それ以外であろうと、関係ありません。フラクタルのダンサーたちは皆、同じビートに合わせて行進するのです。

彼らが否定したもの（そして否定しなかったもの）

著者たちは、この普遍性の限界をテストするために細心の注意を払いました。

1. 「ジャンプのルール」の特定の分布は重要ではないことを彼らは否定しました。 ジャンプがガウス分布であっても、ラデマッハー（単に +1 または -1）であっても、スパースであっても、あるいは三角形の分布であっても、結果は同じです。
2. 彼らは別のモデル、量子ランダムエネルギーモデル（QREM）を明示的にテストしました。 これは、無秩序な量子系（スピングラスのようなもの）のより現実的なモデルです。これは、紙の上ではRPモデルと非常によく似ています。
   • 結果： ここでは普遍性が崩壊します。
   • QREMのダンサーたちは、混沌から凍結への遷移を示しますが、その曲線は異なります。その曲線は、より広く、より緩やかです。
   • なぜか？ 著者らは、これはQREMのダンサーが単なるフラクタルではなく、マルチフラクタルであるためだと示唆しています。RPのダンサーは単一の「フラクタル次元」（広がりを測る一つの方法）を持っていますが、QREMのダンサーは、波動関数の異なる部分が異なるスケールで変化する、複雑で多層的な構造を持っています。単純なRPのレシピでは、この追加の複雑さを捉えることはできません。

彼らの確信度はどの程度か？

• 数学： 一般化されたRPモデル（ウィグナーおよび直交不変モデル）については、結果は解析的に導出されています。彼らは厳密な数学的証明を用いて、スケーリング関数が普遍的であることを示しました。
• シミュレーション： 彼らはこれを、巨大な行列（最大 $N=30,000$）の厳密な数値対角化によって裏付けました。コンピュータ・シミュレーションからのデータポイントは、彼らの数学的曲線と完璧に一致しました。
• 稀なイベント： 彼らはまた、極めて稀なイベント（確率が $10^{-40}$ まで下がるような事象）を調べるための特別な「大偏差」アルゴリズムも使用しました。ここでのデータは、コンピュータのサイズの限界により多少ノイズが混じりましたが、中心付近では依然として理論と一致していました。
• QREM： 量子ランダムエネルギーモデルについては、結果をシミュレーションしました。彼らはこれを数学的にまだ証明していません。彼らは数値的に観察したのです。論文は、この差異はマルチフラクタリティによるものであると示唆していますが、これはデータに基づいた強力な仮説として提示されており、最終的な数学的証明ではありません。

まとめ

もし、あるシステムが特定の単純な方法で「フラクタル」であるならば（一般化されたローゼンツヴァイク・ポーター・モデルのように）、そのエネルギー準位の間隔は普遍的な法則に従います。それはノイズの微視的な詳細には関心がなく、スケールのみに関心を持ちます。

しかし、もしシステムが「マルチフラクタル」であるならば（量子ランダムエネルギーモデルのように）、その普遍的な法則は崩れます。ダンスフロアはより複雑になり、RPモデルの単純でエレガントな曲線はもはや適合しなくなります。著者らは、RPモデルは混沌を理解するための素晴らしい「おもちゃ」ではあるものの、限界があり、マルチフラクタルな量子系の完全で混沌とした現実を捉えるには、新しいモデルが必要であると結論づけています。