理想純物質とその物質状態のための統計熱力学

<a id="gesamt"></a>

全体要約

モデル3は、固体・液体・気体という三つの物質状態とその相転移を、分子結合による単一の結合像から導出する 定性的な統計力学的理想物質モデルである。基本原理は、内部エネルギーに対して三つのエネルギー貯蔵を考慮する ことにある。すなわち、分子結合中の振動エネルギー、分子の並進エネルギー、そして内部の密度ゆらぎを均すための 内部物質流の理想化された形としてのミセル回転における協同的運動エネルギーである。この文脈における重要な基本量 は平均配位数 z̄ = z_max·x であり、これは協同的二状態結合(イジング型、協同性 K = J·z_max、破断状態の縮退 G = z_max)を通じて自己無撞着に定まる。ここから、ただ一つのエネルギーパラメータ ε で較正された共通の τ 尺度上に 三つの転移が導かれる。融解はフォノン/剛性のしきい値として(z̄ が剛性に達する)、沸騰は臨界点における 連結性(パーコレーション)の崩壊として、そして任意選択の固体(1)→固体(2)転移(FCC↔BCC)が 1.0 k_B というパラメータ不要の c_v 段差を伴って現れる。転移は ln G > 2 + 4ε/K となった時点で一次転移になる。 ここで破断状態の縮退 G は破断しうる隣接接触の数(G = z_max、最大配位数)、結合エネルギー ε は隣接接触ごとに 得られるエネルギー、協同性 K = J·z_max はイジング型結合強度 J と z_max の積である。したがって、接触の多様性 G が大きいほど、また結合エネルギー ε に対して協同性 K が強いほど、転移はより早く一次転移になる。

本稿は七つの公理(第1章)から矛盾分析(第2章)と静的な漫画コマによる絵画的記述(第3章)を経て、統計力学的 骨組み — 協同的結合、分配関数、熱容量(第4–5章)— へ、さらに三次元的な数学的もっともらしさ(第6章)へ、 そして最後に較正と状態方程式(第7章)へと導く。これらの考えは、一次融解の完全な三次元 c_v 走査(Δz̄、 潜熱 ΔU、ピークの移動)、フォノン剛性 θ(格子) に対するアインシュタインばねモデル、および クラウジウス–クラペイロンによる固体(1)→固体(2)変態線によって補完される。さらに θ(格子) は デバイ型フォノン状態密度へと精密化され、次いで固体(1)→固体(2)変態が(鉄で較正された) (τ, p) 図における線としても導入され、ファン・デル・ワールス既定値の隣に分離されたレドリッヒ–クォン枝が 追加される。その際 θ(格子) はレナード–ジョーンズ完全格子(FCC は安定、理想 BCC はせん断不安定)で、 固体(1)→固体(2)のトポロジーはチタン(第二の確認例)とアルゴン(陰性対照)で検証される。第8章は 本モデルを静的熱力学の十一の確立したモデル観の中に位置づける。付録は文献出典(著者順)と AI 支援 (推定革新寄与率、推定費用、各開発段階の所要時間)を記録し、補助部は読み方、残された批判点、展望をまとめる。

核心的主張。 ただ一つの協同的結合尺度 z̄ が、本来別々の三つの伝統の系統 — 格子気体/イジング像(協同的二状態転移)、パーコレーション像(連結性 → 沸騰)、 アインシュタイン/デバイのフォノン像(剛性 → 融解)— を、意識的にファン・デル・ワールス型の理想像の 内部で、整合的で定性的に完全な相図へと連結する。しかも理論の建物は、容易に思い描ける動的な像から自然に 導出できる。精度ではなく絵画性こそがこのモデルの本来の強みであり、その科学的モデル価値を特徴づける。

目次

• 全体要約
• 1 — 標語名を付した公理(統合版)
• 2 — 矛盾と欠落の分析
• 3 — 絵画的記述
• 3.1 三つの個別相
• 3.2 相転移
• 3.3 なぜ液体中の音速は T とともに下がるのか
• 3.4 ブラウン運動はいかにして生じるか — 鋭くしたミセル概念とともに
• 3.5 なぜミセル境界層で反応が速くなるのか
• 3.6 三次元の液体窓
• 4 — 統計力学的モデル
• 4.1 自由度と変数
• 4.2 協同的二状態結合(イジング型、パラメータ J)
• 4.3 エネルギー収支 — 協同項を含む
• 5 — 分配関数と熱容量
• 5.1 因子化された分配関数
• 5.2 相ごとの熱容量(解析的)
• 5.3 協同項を**含む**数値的 c_v(τ) 曲線 — 潜熱の定量化
• 5.4 結果としての連結性パーコレーションしきい値(変更なく確認) - 5.4.a 二つの異なるしきい値を明確に分離
• 5.5 融解転移の次数
• 5.6 パーコレーションしきい値の有限サイズ研究
• 5.7 一次融解の完全な三次元 c_v(τ) 走査
• 6 — 三次元解釈の数学的もっともらしさ
• 6.1 回転自由度の 2D → 3D 変化
• 6.2 各種結晶格子の影響 — 識別特徴としての連結性
• 6.3 追加の転移 固体(1) → 固体(2)
• 6.4 考察:液体の水の音響異常
• 6.5 代替案:ミセル回転に代わるフラクタル渦場
• 6.6 微視的フォノン剛性 θ(格子):アインシュタインばねモデル
• 6.7 準安定性の数学的見積り(2D 三角 ↔ 正方、16b)
• 6.8 θ(格子) のデバイ型精密化
• 6.9 レナード–ジョーンズ完全格子から得る θ(格子)
• 7 — 較正と状態方程式
• 7.1 較正法(エネルギーパラメータ ε 一つ)— 調整済み
• 7.2 較正表
• 7.3 密度/体積軸:硬核ファン・デル・ワールス状態方程式
• 7.4 定量的な液/気較正:a → ε, b → D²
• 7.5 三重点を含む共通の (τ, v) 面
• 7.6 潜熱と融解エンタルピーの対比
• 7.7 エネルギー尺度の統一:ε ↔ ε_lg
• 7.8 z̄_c から得るパラメータ不要の三重点
• 7.9 三次元検算:単一の ε 尺度上の明瞭な液体窓
• 7.10 真の三次元ボンドパーコレーション走査(FCC/BCC)と窓端の感度
• 7.11 結合速度論の真の三次元 τ 走査:τ_b(3D) と単一 ε 尺度上の融解線
• 7.12 (τ, v) 面における固体(1)→固体(2)変態線
• 7.13 鉄で較正した (τ, p) 図における固体(1)→固体(2)変態線
• 7.14 分離した理想像の枝としてのレドリッヒ–クォン
• 7.15 他の多形における固体(1)→固体(2):チタンと希ガス対照
• 8 — 文献比較:静的熱力学のモデル観
• 8.1 一覧表
• 8.2 個々のモデル
• 8.3 モデル3との類似点と相違点
• 9 — 付録
• 9.1 文献
• 9.2 AI 支援
• 9.3 略号・記号集
• 9.4 計算スクリプトとテストカバレッジ

<a id="kap1"></a>

1 — 標語名を付した公理(統合版)

要約。 標語名を付した七つの公理が理想物質を定める。保存される粒子数、純粋に古典的な運動、協同的二状態結合、 そして秩序パラメータとしての平均配位数 z̄ である。これらは以後のすべての導出を支える。ドイツ語原文では標語に 同形異字の透かしが施されている(鍵は本文末)。

精密化した公理:

1. 自由飛行(補助公理)。 結合していないまま互いの傍を飛び過ぎる二分子の間には、摩擦も静電気力も 働かない。運動量とエネルギーは、衝突時の硬い 硬核 斥力と 結合 の生成・破断によってのみ変化する。この公理は もっぱら古典力学とエネルギー収支の観点で重要である。すなわち、分子レベルに散逸的な経路が存在しないことを 保証する — 運動エネルギーは決して「失われ」ず、三和音 の三つの貯蔵か結合エネルギーへ移るだけである。 (巨視的な摩擦・粘性は微視的な基本力としてではなく、スライドパズルと 亀裂 から*創発的に*生じる。)
2. 多数性。 理想物質はきわめて多数の同種分子から成る。巨視的性質は熱力学的極限(N → ∞)で統計的に記述される。
3. 結合。 幾何学的に隣り合う任意の二分子は一本の結合を担いうる。健全な結合はエネルギーを有限量 ε だけ下げ、 その破断はちょうど ε を要する。結合は {健全, 破断} の二状態系である。
4. 三和音。 内部エネルギーは古典的に三つの貯蔵へ配分される:(a) 分子の並進、(b) 健全な結合中の振動 (結合ごとの調和ばね)、(c) 剛体分子ブロックの回転(連結領域の集団的回転自由度)。
5. 硬核。 分子間距離が最小値 σ を下回ると硬い斥力が働く(固有体積)。結合 と合わせて実効的な対ポテンシャル (硬い殻 + 有限の結合井戸)が得られる。
6. 保存。 全運動量と全エネルギーは保存する。個々の分子はエネルギー的に中立であり、内部にエネルギーを溜め 込まない。エネルギーは結合(結合)と運動学的自由度(三和音)に宿る。
7. 亀裂。 結合は亀裂/治癒の原理に従って崩壊・生成する。欠けた隣接結合が多いほど → 破断確率は高く、結合した 隣接が多いほど → 治癒確率は高い。この速度論は協同的である。
8. 敏捷。 破断と治癒は速く低活性化で進む。結合部分系はどの温度でも局所平衡にあり、分子運動に対して断熱的に 分離してよい(因子化された分配関数の基礎)。

モデル化の姿勢。 幾何学的配位は理想化されたまま保たれ、変わるのは破断による配位数 z だけである。運動は スライドパズルのように生じる。空隙(硬核 の固有体積と 亀裂 の破断に由来)が移動し、ブロックが後を追って ずれる。回転しながら移動する小結晶という意識的に「誤った」像は、直観的基礎として保持される。

<a id="kap2"></a>

2 — 矛盾と欠落の分析

要約。 この分析はモデルを内的矛盾と欠落について体系的に検査し、作業状況を書き継ぐ。認識された弱点は名指し され、以後の章で解消されるか、あるいは正直にモデルの限界として印される。

批判点の一覧

1. 並進の二重計上(公理 3a 対 3c)。 原文は並進を二度挙げている(分子の並進*と*ブロックの並進)。独立した 並進貯蔵が二つあることはありえない。解消: ブロックの*並進*は、すでに数えた分子並進の集団モードにすぎない。 真に新しいのはブロックの回転自由度だけである。自由度収支(第4章)では並進はちょうど一度だけ数えられる。
2. 「エネルギー的に中立」(公理 5)対「結合がエネルギーを受け取る」(公理 2/3b)。 見かけの矛盾。 解消: 中立性は*物体としての分子*に関わる。エネルギーは分子と分子の間(結合 の中)と集団的運動モードに 宿り、分子内の緩衝器に宿ることは決してない。分子 ≠ 結合 — これで矛盾は解消する。
3. 引力の欠如(公理 4 は斥力しか挙げない)。 引力がなければ凝縮相は存在しない。解消: 引力は 結合(深さ ε の井戸)が、斥力は 硬核(σ における硬い殻)が担う。合わせて実効的な対ポテンシャルとなる。
4. 対抗項のないフィードバック(公理 6)。 「欠けた隣接が多いほど破断が増える」だけでは正のフィードバック (全崩壊まで亀裂が進行)である。解消: 治癒項(隣接が多いほど治癒が増える)と有限結合エネルギー ε(結合)の ボルツマン因子が 亀裂 を釣り合わせ、暴走ではなく*明確な*転移を生む。
5. 語られなかった時間尺度の分離(公理 7)。 「速く低活性化」は、そこから結合動力学と並進動力学の断熱分離を 導いて初めて使える。解消: 因子化 Z = Z_trans·Z_vib·Z_rot·Z_bond の前提として明示的に記した。
6. 幾何学の欠如と運動機構の欠如。 公理は配位格子も、流動を可能にするものも挙げていない。 解消(明示した追加仮定として): 結合占有が可変な固定参照配位幾何。正準的な二次元の場合として三角配位 (z_max = 6)。可動性はスライドパズル的空隙機構による。この二点だけが七公理を超える真の追加であり、第6章で 戦略的問題として印されている。
7. 潜在的な散逸の隙間。 版07は、結合していない分子が通過時にエネルギーを失うか否かを未決のままにしていた。 自由飛行 による解消: 微視的な摩擦・静電気力は存在しない。すべてのエネルギー流と運動量流はもっぱら 硬核 衝突と 結合 速度論を経る。これによりエネルギー収支(第4.3章)は閉じ、気体極限でのマクスウェル–ボルツマン統計は 厳密になる(弾性衝突の合間は自由な理想気体)。
8. 秩序パラメータのエネルギー的不整合。 版08は協同項 −(K/2)x² を*自由*エネルギー g(x) では(x(τ) を生成する ために)用いたが、内部エネルギー U/N では用いていなかった。これは熱力学的に不整合である。同じ分子場項は 両方に現れねばならない。解消(第4.3章): −(K/2)x² を U/N に反映した。これにより潜熱は測定可能なほど増加し (第5.3章)、固体プラトーと気体プラトーは変わらない。
9. 液/気較正が定性的にとどまる。 版09は密度軸(ファン・デル・ワールス状態方程式)を換算形でしか示していな かった。解消(第7.4章): 臨界点 (T_c, p_c) を通じて絶対的な沸騰線を計算し、T_b(p_atm) を実測値と対置する。 系統的なずれ(グルドベルク則)を名指しする。
10. 三重点のない分離した軸。 融解軸と沸騰軸は別々の計算であった。解消(第7.5章): 両者を一つの (τ, v) 平面に統合し、三重点を構成し、臨界点を印す。

<a id="kap3"></a>

3 — 絵画的記述

要約。 第3章は固体・液体・気体とその転移を(静的な漫画コマで)直観的に説明する。三つの個別相(3.1)、 再結晶/融解/沸騰/昇華の相転移と三重点・臨界点(3.2)、音響異常(3.3)、ブラウン運動と鋭くしたミセル概念 (3.4–3.5)、そして三次元の液体窓(3.6)である。 これらの記述に基づいて、異なる相における動的な分子的関係を心の目に思い描く助けとなるべき図解(漫画コマ)が 作られた。この思い描けることこそがモデルの本来の核心である。統計力学を通じて、このモデルが大がかりな数学なしでも 定性的な経験的傾向・観測を予測可能にすることをもっともらしく示せたとき、それは定性的かつ科学的に価値をもつ。

<a id="kap3-1"></a>

3.1 三つの個別相

