スピングラス動力学の統一的理解：実験・理論・シミュレーションの定量的収束

---

• 執筆日: 2026-03-24
• トピック: スピングラスのオフ平衡動力学 ── コヒーレンス長を軸にした実験・理論・シミュレーションの収束
• 中心論文 (main-1): J.A. Mydosh, "Spin glasses: redux: an updated experimental/materials survey", Rep. Prog. Phys. 78, 052501 (2015). DOI: 10.1088/0034-4885/78/5/052501
• 中心論文 (main-2): E.D. Dahlberg, I. González-Adalid Pemartín, E. Marinari, G. Parisi, F. Ricci-Tersenghi, V. Martin-Mayor, J. Moreno-Gordo, R.L. Orbach, I. Paga, J.J. Ruiz-Lorenzo, D. Yllanes, "Spin-glass dynamics: Experiment, theory, and simulation", Rev. Mod. Phys. 97, 045005 (2025). arXiv:2412.08381
• 参照した references 内論文数: 0（references フォルダに sub 論文なし）
• 追加検索した arXiv 論文数: 約 15 本
• Primary area: Condensed Matter Physics — Disordered Systems and Spin Glasses
• Secondary area: Statistical Mechanics, Computational Physics

---

図についての注記: 本記事で引用した論文はすべて arXiv の "nonexclusive-distrib/1.0" ライセンス、あるいは IOP Publishing / APS の通常著作権下にあり、図の再利用許可が明確でない。したがって原論文の図は使用せず、文章による詳細な説明、および概念図の案（挿入推奨箇所）で代替する。

---

1. 導入：なぜスピングラス動力学の問いは重要か

「ランダムに配置された磁気モーメントが互いに競合する相互作用を通じて凍りつく」という、一見シンプルな設定が、20 世紀物理学が生み出した最も豊かで難解な問題の一つを生む。それがスピングラス（spin glass; SG）である。

スピングラスが注目される理由は二重の意味で深い。第一に、スピングラスは「複雑系の模範例」である。多数の自由度がランダムな相互作用で結びついた系の統計力学は、磁性体を超えて、ニューラルネットワーク（Hopfield モデル）、タンパク質折り畳み、組合せ最適化問題、進化生態学、金融市場のモデルとして直接応用されてきた。2024 年のノーベル物理学賞が J.J. Hopfield と G. Hinton に与えられた文脈で、Hopfield が強調したのはまさにスピングラスの統計力学からの着想であった。第二に、スピングラスは「理論の実験場」である。レプリカ法、くりこみ群、数値シミュレーション、そして専用目的計算機の開発まで、スピングラス研究は物性物理・数理物理・計算科学の最先端を牽引してきた。

しかし約 50 年のあいだ、スピングラス研究には一つの根本的ギャップが存在した。実験と理論・シミュレーションの定量的対応が取れていなかったのである。理論は互いに競合する二つの立場（レプリカ対称性の破れ理論 vs. ドロップレット理論）が拮抗し、実験は多結晶試料の複雑さに悩まされ、シミュレーションは実験タイムスケールに遠く届かなかった。

この状況を劇的に変えた二つのブレークスルーが本記事の核心をなす。一つは CuMn 単結晶試料の育成（Ames Laboratory、D.L. Schlagel のグループ）によって、150 nm にも達するスピングラスコヒーレンス長の実験的計測が可能になったこと。もう一つは Janus II 専用スーパーコンピュータ（Janus Collaboration）の稼働により、実験と同等のタイムスケール・長さスケールでの数値シミュレーションが初めて実現したことである。

本記事は二本の主論文を軸に据える。Mydosh (2015) による実験・材料サーベイ（Rep. Prog. Phys.）は、スピングラスの正典的判別基準、理論モデルの現状、および新材料の広がりを俯瞰する。Dahlberg ら (2025) による包括的レビュー（Rev. Mod. Phys.）は、このギャップが埋まった瞬間の記録であり、コヒーレンス長というただ一つの長さスケールを通じて、エイジング・温度カオス・リジュビネーション・記憶効果がすべて統一的に理解できることを示す。

---

2. スピングラスは何を問うているのか

2.1 問いの本質：秩序とは何か

スピングラスとは何か。もっとも簡潔な定義は、Mydosh (2015) の言葉を借りれば「ガラス転移温度 $T_g$ において凍りつく、スピンの相転移」である。しかしこの「凍りつき」は、強磁性（全スピンが同方向に揃う）でも反強磁性（二副格子に交互に揃う）でもなく、各格子点ごとに異なるランダムな方向への凍結である。

秩序変数としての長距離磁化はゼロでありながら、局所的な自発磁化（自己相関 $q_i=⟨{\sigma }_i{⟩}^2$）は有限の値を持つ。この問いを定式化したのが Edwards-Anderson (EA) 秩序変数

$$q_{\mathrm{E}\mathrm{A}}=\frac{1}{N}\sum _i⟨{\sigma }_i{⟩}^2$$

であり、$T<T_g$ では $q_{\mathrm{E}\mathrm{A}}>0$ となる。${\sigma }_i$ は格子点 $i$ のスピン（イジングモデルでは $\pm 1$）、$⟨\cdots ⟩$ は熱平均を表す。

2.2 スピングラスを作る三つの要素

Mydosh (2015) は、スピングラスが生まれるための三つの本質的要素を明快に整理している。

（1）ランダム性（disorder）: スピンの位置、あるいは相互作用の符号と強さがランダムに分布していること。例えば CuMn では、Cu 母格子の中に Mn 原子がランダムに溶け込む（希薄磁性合金）。

（2）フラストレーション（frustration）: すべての相互作用を同時に満足する配置が存在しないこと。三角格子上の反強磁性体がその典型例だが、スピングラスではランダムな F（強磁性）と AF（反強磁性）結合が生み出すフラストレーションが本質的である。

（3）競合する相互作用: 金属スピングラス（CuMn, AuFe など）では、RKKY（Ruderman-Kittel-Kasuya-Yosida）相互作用が主役である。伝導電子を媒介としたこの相互作用は距離 $r_{ij}$ に対して

$$J_{ij}\propto \frac{\cos (2k_{\mathrm{F}}{r}_{ij})}{r_{ij}^3}$$

のように振動するため（$k_{\mathrm{F}}$: フェルミ波数）、磁気モーメント間の距離に応じて F にも AF にもなりうる。これがランダム配置と組み合わさることで、フラストレーションが必然的に生じる。

