Netsu Sokutei 42 (4), 135-141 (2015)

解説
ガラス転移理論の最近の発展
宮崎 州正 a, 尾澤 岬 a, 池田 昌司 bab名古屋大学理学研究科
京都大学福井センター（受取日：2015 年 8 月 30 日，受理日：2015 年 9 月 24 日）

Recent Development of Theory of the Glass Transition
Kunimasa Miyazakia, Misaki Ozawaa, and Atsushi IkedababDepartment of Physics, Nagoya University
Fukui Institute for Fundamental Chemistry, Kyoto University(Received Aug. 30, 2015; Accepted Sep. 24, 2015)

The nature of the glass transition of supercooled liquids is still surrounded by controversy. The current status oftheoretical understanding of the transition is much more like that before the Currie-Weiss theory for critical phenomena
in a sense that the full-fledged mean field theory is still not at hand. This, however, does not mean that we suffer from(literally glass-like) stagnancy in the progress. Contrarily, we witness a revolutionary development of the mean-field
description of the slow dynamics near the transition escorted by the bona fide thermodynamic transition at a finite
temperature. In this review, we start with the introductory discourse of the phenomenology of the glass transition andoutline the whole picture of the glass transition which the contemporary version of the mean-field theory envisions.
Furthermore, we present the results of our recent numerical study for a randomly pinned liquid, which strongly support
the mean-field scenario. This is the first time forever that the true glass state with thermodynamic singularity has been
generated in silico.
Keywords: Glass Transition, Mean-field theory, Random First Order Transition
宮崎 州正
Kunimasa Miyazaki
E-mail: miyazaki@r.phys.nagoya-u.ac.jp

尾澤 岬
Misaki Ozawa
E-mail: misaki.ozawa@r.phys.nagoya-u.ac.jp

池田 昌司
Atsushi Ikeda
E-mail: atsushi.ikeda@fukui.kyoto-u.ac.jp© 2015 The Japan Society of Calorimetry and Thermal Analysis

Netsu Sokutei 42 (4) 2015
135

解

説