固体。 ほとんどすべての隣接分子が分子結合で結ばれており(結合)、分子の配位数の平均 z は幾何学的最大値の 近くにある。分子網は空隙のない単一の剛体ブロックである(硬核 が密に詰まっている)— 液体のような内部スライド パズルは、この緻密な版では不可能である。唯一の運動は静止位置のまわりの震え(結合振動、三和音-b)であり、それは 隣接分子との分子結合によって決まる。 運動量とエネルギーは網の中で弾性的に往復する(保存、自由飛行:損失なし)。極限の場合(きわめて低い温度、理想的 充填)には理想結晶という理想状態が現れる。

液体。 協同的な 亀裂 が網を貫いた。結合が斑状に欠け、分子の配位数の平均 z は下がり、剛体ブロックは内部で多数の 小ブロックへ砕けた。一枚岩は瓦礫の山になった。これらのブロックは並進し、また回転する(三和音-c)。その間には 空隙が開き、そこを通って互いにずれ違う — ブラウン運動のスライドパズルが働いている。局所的な配位幾何は理想化された まま保たれ、変わるのは結合が実際に*いくつ*存在し形成されているかだけである。

気体。 亀裂 が勝った。結合はほぼ至る所で破断している(z → 0)。極端な場合、どの分子ももはや他の分子への分子 結合をもたない。ブロックは個々の分子へ分解した。分子は自由に飛び、硬核 衝突によってのみ結び付く(その合間は 自由飛行)。速度はマクスウェル–ボルツマンに従う(保存)。結合振動とブロック回転は消え、純粋な並進だけが残る (三和音-a)。極限の場合(高温、低圧)には理想気体という理想状態が現れる。

<a id="kap3-2"></a>

3.2 相転移

再結晶(固体(1)/固体(2))。 融解より下でなお、固体は剛性を失うことなくその格子を替えうる。異なる 配位数をもつ二つの理想充填(FCC z_max = 12 ↔ BCC z_max = 8)は異なる自由エンタルピー μ(τ, p) をもつ。 μ_FCC = μ_BCC のところで配位幾何が跳躍的に切り替わる(融解のように占有 z̄ が変わるだけではない)。見えるのは 小さな c_v 段差(Δc_v = z_max/4 の変化 = 1.0 k_B)と弱い潜在変態熱である。AI の主張:これは三次元でのみ 考えられる(離散的で量子化された連結性、第6.3章)。 二次元の場合には、真の再結晶、すなわち二つの固相間の転移をより容易に思い描ける。平面には穴のない正則な敷き詰めは 三種類しかない。正三角形、正方形、正六角形で平面を穴なく完全に覆う(敷き詰める)ことができる。これらの格子で辺が 分子結合を、頂点が分子を表すとすれば、固体(三角)→固体(正方)という相転移が考えられる。その際、三角敷き詰めの 分子は隣接分子への結合を六本もつのに対し、正方格子の分子は四本しかもたない。

融解(固/液)。 温度が上がると網が十分に多くの箇所で裂け、剛体ブロック全体がその剛性を失う(協同的な 亀裂 が マクスウェル剛性しきい値 z̄ = 2d を下回る — 網はフロッピーになる)。熱力学的ミセル、すなわちより大きな隣接領域と もはや結合で結ばれていない領域が形成される。一枚岩は分子的視点では瓦礫の山となり、条件付きで流動できる。最初の 空隙が、スライドパズルのように分子ブロックの移動を許す。隣接への結合の欠如がミセル内部の回転を許す。ミセル回転 (三和音-c)が始動する。多数の結合(結合)の破断が潜熱を呑み込む。

沸騰(液/気)。 さらに加熱すると空隙そのものがパーコレートする。空隙は連続した空虚へと結び付き、その中で ブロックは解け、個々の分子が飛び去る。沸騰は「物質の中の空隙」から「空隙の中の物質」への転移である。

昇華(固/気)。 低圧では液体の中間領域が欠ける。網はミセルのようなブロック形成を経ずに、個々の分子を直接 空虚へ放出する。ここでは、幾何学的側面(長い分子、絡み合う分子)、静電的側面(磁気モーメントや双極子モーメントに よる配向)、種々の分子結合型(ファン・デル・ワールス、メソメリー系間の相互作用、水素結合)が熱力学的ミセルの生成 あるいは非生成にどう影響しうるかを考えることができる。

三重点。 まさにそこで二つのしきい値が一致する。マクスウェル剛性(固体 ↔ 液体、剛性の喪失)と 空隙パーコレーション(液体 ↔ 気体)が*同じ* (τ, p) で生じる。三つの像(剛体ブロック[一枚岩]、ずれるブロック [瓦礫の山]、自由な分子[塵の雲])が共存する。

臨界点。 圧力が上がるにつれて液体→気体の跳びはますます小さくなり、空隙割合の変化はますます連続的になる。臨界点 では差異が完全に消える —「空虚の少ない密なブロック」と「分子の少ない多くの空虚」とが連続的に移り合う。 (「液体の視点ではグラスは半分空、気体の視点では半分満ちている。」)

二次元像における再結晶。 漫画は再結晶をいまや二次元の場合として上下に示す。下は密な三角格子 (配位 z = 6、二次元で最もエネルギーの低い理想)、上は緩い正方格子(z = 4)である。三角 → 正方の転移は 体積を増やし(圧力に抗する仕事)、結合数を下げる。配位数が替わっても結晶格子の剛性は保たれることが見て取れる。

エネルギー転移と熱容量の定性的予測。 連結性 z だけから両指標を見積もれる。(i) 粒子あたりの結合エネルギーは U/N ≈ −(ε/2)·z(各結合は半分ずつ数える)であり、z が高いほどエネルギーは低い。(ii) 格子交替時の c_v の跳びは モデル則 Δc_v = z_max/4 の変化に従う。二次元の三角(z = 6)↔ 正方(z = 4)の場合、これは粒子あたり ΔU ≈ εの エネルギー差と Δc_v = (6 − 4)/4 = 0.5 k_B の c_v 段差を与える — 三次元値(1.0 k_B)の半分である。標語: 密なほどエネルギーは低いが、転移は高くつく(潜熱が大きく、c_v 段差も大きい)。広がる際により多くの結合と より多くの凍結自由度が解放されるからである。

格子構造が準安定性を決めるのか。 部分的には。格子構造はエネルギーの順序を定める — どの形が基底状態(三角、 最低の μ)で、どれが励起状態(正方)かである。だが緩い形が低温で準安定に生き延びるかどうかを決めるのは、格子 構造だけではない。厳密に理想的なモデル(障壁なし)では、物質はどの (τ, p) でも直ちに最低の μをもつ相をとる — そこに準安定性はない。それは有限の活性化エネルギー(核形成障壁)によって初めて生じ、これが組み替えを速度論的に 妨げる。したがって格子幾何は熱力学的な序列を、障壁は速度論的な持続性を与える — 過飽和相を説明するのと同じ 成分である(下記参照。定量的には第6.7章)。

過飽和相 — 活性化エネルギー。 補足図(3.2、第2列)は液体の上の過飽和気相を示す。単独分子のほかに、2–3 個の 寄り集まった粒子から成る小集団(灰色の棒)が全体として飛んでいる。モデルはこの過飽和を分子結合に対する明確な 活性化エネルギーによって説明する。小さなクラスター(2–3)は縁と面積の比が不利で、その縁の結合は不完全である — 臨界核になって初めて自発的に成長を続ける。超臨界核が生じないかぎり、μ_気体 > μ_液体 がすでに凝縮を許していても 気体は過飽和のままである。まさにこの障壁が理想平衡像には欠けており、ここでは「理想像の外に位置する」独立の寄与 として名指しされる(定量的見積りは第6.7章)。

<a id="kap3-3"></a>

3.3 なぜ液体中の音速は T とともに下がるのか

音は弾性波である:v = √(K_eff/ρ)。K_eff はもっぱら健全な結合(結合)から生じる — それがばねである。T が上がると 協同的にますます多くの結合が破断し(亀裂)、ばね網は薄くなって K_eff は密度 ρ よりも急に落ちる。破断が支配する のは分母ではなく分子である — だから v は下がる。分子的解釈: 分子結合を介する音の伝達は、力学的衝突による機械的な 音の伝達より速い。温度上昇にともなって液体中の亀裂が増えるほど、音は分子結合ではなく力学的衝突を介して運ばれる回数が 増える。それに応じて音速は温度上昇とともに遅くなる。(水の*異常* — 約 74 °C における音速の*極大* — は第6.4章で扱う。)

<a id="kap3-4"></a>

3.4 ブラウン運動はいかにして生じるか — 鋭くしたミセル概念とともに

統計熱力学的な意味で、ここでのミセルは界面活性剤の集合体ではなく、亀裂 に縁取られたブロックである。欠けた 結合(ブロック境界における非結合)が、閉じ込められた小結晶(健全に結合した核)のまわりに「殻」を作る。この核は 内部の協同的な回転運動を行いうる(三和音-c)— 境界殻が隣接から切り離している間、ブロックは全体として回る。

ブラウン運動はこうして生じる。分子結合の交換(亀裂、敏捷)によって、個々の分子はあるミセルの殻を離れ、別の ミセルへ移りうる。そうした移行のたびに運動量が伝わり(保存。自由飛行 のおかげで損失なし)、関与する二つのミセルの 形と角運動量が跳躍的に変わる。大きな試験粒子(塵粒)は、そのような相関のないミセル再編の多数に洗われる。運動量 の突きの総和は平均でゼロだが有限の分散をもつ → ランダムウォーク。T が高いほど 亀裂 事象とミセル交替は頻繁になり、 震えは強くなる。ブラウン運動はこうして微視的なミセル交換過程の巨視的な影あるいは帰結である。

<a id="kap3-5"></a>

3.5 なぜミセル境界層で反応が速くなるのか

きまり悪い前置き。 以下の議論は像からの導出であり、私の知るかぎりこれまで実証的に裏づけられても研究されても いない。もっとも、それについて詳しく調べる労を取らなかったことは白状しなければならない。新規性の手がかりは、AI システムがこの帰結に触れないという一点にすぎない。この命題は、絵画的な想像がいかに新しい問いへ導きうるかの一例である。

中程度にしか溶かさない溶媒では、反応物は好んでミセルの境界殻に付く — モデルではまさに亀裂/ブロック境界(亀裂)で ある。二つの効果が反応を速める。(i) 次元の低下 — 反応物はもはや全体積の中で互いを探すのではなく、擬二次元の界面に 沿って探す。二次元の探索は三次元の探索より統計的に速く当たる。(ii) 濃縮 — 境界殻が反応物を集め、その局所密度が 上がる。絵画的には、亀裂境界は相手どうしを導き寄せるガイドレールのように働く。

次元効果には厳密な核がある:ランダムウォーク(ポリア)。 二次元の探索が三次元の探索より統計的によく当たる理由は、 直観だけでなく確率論の定理から従う。ジョージ・ポリア(1887–1985)は 1921 年に格子上のランダムウォークについて次を示した: 一次元と二次元では歩行者は確率 1 で(すなわち確実に)任意の点へ無限回戻る。平面格子 ℤ² でも直線 ℤ でも同様である。 三次元以上ではこれはもはや成り立たない。回帰確率は真に 1 より小さい(単純立方格子 ℤ³ ではおよそ 0.34 にすぎない)。 直観的には、ポリアの言い回しで:「酔った人はいつか確実に家へ帰り着くが、酔った魚はそうならない。」成功する衝突のために 何度も出会わねばならない二つの反応物に当てはめると、こうなる:自由な三次元ランダムウォークをもつ非常によく溶かす 溶媒では、いったん離れた相手は有限の確率で永遠に離れていく — 反復する反応の試みは失敗し、収率はこの純粋に幾何学的・ 統計的な理由から低いままになる。これに対し(擬)二次元のミセル境界層では再会は確実である(確率 1)。相手どうしは 反応が成功するまで何度でも巡り合う。したがって境界層は濃縮するだけでなく、同時に再会に有利な次元を用意する。 (意識的に尖らせた)主張はこうである:非常によく溶かす溶媒での収率は、ランダムウォーク的理由から総じて悪くなりがちである。

典拠: G. Pólya (1921), Über eine Aufgabe der Wahrscheinlichkeitsrechnung betreffend die Irrfahrt im Straßennetz, Math. Ann. 84, 149–160 — ここでは M. Eisermann, Höhere Mathematik 3(シュトゥットガルト大学)、問題 W-2×2 に よる引用。(原典閲覧 2026-09-05)

<a id="kap3-6"></a>

3.6 三次元の液体窓

3.1–3.5 の図とは異なり、この図は分子の場面を示すのではなく、第7.9/7.10章の計算結果を τ 軸上の模式図として 示す(結合エネルギー ε の単位)。 目的は、単一の二次元格子では液相のための領域がほとんど許されないのに対し、三次元格子では問題なく可能であることを 見えるようにすることであった。しかも両者には同じ剛性則(z̄_c = 2d)が当てはまる。

図① — 二次元(窓なし)。 τ 軸上で τ_m = 1.20 と τ_triple = 1.16 は実質的に重なり(τ_triple は τ_m の下)、 τ_c = 1.50 がそのすぐ上にある。液相のための領域は隙間に押し潰されている — 固体はほぼ直接に気体へ移る(昇華に近い)。

図② — 三次元(明瞭な窓)。 同じ軸だが引き伸ばされている:τ_m = 1.20 < τ_triple = 2.83 < τ_c = 3.00。τ_triple と τ_c の間に独立した液体の領域がある(下方には広く、τ_triple − τ_m = 1.63。上方には狭く、τ_c − τ_triple = 0.18)。 狭い上端はそれとして印されている。

図③ — その原因。 z̄_c = 2d は次元とともに増え(4 → 6)、それにより v_solid は下がり(換算体積で 0.50 → 0.667、 すなわち固体は液体に対して「より緩く」なる)、τ_c = z_max/4 は倍になる(1.50 → 3.00)。臨界点は融点のはるか上へ 移る。これによって初めて液相の場所が生まれる。これは純粋な次元効果である。

<a id="kap4"></a>

4 — 統計力学的モデル(版09から変更なし)

要約

第4章はモデルの統計力学的骨組みを二次元で立てる(三角格子、z_max = 6):

• 4.1 自由度と変数: 秩序量は結合割合 x ∈ [0,1] であり、平均配位数は z̄ = z_max·x。ブロック(≥ 2)ごとに回転 自由度が一つ加わる。温度 τ は結合エネルギー ε の単位で数える。
• 4.2 協同的二状態結合: 自己無撞着な結合方程式 x = 1/(1 + G·exp[−(ε + K·x)/τ])(K = J·z_max、G = z_max)は ブラッグ–ウィリアムズ自由エネルギーの停留条件として従う。その大域的最小が物理的な枝を与え、その枝が連続か 跳びを伴うか(一次)は第5.5章が決める。
• 4.3 協同項を含むエネルギー収支: U/N = τ·[1 + z̄/4 + ½ n_blk] − (ε/2)·z̄ − (K/2)·x²(二次元)は並進・ 結合振動・ブロック回転、結合の配置エネルギー、そして新しい協同的分子場エネルギー −(K/2)x² を統合する。対応する 熱容量 c_v は正の追加項 −K·x·ẋ を得る。これは融解ピークを鋭くするが、固体プラトーと気体枝には手を触れない。