2.3 正典的スピングラスの四つの判別基準

Mydosh (1993, 2015) は、「本物の」スピングラス相転移と単なる超常磁性体の凍結を区別するための四つの実験的基準を提示している。

1. 交流帯磁率のカスプと周波数依存性: ${\chi }_{\mathrm{A}\mathrm{C}}(T)$ は $T_g$ にカスプを示し、その温度が測定周波数の対数に弱く依存する（$\mathrm{\Delta }{T}_g/\mathrm{\Delta }\ln \nu \sim 0.2$–$0.3$ K）。
2. ZFC/FC 磁化の既約性: ゼロ磁場冷却（ZFC）磁化と磁場中冷却（FC）磁化の乖離が $T_g$ 以下で生じる。
3. 比熱の線形性からの逸脱: $T_g$ 以下で比熱は線形温度依存性から外れ、$T_g$ 以上で幅広いピークを示すが、明確なλ型異常はない。
4. 待機時間効果（エイジング）: $T<T_g$ での磁化緩和が、試料が凍りついてから経過した時間（待機時間 $t_w$）に依存する。

概念図の挿入推奨（1）: ZFC/FC 磁化の温度依存性と熱残留磁化（TRM）の模式図。横軸: 温度 $T$、縦軸: 磁化 $M$。$T>T_g$ で3本が合流し、$T<T_g$ で FC がほぼ一定、ZFC が急落する様子。

---

3. main論文は何を新しく示したのか

3.1 Mydosh (2015) の貢献：材料世界の広がり

Mydosh (2015) のレビューは、2015 年時点でのスピングラス材料世界の全景を描く。正典的スピングラス（AuFe, CuMn, AgMn, CuFe など）から始まり、以下の多様な新材料群がスピングラス的ふるまいを示すことを整理している。

• 超伝導体: La${}_{1-x}$Sr${}_x$CuO${}_4$ などの酸化物高温超伝導体でのストライプ・ガラス状態
• 重フェルミオン系: URh${}_2$Ge${}_2$（ランダムボンド型）、CeNi${}_{1-x}$Cu${}_x$（ランダムサイト型）
• 金属間化合物・ホイスラー合金: Fe${}_2$MnAl、Pd${}_2$MnIn など
• パイロクロア・スピネル: Y${}_2$Mo${}_2$O${}_7$（幾何学的フラストレーション主体）
• 酸化物・カルコゲナイド: CdCr${}_{1.7}$In${}_{0.3}$S${}_4$（絶縁体スピングラスの模範例）
• 準結晶: Al-Mn 系準結晶でのスピングラス的ふるまい

Mydosh は特に**スピン密度波（SDW）**の概念を使って、正典的スピングラスと従来の長距離磁気秩序との違いを整理する。SDW が短距離で乱れた（short-range disordered）系がスピングラスに対応し、長距離整合 SDW が通常の AF 秩序に対応するという視点は、二つを連続的に理解する一助を与える。

さらに Mydosh はカイラルガラス（chiral glass）と量子スピングラス（quantum SG）の現状を展望する。カイラルガラスでは、スピン自由度とカイラリティ（隣接スピンの渦巻き度）が独立にあるいは同時に秩序化するという Kawamura の提案が議論される。量子スピングラスでは、低温での量子ゆらぎが古典的な凍りつきを修正し、量子臨界点近傍での異常な動力学が問題となる。

3.2 Dahlberg ら (2025) の核心的貢献：定量的収束

Dahlberg ら (2025) のレビューは、スピングラス研究史における転換点を記録する。その中心的主張は次の一文に集約できる。

単結晶試料の到達と Janus II シミュレーションによって、スピングラスのオフ平衡動力学は初めて実験と理論の定量的一致に達した。

この収束を可能にした鍵がスピングラスコヒーレンス長 $\xi (t_w,T)$ である。$t_w$ は待機時間（温度クエンチ後の経過時間）、$T$ は測定温度である。コヒーレンス長はスピングラスドメイン（局所的に秩序を持つ領域）の典型的なサイズを表し、エイジングとともに緩やかに成長する。

実験側の進展: CuMn 単結晶（6 at.%、$T_g=31.5$ K）を用いた Zhai ら (2019) の実験では、${\xi }_{\mathrm{Z}\mathrm{e}\mathrm{e}\mathrm{m}\mathrm{a}\mathrm{n}}$ が 150 nm（格子定数の約 240 倍）にまで到達した。これは多結晶試料では実現不可能な値であり、グレイン境界による RKKY 相互作用の変質を排除した純粋な測定を可能にした。

シミュレーション側の進展: Janus II は 3D イジング EA モデルを対象に設計された専用 FPGA ベースのスーパーコンピュータであり、1 モンテカルロ・ステップを 1 ps と対応付けると、実験タイムスケール（数時間〜数日）に相当するシミュレーションを実現した。シミュレーションから抽出される ${\xi }_{\mathrm{m}\mathrm{i}\mathrm{c}\mathrm{r}\mathrm{o}}$（$C_4$ 相関関数から得られる）と実験の ${\xi }_{\mathrm{Z}\mathrm{e}\mathrm{e}\mathrm{m}\mathrm{a}\mathrm{n}}$ が、スケーリング則を通じて定量的に対応付けられた（Baity-Jesi ら 2017a; Zhai ら 2020; Paga ら 2021）。

統一スケーリング則: Zhai ら (2020) および Paga ら (2021) が確立した普遍スケーリング則

$$\ln \frac{t_H^{\mathrm{e}\mathrm{f}\mathrm{f}}}{t_{H\to 0}^{\mathrm{e}\mathrm{f}\mathrm{f}}}=\frac{\hat{S}}{\hat{T}}{\xi }^{D-\theta /2}{H}^2+{\xi }^{-\theta /2}\mathcal{G}({\xi }^{D-\theta /2}{H}^2;T)$$

は、実験・シミュレーション双方のデータを、磁場 $H$、コヒーレンス長 $\xi$、温度 $T$ の全域にわたって記述する。ここで $\theta$ はレプリカン指数、$D=3$ は空間次元、$\mathcal{G}$ はスケーリング関数である。

