Local micromechanics in a mean-field model of glasses reveal key properties of its non-equilibrium RSB phase

技術要約：平均場モデルにおけるガラスの局所的な微視力学

問題と動機
ガラスは、無秩序さと長距離秩序の欠如を特徴とする非平衡固体である。連続体スケールの特性（例：弾性）はよく理解されているが、非平衡現象（例：非デバイ型振動スペクトル、イールド転移、機械的メモリ）の微視的な起源は依然として困難な課題である。有限次元のコンピュータ・グラスにおける最近の進展により、「準局在化」した低周波モードと、振動状態密度（VDoS）$D_G(\omega)$における普遍的な$\sim\omega^4$のテイルが特定されている。

最近定式化された無秩序なハミルトニアン平均場モデル[47, 48]は、これらのガラス特性を統一する有望な手法を示している。このモデルは、特定のパラメータ範囲において、平衡状態のゼロ温度におけるレプリカ対称性の破れ（RSB）転移を特徴とする。レプリカ対称（RS）相における解は、結合していない自由度に対する有効な自己整合的ランダムポテンシャルによって記述される。しかし、このモデルが非平衡RSB相の深部（瞬時クエンチ後）にある場合の挙動、および必要とされる有効ポテンシャル記述との関連については、さらなる調査が必要であった。具体的には、空間構造を欠く完全結合モデルにおいて、「微視力学的（micromechanical）」な応答関数を定義し、局所的な観測量とグローバルな感受性の間の厳密な関係を確立する必要があった。

手法
著者らは、以下のハミルトニアンで定義される平均場ガラスモデルを解析している：
$$H = \frac{1}{2}\sum_{i=1}^N k_i x_i^2 + \frac{1}{4!}\sum_{i=1}^N x_i^4 + \sum_{i<j} J_{ij}x_i x_j - f(t)\sum_{i=1}^N x_i$$
ここで、$k_i$は初期剛性であり、分布（検討されているケースでは$[0,1]$上の$p(k) \propto k^{-1}$）から抽出される。$J_{ij}$はガウス型ランダム結合である。システムは、熱平衡化を行わずに$H$を最小化すること（瞬時クエンチ）によって、非平衡RSB相において研究される。

核心となる手法上の革新は、この平均場コンテキストにおける力モノポール応答（force monopole response）の定義である。有限次元のガラスにおける力双極子（長波長フォノンを遮蔽するもの）に類似して、著者らは局所的な力モノポール $f_i$（$i$番目の振動子への単位力）を印加し、線形応答 $z_i = M^{-1} \cdot f_i$ を測定する。ここで $M$ はヘッセ行列である。彼らは、局所的なモノポール剛性 $\kappa_i$ を以下のように定義する：
$$\kappa_i = \frac{f_i \cdot M^{-1} \cdot f_i}{f_i \cdot M^{-2} \cdot f_i} = \frac{(M^{-1})_{ii}}{(M^{-2})_{ii}}$$
ランダム行列理論（RMT）の手法を用いて、著者らはこれらの局所的な剛性と、グローバルな線形（$\chi$）およびスピングラス（$\chi_{SG}$）感受性を結びつける厳密な平均場恒等式を導出する。その後、数値シミュレーションを用いてこれらの関係を検証し、モノポール統計から完全なVDoSを再構成し、限界安定性基準に基づいて振動子の「安定化セクター（stabilized sector）」を特定する。

主な貢献と結果

1. 厳密な平均場恒等式： 著者らは、局所的なモノポール剛性 $\kappa_i$ をグローバルな感受性に結びつける厳密な関係を導出した：
   $$\kappa_i = \frac{(k_i - J^2\chi) + \frac{1}{2}x_i^2}{1 + J^2\chi_{SG}}$$
   この恒等式は、完全に結合されたモデルの性質を反映して、完全に局所的な特性がグローバルな量によって決定されることを示している。決定的なことに、$\kappa_i$ が有限であるためには、スピングラス感受性 $\chi_{SG} = \int \omega^{-4} D_G(\omega) d\omega$ が有限でなければならない。これは、モデルで観察される $\sim\omega^4$ のテイルと一致する、擬似ギャップを持つVDoSを必要とする。

2. 微視力学的局在化と擬似ギャップ： 数値解析により、モノポール応答 $z_i$ は強く局在しており（参加数 $\langle N_e \rangle \approx 4$、システムサイズ $N$ に依存しない）、これは有限次元のガラスにおける準局在モードを反映している。さらに、モノポール剛性の分布 $p(\kappa)$ は擬似ギャップ（$\kappa \to 0$ で消失）を有しており、限界安定性を示している。

3. 振動スペクトルの再構成： 著者らは、モノポール剛性 $\{\kappa_i\}$ の統計のみから、完全な非フォノン的VDoS $D_G(\omega)$ を再構成できることを示している。レゾルベントに関する自己整合方程式を、導出された剛性統計を用いて解くことにより、直接的なヘッセ行列の対角化に対して優れた定量的一致をもって $D_G(\omega)$ を復元している。

4. 特性スケールとボゾンピーク： 本研究は、非フォノン的VDoSを制御する単一の特性周波数スケールを特定している。

   • 著者らは、初期剛性が $k_i < J^2\chi$ である振動子の「安定化セクター」を定義している。これらの振動子は、正の剛性を維持するために非調和変位項（$\frac{1}{2}x_i^2$）に依存している。
   • この安定化セクターの平均剛性 $\langle \kappa \rangle_{st}$ は、固有のボゾンピーク周波数の二乗に比例する：$\langle \kappa \rangle_{st} \sim \omega_p^2$。
   • さらに、力変動 $\Delta f$ から導かれる有効ランダムポテンシャルの特性周波数スケールは、$\omega_p \sim (\Delta f)^{1/3}$ とスケールする。

5. メゾスコピック vs グローバルな弾性： メゾスコピックな剛性（$\langle \kappa \rangle_{st}$ によって特徴付けられる）とグローバルな剪断弾性率 $G$ の比は、無秩序度 $J$ によって大きく変化する。具体的には、$\langle \kappa \rangle_{st}$ は $G$ よりも広い変動範囲を示しており、これは有限次元のコンピュータ・グラスでの観察結果と一致する挙動である。

意義と主張
本論文は、平均場モデルの統計的な微視力学、特に力モノポール応答が、平衡RS相における有効ランダムポテンシャルによって記述される集団的自由度の具体的な実現を提供することを主張している。知見は以下を示唆している：

• 有効ポテンシャルの数学的構造（平衡RS相のために導出されたもの）は、非平衡RSB相の深部においても有効である。
• 「安定化セクター」の振動子（$k_i < J^2\chi$ であるもの）は、$\omega^4$ のテイルとボゾンピークの原因となるソフトなガラスモードの物理的な実体である。
• 本モデルは、平均場理論と有限次元のガラス物理学を橋渡しする統一的な枠組みを提供する。そこでは、局所的な微視力学的観測量（モノポール剛性）が、グローバルなスペクトル特性と感受性の限界を直接的に符号化している。

著者らは、これらの結果が、このモデルのエネルギーランドスケープが本質的にガラス的であるという理論的示唆を強化するものであると結論付けている。また、有効ポテンシャルの描像がその非平衡特性に対する統一的な記述を提供すると述べている。彼らは、今後の課題として、疎な結合性、非調和特性、および強い駆動ダイナミクスの探索を挙げている。