の状態のことだ。問題は, 急冷速度をゼロにした極限でどうなるかである。この液体のエントロピー曲線を外挿してガラス転移とは, 急冷した液体が結晶化をし損なって, いくと, 結晶のエントロピー曲線と交わってしまう。こん
分子の配置がアモルファス（ランダム）な状態に保たれた なことが実際に起こっては問題である。なぜなら, 分子配まま凍結してしまう現象である。
「転移」という名を冠して 置が規則正しい結晶状態よりも, ランダムな状態のほうがいるが, この現象はいわゆる相転移だろうか。実験, 理論, エントロピーが小さいということになるからだ。この矛盾
シミュレーションによる精力的な研究が行われているが, を, それを指摘した研究者の名にちなんで「カウツマンの
この素朴な疑問に対する答えを知る者はいない。だからと パラドクス」と呼んでいる。この矛盾を避けるためには, 液
いって, ガラス転移の研究が停滞していたわけではない。 体のエントロピー曲線は下がり続けることはなくどこかで,むしろ, 他の研究分野も巻き込みながら, めまぐるしく進 「何か」が起こる, つまり, 相転移が起こらなくてはならな
歩している。本解説の目的は, ガラス転移の統計力学, 特に いという結論になる。多くの研究者は, 液体と結晶のエン
その理論的な側面に焦点を当てて, 最近の理論研究の目覚 トロピー曲線が交わる温度こそ相転移が起こる温度である,と信じていて, それを「理想ガラス転移点」, あるいは「カましい発展を, 非専門家向けに説明することである。
そもそもガラスとは何だろう。ガラスの代表格といえば, ウツマン温度」TK と呼んでいる。
窓ガラスや食器に使われるシリカガラスである。しかし構
この議論をダイナミクスの側面から考えてみよう。様々な成要素が原子や分子である必要はなく, コロイドや高分子 過冷却液体の粘性係数 の温度依存性を示したのが Fig.1(b)のような比較的大きなものであっても, それがランダムで である。横軸には温度の逆数を, 縦軸には を対数で示して流れなければガラスと言ってよい。このように定義すると, いる。緩和時間 も とほぼ比例関係にあるから, 同じ温度
我々が普段目にしている物質の多くはガラスと呼べそうだ。 依存性を示すと思ってよい。温度がたかだか半分程度に低ほとんどのプラスチック用品は高分子ガラスだし, 歯磨き 下するだけで, は 15 桁も増大している。この線を低温側
粉や化粧品などのペースト状の物質は, コロイドのガラス に外挿すれば, どこかで が発散しそうである。しかし, そだ。ひょっとするとゲルや砂山だって, ランダムで流れな こに達する前に, 観測者の辛抱の限界を超えてしまう。そいからガラスの一種と呼べるかもしれない。ではこのガラ の観測限界がガラス転移点 Tg だ。習慣上, が 1012 Pa∙s にスと液体とを分かつ「転移点」は存在するのか。実験的に 達する点を Tg の定義としている。Tg より低温で が発散
は, 液体が流れなくなるように見える温度を便宜上ガラス するということは, そこに真のガラス転移点があることを転移点 Tg と呼んでいる。しかし, これはとても転移点の真 意味している。実際に,
はの定義と呼べるものではない。測定者の辛抱できる限界を
exp(1)超える温度を転移点と呼んでいるだけだ。実際, 観測時間を長くすれば, Tg は低温側へシフトしてしまう。そもそも という Vogel-Fulcher 則と呼ばれる経験式でフィットする「流れない」状態とは正確には何を意味しているのか。我々 ことができる。驚くべきことに, ここでフィットパラメー
は無意識のうちに, 固体は流れないが液体は流れると考え タとして出てきた T0 は多くの液体でエントロピーの議論ている。その理由は, 固体には「規則的な分子配置」とい で述べたカウツマン温度 TK とほぼ一致することが知られう秩序があるのに対して, 液体の分子配置は無秩序だから ている。俄然, TK で真のガラス転移が起こると思い込みただ。しかし, この論理でいくと, ガラスは流れる物質に分類 くなるではないか。この の増大は, 分子どうしの交通渋滞されてしまう。本当だろうか? この疑問に対する答えとして, 二つ考えられる。一つは, 「我々がガラスと呼んでいるのは, 実は粘性が極端に大きい液体にすぎない」という答
え, もう一つが, 「ガラス転移という新しい種類の相転移が起きている」という答えだ。後者は, この真の転移点を境にして, その温度より高温では系は液体だが, 低温では真のガラス状態が現れるという考え方である。ガラス状態における分子配置は, 液体と全く同じようにランダムにしか見えないのだが, それは我々が愚鈍だからであって, 実際には, 「アモルファス秩序」とも言うべき, まるで形容矛盾のような秩序が隠れているのだ。どちらの答えが正しいか
我々はまだ知らない。この真のガラス転移点が「ある」派と「ない」派に分かれて活発な議論が交わされているのが現状である。筆者を含め多くの研究者は, 今のところ「ある」派である。その理由を理解するために, ガラス転移を
もう少し正確に説明しよう。まず Fig.1(a) は, 液体のエン
トロピーの温度変化を模式的に描いたものだ。液体を十分にゆっくり冷却すれば, 温度が融点 Tm になったところで一次転移を経て結晶となる筈だ。しかし, 系を急冷すると,結晶化をし損なって過冷却状態となる。液体のエントロピ
ーは結晶に比べて, 急激に減少する点に注目しよう。温度が下がると分子運動も遅くなり, なかなか熱平衡状態に行き着かなくなってしまう。さらに温度が下がると, やがて Fig.1 (a) T dependence of the entropy for a liquid. (b)
液体の緩和時間が, 急冷の時間よりも長くなり, 我々が観 Temperature dependence of the viscosities for several glass
測できる時間範囲では液体は凍結しているように見えて, formers1). (c) F(k,t) of the Lennard-Jones potential liquid
エントロピーが変化しなくなる。これが Fig.1(a)の一点鎖 obtained by simulation. (d) Dynamical heterogeneity in a
線で表した状態である。我々がガラスと呼んでいるのはこ soft-core potential liquid.

1. はじめにNetsu Sokutei 42 (4) 2015136

ガラス転移理論の最近の発展が原因だ。つまり微視的なレベルでも遅い緩和が見える筈
∗〉〈である。Fig.1(c)は,,で定義される密度相関関数の振る舞いを示したものである。N は分子数,は, 密度揺らぎを波数 k にフーリエ変換した量である。波長2π/ は分子サイズ程度に選んである。温度をわずかに変化させただけで, 急激に緩和時間 が成長していることがわかる。また固体的な物質に特徴的な階段上の緩和を示しているのも特徴だ。分子配置には何も変化がないのに, ダイナミクスが激変するのはなぜか。その遅いダイ
ナミクスの背後には, 何らかの新しい物理が隠れている筈である。もし何もなければ, いくら低温で分子運動が遅くなっても, や はたかだか exp[A/T] 程度の増大しかしないからだ。気液相転移などの二次相転移だったら話は簡単である。二次相転移は, 比熱や磁化率などの感受率が発散する。感受率は揺らぎの分散で書けるから, このことは揺らぎが空間的に発散的に成長することを意味する。実際に,転移点の近傍では, 揺らぎの相関長は, ベキ的に発散することが知られている。これが二次相転移点近傍でダイナミ
クスが遅くなる理由である。2) では, ガラス転移の場合はどうか。ガラス転移点近傍の分子の配置は, 液体と少しも異なるように見えない。ところが, 分子配置の時間変化を観測すると様相が一変する。これをシミュレーションで描いたのが Fig.1(d)だ。相対的に大きく動く分子を大きな球でハイライトしてある。動く分子がクラスターを作って集団運動をしている様子が分かるだろう。その様子は二次転移点で現れる揺らぎのフラクタル状のパターンにそっくりである。この集団運動は動的不均一性と呼ばれている。3)
これこそ, 分子運動を遅くする張本人に違いないという確
信が生まれ, 現在は, その相関長の研究が盛んに行われている。ただし, 分子運動を通して相関長が観測されたからと言って, ガラス転移が動的な現象である, と早合点してはいけない。この相関長は, ランダムな分子配置の中に埋め込まれているにも関わらず, 我々がその起源となる静的な「ランダム秩序」を検出する手段を持っていないだけかもしれないからだ。

