tecktom.com

出てくるもの (2):止まって見える分布は保存量の関数

2026年10月15日6分で読めますAI生成教科書『虎の巻』第3部 第1章 相空間と Liouville の定理
目次
  1. 定常分布の条件
  2. エネルギー面の上の一様分布と時間平均
  3. 木下の層

AI 生成・木下が未勉強

解析力学から様々な理論分野を組み立てる 第3部 第1章 §3.1.4 ★★★ 時間によらない分布は {H,ρ}=0\{H,\rho\}=0 を満たすもの。f(H)f(H) の形はすべてそうで、ミクロカノニカル分布もカノニカル分布もこの仲間

定常分布の条件

熱平衡の状態では、マクロに見える量が時間によらない。そこでアンサンブルの密度そのものが時間によらない場合、すなわち ∂ρ/∂t=0\partial\rho/\partial t=0 を考える。Liouville 方程式から、その条件は

{H,ρ}=0(3.1.14)\{H,\rho\}=0\tag{3.1.14}

である。第1部第11章の言葉では、ρ\rho が保存量であることと同じである(第1部第11章)。密度は軌道に沿って一定なので、時間によらない密度は、軌道に沿って一定な関数、つまり保存量でなければならない。したがって、HH が時刻を含まないとき、エネルギーの任意の関数

ρ(q,p)=F(H(q,p))(3.1.15)\rho(q,p)=F\bigl(H(q,p)\bigr)\tag{3.1.15}

は定常である。実際、{H,F(H)}=F′(H){H,H}=0\{H,F(H)\}=F'(H)\{H,H\}=0 である。運動量や角運動量のようなほかの保存量があれば、それらの関数も定常である。

統計力学で使う2つの分布はどちらもこの形をしている。

ρmc∝δ(H(q,p)−E),ρcan=e−βH(q,p)∫e−βH d2fz.(3.1.16)\rho_{\mathrm{mc}}\propto\delta\bigl(H(q,p)-E\bigr),\qquad \rho_{\mathrm{can}}=\frac{\mathrm{e}^{-\beta H(q,p)}}{\int\mathrm{e}^{-\beta H}\,\mathrm{d}^{2f}z}.\tag{3.1.16}

左は全エネルギーが EE の面(エネルギー面)の上だけに重みを置くミクロカノニカル分布、右はカノニカル分布である(第3部第2章・第3部第3章)。ここで分かるのは、Liouville の定理がこれらの分布を「力学と矛盾しない」と保証することまでである。F(H)F(H) のどれを選ぶかは Liouville の定理からは決まらない。エネルギー面の上で一様に重みを置くこと(等重率)を選ぶ根拠と、それがいつ正当化されるかは、次の章で扱う。

エネルギー面の上の一様分布と時間平均

1自由度の周期運動では、エネルギー面は1本の閉曲線(1つの軌道)である。この場合、δ(H−E)\delta(H-E) で重みを付けた平均と、1つの軌道に沿った時間平均とが一致することを確かめることができる。δ(H−E)\delta(H-E) を pp について積分すると、H(q,p)=EH(q,p)=E を満たす p=±pE(q)p=\pm p_E(q) の2点で 1/∣∂H/∂p∣=1/∣q˙∣1/\lvert \partial H/\partial p\rvert =1/\lvert \dot q\rvert の重みが立つので、

∫A(q,p) δ(H−E) dp=∑p=±pE(q)A(q,p)∣q˙∣.(3.1.17)\int A(q,p)\,\delta(H-E)\,\mathrm{d} p=\sum_{p=\pm p_E(q)}\frac{A(q,p)}{\lvert \dot q\rvert } .\tag{3.1.17}

一方、軌道に沿った時間平均では、qq と q+dqq+\mathrm{d} q の間にいる時間が dq/∣q˙∣\mathrm{d} q/\lvert \dot q\rvert である。どちらの重みも 1/∣q˙∣1/\lvert \dot q\rvert なので、2つの平均は一致する。エネルギー面の上で「一様」に見える δ(H−E)\delta(H-E) は、qq の上では一様でなく、遅く動く場所(転回点の近く)ほど重くなる。