概念図の挿入推奨（2）: ${\xi }_{\mathrm{m}\mathrm{i}\mathrm{c}\mathrm{r}\mathrm{o}}$（シミュレーション）と ${\xi }_{\mathrm{Z}\mathrm{e}\mathrm{e}\mathrm{m}\mathrm{a}\mathrm{n}}$（実験）の時間依存成長の比較図。縦軸: $\xi$（格子定数単位）、横軸: $T\ln (t_w/{\tau }_0)$。両者が同一スケーリング曲線に乗ることを示す模式図。

---

4. 背景と文脈：スピングラス理論の地形

4.1 二つの理論的立場の対立

スピングラスの低温相の性質をめぐって、半世紀にわたる理論的論争が続いてきた。主要な立場は以下の三つである。

（1）レプリカ対称性の破れ（RSB）理論: Parisi (1979) が確立し、Sherrington-Kirkpatrick (SK) モデル（無限次元全結合イジングスピングラス）の厳密解を与える。低温相には無限個の「純粋状態」が存在し、これらの状態間の重複分布 $P(q)$ が非自明な形を持つ（ultrametric 構造）。$P(q)$ は

$$q=\frac{1}{N}\sum _{i=1}^N{\sigma }_i^{(1)}{\sigma }_i^{(2)}$$

と定義される二つのレプリカ間の重複であり、RSB 相では $P(q)$ が広がりを持つ。この理論では弱い外部磁場をかけても SG 相が安定に存在し（de Almeida-Thouless 線）、外部磁場中でも相転移がある。

（2）ドロップレット理論: Fisher と Huse (1986, 1988) が提案した現象論的理論。低温相には（外部磁場なしの場合）たった二つの純粋状態（時間反転対称で関係する）しか存在しない。低エネルギー励起は「コンパクトドロップレット」（スピン反転のコンパクトな島）で記述され、スケール $L$ のドロップレットのエネルギー障壁は $\mathrm{{\rm Y}}{L}^y$ ($y>0$) のように成長する。この理論では外部磁場は SG 相を破壊し、$D=3$ では SG 相転移が存在しない。

（3）trivial-nontrivial（TNT）理論: Krzakala と Martin (2000)、Palassini と Young (2000) が提案。重複分布 $P(q)$ は RSB 的に非自明だが、リンク重複（link overlap）$P(q_l)$ はドロップレット理論的に自明（デルタ関数）に振る舞う。これは RSB とドロップレットの中間的立場である。

Dahlberg ら (2025) の重要な指摘は、この論争がオフ平衡動力学の理解には大きな障壁にならないという点である。なぜなら、動力学的観測は常に有限時間・有限スケールの問題であり、二つの理論の予測が実験タイムスケールで区別できないほど近いからである。実際、レプリカン指数 $\theta$ の実験的抽出値（$\theta \approx 0.4$–$0.5$）は RSB 理論とドロップレット理論の両方に整合し得る（Baity-Jesi ら 2018）。

4.2 Edwards-Anderson モデルとそのハミルトニアン

3D 有限次元 SG の標準的なモデルは Edwards-Anderson (EA) モデルである。ハミルトニアンは

$${\mathcal{H}}_{\mathrm{E}\mathrm{A}}=-{\sum _{⟨i,j⟩}}^{\prime }{J}_{ij}{\sigma }_i{\sigma }_j+H\sum _i{\sigma }_i$$

ここで ${\sum }^{\prime }$ は最近接格子対の和、$J_{ij}$ はランダム結合（ガウス分布または $\pm J$ 分布）、$H$ は外部磁場である。クエンチされた（時間不変の）ランダム性が本質的で、$J_{ij}$ は系の時間発展とともに変化しない。