<a id="kap4-1"></a>

4.1 自由度と変数

• 固定された参照配位幾何の上に N 個の分子。正準的な二次元の場合:三角格子、z_max = 6、ブロックあたり回転自由度 一つ(多数性、結合)。本章は明示的に二次元解釈に当てはまる。三次元への一般化は第6章にある。
• 幾何学的隣接対ごとの結合変数 s_ij ∈ {0, 1}(1 = 健全)。
• 配位数 z_i = Σ_j s_ij、その平均 z̄ = 2 N_bond / N が秩序パラメータである。結合割合 x = p_b ∈ [0, 1] を 用いて z̄ = z_max · x。
• エネルギー単位 τ ≡ k_B T(古典的、二次形式の自由度あたり ½τ)。すべてのエネルギーは結合エネルギー ε の単位で表す。

<a id="kap4-2"></a>

4.2 協同的二状態結合(イジング型、パラメータ J)

孤立した結合(健全で −ε、破断で 0)はマクスウェル–ボルツマン平衡で p_b⁰(τ) = 1/(1 + e^(−ε/τ)) をもつ。協同的な 亀裂 (公理 6)は、実効的な結合安定性を結合割合 x に比例する分子場の分だけずらす。分子場(イジング)近似では自己無撞着な 結合方程式は次となる:

x = 1 / (1 + G · exp[−(ε + K·x) / τ]), K = J·z_max, z̄ = z_max·x

ここで:

• K = J·z_max は協同的結合強度である(亀裂 のフィードバック:隣接結合が多いほど → 実効的な結合安定性が高い)。 J > 0 は転移を鋭くする。臨界の J を超えると z̄(τ) は不連続になる — それが融解の潜熱である。純粋な二状態の 場合は極限 J = 0(分子場なし、連続的クロスオーバー)である。
• G > 1 は破断状態の縮退である。自由になった分子は位相空間(並進の自由)を得る。G は高い τ で結合の破断を駆動し、 x が ½ に留まらず x → 1/(1+G) となるようにする(気体極限)。正準的には G = z_max。

この方程式は分子あたりのブラッグ–ウィリアムズ自由エネルギーの停留条件として従う(n_slots = z_max/2 の結合席):

g(x)/n_slots = −ε·x − (K/2)·x² + τ·[x ln x + (1−x) ln(1−x)] − τ·ln(G)·(1−x)

項 −(K/2)x² は協同的(強磁性的)結合であり、最後の項は破断の縮退である。x にわたる g の大域的最小が物理的な枝を 与える。二つの最小が同じ深さになるところに一次転移がある。敏捷(公理 7)は x(τ) がどの温度でも平衡値をとることを 正当化する。

<a id="kap4-3"></a>

4.3 エネルギー収支 — 協同項を含む

分子場項 −(K/2)x² は真のエネルギー(τ に依らない)であり、U/N に現れねばならない。整合的な収支はこうなる(二次元):

U/N = τ·[1 + z̄/4 + ½ n_blk] − (ε/2)·z̄ − (K/2)·x²

• τ·[…]:並進(1)+ 結合振動(z̄/4)+ ブロック回転(½ n_blk)— 三和音。
• −(ε/2)·z̄:結合の配置エネルギー — 結合。
• −(K/2)·x²:協同的分子場エネルギー — 亀裂。 U/N において新しい。

対応する熱容量 c_v = d(U/N)/dτ(k_B 単位)は追加の寄与を得る:

c_v = [1 + z̄/4 + ½ n_blk] + τ·[¼ ż̄ + ½ ṅ_blk] − (ε/2)·ż̄ − K·x·ẋ, (ẋ = dx/dτ)

ここで新項 −K·x·ẋ は正であり(ẋ < 0)、多数の結合が同時に破断するまさにその場所で効く。これは融解ピークを 鋭くするが、固体プラトー(ẋ ≈ 0)も気体枝も変えない。

<a id="kap5"></a>

5 — 分配関数と熱容量

要約

第5章は分配関数と熱容量を計算し、較正(第7章)のためのモデルデータを提供する:

• 5.1 因子化された分配関数: 時間尺度の分離により Z = Z_trans·Z_vib·Z_rot·Z_bond と因子化する。
• 5.2 相ごとの熱容量: c_v = d(U/N)/dτ。正の項 −(ε/2)(dz̄/dτ) が潜熱の類似物として ショットキー型のピークを生む。
• 5.3 数値的 c_v(τ) 曲線: τ_m = 1.20 で c_v ピーク 5.51、固体プラトー 2.50、気体枝 → 1。 潜熱 L ≈ 1.41 ε/分子(協同項あり)は 0.48 ε(なし)の三倍である。
• 5.4 / 5.4a 二つの分離したしきい値: 連結性パーコレーション z̄_c = z_max·p_c ≈ 2.08 (τ_b ≈ 1.85 において)はクラスター崩壊(沸騰近傍)を印し、結果である。これは 一般的マクスウェル剛性 z̄_c = 2d(融解密度、剛性の喪失)と明確に区別されねばならない。
• 5.5 融解転移の次数: 判定条件 ln G > 2 + 4ε/K。2D(G = 6)は*どの* J でも連続、3D(G = 12)は K > ~8.2 から一次(実演:Δz̄ ≈ 5.3, ΔU ≈ 1.75 ε)。
• 5.6 有限サイズ研究: p*(L) は下から単調に p_c = 0.3473 へ収束する。z̄_c ≈ 2.08 は有限サイズに 対して頑健である。
• 5.7 完全な三次元 c_v 走査: K > K_臨界 ≈ 8.25(J > 0.68)から真の跳び Δz̄ が潜熱 ΔU とともに開く。 融解ピークは J とともにより高い τ へ移動する(連続な二次元ピークに対する正の移動 +0.5 … +2.2)。

<a id="kap5-1"></a>

5.1 因子化された分配関数

時間尺度の分離(敏捷)により、正準分配関数は因子化する:

Z = Z_trans · Z_vib · Z_rot · Z_bond.

• Z_trans = (1/N!)·(A/λ_T²)^N, λ_T = h/√(2π m τ) (古典理想の二次元並進。多数性、保存、自由飛行)。
• Z_vib = (τ/ħω₀)^(N_bond) (健全な結合ごとの古典調和振動子。三和音-b)。
• Z_rot = ∏_(ブロック≥2) (2π I_b τ)^(1/2)/ħ (ブロックごとの古典二次元剛体回転子。三和音-c)。
• Z_bond = Σ_{s} exp[(ε/τ)·Σ s_ij + (K/τ)·協同項] (イジング型の結合和。結合、亀裂)。分子場評価が 4.2 の x(τ), z̄(τ) を与える。

<a id="kap5-2"></a>

5.2 相ごとの熱容量(解析的)

U/N(第4.3章)から、c_v = d(U/N)/dτ(k_B 単位):

c_v = [1 + z̄/4 + ½ n_blk] + τ·[¼ (dz̄/dτ) + ½ (dn_blk/dτ)] − (ε/2)·(dz̄/dτ).

項 −(ε/2)·(dz̄/dτ) は正であり(dz̄/dτ < 0:加熱は結合を破断する)、多数の結合が同時に破断するところで 支配的になる — それがショットキー型のピーク、潜熱の類似物である。

<a id="kap5-3"></a>

5.3 協同項を含む数値的 c_v(τ) 曲線 — 潜熱の定量化

自己無撞着な結合(大域的最小)、n_blk(τ) はボンドパーコレーションのモンテカルロ(Union-Find, L = 100, 5 回の実現)から、c_v = dU/dτ(中心差分)、U/N は −(K/2)x² を含む(cv_scan_09.py, ε = 1, J = 0.55, G = z_max = 6)。選んだ支点:

結果(数値的)。 三相の順序は引き続き確認され、いまやより鋭い:

c_v(固) ≈ 2.50 < c_v(液、最大) ≈ 5.51 > c_v(気体枝) ↓ → 1.

潜熱(ベースラインより上のピーク面積)。 (c_v − ベースライン)を融解窓 τ ∈ [0.70 ; 1.70] にわたって 積分すると(ベースライン = c_v(0.70) と c_v(1.70) を通る直線)、

L ≈ 1.41 ε/分子 (協同項あり)、 L ≈ 0.48 ε/分子 (協同項なし、旧収支)、 ΔL ≈ +0.92 ε/分子。

協同項はピークを約 1.7 倍に持ち上げ(3.25 → 5.51)、潜熱を三倍にする。物理的には、多数の隣接結合が 同時に破断すること(亀裂)は配置エネルギー(−ε/2·z̄)だけでなく分子場エネルギー(−K/2·x²)も蓄えて おり、融解時に両者が解放される。

<a id="kap5-4"></a>

5.4 結果としての連結性パーコレーションしきい値

巨大クラスター割合 P_∞(τ) は τ ≈ 1.85 で崩壊する。そこでは x ≈ 0.347 = p_c(三角) = 2·sin(π/18) = 0.3473、すなわち z̄_c = z_max·p_c ≈ 2.08 — 設定ではなく*結果*である。これは 連結性のしきい値(貫通クラスター)であり、沸騰(τ_b)近傍のクラスター崩壊を印す。マクスウェル剛性 (z̄_c = 2d、融解)と混同してはならない — 第5.4a章の明確な分離を参照。

<a id="kap5-4a"></a>

5.4.a 二つの異なるしきい値を明確に分離

これまで曖昧だった「パーコレーション/剛性」という二重呼称をここで解消する。混同してはならない 物理的に異なる二つのしきい値である:

理由。 融解は密なブロックの機械的剛性の喪失であって、連結の崩壊ではない。したがって関連する量は マクスウェルの拘束数え上げ(平均場/一般的剛性)である。中心力から成る D 次元の網は、平均配位が z̄ = 2D に達したときちょうど等静的になる(2D → 4, 3D → 6)。ボンド連結性(z̄_c = 2.08)はそれより はるかに下にある。z̄ = 2.08 でかろうじて連結している網はとうにフロッピー(軟らかいモード)であり、 すでに液体である。だから連結性は沸騰近傍のより遅いクラスター崩壊を印し、剛性はより早い 固化/液化(融解密度)を定める。この分離は第7.8章の三重点をパラメータ不要かつ現実的に固定する。 5.4/5.6 で確認された連結性しきい値 2.08 を変えるのではなく、その解釈だけを変える(τ_b = 沸騰近傍 であって融解ではない)。

<a id="kap5-5"></a>

5.5 融解転移の次数

解析的判定条件。 スピノーダル領域(g の二つの最小)は τ < K/4 で存在し、転移はおよそ τ_t ≈ (ε + K/2)/ln G にある。真の跳び(一次転移)は τ_t < K/4 を要求し、すなわち K → ∞ の極限で:

ln G > 2 + 4ε/K.

ここから破断縮退 G = z_max についての注目すべき結果が従う:

数値的確認(τ 刻みごとの max |Δx|、detect_first_order):

解釈。 正準的な二次元モデル(G = 6)では融解は原理的に連続である(「軟らかい」転移/クロスオーバー) — ln 6 < 2 であるから、どの J でも真の跳びは生じない。三次元充填のより高い配位/縮退(G = 12)になって 初めて、十分な協同性(K > 8.2)から有限の潜熱を伴う一次転移が許される。これは現実と符合する。二次元 融解はしばしば連続的(KTHNY 型)であり、三次元融解は一次である。J = 0(分子場なし)は常に連続な参照 の場合である。— 物理的には J は協同エネルギーと結合エネルギーの比(K/ε = J·z_max/ε)を意味する: 健全な隣接結合が自身の結合をどれだけ安定化するか、である。

一次の実演(三次元に動機づけられた、cv_scan_09.py, J = 1.0, G = 12, z_max = 12): z̄ は τ ≈ 2.82 で Δz̄ ≈ 5.30 だけ跳ぶ。潜熱はエネルギーの跳びとして ΔU ≈ 1.75 ε/分子。同じ K を G = 6 の格子 (2D)で用いても連続のままである — 次数の性格はひとえに G = z_max にかかっている。

<a id="kap5-6"></a>

5.6 パーコレーションしきい値の有限サイズ研究

L × L 三角格子(周期的)上のボンドパーコレーション、巨大クラスター割合 P_∞(p) を L にわたって スケールした(cv_scan_09.py; L = 50/100/200/400、それぞれ 8/5/3/2 回の実現)。厳密なしきい値 p_c = 2·sin(π/18) = 0.34730 → z̄_c = z_max·p_c = 2.0838。

P_∞(p)(平均 ± ばらつき):

しきい値推定 p*(L)(P_∞ が 0.5 を横切る):

結果。 p*(L) は単調に下から p_c = 0.3473 へ収束する。p_c におけるばらつきは ±0.103(L = 50)から ±0.017(L = 400)へ下がる — 転移は L とともに鋭くなるが、しきい値を移動させない(普遍的パーコレー ション値)。これにより z̄_c ≈ 2.08 = z_max·p_c は有限サイズに頑健なモデル値として確認された。

<a id="kap5-7"></a>

5.7 一次融解の完全な三次元 c_v(τ) 走査

第5.5章の単一の一次実演走行(J = 1.0, G = 12)を、いまや協同性 J にわたる完全な c_v(τ) 走査へと 拡張する(firstorder_3d_16.py、n_blk には FCC 連結性、三次元エネルギー収支 U/N = τ·[3/2 + z̄/4 + (3/2)n_blk] − (ε/2)z̄ − (K/2)x²)。第5.5章の解析的次数判定条件(一次 ⇔ K > K_臨界 = 4ε/(ln 12 − 2) = 8.25)が場合を分ける。連続な二次元参照ピークは τ_2D = 1.20 (c_v ≈ 5.5)にある。

*連続の場合:Δz̄/ΔU は急峻なクロスオーバー(真の跳びではない)。τ_転移 はそこでは c_v ピークである。

結果。 (i) 連続な三次元の場合 J = 0.55 は τ_転移 = τ_m ≈ 1.72 を与える — 独立な結合速度論 (bondkinetics_3d_15.py: τ_m(FCC) = 1.70)との照合は一致する。(ii) K > K_臨界 ≈ 8.25 (J > 0.68) から真の跳び Δz̄ が有限の潜熱 ΔU(不連続をまたぐエネルギーの跳び)とともに開く — 融解は一次になる。まさに ln 12 > 2 という判定条件が要求するとおりである。(iii) 照合 J = 1.0 は 第5.5章の実演値 Δz̄ ≈ 5.4 を再現する。(iv) 跳びの高さ Δz̄ と位置 τ_転移 は J とともに増す:転移は 上方へ移動する — 三次元充填のより高い配位*と*より強い協同性が融解温度を押し上げる。連続な二次元 ピーク(τ_2D = 1.20)に対して、こうして一貫して正のピーク移動(+0.5 … +2.2)が生じ、それは K とともに 単調に増加する。三次元融解の一次的性格は、これでもはや*一つの*点だけでなく、協同性の全域にわたって 裏づけられた。