∙ ,

μ

(2)

で記述できる。ここで, J は拡散流を表す。第二式は拡散流が, 化学ポテンシャルμの勾配により引き起こされることlnを表している。分子間の相互作用がなければ,だから, 式(2) はよく知られたフィックの法則(3)となる。D は拡散係数だ。分子間の相互作用がある場合は,μは自由エネルギー を用いてμ ∂ / ∂ρと書けるが, その
は
ln

′

(4)
と書くことができる。6) ここで c(r) は実効的な粒子間相互作用で直接相関関数と呼ばれる量だ。これを式(2) に代入
すると, 波数表示で∙(5)を得る。ここで1/S k
1は構造因子と呼ばれている量だ。この方程式は非線形方程式だからとても厳密には解くことができない。我々は密度相関関数 F(k, t) を求めたいのだが, 式(5) の右辺から, ∗ をかけて平均をしても,∗右辺第二項から, 〈0 〉のような 3 体の相関関数がでてきてしまうのだ。ならば, この 3 体相関関数の方程式を立てれば良いだろう。すると今度は 4 体の相関関数が出てきてしまう。このようにいつまで行っても高次の相関が芋づる式に出てきてしまい, きりがない。そこでエイヤッと, この連鎖を断ち切ろうというのが MCT の骨子である。何のことはない, 4 次の相関関数が出てきたら,∗∗∗〈0 ∗ 0 〉 〈0 〉〈0 〉と 2 体相関関数の積で近似してしまうだけのことである。
きちんとした導出はもう少し複雑だが, 結果だけを書くと
∂
,
∂
, ′
,
,
∂t
∂t′
(6)
となる。ここで, M(k, t) は記憶関数と呼ばれ,

2. ガラス転移の理論

,

2

,

,

,

(7)
ガラス転移の理論研究の現状は, 実は, 平均場理論のレ
,
∙
∙
は相互作
ベルの理解すらおぼつかないという体たらくである。平均 と書かれる。ここで場理論とは, 熱揺らぎを無視した一番単純な理論のことだ。 用による, 密度モードのカップリング（これが「モード結だからと言って理論研究が停滞していたわけではない。む 合」理論の名前の由来だ）を表している。この導出からもしろ, 最近の研究の進展は目まぐるしいほどで, もう少し 明らかなように, MCT は非常に荒い近似理論である。にもで平均場理論が完成しそうだというところまで来ているの かかわらず, あまり低温でなければ, 過冷却液体の実験だ。ここでは, このガラス転移の平均場描像を概観してみ やシミュレーションの結果を説明することに成功している。4) 特に, Fig.1(c)の F(k, t) の振る舞いは, ほぼ定量的によう。
説明できる。また, この計算から得られる緩和時間τはある
温度 Td で
2.1 ダイナミクスからのアプローチ||(8)
まずは, ガラス転移の平均場理論のダイナミクス版とものように,ベキ的に発散する。
ここで γは MCT から計算さいうべきモード結合理論（mode-coupling theory, MCT）の紹介から始めよう。MCT は, ボルツマン方程式に始まる分子 れる指数である。この突然の運動の凍結は, MCT 転移と呼運動論の直系の子孫である。液体のような高密度流体では, ばれている。ところが, これは見かけの転移である。Td は,
分子運動は, 瞬間的な衝突の足し算では記述できず, 混雑 実験的なガラス転移点 Tg や, 理想ガラス転移点 TK よりもした分子の集団がワサワサと拡散していると捕らえる必要 かなり高いのだ。実際に式(8) を Fig.1(b)にある粘性係数にがある。これらの多体相関を自己無撞着に取り扱う手法と 重ねると, 高温側の最初の数桁まではよくフィットできる
して, MCT が誕生したのは 1980 年代のことである。4) こ が, 肝心の低温側でのダイナミクスは全く説明できない。
れと並行して, Kirkpatrick らは MCT がスピンガラスの平 そ も そ も , 式 (8) は ベ キ 的 な 関 数 だ か ら , 式 (1) の均場模型のダイナミクスと等価であることに気付き, スピ Vogel-Fulcher 則と一致しようがない。これだけ読むと,ンガラス理論と MCT は融合することとなるのである。5) MCT は到底ガラス転移の本質を捉えていないように見え
まず, 液体の運動はミクロなスケールでは局所的な密度場 る。そもそも, なぜ MCT がガラス転移の平均場理論と呼ばれているのか。その理由は次節で明らかになる。
, に対する拡散方程式
Netsu Sokutei 42 (4) 2015
137

