解析力学から様々な理論分野を組み立てる 第3部 第1章 §3.1.4 ★★★ 時間によらない分布は {H,ρ}=0 を満たすもの。f(H) の形はすべてそうで、ミクロカノニカル分布もカノニカル分布もこの仲間
定常分布の条件
熱平衡の状態では、マクロに見える量が時間によらない。そこでアンサンブルの密度そのものが時間によらない場合、すなわち ∂ρ/∂t=0 を考える。Liouville 方程式から、その条件は
{H,ρ}=0(3.1.14)
である。第1部第11章の言葉では、ρ が保存量であることと同じである(第1部第11章)。密度は軌道に沿って一定なので、時間によらない密度は、軌道に沿って一定な関数、つまり保存量でなければならない。したがって、H が時刻を含まないとき、エネルギーの任意の関数
ρ(q,p)=F(H(q,p))(3.1.15)
は定常である。実際、{H,F(H)}=F′(H){H,H}=0 である。運動量や角運動量のようなほかの保存量があれば、それらの関数も定常である。
統計力学で使う2つの分布はどちらもこの形をしている。
ρmc∝δ(H(q,p)−E),ρcan=∫e−βHd2fze−βH(q,p).(3.1.16)
左は全エネルギーが E の面(エネルギー面)の上だけに重みを置くミクロカノニカル分布、右はカノニカル分布である(第3部第2章・第3部第3章)。ここで分かるのは、Liouville の定理がこれらの分布を「力学と矛盾しない」と保証することまでである。F(H) のどれを選ぶかは Liouville の定理からは決まらない。エネルギー面の上で一様に重みを置くこと(等重率)を選ぶ根拠と、それがいつ正当化されるかは、次の章で扱う。
エネルギー面の上の一様分布と時間平均
1自由度の周期運動では、エネルギー面は1本の閉曲線(1つの軌道)である。この場合、δ(H−E) で重みを付けた平均と、1つの軌道に沿った時間平均とが一致することを確かめることができる。δ(H−E) を p について積分すると、H(q,p)=E を満たす p=±pE(q) の2点で 1/∣∂H/∂p∣=1/∣q˙∣ の重みが立つので、
∫A(q,p)δ(H−E)dp=p=±pE(q)∑∣q˙∣A(q,p).(3.1.17)
一方、軌道に沿った時間平均では、q と q+dq の間にいる時間が dq/∣q˙∣ である。どちらの重みも 1/∣q˙∣ なので、2つの平均は一致する。エネルギー面の上で「一様」に見える δ(H−E) は、q の上では一様でなく、遅く動く場所(転回点の近く)ほど重くなる。
自由度が2以上では、エネルギー面は 2f−1 次元で、1本の軌道がそれを埋め尽くすとは限らない。保存量がほかにあれば軌道はもっと狭い面に閉じ込められる。軌道がエネルギー面をくまなく回るかどうか(エルゴード性)は、Liouville の定理とは別の問題である。
例題 3.1.4 調和振動子のエネルギー面の上の分布 ★★
オリジナル調和振動子 H=p2/(2m)+mω2q2/2 の、エネルギー E(振幅 A=2E/(mω2))の軌道を考える。
- (1) 軌道に沿った時間平均で、位置の分布 P(q) を求めよ。
- (2) ⟨q2⟩ を求め、E/(mω2) に等しいことを示せ。ミクロカノニカル分布 ρmc∝δ(H−E) で計算しても同じになることを確かめよ。
- (3) ∣q∣<A/2 にいる時間の割合を求めよ。
解答
(1) q=Asinωt とすると ∣q˙∣=ωA2−q2。1周期 2π/ω のうち dq にいる時間は往きと帰りで 2dq/∣q˙∣ なので、
P(q)dq=2π/ω2dq/(ωA2−q2)=πA2−q2dq(∣q∣<A).
∫−AAPdq=π1[arcsin(q/A)]−AA=1 で規格化されている。転回点 q=±A で発散するのは、そこで速さが 0 になり、長くとどまるからである。
(2) q=Asinu と置くと ⟨q2⟩=π1∫−π/2π/2A2sin2udu=A2/2=E/(mω2)。ミクロカノニカル分布では、式 (3.1.17) から q の重みが 2/∣q˙∣∝1/A2−q2 になり、規格化すると (1) と同じ P(q) になるので、⟨q2⟩ も同じ E/(mω2) である。結果の 21mω2⟨q2⟩=E/2 は、エネルギーの半分が平均として位置エネルギーにあることを表す。第6章の等分配則(第3部第6章)の1自由度の形である。
(3) ∫−A/2A/2Pdq=π2arcsin21=π2⋅6π=31。q の範囲の半分を占める中央部にいる時間は3分の1しかなく、残りの3分の2は外側の半分にいる。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com