<a id="kap6"></a>

6 — 三次元解釈の数学的もっともらしさ

要約

第6章は結合像の 2D から 3D への一般化を純粋に数学的にもっともらしくする(独自の HTML 可視化なし):

• 6.1 回転自由度 2D → 3D: ブロック回転は一軸から三軸へ、並進は二方向から三方向へ、z_max は 6 から 12 へと増える。c_v の基礎値は上がり(c_v,固 2.5 → 4.5, c_v,気 1 → 3/2)、c_v(τ) 曲線の形は保たれる。
• 6.2 識別特徴としての連結性: z_max は量子化された量であり、最密充填(FCC/HCP, z_max = 12)を BCC(z_max = 8)から中間状態なしに跳躍的に分ける。
• 6.3 追加の転移 固体(1)→固体(2): 離散性から FCC↔BCC の多形がパラメータ不要の c_v 段差 1.0 k_B とともに従う。厳密な理想像(θ_BCC = θ_FCC)では生じず、より軟らかい BCC(開口しきい値 0.61 を下回る θ_BCC < θ_FCC)になって初めて τ_ss < τ_m が開く。
• 6.4 / 6.5 考察: 液体の水の音響異常を二つの格子型の不明瞭な近距離秩序交替として、また、ミセル回転 に対するフラクタル渦場の代替案(純粋に定性的)。
• 6.6 / 6.8 / 6.9 三段階の微視的 θ(格子): アインシュタインばねモデル(θ_BCC/θ_FCC ≈ 0.82…0.94)、 デバイ型フォノン状態密度(≈ 0.87…1.00)、レナード–ジョーンズ完全格子(≈ 1.0;理想 BCC はせん断不安定、 FCC は安定な基底状態)— すべてしきい値 0.61 の上にあり、したがって現実の多形は調和的理想像の外に ある。
• 6.7 準安定性(2D): 三角 ↔ 正方の組み替えについての桁数計算(ΔU = +1 ε, ΔA ≈ +15.5 %)と、速度論的 追加としての核形成障壁 ΔG\* = πλ²/Δμ_A。

本章は三次元への一般化を純粋に数学的にもっともらしくする。三次元は HTML 版では決して可視化され ない。 目的はひとえに、回転自由度と連結性が 2D から 3D へどう変わるか、そしてそれによってどんな新しい 転移が考えられるようになるかを示すことにある。

<a id="kap6-1"></a>

6.1 回転自由度の 2D → 3D 変化

2D に比べて質的に新しい唯一の点は、ブロックの回転である:

• 2D: 剛体ブロックは回転自由度を 1 もつ(面内回転)→ U_rot/N = ½ n_blk τ。
• 3D: 剛体ブロックは回転自由度を 3 もつ(三つのオイラー軸)→ U_rot/N = (3/2) n_blk τ。

同時に並進は 2 から 3 自由度へ増え(U_trans/N = (3/2)τ)、最大配位数は z_max = 6(三角)から z_max = 12 (最密球充填)へ増える。三次元のエネルギー収支はこうなる:

U/N = τ·[3/2 + z̄/4 + (3/2) n_blk] − (ε/2)·z̄,

熱容量の基礎値も整合的に移動する:

曲線の構造(固体プラトー、−(ε/2)dz̄/dτ による液体最大、気体値への低下)は同一のままであり、次元と z_max とともにスケールするのは水準だけである。

<a id="kap6-2"></a>

6.2 各種結晶格子の影響 — 識別特徴としての連結性

3D では少なくとも二つの理想格子を区別せねばならず、それらは連結性(配位数 z_max)で分かれる:

連結性 z_max こそ決定的な識別特徴である: それは固体の熱容量水準(z_max/4)もパーコレーションしきい値 z̄_c も定める。どの格子も*一つの*固定した配位数しか許さないから、連結性は量子化された有限の量である — z_max = 12 と z_max = 8 の間に連続的な移行はなく、ある格子型から別の格子型への離散的な跳びがあるだけ である。

<a id="kap6-3"></a>

6.3 追加の転移 固体(1) → 固体(2)

連結性の有限性/離散性から、二次元モデルには存在しない転移が従う:固体(1) → 固体(2)(多形/同素体)。 異なる z_max をもつ二つの固相(たとえば FCC ↔ BCC)は異なる自由エンタルピー μ(τ, p) をもつ。 μ_{固(1)} = μ_{固(2)} となる温度・圧力の組み合わせで、物質は融けることなく格子を替える。モデル像では その際配位幾何が跳躍的に変わる(融解のように占有 z̄ だけが変わるのではない)— 通常は固定される参照 幾何そのものが跳ぶ唯一の状況である。観測されるのは小さな c_v 段差(z_max/4 の変化)と潜在変態熱で あり、融解ピークよりはるかに弱い。

一次判定条件 ln G > 2 + 4ε/K(第5.5章)は、なぜまさに三次元解釈が定性的により豊かなのかについて 具体的な理由を与える。そこでのみ(G = z_max = 12, ln 12 = 2.485 > 2)融解が一次になりうる — そして連結性が量子化/有限であるがゆえに、そこでのみ離散的な格子交替 固体(1) → 固体(2)(FCC ↔ BCC) が存在する。両現象は同じ原因の二つの装いである:三次元充填のより高く離散的な配位が、鋭い融解の跳びも 多形も生む。二次元単一格子(G = 6)はどちらももたない — 第5.5章と整合的である。

定量化(bondkinetics_3d_15.py)。 この転移をいま簡潔に計算する。両格子はそれぞれ古典調和像での 自由エンタルピーを受け取る(共通の p·v と並進の寄与は比較で相殺する):

μ(τ) = h₀ − c_v·τ·ln(τ/θ), h₀ = −(ε/2)·z_max, θ = ħω/k_B(ε 単位のフォノン剛性)、

第6.1/6.2章からのパラメータ不要の特性値とともに:

ここから c_v 段差がパラメータ不要で:|Δc_v| = z_max/4 の変化 = (12−8)/4 = 1.0 k_B(固体は FCC→BCC の転移で熱容量が 1 k_B「軽く」なる)— 融解ピーク(数 k_B、第5.3/7.11章)よりはるかに小さい。

変態温度 τ_ss(μ_FCC = μ_BCC)と潜在変態熱 L = τ_ss·(S_BCC − S_FCC):

結果(正直な所見)。 厳密な理想像(すべての結合に同じフォノン剛性、θ_BCC = θ_FCC)では固体(1)→ 固体(2)転移は存在しない。FCC はエンタルピー(より深い、−6 対 −4 ε)でも、結合数がより多いために (c_v = 4.5 > 3.5)振動エントロピーでも勝つ。μ_FCC はどの温度でも μ_BCC の下にとどまる。格子依存の フォノン剛性 θ_BCC < θ_FCC(より開いた BCC は軟らかく振動し、したがってエントロピーが大きい)になって 初めて、高温で μ_BCC が μ_FCC の下に来る — そのとき変態が開き、それは融解ピーク τ_m の下に位置し、 固体(1)→固体(2)→液体→気体という順序が生じる(τ_ss < τ_m < τ_c ≤ τ_b(3D)、第7.11章)。軟らかい BCC 剛性は厳密な理想像の外への穏当な拡張である(現実には鉄の γ/δ 多形として知られる:高温で BCC、その 間に FCC)。これに対し c_v 段差 1.0 k_B と小さな正の変態熱はモデル内在的かつパラメータ不要である。

θ_BCC/θ_FCC は実際どれほど小さいか?(theta_springs_16.py) 変態を自由な θ パラメータから解放する ため、フォノン剛性をいま微視的に見積もる(第6.6章)。中心対ばねによるアインシュタインばねモデルは θ_BCC/θ_FCC ≈ 0.82(最近接のみ)から ≈ 0.94(BCC 第二殻を含む)を与える。これに対し μ 比較から 必要な開口しきい値は θ_BCC/θ_FCC < 0.61 にある(それを下回れば τ_ss が存在し τ_m の下にくる)。 すなわち幾何学的な剛性比は大きすぎる:厳密な理想像では FCC はどの温度でも安定であり、純粋に幾何 学的に駆動される固体(1)→固体(2)転移は存在しない。現実の多形(鉄 γ↔δ)は理想的な調和像の外にある 追加の(電子的/磁気的/非調和的)寄与を必要とする — まさに上で置いた「穏当な拡張」である。これは正直な 所見を裏づけ、θ 比を自由な仮定から見積もられ、名指しされたモデルの限界へと変える。

<a id="kap6-4"></a>

6.4 考察:液体の水の音響異常

普通の液体では v(T) は単調に下がる(第3.3章)。水は ≈ 74 °C に音速の*極大*をもつ:v はまず T ととも に上がり、その後下がる。モデルではこれを二つの局所格子型の競合(第6.2章)として解釈できる: 開いた、低配位の氷に似た近距離秩序(小さい z_max、多くの「穴」)と、密で、より高配位の近距離秩序 である。低い T では開いた構造が支配する。加熱はまず*開いた*近距離秩序を、より密なものに有利なように 壊す — 実効的な連結性、したがって K_eff が上がる(v が上がる)。極大を超えて初めて通常の亀裂による 結合破断が優勢になり、K_eff が下がり、v が減る。この異常はこうして液相内部での、不明瞭で不完全な 固体(1)↔固体(2)的な近距離秩序交替の表れである — 二つの格子型をもつ三次元モデルでのみ考えられ、二次元 単一格子では考えられない。

<a id="kap6-5"></a>

6.5 代替案:ミセル回転に代わるフラクタル渦場

局所的なブラウン運動(第3.4章)はミセル回転と交換によって説明された。代替として、それをフラクタル 渦場の効果として思い描くこともできる:亀裂によって絶えず組み替えられるブロック境界が局所的な運動量流の 場を生み、その場は — 乱流に似て — あらゆる長さスケールで自己相似な渦をもつ。すると試験粒子は離散的な ミセルの突きではなく、スケール不変な渦の連続体によって駆動される。両方の像(離散的ミセル ↔ フラク タル渦)は同じ亀裂速度論の相補的な極限であり、どちらが効くかは試験粒子の大きさと平均ミセル寸法の比に 依存する。この試論は純粋に定性的であり、ここでは数学的に展開しない。

<a id="kap6-6"></a>

6.6 微視的フォノン剛性 θ(格子):アインシュタインばねモデル

第6.3章の固体(1)→固体(2)比較は相対的なフォノン剛性 θ_BCC/θ_FCC にかかっている。それを自由に選ぶ 代わりに、ここでは配位と幾何から見積もる(theta_springs_16.py)。

モデル。 中心原子はその z 個の隣接がつくる調和的な椀の中に座る。各結合は径方向の対ばね κ である。 変位 u は結合方向へ ⟨cos²⟩ = 1/d = 1/3(d = 3)で等方的に射影されるので、実効ばね定数は k_eff = (z/3)·κ となる。アインシュタイン振動数は ω = √(k_eff/m)、すなわち θ = ħω/k_B ∝ √k_eff。 同じ質量と同じ対ばねなら、

θ_BCC/θ_FCC = √(z_BCC/z_FCC) = √(8/12) = 0.816.

精密化。 (B) BCC の第二殻(r₂/r₁ = 2/√3 = 1.155 に 6 原子)を距離でスケールしたばね κ(r) ∝ (r₁/r₂)ⁿ(n = 6…8、レナード–ジョーンズ型)とともに考えると比は 0.91…0.94 に上がる(近い 第二殻によって BCC は実効的に硬くなる)。(C) 同じ密度における BCC のより短い結合 (r_nn(BCC)/r_nn(FCC) = 0.972, −2.8 %)はばねをわずかに硬くし、0.89…0.91 を与える。総じて 幾何学的な窓 θ_BCC/θ_FCC ≈ 0.82…0.94。

開口しきい値。 μ 比較(第6.3章)が変態を開く(τ_ss が存在し τ_m = 1.70 の下にくる)のは θ_BCC/θ_FCC < 0.61 になってからである。幾何学的な窓はそれよりはるかに上にある。

所見。 純粋なばねモデルは BCC を十分に軟らかくしない:厳密な理想像では FCC が安定なままであり、 幾何学的に駆動される固体(1)→固体(2)転移は生じない。これにより θ_BCC/θ_FCC は自由な仮定から 見積もられた量になり、その自然な値(≈ 0.8–0.9)は変態を閉じたままに保つ — 第6.3章の正直な所見と整合 的である。したがって現実の多形は、単なる格子幾何ではなく、理想的な調和像の外(電子的/磁気的寄与、 非調和性)にあることが立証される。

<a id="kap6-7"></a>

6.7 準安定性の数学的見積り(2D 三角 ↔ 正方)

第3.2章の定性的議論に対し、ここでは二次元像での対応する桁数計算を示す(硬い円板、直径 D、接触距離 a = D)。

(1) エネルギー。 粒子あたり U/N = −(ε/2)·z(各結合は半分ずつ):三角 z = 6 → −3 ε、正方 z = 4 → −2 ε。したがって三角 → 正方の組み替えは粒子あたり ΔU = +1 ε を要する(吸熱的。緩い形のほうが エネルギーが高い)。

(2) 体積(面積)と仕事。 粒子あたりの面積:三角 A_tri = (√3/2)·D² ≈ 0.866 D²、正方 A_quad = D²。 広がりは ΔA = A_quad − A_tri ≈ 0.134 D²(+15.5 %) — 漫画で示された体積増大である。(換算)圧力 p に対してそのために仕事 W = p·ΔA を要する。自由エンタルピーは Δμ = ΔU + p·ΔA − τ·ΔS だけ増える。 エントロピー項 ΔS > 0(緩い形は軟らかいモードをもち、c_v 段差も大きい Δc_v = (6−4)/4 = 0.5 k_B)は 高い τ になって初めて μ_quad を μ_tri の下へ下げる — だから三角が低温の基底状態で、正方が高温/準安定 の形なのである。

(3) 核形成障壁(2D)。 準安定性は障壁を前提とする。2D では新相の核(半径 r)が面積利得 −Δμ_A·π r²(Δμ_A = 面積あたりの駆動力)で、縁の寄与 +λ·2π r(λ = 長さあたりの縁エネルギー)に抗して 成長する: ΔG(r) = −Δμ_A·π r² + 2π λ r. 最大は r\* = λ / Δμ_A にあり、障壁は ΔG\* = π λ² / Δμ_A. 臨界核 r\ は駆動力 Δμ_A の増大とともに縮む。熱的スケール τ = k_BT が障壁に達しないかぎり (**ΔG\ ≳ 数 k_BT)、緩い(あるいは沸騰では過飽和の)相は速度論的に凍結したままである。小さな クラスター(2–3 粒子、過飽和のコマのように)では縁と面積の比が最も不利であり、それらは r\* の左に あって平均としては再び解ける。こうして漫画で示された活性化エネルギーは、準安定な結晶形も過飽和気相も 支える同じ量 ΔG\* である。それは理想的な(障壁のない)平衡モデルへの速度論的**な追加である。