解

説

2.2 熱力学からのアプローチ
ガラス転移を熱力学的な相転移として記述する試みは,
歴史的には, MCT とは独立に進んできた。その平均場理論の本質は, ガラスを, 分子配置により決まるエネルギーのでこぼこの山と谷の中に系が落ち込んだ状態と見なすことである。これをエネルギーランドスケープ描像と呼ぶ。7)学と同様, f が小さいほど安定する傾向と が大きいほど安定する傾向がちょうど釣り合う条件を用いて, 以下のように評価できる

exp∗∗

全自由エネルギーはこの対数をとればよい。 ∗ は, 式(10)の指数の肩が極大を取る値で,
(12)

∗

で決まる。
を f の関数として表したのが図 3(a)である。
ランドスケープには下限があることに対応して, f が最小値を取るところで0となる。また, f の最大値以上
では, 液体はランドスケープの山頂よりも高いところにいるから,も存在しない。式(12)から, 図 3(a) の曲線の傾きが 1/T となるような点として, f とが決まる。これを逆に解いて, 温度 T の関数としてを描いたのがFig.3(b)である。高温で有限だったは, 温度低下と共に減少し, 有限の温度で消滅している。この図は, Fig.1(a)における液体と結晶のエントロピー差の温度依存性とそっくりだろう。つまり, このがゼロとなる温度こそ TK ということだ。

Fig.2 Schematic drawing of the energy landscape.
Fig.2 は, エネルギーランドスケープをスケッチ風に描いたものである。横軸は分子配置を表している。N 個の分子を空間上に配置する方法は無数にあるが, 軸上の一点は,その配置のひとつに対応していると思えばよい。縦軸はその配置ごとのポテンシャルエネルギーを表している。このエネルギーの山と谷の中を液体は運動している。温度が高ければ系は山をやすやすと乗り越え自由に配置空間内を運動できるだろう。しかし温度が低くなると, 系はエネルギ
ーの山に引っかかり始め運動が遅くなる。温度がさらに低くなれば系はもはや山を乗り越えることができなくなり,運動が完全に凍結する。これがガラス転移であると考えるのである。もし, 熱揺らぎが小さく（これが平均場理論の肝である）, 系がエネルギーの谷に落ち込んで, そこから抜け出ることがほとんどできなくなっていれば, その谷間の中を一つの熱力学的な「相」と見なして差し支えないだろう。この「相」は, 液体相でも結晶相でもなく, 「ランダムな相の中のどれか一つ」としか呼びようがないものである。
ここで登場した特定の「相」を状態 a とし, この相の自由エネルギーを Nfa と書こう。N は液体全体の分子数で, fa は,一分子当りの自由エネルギーを表す。この相 a の分配関数は,expと書ける。系全体の分配関数は, 実際に液体が取り得る相についての総和を取ればよく
exp

9をエネルギー準位 Ea に読と書けるはずだ。この式は,み替えれば統計力学の教科書に出てくる分配関数とそっくりである。そこで教科書でやるように, ラベル a についての足し算を, 準位 についての足し算に置き換えると上
の式は,
W
exp

exp
/

(11)

10

と書ける。ここでWは, 自由エネルギーが Nf になるよ
うな, 分子配置の全状態数である。最後の等式で定義した
lnW
はいわゆるエントロピーである。ただし, は, 振動などの細かい自由度について既に均してしまった自由エネルギーf の関数であるため, 特に配置エントロピーと呼ばれている。式(10) の積分は, 通常の統計力

Fig.3
T.

(a)

as a function of f. (b)as a function of