EA モデルは $D=3$ では $T_g\approx 1.1019(29)J$ ($J$ は結合定数の標準偏差、$k_{\mathrm{B}}=1$ 単位）で相転移を示し（Baity-Jesi ら 2013）、臨界指数は $\nu \approx 2.56$, $\eta \approx -0.39$ が数値的に確立されている。

4.3 Janus コラボレーション：実験と数値の橋架け

Janus I (2008) および Janus II (2014) は、3D EA イジングスピングラスのモンテカルロシミュレーション専用に設計された FPGA ベースのコンピュータクラスターである。通常の汎用スパコンと比べて数百倍の速度でスピングラスシミュレーションを実行でき、格子サイズ $L=32$（約 $3\times {10}^4$ スピン）を $T\approx 0.64T_g$ 程度まで熱平衡化できる。

動力学シミュレーションでは 1 モンテカルロステップ $\approx$ 1 ps の換算を使い、実験で数時間〜数日かかる待機時間 $t_w\sim {10}^3$–${10}^5$ s に相当するシミュレーションが可能になった。この対応があって初めて、同一スケーリング関数上に実験点とシミュレーション点を乗せることができた（Zhai ら 2019; Baity-Jesi ら 2017a, 2018）。

---

5. 比較と解釈：何がわかり、何が揺れているか

5.1 エイジング：コヒーレンス長の成長と障壁の増大

スピングラスのエイジングとは、$T<T_g$ に急冷したあと、系が真の熱平衡から外れたまま時間とともにゆっくりと変化し続ける現象である。その本質はコヒーレンス長 ${\xi }_{\mathrm{m}\mathrm{i}\mathrm{c}\mathrm{r}\mathrm{o}}(t_w)$ の時間的成長にある。

待機時間 $t_w$ 後のコヒーレンス長は、動的臨界指数 $z(T)$ を使って

$${\xi }_{\mathrm{m}\mathrm{i}\mathrm{c}\mathrm{r}\mathrm{o}}(t_w,T)\simeq A(T)t_w^{1/z(T)}$$

と成長する。ここで $z(T)=z_c\cdot T_g/T$ の形で温度依存し、$z_c\equiv z(T_g)\approx 6.86$ が数値的に確立されている（Belletti ら 2008; Lulli ら 2016）。実験的には $z_c{T}_g/T\approx 9.62$ が CuMn で測定されており（Zhai ら 2017）、数値との差が問題となっていたが、より精密な理論式（Baity-Jesi ら 2018 の式 (54)）を用いることで一致が改善された。

コヒーレンス長が増大するにつれ、スピングラスが超えなければならない自由エネルギー障壁 ${\mathrm{\Delta }}_{\mathrm{m}\mathrm{a}\mathrm{x}}$ も増大する。Joh ら (1999) の実験と理論の組み合わせでは

$${\mathrm{\Delta }}_{\mathrm{m}\mathrm{a}\mathrm{x}}={\chi }_{\mathrm{F}\mathrm{C}}{N}_s{H}^2=k_{\mathrm{B}}T\ln t_w^{\mathrm{e}\mathrm{f}\mathrm{f}}-k_{\mathrm{B}}T\ln {\tau }_0$$

が成立し、$N_s=\xi (t_w{)}^{D-\theta /2}$ は相関スピン数、${\tau }_0\sim {10}^{-12}$ s は試みの時間スケールである。これによってコヒーレンス長の磁場依存実験（$\ln t_H^{\mathrm{e}\mathrm{f}\mathrm{f}}$ の $H^2$ 依存性）からコヒーレンス長が抽出できる。

概念図の挿入推奨（3）: 待機時間に対するコヒーレンス長の成長。縦軸: $\xi$（対数スケール）、横軸: $t_w$（対数スケール）。温度パラメータを複数取り、$T$ が低いほど成長が速い傾向を示す模式図。シミュレーション点と実験点が同一直線上に乗る様子を描く。

5.2 温度カオス：自由エネルギー景観の敏感さ

温度カオスとは、スピングラスの低温相（$T<T_g$）における平衡配置が、ごくわずかな温度変化 $\mathrm{\Delta }T$ に対して急激に変化する現象である。温度 $T_1$ での平衡配置と $T_2=T_1+\mathrm{\Delta }T$ での平衡配置の類似度は、カオスパラメータ

$$X_{T_1,T_2}^J=\frac{⟨q_{T_1,T_2}^2{⟩}_J}{(⟨q_{T_1,T_1}^2{⟩}_J⟨q_{T_2,T_2}^2{⟩}_J{)}^{1/2}}$$

で測定できる。カオスがある場合 $X\to 0$、ない場合 $X\to 1$ である。

温度カオスはドロップレット理論では長スケール ${\ell }_c\sim (\mathrm{\Delta }T{)}^{-1/(D_s/2-y)}$ のみで現れ、RSB 理論では平均場極限でも存在するが弱い。実験と数値シミュレーションの両方で、温度カオスが実際に存在することが確認された（Zhai, Orbach, Schlagel 2022）が、その現れ方は均一ではなく、まれな事象（rare event）として局所的な「カオス球」が存在する（Baity-Jesi ら 2021; Fernandez ら 2016）。

Zhai ら (2022) の実験では、CuMn 単結晶（$T_g=31.6$ K）で $T_1=18$ K における FC 磁化のサイクリックプロトコルを用い、$\mathrm{\Delta }T\approx 0.45$–$0.46$ K でカオスの開始を 10 mK の精度で検出した。このことから温度カオスの指数 $1/(D_s/2-y)\approx 1.1$ が得られており、平衡状態での理論値と整合する。

概念図の挿入推奨（4）: 温度サイクルプロトコルの模式図。FC 磁化 $M_{\mathrm{F}\mathrm{C}}(t)$ の「参照曲線」と「サイクル曲線」の重ね合わせ。可逆温度変化ではシフトにより完全重合、カオス領域では部分的にしか重合しない様子を示す。

5.3 リジュビネーション（若返り）と記憶効果

リジュビネーションとは、温度 $T_1$ でエイジングしたスピングラスを低温 $T_2<T_1$ に冷却すると、あたかも「新鮮な」スピングラスのように振る舞いを再開する現象である。対照的に記憶効果（memory）は、$T_1$ に温度を戻したとき、以前のエイジング状態の記憶が保持されている現象である。この二つの現象が「同時に」起きることが、1998 年の Jonason らの実験（CdCr${}_{1.7}$In${}_{0.3}$S${}_4$ の交流帯磁率測定）で発見された。

リジュビネーションの機構として長年議論されてきたのが温度カオスである。$T_2$ での自由エネルギー景観が $T_1$ のそれとまったく異なれば、$T_1$ でのエイジング構造は $T_2$ では不安定となり（リジュビネーション）、$T_1$ に戻ったときに元の景観が復元される（記憶）。

Dahlberg ら (2025) は、Janus Collaboration による一連の数値実験（Baity-Jesi ら 2023）と Freedberg ら (2024) の CuMn 単結晶実験を結合して、記憶効果の定量的モデルを提示する。記憶の度合いは

$$\mathcal{M}=\mathrm{\Delta }{\chi }_H^{″}/\mathrm{\Delta }{\chi }_C^{″}$$

で定量化され（${\chi }_H^{″}$: 加熱時の虚部磁化率変化、${\chi }_C^{″}$: 冷却時の変化）、この量がコヒーレンス長の比 ${\xi }_1({\xi }_2/{\xi }_1{)}^3$ の関数として普遍的にスケールすることが示された（Paga ら 2024）。

特に重要な知見は三つの長さスケールが必要という点である。${\xi }_{\mathrm{m}\mathrm{i}\mathrm{c}\mathrm{r}\mathrm{o}}(t_w,T)$（成長するコヒーレンス長）、${\xi }_{\mathrm{Z}\mathrm{e}\mathrm{e}\mathrm{m}\mathrm{a}\mathrm{n}}$（外部磁場でアクセスできるコヒーレンス長）、$\zeta (t_1,t_2)$（二時刻の相関関数から得られる「再配置長さ」）の三者を組み合わせることで、リジュビネーション・記憶の全体像が記述できる（Baity-Jesi ら 2023）。

5.4 揺れている問題：未解決点

Dahlberg ら (2025) が明示する未解決問題を整理する。

（a）低温相の場の理論: EA モデルの低温相を繰り込み群（RG）で厳密に記述する場の理論はまだ存在しない。平均場 RSB 理論は $D\ge 6$ でのみ厳密だが、$D=3$ への拡張は技術的に困難である。

（b）臨界指数の矛盾: 数値シミュレーション（イジング EA モデル）で得られる臨界指数 $\gamma =6.33\pm 0.11$ は、実験値（AgMn: $\gamma =2.3\pm 0.2$、Lévy 1988）と大きく異なる。実材料の CuMn は Heisenberg スピンガラスであり、イジングモデルとの差異が問題の核心にある。

（c）外部磁場中の相転移: de Almeida-Thouless (AT) 転移線（RSB 理論が予言する外部磁場中の SG 相転移）の存在が $D=3$ で成立するかどうかは未確定である。Miranda-Filho と Jin (2026, arXiv:2601.19192) による最新の数値的証拠は AT 転移の存在を支持するが、決着には至っていない。

（d）薄膜ダイナミクスの解釈: 薄膜スピングラス（CuMn 多層膜）では、コヒーレンス長が膜厚で制限されるため、2D 的な動力学への次元クロスオーバーが議論される。シミュレーション（Fernandez ら 2019）は 4 つの成長レジームを予言するが、実験との対応はまだ完全ではない。

---

6. 材料・手法・応用への広がり

6.1 スピングラス類似現象を示す新材料群

Mydosh (2015) が詳述するように、スピングラスはもはや希薄磁性合金（CuMn, AuFe, AgMn）に限定されない。以下の材料群でもスピングラス的ふるまいが観測されており、各々が独自の物理的問いを提供する。

高温超伝導関連物質: La${}_{2-x}$Ba${}_x$CuO${}_4$（LBCO）や La${}_{2-x}$Sr${}_x$CuO${}_4$（LSCO）の特定のドーピング領域では、スピンストライプとスピングラスが共存するスピンガラス状態が現れる。磁性と超伝導の競合を理解する文脈で注目される。

重フェルミオン化合物: URh${}_2$Ge${}_2$ は Rh サイトの非磁性不純物がランダムボンド型フラストレーションをもたらすスピングラスの模範例である。強相関電子系のコヒーレンス温度（$T^{\ast }$）以下でのスピングラス形成は、重フェルミオン状態とガラス転移の相互作用として興味深い。

幾何学的フラストレーション系（パイロクロア・スピネル）: Y${}_2$Mo${}_2$O${}_7$（パイロクロア酸化物）は反強磁性的超交換相互作用にもかかわらず、幾何学的フラストレーションにより長距離秩序を持たず SG 転移を示す（$T_g\approx 22$ K）。ZnCr${}_2$O${}_4$ などのスピネルでも類似した現象が報告されている。

ホイスラー合金: Fe${}_2$MnAl、Ni-Mn-Ga 系などのホイスラー型磁性体では、秩序-不規則変態に伴う無秩序状態でスピングラス的ふるまいが現れる。これらは形状記憶合金・磁気形状記憶合金との接点でも重要である。

6.2 実験手法の進化：単結晶育成と多様な測定手法

Dahlberg ら (2025) が強調する実験的進歩の中心は CuMn 単結晶の育成である。Ames Laboratory の Schlagel らが育成した CuMn 6 at.%（$T_g=31.5$ K）の大型単結晶は、グレイン境界によるRKKY の変質を排除し、スピングラスコヒーレンス長の成長が初めてマクロスケール（150 nm）で観測できる試料を提供した。

これに加えて、コヒーレンス長を実験的に抽出するための複数の手法が比較されている。

• 磁場依存 TRM 測定（Joh ら 1999 の方法）: $S(t,t_w;H)$ のピーク時刻 $t_H^{\mathrm{e}\mathrm{f}\mathrm{f}}$ の $H$ 依存性から $N_s\propto {\xi }^{D-\theta /2}$ を経由して $\xi$ を抽出。これが主流の方法。
• 非線形帯磁率 ${\chi }_3$: ${\chi }_3(T)\propto {\xi }^{2-\eta }$ が相転移点付近で発散し、臨界指数の決定に使われる。Gunnarsson ら (1991) による CuMn での測定が相転移の熱力学的証拠を与えた。
• 非共鳴ホールバーニング（NRHB）: Chamberlin (1999) が提案した光学的手法の磁気版。コヒーレントな局所ドメインの証拠を与える。
• メゾスコピック電気抵抗ノイズ: Weissman ら (1988–1992) の実験。伝導ゆらぎ（UCF）がスピングラス動力学の鋭敏なプローブとなる。

6.3 量子情報・機械学習との接点

スピングラスの理論的遺産は、現代の量子情報・機械学習にも直接的な影響を持つ。

ニューラルネットワーク（Hopfield モデル）: 1982 年に Hopfield が提案した連想記憶モデルは、SK スピングラスのエネルギー関数と数学的に同型である。$p$ 個のパターン $\{{\xi }_i^{\mu }\}$ を記憶させた Hopfield ネットワークのハミルトニアンは $\mathcal{H}=-\sum _{i<j}{J}_{ij}{S}_i{S}_j$（$J_{ij}=N^{-1}\sum _{\mu }{\xi }_i^{\mu }{\xi }_j^{\mu }$）であり、SK モデルと同じ形を持つ。パターン数 $p$ が $p/N\to \alpha$ の極限でスピングラス転移（連想記憶の破綻）が起きることが RSB 理論によって正確に記述される（${\alpha }_c\approx 0.138$）。2024 年のノーベル物理学賞に向けた評価委員会のコメントでも、この接続が強調されている（Dahlberg ら 2025）。

量子アニーリングと組合せ最適化: D-Wave のような量子アニーリングマシンは、その動作原理においてスピングラスエネルギー景観の探索問題と密接に関連する。スピングラスの基底状態探索は NP 完全問題であり（Barahona 1982; Istrail 2000）、効率的なアルゴリズムの開発は計算科学的にも重要である。

機械学習への応用: スピングラス的複雑な自由エネルギー景観の解析手法（レプリカ法、TAP 方程式）は、ディープラーニングの損失関数のランドスケープの理解にも応用されている（arXiv:2508.07397）。

---

7. 基礎から理解する

7.1 ランダム性とフラストレーション：スピングラスの二つの魂

スピングラスを理解するためにまず必要なのは、「ランダム性」と「フラストレーション」という二つの概念を直感的に掴むことである。

ランダム性を理解するには、銅（Cu）の結晶の中にマンガン（Mn）原子がランダムに散らばっている状況を想像せよ。Mn 原子は局所的な磁気モーメントを持つが、その位置は不規則である。Mn 原子間の交換相互作用は、間にある Cu の伝導電子を媒介として間接的に伝わる（RKKY 相互作用）。この相互作用は二つの Mn 原子の間の距離に依存して符号が振動する：近ければ強磁性的（互いに同じ方向を好む）、少し離れると反強磁性的（逆向きを好む）、さらに離れるとまた強磁性的……と続く。結果として、あるペアは「揃いたい」、別のペアは「逆になりたい」という相反する要求を持つ。

フラストレーションを理解するために、三角格子上のイジングスピンを考えよう。三つのスピン $S_1,S_2,S_3$ のすべてのペアが反強磁性（逆向きを好む）とする。$S_1=\uparrow$ とおくと $S_2=\downarrow$ が望ましいが、$S_3$ は $S_1$ に対して $\downarrow$、$S_2$ に対して $\uparrow$ を同時には取れない。「すべての相互作用を同時に満足する配置が存在しない」状態がフラストレーションである。

スピングラスではこの両者が組み合わさる。ランダムな配置のため多数のフラストレーションが生まれ、系は「すべてを妥協した」多数のほぼ等エネルギーな配置（メタ安定状態）に捕まりやすくなる。これが凍りつきの本質である。

7.2 Edwards-Anderson モデル：スピングラスの最小模型

EA モデルのハミルトニアンを再掲すると

$${\mathcal{H}}_{\mathrm{E}\mathrm{A}}=-{\sum _{⟨i,j⟩}}^{\prime }{J}_{ij}{\sigma }_i{\sigma }_j$$

ここで $J_{ij}$ はガウス分布（平均 0、分散 $J^2$）または $\pm J$ 分布に従うランダム変数であり、シミュレーション開始前に固定される（これが「クエンチされた disorder」）。この「固定されたランダム性」が、アニール（熱的に平均化される）disorder とは本質的に異なる。

系を高温から $T<T_g$ にクエンチすると、スピングラス秩序変数 $q_{\mathrm{E}\mathrm{A}}=N^{-1}\sum _i⟨{\sigma }_i{⟩}^2$ がゼロから有限の値に成長し始める。しかし真の平衡 $q_{\mathrm{E}\mathrm{A}}$ に到達するためには、成長し続けるコヒーレンス長（$\xi \propto t_w^{1/z}$）が系のサイズを超える必要があり、有限サイズ系では非常に長い時間がかかる。実際、$T\approx 0.64T_g$ で $L=32$ 格子点の系を熱平衡化するには、現状最高速のシミュレーションでも数ヶ月のコンピュータ時間を要する。

7.3 コヒーレンス長：スピングラスのルーラー

スピングラス動力学を理解する「物差し」がコヒーレンス長 $\xi (t_w,T)$ である。直感的には、「局所的にスピングラス秩序が形成されている領域のサイズ」である。

$t_w=0$ でのクエンチ直後、$\xi$ は格子定数程度の微小な値から始まる。時間が経つにつれて $\xi$ が成長し、その成長に伴ってエイジング効果（待機時間依存性）が現れる。$\xi$ の成長が遅い理由は、コヒーレント領域が成長するにつれ越えなければならない自由エネルギー障壁が大きくなるためである（障壁は $\sim {\xi }^y$ で増大）。

実験的に $\xi$ を「見る」ためには外部磁場の助けが必要である。弱い磁場 $H$ を印加すると、コヒーレント領域内のスピン全体が磁場方向に弱く磁化する（ゼーマン効果）。この磁化の $H$ 依存性を精密に測定することで、コヒーレント領域のサイズ（スピン数 $N_s\propto {\xi }^{D-\theta /2}$）が推定できる。これが ${\xi }_{\mathrm{Z}\mathrm{e}\mathrm{e}\mathrm{m}\mathrm{a}\mathrm{n}}$ の測定原理である。

数式補足: コヒーレンス長抽出の式

$$\ln \frac{t_H^{\mathrm{e}\mathrm{f}\mathrm{f}}}{t_{H=0}^{\mathrm{e}\mathrm{f}\mathrm{f}}}\propto N_s{H}^2={\xi }^{D-\theta /2}{H}^2$$

ここで $t_H^{\mathrm{e}\mathrm{f}\mathrm{f}}$ は $H$ 中での効果的待機時間、$D=3$、$\theta \approx 0.4$–$0.5$ (レプリカン指数) である。$\ln t_H^{\mathrm{e}\mathrm{f}\mathrm{f}}$ の $H^2$ 依存性をプロットすることで $\xi$ が決まる。

7.4 揺動散逸関係と静力学・動力学の対応

有限温度熱平衡系では、線形応答（磁化率）と自発ゆらぎ（自己相関関数）が揺動散逸定理

$$\chi (t,t_w)=\beta [C(t_w,t_w)-C(t,t_w)]$$

で結ばれる（$\beta =1/(k_{\mathrm{B}}T)$）。この関係はスピングラスの非平衡エイジング領域では破れ、

$$\mathcal{R}(\lambda )=\beta X[C(\lambda )]C^{\prime }(\lambda )$$

という一般化揺動散逸関係（GFDR）となる（$0\le X\le 1$ は「揺動散逸比」）。RSB 理論では GFDR が厳密に成り立ち、$X$ と平衡重複分布 $P(q)$ が $X[C]={\int }_0^{C}dqP(q)$ で関係付けられる。Hérisson と Ocio (2002) の実験的測定は、この GFDR が $D=3$ の実験系でも近似的に成立することを示した画期的な結果である。

Dahlberg ら (2025) の「静力学・動力学辞書（SDD）」は、この対応を有限時間・有限サイズの実験条件に対して明示的に拡張したもので、有効系サイズ $L_{\mathrm{e}\mathrm{f}\mathrm{f}}(t,t_w)=\xi (t+t_w)h[\xi (t+t_w)/\xi (t_w)]$ という二つのコヒーレンス長を用いたスケーリングが必要であることを示した（Baity-Jesi ら 2017b）。

---

8. 重要キーワード10個の解説

キーワード 1: スピングラス（Spin Glass, SG）

定義: ランダム配置された磁気モーメントが競合する交換相互作用によりフラストレーションを受け、特定の凍結温度 $T_g$（ガラス転移温度）以下でランダムな方向に凍りつく磁気相。長距離磁化はゼロだが、局所的な自発磁化が有限。 main 論文との接続: 両論文の中心テーマ。Mydosh (2015) は正典的判別基準を整理し、Dahlberg ら (2025) はその動力学を詳細に分析する。

キーワード 2: Edwards-Anderson 秩序変数（EA order parameter）

定義: $q_{\mathrm{E}\mathrm{A}}=N^{-1}\sum _i[⟨{\sigma }_i{⟩}_{\mathrm{t}\mathrm{h}}{]}^2$（$⟨\cdots {⟩}_{\mathrm{t}\mathrm{h}}$: 熱平均、disorder 平均の後）。SG 相では $q_{\mathrm{E}\mathrm{A}}>0$ となり、SG 相の秩序変数として機能する。平衡でのスピングラス秩序の強さを表す。 main 論文との接続: Dahlberg ら (2025) §II.A で詳説。$q_{\mathrm{E}\mathrm{A}}$ は SG 相の order parameter であり、クエンチ後の動力学の「ターゲット」として機能する。

キーワード 3: レプリカ対称性の破れ（Replica Symmetry Breaking, RSB）

定義: Parisi (1979) が SK モデルで発見した低温相の性質。自由エネルギーを計算するために $n$ 個のレプリカ（同じランダム結合を持つコピー）の平均を取るが、$n\to 0$ 極限でこれらのレプリカ間の重複分布 $P(q)$ が非自明（広がりを持つ）になる。すなわちスピングラスの低温相は無限個の純粋状態を持ち、それらの間に ultrametric な階層構造が存在する。 main 論文との接続: Dahlberg ら (2025) §II.B, III で RSB vs. ドロップレットの比較として議論。実験タイムスケールでの両者の差は限定的。

キーワード 4: コヒーレンス長（Coherence Length, $\xi$）

定義: スピングラスドメイン（局所的に凍ったスピン秩序の領域）の典型的サイズ。$C_4(r,t_w)=N^{-1}\sum _x⟨q_x(t_w)q_{x+r}(t_w)⟩$ の空間的減衰から ${\xi }_{\mathrm{m}\mathrm{i}\mathrm{c}\mathrm{r}\mathrm{o}}$ として定義。実験的には Zeeman エネルギーを使って ${\xi }_{\mathrm{Z}\mathrm{e}\mathrm{e}\mathrm{m}\mathrm{a}\mathrm{n}}$ として測定。 main 論文との接続: Dahlberg ら (2025) 全体を貫く中心的概念。エイジング・温度カオス・記憶・リジュビネーションはすべて $\xi$ の成長によって統一的に記述される。

キーワード 5: エイジング（Aging）

定義: $T<T_g$ に急冷した後、系の緩和特性（例：TRM 磁化の時間依存性）が待機時間 $t_w$ に依存して変化し続ける現象。真の平衡状態には到達せず、系が「経年」するにつれてダイナミクスがより遅くなる。 main 論文との接続: Dahlberg ら (2025) §IV.A に詳しい。コヒーレンス長の成長 $\xi \propto t_w^{1/z(T)}$ がエイジングの直接的な原因。CuMn 単結晶での系統的測定が理論を検証。

キーワード 6: 温度カオス（Temperature Chaos）

定義: スピングラスの低温相での平衡配置（Gibbs 状態）が、温度のわずかな変化 $\mathrm{\Delta }T$ に対して指数関数的に変化する性質。長さスケール ${\ell }_c\sim (\mathrm{\Delta }T{)}^{-1/(D_s/2-y)}$ より大きな系では配置が「カオス的」に変わる。 main 論文との接続: Dahlberg ら (2025) §II.D, VI.A で詳説。リジュビネーション現象の機構として議論。Zhai ら (2022) の実験で直接観測、$1/(D_s/2-y)\approx 1.1$ が確認された。

キーワード 7: リジュビネーション（Rejuvenation）

定義: $T_1$ でエイジングしたスピングラスを $T_2<T_1$ に冷却すると、$T_2$ での緩和曲線（$S(t,t_w;H)$）が $t_w$ 依存性を失い、あたかも「若い」系のように振る舞う現象。温度カオスにより $T_1$ でのコヒーレント構造が $T_2$ では無意味になるため生じる（という解釈）。 main 論文との接続: Dahlberg ら (2025) §VI.B。Janus II のシミュレーション（Baity-Jesi ら 2023）で $S_{\mathrm{Z}\mathrm{F}\mathrm{C}}(t,t_w;H)$ のピークが $t_w^{\mathrm{e}\mathrm{f}\mathrm{f}}\ll t_{w1}+t_{w2}$ に現れることが「リジュビネーション」の数値的証拠として示された。

キーワード 8: 記憶効果（Memory Effect）

定義: リジュビネーションが起きていても、元の温度 $T_1$ に戻ったとき、$T_1$ でのエイジング構造が部分的に回復される現象。スピングラスが $T_1$ での情報を「記憶」している。 main 論文との接続: Dahlberg ら (2025) §VI.B. Freedberg ら (2024) の定量化により、記憶の大きさが ${\xi }_1({\xi }_2/{\xi }_1{)}^3$ でスケールすることが判明。三つの長さスケールが必要。

キーワード 9: 静力学・動力学等価性（Statics-Dynamics Equivalence, SDE）

定義: スピングラスのオフ平衡動力学（待機時間 $t_w$ 後の系）が、同じ disorder 実現の下での有効サイズ $L\sim \xi (t_w)$ の平衡系に等価であるとする対応。揺動散逸定理の非平衡版（GFDR）で保証される。 main 論文との接続: Dahlberg ら (2025) §IV.C, IV.D で詳述。この等価性が実験とシミュレーションの橋渡しを可能にし、Janus II との定量比較の理論的基盤となった。

キーワード 10: Janus コラボレーションと専用計算機

定義: スピングラスシミュレーション専用に設計された FPGA ベースのスーパーコンピュータ群（Janus I, Janus II）を運用する国際研究コラボレーション。3D EA イジングスピングラスを対象として、格子サイズ $L$ 最大 32、温度 $T\approx 0.64T_g$ での大規模な平衡・非平衡シミュレーションを実行。 main 論文との接続: Dahlberg ら (2025) §II.H および全体。単結晶実験との定量的合致（Zhai ら 2019, 2020; Paga ら 2021, 2024）を実現した中核的なシミュレーション基盤。

---

9. まとめと今後の論点

スピングラス研究は 2020 年代に入り、「定性的描像の競争」から「定量的比較の時代」へと移行した。この転換を可能にしたのは二つの技術革新 ── CuMn 単結晶育成と Janus II 専用計算機 ── であり、両者が初めて同一のスケーリング関数上で結びついた。

確立された核心的知見: スピングラスコヒーレンス長 $\xi (t_w,T)$ は、エイジング・温度カオス・リジュビネーション・記憶効果という多様な現象を統一的に記述する物差しである。温度カオスは実験的に観測され、リジュビネーションとの因果関係が定量的に確かめられた。静力学・動力学等価性は $D=3$ の EA モデルで数値的に確立され、実験との橋渡し辞書を提供した。

残された主要課題:

理論側: $D=3$ における低温相の場の理論、臨界指数（特に $\gamma$）の実験との矛盾解消、外部磁場中の de Almeida-Thouless 転移の存否、カイラリティとスピンの転移温度の関係（特に等方的 Heisenberg スピングラスにおけるカイラルガラス転移）。

実験側: $\xi$ が大きくなった領域での GFDR のより精密な検証、薄膜動力学の 4 段階成長レジーム（シミュレーション予言）の実験的観測、より高磁場（〜数十テスラ）でのスピン磁化曲線を使った希薄スピン間相互作用の直接測定（Dahlberg ら 2025 §VII.B で予告）。

波及分野への展開: スピングラス動力学の普遍的側面 ── エイジング、カオス、記憶 ── は、ニューラルネットワークの訓練力学、量子アニーラーの最適化性能、強相関電子系のグラッシーダイナミクスとの深い接続を示唆する。特に Mpemba 効果（熱い系が冷たい系より速く冷える）のスピングラスにおける理論的予言（Baity-Jesi ら 2019）と実験での確認は、エントロピー論的な「記憶」の新しい形を示している。

最後に、2021 年のノーベル物理学賞（Parisi）と 2024 年の同賞（Hopfield）という連続したノーベル賞が、スピングラス研究の核心的な知的遺産 ── RSB 理論と Hopfield ネットワーク ── にそれぞれ与えられた事実は、この「役に立たない材料」の普遍的な重要性を雄弁に物語る。スピングラスはいまもなお、複雑系科学の「コンパス」であり続けている。

---

10. 参考にした論文一覧

中心論文

1. J.A. Mydosh, "Spin glasses: redux: an updated experimental/materials survey", Rep. Prog. Phys. 78, 052501 (2015). DOI: 10.1088/0034-4885/78/5/052501

2. E.D. Dahlberg, I. González-Adalid Pemartín, E. Marinari, G. Parisi, F. Ricci-Tersenghi, V. Martin-Mayor, J. Moreno-Gordo, R.L. Orbach, I. Paga, J.J. Ruiz-Lorenzo, D. Yllanes, "Spin-glass dynamics: Experiment, theory, and simulation", Rev. Mod. Phys. 97, 045005 (2025). arXiv:2412.08381

追加検索論文（arXiv）

3. J. Freedberg, D.L. Schlagel, R.L. Orbach, E.D. Dahlberg, "Interplay between aging and rejuvenation in spin glasses", Phys. Rev. B 112, 064404 (2025). arXiv:2406.07628

4. J. Freedberg, W.J. Meese, J. He, D.L. Schlagel, E.D. Dahlberg, R.L. Orbach, "Memory and Rejuvenation in Glassy Systems", Phys. Rev. B 110, L060411 (2024). arXiv:2305.17296

5. Janus Collaboration (M. Baity-Jesi et al.), "Temperature chaos is present in off-equilibrium spin-glass dynamics", Commun. Phys. 4, 74 (2021). arXiv:2011.09419

6. I. Paga et al. (Janus Collaboration), "Spin-glass dynamics in the presence of a magnetic field: exploration of microscopic properties", Phys. Rev. B 103, 174421 (2021). arXiv:2101.00821

7. Q. Zhai, V. Martin-Mayor, D.L. Schlagel, G.G. Kenning, R.L. Orbach, "Slowing down of spin glass correlation length growth: simulations meet experiments", Phys. Rev. B 100, 094202 (2019). arXiv:2002.07478

8. Q. Zhai, R.L. Orbach, D.L. Schlagel, "Scaling law describes the spin-glass response in theory, experiments and simulations", Phys. Rev. B 105, 014439 (2020). arXiv:2007.03871

9. M. Baity-Jesi et al. (Janus Collaboration), "The Mpemba effect in spin glasses is a persistent memory effect", Proc. Natl. Acad. Sci. 116, 15350 (2019). arXiv:1804.07569

10. M. Baity-Jesi et al. (Janus Collaboration), "A statics-dynamics equivalence through the fluctuation-dissipation ratio provides a window into the spin-glass phase from nonequilibrium measurements", Proc. Natl. Acad. Sci. 114, 1838 (2017). arXiv:1610.01418

11. M. Baity-Jesi et al. (Janus Collaboration), "Temperature chaos in 3D Ising spin glasses is present in equilibrium simulations", Phys. Rev. Lett. (2021). arXiv:2011.09419

12. R.A. Baños et al. (Janus Collaboration), "The three-dimensional Ising spin glass in an external magnetic field: the role of the silent majority", J. Stat. Mech. P07026 (2014). arXiv:1403.2622

13. M. Baity-Jesi et al. (Janus Collaboration), "Critical parameters of the three-dimensional Ising spin glass", Phys. Rev. B 88, 224416 (2013). arXiv:1310.2910

14. M. Baity-Jesi et al. (Janus Collaboration), "Dynamical transition in the D=3 Edwards-Anderson spin glass in an external magnetic field", Phys. Rev. E 89, 032140 (2014). arXiv:1307.4998

15. L. Miranda-Filho, Y. Jin, "Evidence of the de Almeida-Thouless transition in three-dimensional spin glasses", arXiv:2601.19192 (2026)

16. M. Aguilar-Janita et al., "Small field chaos in spin glasses: universal predictions from the ultrametric tree", arXiv:2403.08503 (2024)

17. G. Caprioti, S. Haas, H. Vasconcelos, M. Girardi-Schappo, "Learning About Learning: A Physics Path from Spin Glasses to Artificial Intelligence", arXiv:2601.07635 (2026)

18. I. Gonzalez-Adalid Pemartín, V. Martín-Mayor, G. Parisi, J.J. Ruiz-Lorenzo, "Finite size scaling analysis of the modified Villain's model and the 4D spin glass transition", arXiv:2307.02224 (2023) [Janus: "Quantifying memory in spin glasses"]

---

本記事は automated scheduled task として生成されました（2026-03-24）。記事の科学的内容は main-1.pdf（Mydosh 2015）および main-2.pdf（Dahlberg ら 2025）を中心に、arXiv 上の追加論文を参照して執筆されています。すべての図はライセンスの制約（nonexclusive-distrib/1.0 または IOP/APS の通常著作権）により転載不可のため、文章説明および概念図挿入推奨箇所として記載しています。