解析力学から様々な理論分野を組み立てる 第3部 第1章 §3.1.5 ★★ −∫ρlnρ は Liouville の流れで厳密に一定。分布は細い筋に引き伸ばされ、粗く見たエントロピーだけが増える
Gibbs のエントロピーは一定
相空間の密度 ρ に対して、Gibbs のエントロピーを
SG[ρ]=−kB∫ρln(hfρ)d2fz(3.1.18)
で定義する。hfρ は「升目1つあたりの確率」で、対数の中を無次元にするために hf を掛けた。離散の分布の −kB∑ipilnpi(記法の正本の Gibbs のエントロピー)を連続の相空間に移したものである。このエントロピーは、どんな H のもとでも時間によらない。実際、任意の関数 G について
∫G(ρ(z,t))d2fz=∫G(ρ0(φ−t(z)))d2fz=∫G(ρ0(z′))d2fz′(3.1.19)
で、2つ目の等号では z′=φ−t(z) と変数を変え、Liouville の定理 d2fz=d2fz′ を使った。G(ρ)=−kBρln(hfρ) と置けば SG は一定である。3.1.3節で見たとおり、密度の値の分布が変わらないので、密度の値だけで決まる量はすべて変わらない。
一方、孤立した気体を仕切りの片側に入れて仕切りを外すと、気体は広がってもとに戻らず、熱力学のエントロピーは増える。この2つはどう両立するのか。答えの手がかりは、例題 3.1.3 で見た「相空間の密度は薄まらないが、射影すると薄まって見える」ことにある。
分布は細い筋になる
もっとも簡単な例として、円周の上を自由に回る回転子の集まりを考える。角度を q(2π で同じ点に戻る)、角運動量を p とし、慣性モーメントを 1 にとって H=p2/2 とする。流れは φt(q,p)=(q+ptmod2π, p) で、ずりなので面積を保つ。時刻 0 に、0≤q≤1, 1≤p≤1.5 の長方形に一様な密度 ρ0=1/0.5=2 を置く。角運動量の大きい回転子ほど速く回るので、長方形は傾いた平行四辺形になり、円周を何度も巻き付く(図 3.1.2)。t=16 では、平行四辺形が4本の細い筋に分かれて q 方向の全体に広がっている。筋の1本の幅は時間に反比例して細くなり、筋の数は時間に比例して増える。面積は 0.5 のままで、密度も筋の中では 2 のままである。

図 3.1.2 回転子 H=p2/2 の集まり。0≤q≤1, 1≤p≤1.5 の長方形(左)が、t=4(中)・t=16(右)に運ばれた先。q は 2π で同じ点に戻るので、はみ出した部分は反対の端から入り直す。面積はどの時刻でも 0.5。
位置の分布は一様に近づく
この筋を p について積分し、角度の分布 n(q,t)=∫ρdp だけを見る。n の Fourier 係数 ck(t)=∫02πn(q,t)e−ikqdq(k=1,2,…)は、回転子1つ1つの e−ikq の平均なので、
ck(t)=∫ρ0(q0,p)e−ik(q0+pt)dq0dp=−ike−ik−1⋅−ikΔpte−ikp2t−e−ikp1t,p1=1, p2=1.5, Δp=0.5(3.1.20)
となる。第2の因子の絶対値は 2/(kΔpt) 以下なので、t→∞ で ck(t)→0 となる。k≥1 の Fourier 係数がすべて 0 に近づくことは、n(q,t) が一様分布 1/(2π) に近づくことを意味する。初期の密度が長方形でなくても、p について滑らかに広がっていれば同じ結論になる(Riemann–Lebesgue の補題)。
相空間の密度 ρ そのものは一様に近づかない(筋の中では 2、外では 0 のまま)。一様に近づくのは、q の小さな区間に入る確率のような、粗い量だけである。このように、細かく見た分布は変わらないまま、粗い量だけが平衡の値に近づくことを混合と呼ぶ。混合は運動の性質で、Liouville の定理とは独立した性質である。調和振動子のようにすべての点が同じ周期で回る系では、分布はただ回るだけで筋にならず、混合しない(章末の総合問題)。
粗く見たエントロピー
q と p の相関を忘れ、ρ を周辺分布の積 ρˉ(q,p,t)=n(q,t)g(p) に置き換える操作を考える(g(p) は p の分布で、回転子では時間によらない)。この ρˉ の Gibbs のエントロピー SG[ρˉ] を粗視化したエントロピーと呼ぶ。対数が和に分かれるので SG[ρˉ] は q の部分と p の部分の和になり、p の部分は一定なので、増えるのは q の部分 −kB∫nlnndq だけである。時刻 0 では ρ0 がもともと積の形なので SG[ρˉ]=SG[ρ] で、そこからの増え ΔS を図 3.1.3 に示す。t=2 で ΔS/kB=0.50、t=4 で 0.94、t=16 で 1.80 と増え、上限 ln(2π/1)=1.84(n が幅 1 の一様分布から 2π の一様分布になったときの増え)に近づく。Δpt がちょうど 2π の倍数のとき(t=4π≈12.6 など)は n が厳密に一様になって上限に届き、その後わずかに下がってまた戻るので、増え方は単調ではない。細かく見たエントロピー SG[ρ] は、図の点線のとおり 0 のままである。