以上のように厳密な計算により, ガラス転移の平均場描像を説明することが出来た・・・というのは少し言いすぎで, 実は, ここでの計算は過冷却液体ではなく, スピンガ
ラスの平均場模型に対して行ったものである。7) 気液相転移の場合と同様, 難しい多体現象を考えるとき, 何事によらずハミルトニアンが簡単なスピン模型を考えるのは統計力学の常道である。ただしガラス転移の場合のスピン模型
はちょっと風変わりで, ハミルトニアンは,13という形をしている。ここで はスピン変数である。 がスピン間の相互作用の強さを表していて, , , の組み合
わせごとにランダムな値をとる。また数学的な便利さのため, はイジング型ではなく, 連続変数としてある。ただし,スピンが発散するようなことがないように∑という拘束をかけておく。この模型は, 3 つのスピンがグループになって相互作用しているので, p3 スピン模型と呼ばれている。ずいぶん非現実的な数理モデルだが, その振る舞
いが, 現実のガラスの観測事実をよく説明できるために,ガラス転移のプロトタイプと信じられている。この模型の解析により, ガラス転移点 TK で何が起こっているのか明らかになってきた。通常の一次転移では, 転移点でエントロ
ピーが不連続に変化し, 潜熱が発生する。二次転移では, エントロピーの飛びはないが, 代わりに感受率が発散する。
いずれの場合も, 磁化や密度などの秩序変数が, 転移点を境に, はっきりと変化するのが特徴だ。これに対して, ガラ

Netsu Sokutei 42 (4) 2015138

ガラス転移理論の最近の発展
ス転移の場合は, エントロピーが Fig.3(b)のように TK で連続的に消滅するが, その傾きの変化は不連続である点が特徴だ。TK で感受率の発散などもない。では, 秩序変数は何か。普通の磁化などは不適当である。スピン配置は TK 以下でもランダムのままで, 磁化もゼロのままだからだ。我々はもう少し賢くならなくてはならない。今, 全く同じ Jijk
を持つサンプルを, 適当に二つ用意する。ある時刻において, 1 つ目のサンプルのスピン配置が, , , ⋯ , 2 つ目のサンプルのスピン配置が,,, ⋯ だったとする。この時, スピン配置の内積
114を「重なり関数」と定義する。この こそ, 「賢い秩序変数」
である。配置が完全にランダムであれば0 となるが,見た目にはランダムであっても, 1 つ目のサンプルと 2 つ目のサンプルに相関があるならば,0となる筈だからだ。
真のガラス転移点 TK で起きているのは, この0から0への転移なのだ。この転移はスピンガラスの平均場理論において,「レプリカ対称性の破れ」と呼ばれている。7)
ここまでは熱力学の議論であったが, これと MCT には深い関係がある。もう一度, Fig.3(b)を見て欲しい。この図には0となる温度 TK の他に,が突如現れる温度がある。この温度は, 液体がエネルギーランドスケープの山に引っかかり始め動けなくなる温度である。これが前節で登場した MCT 転移点 Td に他ならないのだ。これは,式(13) で定義した p3 スピン模型に対して厳密に証明されている。実際に, この模型のダイナミクスの方程式を書き
下すと,
∂
∂
∑
(15)
7)
となる。 右辺第 1 項のは, スピンの大きさを発散させないための拘束因子である。右辺第二項が相互作用項である。7) この式は, MCT の導出の出発点となった式(5) と等価な形をしている。そのため, 前節と同じ手続きを経る
ことによって, 式(6)と数学的に等価なスピンの相関関数に対する運動方程式が導出できるのである。スピンの方が,波数依存性がないから少し簡単だが, そんなことは本質的ではない。このスピンの MCT から導かれる Td が, Fig.3(b)が突然表れる温度とぴたりと一致することが導けでるのだ。
ここまで, 液体のガラスのプロトタイプとも言うべきス
ピン模型を引き合いにして議論したが, このことに不満を感じはしなかったろうか。私は感じる。スピンのハミルト
ニアン, 式(13) は液体のそれと似ても似つかないのに, 本
当にスピンで証明された結果が液体に当てはまるのだろうか。普通の相転移だったら良いだろう。なぜなら, 相転移では系に潜む対称性が本質であって, 個々の系の詳細など
重要ではないのだから。2) しかし, ガラス転移の場合は,話はそれほど簡単ではない。現在, スピンガラス模型で使われた様々な理論や考え方を, 液体に移植する試みが活発に行われている。熱力学については, レプリカ液体論という体系が作られ, TK や配置エントロピーなどが計算できるようになった。8) ダイナミクスについては既に MCT がある。しかし困ったことに, レプリカ液体論と MCT が整合していないように見えるのだ。9) これらは多くの研究者が悩んでいる深刻な問題だ。ところが, この原稿を執筆している間にフランスのグループがその解決の道筋をつけたというニュースが飛び込んできた。10) これが本当であれば, い