<a id="kap6-8"></a>

6.8 θ(格子) のデバイ型精密化

アインシュタインばねモデル(第6.6章)はフォノン分布全体を一つの振動数で置き換える。ここでの手法は 格子ごとの簡潔なフォノン状態密度によって θ(格子) を精密化する(theta_debye_17.py):同じ中心対ばね κ(r) = r^(−n) を最近接および第二隣接殻に置いて動力学行列 D(k) = Σ_j κ_j·(R̂_j⊗R̂_j)·(1 − cos(k·R_j)) をつくり、これを原始ブリルアンゾーンの k 格子(モンクホルスト –パック、12³ 点)で対角化する(解析的な 3×3 固有値ソルバ、標準ライブラリのみ)。三つのフォノン枝から、 振動数モーメントとして二つのデバイ型剛性尺度が従う:

• 2 次モーメント θ(2) ∝ √⟨ω²⟩ — アインシュタイン極限(1 − cos の BZ 平均は 1 となるので ⟨ω²⟩ = z·κ/3)。最近接のみの対照、純粋な幾何:√⟨ω²⟩_BCC / √⟨ω²⟩_FCC = 0.816 = √(8/12) — まさに第6.6章のアインシュタイン値である。
• 対数モーメント θ(0) ∝ exp⟨ln ω⟩ — 本来のデバイ型の量であり、アインシュタイン値が無視する軟らか い音響枝に重みを与える。

スケールをもつ定数(ħ/k_B、κ のスケール、デバイの前因子)はすべて θ_BCC/θ_FCC という比で相殺し、 DOS の形だけが残る。BCC をそもそも機械的に安定化する第二隣接殻とばね則 n = 0…8 のもとで、

θ_BCC/θ_FCC ≈ 0.87…1.00 (モーメントの定義とばね則にわたる誤差帯)。

これは純粋なアインシュタイン値(0.82…0.94)よりやや高い。同じ密度でより短い BCC 結合がばねを 硬くするからである。μ 比較(第6.3章)から必要な開口しきい値 0.61 は依然としてはるかに下にある: フォノン DOS もまた BCC を十分軟らかくしない。第6.6章の所見はこうしてより細かい近似によって裏づけられ、 精密化される — 厳密な理想像では FCC はどの温度でも安定であり、純粋に幾何学的に駆動される固体(1)→ 固体(2)転移は生じない。現実の多形(鉄、第7.13章)は理想的な調和像の外の寄与を必要とする。

<a id="kap6-9"></a>

6.9 レナード–ジョーンズ完全格子から得る θ(格子)

第6.8章は nn + 2nn 上で打ち切った中心ばね κ(r) = r^(−n) を用いる。版17は同じ問いを、対ポテンシャル から得た本物の力定数で完全格子にわたって検証する(theta_ljlattice_18.py):各結合はレナード– ジョーンズポテンシャル φ(r) = 4ε[(σ/r)¹² − (σ/r)⁶] から得られるばねとなり、縦成分 κ_L = φ″(r) および 横成分 κ_T = φ′(r)/r をもつ。動力学行列

D_ab(k) = Σ_j [ (κ_L − κ_T)·û_a û_b + κ_T·δ_ab ]·(1 − cos(k·R_j))

はカットオフ(2.0…3.0 nn)までのすべての隣接殻にわたって和をとる。各格子はあらかじめそれ自身の LJ 平衡へ緩和される(自由な格子定数。同じ密度を強制しない)。これは簡潔な nn+2nn 近似(第6.8章)が 剛性の順序を正しく捉えているかどうかの正直な検証である。

結果(決定論的、標準ライブラリのみ):

二つの所見、いずれも第6.6/6.8章を強める:

1. FCC は安定な LJ 基底状態である(カットオフ 3.0 nn で E ≈ −8.4 ε/原子、BCC の −8.0 ε に対して)— 対ポテンシャルにおける密充填の既知の優位である。
2. 理想 LJ BCC はせん断不安定である:フォノンモードのおよそ 3…3.5 % が虚(ω² < 0、ボルン不安定)。 純粋な対ポテンシャルの理想像は機械的に安定な BCC をまったく支えない。残る実モードにわたって θ_BCC/θ_FCC ≈ 1.0 であり、μ 比較から必要な開口しきい値 0.61 はなおさら下回られない。

こうして第6.6/6.8章の主張は、打ち切りばねから完全格子へ、また純粋な縦モデルから完全な動力学行列へという 精密化に対して頑健である:古典的な理想像では純粋に幾何学的な固体(1)→固体(2)転移は生じない — 現実の 多形(鉄 第7.13章、チタン 第7.15章)は対ポテンシャルの外の寄与(電子的/磁気的、非調和的)を必要とする。 本格的な第一原理フォノン計算は意識的に対象外とする(簡潔な理想像の性格から外れてしまう)。

<a id="kap7"></a>

7 — 較正と状態方程式

要約。 第7章はモデルをただ一つのエネルギーパラメータ ε によって実在の物質へ較正し、液体/気体の枝を 状態方程式へと拡張する:

• 7.1 / 7.2 較正法と較正表: 融点に一つの錨点をもつ線形の温度対応 T = (ε/k_B)·τ。そこから精密化された 普遍比 R = τ_b/τ_m ≈ 1.54。
• 7.3 / 7.4 硬核ファン・デル・ワールス状態方程式: 換算状態方程式 p_r = 8τ_r/(3v_r − 1) − 3/v_r²、 臨界点 (1, 1, 1) とマクスウェル共存。a ∝ ε, b ∝ D²。沸点は系統的に T_b/T_c ≈ 0.48 にあり (グルドベルクとの隔たり)、レドリッヒ–クォンによる比較計算はこれを ≈ 0.60 へ引き上げる — 状態方程式を 替えることなく、名指しされた誤差予算として。
• 7.5 共通の (τ, v) 面: 共存ドーム(沸騰)、融解線、硬い壁が三重点とともに一つに組み合わさる (v_solid = 0.40 に対し τ_r,triple ≈ 0.47)。
• 7.6 潜熱: L ≈ 1.41 ε は ΔS_fus ≈ 1.175 R を与える(リチャーズ則に近い)。金属は見事に当たり (銅 ×1.00)、イオン結晶・水素結合結晶・分子結晶は理想像の外にある。
• 7.7 ε ↔ ε_lg の統一: k_B T_c = z_max·ε/4 が両尺度を結ぶ(ε_lg = (27/32)·z_max·ε)。三次元の予測 τ_c/τ_m = 2.50 は単純な分子性物質に当たり、二次元は特別な場合のままである。
• 7.8 / 7.9 パラメータ不要の三重点: v_solid は連結性からではなく剛性しきい値 z̄_c = 2d から従う。 3D では清潔な液体窓が開き(τ_m = 1.20 < τ_triple = 2.83 < τ_c = 3.00)、2D では閉じたままである (三重に裏づけられた二次元の特殊な役割)。
• 7.10 / 7.11 本物の三次元走査: FCC/BCC 上のボンドパーコレーションは z̄_c(3D) ≈ 1.5–1.6 を与える。 結合速度論の τ 走査は τ_m ≈ τ_rig(融解 = 剛性、低い)と τ_b(3D) ≈ τ_c(沸騰 = 臨界点における連結性 の崩壊)を示す。
• 7.12 / 7.13 / 7.15 固体(1)→固体(2)線: 点 τ_ss から圧力を通じて一本の線になり(dp/dτ = ΔS/Δv > 0, Δv = +8.9 %)、(τ, v) 面における第四の境界線として、また 鉄(γ→δ)に較正された (τ, p) 線として、 チタン(α→β)で確認され、アルゴンを陰性対照として。
• 7.14 レドリッヒ–クォン枝: ファン・デル・ワールスの既定値の隣に、独自の較正をもつ自立した理想像の 変種として(T_b/T_c ≈ 0.60、誤差 −1 %)。

<a id="kap7-1"></a>

7.1 較正法(エネルギーパラメータ ε 一つ)— 調整済み

線形に:T = (ε/k_B)·τ、物質ごとに一つの錨点。融点での錨づけ:ε/k_B = T_m,実験/τ_m、予測 T_b,pred = T_m,実験·(τ_b/τ_m)。精密化された融解ピーク(協同項)により、c_v 最大はいま τ_m = 1.20 にあり(1.15 ではなく)、連結性パーコレーション(沸騰近傍)は変わらず τ_b = 1.85 にある。そこから 精密化された普遍比:

R = τ_b / τ_m = 1.85 / 1.20 ≈ 1.54 (版08:1.61)。

版08の較正表(第7.2章)は R とともに線形にスケールする:すべての T_b,pred は係数 1.54/1.61 ≈ 0.957 だけ下がる(≈ −4.3 %)。これは希ガスの過大評価をわずかに改善し(アルゴン +54 % → +48 %)、NaCl と鉄の 当たりをわずかに悪化させる(NaCl −0.6 % → −4.8 %)。基本的な主張は変わらない:R は等方的に密充填し、 無指向的に結合する物質に合う。希ガスや鎖状物質は理想像の外にある。

<a id="kap7-2"></a>

7.2 較正表

modell_08.md 第7.2章から変更なく引き継いだ(7.1 の線形 R スケーリングの注記つき)。cv_scan_08.py(表) あるいは cv_scan_09.py(τ_m, τ_b, R)で再現できる。

<a id="kap7-3"></a>

7.3 密度/体積軸:硬核ファン・デル・ワールス状態方程式

圧力は「パラメータ」から真の変数へと格上げされる。硬核 が固有体積 b を、結合の引力(ε·z̄)が凝集圧 a を 与える。換算単位では(臨界点が (v_r, τ_r, p_r) = (1, 1, 1)):

p_r = 8 τ_r / (3 v_r − 1) − 3 / v_r², v_r > 1/3(硬い壁)。

マクスウェル作図による共存(cv_scan_09.py、等面積の切断):

結果。 τ_c の下では液/気の密度の跳びが開く(v_gas/v_liq は 2.5 から 11.5 へ増え、冷たいほど大きい)。 これが沸騰線であり — いまや*固定*密度の格子の上に生きる融解の c_v ピーク(第5.3章)から明確に分離 された。これで版08の未決点3が閉じる:融解(固定密度における剛性/結合の秩序)と沸騰(密度/体積の跳び)は 同じモデルの異なる軸である。

<a id="kap7-4"></a>

7.4 定量的な液/気較正:a → ε, b → D²

換算量の実在物質への写像。 ファン・デル・ワールス定数は臨界点に依存する:

b = k_B·T_c / (8 p_c) (固有体積、硬核), a = 27 k_B²·T_c² / (64 p_c) (凝集、結合)。

したがって a ∝ ε(凝集は結合エネルギーに由来する)であり、b ∝ D²(二次元の固有体積は硬核直径から) である。液/気の枝の自然なエネルギー尺度は

ε_lg / k_B = (a/b) / k_B = (27/8)·T_c = 3.375·T_c.

錨物質アルゴン(T_c = 150.87 K, p_c = 48.63 bar):b = 5.35·10⁻²⁹ m³/分子、 a = 3.76·10⁻⁴⁹ J·m³/分子²、ε_lg/k_B = 509 K、分子直径 D ≈ b^{1/3} = 3.8·10⁻¹⁰ m — アルゴンのファン・ デル・ワールス直径(~3.4·10⁻¹⁰ m)に対して 11 % の精度である。

沸点の予測。 一定の外圧のもとで沸点は、換算飽和圧が p_sat(τ_r)·p_c = p_atm を満たすところにある。 そこから T_b = T_c·τ_r,b(lg_calib_10.py のマクスウェル飽和曲線):

*CO₂ は 1 atm で昇華する(T = 194.7 K = 昇華点)。

結果。 分子性物質(p_c ≈ 30–70 bar)についてモデルは一貫して T_b/T_c ≈ 0.48–0.51 を与え、測定値より 系統的に ~17–24 % 低い。これは狭い意味でのモデルの誤りではなく、ファン・デル・ワールス方程式が経験的な グルドベルク則(T_b/T_c ≈ 0.6)から隔たっているという既知の事実である:ファン・デル・ワールス等温線は 低温で飽和圧を過小評価する。水(および T_c/p_c がもともと粗く見積もられているだけの金属)は p_r,atm < 0.005 で意味のあるファン・デル・ワールス領域の外にある(τ_r,b < 0.40)— そこでは凝集がもはや ファン・デル・ワールス的ではない(水素結合、金属結合)。沸点の桁と順序は正しく当たっている。

誤差表としてのグルドベルク位置づけ。 この隔たりはここでは名指しされるだけで、状態方程式の交換によって 是正されはしない(理想像はファン・デル・ワールス的なままである)。経験的なグルドベルク則(Guldberg 1890) は、沸点がおよそ T_b/T_c ≈ 0.6 にあると述べる(物質依存はわずかで ~0.58–0.64)。換算ファン・デル・ワールス モデルは代わりに一定の値 T_b/T_c = τ_r,b ≈ 0.48–0.51 を与える:

解釈。 モデルの平均値 0.49 はグルドベルク定数 0.60 におよそ −18 % で当たる。ずれは系統的かつ一方向的 (常に低すぎる)であり、したがって分子像ではなくファン・デル・ワールス等温線の性質である(低い τ での飽和圧が 小さすぎる)。より良い凝集の状態方程式(たとえばレドリッヒ–クォン)は τ_r,b を上げるであろう — が、ここでは 意識的に理想像の外にとどめる。

レドリッヒ–クォンによる比較計算を用いた誤差予算(guldberg_rk_15.py)。 ファン・デル・ワールスの グルドベルク則からの隔たりを数値化するため(治すためではない)、同じ沸点テストをレドリッヒ–クォン型の、 温度依存の凝集 a/√T で繰り返す。換算単位で(臨界点 (v_r, τ_r, p_r) = (1, 1, 1)、Z_c = 1/3):

p_r = 3 τ_r/(v_r − β) − A/(√τ_r · v_r(v_r + β)), β = 3·Ω_b = 0.25992, A = 9·Ω_a = 3.84732,

マクスウェルの等面積作図と、上と同じ換算外圧 p_r,atm を用いる。結果(錨物質ごとに一つの数値、 T_b/T_c = τ_r,b):

結果。 温度依存の RK 凝集は T_b/T_c を vdW ≈ 0.48 から ≈ 0.60 へ引き上げる — 実質的に グルドベルク定数であり実測値(≈ 0.58–0.63)である。測定値に対する平均誤差は −21 %(vdW)から −1 %(RK)へ下がる。これでグルドベルクとの隔たりの原因は明確に数値化された:それはまさに ファン・デル・ワールス凝集の温度依存性の欠如(低い τ での飽和圧が小さすぎる)であって、分子像ではない。 この比較計算は意識的に理想像の外にある。モデルはファン・デル・ワールス的なままである(状態方程式の 交換なし)— 誤差予算は名指しされただけで、組み込まれてはいない。