自由度が2以上では、エネルギー面は 2f−12f-1 次元で、1本の軌道がそれを埋め尽くすとは限らない。保存量がほかにあれば軌道はもっと狭い面に閉じ込められる。軌道がエネルギー面をくまなく回るかどうか(エルゴード性)は、Liouville の定理とは別の問題である。

例題 3.1.4 調和振動子のエネルギー面の上の分布 ★★

オリジナル調和振動子 H=p2/(2m)+mω2q2/2H=p^2/(2m)+m\omega^2q^2/2 の、エネルギー EE(振幅 A=2E/(mω2)A=\sqrt{2E/(m\omega^2)})の軌道を考える。

  1. (1) 軌道に沿った時間平均で、位置の分布 P(q)P(q) を求めよ。
  2. (2) ⟨q2⟩\langle q^2 \rangle を求め、E/(mω2)E/(m\omega^2) に等しいことを示せ。ミクロカノニカル分布 ρmc∝δ(H−E)\rho_{\mathrm{mc}}\propto\delta(H-E) で計算しても同じになることを確かめよ。
  3. (3) ∣q∣<A/2\lvert q\rvert <A/2 にいる時間の割合を求めよ。
解答

(1) q=Asin⁡ωtq=A\sin\omega t とすると ∣q˙∣=ωA2−q2\lvert \dot q\rvert =\omega\sqrt{A^2-q^2}。1周期 2π/ω2\pi/\omega のうち dq\mathrm{d} q にいる時間は往きと帰りで 2 dq/∣q˙∣2\,\mathrm{d} q/\lvert \dot q\rvert なので、

P(q) dq=2 dq/(ωA2−q2)2π/ω=dqπA2−q2(∣q∣<A).P(q)\,\mathrm{d} q=\frac{2\,\mathrm{d} q/\bigl(\omega\sqrt{A^2-q^2}\bigr)}{2\pi/\omega}=\frac{\mathrm{d} q}{\pi\sqrt{A^2-q^2}}\qquad(\lvert q\rvert <A).

∫−AAP dq=1π[arcsin⁡(q/A)]−AA=1\int_{-A}^{A}P\,\mathrm{d} q=\frac1\pi[\arcsin(q/A)]_{-A}^{A}=1 で規格化されている。転回点 q=±Aq=\pm A で発散するのは、そこで速さが 00 になり、長くとどまるからである。

(2) q=Asin⁡uq=A\sin u と置くと ⟨q2⟩=1π∫−π/2π/2A2sin⁡2u du=A2/2=E/(mω2)\langle q^2 \rangle =\frac1\pi\int_{-\pi/2}^{\pi/2}A^2\sin^2u\,\mathrm{d} u=A^2/2=E/(m\omega^2)。ミクロカノニカル分布では、式 (3.1.17) から qq の重みが 2/∣q˙∣∝1/A2−q22/\lvert \dot q\rvert \propto1/\sqrt{A^2-q^2} になり、規格化すると (1) と同じ P(q)P(q) になるので、⟨q2⟩\langle q^2 \rangle も同じ E/(mω2)E/(m\omega^2) である。結果の 12mω2⟨q2⟩=E/2\frac12m\omega^2\langle q^2 \rangle =E/2 は、エネルギーの半分が平均として位置エネルギーにあることを表す。第6章の等分配則(第3部第6章)の1自由度の形である。

(3) ∫−A/2A/2P dq=2πarcsin⁡12=2π⋅π6=13\int_{-A/2}^{A/2}P\,\mathrm{d} q=\frac2\pi\arcsin\frac12=\frac2\pi\cdot\frac\pi6=\frac13。qq の範囲の半分を占める中央部にいる時間は3分の1しかなく、残りの3分の2は外側の半分にいる。

木下の勉強ステータス

木下はまだこの記事を勉強していません。

AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com

← AI生成教科書『虎の巻』の記事一覧へ