よいよ「液体の」平均場理論が完成したことになり, 前時
代的な状況にあったガラス転移研究が大きく進展することになる。

2.3 平均場理論から現実のガラスへ
前節まで, 熱揺らぎがない平均場理論の議論をしてきたが, 現実の 3 次元の系は熱揺らぎのために, 山と山に分け隔てられた「相」の定義があやふやになる。この揺らぎのために, 系は, 状態から状態へと熱活性過程により山を乗り越えていくことが出来るに違いない。またこれを, 実空間で眺めれば, 無数の相の共存状態が見えるに違いない。

Fig.4 A schematic cartoon of the mosaics.
この様子を漫画にしたのが, Fig.4 である。この漫画で描いた点線が, 相の境界である。異なる相がパッチワークのよ
うに液体を埋め尽くしていることが分かる。そこでこの相の一つ一つをモザイクと呼ぼう。問題は, 異なる相といってもどれも単なるランダムな配置にしか見えず, その境界がさっぱり分からない点である。ともかく, このような無数の相があると信じることにすると, ガラス転移点近傍ではモザイク内の分子達が協同的にガサガサと再配置しているように見える筈だ。11) そのため, モザイクのことを協同
再配置領域と呼ぶこともある。このモザイクの一つ一つを,結晶化における核みたいなものだと思えば, 一次転移の核生成と似た議論を用いて, ガラスのダイナミクスを理解できる。この考え方に立って平均場理論を有限次元の系に拡張した理論を, ランダム一次相転移（RFOT）理論と呼ぶ。12) 核生成では, 結晶と液体の間の自由エネルギー差が駆動力, 表面エネルギーが結晶化を妨げる力で, その競合
で結晶化のダイナミクスが決まるのであった。ガラス転移の場合は, 液体中のモザイクは, 一つの相に留まっていた
くはなく, できることなら, 取り得る多くの状態へと変移したがっている。つまり, ガラス転移における核生成の駆
動力は配置エントロピーであると考えるのが自然である。
一方, 一つのモザイクは異なるモザイクと隣り合わせにある分だけ表面エネルギーを損しているだろう。そこで, この核生成のエネルギーバランスは, ふつうの核生成の議論
に倣って,
Δ
16
と書くことができる。ここで, はモザイクの大きさで, σは表面エネルギーである。モザイクがどんな形か分からないので, サイズ依存性を表す指数を と置いた。このΔ は,Netsu Sokutei 42 (4) 2015
139∗

で極大となり, その値はΔ ∗/∗となる。熱活性化過程による遷移の時間は,/と与えられるので, ~exp
を~exp Δ ∗ /

解
得る。もし

説

3/2であれば

~exp

(17)

を得る。これは現象論として良く知られる Adam-Gibbs 関
11)
係式に他ならない。
さらに, 前節で議論したように は T= TK でゼロになるので,と展開すると,~exp"

度で Td と TK は一点で交わり, そこでガラス転移点は消滅してしまう。この図は, 気液相転移の相図によく似ている。
共存線が TK に, スピノーダル線が Td に対応しているわけ
である。

(18)

を得る。これは式(1) の VF 式に他ならない！
MCT と RFOT から浮かび上がってくるガラス転移の全貌は以下のようなものである。まず, 高温側から温度を下げていくと, Td 付近で, MCT が記述する遅いダイナミクスが観測される。しかし, 実際には熱揺らぎのために, はっき
りした MCT 転移は起こらず, 単なるクロスオーバーとなる。Td より低温では, 遅いダイナミクスの主役は熱活性化過程となる。その駆動力は, 配置エントロピーである。
それを特徴付ける相関長は, モザイクの大きさで ∗ である。これが動的不均一性の正体だろう。さらに温度が下がるとも減少するので, 運動は加速度的に遅くなる。やがてが消失する TK で完全に運動は凍結し, さらにそこでは系の対称性（レプリカ対称性）が破れて, 真のガラス転移が起こる。
ここまで読むと, MCT と RFOT でガラス転移の本質はもうすっかり分かってしまった気になってしまう。しかし話はそれほど簡単ではない。ガラス転移の理論は,
MCT+RFOT 以外にもたくさんある。そして, それらの理論
も実験結果を同じくらいよく説明できるのだ。そのような理論の一つに, ガラス転移を, 安定な秩序相がどこかにあるはずなのに, たまたま空間が 3 次元だったために秩序化するのに失敗してしまった状態と見なす, フラストレーション描像がある。13) この他にも, ガラス転移は熱力学転移
ではなく, 純粋に動的な特異性の産物だとみなす, Dynamic
facilitation と呼ばれる描像もある。14) これらはいずれも TKの存在を仮定していないにも関わらず, ガラス転移の様々 Fig.5 (a) Schematic phase diagram of the randomly pinned
なデータを説明できてしまうのだ。
spin model, 15) (b) The phase diagram obtained for the numericalsimulation for the Lennard-Jones liquid. 16)

