院試 統計力学の出題傾向と対策

統計力学は、物理系院試において量子力学と並ぶ「逃げ場のない必答科目」です。試験中の本当の負荷は計算量ではなく、熱力学量とミクロな状態和をつなぐ論理を毎問ごとに自分で組み立てる必要があることにあります。学部の熱力学で習った F・G・Ω といった熱力学関数がどの自然変数に対して何の役割を果たすか、それが分配関数の対数とどう対応するか——この接続が頭の中で自動化できていないと、設問ごとに「どのアンサンブルを取り、何を変数にして偏微分すべきか」を毎回ゼロから考え直すことになり、時間切れで失点します。逆に、この接続さえ自動化できていれば、計算量の多い問題でも答案の流れは一定です。本記事は、外部生が統計力学をどう「概念と計算の橋渡し」として整理し、対策していくかを実務的に整理した手引きです。

学部段階で固めておくべき基礎

• 熱力学の関数関係：内部エネルギー U(S, V, N)、ヘルムホルツ自由エネルギー F(T, V, N)、ギブズ自由エネルギー G(T, P, N)、グランドポテンシャル Ω(T, V, μ) の自然変数とルジャンドル変換の対応。これが脳内表として即座に出ないと、どの分配関数からどの熱力学量を取り出すかで毎回混乱する
• 組合せ論と Stirling 近似：ln N! ≈ N ln N − N の使いどころと、N が大きいときに対数を取って良い理由。ミクロカノニカルの状態数計算で必須
• 量子力学の状態数：1次元井戸・調和振動子・自由粒子のエネルギー準位と縮退度。古典極限と高温極限の取り方
• 確率論の基礎：大数の法則と中心極限定理、揺らぎが N^(−1/2) で減衰するという定性的感覚。揺らぎと応答の問題で配点に直結
• 偏微分と積分のテクニック：ガウス積分、ゾンマーフェルト展開、リーマンゼータ関数の値（ζ(2), ζ(3/2), ζ(3) 程度の数値）。電子比熱・BEC 転移温度で必須
• ラグランジュの未定乗数法：最大エントロピー原理から分配関数を導出する流れで、温度 β と化学ポテンシャル μ が乗数として現れる構造の理解

特に熱力学関数の自然変数の対応は、紙に手書きで一覧表を作り、机に貼っておくくらいで丁度よい。

論点別に見る出題と失点

統計力学の論点ごとに、頻出パターンと答案で失点する典型例を独立に整理しておきます。論点が違えば失点の原因も違うため、まとめて覚えるより1つずつ理解した方が頭に残ります。

ミクロカノニカル分布

頻出は、孤立系の状態数 W(E, V, N) の計算とエントロピー S = k ln W の導出、そこから温度 1/T = (∂S/∂E)_(V,N)、圧力 P/T = (∂S/∂V)_(E,N)、化学ポテンシャル −μ/T = (∂S/∂N)_(E,V) を順に取り出す流れです。失点の典型例は3つあります。第一に、状態数 W と状態密度 Ω(E) を混同して、対数を取った後の処理がずれるケース。理想気体の状態数を体積 V^N と運動量空間の球殻体積で書くとき、エネルギー殻の厚みを δE のまま残すか積分して状態数に直すかを答案で明示しないと、エントロピーの定義が一意に定まりません。第二に、Gibbs のパラドックスを意識せずに N! で割らない答案。同種粒子の不可弁別性から出てくる 1/N! を最初の状態数の段階で書いておかないと、エントロピーの示量性が壊れます。第三に、温度や化学ポテンシャルを偏微分で取り出す際、どの変数を固定しているかを明示しない答案。これは熱力学関数の自然変数の理解と直結する論点で、減点されやすい。

カノニカル分布

頻出は、分配関数 Z(T, V, N) = ∑ exp(−β E_n) の計算、自由エネルギー F = −kT ln Z、内部エネルギー U = −(∂ ln Z / ∂β)_(V,N)、比熱 C_V = (∂U/∂T)_V、エントロピー S = −(∂F/∂T)_(V,N) の導出。古典統計では理想気体・調和振動子集団・常磁性体（パウリスピン・ランジュバン）の3類型がほぼ毎年どこかで出ます。失点の典型は、分配関数の対数から内部エネルギーを算出する際の偏微分変数の混乱です。U = −(∂ ln Z / ∂β) と書くか U = kT² (∂ ln Z / ∂T) と書くかで符号と係数が変わるため、最初に β = 1/(kT) を明記してから一貫した変数で押し通さないと、調和振動子集団の問題で内部エネルギーの係数を取り違えます。古典極限での因子 h^(3N) の起源（プランク定数による位相空間の量子化）を答案で1行触れておくのも、上位校では効きます。常磁性体の問題では、磁化 M = −(∂F/∂B)_(T,N) として磁場 B についての偏微分を取る流れを、必ず F を経由して書いてください。直接 ⟨S_z⟩ から平均を取る方法は計算は速いが、磁場と熱力学量の対応関係を採点者に示せないため部分点が薄くなります。

