第 3 章 蒸気機関の効率を上げたい — 熱力学とエントロピーの誕生¶
前回までのあらすじ: 第 1 章で Newton が重力のモデルを作り、第 2 章で Maxwell が電気と磁気を統一した。どちらも「好奇心」が主な動機だった。この章では、動機が全く異なる——実用的な必要性から生まれたモデルの物語を見る。そして、蒸気機関の効率という実用的な問いが、宇宙の根本的な性質——エントロピー——に到達する。
この章のゴール
- 「必要性」から始まった探求が、宇宙の根本的な性質に到達する物語を追う
- Carnot の効率限界と Boltzmann のエントロピーを理解し、「ミクロな世界の状態の数え上げ」がマクロな熱現象を説明するという統計力学の核心を把握する
3.1 動機:蒸気機関の効率をもっと上げられないか?¶
🟡 リナ: 第 1 章と第 2 章では、「惑星はなぜ動くのか」「電気と磁気は統一できないか」という好奇心がモデルを生んだ話をしたわね。今日は全く違う動機——お金よ。
🔵 カイ: お金?
🟡 リナ: 18 世紀のイギリス。産業革命の真っ只中。石炭を燃やして蒸気を作り、蒸気の力で機械を動かす——蒸気機関が経済の心臓だった。でも、蒸気機関の効率は低かった。投入した熱エネルギーのうち、実際に仕事に変換できるのはほんの一部。残りは排熱として捨てられる。
🔵 カイ: それってもったいないですね。どのくらい無駄になってたんですか?
🟡 リナ: 初期の蒸気機関は効率がわずか数%——投入した熱の 90%以上が無駄になっていたの。だから「もっと効率を上げられないか? 理論的な限界はあるのか?」という問いが切実だった。これに答えたのが、フランスの Sadi Carnot (サディ・カルノー)。1824 年の論文『火の動力についての考察』。
🔭 科学哲学メモ: 物理学のモデルは「好奇心」だけでなく「実用的必要性」からも生まれる。しかし、どちらの動機から始まっても、深い場所に到達できる。動機の種類はモデルの価値を決めない。
✅ 理解度チェック: 蒸気機関の効率の理論的限界を最初に示したのは誰でしょうか?
答え
Sadi Carnot(サディ・カルノー)。1824 年の論文『火の動力についての考察』で示した。
✅ 理解度チェック: この章のモデルが生まれた動機は、第 1 章・第 2 章とどう異なるでしょうか?
3.2 Carnot の問い — 効率に限界はあるか?¶
🟡 リナ: Carnot の問いはシンプルだった。熱を仕事に変換する効率に、理論的な上限はあるか?
🔵 カイ: 上限があるんですか? 技術が進歩すれば、いくらでも効率を上げられそうですけど。
🟡 リナ: それが、上限があるの。どんなに技術が進歩しても、原理的に超えられない壁がある。Carnot は実験ではなく、思考実験でこれを導いた。まず、熱機関が何をしているかを整理しましょう。図 3.1「熱機関のエネルギーフロー」 を見て。図では高温源の温度を \(T_\text{hot}\)、低温源を \(T_\text{cold}\) と書いてあるけど、これは「高温 (hot)」「低温 (cold)」を明示するためで、これから先の本文では \(T_H\)(\(H\) = hot)、\(T_C\)(\(C\) = cold)と略記するわね。熱も同様に \(Q_H\)(高温源から受け取る熱)、\(Q_C\)(低温源に排出する熱)と書くわ。
%%{init: {"theme": "default", "themeCSS": ".edgePath .path, .flowchart-link { stroke-width: 2px !important; }"}}%%
flowchart LR
Hot["高温源<br/>温度 T_hot"] -->|"熱 Q_hot"| Engine["熱機関"]
Engine -->|"仕事 W"| Work["仕事<br/>(有用な出力)"]
Engine -->|"排熱 Q_cold"| Cold["低温源<br/>温度 T_cold"]
style Hot fill:#f66,color:#fff
style Cold fill:#66f,color:#fff
style Engine fill:#ff9
図 3.1: 熱機関のエネルギーフロー
熱機関は高温源(温度 \(T_H = T_\text{hot}\))から熱 \(Q_H = Q_\text{hot}\) を受け取り、一部を仕事 \(W\) に変換し、残り \(Q_C = Q_\text{cold}\) を低温源(温度 \(T_C = T_\text{cold}\))に排出する。
Carnot サイクルの 4 つのステップ¶
🟡 リナ: Carnot は「もし完璧に無駄のない理想的な熱機関があったとしたら」と考えた。全ての過程が可逆(元に戻せる)な理想的なサイクルを構成したの。4 つのステップからなるわ。
- 等温膨張(温度 \(T_H\) で):「等温」は温度一定という意味。高温源に接触したまま気体を膨張させる。温度が変わらないよう熱源から熱 \(Q_H\) を吸収しながら、仕事をする。
- 断熱膨張:「断熱」は熱の出入りがないという意味。高温源から切り離し、さらに膨張。熱の出入りなしに温度が \(T_H\) から \(T_C\) に下がる。
- 等温圧縮(温度 \(T_C\) で):低温源に接触したまま気体を圧縮する。温度が変わらないよう気体は熱 \(Q_C\) を低温源に放出する。
- 断熱圧縮:低温源から切り離し、さらに圧縮。熱の出入りなしに温度が \(T_C\) から \(T_H\) に戻る。
⚪ メイ: 温度の上げ下げを「断熱」で、熱のやり取りを「等温」で——役割がきれいに分かれているのね。
🟡 リナ: そう。各ステップで何が起きているかを物理的に図にしたから見て(図 3.2「Carnotサイクルの4ステップ」)。
図 3.2: Carnotサイクルの4ステップ。シリンダーとピストンで各ステップを物理的に図示。① 高温源から熱 \(Q_H\) を吸収して膨張、② 断熱膨張で温度低下、③ 低温源へ熱 \(Q_C\) を放出して圧縮、④ 断熱圧縮で温度回復。
🟡 リナ: この 4 つのステップを圧力-体積図(PV 図)に描くと 図 3.3「CarnotサイクルのPV図」 のようになるわ。
図 3.3: CarnotサイクルのPV図。Carnot サイクルの 4 つのステップ(等温膨張→断熱膨張→等温圧縮→断熱圧縮)を圧力-体積図上に示す。囲まれた面積が正味の仕事 \(W\) に対応する。
🟡 リナ: 4 つのステップを経て、気体は元の状態に戻る——だから「サイクル」と呼ぶの。
🔵 カイ: 元に戻るってことは……内部エネルギーも元に戻る?
🟡 リナ: その通り。状態が元に戻るなら、内部エネルギーの変化はゼロ。エネルギー保存から:
つまり、仕事は「吸収した熱」から「排出した熱」を引いたもの。
🔵 カイ: なるほど、差し引きで残った分だけが仕事になるんですね。でも、\(Q_C\) をゼロにできれば全部仕事にできるんじゃ……?
効率の定義と Carnot の結論¶
🟡 リナ: 効率 \(\eta\) は「投入した熱のうち、どれだけが仕事に変わったか」で定義される:
🔵 カイ: 効率 100%にするには \(Q_C = 0\)、つまり排熱ゼロにすればいい? でも、排熱をゼロにできない物理的な理由があるんですか?
🟡 リナ: まさにそこが核心。Carnot は、可逆サイクルでは \(Q_H/T_H = Q_C/T_C\) が成り立つことを示した——つまり \(T_C > 0\) なら必ず \(Q_C > 0\) で、排熱をゼロにはできないの。
🔵 カイ: えっ、なんで \(Q_H/T_H = Q_C/T_C\) が成り立つんですか? 温度で割ると等しくなるって、直感的にはピンと来ないんですけど。
🟡 リナ: いい疑問。今から具体的に導出するわ。
Carnot 効率の導出¶
🟡 リナ: 理想気体を作業物質として使うと、各ステップを具体的に計算できる。理想気体の状態方程式は \(pV = Nk_B T\) よ。高校では \(pV = nRT\)(\(n\) は物質量、\(R\) は気体定数)と習ったかもしれないけど、\(N\) 個の分子で書き直すと \(Nk_B = nR\) だから同じ式。ここで \(k_B \approx 1.381 \times 10^{-23}\;\text{J/K}\) は Boltzmann 定数——1 個の分子レベルでエネルギーと温度を結びつける換算係数よ。高校で習った気体定数 \(R\) とは \(R = N_A k_B\)(\(N_A\) はアボガドロ数)の関係にあるの。この章では分子の数 \(N\) で書く方が、後で統計力学に繋げやすいの。
ステップ 1(等温膨張 \(A \to B\)、温度 \(T_H\)):
まず、内部エネルギーとは、系を構成する全ての粒子が持つエネルギーの合計のこと——粒子の運動エネルギーと、粒子間の相互作用のポテンシャルエネルギーを全部足し合わせたものよ。理想気体では粒子間に力が働かないから、ポテンシャルエネルギーはゼロで、内部エネルギーは全て運動エネルギー——つまり粒子の速さだけで決まる。等温過程では温度が一定だから、粒子の平均的な速さも一定——したがって内部エネルギーも変わらない。
🔵 カイ: なるほど、温度が変わらなければ粒子の動きも変わらないから、内部エネルギーも一定なんですね。
🟡 リナ: その通り。次に、エネルギー保存を使うわ。系に加えた熱は、内部エネルギーの増加と系がした仕事に分配される——つまり「加えた熱 = 内部エネルギーの変化 + した仕事」。式で書くと:
ここで \(dU\) は内部エネルギーの微小変化、\(\delta Q\) は加えた微小な熱量、\(p\,dV\) は体積が \(dV\) だけ変化したときに系がする微小な仕事よ。\(\delta\) と \(d\) で記号を使い分けているのが気になるかもしれないけど、今は「どちらも微小な量を表す記号」と思っておいて。違いの理由は 「第 1 法則(エネルギー保存)」 で説明するわ。これはエネルギー保存そのもので、後のセクションで第 1 法則として正式に名前をつけるわ。
今は「加えた熱=内部エネルギー変化+した仕事」と覚えておいて。等温過程では内部エネルギー変化がゼロだから \(\delta Q = p\,dV\)。つまり吸収する熱は \(\delta Q = p\,dV\) を全過程にわたって足し上げた(積分した)もの。状態方程式 \(pV = Nk_BT_H\) から \(p = Nk_BT_H/V\) を代入すると:
ここで \(\int_{V_A}^{V_B} \frac{dV}{V} = [\ln V]_{V_A}^{V_B} = \ln V_B - \ln V_A = \ln\frac{V_B}{V_A}\) を使ったわ。\(\int 1/x\,dx = \ln x\) は高校で習ったわね。
⚪ メイ: 吸収した熱が温度と体積比の対数で決まるのね。体積が大きく膨張するほど、たくさん熱を吸い込むと。
ステップ 3(等温圧縮 \(C \to D\)、温度 \(T_C\)):
🟡 リナ: 等温過程だから、ステップ 1 と同じく内部エネルギー変化がゼロで \(\delta Q = p\,dV\) が成り立つ。ただし今度は温度が \(T_C\) だから、状態方程式は \(p = Nk_BT_C/V\) よ。ステップ 1 では膨張(\(dV > 0\))だったから \(\delta Q > 0\)——系が熱を吸収していたわね。今度は圧縮なので \(dV < 0\)——つまり \(p\,dV < 0\) で、\(\delta Q = p\,dV < 0\) になる。\(\delta Q\) は「系に加えた熱」だから、これが負ということは「系から熱が出ていった」——つまり系が熱を放出しているということよ。過程全体で系に加えた熱の合計は、体積が \(V_C\) から \(V_D\) に変化する間の積分で求まる(圧縮だから \(V_D < V_C\)):
\(\int_{V_C}^{V_D} \frac{dV}{V} = [\ln V]_{V_C}^{V_D} = \ln V_D - \ln V_C = \ln(V_D/V_C)\) と機械的に計算すれば正しい結果が出るわ。
ここで \(V_D < V_C\)(圧縮だから体積が減る)なので \(\ln(V_D/V_C) < 0\)——つまり積分の結果 \(Nk_B T_C \ln(V_D/V_C)\) は負の値になる。確かに系は熱を放出している。
ここで \(Q_C\) を「系が低温源に放出した熱の大きさ」として正の値で定義したいの。積分結果が負なのは「系が熱を失った」ことを表しているから、その絶対値を取ればよい。つまり:
最後の等号は \(-\ln(a/b) = \ln(b/a)\) を使ったわ。
🔵 カイ: 2 つの断熱過程はどう使うんですか?
🟡 リナ: 断熱過程では \(\delta Q = 0\) だから、第 1 法則は \(dU = -p\,dV\) になる。単原子理想気体(ヘリウムのような 1 原子の気体)では、粒子は \(x, y, z\) の 3 方向に運動するだけで、各方向あたり \(\frac{1}{2}k_BT\) のエネルギーを持つことが実験的に知られているの(\(k_B\) は先ほど状態方程式で出てきた Boltzmann 定数よ)。つまり 1 粒子あたり \(\frac{3}{2}k_BT\)、\(N\) 粒子で \(U = \frac{3}{2}Nk_BT\)。これがなぜ成り立つかは、この章の後半(3.7「温度の統計力学的な意味」)で統計力学から導出するわ。今は実験事実として使わせてね。微小変化は \(dU = \frac{3}{2} N k_B\,dT\)。これと状態方程式 \(p = Nk_BT/V\) を組み合わせて、断熱過程で温度と体積がどう関係するか導出してみましょう。\(dU = -p\,dV\) に代入すると:
両辺を \(Nk_B T\) で割ると \(\frac{3}{2}\frac{dT}{T} = -\frac{dV}{V}\)。これは「左辺は \(T\) だけの式、右辺は \(V\) だけの式」という形になっているわね。こういうとき、両辺をそれぞれの変数で積分できるの(変数分離法と呼ばれる手法よ)。\(\int dT/T = \ln T\) と \(\int dV/V = \ln V\) を使うと、\(\frac{3}{2}\ln T = -\ln V + \text{const}\)。左辺に \(\ln V\) を移項して \(\frac{3}{2}\ln T + \ln V = \text{const}\)。対数の性質 \(a\ln x = \ln x^a\) と \(\ln x + \ln y = \ln(xy)\) を使うと \(\ln(T^{3/2} V) = \text{const}\)。両辺の指数を取れば \(T^{3/2} V = \text{一定}\)。
🔵 カイ: おお、温度と体積がセットで決まるんですね。膨張すれば温度が下がる。
🟡 リナ: そう。このままでも使えるけど、\(T\) の指数を 1 にして「\(T\) と \(V\) の関係」を見やすくしたいから、両辺を \(2/3\) 乗するわ。\((T^{3/2} V)^{2/3} = T^{(3/2)(2/3)} \cdot V^{2/3} = T \cdot V^{2/3}\) だから、\(T V^{2/3} = \text{一定}\) と書ける(右辺の定数の具体的な値は変わるけど「一定」であることは同じ)。高校で \(pV^\gamma = \text{一定}\)(\(\gamma\) は比熱比)と習ったかもしれないけど、単原子理想気体では \(\gamma = 5/3\) で、\(TV^{\gamma-1} = TV^{2/3} = \text{一定}\) と同じ式になるの。大事なのは「断熱過程では温度と体積が独立に変化できない」ということ。これを断熱膨張 \(B \to C\) と断熱圧縮 \(D \to A\) に適用すると:
1 つ目の式を 2 つ目の式で割ってみましょう。左辺は \(\frac{T_H V_B^{2/3}}{T_C V_D^{2/3}}\)、右辺は \(\frac{T_C V_C^{2/3}}{T_H V_A^{2/3}}\)——あ、ちょっと複雑ね。もっと簡単にやりましょう。1 つ目の式から \(V_B^{2/3} = (T_C/T_H) V_C^{2/3}\)、2 つ目の式から \(V_A^{2/3} = (T_C/T_H) V_D^{2/3}\)。この 2 つを割ると \(T_C/T_H\) がキャンセルして:
\(2/3\) 乗を外す(両辺を \(3/2\) 乗する)と:
🔵 カイ: あ、対数の中身が同じになる! ということは \(Q_C/Q_H\) の比を取ったら……対数が消えるんじゃないですか?
🟡 リナ: その通り。\(\ln(V_B/V_A) = \ln(V_C/V_D)\) だから、\(Q_H\) と \(Q_C\) の比を取ると対数因子がキャンセルして温度の比だけが残る:
したがって、Carnot サイクルの効率は:
🔵 カイ: これって単原子理想気体で計算したから成り立つ結果ですよね? 他の気体だと変わりませんか?
🟡 リナ: いい疑問。実は Carnot の効率は作業物質に依存しないの。ここでは単原子理想気体で具体的に計算したけど、2 原子分子でも液体でも、可逆サイクルなら同じ結果になる。これは Carnot の定理の帰結——もし作業物質によって効率が変わるなら、効率の高い方と低い方を組み合わせて第 2 法則に矛盾する装置が作れてしまうから。
⚪ メイ: つまり、低温源の温度がゼロでない限り、効率は 100%にならない。
🟡 リナ: さらに重要なのは、Carnot の定理——同じ温度の間で動作する全ての熱機関の中で、可逆機関(Carnot 機関)の効率が最大だということ。不可逆な機関は必ずこれより低い効率しか出せない。
🔵 カイ: なんで?
🟡 リナ: もし Carnot 機関より効率の高い機関があったとすると、その機関と Carnot 機関の逆回し(ヒートポンプ)を組み合わせることで、「低温源から高温源へ、何の代償もなく熱を移す」ことが可能になる。これは後で述べる第 2 法則に矛盾するの。
📝 練習問題:
- Carnot 効率の具体的な計算 → 問題 B-1. Carnot 効率の計算
✅ 理解度チェック: Carnot サイクルの効率の上限を表す式を書いてください。
答え
\(\eta_{\text{Carnot}} = 1 - \frac{T_C}{T_H}\)
✅ 理解度チェック: Carnot の効率限界が 100%にならないのはどのような場合でしょうか?
答え
低温源の温度 \(T_C\) がゼロでない限り、効率は 100%にならない。
3.3 熱力学の法則 — エネルギー保存と方向性¶
🟡 リナ: Carnot の仕事を受けて、19 世紀に熱力学の法則が整備された。全部で第 0 法則から第 3 法則まであるけど、特に重要なのは第 0 法則、第 1 法則、第 2 法則。
第 0 法則(熱平衡の推移律)¶
🟡 リナ: 第 0 法則は温度の存在を保証する法則よ。
系 \(A\) と系 \(C\) が熱平衡にあり、系 \(B\) と系 \(C\) も熱平衡にあるなら、系 \(A\) と系 \(B\) も互いに熱平衡にある。
⚪ メイ: 当たり前に聞こえますけど……。
🟡 リナ: 「当たり前」に見えるけど、この推移律が成り立つからこそ、全ての系を「温度」という一つの数値でラベル付けできるの。もしこの法則が成り立たなければ、「\(A\) と \(C\) は平衡だけど \(A\) と \(B\) は平衡でない」なんてことが起きて、温度という概念自体が意味をなさなくなる。
✅ 理解度チェック: 熱力学の第 0 法則はなぜ重要でしょうか? この法則が成り立たなければ何が困るでしょうか?
答え
第 0 法則(熱平衡の推移律)が成り立つからこそ、全ての系を「温度」という一つの数値でラベル付けできる。成り立たなければ温度という概念自体が意味をなさなくなる。
第 1 法則(エネルギー保存)¶
🟡 リナ: 第 1 法則は、エネルギー保存。微小変化の形で書くと:
ここで \(dU\) は内部エネルギーの変化、\(\delta Q\) は系に加えられた熱、\(\delta W\) は系がした仕事。気体の膨張・圧縮を考えるなら \(\delta W = p\,dV\) だから:
これは Carnot サイクルの導出で使った「加えた熱=内部エネルギー変化+した仕事」(\(\delta Q = dU + p\,dV\))と全く同じ式を、\(dU\) について解いた形よ。
🔵 カイ: \(d\) と \(\delta\) が違うのはなぜですか?
🟡 リナ: いい質問。\(dU\) は状態量の微小変化——系の状態が決まれば \(U\) の値は一意に決まる。一方、\(\delta Q\) と \(\delta W\) は経路に依存する量——どういう過程を経たかによって値が変わる。銀行口座の残高(状態量)と、現金入金額(経路に依存)の違いと同じよ。
⚪ メイ: 整理するとこうなるわね。
表 3.1: 状態量と経路依存量の比較
| 分類 | 記号 | 例 | 特徴 |
|---|---|---|---|
| 状態量(\(d\)) | \(dU\), \(dS\), \(dV\) | 内部エネルギー、エントロピー、体積 | 状態が決まれば値が一意に決まる。経路に依存しない |
| 経路依存量(\(\delta\)) | \(\delta Q\), \(\delta W\) | 熱、仕事 | どういう過程を経たかで値が変わる |
✅ 理解度チェック: 第 1 法則の式 \(dU = \delta Q - \delta W\) で、\(d\) と \(\delta\) の記号が使い分けられている理由は何でしょうか?
答え
\(dU\) は状態量(系の状態が決まれば値が一意に決まる量)の微小変化を表し、\(\delta Q\) と \(\delta W\) は経路に依存する量(どういう過程を経たかで値が変わる量)を表すため。
第 2 法則(エントロピー増大・不可逆性)¶
🟡 リナ: 第 2 法則は「方向性」の法則。いくつかの同値な表現があるけど:
Clausius の表現: 熱は、自発的には低温の物体から高温の物体へ流れない。
Kelvin の表現: ある過程の唯一の結果が、熱源から熱を取り出してそれを完全に仕事に変換することであるような過程は存在しない。
🔵 カイ: それは当たり前じゃないですか? 熱いコーヒーは冷めるけど、冷めたコーヒーが勝手に熱くなることはない。
🟡 リナ: 「当たり前」に見えるでしょう? でも、実はこれは深い謎なの。Newton 力学の方程式 \(F = ma\) は時間反転対称——つまり、時間を逆回しにしても方程式の形が変わらない。ミクロな粒子の運動は、前向きでも逆向きでも物理法則に矛盾しない。
🔵 カイ: え、ちょっと待ってください。ミクロな法則が時間反転対称なら……「熱が低温から高温に流れる」過程もミクロでは許されてるってことですか? でもそれって矛盾しません? ミクロでは OK なのに、マクロでは絶対に起きないなんて。じゃあ第 2 法則は Newton の法則みたいな「絶対の法則」じゃないってことですか?
🟡 リナ: いい直感ね。実はその通り——第 2 法則の性格は Newton の運動方程式とは本質的に異なるの。この矛盾を解いたのが Boltzmann——でもその前に、エントロピーの熱力学的な定義を見ておきましょう。
✅ 理解度チェック: 熱力学の第 1 法則を微小変化の形で表してみましょう。
答え
\(dU = \delta Q - p\,dV\)(系に加えた熱はエネルギー保存により、内部エネルギーの変化と系がした仕事の和に等しい)。
✅ 理解度チェック: 熱力学の第 2 法則の Clausius による表現は何でしょうか?
答え
熱は、自発的には低温の物体から高温の物体へ流れない。
3.4 エントロピーの熱力学的定義¶
🟡 リナ: Carnot サイクルの結果から、驚くべきことが導ける。可逆サイクルでは:
ここで少し記法を整理させて。さっきまで \(Q_H\) も \(Q_C\) も「熱の大きさ」として正の値で扱ってきたわね。でもこれから一般化するために、系が受け取る熱を正、系が放出する熱を負という符号の約束に切り替えるわ。この約束だと、等温圧縮で系が放出する熱は \(-Q_C\) と書ける。すると上の式は:
🔵 カイ: 一周して足し合わせたらゼロ……これって偶然ですか?
🟡 リナ: 偶然じゃないの。実はこれを一般化できる。Carnot サイクルは特殊な形の閉曲線だけど、任意の形の可逆サイクルでも同じことが成り立つことを示したいの。PV 図上で任意の閉曲線を描いたとき、その内部を細かい格子状に区切ることを考えて。なぜ格子が作れるかというと、PV 図の各点には「その点を通る温度一定の曲線(等温線)」と「その点を通る熱の出入りがない曲線(断熱線)」がそれぞれ 1 本ずつ引けるの。等温線は「\(T\) を固定して状態方程式 \(pV = Nk_BT\) を満たす \((p, V)\) の集まり」だから各温度に 1 本ずつ存在する。断熱線も同様に「\(TV^{2/3} = \text{一定}\)」を満たす曲線として各点を 1 本ずつ通る。ちょうど地図上のどの点にも「経線」と「緯線」が 1 本ずつ通るのと同じ。
🔵 カイ: あ、地図の経線・緯線みたいに、PV 図のどこでも等温線と断熱線が交差して格子を作れるんですね。
🟡 リナ: その通り。等温線と断熱線を格子線として使えば、PV 図を細かいマス目に分割できる。各マス目は「等温→断熱→等温→断熱」の 4 辺を持つ——つまり微小な Carnot サイクルになっているの。先ほど理想気体で示した \(Q_H/T_H = Q_C/T_C\) は、実は Carnot の定理により作業物質に依存しない——どんな物質を使っても可逆 Carnot サイクルでは同じ関係が成り立つの(もし成り立たない物質があれば、それと理想気体を組み合わせて第 2 法則に矛盾する装置が作れてしまう)。だから各微小 Carnot サイクルでも \(\delta Q_H/T_H + \delta Q_C/T_C = 0\) が成り立つ。任意の閉曲線の内部をこのマス目で埋め尽くせば、曲線を微小な Carnot サイクルの集まりとして近似できる——ちょうど曲線を階段状に近似するようなものよ。
🔵 カイ: 隣り合うサイクルの共有辺で打ち消し合うっていうのが、ちょっとイメージしにくいんですけど……。
🟡 リナ: こう考えて。タイルを敷き詰めるイメージよ。床全体を小さな正方形のタイルで埋めたとき、隣り合うタイルの境界線は 2 枚のタイルに共有されているわね。一方のタイルから見ると「右の辺」だけど、隣のタイルから見ると「左の辺」——同じ辺を逆向きに通ることになる。
物理的に言うと、ある微小サイクルの等温膨張の辺は、隣の微小サイクルから見ると等温圧縮の辺になっている。同じ等温線上の同じ区間を、一方は膨張方向に(\(\delta Q > 0\))、他方は圧縮方向に(\(\delta Q < 0\))通る。温度が同じで変化量の大きさも同じだから、\(\delta Q/T\) の寄与が正確に打ち消し合うの。
🔵 カイ: あ、つまり隣のサイクルの辺を足すと \(+\delta Q/T + (-\delta Q/T) = 0\) になるってことですね。
🟡 リナ: その通り。そして断熱線の辺はどうかというと——断熱とは「熱の出入りがない」という意味だったわね(ステップ 2 と 4 で確認したでしょう?)。だから断熱線上では \(\delta Q = 0\) で、\(\delta Q/T = 0\) だから寄与しない。結局、内部の辺は全て「等温線上の辺が隣同士で相殺」か「断熱線上の辺でそもそもゼロ」のどちらかで消えて、打ち消されずに残るのは外周の辺だけ。
🔵 カイ: あ、タイルの比喩で分かりました! 内部は全部相殺して、外周だけ残るんですね。……でもそれって、外周の形がどんなにグニャグニャでも成り立つんですか?
🟡 リナ: その通り。各微小サイクルで \(\delta Q_H/T_H + \delta Q_C/T_C = 0\) が成り立つから、全体を足し上げると外周に沿った積分だけが残って:
が成り立つ。これは「\(\delta Q_{\text{rev}}/T\) が完全微分のように振る舞う」ことを意味するの。
🔵 カイ: 完全微分って何ですか?
🟡 リナ: 山登りで例えるわね。標高差は出発点と到着点だけで決まる——どの登山道を選んでも \(h(B) - h(A)\) は同じ値になるでしょう? こういう「経路に依存しない微小変化」を完全微分と呼ぶの。数学的に言えば、「ある関数 \(S\) が存在して、その微小変化 \(dS\) として書ける量」のこと。一方、歩いた距離は経路によって変わる——こちらは不完全微分。\(\delta Q\) がまさにこれで、同じ始点と終点でも途中の過程によって値が変わるの。
⚪ メイ: つまり \(\delta Q\) 自体は「歩いた距離」で経路に依存するけど、\(T\) で割ると「標高差」のように経路に依存しなくなるのね。
🟡 リナ: いい例えね。\(\oint \delta Q_{\text{rev}}/T = 0\) が成り立つということは、\(\delta Q_{\text{rev}}/T\) を \(A\) から \(B\) まで足し上げた結果が経路に依存しないことを意味するの。なぜかを直感的に説明するわね。\(A\) から \(B\) への行き方が経路 I と経路 II の 2 通りあったとして、もし結果が違ったら——経路 I で行って経路 II の逆で戻ると一周の値がゼロにならない。でも \(\oint = 0\) と矛盾する。だからどの経路でも同じ値になるはず。次の小節で厳密に証明するわ。つまり、ある関数 \(S\) が存在して \(dS = \delta Q_{\text{rev}}/T\) と書ける——\(\delta Q_{\text{rev}}/T\) は完全微分なの。
思い出して——\(\delta Q\) 自体は経路に依存する不完全微分だった。でも \(T\) で割ると、経路によらない量に変わる。こういう「割ることで不完全微分を完全微分に変換してくれる量」を数学では積分因子と呼ぶの。ここでは \(1/T\) が積分因子の役割を果たしているわ。
✅ 理解度チェック: \(\delta Q\) は経路に依存する量だが、\(\delta Q_{\text{rev}}/T\) はどのような性質を持つでしょうか?
答え
可逆過程に沿った \(\delta Q_{\text{rev}}/T\) の積分は経路に依存しない(完全微分のように振る舞う)。これにより状態量であるエントロピー \(S\) を定義できる。
状態量としてのエントロピー¶
🟡 リナ: 可逆サイクルで一周積分がゼロということは、可逆過程に沿った \(\int \delta Q_{\text{rev}}/T\) の値が経路に依存しないことを意味する。
証明しましょう。2 つの異なる可逆経路 I と II で状態 \(A\) から状態 \(B\) へ行くとする。「経路 I で \(A \to B\)、経路 II の逆で \(B \to A\)」というサイクルを考えると:
したがって:
🔵 カイ: おお、どの経路で行っても同じ値になる! だから「状態量」と言えるんですね。
🟡 リナ: 経路によらないから、これは状態量を定義する。基準状態 \(O\) を固定して:
微小変化の形で書けば:
これがエントロピー \(S\) の熱力学的定義。\(\delta Q\) は経路に依存する不完全微分だったのに、\(T\) で割ることで経路に依存しない状態量に変身した。\(T\) は積分因子の役割を果たしているの。
熱力学の基本関係式¶
🟡 リナ: \(\delta Q_{\text{rev}} = T\,dS\) を第 1 法則に代入すると:
これが熱力学の基本関係式。\(dU\)、\(dS\)、\(dV\) は全て状態量の全微分で、不完全微分は含まれていない。
🔵 カイ: きれいですね。でも、この式は可逆過程でしか使えないんですか?
🟡 リナ: いい質問。\(U\)、\(S\)、\(V\) は全て状態量だから、この式は任意の平衡状態間の関係として成り立つ。過程が可逆か不可逆かに関係なく、2 つの平衡状態の間の差を表す式として普遍的に使えるの。
Carnot サイクルでの確認¶
🟡 リナ: Carnot サイクルでエントロピーの変化を確認しましょう。
- 等温膨張(\(T_H\)):\(\Delta S_1 = Q_H / T_H\)
- 断熱膨張:\(\Delta S_2 = 0\)(\(\delta Q = 0\) だから)
- 等温圧縮(\(T_C\)):\(\Delta S_3 = -Q_C / T_C\)
- 断熱圧縮:\(\Delta S_4 = 0\)
一周の変化:
⚪ メイ: 確かにゼロ。エントロピーは状態量だから一周すれば元に戻る。
🔵 カイ: でも、不可逆な過程だったらどうなるんですか? 一周してもゼロにならない?
🟡 リナ: いい質問。不可逆過程では \(\delta Q/T < dS\)、つまり:
等号は可逆過程のとき。断熱系(\(\delta Q = 0\))なら:
これがエントロピー増大の法則——第 2 法則のエントロピーによる表現よ。
✅ 理解度チェック: エントロピーの熱力学的定義を微小変化の形で書いてください。
答え
\(dS = \delta Q_{\text{rev}} / T\)(可逆過程で系に加えられた微小熱量を温度で割ったもの)。
3.5 Boltzmann のエントロピー — ミクロからマクロへ¶
🟡 リナ: ここまでの熱力学は、ミクロな世界の情報を一切使っていない。「熱」「温度」「エントロピー」を全てマクロな量として定義してきた。Boltzmann (ボルツマン) のアイデアは革命的だった——マクロな熱現象は、ミクロな粒子の統計的な振る舞いから説明できる。
ミクロ状態とマクロ状態¶
🔵 カイ: 統計的?
🟡 リナ: 例えば、部屋の中の空気を考えて。空気分子は約 \(10^{23}\) 個ある。一つ一つの分子の運動を追跡するのは不可能。でも、「全体としてどう振る舞うか」は確率と統計で予測できる。
ここで 2 つの概念を区別するわ。
- マクロ状態: 温度、圧力、体積など、マクロに測定可能な量で指定される状態
- ミクロ状態: 全ての粒子の位置と速度(量子力学なら量子状態)を完全に指定した状態
同じマクロ状態に対応するミクロ状態は、一般に膨大な数ある。
🔵 カイ: つまり、1 個の分子がどう動くかは分からなくても、\(10^{23}\) 個まとめて見れば予測できるってことですか?
🟡 リナ: その通り。サイコロ 1 回の結果は予測できないけど、100 万回振れば各目がほぼ均等に出ることは予測できるでしょう? 同じ原理よ。
Boltzmann のエントロピー¶
🟡 リナ: Boltzmann はエントロピーを次のように定義した:
ここで: - \(\Omega\) = あるマクロ状態に対応するミクロ状態の数 - \(k_B \approx 1.381 \times 10^{-23}\;\text{J/K}\) = Boltzmann 定数。単位が J/K(エネルギー÷温度)であることからわかるように、ミクロな世界(エネルギー)とマクロな世界(温度)を繋ぐ換算係数の役割を果たす - \(\ln\) = 自然対数
🔵 カイ: なぜ対数を取るんですか?
🟡 リナ: 加法性のため。2 つの独立な系を並べたとき、合成系のミクロ状態数は積になる:
なぜ積かというと、系 1 のどのミクロ状態に対しても系 2 の全てのミクロ状態が可能だから。サイコロ 2 つの出目の組合せが \(6 \times 6 = 36\) 通りになるのと同じ。
対数を取れば積が和に変わる:
エントロピーは足し算できる量になる。エネルギーや体積と同じ性質ね。
⚪ メイ: 対数を取るのは「掛け算を足し算に変える」ためなのね。数学的なトリックじゃなくて物理的な要請から来ているのか。
コインの例で理解する¶
🟡 リナ: 具体例で考えましょう。コインを 4 枚投げる。
表 3.2: コイン4枚のマクロ状態とエントロピー
| マクロ状態(表の枚数) | ミクロ状態の数 \(\Omega\) | エントロピー \(S/k_B = \ln\Omega\) |
|---|---|---|
| 0 枚(全部裏) | 1 | 0 |
| 1 枚 | 4 | 1.39 |
| 2 枚 | 6 | 1.79 |
| 3 枚 | 4 | 1.39 |
| 4 枚(全部表) | 1 | 0 |
🔵 カイ: 全部表とか全部裏は 1 通りしかないのに、半々だと 6 通りもあるんですね。
⚪ メイ: つまり「表 2 枚・裏 2 枚」が最も \(\Omega\) が大きくて、エントロピーが最大になるのね。
🟡 リナ: そう。グラフにすると一目瞭然よ(図 3.4「コインのミクロ状態数とエントロピー」)。真ん中がピークになっているでしょう?
図 3.4: コインのミクロ状態数とエントロピー。コイン 4 枚の場合のマクロ状態(表の枚数)ごとのミクロ状態数 \(\Omega\) とエントロピー \(S/k_B = \ln\Omega\)。最も均等な分配(表 2 枚)でエントロピーが最大になる。
✅ 理解度チェック: コイン 4 枚の例で、エントロピーが最大になるマクロ状態はどれでしょうか? その理由は?
答え
表 2 枚・裏 2 枚のマクロ状態。対応するミクロ状態の数 \(\Omega = 6\) が最大だから(\(S = k_B \ln \Omega\) が最大)。
🟡 リナ: そう。そして \(10^{23}\) 個の粒子になると、「均等に分布している状態」の \(\Omega\) は「一箇所に偏っている状態」の \(\Omega\) に比べて圧倒的に大きい。だから、系は自然に \(\Omega\) が大きい状態——つまりエントロピーが高い状態——に向かう。図 3.5「自由膨張とエントロピー増大」 で具体的なイメージを見て。箱の仕切りを外すと気体は自発的に全体に広がるわね。「各粒子が左右どちらの半分にいるか」だけに注目すると、初期状態(全粒子が左半分)の配置は 1 通りだけど、平衡状態では各粒子が左右どちらにもいられるから配置数は \(2^N\) 通りに爆発的に増える。
図 3.5: 自由膨張とエントロピー増大。仕切りを外すと気体は自発的に全体に広がる。「各粒子が左右どちらの半分にいるか」だけに注目すると、初期状態(全粒子が左半分)の配置は 1 通りだが、平衡状態(均一分布)では各粒子が左右どちらにもいられるため配置数は \(2^N\) 通りに増大する。
なぜ第 2 法則が成り立つか¶
🔵 カイ: それが「熱は高温から低温に流れる」の理由?
🟡 リナ: そう。熱が高温から低温に流れるのは、そうなった方が全体のミクロ状態の数 \(\Omega\) が増えるから。禁止されているのではなく、逆方向に起きる確率が天文学的に小さいだけ。
🔵 カイ: つまり「禁止」じゃなくて「ほぼ確実にそうなる」ってことか……。でも「ほぼ確実」って、具体的にどのくらいの確率なんですか? 宝くじに当たるくらい? それとも桁違いに小さい?
🟡 リナ: 具体的に見積もってみましょう。\(N\) 個の気体分子が箱の中にいるとして、全ての分子が箱の左半分に集まる確率は:
\(N = 10^{23}\) なら:
🔵 カイ: \(10^{-3 \times 10^{22}}\)……宝くじどころの話じゃない。ゼロが \(10^{22}\) 個並ぶ確率って、もはや「起こらない」と同じですね。
🟡 リナ: 宇宙の年齢(約 \(10^{10}\) 年 \(\approx 10^{17}\) 秒)の間、毎秒 1 回試行しても、一度も起こらない。
🔭 科学哲学メモ: 第 2 法則は「絶対に破れない法則」ではなく、「破れる確率が天文学的に小さい統計的法則」である。これは Newton の運動方程式のような決定論的法則とは本質的に異なる性格を持つ。自分で判断してほしい——「確率が \(10^{-10^{22}}\) の事象は起こらない」と言い切ってよいか?
📝 練習問題:
- コインのミクロ状態数とエントロピー → 問題 M-1. コインのエントロピー
✅ 理解度チェック: Boltzmann のエントロピーの定義式を書き、\(\Omega\) が何を表すか答えてください。
答え
\(S = k_B \ln \Omega\)。\(\Omega\) はミクロ状態の数(マクロな状態が同じになるようなミクロな配置の総数)。
✅ 理解度チェック: 熱が高温から低温に流れる理由を、ミクロ状態の数を用いて説明してみましょう。
答え
熱が高温から低温に流れた方が全体のミクロ状態の数 \(\Omega\) が増えるため。逆方向に起きる確率が天文学的に小さいだけで、禁止されているわけではない。
3.6 理想気体のミクロ状態数 — 具体的な計算¶
🟡 リナ: \(S = k_B \ln \Omega\) が抽象的に見えるかもしれないから、理想気体で具体的に \(\Omega\) を計算してみましょう。
問題設定¶
🟡 リナ: \(N\) 個の同種粒子(質量 \(m\))が体積 \(V\) の箱に閉じ込められていて、全エネルギーが \(E\) の場合を考える。
古典力学では、各粒子の状態は位置 \((x, y, z)\) と運動量 \((p_x, p_y, p_z)\) の 6 つの数で完全に指定される。運動量は「質量×速度」で、各方向の成分 \(p_x = mv_x\), \(p_y = mv_y\), \(p_z = mv_z\) を持つベクトル量よ。
🔵 カイ: 位置だけじゃダメなんですか?
🟡 リナ: ダメなの。同じ位置にいても、速く動いている粒子と止まっている粒子では状態が違うでしょう? だから「どこにいるか」と「どう動いているか」の両方が必要。1 粒子なら \((x, y, z, p_x, p_y, p_z)\) の 6 個——これを 6 本の座標軸とする 6 次元空間を考えると、その中の 1 つの点が「粒子がどこにいて、どう動いているか」を完全に表す。粒子が 2 個なら \(6 \times 2 = 12\) 個の数で系全体の状態が決まる。\(N\) 個の粒子なら、全体で \(6N\) 個の変数が必要ね。この \(6N\) 個の変数を座標軸とする \(6N\) 次元の空間を位相空間(phase space)と呼ぶの。「位相」は波の位相とは無関係よ——名前の由来には諸説あるけど、ここでは「系の状態を完全に表す空間」という意味だと思っておいて。位相空間の 1 つの点が、系全体の 1 つのミクロ状態に対応するわ。
⚪ メイ: つまり、\(10^{23}\) 個の粒子なら \(6 \times 10^{23}\) 次元の空間の中の 1 点で、系全体が完全に記述されるのね。途方もない次元だけど概念はシンプル。
エネルギーの制約¶
🟡 リナ: 理想気体では粒子間の相互作用がないから、全エネルギーは運動エネルギーの和:
この式の構造を見やすくするために、\(3N\) 個の運動量成分を \(\xi_1, \xi_2, \ldots, \xi_{3N}\)(\(\xi_1 = p_{x,1}\), \(\xi_2 = p_{y,1}\), \(\xi_3 = p_{z,1}\), \(\xi_4 = p_{x,2}\), …)と一列に並べ直すわ。すると全エネルギーの条件は:
🔵 カイ: \(\xi_1^2 + \xi_2^2 + \cdots = 2mE\) って、何かの球みたいな形ですね? 3 次元なら \(x^2 + y^2 + z^2 = R^2\) が球面だから……。
🟡 リナ: いい直感ね。まさにその通りで、これは \(3N\) 次元版の球面の方程式——半径 \(\sqrt{2mE}\) の球面なの。つまり、エネルギーが決まれば運動量の「配り方」が球面上に制限されて、エネルギーが大きいほど球面が大きくなるから状態数が増える。エネルギーの制約が幾何学的な制約——球面の大きさ——に翻訳されるわけ。
⚪ メイ: つまり、エネルギーの制約が「球面の半径」という幾何学的な条件に置き換わるのね。
🟡 リナ: そう。図 3.6「運動量空間の球面とミクロ状態数」 にイメージを描いたわ。
図 3.6: 運動量空間の球面とミクロ状態数。理想気体のミクロ状態は運動量空間で半径 \(\sqrt{2mE}\) の球面(の内部)に対応する。エネルギーが増えると球が大きくなり、ミクロ状態数が増える。
エネルギーがちょうど \(E\) の状態は球面上にあるけど、ミクロ状態の数を数えるときは球の内部(エネルギーが \(E\) 以下の全状態)の体積を使うの。実は高次元では球の体積のほとんどが表面付近に集中するの。直感的に言うと、\(d\) 次元球の体積は \(R^d\) に比例するから、半径 \(R\) の球の体積に対して半径 \(0.99R\) の球の体積の割合は \((0.99)^d\)——\(d = 3\) なら約 97%だけど、\(d = 10^{23}\) なら \((0.99)^{10^{23}} \approx 0\)、つまり体積のほぼ全てが外側 1%の殻に集中する。だから球の内部の体積で数えても球殻の表面積で数えても結果は同じになる。ここでは球の内部の体積に比例するとして計算を進めるわ。
🔵 カイ: 高次元だと体積が表面に集中する……直感とは違いますね。でも \((0.99)^{10^{23}} \approx 0\) って言われると納得です。
\(3N\) 次元球の体積¶
🔵 カイ: \(3N\) 次元の球の体積ってどうやって求めるんですか? 3 次元の球なら \(\frac{4}{3}\pi R^3\) ですけど……。
🟡 リナ: いい質問。ミクロ状態の数を求めるには、運動量空間で半径 \(\sqrt{2mE}\) の球の体積を計算する必要がある。3 次元なら球の体積は \(\frac{4}{3}\pi R^3\) よね——\(R^3\) に比例する。同じ発想で \(d\) 次元に拡張すると、多重積分を使って計算でき(導出は省略するけど):
という公式が得られるの。ここで \(\Gamma\) はガンマ関数——階乗 \(n!\) を整数以外にも拡張した関数で、整数 \(n\) に対しては \(\Gamma(n+1) = n!\) を満たすもの。なぜ \(\pi\) が出てくるかというと、2 次元の円の面積 \(\pi R^2\) を思い出して——\(\pi\) は「角度方向をぐるっと一周する」ことから来ているわね。高次元でも同じで、球の体積を多重積分で計算するとき、角度方向を全方向にわたって積分する部分から \(\pi\) のべきが現れ、動径方向(中心からの距離)の積分から階乗(の一般化であるガンマ関数)が自然に現れるの。
⚪ メイ: \(\pi\) は「ぐるっと一周」の幾何から、\(\Gamma\) は「半径方向に積み上げる」計算から——それぞれ別の由来で出てくるのね。
🟡 リナ: その通り。公式を丸暗記する必要はないわ——大事なのは \(R^d\) に比例するという部分。つまり次元 \(d\) が大きいほど、半径 \(R\) が少し増えただけで体積が爆発的に増えるということ。低次元で確認してみましょう。\(d = 1\) なら \(V_1 = \frac{\pi^{1/2}}{\Gamma(3/2)} R = \frac{\sqrt{\pi}}{\frac{1}{2}\sqrt{\pi}} R = 2R\)——これは長さ \(2R\) の線分で正しい。\(d = 2\) なら \(V_2 = \frac{\pi}{\Gamma(2)} R^2 = \pi R^2\)——円の面積。\(d = 3\) なら \(V_3 = \frac{\pi^{3/2}}{\Gamma(5/2)} R^3 = \frac{\pi^{3/2}}{\frac{3}{4}\sqrt{\pi}} R^3 = \frac{4}{3}\pi R^3\)——球の体積。ちゃんと合うわね。
🔵 カイ: \(\Gamma(d/2 + 1)\) って、\(d\) が奇数だったら \(\Gamma(3/2)\) とか \(\Gamma(5/2)\) みたいに整数じゃない値が入りますよね? それって計算できるんですか?
🟡 リナ: いい指摘。ガンマ関数は階乗 \(n!\) を整数以外にも拡張した関数で、漸化式 \(\Gamma(n+1) = n \cdot \Gamma(n)\) を満たすの。\(\Gamma(1/2) = \sqrt{\pi}\) という値が知られていて、そこから \(\Gamma(3/2) = \frac{1}{2}\sqrt{\pi}\)、\(\Gamma(5/2) = \frac{3}{4}\sqrt{\pi}\) と順に計算できる。でも今大事なのはこの係数の詳細じゃなくて、体積が \(R^d\) に比例するという部分——つまり次元 \(d\) が大きいほど、半径が少し増えただけで体積が爆発的に増えるということ。ガンマ関数の詳しい性質は必要になったときに改めて学ぶわ。\(R = \sqrt{2mE}\) を代入すると:
ミクロ状態数の組み立て¶
🟡 リナ: 次に位置の自由度を考えましょう。1 個の粒子は体積 \(V\) のどこにでもいられるから、位置の「選択肢」は \(V\) に比例する。2 個なら、粒子 1 が \(V\) のどこか、粒子 2 も独立に \(V\) のどこか——だから \(V \times V = V^2\)。\(N\) 個なら \(V^N\)。つまり位置の自由度は \(V^N\) の因子を与えるの。
🔵 カイ: 運動量空間の球の体積と \(V^N\) を掛ければミクロ状態数になるんですか?
🟡 リナ: もう少し補正が必要なの。まず、位相空間の「体積」は連続的だから、そのままでは状態の「数」にならない——体積は実数値で、「何通り」とは数えられないわ。
🔵 カイ: 確かに。「面積が 10 平方メートル」と言われても、「何マス分」とは言えないですよね。マス目のサイズを決めないと。
🟡 リナ: その通り。状態を数えるには、位相空間を何らかの最小単位で区切って「マス目の数」にする必要がある。ちょうど方眼紙のマス目のサイズを決めるようなものね。では、そのマス目のサイズは何で決まるか? ここで量子力学(次の章以降で詳しく学ぶ)が登場するの。量子力学には「位置と運動量を同時に完全に正確に決めることはできない」という原理があって——たとえば粒子の位置を精密に測ろうとすると、運動量の情報がぼやけてしまう。この原理の帰結として、位相空間にはある最小の面積より細かい区別ができない「最小マス目」があるの。その最小面積を決める定数が \(h\)(\(h \approx 6.626 \times 10^{-34}\;\text{J·s}\)、Planck 定数と呼ばれる)——自然界が「これ以上細かくは区別しない」と決めている基本定数よ。なぜこの値なのか、なぜこんな原理があるのかは次の章以降で詳しく学ぶわ。
🔵 カイ: つまり、\(h\) の具体的な値は実験で決まっていて、今は「自然界が決めたマス目のサイズ」として受け入れればいいんですね?
🟡 リナ: その通り。1 方向あたり \(\Delta x \cdot \Delta p_x\) の最小値が \(h\) 程度で決まる。1 粒子は \(x, y, z\) の 3 方向それぞれに位置と運動量を持つから、最小区画は \(h \times h \times h = h^3\)(各方向が独立なので積になる)。\(N\) 粒子なら \(h^{3N}\)。だから位相空間の体積を \(h^{3N}\) で割って「状態の数」に変換するの。
さらに、同じ種類の粒子は入れ替えても区別できない——これは量子力学から来る要請で、次の章以降で詳しく学ぶわ。直感的に言うと、ヘリウム原子 A とヘリウム原子 B に「名札」を貼ることは実験的にできない——どちらも全く同じ性質を持つから、交換しても物理的に区別がつかないの。今は「同じ種類の粒子にラベルを貼って区別することはできない」と受け入れて。粒子 1 と粒子 2 を交換しても同じミクロ状態。\(N\) 個の粒子の並べ替えは \(N!\) 通りあるから、重複を除くために \(N!\) で割る。全部合わせると:
(\(N\) が大きいとき、\(\Gamma(3N/2+1)\) は次の小節で紹介する Stirling の近似で扱えるので、整数かどうかを気にする必要はないわ。)
🔵 カイ: 複雑ですね……。
🟡 リナ: 大事なのは \(E\) と \(V\) への依存性。\(N\) を固定して \(E\) と \(V\) だけの関数として見ると(\(N\) に依存する因子は全て定数扱い):
⚪ メイ: 体積を広げるか、エネルギーを増やすか——どちらでもミクロ状態の数が増えるのね。シンプルな結果で安心するわ。
エントロピーの計算¶
🟡 リナ: 対数を取ると:
このままだと \(\ln(N!)\) や \(\ln\Gamma(3N/2+1)\) が扱いにくいわね。ここで Stirling の近似を使うわ。\(n\) が大きいとき、\(n!\) の対数は:
と近似できる。\(\ln n! = \ln 1 + \ln 2 + \cdots + \ln n\) は \(\ln x\) の積分 \(\int_1^n \ln x\,dx = n\ln n - n + 1 \approx n\ln n - n\) で近似できるからよ。
🔵 カイ: 階乗の対数が \(n \ln n - n\) で近似できるのは便利ですね。\(10^{23}!\) なんて直接計算できないですもんね。
🟡 リナ: これを適用すると、\(\ln(N!) \approx N\ln N - N\)、\(\ln\Gamma(3N/2+1) \approx \frac{3N}{2}\ln\frac{3N}{2} - \frac{3N}{2}\)。整理すると \(N\ln V - N\ln N = N\ln(V/N)\) や \(\frac{3N}{2}\ln(2mE) - \frac{3N}{2}\ln\frac{3N}{2} = \frac{3N}{2}\ln\frac{4mE}{3N}\) のように変形できる。さらに \(\ln\frac{4mE}{3N} = \ln\frac{E}{N} + \ln\frac{4m}{3}\) と分解すると、\(\ln\frac{4m}{3}\) の部分は \(E\) にも \(V\) にも依存しない定数だから、\(m\), \(h\) を含む他の定数項と一緒にまとめられる(対数の引数が無次元になるのはこれらの定数と組み合わせたときであることに注意)。最終的に Sackur-Tetrode の公式と呼ばれる結果が得られる:
ここで (定数) は \(m\), \(h\), \(k_B\) などから構成される項を全て含んでいるの。「\(\ln(E/N)\) の中身はエネルギーの次元を持つのに、対数に入れていいの?」と思うかもしれないけど、実際には定数の中に \(\frac{3}{2}\ln(4\pi m / (3h^2))\) のような項が含まれていて、合わせると \(\ln\) の引数は無次元になる。ただし \(E\) や \(V\) で偏微分する際にはこれらの定数項は消えるので、以下の計算には影響しない。だから今は「定数部分は気にしなくてよい」と思って大丈夫。
⚪ メイ: \(E\) が大きいほど、\(V\) が大きいほど、エントロピーが大きい。直感と合いますね。
🟡 リナ: そう。エネルギーが増えれば運動量空間の球が大きくなり、体積が増えれば位置空間が広がる。どちらもミクロ状態の数を増やすの。
✅ 理解度チェック: 理想気体のミクロ状態数 \(\Omega\) はエネルギー \(E\) にどのように依存するでしょうか?
答え
\(\Omega \propto E^{3N/2}\)。エネルギーが増えると運動量空間の球の体積が増えるため。
3.7 温度の統計力学的な意味¶
🟡 リナ: Boltzmann のエントロピーを使うと、温度にも統計力学的な意味が与えられる。
2 つの系の熱平衡条件からの導出¶
🟡 リナ: 2 つの系(系 1 と系 2)がエネルギーを交換できる状況を考えましょう。全エネルギーは保存される:
合成系のミクロ状態数は:
ここで一つ重要な仮定を置くわ。等重率の原理——孤立系では、エネルギーなどの条件を満たす全てのミクロ状態が等しい確率で実現される、という仮定よ。
🔵 カイ: なんで等確率って言えるんですか?
🟡 リナ: 深い問いね。厳密な証明はないの。でも「特定のミクロ状態を優遇する理由がない」という対称性の議論と、この仮定から導かれる結果が実験と一致するという事実が、この原理を支えている。統計力学の出発点となる仮定だと思って。
🔵 カイ: それがどう \(\Omega\) の最大化に繋がるんですか?
🟡 リナ: 各ミクロ状態が等確率なら、あるマクロ状態が実現する確率は、そのマクロ状態に対応するミクロ状態の数 \(\Omega\) に比例する。だから \(\Omega\) が最大のマクロ状態が最も高い確率で実現される——つまり系は \(\Omega_{\text{total}}\) が最大になるエネルギー分配に向かうの。
図 3.7: 熱平衡とエントロピー最大化。2 つの系がエネルギーを交換するとき、全エントロピー \(S_{\text{total}} = S_1(E_1) + S_2(E_{\text{total}} - E_1)\) が最大になる点で平衡に達する。この点で \(\partial S_1/\partial E_1 = \partial S_2/\partial E_2\)、すなわち \(T_1 = T_2\)。
🟡 リナ: \(\Omega_{\text{total}}\) を最大化する条件を求めればいいの。図 3.7「熱平衡とエントロピー最大化」 にイメージを描いたわ。対数を取った方が計算しやすいから、エントロピーで書くと:
これを \(E_1\) で微分してゼロとおく:
\(E_2 = E_{\text{total}} - E_1\) だから \(\partial E_2 / \partial E_1 = -1\)。したがって:
🔵 カイ: 平衡では、両方の系で \(\partial S / \partial E\) が等しくなる!
🟡 リナ: そう。そして私たちは「熱平衡では温度が等しくなる」ことを知っている。だから、\(\partial S / \partial E\) は温度に関係する量のはず。実際:
これが温度の統計力学的定義。
✅ 理解度チェック: 2つの系がエネルギーを交換して熱平衡に達する条件を、エントロピーを用いて表せ。
答え
\(\frac{\partial S_1}{\partial E_1} = \frac{\partial S_2}{\partial E_2}\)(両系でエントロピーのエネルギー微分が等しくなる)。これは温度が等しくなることに対応する。
理想気体での確認¶
🟡 リナ: 先ほど求めた理想気体のエントロピーは \(S = k_B N \left[\frac{3}{2}\ln \frac{E}{N} + \ln\frac{V}{N} + (\text{定数})\right]\) だったわね。\(E\) で偏微分するとき、\(\ln(E/N) = \ln E - \ln N\) だから \(\ln N\) の部分は定数として消える。微分してみましょう:
整理すると:
🔵 カイ: \(E = \frac{3}{2}Nk_BT\) って、さっき断熱過程の計算で使った式と同じですね! あのとき天下りで使ったものが、ここで導出されたってことですか?
🟡 リナ: その通り。さっきは「知られている事実」として使ったけど、今まさに統計力学から導出できたわ。\(3\) は各粒子が持つ自由度の数(\(x, y, z\) 方向の運動)で、各自由度あたり \(\frac{1}{2}k_B T\) のエネルギーを持つ。これをエネルギー等分配則と呼ぶの。
🔵 カイ: じゃあ、もし分子が回転もできたら自由度が増えて、同じ温度でもエネルギーが大きくなるってことですか?
🟡 リナ: その通り。2 原子分子なら回転の自由度が 2 つ加わって、\(E = \frac{5}{2}Nk_BT\) になる。等分配則の威力ね。さらに言えば、統計力学の定義 \(1/T = \partial S/\partial E\) から出発するだけで、こうしたおなじみの公式が全て導出できるの。同様に、エントロピーの \(V\) 微分から圧力も出る:
整理すると:
🔵 カイ: 理想気体の状態方程式だ! でも待ってください。\(\partial S/\partial V\) を計算するとき、\(E\) を固定してますよね。温度じゃなくてエネルギーを固定するのはなぜですか?
🟡 リナ: いい質問。ここでは孤立系(外部とエネルギーのやり取りがない系)を考えているから、\(E\) が一定なの。温度は \(E\) から導出される量だから、出発点では \(E\) を固定するのが自然なのよ。
⚪ メイ: つまり、\(S = k_B \ln \Omega\) からエントロピーを求めて、\(E\) で偏微分すれば温度、\(V\) で偏微分すれば圧力が出てくるのね。
🟡 リナ: その通り。ミクロ状態を数えてエントロピーを求め、偏微分するだけで、温度も圧力も状態方程式も全て導出できる。マクロな量が全てミクロ状態の数え上げから出てくる——これが統計力学の威力よ。
✅ 理解度チェック: 理想気体のエントロピーから \(E = \frac{3}{2}Nk_BT\) を導出する際、どのような操作を行うでしょうか?
答え
エントロピー \(S = k_B N [\frac{3}{2}\ln E + \cdots]\) をエネルギー \(E\) で偏微分して \(1/T = \partial S/\partial E\) を計算し、\(E\) について解く。
温度の直感的な意味¶
🟡 リナ: 温度の定義 \(1/T = \partial S / \partial E\) を言葉にすると:
温度が低い系ほど、エネルギーを少し加えたときにエントロピー(状態数の対数)が大きく増える。
高温の物体はすでに「乱雑」だから、エネルギーを追加しても新しい状態がそれほど開かない。低温の物体はまだ「整然」としているから、少しのエネルギーで多くの新しい状態が開く。
🔵 カイ: 温度が低いほど「伸びしろ」が大きいってことですね。エネルギーを渡す効果が大きい。
🟡 リナ: そうね。図 3.8「温度の統計力学的意味」 のグラフで視覚的に確認しよう。\(S\)-\(E\) グラフでは、同じ系でもエネルギーが小さい領域(低温)では傾き \(\partial S/\partial E\) が急で、エネルギーが大きい領域(高温)では傾きが緩やかになるの。
図 3.8: 温度の統計力学的意味。\(S\)-\(E\) グラフの傾き \(\partial S/\partial E = 1/T\) が温度の逆数を与える。傾きが急な領域(青)は低温で、同じ \(\Delta E\) に対してエントロピーが大きく増える。傾きが緩やかな領域(赤)は高温で、エネルギーを加えてもエントロピーがあまり増えない。
だから、エネルギーが高温の系から低温の系に流れると、高温側でのエントロピー減少よりも低温側でのエントロピー増加の方が大きく、全体としてエントロピーが増える。
🔵 カイ: だから熱は高温から低温に流れるんですね! でも、もし系が小さくて粒子が 10 個くらいしかなかったら、逆方向に流れることも観測できるんですか?
🟡 リナ: 鋭いわね。実際、ナノスケールの系では「揺らぎ」として第 2 法則からのずれが観測されている。粒子数が少ないと \(\Omega\) の比が圧倒的ではなくなるから。でもマクロな系(\(N \sim 10^{23}\))では事実上起こらない。
🔵 カイ: へえ、じゃあ第 2 法則って、粒子数が多いからこそ成り立つ「近似的な法則」なんですね。
🟡 リナ: その通り。Newton の法則は「個々の粒子に対して決定論的に成り立つ」けど、第 2 法則は「\(N\) が大きいから事実上成り立つ統計的な法則」——法則の性格が根本的に違うの。物理法則には少なくとも 2 種類の「成り立ち方」がある——粒子 1 個でも厳密に成り立つ決定論的法則と、\(N\) が大きいときにだけ事実上成り立つ統計的法則。第 2 法則は後者だから、ナノスケールでは破れが見えるけど、マクロでは絶対的に見えるの。
⚪ メイ: つまり、同じ「物理法則」でも成り立つ仕組みが全然違うのね。決定論的な法則と統計的な法則——区別して考える必要があるということか。
📝 練習問題:
- 熱平衡条件から温度の等しさを導出する → 問題 M-2. 温度の統計力学的定義
✅ 理解度チェック: 温度の統計力学的な定義式を書いてください。
答え
\(\frac{1}{T} = \frac{\partial S}{\partial E}\bigg|_{V,N}\)(エントロピーをエネルギーで偏微分したものが温度の逆数)。
✅ 理解度チェック: 低温の物体にエネルギーを加えると、ミクロ状態の数はどう変化するでしょうか?
答え
大きく増える(\(\partial S / \partial E\) が大きいので \(T\) が小さい)。
3.8 自由エネルギー \(F = U - TS\)¶
🟡 リナ: 最後にもう一つ重要な概念を導入しておくわ。Helmholtz の自由エネルギー。
動機:温度一定の系¶
🟡 リナ: 実験室では、系を大きな熱浴に接触させて温度を一定に保つことが多い。このとき、系のエントロピーは一定ではない(熱浴とエネルギーをやり取りするから)。代わりに温度 \(T\) が一定。
こういう状況で「系はどの状態に落ち着くか?」を判定する量が自由エネルギー:
\(F\) の最小化原理¶
🟡 リナ: 系と熱浴の全系を考えると、全系のエントロピーは増大する(第 2 法則)。ここで系の体積は一定(仕事をしない)と仮定するわ——実験室で容器に入れた系を熱浴に浸す典型的な状況ね。この条件のもとで、全系のエントロピー変化から系だけの量で書かれた条件を導出しましょう。
🟡 リナ: 熱浴のエントロピー変化を求めましょう。ポイントは 2 つ。
- 熱浴は非常に大きいので、系からエネルギーを受け取っても温度がほとんど変わらない。巨大なプールに一滴の熱湯を垂らすイメージね——プールの温度はほとんど変わらないでしょう? つまり熱浴は常にほぼ平衡状態にあり、系とのエネルギーのやり取りは熱浴にとって「ほんの微小な変化」。だから熱浴の側では \(dS = \delta Q_{\text{rev}}/T\) がそのまま使える——熱浴が受け取った熱 \(\Delta Q_{\text{bath}}\) に対して \(\Delta S_{\text{bath}} = \Delta Q_{\text{bath}}/T\) と書けるの。(厳密に言えば「準静的で散逸がない過程は可逆」という条件が必要だけど、巨大な熱浴にとっては微小な変化だから、この条件は自動的に満たされると考えてよいわ。)
- 体積一定だから系は仕事をしない(\(\Delta W = p\Delta V = 0\))。第 1 法則 \(\Delta U = \Delta Q - \Delta W\) で \(\Delta W = 0\) とすると、系のエネルギー変化は \(\Delta U_{\text{sys}} = \Delta Q_{\text{sys}}\)(系が受け取った熱がそのまま内部エネルギーの変化になる)。エネルギー保存から、系が受け取った熱は熱浴が失った熱に等しいので \(\Delta Q_{\text{bath}} = -\Delta Q_{\text{sys}} = -\Delta U_{\text{sys}}\)
🔵 カイ: なるほど、系がエネルギーをもらった分だけ熱浴が失う——エネルギー保存ですね。
🟡 リナ: まとめると:
\(\Delta S_{\text{total}} \geq 0\) を要求すると:
両辺に \(-T\) をかけて(\(T > 0\) だから不等号が反転):
ここで温度 \(T\) は一定(熱浴が温度を固定している)だから、\(\Delta U - T\Delta S = \Delta(U - TS)\) と書ける。つまり:
🟡 リナ: つまり、温度一定・体積一定では、系は自由エネルギー \(F\) を最小化する方向に変化するの。
⚪ メイ: エントロピー増大の法則を系だけの量で書き直すと、\(F\) の減少になるのね。
\(F\) の物理的意味¶
🟡 リナ: \(F = U - TS\) の意味を考えてみて。
- \(U\) はエネルギーを下げたい(安定性)
- \(-TS\) はエントロピーを上げたい(乱雑さ)
\(F\) はこの 2 つの競合を一つの量にまとめたもの。低温(\(T\) が小さい)ではエネルギー最小化が支配的で、系は秩序状態を好む。高温(\(T\) が大きい)ではエントロピー最大化が支配的で、系は乱雑な状態を好む。
✅ 理解度チェック: 自由エネルギー \(F = U - TS\) において、低温と高温でそれぞれどちらの項が支配的になるでしょうか?
答え
低温では \(U\)(エネルギー最小化)が支配的で系は秩序状態を好み、高温では \(-TS\)(エントロピー最大化)が支配的で系は乱雑な状態を好む。
🔵 カイ: 水が低温で氷(秩序)になり、高温で水蒸気(乱雑)になるのは、この競合の結果なんですね。でも、ちょうど境目——氷が溶ける温度とか——はどうやって決まるんですか?
🟡 リナ: いい質問。それは \(F\) が最小になる状態が「秩序状態」から「乱雑な状態」に切り替わる温度で、これが相転移温度よ。温度を上げていくと、ある温度で \(-TS\) の項が \(U\) の項に勝ち始めて、乱雑な状態の方が \(F\) が小さくなる——その切り替わりの瞬間が相転移。図 3.9「自由エネルギーにおけるエネルギーとエントロピーの競合」 にこの競合のイメージを示したわ。
図 3.9: 自由エネルギーにおけるエネルギーとエントロピーの競合。\(F = U - TS\) において、低温ではエネルギー最小化(秩序)が支配し、高温ではエントロピー最大化(乱雑)が支配する。
🟡 リナ: 相転移の本質は \(F\) の最小化条件が温度によって変わることなの。
微分形¶
🟡 リナ: \(F\) の微小変化を計算すると:
\(dU = T\,dS - p\,dV\) を代入すると:
🔵 カイ: おお、\(T\,dS\) がきれいに消えた!
🟡 リナ: \(F\) の自然な変数は \((T, V)\)。偏微分から:
⚪ メイ: \(F\) を一つ知っていれば、\(T\) で微分してエントロピー、\(V\) で微分して圧力——全部出てくるのね。
✅ 理解度チェック: 自由エネルギー \(F = U - TS\) が最小化されるのはどのような条件のときでしょうか?
答え
温度一定・体積一定の条件で、系は自由エネルギー \(F\) を最小化する方向に変化する。
3.9 「必要性」から「宇宙の根本法則」へ¶
🟡 リナ: ここで一歩引いて、この章の物語を振り返りましょう。
🔵 カイ: 蒸気機関の効率を上げたいという実用的な問いから始まった。
⚪ メイ: そこから Carnot の効率限界が出て、第 2 法則、Boltzmann の \(S = k_B \ln \Omega\) まで繋がったわね。
🔵 カイ: でも不思議ですよね。Carnot は蒸気機関の効率を考えていただけなのに、なぜ宇宙全体に通用する法則に辿り着いたんですか?
🟡 リナ: いい問いね。Carnot が「理想的な限界」を追求したからよ。具体的な機械の細部を捨てて、「原理的に何が可能か」を問うた瞬間、問題は特定の機械を超えて普遍的になった。蒸気機関という実用的な必要性から始まった探求が、最終的に宇宙の根本的な性質——エントロピーは増大する——に到達した。プロローグで「動機は必要性と好奇心の 2 種類」と言ったけど、必要性から始まっても、好奇心と同じくらい深い場所に到達できるの。
✅ 理解度チェック: この章の物語は、どのような実用的問いから始まり、最終的にどのような宇宙の根本的性質に到達したでしょうか?
答え
蒸気機関の効率を上げたいという実用的な問いから始まり、最終的に「エントロピーは増大する」という宇宙の根本的な性質に到達した。
3.10 伏線:エントロピーとブラックホール¶
%%{init: {"theme": "default", "themeCSS": ".edgePath .path, .flowchart-link { stroke-width: 2px !important; }"}}%%
flowchart TD
A["蒸気機関の効率<br/>(実用的動機)"] --> B["Carnot サイクル<br/>η = 1 − T_C/T_H"]
B --> C["熱力学の第2法則<br/>エントロピー増大"]
C --> D["Boltzmann<br/>S = k_B ln Ω"]
D --> E["統計力学<br/>ミクロ状態の数え上げ"]
E --> F["Bekenstein-Hawking<br/>ブラックホールエントロピー<br/>S_BH = A/(4ℓ_P²) k_B"]
F --> G["Strominger-Vafa (1996)<br/>弦理論でミクロ状態を計算<br/>(第20章)"]
style A fill:#ffa,stroke:#333
style G fill:#afa,stroke:#333
図 3.10: エントロピー概念の発展系譜
🟡 リナ: 最後に一つ、伏線を張っておくわね。
🔵 カイ: 伏線?
🟡 リナ: この章で学んだエントロピーの概念は、後の章で再登場する。1970 年代に Bekenstein (ベケンシュタイン) と Hawking (ホーキング) が、ブラックホールにもエントロピーがあることを示した:
ここで \(A\) は事象の地平面の面積、\(\hbar = h/(2\pi)\) はディラック定数(Planck 定数 \(h\) を \(2\pi\) で割ったもの。量子力学では \(h\) よりこちらの方が頻繁に現れる——なぜかは次の章以降で分かるわ。今は「\(h\) の仲間の定数」と思っておけば大丈夫)、\(\ell_P = \sqrt{G\hbar/c^3}\) は Planck 長さ(後の章で再び詳しく扱う)。
🔵 カイ: ブラックホールにエントロピーがある……ということは、\(S = k_B \ln \Omega\) から考えると、ブラックホールにもミクロ状態がある?
🟡 リナ: まさにその問いが、弦理論の最大の成功に繋がるの。1996 年に Strominger と Vafa が、弦理論を使ってブラックホールのミクロ状態を数え上げ、\(S = k_B \ln \Omega\) から Bekenstein-Hawking のエントロピーを導出した(第 20 章)。
⚪ メイ: 蒸気機関から始まった概念が、ブラックホールに繋がるなんて!
🟡 リナ: 物理学の面白いところよね。一見無関係に見える分野が、深いところで繋がっている。でも注意して——Bekenstein-Hawking エントロピーは一般相対論と量子力学の組み合わせから出てくるけど、そのミクロ的な起源を弦理論で説明できたのは特定のブラックホール(BPS ブラックホール)に限られている。一般のブラックホールへの拡張は未解決問題よ。
🔭 科学哲学メモ: 弦理論によるブラックホールエントロピーの導出は美しい成果だが、これだけでは弦理論の「実験的検証」にはならない。なぜなら、同じ結果を別のモデル(ループ量子重力など)でも導出できる可能性があるから。「正しい答えを出すモデル」と「唯一正しいモデル」は異なる。反証可能性の観点から、自分で判断してほしい。
✅ 理解度チェック: ブラックホールにもエントロピーがあることを示したのは誰でしょうか?
答え
Bekenstein(ベケンシュタイン)と Hawking(ホーキング)。
次章予告¶
第 4 章 ——Newton 力学、電磁気学、熱力学という 3 つの成功したモデルが、19 世紀末に次々と破綻する。黒体輻射、光電効果、水星の近日点移動。この「危機」が、20 世紀の二大革命——相対論と量子力学——への扉を開く。特に黒体輻射の問題は、この章で学んだ熱力学と量子論の接点として重要な役割を果たす(詳しくは「量子力学」編の「量子力学」編 第 1 章を参照)。
参考文献¶
この章の内容は以下の文献を参考に構成した。
- David Tong, Lectures on Statistical Physics, Ch.1: "The Fundamentals of Statistical Mechanics" — ミクロカノニカルアンサンブル、等重率の原理
- David Tong, Lectures on Statistical Physics, Ch.2: "Classical Gases" — Boltzmann エントロピー、温度の統計力学的定義
- David Tong, Lectures on Statistical Physics, Ch.4: "Classical Thermodynamics" — 熱力学の法則、Carnot サイクル
- David Tong, Lectures on Statistical Physics, Ch.5: "Phase Transitions" — 自由エネルギー
このページについてフィードバック
分からなかった箇所、誤りの指摘、改善提案などをお寄せください。