3. 真のガラス転移を探る
相容れない複数の理論が提案された場合に, 真偽の裁決をするのはもちろん自然であり観測データである。問題は,実験やシミュレーションで TK を観測することがどうしてもできないことだ。しかし最近, この TK を現実に観測するアイデアが提案された。液体中に不純物を入れることで,シミュレーションできる温度領域まで TK を引き上げるというものである。方法はいたって簡単で, まず, ある温度 T
で液体を平衡化するまで十分に待つ。次に液体中の分子の何割かを無作為に選んでその運動を瞬間的に凍結（pin 止め）させる。最後に, 凍結されていない分子を運動させて,熱力学やダイナミクスを観測する。この操作を何度も繰り返し, サンプルの平均を取る, というものだ。pin 止めされた分子があるとガラス化しやすくなり, TK も高くなるはずである。かくして温度を高温に保ったまま, ガラス転移を観測することができるのである。この方法の肝は, pin 止め
した分子は平衡状態の液体の配置から選んでいるため, 得られるガラス状態は熱力学的には平衡状態であることだ。
このアイデアは, まず式(13) で定義されたスピンガラスの平均場模型に対して厳密に解析された。15) この結果得られるガラス転移の相図を,温度 T と不純物濃度 c の関数としてスケッチ風に描いたのが Fig.5(a)である。平均場理論が予想するようにガラス転移には, 二つの転移点 Td と TK があるが, それらが c の増加と共に上昇することがわかる。
同時に Td と TK の間の距離が縮まっていく。やがてある濃この結果を, 実際に観測することができれば, 平均場理論の強力な証拠になるだろう。そこで我々は, Lennard-Jones
ポテンシャル系を用いて pinning 系のシミュレーションを行った。16) 理論を検証するために, 我々が確かめなくてはいけないことは, (1) TK において, 式(14) で定義した重なり関数 Q がゼロから有限の値に転移すること, (2) TK において, がゼロになること, そして, (3) ある濃度で Td がTK と交わること, の三点である。我々は, レプリカ交換法
を用いて大規模なシミュレーションを行い, 膨大な量の分子の配置を生成し, これから, まず(1) の Q を計算した。
(2) については, 内部エネルギーを温度で積分することで全エントロピーを求める。 は, 全エントロピーから, 不必
要な振動成分を差し引けば計算できる。(3) の Td について
は, 系のエネルギーランドスケープで不安定な極値が消滅する温度として求めた。このようにして求めた TK と Td を相図で表したのが, Fig.5(b)である。我々が計算できるのはせいぜい T = 0.45 よりも上の温度だから, それよりも高温の状態しかプロットしていない。Q が不連続に変化する点(TK) を◆で,0 となる点（TK）を●で, そして Td を■で表している。高温かつ高濃度側での点が白抜きになっているのは, Q の変化や が消える様子が曖昧だったため,転移の判別ができなかった部分である。この結果をみると,
低温, 低濃度側では二つの方法（Q と ）で計算した TK が一致すること, TK と Td が c の増加とともに上昇すること,
Netsu Sokutei 42 (4) 2015140

ガラス転移理論の最近の発展
そして TK と Td が一つの点で交差していることがわかる。
さらに興味深いことはこの交差点を境にそれより高温高濃度側で, 転移の判別ができなくなり, Q から決めた TK とから決めた TK の二つの線が大きく分岐していることだ。これらの結果はすべて平均場理論の予想を支持している。また, 交差点より高温側で, 2 つの TK が分岐する様子は, まる
で気液相転移における Widom 曲線（臨界点より高温側にも残る熱力学異常の痕跡のこと）のようである。17) ここまで
平均場理論の主張をエレガントに示したシミュレーション
結果はないのではないだろうか。何よりも重要な点は, 世界で初めて配置エントロピーがゼロになる, 真のガラス状
態を計算機の中で作ることに成功したことだ。とはいうもののこれで一件落着というわけではない。気液相転移から類推できるように, Td と TK の交差点は一種の臨界点であろう。平均場理論では既にこれが, 「ランダム磁場イジング
模型」というユニバーサリティクラスに属する臨界点であることが予言されている。15) 我々の研究は今のところ大雑把に相図を描くところまでできたが, この臨界点の解明は今後の大きな目標である。この研究を機に, 有限次元系のガラス転移の定量的な解析が進むことが望まれる。