<a id="kap7-5"></a>

7.5 三重点を含む共通の (τ, v) 面

融解軸と沸騰軸はこれまで別々に生きていた:融解(剛性/結合の秩序)は固定密度で、沸騰(ファン・デル・ ワールス)は密度の跳びを通じて。共通の (τ_r, v_r) 平面では両者が組み合わさる:

• 共存ドーム(沸騰):マクスウェル作図から得られる v_liq(τ_r)(左枝)と v_gas(τ_r)(右枝)。両者は 臨界点 (τ_r, v_r) = (1, 1) で出会う。
• 固体(融解):硬い壁 v_r > 1/3 のすぐ上の狭い帯。高く、ほぼ固定した密度で(換算体積 v_solid ≈ 0.40)。
• 融解線:v_solid でほぼ垂直(固定密度)、すなわちマクスウェル剛性の温度(z̄_c = 2d。第7.8章は そこから v_solid をパラメータ不要で導く — 旧版11の対応づけ「z̄_c ≈ 2.08」はこれで訂正された)。

三重点は、液体の共存境界 v_liq(τ_r) が固体の充填密度 v_solid に達する点である — その下ではもはや液体の 枝は存在しない(直接昇華):

結果。 v_solid = 0.40 に対して三重点は τ_r,triple ≈ 0.47、p_sat ≈ 0.018 にあり、蒸気体積は v_gas ≈ 68(きわめて希薄な昇華蒸気)である。比 τ_triple/τ_c ≈ 0.47 は現実の窓の内にある(水 0.42、 アルゴン 0.56、CO₂ 0.71)。固体が密であるほど(v_solid が小さいほど)換算三重点は低く、共存する蒸気は 希薄になる — 定性的に正しい。これで両方の転移と三重点は一つの (τ, v) 面に錨づけられた。臨界点は (1, 1, 1) に座る。版15は第7.12章に第四の境界線を加える — 固体(1)/固体(2)変態(FCC↔BCC)を、 τ_ss で互いに切り替わる二本のほぼ垂直な固体密度帯(v ≈ 0.667 および 0.726)として。

これ以前の開発段階: ここでは v_solid = 0.40 はまだ*置かれた*ものだった。第7.8章は v_solid を パーコレーションしきい値 z̄_c からパラメータ不要で導き、三重点を真の予測にする。

<a id="kap7-6"></a>

7.6 潜熱と融解エンタルピーの対比

モデルの潜熱 L ≈ 1.41 ε(第5.3章)は融点較正 ε/k_B = T_m/τ_m(τ_m = 1.20)を通じて実在の単位に翻訳される:

ΔH_m,pred = L·R·T_m/τ_m = 1.175·R·T_m, これに応じて ΔS_fus,pred = L·R/τ_m = 1.175·R。

これは — 注目すべきことに — ほぼ正確にリチャーズ則である(単純な構造の固体では ΔS_fus ≈ R)。測定された 融解エンタルピーとの比較(lg_calib_10.py):

結果。 モデルは桁(kJ/mol)を一貫して当て、単純な金属についてリチャーズ則を見事に当てる (銅 ×1.00、鉄 ×1.28)。液化ガスはそのすぐ下にある(メタン ×0.94、窒素 ×0.87、アルゴン ×0.69)。 イオン結晶、水素結合結晶、分子結晶(NaCl ×0.37、水 ×0.44、ベンゼン ×0.28、CO₂ ×0.23)は現実の融解 エントロピーが明らかに高い(ΔS_exp/R = 2.6–5.0)。融解時に追加の秩序自由度(イオン対形成、多原子分子の 配向)が解放され、それを一つの秩序パラメータのモデルは捉えないからである。限界はこうして明確に名指し された:理想像は等方的で無指向的に結合する物質に妥当する(フィッティングなし)。

<a id="kap7-7"></a>

7.7 エネルギー尺度の統一:ε ↔ ε_lg

これまで融解の枝は結合エネルギー ε(τ_m = 1.20)を担い、液/気の枝は別個の凝集尺度 ε_lg/k_B = 3.375·T_c(第7.4章)を担っていた。両者は一つのモデルパラメータから導ける。ファン・デル・ ワールス凝集 a はモデルでは分子あたりの結合エネルギーの総和から生じる(結合):各分子はエネルギー ε の 結合を z_max/2 本担う。隣接引力 ε をもつ平均場ブラッグ–ウィリアムズ格子気体は、臨界温度についての 古典的な関係を与える:

k_B·T_c = z_max·ε / 4, すなわち τ_c ≡ k_B T_c/ε = z_max/4.

これを ε_lg = a/b = (27/8)·k_B·T_c(第7.4章)に代入すると、両尺度の結合が従う:

ε_lg = (27/8)·(z_max/4)·ε = (27/32)·z_max·ε, 定数 c = 27/32 = 0.844.

これで融解ピーク、パーコレーション、臨界点は一つの ε 尺度上に乗る(unify_scale_12.py):

予測 τ_c/τ_m = z_max/(4·τ_m): 2D → 1.25、3D → 2.50。実在物質との対比:

結果。 三次元の予測 τ_c/τ_m = 2.50 は単純な分子性物質の現実の指標(T_c/T_m ≈ 2.0–2.1)によく当たる。 金属はより高い(T_c は粗い見積もり、結合の種類が異なる)。二次元の予測(1.25)は低すぎる — しかも τ_c(2D) = 1.50 は τ_b = 1.85 よりさらに下にある(臨界点がパーコレーションしきい値の下に来てしまう)。 これは計算違いではなく、既知の二次元の特殊な役割である:そこでは凝集尺度と結合尺度が明確に並ばず、 頑健な液体窓が欠ける(連続的な二次元融解、第5.5章を参照)。3D になって初めて(τ_m = 1.20 < τ_c = 3.00) 尺度が物理的に意味のある形で分離する。

<a id="kap7-8"></a>

7.8 z̄_c から得るパラメータ不要の三重点

第7.5章では v_solid = 0.40 をまだ置いていた。これ以前の開発段階は固体の充填密度をパーコレーション しきい値から導く。換算体積 v_r における格子の配位数は密度–配位関係に従う(硬い壁 v_r = 1/3 ↔ 最密充填、 z̄ = z_max):

z̄(v_r) = z_max·(1/3)/v_r = z_max/(3 v_r), 逆に v_solid = z_max/(3·z̄_c).

固体は z̄ がしきい値 z̄_c を下回るとただちに剛性を失う。それに対応する密度が v_solid である。三重点は そのとき、液体の共存境界 v_liq(τ_r) がこの密度に達する点である(unify_scale_12.py、2D、z_max = 6):

改善版。 融解密度 v_solid は一般的マクスウェル剛性 z̄_c = 2d が定める(2D → 4、3D → 6)。ボンド 連結性が定めるのではない(理由は第5.4a章:融解は*剛性*の喪失であって*連結*の喪失ではない)。したがって 信頼できる三重点はこうなる(rigidity_3d_13.py、2D):

二次元の結果。(物理的に正しい)剛性しきい値をもってしても二次元単一格子は清潔な液体窓を与えない: v_solid = 0.50 はたしかに現実の指標 τ_triple/τ_c ≈ 0.77 を与えるが、絶対的な ε 単位では τ_triple = 0.77·1.50 = 1.16 であり、こうして τ_m = 1.20 の下にある。連結性しきい値(z̄_c = 2.08)は v_solid ≈ 0.96 ≈ v_c を共存ドームの外へ押し出しさえする — 三重点と臨界点が一致してしまう。両方の計算が 二次元の特殊な役割(第5.5章、7.7章)を裏づける:二次元単一格子は自立した液相をほとんど担わない。これは 正直な所見であり、狭い意味での計算やモデルの弱さではない。

したがって決定的なテストは三次元の検算(第7.9章)である:そこでは同じ剛性則(z̄_c = 6 → v_solid = 0.667) が清潔な窓を解き放つ。いずれの場合も三重点はいまやパラメータ不要で z̄_c に結びついている — 自由 パラメータ v_solid は消去された。

<a id="kap7-9"></a>

7.9 三次元検算:単一の ε 尺度上の明瞭な液体窓

すべての目印が一つの ε 尺度上に乗るとき、二次元の崩壊が 3D で解消されるかどうかの検証 (第7.7章からの τ_c = z_max/4;一般的剛性 z̄_c = 2d からの v_solid;三重点は、換算ファン・デル・ワールス 液体密度 v_liq(τ_r) が固体密度に達するところ)。結果(rigidity_3d_13.py):

結果。 3D では清潔な液体窓が開く:τ_m = 1.20 < τ_triple = 2.83 < τ_c = 3.00(窓の幅: τ_triple − τ_m = 1.63;τ_c − τ_triple = 0.18)。上端は狭い — 三重点は臨界点のすぐ下にある — が、順序は 正しく、液相は自立した領域として存在する。2D では窓は閉じたままである(τ_triple = 1.16 < τ_m)。これで 二次元の特殊な役割は三度目に独立して裏づけられ、同時にそれが次元効果であることが示された:同じ剛性則 が 3D ではモデルを癒す。狭い上端は第7.10章で本物の 3D パーコレーション走査と感度解析によって補強される。

モデル限界の囲み — 二次元の特殊な役割。 三つの独立な計算が同じ所見を示す:二次元単一格子は自立した液相をほとんど担わない。 1. τ_c < τ_b(第7.7章):臨界点(1.50)が連結性パーコレーション(1.85)の下に来てしまう — 凝集尺度と結合尺度が 2D では明確に並ばない。 2. v_solid ≈ v_c(第7.8章):連結性による融解密度(0.96)が臨界密度とほぼ一致する — 三重点と臨界点が 崩れて重なる。 3. τ_triple < τ_m(第7.9章):正しいマクスウェル剛性を用いてすら三重点(1.16)が融解ピーク(1.20)の 下にある — 窓がない。

これは正直な所見であり、計算やモデルの弱さではない:既知の次元効果である(連続的な二次元融解、 KTHNY 型;第5.5章)。3D では同じ剛性則がモデルを癒す(第7.9/7.10章)。第二の近距離秩序(第6.4章)が 二次元像をさらに改善しうるかどうかは拡張の問いであって、理想像の欠陥ではない。モデルは可視化に資するもの であるから、私はこの弱さを受け入れ、モデルの公理的な拡張を断念する。

<a id="kap7-10"></a>

7.10 本物の三次元ボンドパーコレーション走査(FCC/BCC)と窓端の感度

第7.8/7.9章は二次元の連結性値 z̄_c = 2.08 を観念的に 3D へ移し(そこでは v_solid = 1.92 で非物理的)、3D の 主張を一般的剛性 z̄_c = 2d に支えていた。版13は連結性を 3D でいま直接検証する:FCC(z_max = 12) と BCC(z_max = 8)上のボンドパーコレーションを Union-Find によるモンテカルロで、周期境界、決定論的な LCG シードで(rigidity_3d_14.py)。しきい値 p_c = P_∞(p) = 0.5 の交差(第5.6章の二次元有限サイズ研究と 同じ方法論):

MC 値は控えめな箱の辺長(L = 12、純 Python)のために文献のしきい値よりわずかに上にある(有限サイズ バイアス。2D と同様、推定量に応じて下からまたは上から);順序と桁は合っている(FCC はより密 → BCC より 小さい p_c)。

解釈(独立な三次元的裏づけ)。 三次元の連結性しきい値は z̄_c(3D) ≈ 1.5–1.6 にあり — 二次元の値 2.08 よりさらに低く、マクスウェル剛性 2d = 6 よりはるかに下である。そこから従う「連結性の融解密度」 v_solid = 1/(3 p_c) ≈ 1.7(BCC)ないし 2.6(FCC)は共存ドーム(v_liq,max ≈ 0.94)のはるか外にある。 これで — 二次元の値を 3D に移すことなく — 直接示された:3D では連結性はなおさら融解密度を定めえない。 それは剛性から来なければならない(z̄_c = 6 → v_solid = 0.667)。二つのしきい値は 3D では 2D よりさらに 明確に分離している。本物の τ_b(3D)(z̄(τ) が低い連結性しきい値 1.5 に達する温度)は、したがって二次元の 値 τ_b より上、臨界点のごく近くに来るであろう — 「連結性の崩壊 = 沸騰」という像と整合的である。

狭い上側の窓端の感度(τ_c − τ_triple ≈ 0.18)。v_solid = 0.667 のまわり(±2 %)で数値的に dτ_triple/dv_solid ≈ 1.55 が得られる。窓が閉じる(τ_triple = τ_c = 3.00)のは v_solid が ≈ 0.11 上がったとき(≈ 0.78 へ)であり、すなわち実効的な三次元剛性しきい値が 2d = 6 より明らかに上にある 場合(z̄_c ≳ 7.7)である。一般的マクスウェル数え上げ(z̄_c = 2d = 6)に対して、順序 τ_m < τ_triple < τ_c は したがって頑健である。上端は実際に狭いが、危うくはない。離散的な格子交替 固体(1)→固体(2)(FCC↔BCC)は 版14で定量化され(第6.3章)、結合速度論の本物の三次元 τ 走査は第7.11章に続く。

<a id="kap7-11"></a>

7.11 結合速度論の真の三次元 τ 走査:τ_b(3D) と単一 ε 尺度上の融解線

第7.8–7.10章は 3D の主張を一般的剛性(z̄_c = 2d = 6)と静的な連結性走査(z̄_c(3D) ≈ 1.55)に支えていた。 未決だったのは、3D における結合速度論そのものの τ 走査である:自己無撞着な z̄(τ) は*温度軸のどこで* この低い連結性しきい値に達するのか(= 本物の τ_b(3D))、そして融解ピークはどこにあるのか。版14は自己無撞着な z̄(τ)(第4.2章、x = 1/(1 + G·e^{−(ε+K·x)/τ}])を直接 FCC(z_max = 12)と BCC(z_max = 8)上で、 G = z_max と二次元参照の場合と同じ協同性 J = 0.55 で計算する。n_blk(τ) は FCC/BCC のボンド モンテカルロ(Union-Find)から得られ、

U/N = τ·[3/2 + z̄/4 + (3/2) n_blk] − (ε/2)·z̄ − (K/2)·x² (3D、協同項つき)、

c_v = dU/dτ(bondkinetics_3d_15.py、決定論的な LCG シード。n_blk は平滑化してある。それに依存するのは 弱いブロック回転項だけだからである)。三つの目印が読み取られる:c_v ピーク τ_m(融解)、剛性の目印 τ_rig(z̄(τ) が 2d = 6 を横切る)、および本物の τ_b(3D)(z̄(τ) が第7.10章の連結性しきい値 z̄_c(3D) を横切る):

結果。 二つの所見。いずれもこれまでの像と整合的である:

1. 融解は低い位置にあり、二重に裏づけられている。 熱力学的な融解ピーク τ_m(FCC 1.70、BCC 1.40)は 実質的に剛性の目印 τ_rig(z̄ = 6;FCC 1.73、BCC 1.35)と一致する — 網が剛性を失うのと同じ温度である。 これが三次元の融解線を ε 尺度に錨づける:T_m = (ε/k_B)·τ_m、たとえばアルゴン(T_m = 83.8 K)なら ε/k_B = 83.8/1.70 ≈ 49 K。K = J·z_max は 3D で 2D より大きい(3.30 ではなく 6.60)から、三次元の融解 ピークは二次元の値 1.20 より高い — 矛盾ではなく、より密な充填のより強い協同性である。
2. 本物の τ_b(3D) は臨界点に座る。 三次元の連結性しきい値がきわめて低い(≈ 1.55)ため、z̄(τ) が それに達するのはごく上である:τ_b(3D) = 3.18(FCC)ないし 2.85(BCC)— τ_c = z_max/4(3.00 ないし 2.00)と同じか、わずかに上である。すなわち結合の巨大クラスターは液相全体を生き延び、臨界点に なって初めて崩壊する。これはまさに第7.10章の予測(「本物の τ_b(3D) は二次元の値より上、臨界点の ごく近くに来るであろう」)であり、「沸騰 = 連結性の崩壊」という像を独立に裏づける:3D では物質は 臨界点の間近まで貫通して結合したままである。

両方の所見は第5.4a章のしきい値の分離を強める:融解 ↔ 剛性(低い、τ_m ≈ τ_rig)、沸騰 ↔ 連結性 (高い、τ_b(3D) ≈ τ_c)。融解密度はしたがって剛性から来る(v_solid = 0.667)のであって連結性からではない。 二つのしきい値は 3D では温度軸上で、このモデルでかつてないほど遠く引き離されている。

<a id="kap7-12"></a>

7.12 (τ, v) 面における固体(1)→固体(2)変態線

第6.3章は FCC↔BCC の格子変態を点 τ_ss(p = 0 において)として与えた。圧力を通じてそこから一本の 線が生じ(polymorph_line_16.py)、それは第7.5/7.9章の共通の (τ, v) 面に第四の境界線として属する — 共存ドーム(沸騰)、融解線、硬い壁と並んで。

体積の跳び。 二つの固体はほぼ固定した密度に座る。両者の体積差は充填密度から従う(v ∝ 1/η): η_FCC = 0.7405、η_BCC = 0.6802。v_FCC = 0.667(第7.9章)とすると v_BCC = v_FCC·η_FCC/η_BCC = 0.726 となり、すなわち

Δv = v_BCC − v_FCC = +0.059(+8.9 %)。

より開いた BCC のほうが大きな体積をもつ — 再結晶 FCC→BCC は体積の増大を伴う(漫画「再結晶」で 示したのと同じ主張である)。

クラウジウス–クラペイロン。 ΔS = S_BCC − S_FCC(高温相引く低温相)とすると、傾きは

dp/dτ = ΔS/Δv.

両方の量は正である(BCC は高温相であるから、より多くのエントロピーを*かつ*より大きな体積をもつ)。したがって dp/dτ > 0:より高い圧力は変態温度を上げる。より密な FCC を有利にするからである。参照 θ_BCC/θ_FCC = 0.60(τ_ss = 1.333、ΔS = 0.50 k_B、L = 0.667 ε)については:

参照に対する (τ, p) の線(上昇):p = 0 → τ_ss = 1.33;p = 2 → 1.57;p = 4 → 1.81;p = 6 → 2.04 (換算単位)。

(τ, v) 面では変態はほぼ垂直な二本の固体密度帯として現れる — 固体(1) = FCC が v ≈ 0.667(より密、 低い T)、固体(2) = BCC が v ≈ 0.726(より開いた、高い T)であり、τ_ss で互いに切り替わる。その上、τ_m で 融解が続く。τ 軸上の完全な順序はこうして τ_ss < τ_m < τ_c ≤ τ_b(3D) である(第7.11章)。

実在の単位。 融点較正 ε/k_B = T_m/τ_m(τ_m = 1.70、3D-FCC)を通じて c_v 段差と変態熱が手に取れる ものになる。c_v 段差は較正に依存しない:|Δc_v| = 1.0 k_B = 1 R = 8.31 J/(mol·K) = 0.0083 kJ/(mol·K)。変態熱 ΔH_ss = L·R·(ε/k_B)(参照 L = 0.667 ε):

変態熱は小さく正であり(0.2–0.9 kJ/mol の桁)、したがって融解エンタルピー(第7.6章、数 kJ/mol)より 明らかに下にある — 小さな c_v 段差を伴う弱い多形変態として期待されるとおりの像である。(θ 比 0.60 は 意識的に厳密な理想像の外にある。第6.6章は、純粋な幾何なら θ ≈ 0.82–0.94 となり変態は閉じたままであろうことを 示す — この線は軟らかい BCC 相の*もし〜ならば*の帰結である。)

<a id="kap7-13"></a>

7.13 鉄で較正した (τ, p) 図における固体(1)→固体(2)変態線

第7.12章は変態を (τ, v) 線として導く。版16はそれを (τ, p) 図における線として補い、実在の FCC↔BCC 多形で較正する(solidsolid_pt_17.py)。

硬核ファン・デル・ワールス状態方程式からの圧力場。 二つの固体密度 v_FCC = 0.667 と v_BCC = 0.726 (第7.12章)において、換算状態方程式 p_r = 8τ_r/(3v_r − 1) − 3/v_r²(第7.3章)は一貫して正の圧力を 与える。同じ τ_r ではより密な FCC がより高い圧力を担う(τ_r = 1:p_r = 1.25 対 1.10)。すなわち圧縮は FCC を有利にする — 正のクラウジウス–クラペイロン傾き dp/dτ = ΔS/Δv > 0 と整合的である。

線と変態点。 参照 θ_BCC/θ_FCC = 0.60(τ_ss = 1.333、ΔS = 0.50 k_B、dp\/dτ = 8.47、換算単位 p\ = p·D³/ε)に対して線 p\(τ) = 8.47·(τ − 1.333) は上昇する:p\ = 0 → τ = 1.33;p\ = 2 → 1.57; p\ = 4 → 1.81;p\ = 6 → 2.04。**p\_tri ≈ 3.1 で固体–固体線が融解の目印 τ_m = 1.70 に出会う — (τ, p) 平面における固体(1)–固体(2)–液体変態点**(三重点の類似物)である。

鉄での較正(ε/k_B = T_m/τ_m = 1811 K / 1.70 = 1065 K)。鉄には逆符号の FCC↔BCC 転移が二つある — まさに試験例である:

γ→δ について(τ_ss = 1667/1065 = 1.565 → θ_BCC/θ_FCC ≈ 0.61)モデルは符号と桁を当てる: ΔH_ss ≈ 4.2 kJ/mol(実際は ≈ 0.84 kJ/mol — 桁は kJ/mol、鋭いしきい値近似のため値は ~5 倍高い)、c_v 段差は パラメータ不要で 1 R = 8.31 J/(mol·K)、dT/dp ≈ +245 K/GPa。これに対し α→γ 転移は密な相(γ)が 高温相である — これは純粋な充填幾何では与えられない(密 = 低エネルギー = 低い T)。α 鉄は磁気的に 安定化されている。鉄はこうして正直な所見の両側を裏づける:幾何学的トポロジー(γ→δ)は理想像の内にあり、 磁気的異常(α→γ)はその外にある — 第6.6/6.8章が θ ≈ 0.87…1.00 > 0.61 として定量化するのと同じ「外」である。

<a id="kap7-14"></a>

7.14 分離した理想像の枝としてのレドリッヒ–クォン

ファン・デル・ワールス枝のグルドベルクとの隔たり(第7.4章、T_b/T_c ≈ 0.60 ではなく ≈ 0.48)は、そこでは 名指しされただけであった。明示的な要望により、改善版はレドリッヒ–クォン凝集を独自の較正、独自の共存、 独自の誤差節をもつ自立した理想像の変種として導入する(rk_branch_17.py)— ファン・デル・ワールス枝が 既定値のままである。モデルのそれ以外の物理(融解、パーコレーション、固体(1)→固体(2))は凝集の状態方程式の 選択に左右されない。

独自のエネルギー尺度。 a_RK = Ω_a·R²T_c^{2,5}/p_c と b_RK = Ω_b·R T_c/p_c から、この枝の凝集エネルギー は次となる:

ε_lg,RK/k_B = (Ω_a/Ω_b)·T_c = 4.93·T_c (vdW の ε_lg/k_B = 27/8·T_c = 3.375·T_c に対して)、

より大きい。温度依存の凝集 a/√T が因子 √T_c を畳み込むからである。b_RK から従う直径は、ファン・デル・ ワールス枝と同程度によく実測値を当てる(アルゴン D ≈ 3.34·10⁻¹⁰ m)。

独自の共存(換算、臨界点 (1, 1, 1)、Z_c = 1/3; p_r = 3τ_r/(v_r − β) − A/(√τ_r·v_r(v_r + β))、β = 0.25992、A = 3.84732):

独自の誤差節(p_r,atm = p_atm/p_c における沸点、T_b/T_c = τ_r,b):

結果。 RK 枝は T_b/T_c を vdW ≈ 0.48 から ≈ 0.60 へ引き上げ、こうしてグルドベルク定数と実測値を実質的に 当てる(平均誤差は vdW の −21 % に対し −1 %)。それはファン・デル・ワールス既定値の隣に第二の、 明確に分離した理想像の変種として立つ。どちらの枝が「より正しい」かは要求される精度の問題であって、モデル 構造の問題ではない — モデル3の核(秩序パラメータとしての z̄)は両方の枝で同一である。

<a id="kap7-15"></a>

7.15 他の多形における固体(1)→固体(2):チタンと希ガス対照

第7.13章は固体(1)→固体(2)の (τ, p) 線を鉄で較正する。改善版は符号と桁を第二の、磁気的に目立たない 多形で確保し、陰性対照を補う(solidsolid_titanium_18.py)。較正と体積の跳びは第7.13章と同じである (ε/k_B = T_m/τ_m、τ_m = 1.70;充填幾何から Δv = +8.9 % 密→開)。hcp は FCC と同様に z = 12 で近似する (fcc/hcp の差は名指しするだけで、解消しない)。

チタンは 1155 K における α(hcp) → β(bcc) でまさにモデルのトポロジーをもつ:密な相(hcp、z = 12) が低温相であり、加熱すると開いた相(bcc、z = 8)になる。すなわち ΔS > 0 かつ Δv > 0 → dp/dτ > 0 — 鉄の γ→δ と同じ符号であるが、磁気的異常はない。二つの独立な密→開の多形の比較:

いずれも符号(dp/dτ > 0、圧力が密な相を有利にする)と桁(ΔH ~ kJ/mol、c_v 段差はパラメータ不要で 1 R = 8.31 J/(mol·K)、dT/dp ~ 10²·K/GPa)を当てる。ΔH_ss の値は両者とも系統的に~一因子だけ高い — μ 比較の鋭いしきい値近似が潜在変態熱を過大評価している。これは名指しされたモデルの限界であって、符号の 誤りではない。

アルゴンは陰性対照である:固体アルゴンは FCC であり、融解(83.8 K)まで FCC のままである — 常圧では安定な bcc は存在しない。まさにそれを純粋な対ポテンシャルの理想像は予測する: θ_BCC/θ_FCC ≈ 0.9…1.0 > 0.61 かつせん断不安定な LJ-BCC(第6.9章)⇒ 固体(1)→固体(2)転移はない。 この対照実験は所見の構造を裏づける:理想像が固体(1)→固体(2)線を生むのは、非幾何学的な寄与(鉄やチタンに おけるような電子的/磁気的なもの)が BCC をそもそも安定化するときだけである。チタンは肯定的な、 アルゴンは否定的な証拠を、第6.8/6.9章の同じ正直な枠組みに与える。

<a id="kap8"></a>

8 — 文献比較:静的熱力学のモデル観

要約。 第8章はモデル3を静的熱力学の確立したモデル観の中に位置づける:

• 8.1 一覧表: 十一のモデル(理想気体、ファン・デル・ワールス、レドリッヒ–クォン、ビリアル展開、 アインシュタイン、デバイ、レナード–ジョーンズ/剛体球、セル理論、イジング/格子気体、パーコレーション、 フローリー–ハギンズ)を相と統計熱力学的な基本着想とともに示す。
• 8.2 個々のモデル: それぞれ特徴・応用・限界を一段落で、典拠へのリンクから始めて述べる。自由に アクセスできる出典は model3/links/ 以下に保存してあり、セル理論(有料)は引用のみである。
• 8.3 対比: モデル3は秩序パラメータ z̄ をもつ協同的二状態結合モデルである。最も近い親戚は イジング/格子気体であり、沸騰はパーコレーションから直接受け継ぎ、フォノン像(アインシュタイン/デバイ)が 融解の剛性を担う。モデル3はこれら三つの伝統の系統を共通の ε 尺度上で連結する。

本章はモデル3を、固体・液体・気体の平衡(「静的」)熱力学の確立したモデル像の中に位置づける。まず表に よる一覧(8.1)、次に各モデルについて典型的な特徴・応用範囲・限界を述べる短い散文の段落 — 各段落は典拠への リンクで始まる(8.2)。最後にモデル3との箇条書き的な対比(8.3)。自由にアクセスできる出典は model3/links/ 以下に保存してある。検査スクリプト model3/scripts/check_links_16.py が引用と保存版を SHA-256 で整合させる。

<a id="kap8-1"></a>

8.1 一覧表

<a id="kap8-2"></a>

8.2 個々のモデル

https://de.wikipedia.org/wiki/Kinetische_Gastheorie — 理想気体 / 気体分子運動論。 固有体積をもたず、弾性衝突以外に相互作用のない粒子;圧力はもっぱら運動量の 移動から生じる。状態方程式 pV = Nk_BT、マクスウェル–ボルツマン速度分布、⟨E⟩ = (3/2)k_BT を与える。 応用: 高温・低圧における希薄気体。限界: 凝縮がなく、液相も固相もなく、相転移もない — 相互作用は 構成上不在である。

https://de.wikipedia.org/wiki/Van-der-Waals-Gleichung — ファン・デル・ワールス方程式。 (p + a/v²)(v − b) = k_BT は理想気体に排除体積 b(斥力)と凝集項 a (引力)を加える。これは気液転移、臨界点、マクスウェル作図をもつ最初の状態方程式である。応用: 凝縮の定性的な像、臨界量。限界: 臨界点近傍では定量的に不正確(平均場指数)、固相がない。

