解析力学から様々な理論分野を組み立てる 第3部 第2章 §3.2.1 ★★★ 分かっているのは E,V,N だけで、相空間のどの点にいるかは知らない。その知らなさを分布 ρ で表す
統計力学の出発点として、多くの教科書は「孤立系では、エネルギーが等しい微視的状態はすべて同じ確率で実現する」という等重率の原理を置く。この一文は短いが、読み返すと分からないことが2つ残る。1つは、数えている「微視的状態」とは何か、という問いである。古典力学の状態は相空間の点で、点は連続に分布しているので、そのままでは数えられない。もう1つは、なぜ「同じ確率」なのか、という問いである。同じ確率といっても、どの変数で測って一様なのか(位置と速度で一様なのか、位置と運動量で一様なのか)で答えが変わる。
この章では、選ぶものは前の章と同じで、N 個の粒子からなる孤立系の Hamiltonian H(Γ) である。当てる道具は前の章(第3部第1章)の Liouville の定理と、第1部の作用変数(第1部第14章)である。Liouville の定理は、相空間の体積 dΓ が時間発展で保たれることを言う。そのため「相空間の体積で測って一様」という分布だけが、時間がたっても形を変えない。等重率の「何で測って一様か」への答えはこれで決まる。数える単位は、量子化の条件から1状態あたり h の面積と決まり、同じ種類の粒子を入れ替えた状態を数え直さないために N! で割る。こうして数えた体積の対数がエントロピーで、そのエネルギー微分が温度の逆数になる。
3.2.1節で孤立系の H と、この章で定義する分布(ミクロカノニカル分布)を書く。3.2.2節では Liouville の定理から、エネルギーの殻の上で一様な分布が時間で変わらないこと、自由度1ではその平均が時間平均と厳密に一致することを示す。3.2.3節から3.2.5節が「出てくるもの」で、理想気体の状態の数え方(h と N!)、2つの系を触れさせたときの温度、殻で数えても体積で数えても結果が変わらない理由を順に扱う。3.2.6節で量子の状態を数える授業の形と並べる。
この章では温度 T を使うので、運動エネルギーには T の字を使わず ∑i∣pi∣2/(2m) と書き、周期は tp と書く。ポテンシャルは体積 V と区別するため、必ず引数をつけて V(r1,…,rN) と書く。
孤立系の Hamiltonian
体積 V の容器に閉じ込めた N 個の同じ粒子(質量 m)を考える。Hamiltonian は
H(Γ)=i=1∑N2m∣pi∣2+V(r1,…,rN)(3.2.1)
で、Γ=(r1,…,rN,p1,…,pN) は 6N 次元の相空間の点である。ポテンシャル V(r1,…,rN) には粒子どうしの相互作用と、容器の壁(壁の外では無限大になる)を含める。系は外とエネルギーをやり取りしない(孤立系)ので、H は時刻を含まず、Hamilton 方程式の解の上でH(Γ(t))=E が保たれる(第1部第8章)。相空間の点は、エネルギー E の等エネルギー面 H(Γ)=E の上を動き続ける。
分かっていることと、知らないこと
実験で決められるのは、全エネルギー E、体積 V、粒子数 N といった少数の量だけである。1023 個の粒子の位置と運動量、つまり相空間のどの点 Γ に系がいるかは分からないし、分かっても使いようがない。そこで、どの点にいるかを確率で表す。相空間の上の確率密度 ρ(Γ) を与え、測りたい量 A(Γ) の値を平均
⟨A⟩=∫A(Γ)ρ(Γ)dΓ(3.2.2)
で予言する。この ρ を統計集団(アンサンブル)と呼ぶ。
ここで問題は、どの ρ を選ぶかである。孤立系について分かっているのは「エネルギーが E である」ことだけなので、ρ はエネルギー E の等エネルギー面の外では 0 にする。面そのものは 6N 次元の空間の中の 6N−1 次元の面で体積が 0 なので、厚み ΔE を持たせた殻 E<H(Γ)≤E+ΔE を考え、後で ΔE の大きさは結果に関係しないことを確かめる(3.2.5節)。殻の中でどう分布させるかについて、この章の答えは「相空間の体積 dΓ で測って一様」である。
定義 3.2.1(ミクロカノニカル分布と状態数)
相空間の測度を
dΓ=h3NN!d3Nqd3Np(3.2.3)
とする(h は Planck 定数。h と N! を入れる理由は3.2.3節で述べる)。エネルギーが E 以下の相空間の体積を
Φ(E)=∫H(Γ)≤EdΓ,(3.2.4)
厚み ΔE の殻の体積を状態数 W(E)=Φ(E+ΔE)−Φ(E) と呼ぶ。殻の中で一様な分布
ρmc(Γ)={1/W(E)0(E<H(Γ)≤E+ΔE)(それ以外)(3.2.5)
をミクロカノニカル分布と呼び、状態数の対数に kB を掛けた
S=kBlnW(E)(3.2.6)
をエントロピーと呼ぶ。
式 (3.2.3) の dΓ は次元のない量なので、Φ(E) と W(E) はただの数になり、対数を取ることができる。「殻の中で一様」という仮定が等重率の原理の中身で、エントロピーの式 (3.2.6) は Boltzmann の名で呼ばれる。この時点ではどちらも定義と仮定にすぎない。S が熱力学のエントロピーと同じふるまいをすることは3.2.4節で確かめる。
測定は時間平均である
ρ を導入する理由を、もう1つの面から見ておく。実際の測定(圧力計の読み、温度計の読み)は、粒子の運動の時間尺度(衝突の間隔は気体で 10−10s 程度)よりずっと長い時間にわたる平均である。測っているのは、1つの軌道に沿った時間平均
A=tf→∞limtf1∫0tfA(Γ(t))dt(3.2.7)
に近い。時間平均を計算するには運動方程式を解かなければならないが、1023 個の粒子では解けない。そこで、時間平均を殻の上の平均⟨A⟩mc=∫AρmcdΓ で置き換える。これが統計力学の基本の手続きである。この置き換えが正しいかどうかは力学の問いで、自由度1なら厳密に正しい(3.2.2節)。まず自由度1の例で、両者が一致する様子を手で確かめる。
例題 3.2.1 4次ポテンシャルの時間平均と殻の平均 ★★
オリジナル1次元の粒子 H=2mp2+aq4(a>0)がエネルギー E で往復している。転回点を Q=(E/a)1/4 とする。
- (1) q が dq を通る時間 dt=dq/∣q˙∣ を重みにして、q2 の時間平均 q2 を u=q/Q の積分で表せ。
- (2) 殻 E<H≤E+ΔE の中で (q,p) の面積について一様な分布を考える。q を止めたときの殻の p 方向の幅を ΔE の1次まで求め、q の分布の重みが (1) の dt に比例することを示せ。これから ⟨q2⟩mc=q2 を結論せよ(ΔE→0)。
- (3) 積分をベータ関数で表し、⟨q2⟩mc を (E/a)1/2 の倍数で出せ。調和振動子の A2/2(A は転回点)と比べよ。
解答
(1) エネルギーの保存から ∣q˙∣=2(E−aq4)/m である。1周期の間に q の各値を往きと帰りの2回通り、往きと帰りの時間は等しいので、時間平均は q>0 の片道 [0,Q] の重み dq/∣q˙∣ で取ればよい。q=Qu と置くと E−aq4=E(1−u4) なので
q2=∫0Qdq/∣q˙∣∫0Qq2dq/∣q˙∣=Q2I0I2,Ik=∫011−u4ukdu
となる。分母と分子に共通の m/(2E)Q は約分で消えた。
(2) q を止めると、殻の中の p は 2m(E−aq4)<∣p∣≤2m(E+ΔE−aq4) の範囲にある。ΔE の1次まで展開すると、正の側の幅は dp/dE⋅ΔE=mΔE/∣p∣=ΔE/∣q˙∣ で、負の側も同じ幅がある。殻の面積を q で切った重みは
2ΔE∣q˙∣dq ∝ dt
で、(1) の重みと定数倍しか違わない。q2 は p によらないので、殻の一様な分布での平均は q の分布だけで決まり、⟨q2⟩mc=q2 である。重みが dt になるのは ∣p∣=m∣q˙∣ で、殻の幅 ΔE/∣q˙∣ が「遅いところほど厚い」ためである。この形は V の形によらない(3.2.2節で一般の形にする)。
(3) u4=s と置くと Ik=41∫01s(k+1)/4−1(1−s)−1/2ds=41B(4k+1,21) である。Γ(45)=41Γ(41) と Γ(41)Γ(43)=π2 を使うと
I0I2=Γ(45)Γ(41)Γ(43)2=Γ(41)24Γ(43)2=Γ(41)48π2=0.457,⟨q2⟩mc=q2=0.457(aE)1/2
を得る(I2=0.599、I0=1.311)。転回点で測ると q2=0.457Q2 で、調和振動子の 0.5A2 より小さい。4次の井戸は底が平らで壁が急なので、粒子は転回点の近くで速さを落とす区間が短く、中央に長くいる。値の E1/2 への比例は、q の大きさが Q∝E1/4 で決まることだけから出る。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com