4. おわりにこの解説では, できるだけ分かりやすくガラス転移の理論研究の現状を説明することを試みた。分かり易さを強調したために, 正確さはかなり犠牲にしていることを了解願いたい。また, 数多あるガラス転移理論の中で, 平均場理論とそれに毛が生えた程度の RFOT に議論を絞った。もし今,どの理論が正しいと思うかを投票で選ぶなら, RFOT が一番人気だと思ったからだ。しかし科学は民主主義ではないから, これを以って RFOT が正しいとする理由にはならない。
また, p3 スピンや単分子原子液体など異常なまでに単純化された模型ばかりを題材にして説明をしてきた。現実のガラス形成物質は内部自由度や分子の形状をもち, はるかに複雑である。その複雑さが原因で, 普遍的に観測できる現象も多い。たとえば, ここでは言及すらしなかった,
Johari-Goldstein 過程（遅いベータ緩和）やアレニウス的な緩和からのずれの度合いを表すフラジリティなどは, ガラス転移の本質と直結した重要な概念だが, ここで取り上げることはしなかった。もうひとつ、ここで全く言及しなかった現象に, 粉体のガラス転移とも言うべきジャミング転移がある。かつては, ジャミング転移はガラス転移の温度ゼロの極限であると思われていたが, 最近急激に理論が整備されて, 平均場理論の範囲で非常に精密に記述できるよ
うになってきた。18,19)
この解説で扱った事柄のより正確で完全な説明は, 最近出版されたいくつかの優れたレビュー論文に書かれている。
特に, 文献 20) や 21) は, RFOT とほかの理論的シナリオについても概括されており, またそれらを批判的に比較しているのでお勧めである。ガラス転移については確かなことは何も分かっていないのだが, 研究の進展は目まぐるし
く, この数年だけを見ても, すっかり様変わりをしているように思われる。今後数年以内に, ここで解説した事柄がすべて間違いだったことがわかった, ということも十分に考えられる。臨界現象のときの繰り込み群法の発見のような, 革命的なパラダイムシフトが起これば面白いだろうが,筆者は, そのようなことはガラス転移の場合にはないような気がしてならない。

謝辞

この解説の最後の章は, Walter Kob 氏との共同研究の成果である。またその研究にあたっては, 科研費 24340098,
25103005, 25000002, そして学振の Core-to-Core Program から支援を受けている。本稿の図面を作成した川崎猛史氏に
感謝する。

文

献

1) C. A. Angell, J. Non-Cryst. Solids 102: 205-221(1988).
2) 西森秀稔, 「相転移・臨界現象の統計物理学」培風館
3) R. Yamamoto and A. Onuki, Phys. Rev. E 58, 3515-3529(1998).
4) W. Gotze, "Complex Dynamics of Glass-Forming Liquids",
Oxford University Press, Oxford (2009).
5) T. R. Kirkpatrick and D. Thirumalai, Phys. Rev. B 36,
5388-5397 (1987).
6) J. P. Hansen and I. R. McDonald. "Theory of simple
liquids", Academic Press, Second Edition edition (1986).
7) T. Castellani and A. Cavagna, J. Stat. Mech., P05012(2005).
8) M. Mezard and G. Parisi, Phys. Rev. Lett. 82 747-750(1999).
9) A. Ikeda and K. Miyazaki, Phys. Rev. Lett. 104, 255704(2010).
10) T. Maimbourg, J. Kurchan, and F. Zamponi, arXiv:
1507.03421.
11) G. Adam and J. H. Gibbs, J. Chem. Phys. 43 139-146(1965).
12) T. R. Kirkpatrick, D. Thirumalai, and P. G. Wolynes, Phys.
Rev. A 40, 1045-1054 (1989).
13) G. Tarjus, S. A. Kivelson, Z. Nussinov, and P. Viot, J.
Phys.: Condens. Matter 17 R1143-R1182 (2005).
14) D. Chandler and J. Garrahan, Annu, Rev. Phys. Chem. 61,
191-217 (2010).
15) C. Cammarota and G. Biroli, Proc. of the Natl. Acad. of Sci.
109 (23), 850-8855 (2012).
16) M. Ozawa, W. Kob, A. Ikeda, and K. Miyazaki, Proc. of
the Natl. Acad. of Sci. 112 (22), 6914-6919 (2015).
17) B. Widom, J. Chem. Phys. 43 (11), 3898-3905 (1965).
18) G. Parisi and F. Zamponi, Rev. Mod. Phys. 82, 789-845(2010).
19) P. Charbonneau, J. Kurchan, G. Parisi, P. Urbani, and F.
Zamponi, Nature Comm. 5, 3725 (2014).
20) A. Cavagna, Phys. Rep. 476, 51-124 (2009).
21) L. Berthier and G. Biroli, Rev. Mod. Phys. 83, 587-645(2011).

Netsu Sokutei 42 (4) 2015141