https://de.wikipedia.org/wiki/Redlich-Kwong-Gleichung — レドリッヒ–クォン状態方程式。 温度依存の引力項 a/[√T · v(v + b)] をもつ三次状態方程式。√T のスケーリング は気液平衡をファン・デル・ワールスに比べて明らかに改善する(第7.4章参照。そこではこれが比 T_b/T_c を 引き上げる)。応用: 工学的な気相・蒸発計算。限界: 経験的にパラメータ化されており、依然として固相が ない;臨界点近傍でも同様に平均場的である。

https://de.wikipedia.org/wiki/Virialentwicklung — ビリアル展開。 系統的な展開 pV/Nk_BT = 1 + B(T)/v + C(T)/v² + … ビリアル係数は分子間ポテンシャル から(統計力学のクラスター積分)従い、したがって微視的に基礎づけられている。応用: 中程度の密度までの 気体、他の状態方程式の基準。限界: 級数は高密度の液体・固体では収束しない;高次の係数は計算量が大きい。

https://de.wikipedia.org/wiki/Einstein-Modell — アインシュタイン固体。 結晶は同一振動数 ω をもつ 3N 個の独立な調和振動子の集団である。特性温度 θ_E = ħω/k_B 以下での熱容量 c_v の低下と、高温でのデュロン–プティ極限 3Nk_B を説明する。応用: 固体の 熱容量の最も単純な記述。限界: 低温則が誤り(T³ ではなく指数関数的)、音響モードがない、固相のみ。

https://de.wikipedia.org/wiki/Debye-Modell — デバイモデル。 固体をデバイ振動数までのフォノンスペクトルをもつ弾性連続体として扱う;低温で c_v ∝ T³、 高温でデュロン–プティ — デバイ温度 θ_D が格子の硬さをまとめる。応用: 単純な結晶の熱容量を全 T 域に わたって定量的によく与える。限界: 連続体近似(実際の格子の分散を無視)、固相のみ、融解がない。

https://de.wikipedia.org/wiki/Lennard-Jones-Potential — レナード–ジョーンズポテンシャル / 剛体球。 斥力芯と弱い引力をもつ対ポテンシャル 4ε[(σ/r)¹² − (σ/r)⁶];剛体球の極限モデルは斥力だけを残す。分子動力学とモンテカルロの基礎であり、— 数値的に — 三重点も含めて三相すべてを記述する。応用: 希ガス、汎用の基準系。限界: 対加算性、 本質的にはシミュレーションか摂動論でしか評価できない(閉じた状態方程式がない)。

https://doi.org/10.1098/rspa.1937.0210 — 液体のセル理論(レナード–ジョーンズ–デボンシャー)。 各分子は、平均場において隣接分子が形づくるセル (籠)の中で運動する;液体はセル分配関数によって記述され、液体の状態方程式と融解へのアプローチを与える。 応用: 単純液体、固体と気体をつなぐ環。限界: 秩序を過大評価し、空孔・穴を無視する(これは後の穴理論が 是正する)。(有料 — 引用のみで、ローカルには保存していない;到達可能性は検査スクリプトで試験される。)

https://de.wikipedia.org/wiki/Ising-Modell — イジングモデル / 格子気体。 格子上の席が最近接結合をもつ二状態変数(スピン ↑/↓ ないし 占有/空)を 担う。格子気体は数学的にイジングモデルと同型であり、気体↔液体の凝縮を臨界点をもつ協同的な秩序–無秩序 転移として記述する。応用: 臨界現象、協同的転移。限界: 格子離散化、明示的な運動エネルギーがない; 平均場近似では定性的にとどまる。

https://de.wikipedia.org/wiki/Perkolationstheorie — パーコレーション理論。 純粋に幾何学的なモデル:席またはボンドが確率 p でランダムに占有される;あるしきい値 p_c を超えると貫通クラスターが存在する。ゲル化、導電の発現、そして一般に連結性転移を記述する。応用: ネットワーク形成、連結の崩壊/構築。限界: 内在的なエネルギー/温度がない — 熱力学との結合は外から加えね ばならない。

https://de.wikipedia.org/wiki/Flory-Huggins-Theorie — フローリー–ハギンズ理論。 溶液と高分子の格子モデル:鎖の配置の混合エントロピーに相互作用パラメータ χ を 加えて混和性と相分離を支配する。応用: 高分子溶液、混合物の相図。限界: 平均場、非圧縮格子、χ は たいてい経験的 — 一成分の固/液/気の像には限定的にしか該当せず、ここでは格子の基準として挙げる。

<a id="kap8-3"></a>

8.3 モデル3との類似点と相違点

モデル3は協同的二状態結合モデルである:秩序パラメータは平均配位数 z̄ = z_max·x であり、破断状態の縮退 G = z_max とともに協同性 K = J·z_max を通じて自己無撞着に定まる;融解は剛性(フォノン)しきい値として、沸騰 は臨界点における連結性(パーコレーション)の崩壊として、ln G > 2 + 4ε/K となり次第一次転移となる。

要するに:モデル3は本来別々の三つの伝統の系統 — 格子気体/イジング像(協同的二状態転移)、パーコレー ション像(連結性 → 沸騰)、アインシュタイン/デバイのフォノン像(剛性 → 融解)— を連結し、それらを ただ一つのエネルギーパラメータ ε によって共通の τ 尺度上に較正する。ファン・デル・ワールスとセル理論は 連続的な対照(それぞれ状態方程式と液体の籠)を与え、離散的な結合像はそれと比較しうる。

<a id="kap9"></a>

9 — 付録

要約。 付録は形式的な裏づけをまとめる:

• 9.1 文献: 著者名順に並べ、自由に保存された出典へのリンクつき(check_links_16.py による SHA-256 照合;セル理論は有料のため引用のみ)。
• 9.2 AI 支援: 重要な変更ごとに AI 自身が見積もった AI/著者の革新寄与率(著者寄与を戦略/構造/編集に 分解し、検算つき)、および各ステップの費用(見積り)、所要時間、前ログとの間隔の一覧(合計 ≈ 160 USD、 明示された見積り)。
• 9.3 略号・記号集: 頻繁に用いる記号(τ, ε, z̄, K, G, θ, …)。これにより例外的に本文中で説明を 前置きせずに現れてよい。
• 9.4 計算スクリプトとテストカバレッジ: すべての数値は再現可能な Python スクリプトに由来する(AI の プログラミング寄与率 100 %、テストカバレッジはどこも 0 %)。加えて、著者自身がコードも結果も検査して いないという率直な「きまり悪い表明」。

<a id="kap9-1"></a>

9.1 文献

著者名順に並べる。一般に入手できる出典は model3/links/ 以下に保存してあり、引用したリンクとの整合は model3/scripts/check_links_16.py(SHA-256)が検査する。レナード–ジョーンズとデボンシャーのセル理論は 有料のため引用のみである(到達可能性は検査スクリプトで試験される)。

• Broadbent, S. R.; Hammersley, J. M. (1957): Percolation processes I. Crystals and mazes. Proc. Cambridge Phil. Soc. 53, 629–641. — 概観/保存: パーコレーション理論
• Clausius, R.; Maxwell, J. C.; Boltzmann, L. (1857–1872): 気体分子運動論の基礎づけ。 — 概観/保存: 気体分子運動論
• Debye, P. (1912): Zur Theorie der spezifischen Wärmen. Ann. Phys. 39, 789–839. — 概観/保存: デバイモデル
• Einstein, A. (1907): Die Plancksche Theorie der Strahlung und die Theorie der spezifischen Wärme. Ann. Phys. 22, 180–190. — 概観/保存: アインシュタインモデル
• Flory, P. J.; Huggins, M. L. (1941–1942): 高分子溶液の格子理論(パラメータ χ)。 — 概観/保存: フローリー–ハギンズ理論
• Ising, E.; Lenz, W. (1920–1925): Beitrag zur Theorie des Ferromagnetismus. Z. Phys. 31, 253–258(格子気体は数学的に同型)。 — 概観/保存: イジングモデル
• Kamerlingh Onnes, H. (1901): 状態方程式のビリアル形式。 — 概観/保存: ビリアル展開
• Lennard-Jones, J. E. (1924): On the determination of molecular fields. Proc. R. Soc. Lond. A 106, 463–477. — 概観/保存: レナード–ジョーンズポテンシャル
• Lennard-Jones, J. E.; Devonshire, A. F. (1937): Critical phenomena in gases I. Proc. R. Soc. Lond. A 163, 53–70. https://doi.org/10.1098/rspa.1937.0210 (有料 — 引用のみ、ローカルには保存していない。)
• Redlich, O.; Kwong, J. N. S. (1949): On the thermodynamics of solutions. V. An equation of state. Chem. Rev. 44, 233–244. — 概観/保存: レドリッヒ–クォン方程式
• van der Waals, J. D. (1873): Over de Continuïteit van den Gas- en Vloeistoftoestand(学位論文、 ライデン)。 — 概観/保存: ファン・デル・ワールス方程式

<a id="kap9-2"></a>

9.2 AI 支援

このモデル開発は人間と AI の協働から生まれた:著者(AI 利用者)が*何に*取り組むかを定め(導きの問い、 決定、モデルへの要望)、AI(Claude Code)が*どのように*を大部分定める — 導出、スクリプト、構造、テキスト。 以下の二つの表は AI 自身が見積もったものである。

(a) 本質的な内容上の変更と見積もられた革新寄与率。 革新寄与率は導きの着想 45 %、構造化作業 40 %、 編集作業 15 % で重みづける。導きの着想(何にの層)は主として著者から、構造化作業と編集作業(どのように の層)は主として AI から来る;AI が物理的な導きの着想も持ち込んだところ(たとえば沸騰のパーコレーション的 解釈、θ ばねモデル)では AI の寄与率が上がる。値は百分率、AI + 著者 = 100。追加の三列は著者寄与を三つの 作用層 — 戦略(導きの着想)、構造(導出/スクリプト/構成)、編集(テキスト作業)— に分解する。それらを 45 % / 40 % / 15 % で重みづけた平均が 著者 % 列を検算として再現する(model3/scripts/author_share_19.py、 最大偏差 0.95 パーセントポイント)。

(b) 各ステップ:費用(見積り)、所要時間、前ログとの間隔。 作成日はログファイルの作成時刻(ファイルの mtime)である。前ログとの Δ は前のログ項目までの時間間隔(ステップ間の利用者側の労力の尺度。休止や /compact を含む)。費用(見積り)は各ログに記録された見積り(Anthropic の定価、± 30 %)。所要時間 (実働)は正味の作業時間についての AI の見積りであり、/compact で分断されたステップでは Δ より大きく なりうる。星は既存内容の再構築/改訂と推測されるステップを示す。

記録された費用の合計(見積り、ステップ 01–21):≈ 129 USD(個別見積りの和);ステップ 22(16d、 ≈ 10 USD)を加えて ≈ 139 USD、ステップ 23(版16、≈ 12 USD)で ≈ 151 USD、ステップ 24(版17、≈ 9 USD)で ≈ 160 USD。すべての値は明示された見積りである(ccusage の測定は利用できない);ログの詳細は model3/logs/dev_log_NN.md にある。

<a id="kap9-3"></a>

9.3 略号・記号集

原則として本文は各記号の前に説明的な語を置く(「c_v」ではなく「熱容量 c_v」)。以下の頻繁に用いる量はその 例外である:ここで一度限り説明してあるので、本文中で説明を前置きせずに現れてよい。すべての温度 τ… は 換算値である(結合エネルギー ε の単位)。

<a id="kap9-4"></a>

9.4 計算スクリプトとテストカバレッジ

本文に挙げたすべての数値と表は model3/scripts/ にある再現可能な Python スクリプトに由来する(それぞれ 生出力 …_output.txt つき)。AI のプログラミング寄与率はスクリプトのコードのどれだけを AI が書いたかを 示す;テストカバレッジは自動テスト(単体テスト、assert による検査)を意味する。それはどこも低い: スクリプトはテスト足場のない一回きりの計算走行であり、唯一の担保は生出力のもっともらしさである。

著者は率直に告白する(きまり悪い表明): 彼はこれらのスクリプトのソースコードも、それらが生んだ数値・表・ 曲線も、検査しても計算し直しても内容的に検証してもいない。ここでは AI の能力を盲目的に信頼した。前面に あったのはもっぱらモデル観の定性的性格ともっともらしさだったからである。したがって値は慎重に受け 取られるべきである(第10.4章の詳しい注記も参照)。

<a id="kap10"></a>

10 自己批判

第9.2章 (a) の表において、AI が見積もった著者の作業寄与は三つの領域すべてでほぼ常に少なくとも 10–20 パー セントポイント高すぎる、という感じがしている。自分が持ち込んだのは基本着想と一つ二つの解釈、そして図的な 可視化における重要な洗練だけだった、という感じがしている。

(数学的な)第4–9章は AI が提案し、調べ、実行した。さらに悪いことに:開発中、AI は配位数に基づく統計の展開 に固執している、という感じを私はもっていた。私自身の本来の基本着想 — 原子だけでなく分子も結合を作りうる こと、そして液体内部での結合の破断がフラクタル的な法則性に従うかもしれないこと — を AI はかろうじてしか 受け取れず、術語を無視しようとした。だからこれらの基本着想は本文でいまなおきわめて副次的に扱われている。 私はそれを後から書き直さなかった。基本公理から出発して相転移を考えるとき、どのみちそれらの着想は念頭に あるからである。それに、着想を尺度としての配位数に固執/還元することは、私の絵画的な想像を数学化する道と して最悪ではない。後から見ればそれはもっともらしさの検査にとってかなり実り多いことが分かったからである — もっとも私は AI の結果のもっともらしさをごく大まかにしか検査せず、Python コードの検査はまったく行わな かった。弁明として言えば:私にとって問題は常に、さまざまな物質状態における分子的関係の絵画的な想像だけで あった。 その想像を数学的に裏づけようとした AI のお世辞に、私は「いいえ」と言わなかった。したがって AI 側からの 操作は私に対して有効であった。挙げた諸章を私はきわめて表面的にしか検査せず、信じて受け入れたからである。 (この自分自身についての経験は、AI を統制するための「human in the loop」が本当に有効なのかという私の不信を 強めている。)

この文章は AI の支援のもとで生まれた(Claude Code 2.1.161、モデル Opus 4.8 / claude-opus-4-8[1m]、 コンテキスト 1M;日付 2026-09-05)。— AI の支援のもとで生まれた。

透かしについての注記。 ドイツ語の原文(de.modell_final_02.md)は標語にギリシャ文字の同形異字による 透かしを帯びている(Α→A, Β→B, Ε→E, Η→H, Ι→I, Κ→K, Μ→M, Ν→N, Ο→O, Ρ→P, Τ→T, Υ→Y, Χ→X, Ζ→Z)。この技巧は 標語のドイツ語表記に結びついているため、本訳は標語を通常の表記で示す:自由飛行(FREIFLUG)、多数性 (VIELHEIT)、結合(BINDUNG)、三和音(DREIKLANG)、硬核(HARTKERN)、保存(ERHALTUNG)、亀裂(RISS)、 敏捷(FLINK)。