図 3.1.3 図 3.1.2 の回転子の集まりで、時刻 0 からのエントロピーの増え。実線は q と p の相関を忘れた ρˉ=n(q,t)g(p) のエントロピー(n を厳密に計算した値)、点線は上限 ln2π と、細かく見た Gibbs のエントロピー(増えは 0)。
増えた分 SG[ρˉ]−SG[ρ]≥0 は、q と p のあいだの相関の量(情報理論の相互情報量に kB を掛けたもの)に等しい。流れのもとで、初期に「どの q にいたか」という情報は消えずに、q と p の細かい相関の中へ移る。相関を見ない観測者にとって、その情報は失われたように見える。どの粗視化を選べば熱力学のエントロピーになるのか、それが単調に増えるのはどういう条件のときかは、情報の側から 相対エントロピーから情報と物理を組み立てる 第2部第4章 で、衝突を足した1粒子分布の方程式として 第3部第8章(Boltzmann 方程式と H 定理)で扱う。
体積が保たれると、いつか戻ってくる
体積の保存からは、もう1つ意外な結論が出る。
命題 3.1.2(Poincaré の再帰定理(弱い形) ★)
H が時刻を含まず、エネルギーが E 以下の領域 ΩE={H≤E} の体積が有限であるとする。ΩE の中の体積が正の任意の領域 U と任意の時間 τ>0 に対して、U から出発してある時刻 kτ(k≥1 の整数)に再び U に入る点がある。
証明
U, φτ(U), φ2τ(U),… はどれも ΩE の中にあり(H が保存するため)、Liouville の定理からどれも U と同じ体積を持つ。これらが互いにまったく重ならないとすると、K 個の和の体積は K⋅vol(U) で、K を大きくすると ΩE の体積を超えて矛盾する。したがってある j<l で φjτ(U)∩φlτ(U)=∅。両辺に φ−jτ を当てると U∩φ(l−j)τ(U)=∅ となり、k=l−j とすればよい。
より強い形として、U のほとんどすべての点が U に無限回戻ってくることも示されている。仕切りを外した気体も、十分待てばもとの片側に集まった状態の近くへ戻る。ただし戻るまでの時間は、粒子数が 1023 程度の系では宇宙の年齢と比べものにならないほど長い。再帰定理は1つの系の軌道についての主張で、上の回転子の集まりのように周期の異なる多数の系がそろって戻ることまでは保証しない。
例題 3.1.5 回転子の集まりの位置の分布が一様に近づく速さ ★★
オリジナル図 3.1.2 の回転子の集まり(0≤q0≤1, 1≤p≤1.5 に一様)を考える。
- (1) 長方形の像の面積が時間によらず 0.5 であることを、流れの Jacobi 行列から確かめよ。
- (2) 式 (3.1.20) を導き、∣c1(20)∣ の値を求めよ。
- (3) 長方形が円周を少なくとも1周分覆うのは、時刻がいくら以上のときか。そのとき n(q,t)>0 がすべての q で成り立つか。
解答
(1) φt(q,p)=(q+pt, p) の Jacobi 行列は (10t1) で、行列式は 1。mod2π は平行移動を貼り合わせるだけで面積を変えないので、面積は 1×0.5=0.5 のままである。
(2) ρ0=2 を q0∈[0,1], p∈[1,1.5] で積分する。q0 と p の積分が分かれて
ck(t)=2∫01e−ikq0dq0∫11.5e−ikptdp=2⋅−ike−ik−1⋅−ikte−1.5ikt−e−ikt=−ike−ik−1⋅−ikΔpte−ikp2t−e−ikp1t
(2=1/Δp を使った)。絶対値は ∣e−ix−1∣=2∣sin(x/2)∣ から∣ck(t)∣=k2∣sin(k/2)∣⋅kΔpt2∣sin(kΔpt/2)∣。k=1, t=20 では 2sin0.5=0.959、2∣sin5∣/10=0.192 で、∣c1(20)∣=0.184。時刻 0 の値 0.959 の約5分の1で、1/t に比例して減っていく。
(3) 時刻 t の像の q の範囲は [t, 1+1.5t] で、その長さ 1+0.5t が 2π 以上になればよいので、t≥2(2π−1)=10.6。このとき、どの q に対しても像の中にそれを通る点があるので n(q,t)>0 である。ただし n が一様になるのは別の話で、像が円周をちょうど整数周覆うときなどに限られる。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com