グランドカノニカル分布

頻出は、グランドパーティション関数 Ξ(T, V, μ) = ∑_N e^(βμN) Z_N の計算、グランドポテンシャル Ω = −kT ln Ξ、平均粒子数 ⟨N⟩ = −(∂Ω/∂μ)_(T,V)、粒子数の揺らぎ ⟨(ΔN)²⟩ = kT (∂⟨N⟩/∂μ)_(T,V)、そして圧力 P = −(∂Ω/∂V)_(T,μ) の導出。最大の失点源は、化学ポテンシャル μ の符号取り違えと、開放系で熱浴・粒子浴と系をつなぐ概念の曖昧な扱いです。理想気体のグランドカノニカルでは μ = kT ln(n λ³) と書くか μ = −kT ln(n_Q/n) と書くかで符号が逆転するため、自分が定義した熱的ド・ブロイ波長 λ の式（λ = h / √(2π m kT)）を答案の冒頭で固定してから先に進めることが必須です。フガシティ z = e^(βμ) を導入する問題では、z < 1 がボース系で物理的に課される条件であることを書き添えると、低温極限での処理が一貫します。揺らぎを使う問題では、ゆらぎが kT × 応答（圧縮率や帯磁率）として現れる「揺動応答関係」の枠組みを意識しておくと、設問の見通しが立ちます。

量子統計（ボース・フェルミ）

頻出は、ボース・アインシュタイン分布 ⟨n_k⟩ = 1/(e^(β(ε_k − μ)) − 1)、フェルミ・ディラック分布 ⟨n_k⟩ = 1/(e^(β(ε_k − μ)) + 1) の導出、縮退理想気体の比熱、フェルミエネルギー E_F とゾンマーフェルト展開、自由電子気体のパウリ常磁性。失点の典型は、化学ポテンシャル μ の温度依存性をフェルミエネルギー E_F と区別しないことです。T = 0 で μ = E_F、有限温度ではゾンマーフェルト展開で μ(T) = E_F (1 − (π² / 12)(kT/E_F)² + ...) と書けるところを、混同して E_F そのものを温度の関数にしてしまう答案が後を絶ちません。電子比熱 C_V = (π²/2) N k (kT/E_F) を導出する問題では、ゾンマーフェルト展開を何次まで保持するかを明示せず計算に入ると、係数が合っていてもなぜその係数になるかが示せず減点されます。ボース系では、化学ポテンシャル μ ≤ 0 という符号制約を破ると低温極限で発散が出ることを、最初に1行書いておくのが安全です。

相転移とイジング模型

頻出は、強磁性転移、イジング模型の平均場近似（ワイス分子場近似）、自己無撞着方程式 m = tanh(β(Jzm + h)) からの臨界温度の導出、ボース・アインシュタイン凝縮（BEC）。平均場近似の問題で最も多い失点は、自己無撞着方程式の導出を機械的な式変形だけで済ませ、近似の物理的意味（ある1スピンが見る環境を周囲スピンの平均で置き換える）を答案に1行も書かないケースです。臨界指数の議論に入る問題では、平均場近似での β = 1/2、γ = 1、α = 0、δ = 3 という値が「相互作用の次元によらない値」になってしまうことが平均場近似の限界そのものであり、繰り込み群への動機になるという文脈を1行添えるだけで、上位校では明確に評価が変わります。BEC の問題ではリーマンゼータ関数 ζ(3/2) ≈ 2.612 の数値を答案で示し、転移温度 kT_c = (n / ζ(3/2))^(2/3) × 2π&hbar;²/m を導出する流れを、ボース分布の積分形式から落ち着いて追えるかが見られます。

揺らぎと応答

頻出は、エネルギー揺らぎ ⟨(ΔE)²⟩ = kT² C_V、粒子数揺らぎ ⟨(ΔN)²⟩ = kT × (∂⟨N⟩/∂μ)、磁化揺らぎと帯磁率の対応、揺動散逸定理の入門。失点の典型は、揺らぎが熱力学的応答に比例するという構造を「公式の暗記」として処理し、なぜそうなるかを答案で示せないことです。分配関数の二階微分から揺らぎが出てくる関係（⟨E²⟩ − ⟨E⟩² = ∂²ln Z / ∂β² など）を一度導出しておけば、設問が比熱を聞いてきても帯磁率を聞いてきても、答案の骨格は同じです。上位校では、揺らぎが N^(−1/2) で抑えられる熱力学極限の意味を1行書かせる小問が混じることがあるので、定性的な感覚としても押さえておいてください。