104. 熱力学・統計力学

2019年8月 3日 (土)

ボルツマン方程式と緩和時間

表題のテーマについて詳細に説明すると約束して

いましたが,これについても本ブログの2008年に

関連の過去記事があったので,今回はこれらの記事

を再掲載することから始めます。

 

※以下,まず,再掲記事第1弾:2008年11/2にアップ

した古典統計力学においての考察記事です。

 

表題「ボルツマン方程式とH定理」

 

今日は,不可逆過程の要因解明の関連で.

「ボルツマン(Boltzmann)方程式」と

「ボルツマンのH定理」について述べたい

と思います。

 

 まず,質量がmで全分子数がNの気体が速度

v+dvの間にある粒子数の分布をf()d

します。

 

簡単な考察によって,絶対温度Tで熱平衡状態に

ある場合,f()はMaxwell-Boltzmann分布,つまり,,

f()=N[m/(2πkBT)]3/2exp[-m2/(2kBT)]

に従うことがわかります。

 

これは,∫f()d=Nと規格化しています。

 

次に,同じ気体分子が非平衡状態にあるとして,

その分布関数を位置と速度,および,時刻tの

関数としてf(r,,t)と表記します。

 

つまり,時刻tに+d,+d

ある分子数をf(,,t)dとします。

これも,∫f(,,t)d=Nと規格化して

おきます。

 

 このとき粒子の衝突を無視した自由運動による

各位置の近傍での粒子数の保存を示す連続の方程式

は,∂f/∂t+∇f=0 となります。

 

これはリウビル(Liouville)の方程式を分布関数で

与えたものとなっています。

 

 しかし一般に非平衡状態では衝突による粒子数変化

による「湧き出し項」として.衝突項が存在して連続

の方程式は∂f/∂t+∇f=(∂f/∂t)coll

なるはずです。

これを「ボルツマン(Boltzmann)(輸送)方程式」

と呼びます。

 

そして,ある時刻tに速度’と1’をもつ粒子対が

衝突して単位時間に速度1との粒子対となって,

+d,11+d1領域に

入ってくるプロセスの頻度を,

σ(,1|’,1’) とします。

 

これと全く逆に,+d,11+d1領域

から出て行くプロセスの頻度を

σ(’,1’|,1) とします。

 

すると,衝突(湧き出し)項は,

(∂f/∂t)coll

=∫σ(,1|’,1’)f(r,’,t)

×f(,1’,t)d1’d1

-∫σ(’,1’|,1)f(,,t)

×f(,1,t)d1’d1

になると考えられます。

 

 ところで力学の時間反転に対する対称性

から,1から’と1’に変わる頻度は,

’と-1’から-と-1に変化する頻度

に等しい。ということがいえます。

 

つまり,σ(’,1’|,1)

=σ(-,-1|-’,-1’)

と考えられます。

この等式を「衝突数算定の仮定」といいます。

 

 さらに座標軸の向きを逆転させても,こうした

プロセスの頻度は同じと考えられるので,

σ(-,-1|-’,-1’)

=σ(,1|’,1’) 

も成立するはずです。

 

そこで,結局

σ(,1|’,1’)=σ(’,1’|,1)

としてよいと考えられます。

 

 ここで略記法として,

f=f(r,,t), f’= f(,’,t),

1=f(,1,t),f1’≡ f(,1’,t)

と置くことにすれば,

衝突項は,(∂f/∂t)coll

=∫σ(,1|’,1)( f’f1’- f f1)

1’d1’と簡単になります。

 

気体分子の衝突は,弾性衝突と考えてよいので,

衝突1の前後でエネルギーも運動量も保存されると

考えられるため,1’+1’かつ,

2+v12=v’2+v12 以外の場合には.

σ(,1|’,1)=0 です。

 

 Boltzmann方程式が「不可逆過程」を記述している

ことを示すために,ここで「ボルツマンのH関数」

という関数:Hを.

H(,t)=∫f(r,,t)logf(r,,t) d

で定義します。

ここでlog=ln (自然対数)lnを意味します。

 

 このとき,∂H/∂t=∫(∂f/∂t)(logf+1)d

と書けますが,これにボルツマン方程式:

∂f/∂t+∇f=(∂f/∂t)coll を代入すると,

∂H/∂t=∫(logf+1)[-∇f+(∂f/∂t)coll ]

=-∇[∫∇(flogf)d]

+∫(logf+1)(∂f/∂t)coll

となります。

 

この右辺で,1×(衝突項)の部分の積分は,

∫(∂f/∂t)coll

=∫σ(,1|’,1)( f’f1’- f f1)

1’d1’です。

 

これは,σ(,1|’,1)の,

,1,’,1’の粒子交換に対する対称性

と.( f’f1’- f f’)の交換反対称性から,

ゼロとなります。

 

したがって,Hの流れとして

H∇(flogf)dを定義すると,

∂H/∂t+∇H

=∫[(logf)(∂f/∂t)coll]d

となります。

 

これは,ボルツマンのH関数の.流出入以外の

正味の生成である,dH/dt=∂H/∂t+∇H

が,∫[(logf)(∂f/∂t)coll]dで与えられる

ことを示しています。

 

そして,∫[(logf)(∂f/∂t)coll]d

=∫[(logf)σ(,1|’,1)( f’f1’- ff1)]

1’d1’です。

 

この式の右辺は,

∫[(logf1)σ(1,|1’,’)

×( f1’f’- f1 f)]d1’d1’,

∫[(logf’) σ(’,1’|,1)

×( f f1-f’f1’)]d1’d1’,

∫[(logf1’)σ(1’,’|1,)

×( f1 f-f1’f’)]]d1’d1

の全てと等しいことになります。

 

しかも対称性からこれらのσは全て等しいので

簡略してσ(,1|’,1’)を単にσと

略記することにします。

 

すると,∫[(logf)(∂f/∂t)coll]d

=(1/4)∫[σ( f’f1’- ff1)

×(logf+logf1-logf’-logf1’)]

1’d1

=(1/4)∫[σ( f’f1’- ff1)log(ff1/f’f1’)

1’d1

と書けることになります。

 

ところが,σは衝突頻度ですから,当然σ≧0 であり

しかも,(x-y)log(y/x)≦0 ですから,結局,

dH/dt=∫[(logf)(∂f/∂t)coll]d≦0

が示されたことになります。

 

つまりボルツマンのH関数は時間と共に,常に一定.

または,減少する,ということが示されたわけです。

 

こうして,時間反転不変=可逆な力学法則から,

どういうわけか不可逆変化が導かれました。

これを,「ボルツマンのH定理」といいます。

 

しかし,これに対しては,

「ロシュミット(Loschmidt=の逆行性批判」という

有名な反論があります。

 

すなわち,「ある瞬間に時間的変化を反転する,

つまり全粒子の向きを逆転させると逆にHは過去

に向かって減少する,または未来に向かっては増加

することになる。」という反論です。

これは,まことにもっともな話です。

 

こうしたさまざまな反論に悩んだ末に,とうとう

ボルツマンは自殺に追い込まれてしまったのです。

 

今考えると,実はH定理は確率法則による定理で

あり,例えば衝突頻度σに対して「衝突数算定の仮定」

が導入されています。

 

既に,「速度空間の大きい体積の方には,小さい体積

よりも粒子数が多いはずである。」などの等重率の原理

のような確率的構想が入っていて,単純な可逆的力学

法則からの確率概念的な飛躍があることに気づきます。

 

というわけで,確率法則としてボルツマンのH定理

の主張は正当である,と認めて問題ないと思います。

 

ところで,H関数は非平衡状態に対して与えられたもの

ですが,熱平衡状態は∫[(logf)(∂f/∂t)coll]d=0 ,

つまり,ff1=f’f1’であること,に対応します。

 

そして,このときはfをで積分したものが,速度vに

対するMaxwell-Boltzmann 分布を与えます。

 

平衡統計力学においてのみ定義され,数式的に

与えられるエントロピーSを計算すると,

これは,ボルツマンのH関数と.

S=-kB∫H(,t)d+(定数) という関係

にあることがわかります。

 

そこで、エントロピーの概念を拡張して,非平衡状態

でもエントロピーSをS=-kB∫H(,t)dで定義

すればよいのでは? と考えることができます。

 

「孤立系=(構成粒子の流出入のない系)では.

エントロピーは常に増加する。」という熱力学第2法則

は,平衡状態に対するボルツマンのH定理の言い換えに

過ぎない,ということにもなります。

 

(参考文献):

北原和夫 著「非平衡系の統計力学」

(岩波基礎物理学シリーズ)

エリデ・ランダウ 著「統計物理学」(岩波書店)

テル・ハール「熱統計学Ⅰ」(みすず書房)

 

(以上,再掲記事(1)終わり※)

 

(※次に,再掲記事第2弾 2008年11/8アップ)

表題「量子的ボルツマン方程式」

 

「古典的気体分子」ではなくて「金属内の自由電子」

という量子的”フェルミ気体“にボルツマン方程式;

∂f/∂t+∇f=(∂f/∂t)coll を適用すること

を考えてみます。

 

c=h/(2π)とすると,波数の電子波束の速度

=(群速度)は,k=(1/hc)∂εk/∂であり.の時間変化

は,外力として電場があるときには電子の電荷をe(<0)

として,d/dt=ehhcとなります。

 

そして,波数がの電子波束をフェルミ粒子として,

その分布関数をfkとすると,

「量子的ボルツマン方程式」は,

∂fk/∂t+k∇fk(∂fk/∂)=(∂fk/∂t)coll

となります。この左辺の第3項は,電場による波束の変化、

つまり,運動量;=hcの時間的変化による分布の変化

を示しています。

 

波数:’+d’と1’+d1’の間の電子波束対

が衝突して,単位時間に波数:1の波束対となり,

+d,11+d1領域に入ってくる散乱

の確率を,σ(,1|’,1’)d1’d1

とし,これと全く逆に+d, 11+d1

領域から出て行く散乱の確率を

σ(’,1’|,1)d1’d1

とします。

 

気体分子運動論のときと同様に,f=fk ,f1=fk1,

f’=fk’,f1’=fk1’と置くと,衝突前の波数:,

1の波束の減少は(ff1)にも比例しますが,衝突後

の状態:’と1’が占有されていないことが必要

なので,(1-f’)(1-f1’)にも比例します。

 

したがって,”衝突項”は,

(∂f/∂t)coll=∫σ(,1|’,1’)f’f1

×(1-f)(1-f1) d1’d1

-∫σ(’,1’|,1)ff1(1-f’)(1-f1’)

1’d1’ 

となる,と考えられます。

 

量子力学でも時間に関する可逆性,空間反転対称性

は成り立ち「衝突数算定の仮定」の統計的意味は,

量子力学自体が確率現象なので,さらなる重みを

持ちます。

 

それ故,(∂f/∂t)coll=-∫σ(,1|’,1’)

[ff1(1-f’)(1-f1’)

-f’f1’(1-f)(1-f1)]d1’d1’ 

となるはずです。

 

量子論でのフェルミ粒子系のエントロピーは.

S=-kB∫[flogf+(1-f)log(1-f)]d

ですから,このときのH関数は.

H=∫[flogf+(1-f)log(1-f)]d

と定義すればよい,と思われます。

 

古典論と同じように考察して,

dH/dt=∫[log{f/(1-f)}(∂f/∂t)coll]dk

=(1/4)∫σ[f’f1‘(1-f)(1-f1)

-ff1(1-f’)(1-f1’)]log[ff1 (1-f’)(1-f1’)

/{f’f1’(1-f)(1-f1)}]d1’d1

と書けます。

 

そこで古典的気体分子運動論と同じく,金属内の

「フェルミ気体」である自由電子の運動論でも,H定理:.

dH/dt=∫[log{f/(1-f)}(∂f/∂t)coll]dk≦0

が成立して,ぉれがエントロピーの増大則:dS/dt≧0

を保証します。

 

今度の場合では,熱平衡状態というのは,

∫[(logf)(∂f/∂t)coll]d=0 ,すなわち,

ff1/[(1-f)(1-f1)]

=f'f1'/[(1-f’)(1-f1’)]に対応しており,

fがフェルミ・ディラック分布:

f=1/[exp{(εk-μ)/(kBT)}+1]

に従うならば,この条件は確かに満足されて

います。

 

(参考文献):

北原和夫 著「非平衡系の統計力学」

(岩波基礎物理学シリーズ)

 

(以上,再掲記事(2)終わり※)

 

さて,「ボルツマン方程式」を古典論と量子論の両方

で定義し,考察した過去記事を上に紹介しましたが,

これらの主眼は,「ボルツマンのH定理」の数学的証明

を与えて解釈することにあり,「時間反転不変で可逆的

な1粒子の基本運動方程式から,多体系に移って統計的

に扱うときには,如何にして時間を逆行することができ

ないという巨視的意味での「不可逆性」が導入され,

正当化されるか?を述べた記事たちです。

 

では,最初,固体系内の電子が陽イオン芯であれ,フォノン

であれ,衝突して散乱される際,その衝突間に自由な運動を

する平均時間として導入され,定義された「緩和時間」と

いうのは,ボルツマン方程式の中ではどこに現われている

のでしょうか?

 

fを位置と速度,時間tの関数としての分布関数

とすると,素朴なボルツマン方程式は,

∂f/∂t+∇f=(∂f/∂t)coll ,

あるいは,∂f/∂t+∇f-(∂f/∂t)coll =0

で与えられるものでした。。

 

これに緩和時間の概念が入るとしたら当然,

衝突項:(∂f/∂t)coll の中でしょう。

 

緩和時間τは,たまたま,手にとった参考書によれば.

(∂f/∂t)coll­[f(,v)-f0]/τ

なる式で定義されるとあります。

ただしf0は電場などの外力が無いときの熱平衡状態

の分布関数(古典論ではMaxwell-Boltzmann分布,

量子論ではFermi-Dirac分布)であり,

右辺のf(,v)は電場Eなど外力があるときの,

左辺の分布関数f(r,v,)が,t→大で平衡に

達したときの分布てす。

 

(参考文献):

「量子物理学入門」(東京電機大学出版会),

| | コメント (0)

2019年7月12日 (金)

-熱交換器・ヒートポンプなど

最近,子供のころ疑問に思っていた科学的事象を,今70歳

近くになるまでに知り得た多少の知見,薀蓄で,まとめたもの

をブログ記事に書き留めておこうと思うに到りました。

 

こうした疑問への詳細な回答は,最近はネットにアクセスして

直接質問したり,他人の質問への回答を見たりWikiなどを参照

するなど,いくらでも答を得る手段があるので,今更という感

もありますが,他人を啓蒙しようとする「上から目線」の意図

というより,ある意味自己満足の日記的覚え書きです。

 

さて,今回は冷房,暖房を司るヒートポンプの仕組みに関わる

熱交換器や冷媒についての薀蓄を書こうと思い,それの導入

として,まず,10年以上前にブログを開始して半年くらいの

2006年8/6にアップした記事「接触による2物体の温度交換」

を再掲載するところから始めようと思います。

 

  • 以下は,まず上記の過去記事の丸写し再掲です。を訂正し,途中計算を略してわかりにくいと思った部分を 温度がそれぞれTA,TBであったとし,これらを接触させて放置 にあったと記憶していたものを参考にしています。まず,一連の手順を5つの工程に分けて行います。これらの物体を 
  • 取り囲む環境は断熱で,熱は逃げたり入ったりしないと仮定します。
  •  
  • これは,どこだったか覚えていませんが,ある大学入試の過去問
  • するだけで,AとBの温度を交換する方法を考えてみます。
  • ※同じ質量Mで同じ物質から成る2つの物体(固体:A,Bの
  • 詳細に補足しておきました。
  • ただし,現在読み直して気づいた当時の若干の誤記,間違い
  1. AとBをそれぞれ半分の質量(M/2)の2つの物体:A1,A2,質量が同じですから「熱量保存の法則」により,これらの温度は, 
  2. A1=TB1=(TA+TB)/2 となるはずです。
  3. および,B1,B2に分割します。そしてA1とB1とを接触放置します。
  4. 次にA1とB2とを接触放置します。となるはずです。
  5.  
  6. 温度はTA1=TB2={(TA+TB)/2+TB}/2=(1/4)TA+(3/4)TB
  7. 次にA2とB1とを接触放置します、となるはずです。
  8.  
  9. 温度はTA2=TB1={TA+(TA+TB)/2}/2=(3/4)TA+(1/4)TB
  10. さらにA2とBとを接触放置します。 
  11. 温度はTA2=TB2=(TA+TB)/2 となるはずです。
  12. 最後に、分割していたA1とA2,および,B1とB2を, TA’=(TA1+TA2)/2=(3/8)TA+(5/8)TB,かつ,=(1/2)[{(3/4)TA+(1/4)TB}+(TA+TB)/2]と変わっているはずです。これら一連の手順で,結局,Bの温度はTBからTB’=(5/8)TA+(3/8)TBに当然のことながら,質量が同じで同じ物質なので, であったとすれば,この手順の結果として,TA’=272℃,かつ, すなわち,(1)まず,AとBをそれぞれ4分の1の質量(M/4)します。(2)まず,A1,A2とB1,B2に先の手順を施行します。A,Bの代わりの質量(M/2)の2物体として.上の操作をします。4分割直後の時点ではA1,A2,A3,A4,の温度は全てTAであり,①~⑤を実行して最後に接着したときには,となるはずです。
  13. A1+A2=(3/8)TA+(5/8)TB,B1*B2=(5/8)TA+(3/8)TB
  14. 1,B2,3,B4の温度は全てTBで結果は質量に無関係ですから,
  15.  
  16. つまり,まず,(A1+A2)と(B1+B2)を上記の質量Mの2物体
  17.  
  18. の4つの物体:A1,A2,A3,A4,および,B1,B2,3,B4に分割
  19. ではこのプロセスをもう1段階増やすとどうなるでしょう。
  20. B’=240℃となり,温度の高低は逆転することになります。
  21. ちなみに,最初Aの温度がTA=192℃,Bの温度がTB=320℃
  22. A’+TB’=TA+TBが成立しています。
  23. 変わることになります。
  24. Aの温度はTAからTA’=(3/8)TA+(5/8)TBに,
  25.  
  26. =(5/8)TA+(3/8)TB 
  27. B’= (TB1+TB2)/2
  28. =(1/2)[{(1/4)TA+(3/4)TB}+(TA+TB)/2]
  29. 最後に戻したA,Bの温度は,このとき
  30. それぞれ接着して元のAとBに戻します。

 

(3)続いて,先の操作をさらにA1,A2とB3,B4に施すと

A1+A2=(3/8)TA1+A2+(5/8)TB,かつ,

B3+B4=(5/8)TA1+A2+(3/8)TBとなります。

 

(4)さらにA3, A4とB1, B2で実行すると,

A3+A4=(3/8)TA+(5/8)TB1*B2,かつ,

B1*B2=(5/8)TA+(3/8)TB1*B2 です。

 

(5) さらにA3,A4とB3,B4で実施すると

A3+A4”=(3/8)TA3+A4’+(5/8)TB3+B4,かつ,

B3+B4=(5/8)TA3+A4’+(3/8)TB3+B4です。

 

(6)そして最後にA1,A2,A3,A4を全て接触,接着,

および,B1,B2,3,B4を全て接触,接着させて

このステップは終わりです。

 

この結果,A,Bの最終温度:TA~,TB~は,

A~=(TA1+A2+TA3+A4”)/2

={(3/8)2+(3/8)(5/8)2}TA

+{1-(3/8)2-(3/8)(5/8)2}TB,

および,

B~=(TA+TB)-TA~

={1-(3/8)2-(3/8)(5/8)2}TA

+{(3/8)2+(3/8)(5/8)2}TB

となります。

 

※(注):上式を証明します。

[証明]:簡単のためにa=3/8,b=1-a=5/8とします。

すると,(2)はTA1+A2=aTA+bTB,B1*B2=bTA+aTB

(3)はTA1+A2=aTA1+A2+bTB=a2A+(ab+b)TB,

B3+B4=bTA1+A2+aTB=abTA+(b2+a)TB

(4)はTA3+A4=aTA+bTB1*B2=(a+b2)TA+abTB,

B1*B2=bTA+aTB1*B2=(b+ab)TA+a2Bです。

 

故に(5)は,TA3+A4”=aTA3+A4’+bTB3+B4

=(a2+ab2)TA+a2bTB+ab2A+(b3+ab)TB,

=(a2+2ab2)TA+(a2b+b3+ab)TB,かつ,

B3+B4=bTA3+A4’+aTB3+B4

=(ab+b3)TA+ab2B+a2bTA+(ab2+a2)TB

=(a2b+b3+ab)TA+(a2+2ab2)TBとなります。

 

したがって.TA1+A2+ TA3+A4”

=2(a2+ab2)TA+(a2b+b3+2a+b)TB,

です。ところが,b=1-aなので,

2b+b3+2ab+b=b(a2+b2+2a+1)

=b(2-2ab+2a)=2(1-a)(1-ab+a)

=2(1-ab+a-a+a2b-a2)

=2(1-ab+a2b-a2)=2(1-a2-ab2)

です。

それ故, TA1+A2+ TA3+A4”

=2(a2+ab2)TA+2(1-a2-ab)TB より

A~=(TA1+A2+ TA3+A4”)/2

=(a2+ab2)TA+(1-a2-ab)TB

が得られます。

 

そして,TB~=(TB1+B2+ TB3+B4”)/2ですが

A~+TB~=TA+TBを用いれば

B~=(TA+TB)-TA~~

=(1-a2-ab)TA+(1-a2-ab)TA

を得ます。[証明終わり]  (注終わり※)

 

では,このプロセスをn回繰り返して2個まで分割

し,最後にはn→∞の無限分割まで際限なく繰り返した

場合.極限では,どうなるのでしょうか?

便宜上,n段階の操作後のAの温度をTnとしておきます。

 

すると,ここまでの話から,明らかに,T1=aTA+bTB,

2=(a2+ab2)TA+(1-a2-ab2)TBですから

一般に,Tn=anA+bnBと置くと,bn=1-anであり,

n+1=an2+ann2=an2+an (1-an)2

=an(an2-an+1),bn+1=1-an+1と漸化式で書けます。

 

しかし,an+1=an(an2-an+1)(a1=3/8)は線形な関係式

ではないし,この数列の一般項を陽に導くのは試行しましたが

できませんでした。

 

それでも,an+1/an=an2-an+1=(an-1/2) 2+3/4>3/4

であり,また,an+1/an-1=an2-an=-an(1-an)です。

 

そこで0<an<1なので, 0 <(an+1/an)<1であり,

一般に0<r<1なる定数rが存在して0<(an+1/an)<r

が成立するため,n→∞に対してan=rn-11→ 0 となること

がわかります。したがって,もちろん.limn→∞n → 1です。

 

かくして無限分割の施行では「マクスウェルの悪魔」という

わけではないですが,無限回の接触操作によってAとBの温度

を完全に交換して入れかえることが原理的には可能になります。

(過去記事再掲載終了)※

 

※ここからは今回の新しい記述です。

上述の同一質量Mの固体物質A,Bの同じ単位質量当たり

同一の熱容量=比熱をcとし,TA,TBを絶対温度と考えると,

A,Bの熱量=熱エネルギー:QA,QBは.QA=∫TAcMdT,

B=∫TBcMdTで与えられるため,上記では温度交換と

称していますが実は,熱交換でもあります。

 

比熱cは固体金属の塊りなどでは,常温では剛体と同じく

運動の自由度は6で,振動の自由度や量子効果を無視すれば

古典統計力学の「エネルギー等分配の法則」により1モル

の比熱は,cmol=(R/2)×6=3R~(8.31×3)J/(mol・K)

~(25/4.19)cal/(mol・K)です。例えば鉄(Fe)なら標準の鉄

の分子量は56なので,c~0.11cal/(g・K)です。

これは固体比熱のデュロン・プティ(Dulong-Petit)の法則と

して知られています。

 

さて,上の記事を書いた約13年前にも着目したことですが,

「(何もしなければ)熱は高温部から低温部に向かって移動し,

その逆向きに流れることはない。」という熱力学第2法則に

上記事実は反しているようにも見えますが,実は分解して接触

させるという力学的操作が加えられており「何もしなければ」

とか,「ひとりでに(自然に)」という条件が満足されていない

ので,低温から高温に熱が流れていても,別にこの法則を破って

いるわけではありません。

 

ヒートポンプ(エアコン)や冷蔵庫などでは電気エネルギー

でモーターを回すことで,こうした操作が施行されています。

 

こうした常温の空気を冷やしたり暖めたりする機器には

熱交換器が使用されています。

一般に熱交換器とは性質の異なる2流体が流れている容器

接触させて高温流体から低温流体へと熱を移動させます。

2つの物体を接触させて放置しておくとやがて熱平衡と

呼ばれる平衡状態に到ります。

この熱的平衡にあるかどうか?の指標が温度(経験温度)です。

 

AとBが熱平衡にあることをA~Bと表わすと,これは一種

の同値関係A~B,B~CならA~Cという推移律が

熱力学第0法則と呼ばれています。

そこでBを温度計とみなせば,AとCを接触させなくても

それらの熱平衡の度合いが測れるわけです。

 

熱交換器の2流体として,一方は普通の水,他方は冷媒と

呼ばれ,常温では既に沸点を超えて気体として存在するもの

を,低温の液体状態にして使用します。

元々は,フッ素(F),塩素(Cl).臭素(Br)などのハロゲンを

含むガス状態ではフロンガスと呼ばれるものを冷媒として

多用していましたが,近年,フロンはオゾン層を破壊して

光化学スモッグを発生させたり,有害な紫外線など太陽

から地球に注ぐ放射線(宇宙線)を貫通させて皮膚ガン

を誘発させるなど,さらに温室効果ガスでもあるという

ことで世界的に規制され,何か代替物に取って代かわら

れているらしいです。

 

要するに,常温の温水管が冷媒管に接触する普通に温水

が冷やされ,液体冷媒は水から気化熱を奪って気体(ガス)

になります。

これだけではやがて熱平衡になって冷却という作用は

終わってしまsいますから,通常,熱機関というのはサイクル

と呼ばれるシステムで元の状態に戻し熱の流れを継続する

必要があります。

(何もしないでもサイクルとなる永久機関は無いのです。)

 

気化熱をもらってガスとなった冷媒を再び圧縮して液体に

戻し断熱圧縮で発生する熱は内ではなく外に逃がします。

 

冷蔵庫ではコンプレッサーで圧縮して放熱板で外に逃がし.

ヒートポンプでは室外機です。

具体的には圧縮して熱を持った冷媒を室外で空冷式にファン

で冷却して室内に戻します。

 

日本の都会の夏は,ほとんど全ての建物で冷房が行われ,室内

の熱:トータルでは莫大な熱を屋外に捨てているので,屋外の

気温を上昇させる効果があります。

 

これはヒートアイランド現象とも呼ばれ,田舎に比べ都会は,

地球温暖化による異常気象のほかに,過剰な冷房のためにも

屋外気温が昔より高くなっています。

 

一方,暖房は冷房と逆のプロセスです。冷媒ガスを圧縮して

際の発生熱を室内に与えて空気を暖めます。圧縮で液体と

なった冷媒を室外で冬の冷気からさらに気化熱を奪ってガス

に戻すわけです。いずれにしろ,全体で熱交換機です。

※さて余談ですが万有引力により地球に束縛されている空気中の

酸素:O2は太陽光を吸収してO2+hν→O+Oと分解され,この

酸素原子をO2+O→O3なる反応で吸収してオゾン:O3ができると

考えられていますが,O3は酸素:O2よりエネルギー的に不安定で,

短絡的に書くと,逆反応:O3→O2+O+hνで光を放出して酸素

に戻ります。

 

これらの反応が平衡してオゾン層を形成していれば,太陽から

地上まで直線的に降り注ぐ光の中でも,特に高エネルギーで人体

にも有害な放射線が,この層で一旦吸収されて放出される際には

向きを変えるという意味での散乱により緩和されます。

 

ところが空気中のPAN(ポリアセチルナイトレート)やフロンが

増えると.こうしたものの介在で反応のバランスが崩れ,O3

分解されてO2に戻るという反応の方が過剰になるというのが

オゾン層破壊のシステムと考えられています。

 

日本では,昔,よく晴れた日には光化学スモッグが発生して目が

痛くなったりするということが問題となった時代がありました。

オゾンの平衡が崩れたのが原因ということで対策がなされて

いました。その後,そうした話を聞かなくなったので,今はどう

なっているのかは知らないのですが,これは丁度,私が20,30代

で,まだ環境アセスメントの会社の社員であった頃に,門前小僧

で得た薀蓄です。

 

また,特に南極大陸上空のオゾン層が破壊され,南極に近い南半球

のオーストラリア南部で,メラがニン色素の少ない肌の白人を中心

に日光による赤焼けなど危険で,オゾン層破壊が皮膚ガン増加

に寄与していると聞いています。

(※ 2006年7/12の本ブログ過去記事「オゾンホール」,および,

7/27の「オゾンホール(訂正)」も参照されたい。※)

 

一方,温室効果ガスですが,これは,そもそも太陽から地球表面に

注ぐ総エネルギーからStefan-Boltzmannの法則(T4則)で単純

計算すると地球上平均温度は,ー18℃くらいであるという評価を

過去記事で書きました。

(※2006年11/21の記事:

「地球の平均気温とステファン・ボルツマンの法則」参照※)

実際の地球平均気温が約15℃であるのは,吸収された熱が地上

から放射されて単純に全て宇宙空間に逃げるという計算では

なく,水蒸気の雲などにより反射し返されて実質上地球が温室

の様になるという意味で,これらの気体を温室効果ガスと呼ぶ

わけでです。

 

最大の温室効果ガスは水(水蒸気)ですが米副大統領ゴア氏の

「不都合な真実」以来,二酸化炭素が増えると,さらに温暖化

するとされて温室効果ガスは危険視されており,フロンもこれ

に寄与するらしいです。

 

PS:2月のリニューアル以後,ブログ書きがむずかしくなり,

コピペもうまくできません。

オリジナル発想の遺稿をアップするまでがんばろうと思って

いましたがモチベーションが下がりました。

2006年からぷ愛用していたブログのテノウレートが

リニューアル前のように戻らないなら,もう科学ブログは

やめて,後わずかなとき訪れるだろう永眠に備えようか?

という気になりました。(※実際,上記本文でもワードで

書いた草稿をコピペ゚でアップしただけで改行不具合発生,

文章順不同で訂正する気にもなりません。)

 

 

| | コメント (0)

2014年12月12日 (金)

統計力学の基礎(6)(量子統計力学4)

 統計力学の基礎概念への脱線から離れて,再び,量子統計力学3の

続きです。

 

経路積分と摂動論に戻る最後の関門の第2量子化とGreen関数

について記述します。

 

今日は,まず,第2量子化について述べます。

そのため,最初はBose統計に従う粒子に着目します。

 

粒子数の演算子rを考えます。

添字rは具体的にはBoseスカラー粒子なら波数(運動量)だけ

ですが,Fermi粒子かスピンが1のBose粒子ならスピン変数をσ

としてr=(,σ)などを意味します。

 

便宜上,適宜この添字rを省略して個数演算子を,その固有値

nを持つ状態を,ケットベクトル |n>で記述します。

|n>=n|>です。

 

固有値nは整数値のみを取り,固有ベクトル|n>は

<m|n>=δmn直交規格化されているとします。

 

こうしたの固有状態に対して次のような作用を及ぼす演算子 

,を導入します。

 

すなわち, |n>=n1/2|n-1>,

|n>=(n+1)1/2|n+1> なる作用です。

 

つまり,は粒子数を1つ減らす演算子,は1つ増やす演算子

です。そこでを消滅演算子,を生成演算子と呼びます。

 

これらは状態のケットにのみに作用して,その係数には作用せず

素通りするとします。

 

すると|n>=n1/2|n-1>=n|n>, 

aa|n>=(n+1)1/2|n+1>=(n+1)|n>

が成立します。

 

それ故,が,個数演算子の役割をすることがわかります。

 

そじで添字rを復活させると,rrrです。

 

,の行列要素をつくると,任意のm,nに対して, 

<m||n>=n1/2<m|n-1>=n1/2δm,n-1であり, 

<n||m>=(m+1)1/2δm,n+1=n1/2δm,n-1ですから 

<n||m>=<m||n>* が成立します。

 

したがって,のHermite共役です。

よって,rrのHermite共役です。

 

そこで,上添字+(ダガー)を文字通り演算子のHermite共役を意味

する記号と考えれば,r=(rr)rrrとなるため

rは実数のみを固有値とするHermite演算子なので,この意味では

定義は無矛盾です。 

 

また,|n>=n|n>,aa|n>=(n+1)|n>より,

全ての|n>に対して,(aa|n>=|n>ですから,

aa=1 が成立します。

任意の演算子,に対して交換子[,]を

[,]≡ABBAで定義すれば, 

[,]=1と書けます。

一方,[,]=[,]=0 は自明です。

あるいは,これを一般化して[r,]=δrs,

[r,]=[,]=0 とします。

 

これは,異なる添字r≠sのr,の固有状態なら固有値が

同じnであっても直交し無関係であることを意味しますが,

これは例えば添字rが運動量意味するBose粒子では,異なる

での個数は無関係なので当然の規約です。

 

これらの交換関係は,改めて固有値と固有ベクトルの表現を添字

付きで,|n>=n|n>etc.とすれば,一般化された固有

ベクトルは,|..,n..,n..>=..|n>..|n>..の

ような直積表現と解釈され,

 

r≠sなら.|..,n..,n..>

=(n1/2(1/2|..,n-1 ..,n-1..> 

|..,n..,n..>,

 

|..,n..,n..>

=(n+1)1/2(n+1)1/2|..,n+1 ..,n+1..> 

=a|..,n..,n..>n..,n..>,

 

そして, 

|..,n..,n..>

=n1/2(n+1)1/2|..,n-1 ..,n+1..> 

=a|..,n..,n..>n..,n..>であり,

 

r=sなら,|..,n..,n..>

=(n+1)|..,n..,> 

(+1)|..,n..,n..>となりますから

明らかです。

 

Fermi統計に従うFermi粒子の場合,|n>=n|n

のnは0,1の2つの値しか取り得ないので,消滅,生成演算子

,|n>=n1/2|n-1>,

|n>=(n+1)1/2|n+1>で同じように定義しても,

 

固有状態 |n>としては,固有値nがn=0,1の状態

|0>,|1> の2つしか存在しないため,全ての作用を具体的に

書くことができて,|1>=|0>,|0>=|1>であり,

一方,|0>=|1>=0 でなければなりません。

 

それ故,|1>=0,|0>=0 です。 

また.|0>=0,|1>=0 も明らかですから

演算子として0です。

 

他方,|0>=|0>,|1>=|0>ですから, 

|1>=0, |0>=0の成立と合わせると, 

()|0>=|0>,

()|1>=|1> を得ます。

 

よって,Fermi粒子の場合は,任意の演算子,に対して

反交換子{,}{,}≡ABBAで定義して, 

{r,}=δrs,{r,}={,}=0  

と規約すれば,全ての辻褄があって無矛盾となります。

 

Fermi粒子でも個数演算子はであり,

|n>=n|nが成立します。

 

なお,粒子が存在しない状態 |0>に対しては,Fermi粒子だけ

でなくBose粒子でも|0>=0 が成立して,これ以上は粒子

を消滅させることができないとしています。

 

|0>に対してこのように規約すれば.Bose粒子の

|n>=n|nを満たす固有値nについても許される

のは非負の整数のみ:n=0,1,2..です。

 

そして,全てのrについての粒子がゼロの状態 |0>の直積で

与えられる全く粒子の存在しない状態 

|..,n..,n..>=|0,..,0...,0,..>を真空(vacuum)と

呼ぶことにします。

 

量子論では零点振動とか零点エネルギーと呼ばれる真空の

エネルギ-なるものの存在について論じられることもありますが

そうしたものは無視して,真空はHamiltomian の固有値がゼロ

の固有状態であるとします。

 

つまり,真空 |0,..,0...,0,..>はエネルギ^-がゼロの

基底状態(エネルギーの原点)であるとします。

 

すると,演算子は,εを添字rに対応する1粒子のエネルギー

として,=Σε=Σεと表現されます。

 

特に,粒子が自由粒子でr=(,σ)なら,ε=hc22/(2m)

です。

 

次に,とスピンsの関数としてこの粒子の場の演算子:

φσ(,s)を,φσ(,s)=V-1/2Σexp(ikx)δ(s,σ)

によって導入します。

 

ただし,δ(s,σ)はスピン部分の関数でδ(s,σ)はσ=s

のとき1でそれ以外はゼロの関数です。

 

スピンがゼロのBose粒子なら常にσ=s=0ですから,

δ(s,σ)=δ(s,0)≡1 であって,場はスカラー場

φ()≡φ0(,0)と表現されます。

 

一方,スピンが1/2のFermi粒子の場合,σ=±,で 

φ±(,s)=φ±()δ(±,s)と書くと,これは2成分が

φ±()のスピノール:(φ+(), φ())で表わされます。

 

さらに,スピンが1のBose粒子の場合,やはり,σ=1,-1,0の場

φσ(,s)φσ(,s)=φσ()δ(s,σ)と書いて,3成分

φ±(),φ0()の3行1列の列ベクトル

(φ+(),φ(),φ0())で表現されます。

 

そして,場φσ(,s)のHermite共役 φσ+(,s)を, 

φs+(,s)=V-1/2Σexp(-ikx)δ(s,σ)

とします。

 

改めて,交換子 [,]を[,],反交換子 {,}

を[,]と書くことにします。

つまり,[,]±AB-±BA(複号同順)です。

 

場の演算子に対しては,+の反交換子がFermi粒子,-の交換子

がBose粒子に対応して,交換関係,反交換関係: 

[φσ1(1) φσ2+(2)]±=δ3(12σ1,σ2, 

[φσ1(1) φσ12(2)]±=[φσ1+(1) ψσ2s2(2)]±=0  

が成立しています。

 

また,=Σεは,ポテンシャルUが存在する

場合,=∫d3xφσ+()[-hc22/(2m)+U()]φσ() 

と書けます。

 

何故なら,=Σε=Σεですが, 

簡単のためスピンがゼロのBose粒子を想定しrの添字rが運動量 

(波数)である場合を考えて,r()と書くことにすると, 

ε=ε()=2/(2m)+U() です。

 

ただし,U()はポテンシャルU()のfourier変換,または

位置共役な物理量である運動量表示のポテンシャル

です。

 

()=V-1∫d3U()exp(ikx)

⇔ U()=∫d3U()exp(-ikx) なる式で与えられます。

 

個数演算子もこの表示ではrrr()=()()

と書けます。

 

交換関係 [r,]=δrs,[r,]=[,]=0  

,r=1,s=2として,[(1)(2)]=δ3(12), 

[(1),(2)]=[(1),(2)]=0 となります。

 

故に,=Σε=∫d3ε()()() 

=∫d33ε()()(3 ()ですが, 

δ3 ()=V-1∫d3exp{-i()}なので, 

=V-1∫d3∫d3()exp(-ikx)

∫d3ε()()exp(ipx) 

=V-1/2∫d3xφ()∫d3[2/(2m)+]ε()

()exp(ipx) 

=∫d3xφ()[-hc22/(2m)+U()]φ()

となるからです。

 

(注釈):ここで第2量子化についての薀蓄を少し述べたく

なりました。

 

元々の通常の(第一量子化の)1粒子量子力学では

=-hc22/(2m)+U()であり,多粒子系では

=Σk{-hc22/(2mk)}+U(1,2..,)であった 

のに対して,

 第2量子化で全ての質量mkが同じmの同一粒子

ではH=∫d3φσ+()[-hc22/(2m)+U()]φσ()

と表現されるのには,どういう意味があるのでしょうか?

 

第2量子化(場の量子化)は波動関数を演算子と見なす操作で行う

ことが可能なため,量子論の基礎付けを根底から大転換する新理論

のように見えますが,実はそれほど重大な意味はなくて,単に表示

を位置座標表示表示から個数表示にユニタリ変換するだけの

表示の変換過ぎないと考えられます。

 

すなわち,抽象的なあるHilbert空間の状態ベクトルに作用

する線型作用素(演算子)としてHamiltonianを,運動量を

すると,2/(2m)+U です。

 

このとき,表示での演算子の行列要素は位置座標で対角化

されています。

 

つまり,表示では任意の演算子の行列要素<||>は

対角成分のみがゼロでない対角行列の形をしています。

 

 演算子が運動量なら,特にSchroedinger表現でのについて

表示では,<||>=-ihcδ3()であり,また,

ポテンシャル||>=U()]δ3()で与えられ

ます。

 

 Hについても,

||>=[-hc22/(2m)+U()]δ3() です。

 

多粒子系なら,位置座標が(1,2..,),および,

(1,2..,)の状態を挟むと,行列要素は,

1,2..,||1,2.., 

 =k|-hc2k2/(2mk)]+U(1,2..,)] 

 δ3(113(11)..δ3() です。

 

ただし,U(1,2..,)は,の行列要素が

1,2..,||1,2..,=U(1,2..,)

δ3(113(11)..δ3() なる対角成分

のみを持つとした多体系のポテンシャルを意味します。

 

状態 |Ψ>の波動関数,つまり位置表示での状態ベクトル

Ψ(1,2..,),|Ψ>に含まれる位置の固有状態

|1,2..,>の成分,

あるいは,その位置固有状態で展開した展開係数であって

粒子の存在確率の確率振幅を表わすものですから

Ψ(1,2..,)=<1,2..,|Ψ>なる式で与え

られます。

 

 ここで,便宜上,=(1,2..,)と略記すると,

 Ψ()=<|Ψ>です。

 

 そこで,

<Ψ||Ψ>=∫d<Ψ|><||><|Ψ> 

=∫d<Ψ()<||>Ψ() 

=∫d Ψ()()Ψ() となります。

 

ただし,ここでも<||=H()δ(),と書け

ますが,≡d3132..d3, 

δ()≡δ3(113(11)..δ3() 

と略記しました。

 

つまり,1粒子量子力学の波動関数ψ()で張られる2乗可積分

関数の状態空間に作用するHamiltonian演算子は

()=-hc22/(2m)+U()であり,多粒子系の波動関数

Ψ()で張られる2乗可積分関数の空間に作用するHamiltonian

演算子は()=Σk|-hc2k2/(2m)]+U()です。

 

一方,個数表示では,添字rを省略すると,

であり,|n>=(n!)-1/2()|0>です。

 

任意の演算子の個数表示による行列要素を表示での

完全系条件∫d><|=1を挿入して表示で展開

すると,

<m||n>=∫d<m|><||><|n>

であり,

一方,任意の状態におけるの期待値は

<Ψ||Ψ>=Σm,n<Ψ|m><m||n><n|Ψ> 

=Σn<Ψ|n>n<n|Ψ> と書けます。

 

ただし,<m||n>=Aδmn です。

 

特に,については,

<Ψ||Ψ>=Σm,n<Ψ|m><m||n><n|Ψ> 

=Σn<Ψ|n>n<n|Ψ>であり,

<m||n>=Hδmnです。

 

 この個数状態 |n>全体で張られる空間の状態に作用して

 <m||n>を与えるHamiltonian演算子については, 

 スカラー場のBose粒子の場合に,

 H=Σε

 =∫d3xφ()[-hc22/(2m)+U()]φ() 

 =∫d3xφ()()φ()で与えられることを既に

 見ました。

 

要するに,それぞれのが異なるのは,それが作用すべき状態

の作る状態空間が異なるためです。

第2量子化した場合も,単に状態空間が変わるだけで,本質的 

にそれらは,演算子の期待値,あるいは行列要素がどの空間でも

一致するようにユニタリ変換で結びついていて謂わゆる表示が

異なるだけで量子力学が新しくなったわけではありません。

(注釈終わり※)


 このところ,薀蓄=脱線注釈部分を書いていて参考書も

なく頭の中にある知見をヒケラかそうとする余り,文章の堂々巡りが多くて長時間を費やしてしまいいました。

※本ブログ2007年8/8の過去記事「量子力学の基礎(表示の話)(1)」, 

および,「量子力学の基礎(表示の話)(2)」も参照してください。

 

今日はここまでにします。


(参考文献):阿部龍蔵 著「統計力学(第2版)」(東京大学出版会)


 PS:これに関連して,26歳(1976年)の春4月に一人で電車を

乗り継ぎ,途中東京で東大の友人のところに一泊し翌日青函

連絡船乗って早朝に船酔いでゲロを吐きつつ初めての北海道

に渡り,旅館に一泊後,北大に遠路はるばる博士課程を受けに行

ったのを思い出しました。


(※今は博士課程(大学院後期過程)もあるようですが当時

は私の在籍大学には修士課程しかありませんでしたから。。)


 テストは当然,修士
論文説明であろうと予想していたの

に,思いもかけず,研究室個室朝から晩まで期限無しで参考

書を見たり外出するのも自由の4,5問程度のペーパーテスト

を出されました。

 
 その第1問目が,Diracのテキストにあるような量子論の定式化

の話で,上記と同じような話を如何にして明快な解答にするか?

に夢中になった余り,堂々巡りの悪戦苦闘で,もちろん,参考書

などは持参してないし,まさか無期限で外出自由でも慣れない

札幌の地で大学内や外の図書館や本屋まで行って調べようと

いう気にもなりません。

(※慣れない一人旅で,心身もかなり疲れていました。)


 結局,第1問を考えているだけで夜の8時になり,往復5千円

の切符しかもってなくて,旅館にもう一泊する余裕もなく,

あきらめて,もう帰宅してしまっている出題した教授の部屋

ドアのポストに答案用紙を入れて帰途についたのを昨日

ことのように思い出しました。


 未だに覚えているのは,よほど悔しかったからでしょう。。

これも人生を左右する1つの岐路でしたね。。

| | コメント (0) | トラックバック (0)

2014年11月27日 (木)

統計力学の基礎(8)(統計力学の基礎付け2)

 統計力学の基礎付けの続きです。

 前回の終わりでは,

 もう1つ気になるのが,統計集団(アンサンブル)を考える際の純粋集団(純粋状態)と混合集団(混合状態)の話です。

と書きました。

その続きです。

 

前々回の記事の量子統計力学2では,分配関数ZはN粒子系の多体問題の純粋状態のトレースによりZ=Tr[exp(-β)]で与えられることを示しました。

 

まあ,小正準集団を想定するなら集団を構成する個々の微視状態を

普通の干渉も有り得る観測前のコヒーレントな量子状態(純粋状態)

と考えることも可能てす。

 

しかし,正準集団の考え方では,莫大な個数 M個の同じ系を集めて全体

として孤立系になるような集団とするわけですが,この集団の個々の系

は全て同じ構造を持つとはいっても,それらは既に純粋状態ではなく,

エネルギ「ーEを持つ系(=純粋状態)が率exp(-βE)/Zで混合された

デコヒ-レントな混合状態と考えられます。

 

 統計集団というのは,統計力学の理論を定式化するために便法

として想定されたものに過ぎないので,それが純粋集団か混合集団

か?のようなことについて思い煩う必要はないのですが。。。

私はやはり気になります。

 

これについては,2006年10/23の記事「観測の理論(デコヒーレンス)」

があります。

これもまず全文を再掲載します。

 

(※以下,過去記事の全文です。)

 

今日は,観測に伴なって固有状態の干渉項が消滅すること

=デコヒーレンス (decoherence)の現象を最近の理論に基づいて

述べてみたいと思います。

 

ただし私自身は本質的には多世界解釈の方に傾いています。

まず,"観測可能量(observable)=物理量=線型演算子"その

あらゆる固有値:oiに属する固有状態:|i>の集合,つまり,

|i>=i|i>を満たす|i>の集合があり,これら

完全系を形成している,すなわち,∑i|i><i|=1

が成立しているとします。

"任意の状態=純粋状態":|ψ>は|ψ>=∑ii|i>と展開可能

この同じ状態|ψ>において,物理量を状態を乱すことなく

独立に多数回観測したときにはの固有値以外が観測されること

はなく,観測値がoiである確率が|ci|2で与えられます。

そして∑i|ci|21が成立しているというのが量子力学

の観測に関する枠組みと考えられます。

 

しかも,通常は固有状態|i>は正規直交化されていて,

<i|j>=δijなのでci=<i|ψ>なる式が成立して

います。 

したがって,この純粋状態|ψ>における物理量O^の観測値

の"期待値=平均値"は,

ψ=∑i|ci|2i=∑ii|ψ><ψ|i><i||i>

=∑i<ψ|i><i||i><i|ψ>=<ψ||ψ>

で与えられます。

 

つまり,<ψ=<ψ||ψ>であり.

ψ=∑i<ψ||i><i|ψ>=∑i<i|ψ><ψ||i>

=Tr(ψ)となります。

ここで射影演算子とよばれるψψ|ψ><ψ|で定義され,

物理量の対角和(trace)が,Tr()≡∑i<i||i>と定義

されています。

 

そして対角和の値が,これを定義する完全系{|i>}の選択に依

ないことも簡単にわかります。

ところで,もしもこの体系が,状態間の干渉が存在するような状態

の重ね合わせのみで成り立つ純粋状態ではなく,

 

情報の欠如などによって統計的に純粋状態:ψ,φ,χ,...が

それぞれ確率:W(ψ),W(φ),W(χ),...で混合している混合

状態であるとすれば,

 

の期待値は<>=∑ψ(ψ)<ψ||ψ> 

で与えられます。

 

これも,<>=∑ψ(ψ)<ψ||ψ>

=∑ψiW(ψ)<ψ||i><i|ψ>

=∑iψ(ψ)<i|ψ><ψ||i>=Tr(ρ)

となり,純粋状態の<ψ=Tr(ψ)と同じ形に

書けます。

 

ここでρ^はρ≡∑ψ(ψ)|ψ><ψ|=∑ψ(ψ)ψと定義

されて統計作用素,または,密度演算子(密度行列)と呼ばれます。

対象となる体系のHamiltonianとすると統計作用素ρ

時間に依存する量子力学の線形演算子に相違ないので,

Heisenbergの運動方程式:ic(∂ρ/∂t)=[,ρ]

を満足します。

 

ただし,hc≡h/(2π)でhはPlanck定数です。

 

実は状態|ψ>はSchroedinger表示の時間を含む状態ベクトル

|ψ(t)>で,これがSchroedingerの方程式:

ihc(∂/∂t)|ψ(t)>=|ψ(t)> を満たします。

 

逆に統計作用素ρ≡∑ψ(ψ)|ψ(t)><ψ(t)|が時間

を含むHeisenberg表示の作用素となるため,Heisenbergの

運動方程式:ihc(∂ρ/∂t)=[,ρ]を満たすと考えて

よいわけです。

時間発展の演算子を(t',t)=exp{-i(t'-t)}と

すると,|ψ(t')>=(t',t)|ψ(t)>ですから,

ρ(t)≡∑ψ(ψ)|ψ(t)><ψ(t)|によって

ρ(t')=(t',t)ρ(t)(t',t)-1となります。

 

統計作用素ρの時間発展はユニタリ変換によって行われる

のでρ,ρに関わる関係式は時間発展によって変化しません。

簡単のため,スピンが1/2の区別できる粒子が2個ある体系に

ついて考察します。

 

スピン1/2の1粒子のスピン角運動量の演算子をとすると,

それは2行2列の行列表示では,Pauliのスピン行列σを用いて,

=(c/2)σと表わされます。

 

σz の固有値+1,-1の固有状態を,それぞれ|α>,|β>とします。

 

2つの粒子それぞれのこうした状態を,それぞれ,(i)

と|β(i)>(i=1,2)で指定することにします。

このとき,全系の任意の状態ベクトルは(1)>|α(2)>,

(1)>|β(2)>,|β(1)>|α(2)>,|β(1)>|β(2)>の1次結合

で表わされます。

 

そして,例えばスピンがゼロの状態は,

|0>=(1/21/2)(|α(1)>|β(2)>-|β(1)>|α(2)>)

で与えられます。

 

この状態での統計作用素ρ^0は,

ρ0 =(1/2)(|α(1)>|β(2)>-|β(1)>|α(2)>)

(<α(1)|<β(2)|-<β(1)|<α(2)|)

=(1/2)(|α(1)><α(1)|(2)><β(2)|

-|α(1)<β(1)|(2)><α(2)|

-|β(1)><α(1)|(2)><α(2)|

+|β(1)><β(1)|(2)><α(2)|)

となります。

 

ただし,記号は直積を表わしています。

 

このρ0は確かに,純粋状態を示す"統計作用素=射影演算子"

です。

ここで,一般に粒子1のみに関する物理量(1)を測定する場合

を想定すると,このときも対象としては全体系ですから,

物理量を表わす作用素(1)(2)です。

 

その期待値は,

(1)(2)>=Tr(ρ(1)(2))

=∑ij<i(1)|<j(2)|ρ(1)|j(2)>|i(1)

=Tr(ρ(1)(1))  と書くことができます。

 

ここで,ρ(1)<j(2)|ρ|j(2)>=Tr,2(ρ) です。

そして,部分系である粒子1の物理量(1)の測定の期待値は全て

(1)(2)>=Tr(ρ(1)(1))の形で表わせるので,

 

実質的には,ρ(1)が部分系である粒子1の状態を示す統計作用素

であると見なすことができるでしょう。

ここで,ρρ0 の場合には,

ρ0(1)=<α(2)|ρ0(2)>+<β(2)|ρ0(2)

=(1/2)(|α(1)><α(1)|+|β(1)><β(1)|) です。

 

そこで,全系が純粋状態でも,部分系である粒子1の状態は

z成分のスピンが上向きと下向きが1対1に混合した混合状態

となることがわかります。

話を戻して,体系の状態が|ψ>で物理量O^の固有状態での

展開が,|ψ>=∑ii|i>(∑i|ci|21) で与えられると

します。

 

O^の測定装置はマクロな物体ですが,装置も状態ベクトル

表わすことができると仮想して,その初めの状態を|o>A

とします。

そして,対象が状態|i>にあるとき,それを測定したときの

"対象=体系と装置"の変化を|i>|o>A|i>|i>A

とします。

そこで,|ψ>を測定したときには,

|ψ>|o>A → ∑ii|i>|i>A となります。

 

この最後の状態はもちろん純粋状態であって,物理量の期待値

を取れば当然|,i>|i>A 間の干渉が現われるはずです。

 

最初の状態が純粋状態であって時間発展がユニタリですから当然

それは予想されたことです。

しかし,我々の観測の経験では,測定の最後の状態は

|i>|i>Aの状態がW(oi)=|ci|2の確率で混じり合って

いて,決して干渉作用など起きない混合状態です。

簡単のために,1電子のスピンのz成分を観測するStern-Gerlach

の実験のようなものを考察します。

 

これは,不均一な磁場の中にスピン磁気モーメントを持つ電子

が入射してスピンが上向きか下向きかが検出される実験です。

 

入射電子はある一定のスピン状態にあって,

|ψ>=(c1|α>+2|β>)|φ>,(|c1|2|c2|21)

であるとします。

 

ただし,|φ>は電子線の空間的運動を表わす状態ベクトルです。

 

入射電子が磁場の中を通るとスピンの向きによって空間的運動

は上下に分裂するので,

 

|ψ> → |γ>≡c1|α>|φ>+2|β>|φ

 

となります。

そして,上下にある検出装置の統計作用素=密度行列をそれぞれρAα,ρAβ,対象と装置の全体系の"統計作用素=密度行列"をη0

すると,

 

η0=|γ><γ|ρAαρAβ

=(|c1|2ρ++|c2|2ρ--12*ρ+-21*ρ-+)ρAαρAβ

 

と書けます。

 

ここで,

 

 ρ++=|α><α|><φ|,

 ρ--=|β><β|><φ|,

 ρ+-=|α><β|><φ|,

 ρ-+=|β><α|><φ|

 

です。

測定装置が状態ベクトルで表わされている状況では,ユニタリ性の故,

測定の結果として,干渉項ρ+-,ρ-+が消えることは決して有り得ないことです。

 

そこで装置は初めから混合状態にあると考えます。

 

すなわちマクロな装置はN個~ Avogadro数個程度の粒子の集合系

であり,このN粒子の系の多数の状態ベクトルの混合状態が装置を

表わしていると考えるわけです。

そして,測定にはある時間にわたって全体系の密度行列η0を調べる

必要があります。

 

それぞれ,N,N'粒子系から成る上下の検出装置に対して

η0(N,N')≡|γ><γ|ρAα(N)ρAβ(N') と定義

します。

 

相互作用が起こる直前の時刻をt0 として,時刻tでの全体系の

統計作用素をN,N'を省略してη(t)と書くと,η(t0)=η0

対し, 

η(t)=∑N,NW(N)W(N')(t,t0)η0(N,N')(t,t0)-1

(ただし∑(N)=1)

と書くことができます。

η0(N,N')=|γ><γ|ρAα(N)ρAβ(N')において,

例えばρ+-に関わる部分は,

|α><β|><φ|ρAα(N)ρAβ(N') です。

 

装置との相互作用部分がスピンに依らないとすれば,時間発展

は,(t,t0)ρAα(N)<φ|ρAβ(N')(t,t0)-1

となります。

 

ここで,|φ>はρAα(N)のみ,<φ|はρAβ(N')のみと相互作用

するので左右に分けました。

 

tを相互作用が終わった時刻とし,N個の粒子の個数に比例する

運動長さの単位をL(N)とすると,そのオーダーはL(N)~N1/3

です。

 

そして,比例定数として波数因子kを掛けた位相の変化がある

と考えられるので,

左の(t,t0)>ρAα(N)は因子exp{ikL(N)}を,

右の<φAβ(N')(t,t0)-1は因子exp{-ikL(N')}

を含むはずです。

ここで,η(t)=∑N,N’...を連続化して積分式にすると,

η(t)=∫dL∫dL'W(L-L0)W(L'-L0)(t,t0)

η0(L,L')(t,t0)-1(ただし∫dL(L-L0)=1)

となります。

 

位相部分だけに着目すると,L(N)~N1/3が大きい極限で

密度行列要素は,それぞれ,

 

ρ++ → 1,ρ-- → 1,および, ρ+-→ exp{ik(L-L')},

ρ-+ → exp{-ik(L-L')}

 

となります。

 

ところで,Riemann^Lebesgueの定理によれば,L,L'が無限大の

極限では,

∫dL∫dL'W(L-L0)W(L'-L0)exp{ik(L-L')} → 0

となります。

 

このことから"統計作用素=密度行列"からρ+-ρ-+の干渉項

が消えてρ++ρ--の項のみがそのままの形で残ることになり,

事実上デコヒーレンスが実現されることになると考えられます。

ただし,清水明氏の量子測定の原理とその問題点」に書かれて

いますが,

"測定装置の他に環境も含めたとしても干渉項のオーダー

観測時間をT,光速をcとして,exp[-(正定数)×cT3]が限界

あり決して正確にゼロになって消えるわけではない。"

という問題は残っています。

一方,szの測定によって必ずしもσzの固有状態である

|α>,|β>が観測されると考える必要はないという本質的な

問題もあります。

 

例えば,|χ±>≡(1/21/2)(|α>±|β>)(複号同順)はσx

対してのスピンの+,-の固有状態です。

 

先の統計作用素において非干渉成分として,

ρ++=|α><α|><φ|,ρ--

=|β><β|><φ| の代わりに,

 

ρ'++=|χ><χ||φ'><φ'|,

ρ'--=|χ><χ||φ'><φ'|

 

が残ると考えても何の不都合もないからですね。

こちらの問題は(猫生)か(猫死)のどちらか一方のみの状態が観測

されるとして定式化しても,

 

それらの重ね合わせ状態が観測されるとして定式化しても,

"統計作用素=密度行列"のデコヒーレンスだけからは,

それらは全く同等である,ことから多世界解釈の問題でもあり

超選択則関わる問題ですね。

 

例えば変換群の異なる既約表現にまたがる重ね合わせ状態は

観測されない,とかの原理的問題であると思います。

 

具体的には既約表現の問題とは,ちょっと違うかもしれないです

が,アイソピン(荷電スピン)に関わる2次元特殊ユニタリ群

SU(2)において,

 

陽子と中性子の重ね合わせ状態は決して観測されない,という

のも超選択則の例です。

 

これに対して,φメソンやKメソンにはむしろ混合(mixing)が

ある状態で存在する方が普通なので,自然がどういうメカニズム

になっているのかは不思議なことです。

 

これに関しては,観測を行なう以前の物理系の状態を記述する

"波動関数や密度行列をも実在であると考えるかどうか?"

という哲学的な問題も関連あるかもしれません。

 

参考文献;町田茂 著「基礎量子力学」(丸善),

ボーム 著「量子論」(みすず書房)

(再掲載終了※) 

さて,再び,考察しますが,正準集団,大正準集団が混合集団であることは疑う余地がりません。

 

しかし,小正準集団についてははっきりしません。

 

古典統計力学では,小正準集団に属する微視的状態は相空間の

個々の軌道であり,それらの状態はもちろん干渉するわけではなく

運動方程式のN個の厳密な解の全ての情報から巨視的平均量確率

分布を求める代わりに,統計的原理を導入,それに基づいてそれら

の量を求めるのですから,集団には完全な情報が欠如しているのは

明らかです。

 

しかし,量子論では,古典論と対比すると,N粒子の系の完全な情報

が,そもそも確率波という確率的なものです。

 

そして,個々の微視状態を多体系の波動方程式の解の重ね合わせで

与えられる純粋状態と考えることもできそうです。

 

純粋状態だけで熱平衡状態を表現できれば情報の欠如なく定式化

できると思いますが,これはシミュレーション計算をする計算物理学

の分野でしょうか?。

 

統計力学の原理に基づいて最大確率を与える分布を求めるという

旧来からの統計力学の思想に合致する対象は,小正準集団でも,干渉

などしない混合状態の集団です。

 

実際,系の巨視的量の観測はマクロな測定装置を用いて行なうわけ

ですから,その時点で既にマクリな装置との干渉でデコヒーレント

な状態に移行した後の系を対象としていると考えるのが普通です。

 

ネットを検索してみると,私の種ひゅうの疑念と意図を共有する

 

ものとして,大阪市大,杉田歩氏のPDF

;量子統計力学の基礎付けについて;」, および,

すう理化学化学科学2013年6月号の特集記事として,

量子純粋状態による統計力学の定式化がありました。 

 

今日はこれで終ります。

 

| | コメント (0) | トラックバック (0)

2014年11月24日 (月)

統計力学の基礎(6)(量子統計力学3)

  量子統計力学の続きです。

 

§3.状態和に対する表式 

莫大な数 N個の粒子から構成される系全体のHamiltonian

とし,そのk番目のエネルギー固有値をE,固有関数

Ψ(1,2,..,N)とします。

 

Ψ=EΨです。

 

ここでxの関数f(x)をf(x)=Σn0=で定義します。 

このとき,xをで置き換えたものとして演算子 f()を定義

できます。

 

が線型演算子なのでf()も線型演算子です。

 

Ψ=EΨより,2Ψ=E2Ψ,..,Ψ=EΨ

となるため,(=f(Eとなります。

 

したがって,exp(-β=exp(-βEです。

 

ここで,Diracのブラケット記法に移行します。 

>=E>,exp(-β)|Ψ>=exp(-βE

etc.です。

 

|Φ>=Φ(1,2,..,N),|Ψ>=Ψ(1,2,..,N)に対して,

スカラー積(ユニタリ内積)を,

<Φ|Ψ>=∫Φ(1,2,..,N)Ψ(1,2,..,N)dτ1τ2..dτ 

で定義します。

 

の固有関数>=Ψ(1,2,..,N) (k=1,2,..)は規格化

直交条件<Ψ>=δijを満たし,|Ψ>,|Ψ2>,.. は,

完全系を構成するとします。

 

つまり,任意の|Ψ>は,|Ψ>=Σ>;c=<Ψ|Ψ>

と展開されます。これは,|Ψ>=Σ><Ψ|Ψ>とも書ける

ので,完全系をなすことは><Ψ|=1と表現されます。

 

 一方,任意の演算子をとするとき,<Φ|Ψ>の<Φ|を<Ψ|

とし,|Ψ>を>とすると,<Ψ|>ですが,これを 

i行j列の行列要素とする行列を同じ記号で表わすこtとに

します。

すると,exp(-β)|Ψ>=exp(-βE)|Ψ>によって, 

正準分布の状態和=分配関数Zは,

Z=Σexp(-βE)=Σ<Ψ|exp(-β)|Ψ

と表わされます。

 

行列A=(Aij)のトレース(trace:対角和)を,TrA=Σkk

で定義すると,ij=<Ψ|>の場合は,

Tr=Σkk=Σ<Ψ|> です。

 

以上から,状態和(分配関数)は量子論的期待値の和として,

Z=Tr[exp(-β)と書けることがわかります。

 

これを大きな状態和=大分配関数に拡張するには少し注意を

要します。何故なら,それは粒子数の変動を伴うからです。

 

エネルギーだけなら量子論の物理量としてHamiltonian

いう対応する演算子がありました。

 

そこで,自由粒子の場合,以前Bose分布とFermi分布の導出の際

に与えた,エネルギー固有値がεの1粒子状態を占有する粒子

数n,個数演算子と見なしてと書き,

=Σε,=Σとすれば,大分配関数Zは上記の

分配関数Zのケースと同様にして,Z=Tr[exp{-β(-μ)}]

と書けます。

 

この個数nを演算子と見なする法を第2量子化といいます。 

謂わゆる場の量子化ですね。

 

§4.古典的極限(高温極限)

 

これまでと同じく,莫大な数 N個の粒子が体積Vの箱の中にある

として考察していますが,粒子間には相互作用が働くとします。

そのポテンシャルをUとし,これは粒子の空間座標1,2,..,

のみの関数と仮定します。

 

特に2体力のみならvijを粒子iと粒子jの間に働く力の

ポテンシャルとして,U=Σi<jijと書けます。

 

しかし,以下の論議はより一般的で,UをU=Σi<jijの形に

限定する必要はありません。

 

,対象としている体系の全Hamiltonianは,

0+U(1,2,..,)であり,

0={-hc2/(2m)}Σk=1であるとします。

 

 H0とUは一般には演算子として非可換ですが,古典的極限では

 交換するとしてよいので,

 exp(-β)=exp{-β(0+U)}~ exp(-βU)exp(-β0)

 です。

 

 そこで,状態和は,この古典極限では,

 Z=Trexp(-β)=Tr[exp(-βU)exp(-β0)]

 です。

 

 系がFermi統計に従うとし,上記の状態和のトレース表現

 Z=Trexp(-β)=Σ<Ψ|exp(-β)|Ψ>のトレース

 を取る状態>して,Slater行列式:

 Ψ(1,2,3,..N)=(N!)-1/2Σ(-1)δ(P)

 ψr1(1)ψr2(2)..ψrN(N)を使います。

これは,前に粒子間に相互作用がある場合には自由粒子の集まり

とは異なって,全系の波動関数は1粒子のエネルギー固有関数の

の一次結合の形にはならない,と記述したことに矛盾するよう

ですが。。

2007年6/15の過去記事

ハートリー・フォック(Hartree-Fock)近似(1)」での電子の

集まりに対するHartree近似,すなわち,独立電子近似のような近似

を考えれば相互作用が弱い場合には矛盾ではないと考えられます。

 

したがって,(r1,r2,..,r)の組に対応する1つの,Slater行列式 

 Ψ(1,2,3,..N)=(N!)-1/2Σ(-1)δ(P)

 ψr1(1)ψr2(2)..ψrN(N) r1,r2,..,rN(1,2,3,..N)と書けば,

 Z=Trexp(-β)=Σ<Ψ|exp(-β)|Ψ>のトレース

 を取るべき)|Ψ>の添字kは,今の場合,k=(r1,r2,..,r)

 で与えられます。

 

 ただし,各rj(j=12,..,N)は1粒子波動関数ψrj(j)の添字

 であり,このrjは1粒子エネルギー固有値εjに対応しています。

 

 それ故,トレースの総和Σはあらゆるエネルギー準位の組

 (r1,r2,..,r)Slator行列式 ψr1,r2,..,rN(1,2,3,..N)

 にわたって取られます。

 以上から,Z=Trexp(-β)=Σ<Ψ|exp(-β)|Ψ 

 =Σ(r1,r2,..,rN)exp(-βU)ψr1,r2,..,rN(1,2,3,..N)

 exp(-β0)ψ1,r2,..,rN(1,2,3,..N)dτ1τ2..dτなる

 表式を得ます。

 

この表式のΣ(r1,r2,..,rN)で総和される内容を,F(r1,r2,..,r)

と書いて,Z=Σ(r1,r2,..,rN)F(r1,r2,..,r)と表現すると,

F(r1,r2,..,r)は,次の性質を有します。

 

 (ⅰ)F(r1,r2,..,r)は(r1,r2,..,r)の個々の引数成分の交換

 に対して対称である。 

 (ⅱ)F(r1,r2,..,r)は(r1,r2,..,r)のうち2つ以上のrk

 一致するときにはゼロである。

 

 以上の2点からr1<r2..<r)の制限付きの1つの

(r1,r2,..,r)の組に対して(N!)個の同じF(r1,r2,..,r)

が対応するため,状態和は,Z=Σ(r1,r2,..,rN)F(r1,r2,..,r)

=,(1/N!)Σr1,r2,..,rNF(r1,r2,..,r)と表わされます。

 

ただし,左辺のΣ(r1,r2,..,rN)の添字の組(r1,r2,..,r)は, 

 r1<r2..<r)の制限付きで,右辺のΣ(r1,r2,..,rNの添字

 r1,r2,..,rには,そうした順序の制限は無しです。

 

 そこで,今の場合の,

 Z=Trexp(-β)=Σ<Ψ|exp(-β)|Ψ>を再び,

 陽に書くと,(r1,r2,..,r)をr1<r2<..<r)に制限して

 Z=Σ(r1,r2,..,rN)exp(-βU)ψr1,r2,..,rN(1,2,3,..N)

 exp(-β0)ψ1,r2,..,rN(1,2,3,..N)dτ1τ2..dτ

 ですから,

 

 Z=(1/N!)Σr1,r2,..,rN∫exp(-βU)(N!)-1Σ(-1)δ(P) 

 [ψr1(1)ψr2(2)..ψrN(N)]exp(-β0)

 [ψr1(1)ψr2(2)..ψrN(N)]τ1τ2..dτN です。

 

ここで,ρ(,)≡Σψ()exp{βhc2/(2m)}∇() 

 とおくと,固有関数の完全性 |r><r|=1

 or Σψ(()=δ3()

 =(2π)-3∫d3exp{ik()} が成立しますから 

ρ(,)=(2π)-3∫d3exp{i()-βhc2/(2m)}

ですが,

i()-βhc2/(2m)

=-{βhc2/(2m)}{-im/(βhc2)}2-m//(2βhc2)

なので,

ρ(,)=(2π)-3exp[{-m/(2βhc2)}()2] 

∫d3exp{{βhc2/(2m)}{-im/(βhc2)()}2

 

 結局,ρ(,)

 ={m/{2πβhc2}}3/2 exp{-m//(2βhc2)()2}

 を得ます。

 

 古典的極限, or 高温極限では,βhc2 ~ 0 ですから,上の最後

 に得たGauss誤差関数の表現によると,xのρ(,)の

 寄与は,x=yのρ(,)=ρ(,) の寄与に比してほぼ

 ゼロであり無視できます。

 

 そしてρ(,)

 ={m/{2πβhc2}}3/2 exp{-m//(2βhc2)()2} 

 でx=yとすると,ρ(,)={m/{2πβhc2}}3/2 を得ます。

 

 したがって,βhc2 ~ 0 の極限では, 

 Z=(1/N!)Σ(r1,r2,..,rN∫exp(-βU)(N!)-1

 Σ(-1)δ(P) [ψr1(1)ψr2(2)..ψrN(N)]

  exp(-β0)[ψr1(1)ψr2(2)..ψrN(N)]τ1τ2..dτにおいて,

 Σ(r1,r2,..,rNの各項で,1つの添字rに対するψrjとψrj

 の積因子による関数:ρ(,)

 ≡Σrjψrj()exp{βhc2/(2m)}∇r1()に寄与する

 のは,総和Σ(-1)δ(P)のうちで=1(恒等置換)に対応する

 x=yの寄与のみです。

 

 これは,{m/{2πβhc2}}3/2をj=1,2,..NのN個掛け合わせた

{m/{2πβhc2}}3N/2ですが,ZはこれのΣ(r1,r2,..,rNのあらゆる

順列の(N!)個の総和を2つの(N!)因子で割ったもので

与えられますが,それ故,1つの(N!)は相殺されます。

 

 よって, 

Z=(1/N!){m/{2πβhc2}}3N/2∫d3132..d3exp(-βU) 

が得られます。

 

ここで.Q≡∫d3132..d3exp(-βU)とおき,

β=1/(kT),c=h/(2πを考慮すると,

Z=(1/N!){(2πmkT)3N/2/h3N}Q と書けます。

 

特に,理想気体でU=0ならQ=Vなので,a=h3N

置けば,Z=(1/N!){ V(2πmkT)3N/2/a}となり,前に

古典統計に基づいて求めた理想気体の分配関数をN!で割った

ものに一致します。

 

まだまだ,この項の草稿は続いてるのですが,長いので分ける

ことにして一旦,終わります。

(参考文献):阿部龍蔵 著「統計力学(第2版)」(東京大学出版会)

PS;11/22の土曜日には1回目の忘年会?がありました。

 私は極度の金欠で38円しかなかったので19時池袋待ち合わせには

 障害者無料のフリーパスで巣鴨から都パスで出かけ,池袋,大山,巣鴨

 と3件はしごしてタクシー代も含め全ておごってもらいました。

 もう7年半も前に心臓病手術直後にクビになった会社の先輩1名

 と後輩3名です。

 当時,門前仲町で降りて永大橋のそばに2000年から7年間通って

 て20時から朝7時まで勤務していた職場は既に無く,今は八王子

 で同じ仕事をもっと少人数で引き継いでるらしいです。

 後輩のうち1名には今回はじめて会いました。

 私が辞めた直後くらいに入社したらしいですが,いかにも人が

 好さそうな顔をしていました。

 イヤ,7年も前に辞めたのに,いまだに1年に数回はタダ酒

 さそってくれるのも私の人徳?のセイかな。。

 アリガタイことです。。

 夜中3時に帰宅して翌朝9時まで爆睡して日曜日はサワヤカに

 めざめました。

| | コメント (0) | トラックバック (0)

2014年11月19日 (水)

統計力学の基礎(4)(量子統計力学1)

 古典統計力学3の続きで量子統計力学に入ります。

 

§1.多粒子系の量子力学 

ここまでの古典統計力学における考察と同様,体積Vの容器の中

に莫大な数であるN個の自由粒子が閉じ込められているケースを

考えますが,統計力学を進める前にこうした多粒子系の量子力学

をレビューしておきます。

 

簡単のため,系全体のHamiltonianが系を構成する各粒子それの

単純和で与えられる場合を論じます。

 

N個のうち,k番目の粒子のHamiltonianを(k)とすると,

は,(1)(2)..+(N)と表わされます。

 

この場合,系を支配する定常状態のSchroedingerの波動方程式: 

Ψ=EΨの形式的な解 Ψは, 

Ψ(1,2…,N)=ψr1(1)ψr2(2)..ψrN(N),および, 

固有値 E=ε1+eε2..+εrN の組で与えられます。

 

ただしrk,および,εrkは,それぞれ,1粒子のSchroedinger

波動方程式:(rk) ψrk=εrkψrkを満たす1粒子の波動関数

(固有関数),および,エネルギー固有値です。

 

全系の解:Ψ(1,2,..,N)=ψr1(1)ψr2(2)..ψrN(N)は,

1番目の粒子がr,2番目の粒子がr..という状態,にある

ことを意味します。

 

しかしながら,量子力学の立場では,そのような粒子の個別性は

否定され,より制限された波動関数のみが許されます。

(粒子の分別不可能性)

 

一般に量子力学的粒子は位置座標の他にスピンと呼ばれる

内部自由度を持っており,その大きさは0.1/2,1,3/2,のような

値のみを取り得ます。

 

スピンが0.1.2,..の整数値をとるとき,その粒子はBose粒子

(Boson),1/2,3/2…のような半奇数値をとるときは,粒子は

Fermi粒子(Fermion)と呼ばれます。

 

自由粒子であれば一般にスピンがsの状態は(2s+1)重に縮退

しています。

 

そして,N粒子系が全てBose粒子の集まりであれば,それは

Bose統計(=粒子の交換に対して対称)に,一方,全てFermi

粒子の集まりならFermi統計(=粒子の交換に対して反対称)

に従うことがわかっています。

 

前に,N粒子の定常状態波動関数を

Ψ(1,2…,N)=ψr1(1)ψr2(2)..ψrN(N)と書きました

.例えばψr1(1)は,粒子1の座標変数を1,スピンをs1

とした関数ψ1(1,s1)の引数 (1,s1)をまとめて数1

で代表させたものです。

 

Bose統計では,その波動関数が粒子の交換に対して対称 

すなわち,Ψ(2,1,3,..)=Ψ(1,2,3,..)であり,Fermi統計

では反対称:すなわち.Ψ(2,1,3,..)=-Ψ(1,2,3,..) 

です。

 

そこで,一般にを(1,2,3,..)を(r,r,..,r)に置き換える

置換演算子とすれば,

Bose統計の場合,

Ψ(1,2,3,..N)=Ψ(1,2,3,..N),

Fermi統計の場合,

Ψ(1,2,3,..N)=(-1)δ(P)Ψ(1,2,3,..N)

です。

 

ただし,はδ(P)は,が偶置換のときは偶数,奇遺憾のときは

奇数を表わします。そして,この対称性は自由粒子だけでなく

粒子間相互作用があるときも成立します。

 

しかし特に自由粒子なら,こうした対称性を満たす波動関数を,

1粒子固有関数のN個の積:ψ1(1)ψr2(2)..ψrN(N)の

一次結合として,

Bose粒子ならΨ(1,2,3,..N)

 =(N!)-1/2Σψr1(1)ψr2(2)..ψrN(N),

Fermi粒子ならΨ(1,2,3,..N)

 =(N!)-1/2Σ(-1)δ(P)ψr1(1)ψr2(2)..ψrN(N) 

 と表わすことができます。

 

後者はi行j列の成分がΨ(i)の行列をΨとすると行列式 detΨ

意味します。これは,Slater行列式として知られています。

 また,全体にかかる係数(N!)-1/2は全確率が1となる条件

規格化条件:

∫Ψ(1,2,3,..N)Ψ(1,2,3,..N)dτ1dτ2..dτ=1

を満足させるために必要な規格化定数です。

 

ただし,記号∫dτkはスピンによる総和と積分:Σsk∫d3k

総称しています。

 

そして,Fermi粒子の場合, k≠lなるk,lに対し,もしもr=rl

なら行列式の性質によってΨ=0 です。

 

つまり,Fermi粒子ではk≠lなる異なる粒子k,lが,因子

ψrk(k)ψrl(l)によってr=rlなる同じ1粒子状態を占める

確率はゼロです。

 

したがって,同じ1粒子状態状態に2つ以上のFermi粒子が入ること

はできないことがわかります。

これをPauliの排他律(exclusion principle)といいます。

 

一方,Bose粒子にはこのような制限はなく同じ1粒子状態に原理的

には無限個Bose粒子が入ることができます。

 

この性質を表わすには,1粒子の状態rを占める粒子数nを考察

するのが便利です。

 

そうした粒子数nを考えれば,Bose粒子の場合には,

=0,1,2,...が,Fermi粒子の場合にはn=0,1のみが

許されます。

 

系の粒子の総数がNの場合N=Σです。そして,いずれの場合

も系全体のエネルギーEは,E=Σεで与えられます。

 

特に自由粒子では,[/(2m)]ψ=εψです。Plank定数を陽にh

と書く単位では,=-ihc∇ですからε=hc/(2m)と

書けます。ただし,hc=h/(2π)です。

 

これらの固有関数はスピンには依らず,2重に縮退しています。

 

ところで,粒子間に相互作用が存在すると,系の波動関数

Ψ(1,2,3,..N)は,

Ψ(1,2,3,..N)

=(N!)-1/2Σψr1(1)ψr2(2)..ψrN(N),または,

Ψ(1,2,3,..N)

=(N!)-1/2Σ(-1)δ(P)ψr1(1)ψr2(2)..ψrN(N)

ような形には書けません。

 

しかしながら,このような多粒子系AとBがあって,それぞれの系

内部では粒子が相互作用していても,AとBの間に相互作用がない

なら全体系のエネルギーEはA,BそれぞれのエネルギーE,E

の和としてE=E+Eをと表わすことができます。

 

.量子統計力学における正準集団,大正準集団 

,対象としている非常に多くの粒子から成る多粒子系と全く

同じ構造の体系がM個あるケ-スを想定してそれらを密着させて

並べます。

 

それらM個の系のそれぞれにはごく弱い相互作用がありエネルギー

を交換するものとします。

しかし,M個全体は孤立系で外部とは何の交渉もないとします。

 

古典論での出発点は.莫大な個数M個あると想定した同じ系の

うちでエネルギーがEにあるものの個数をMとするとき.

次の3つの点が成立するということでした。

 

 すなわち,(ⅰ)W=M!/{Π(M!)}, 

(ⅱ)(M1,M2,..)なる組で示されるMの配分が実現される確率は

Wに比例する。(ⅲ)平衡状態ではWが最大になっている。 

 の3つです。

 

 これらは「エルゴード仮設」,なたは「等重率の原理」に基づく

 ものです。

 

量子論でもこの状況は同じです。量子論の古典論との決定的な

違いは系を構成する微視的粒子の区別が不可能なこと:つまり,

1番目の粒子,2番目の粒子..というように粒子に順番を付ける

ことは量子論の世界では不可能であることです。

 

しかし,微視的粒子の集まりでなく巨視的な系のM個の集まりで

ある正準集団や大正準集団では,個々の系を区別して番号を付ける

ことが可能なので,これらの方法においては古典論との違いはない

と考えられます。

 

古典統計との違いが克明に現われるのは小正準集団の方法ですが

これについては別記事で後述する予定です。

 

さて,(ⅰ),(ⅱ),(ⅲ)の仮定が量子論でも正しいと認めて正準集団

の方法を進めると,古典論と同じく,

粒子数がN個と一定な系でエネルギーがEにある確率分布は,

exp(-βE)/Z;Z=Σexp(-βE)で与えられます。

 

さらに,粒子の交換も許す大正準集団では,粒子数Nの系で

エネルギーEk(N)の状態をN粒子の系が取る確率は,次式で

与えられます。

λexp{-βEk(N)}/Z=exp[-β{Ek(N)-μN}]/Z

です。

 

ただし,分母の大分配関数Zは,

=ΣN=0Σexp[-β{Ek(N)-μN}] 

です。

 

E=Σεを用いると,粒子個数がNでΣ=Nという制限

がある場合の自由粒子に対する状態和=分配関数:

Z=Σexp(-βE)は,Z=ΣΣnr=N exp{-β(Σε)}

と表わすことができます。

 

 しかし,大きな状態和=大分配関数では,Σ=Nという粒子数

 制限がないので,は各準位rについて独立に変わるとしてよい

 と考えられるため,

 

 =Σnr exp{-β(Σε-μΣ)}

 =Σnr exp[-β{Σ-μ) }]

 =ΠΣnr exp{-βn-μ)} となります。

 

 そして,粒子がBose粒子なら,n=0,1,2,..を取り獲るので

Σnr exp{-βn-μ)}=Σnr=0 exp{-βn-μ)}  

=1+exp{-β(ε-μ)}+exp{-2β(ε-μ)}+..(無限等比級数) 

[1- exp{-β(ε-μ)}]-1ですから,

=Π[1- exp{-β(ε-μ)}]-1 を得ます。

 

一方,粒子がFermi粒子なら,n=0,1のみなので 

Σnr exp{-βn-μ)}=Σnr=01 exp{-βn-μ)} 

=1+exp{-β(ε-μ)}ですから,

=Π[1+ exp{-β(ε-μ)}] です。

 

そこで,nの平均値をfと書けば, 

≡<n

=Σ,Nexp{-βE-μN}}/[ΣE,N exp{-βE-μN}] 

=Σnrexp{-βn-μ)}/Σnrexp{-βn-μ)} 

=-{∂Z/∂(βε)}/ Z=-∂{lnZ/∂(βε)} より,

 

Bose粒子,Fermi粒子のそれぞれについて, 

exp{-β(ε-μ)}/[1- exp{-β(ε-μ)}]

=1/[exp{β(ε-μ)}-1],

 

exp{-β(ε-μ)}/[1+ exp{-β(ε-μ)}]

=1/[exp{β(ε-μ)}+1] 

が得られます。

この分布則を,それぞれ,Bose分布,Fermi分布と呼びます。

これらは,β(ε-μ)=(ε-μ)/(kT)が大きいときには

=1/[exp{β(ε-μ)}±1]

~ exp{-β(ε-μ)}=exp{-(ε-μ))/(kT)}となって,

古典的なMaxwell-Boltzmann分布になります。

ところで,Z=exp(-βΩ) or Ω=-lnZ/βで熱力学ポテンシャル 

Ωを定義すれば,前と同じくΩ=-PVであり, 

dΩ=-SdT-pdV-<N>dμですから,エントロピー 

はS=-(∂Ω/∂T)V,μで与えられます。 

 

 Z=Π[1±exp{-β(ε-μ)}]±1 より  

Ω=-lnZ/β=-±Σ[ln[1±exp{-β(ε-μ)}]です。

 結局,種々の計算の末,(詳細は略。。)

S=kΣ[-flnf±(1±f)ln(1±f)]なるエントロピーの表現

を得ます。

 

短かいですが,本記事はここで終わります。 

 

すぐ後に小正準集団の方法を含む2007年の過去記事を再掲載する

予定です。

 

参考文献: 阿部流蔵 著「統計力学(第2版)」(東京大学出版会)

| | コメント (0) | トラックバック (0)

2014年11月16日 (日)

統計力学の基礎(3)(古典統計力学3)

 古典統計力学の続きです。

 

エネルギーだけの交換を許し粒子は閉じ込められた系の集団で

あった正準集団よりも一般的で,エネルギーの交換だけでなく,

さらに粒子の交換も許す場合を考察します。

 

こうした系の集団を大正準集団(grand-canonical ensemble)

といいます。

 

簡単のため,対象とする系は1種類の粒子のみから成る体系

とします。正準集団と同じく莫大なM個の同じ構造を持つ

を並べたと想定します。

そのうちのいくつかの系はN個の粒子から成る系で,その

個(巨大な個数)の粒子から成る系全体つくる相空間の

うちエネルギーがEijであるようなj番目の細胞にあるもの

数をMijとします。

 

こうした個数配分の配置数Wは,W=M!/(Πi.jij!)で

与えられます。

 

ここで,M個全体の粒子数をN0,総エネルギーをE0

としてこれらは一定とします。つまり,M個全体では孤立系

とします。

 

前と同じように,熱平衡で実現されるのは,条件:Σi.jij=M

Σi.jijij=E0,Σi.jiij=N0の下で,Wが最大になる

ようなMのMijへの配分であるケースです。

 

そこで,変分δMijに対してΣi.jδMij=0,Σi.jijδMij=0,

Σi.jiδMij=0 の条件下で,

δlnW=Σi.j(lnMij+1)jδMij =0 が成立する条件を

求めます。

 

これは,Lagrangeの未定乗数をα,β,γとして,独立な任意

のδMijに対し,恒等的に

Σi.j(lnMij+α+βEij+γNi)δMij=0 が成立すること

を意味します。

 

よって, lnMij+α+βEij+γNi=0,すなわち, 

ijexp(-α-βEij-γNi)です。

 

それ故,λ≡exp(-γ),Z≡Mexp(α)と定義して 

ij(M/ZNiexp(-βEij)と書くことができます。

 

そしてi.jij=Mから,Z=Σi.jλNiexp(-βEij)

です。

 

 以上から,大正準集団のM個の系うち,1個だけを考察対象の系

とし,残り全部を熱浴,および,粒子の貯槽と考えて,

対象とする系がエネルギーE,粒子数Nの状態を占める確率は,

λexp(-βE)/[Σi.jλNiexp(-βEij)]なる分布で与えられる 

ことがわかりました。

 

ここで,βやλの意味を知るため,特別な場合として粒子の交換

が許されず,i=N=一定として,

λexp(-βE)/[Σi.jλNiexp(-βEij)]の分子,分母を 

共通な定数因子λで割って約分すると,

exp(-βE)/[Σi.jexp(-βEij)]となりますが,これは前回論じた

正準分布に一致するはずです。

 

したがって,今回の場合でもβ=1/(kT)を得ます。

 

 一方,λについて論じるため,粒子交換が許される場合の熱力学

での2相平衡を考察します。

 

すなわち,2つの相が平衡にある条件は化学ポテンシャルμ

(1分子当たりのGibbs自由エネルギー)が互いに等しいこと

ですが,考察しているM個の系の集団では全ての系でλが

共通なので,これが相平衡でのμと同じ役割を果たしていると

考えられます。

 

 まず,前回は,∂(lnZ)/∂β

 =-Σexp(-βε)/Z=-<E>,および

 (∂F/∂V)=-P,F=-kTlnZより,

 ∂(lnZ)/∂V=-βP が成立したので,

 この,粒子数Nが不変でZの係数にかかるべき,λ

 が定数場合に相当するケースでは,

 d(inZ)=-<E>dβ+βPdVという結果

 でした。

 

 今回のinZでは,Nも変わるので.

 d(inZ)=-<E>dβ+βPdV+<N>d(lnλ)

 です。

 ただし,<N>=Σi,jiλNiiexp(-βEij)/[Σi.jexp(-βEij)]

 です。

一方,Gibbs自由エネルギーG=Nμ=F+PV=E-TS+PV

の微分を取ると,dG=dE-TdS-SdT+PdV+VdP

=-SdT+VdP+μdN となります。

 

何故なら,粒子の流出入を許す今の場合,熱力学第一法則は

閉鎖系でのdE=TdS+PdVではなくて,

dE=TdS+PdV+(∂E/∂N)S,VdNとなり,

一方,G=Nμ=E-TS+PVから

E=Nμ+TS-PVなので(∂E/∂N)S,V=μ

ですからdE=TdS+PdV+μdNとかける

からです。

 

このdG=-SdT+VdP+μdNと,G=Nμからの

G=Ndμ+μdNの右辺同士を等置して,

Ndμ+SdT+VdP=0 を得ます

これはGibbs-Duhemの関係として知られています。

 

一方,G=E-TS+PV,Nμ/T=E/T-S+PV/T

から,Nd(μ/T)+μdN/T=-(E/T)dT-dE/T

-dS+d(PV/T)です。

 

dE/T=dS+PdV+μdN/T,より 

結局,d(PV/T)=(E/T)dT+(P/T)dV

+Nd(μ/T)が得られます。

 

ここで,β=1/(kT)より,1/T=kβ,

-dT/T=kdβを代入すると,

(1/k)d(PV/T)=-Edβ+βPdV

+Nd(βμ)ですが,

これを,先に求めた

(inZ)=-<E>dβ+βPdV

+<N>d(lnλ)と比較し,統計平均<E>,<N>

それぞれE,Nと同一視します。

 

 したがって,PV=kTinZ,かつ,

lnλ=βμ=μ/(kT),つまり,λ=exp{μ/(kT)}と置けば

いいことがわかります。

 

λはフガシティー(fugacity)と呼ばれ, Zは大分配関数,

または,大きな状態和と呼ばれています。

 

 先に得た系がエネルギーE,粒子数Nの状態を占める確率分布 

λexp(-βE)/[Σi.jλNiexp(-βEij)は,最終的に

exp{-(E-μN)/(kT)}/[ΣN,Eexp{-(E-μN)/(kT)}]

となることがわかりました。

 

再確認しますが,これは系が温度Tの熱浴に接していて,さらに

化学ポテンシャルμ の貯槽と粒子交換するときの確率分布です。

 

さて,粒子がNに固定されていて,エネルギーのみ交換可能な

正準集団での分配関数Z=Σexp{-E/(kT)}はΣの中に

粒子総数がNであるという情報が含まれているため,V,T,N

の関数ですから,Z=Z(V,T,N)と書くことにします。

 

一方,大分配関数は,Z=ΣΣλexp{-E/(kT)}

で与えられますから,=Σ=0λ(V,T,N)と書くこと

ができます。はV,T,N,λの関数です。

 

そこで,PV=kTinZによってPをV,T,Nの関数として

表現したいときには,λを消去する必要があります。

 

これは<N>=λ[∂(inZ)/∂λ]V.Tを用いて,<N>を

Nと同一視することで可能です。

 

例として,先に求めたN個の単原子分子からなる理想気体の系

を考察します。

 

前回,単原子分子N個の理想気体の分配関数Zが, 

Z=(1/a){∫d∫exp(-2/(2mkT))d}

=V(2πmkT)3N/2/a で与えられることを見ました。

 

ここで6N次元相空間の細胞の測度であるaは6次元空間

の測度bによってa=bで与えられること,そして矛盾の

ない示量関数であるためにはN!で割る必要があることを

考慮して,Z(V,T,N1={V(2πmkT)3/2}/(bN!) 

とすれば,この理想気体の大分配関数:

=Z(V,T,N,λ)は,

 =Σ=0λ(V,T,N)

=ΣN=0(λ/b){V(2πmkT)3/2}/(N!) 

exp{λV(2πmkT)3/2/b} で与えられることが

わかります。

 

そこで,lnZ=λV(2πmkT)3/2/bですから. 

<N>=λ[∂(inZ)/∂λ]V.T=lnZ=PV/(kT)

です。

 

これから,状態方程式:PV=<N>kTが得られます。 

ちなみに内部エネルギーの平均値は粒子数=N=一定の系で 

<E>=-[∂(inZ)/∂β].Vですから,

自由度が3の単原原子分子の理想気体では,

Z=Z(V,T,N1={V(2πm/β)3/2}/(bN!)より 

<E>=(3/2)NkT) です。

 

なお,後述するように,量子論ではZを用いるよりZを用いた

方が便利なことが多いようです。

 

λの代わりにμを用いたとして,Z=exp{―βΩ(V,T,μ)}

によって,熱力学ポテンシャルΩ=Ω(V,T,μ)を定義します。

 

Ω=-lnZ/β=-kTinZなので,Ω=-PVです

から,dΩ=-PdV-VdP,ですが,Gibbs-Duhemの関係:

Ndμ+SdT-VdP=0 を代入し,Nを粒子数平均値

<N>に置き換えると,dΩ==SdT-PdV-<N>dμ

を得ます。

 

一方,G=<N>μ=<E>-TS+PV=<E>-TS-Ω

と表現されるので粒子数Nの変動を許す一般の系での

内部エネルギーの平均値は,粒子数平均値<N>を用いて,

<E>=Ω+μ<N>+TS と表現されます。

 

そこで,T→ 0の極低温の極限では,<E>=Ω+μ<N>

となります。

 

(※閑話休題) 論議に必要な知見である熱力学変数の話をします。

 

1. 示強変数,示量変数,Legendre変換 

熱力学では,(熱平衡にある)現在の状態の巨視的量だけで定まる

物理量を状態量と呼び,これを他の状態量の関数と考えたものを

状態関数,または熱力学関数(熱力学ポテンシャル)といいます。

 

,関数Fを状態変数X,Y,Z,..を独立変数とする状態関数

とします。すなわち,F=F(X,Y,Z,..)とします。

 

すると,dF=xdX+ydY+zdZ+..と書けます。

このとき,x=(∂F/∂X)Y,Z,. etc.ですが,係数x,y,z..

をそれぞれ,X,Y,Z,..に共役な変数と呼びます。

 

新しい関数FをF=F-xXで定義すると,

dF=-XdX+ydY+zdZ+..となります。

こうした変換をLegendre変換と呼びます。

 

これは変換で得られる状態量Fでは, 独立変数が

Fのそれ:X,Y,Z,..から,Xだけがその共役xに代わり,

の独立変数はx,Y,Z,..になること:

=F(x,Y,Z,..)と書けることを意味します。

 

例えば,1種類の物質のみから成る構成分子数,またはモル数

が不変な閉鎖系の内部エネルギーを,慣例に従ってEではなくU

で定義すると,可逆過程(準静的過程)での熱力学第一法則は,

dU=TdS-PdVです。

 

これは,UがS,Vを独立変数とする2変数の関数U=U(S,V)

であることを意味しています。

 

そこで,エンタルピーという状態量 Hを,H=U+PVで定義

すると,dH=TdS+VdPとなり,Helmholtzの自由エネルギー

FをF=U-TSで定義すると.dF=-SdT-PdVです。

 

さらにGibbsの自由エネルギーGをG=F+PVとすると, 

dG=-SdT+VdP となります。

 

これから,H=H(S,P),F=F(T,V),G=G(T,P)と

書けて独立変数の数はどれも2で不変ですが,種々の変数の

関数に変換できることがわかります。

 

こうしたLegendre変換だけではなく,例えば

dU=TdS-PdVから,dS=(1/T)dU+(P/T)dV

により,S=S(U,V)とも書けます。

 

系がN個の分子,またはnモル(n=N/N0)から成る場合,

U,S,V etc.は物質量に比例する値をとる示漁変数ですから

U=nu,S=ns,V=nvで1モル当たりのU,S,Vである

u,s,vを定義すると,第一法則のdU=TdS-PdVは 

 ds=Tds-Pdv となります。

 

このケースでは,全微分式の係数である温度Tと圧力Pは物質量n

には無関係な変数ですから,これらを示強変数と呼びます。

 

一般のdF=xdX+ydY+zdZ+..の場合でも,両辺をnで

割ると(F/n)=xdX/n+ydY/n+zdZ、n+..となり,

Fが示量変数ならdF=xdX+ydY+zdZ+..の右辺各項

のxdX,ydY,zdZ..もdFと同様,それぞれがnに比例

するはずです。

 

それ故,例えば,互いに共役なxとXとは,一方が示量変数で

他方は示強変数です。つまり,示量変数と示強変数は1対1に

対応します。

 

したがって,任意の状態量はLegengre変換などにより,示強変数

のみを独立変数とする関数となるように変換可能であることが

わかります。

 

特に系の物質成分か1種類だけであって,かつ,気体か液体か固体

のどれか1つの相にあるケースであれば,その系での独立変数を

示強変数T.Pのみとして,他の量は全てこの2変数の関数で表

わされる従属変数と考えることができます。

 

2.Gibbsの相律 

対象とする系の温度がT,圧力がPであってc種の成分,

φ種の相から成り,相間の移動反応を除いてr個の区別

できる化学反応が起こっているとします。

 

この場合,系が平衡にあるときには,「幾つの示強状態変数を

任意に定めることができるのか?(独立な変数はいくつなのか?)」

という問題を考察します。

 

任意に定め得る示強変数の数を系の自由度,または可変度

(variance)といいますが,ここでは,系を構成する各物質成分

の量を含む示量変数は問題にせず.ただ.示強変数のみを問題

にします。

 

それらを書き上げると,温度,圧力,および,組成変数が全てです

が,それは,,P,x11,x21,..,x1,x12,x22,..,x2,..,x1φ,

2φ,..,xφと書けます。その数は,(2+cφ)個です。

 

ここで,xα(k=1,2,..,c;α=1,2,..,φ)は,物質成分kの相α

にある分率です。

 

それ故,これらの変数全てが独立というわけではなく,α=1,2,..,φ

の各相αにおいて,恒等的関係:Σk=1α=1があります。

 

 次に,相平衡にある条件として,各成分kの各相への分配を規定する

c個の条件μ1=μ2..=μφ,(k=1,2,..,c)があります。

 

これは,組成分率を=(x11,..,x1,..,x1φ,..,xφ)と略記

すると化学ポテンシャルのそれぞれがT,P,の関数であること

つまり,μα=μα(T,P,)なることが陰に含まれていて,

各成分kについて(φ-1)個の独立な条件式です。

 

 また,系内で生じるr個の化学反応が,ある相α1で生じる際の

 平衡条件として,j=-Σ=1jνijμα1 (j=1,2,..,r)なる

 r個の関係式があります。ただし,jは化学反応jの親和力,

 νjkはその反応jでの成分kに対する化学量論係数です。

 

以上,3種類の条件は,全部で{φ+c(φ-1)+r}個の独立な関係式

となっています。

 

 したがって,最初に挙げた(2+cφ)個の示強変数のうちで独立な

変数の数が任意に取りえる変数,つまり自由度であり,それをfと

すると,f=(2+cφ)-{φ+c(φ-1)+r},すなわち,自由度は

f=2+(c-r)-φで与えられることになります。

 

特に,化学反応が全く無いことが保証されている全く熱力学的な

系ならr=0 であって,この場合はf=2+c-φとなります。

こうした関係を「Gibbsの相律」と呼びます。

 

 この関係は,次のような方法で導くこともできます。

 

まず,Gibbs自由エネルギー G=F+PV=E-TS+PVの

微分を取ると,dG=dE-TdS-SdT+PdV+VdP

ですが,c成分がΦ相で混合している場合,各相αにおいて,

G=Gα=Σαμαであり,

dGα=dEα-TαdSα-SαdT+PαdV+VdPα,

かつ,dEα=TαdSα+PαdV(∂Eα/∂N)S,VdNα

なのでdGα=-SαdTα+VdPαμαdNα

となります。

 

一方, Gα=Σαμαなので,

dGα=Σαdμα+ΣμαdNα 

 多成分に一般化されたGibbs-Duhemの関係式として,

 SαdTα-VαdPα+Σαdμα0 が

 得られます。

 

これは,分率の定義:xα=Nα/Nα; Nα=Σα

により,両辺をNαで割ると示量変数の体積Vαやエントロピ-

αの1分子当たりの量をVα,Sαとして,

Σαdμα+SαdTα-VαdPα=0

と表現されます。

 

 平衡条件:Tα=T,Pα=P,μα=μを考慮すると,

示強変数としては,P,μ(k=1,2,..,c)の(2+c)個だ

けで十分です。

 

これとαdμα+SαdTα-VαdPα=0,または

これにα=T,Pα=P,μα=μを代入した式である 

Σαdμ+SαdT-VαdP=0 (α=1,2,..,Φ)の

Φ個の制約条件から,結局,独立変数の個数がf=2+c-Φである

という相律を再び得ます。

 

ところで,今の(T,P,)系では,組成以外の独立変数をT,P

の2つだけと仮定しましたが,例えば他の示強変数として表面張力

γや磁場(または)などを考慮する必要がある場合もあり,

f=2+c-φの右辺の2は3にも4にも変わることああると

考えられます。 (閑話終わり※)

 

 自分の過去ノートの日付によると,統計力学は1995年から1997年

(45歳から46歳),熱力学は遅れること1年の1997年(47歳)に集中して

やったようです。

 

もちろん私の場合,一応,1年浪人して19歳(1969年)で大学の理学部

物理学科に入学して,27歳(1977年)の4月に普通の会社に就職する

までは,曲りなりにも理科系の学生をやっていたので,

理の基礎分野では門前の小僧的知識があり,全く白紙状態から

独学スタートしたのではなく,そうした40歳代での専念行為は知見

を完全にしたいという趣味的欲求によるものです。

 

年を取ってからも,別に資格取得目的でも外国語会話習得でも

なく,何の現実的利益にもならないだろう趣味としての学問を

やるというのも,オタク的変態のようですが一応,人生最後の

ライフワークです。

 

仕事(Work)であって労働(Labor)じゃないので利益になりません。。

分をわきまえず,脳天気で苦労知らずのオボッチャマの趣味に

興じてるキリギリス,すぐにも7飢え死にするヨ。>自分

本文が短かいと思ったので熱力学や余談を書いたら,いつもより

長くなったので終わります。

 

 相変わらず,金欠病です。昔から,お金を使う趣味が全く不可能

 なときには幸か不幸か?仕方なく筆が進みますね。

 

参考文献:阿部流蔵 著「統計力学(第2版)」(東京大学出版会) 

妹尾学 著「熱力学」(サイエンス社)

| | コメント (0) | トラックバック (0)

2014年11月11日 (火)

統計力学の基礎(2)(古典統計力学2)

古典統計力学の続きです。

前回,粒子数がNで自由度がfの多粒子系で,系のエネルギーが

の孤立系において,エネルギーがε,..,にある粒子数を

それぞれ,n,n,..,とすると,その熱平衡時の分布がM-B分布

つまり,n=fexp(-βε)/Z;ただし,Z=Σexp(-βε)

で与えられることを見ました。

 

しかし,まだ,パラメータβやZがどのような量であるか?が不明

です。

 

そこで,以下では,まず,それらの量を明確にしたいと思います。

 

,異なる種類の粒子からなる2つの系A,Bがあって,これらの系

はどこかの境界で接触していて,それ以外の境界からはエネルギー

の流出入がなく熱平衡にあるとします。

 

このとき,系A,Bそれぞれを構成する粒子数N,Nは,独立に

不変ですが,エネルギーは弱い接触相互作用ですがAとBの間

では自由に交換できます。

 

系AのN個の粒子のうちで,エネルギーがε,..,にある

ものの数が,n,.., で.系BのN個のうちエネルギーが

,e,..,にあるもの数がm,m,..,であるような粒子

微視状態の配分方法の総数をWとすると,

 

W={N!/(n!n!..), }×{N!/(m!m!..)}

と書けます。

 

そこで,Stirlingの公式によって, 

lnW ~ fln(f)-Σln(n)+fln(f)

-Σln(m )

と近似されます。

 

前の考察と同じように=N=N,かつ,

Σε+Σme=E の条件下で,lnWを最大にする

分布をLagrangeの未定係数法で[求めます。

 

すなわち,-δlnW+αΣδn+αΣδm

+β(Σεδn+Σδm=0 ,δn,δm

ついての恒等式になる条件を求めます。

 

そこで,δn,δmの全ての係数がゼロ,つまり, 

lnn+α+βε=0,かつ,lnm+α+βe=0 

です。

 結局,n=Nexp(-βε)/Z,かつ,

=Nexp(-βe)/Zなる表式を得ます

 ただし,Z=Σexp(-βε),Z=Σexp(-βe)です。

 

ここでβは系AとB共通の係数であることに注意します。

 

熱力学によれば,こうした接触した2つの系AとBが熱平衡状態

にある条件は,両者の絶対温度が等しいことです。

 

それ故,βは絶対温度Tのみの関数と考えられます。

 

N粒子から成る孤立系でのM-B分布:n=Nexp(-βε)/Z; 

Z=Σexp(-βε)に戻って,lnZのβによる偏微分を取ると, 

(lnZ)/∂β=∂Z/∂β/Z=Σεexp(-βε)/Z 

=-Σε/N:つまり,∂(lnZ)/∂β=-E/N です。

 ところで,系のHelmholtz自由エネルギーをFとすると,Fの定義

系の内部エネルギーをU,エントロピーをSとしてF=U-TS

ですが粒子の総ネルギーEが系の内部エネルギーUに相当するため,

F=E-TS,/T=E/T-S です。

 

そこで,{∂(F/T)/∂T}

=-E/T+(∂E/∂T)/T-(∂S/∂T) 

ですが,熱力学第一法則により,TdS=dE+PdV

なので,(∂S/∂T)=(∂E/∂T)です。

 

したがって,{∂(F/T)/∂T}=-E/Tを得ます。

 

一方,∂(lnZ)/∂β=-E/fから, 

(flnZ)/∂T={∂(flnZ)/∂β}(dβ/dT) 

=-E(dβ/dT) です。

 

それ故,{∂(F/T)/∂T}

={∂(flnZ)/∂T}{(dβ/dT)-1/T2 ですから,

ある定数をkとして,F/T=-fklnZ,(dβ/dT)T=-1/k 

と置けば,この関係式が成立します。

 

dβ/dT=-1/(kT)を積分すれば,β=1/(kT)+C

を得ますが,特に積分定数Cはゼロでβ=1/(kT)と考えられます。

 

何故なら,非常に高温(T → ∞)では分布は,エネルギーεの値に

関係なく一様になると思われるからです。

 

一方,lnW=flnf-Σlnnに,n=fexp(-βε)/Z,

および,lnn=lnf-βε-ln/Zを代入すると,

f=Σ,E=Σεによって,lnW=βE-flnZ

ですが,F=-fkTlnZより,flnz=-βFなので, 

lnW==βE-βF=βTS=S/kです。

故にS=klnWなる関係式が得られます。

 

ところで,理想気体の系でしかもそのN個の構成分子が構造を

持たない質量mの質点と見なせる場合(例えば単原子分子の場合),

個々の粒子のHamiltonianはH=/(2m)で与えられ,系の自由度

はf=3Nです。

 

このとき,運動量がの粒子数は,n()=Nexp{-/(2mkT)}

で与えられるはずでますが,このN粒子を閉じ込めている容器の体積

Vの増減と共に,取り得る状態の配分数WはVに比例して増減する

はずなので比例定数をAとしてlnW=ln(AV)=Nln(AV)と書け

ます。

 したがって,S=klnW=Nkln(AV)より,

T(∂S/∂V)=NkT/Vを得ます。

ところが,TdS=dE+PdVなので,

P=(∂S/∂V)です。

 

以上から,P=NkT/V,つまりPV=NkTなる関係式が

得られました。

 

理想気体の状態方程式は系の絶対温度をT,圧力をP,体積

をVとするとRを気体定数として,PV=nRT=Nk

で与えられます。

 

ただし,nはこの気体のモル数です。1モル当たりの気体分子

であるAvigadro数(~6.02×1023)をN0とすると,n=N/N0であり

はBoltzmann定数で,k=R/N0で定義されます。

 

この状態方程式PV=NkTを,すぐ前で求めた式:PV=NkT

比較等置して,未知定数 kはk=kB と結論されます。

 

こうして,結局,β=1/(kT)で,N個粒子の孤立系の粒子数分布は 

=Nexp{-ε/(kT)/ZでZ=Σexp{-ε/(kT)であり,

lnZはHelmholtzの自由エネルギーFとF=-NkTlnZなる関係

にあります。

 

そうして,エントロピーはS=klnWで与えられることが

わかります。 

この最後のエントロピーと状態配分の関係式はBoltzmanの原理

(Boltzmannの関係)と呼ばれています。

 

この関係から,熱力学での熱平衡でのエントロピー極大の条件

が正に最大確率の分布に対応していることがわかります。

 

 さて,次に,正準集団の方法を紹介します。

 

体系の全エネルギーが各粒子のエネルギーの総和となることは,

運動エネルギーにおいては常に成立しますが,位置エネルギーに

おいては一般には成立しません。

 

以下ではこれまでと異なって粒子間に微弱とはいえない相互作用

がある場合をも扱える一般的方法について述べます。

 

 前と同じく莫大な個数N個の粒子から成り立つ全粒子系の自由度

がfの体系を考え,その一般化座標を(q1,..,q),共役運動量

=(p1,..,p)とします。

 

この体系の状態を示す.座標が(,)の代表点の作るEuclid空間

は前にΓ空間と呼んだ2f次元の相空間ですが,

今回はH(,)=E(一定)という制約は設けません。

 この系と同じ構造を持つ体系を非常に多く用意できたと想定し

それらを並べるとします。

 そしてその同じ系の総数はM個とします。MもAvogadro数に

匹敵するような巨大な数です。

これらM個の体系はその間に微弱な相互作用があって,相互に

エネルギーを交換ができるとします。

 

しかし,M個全体では外界と何の交渉もない孤立系であるとします。

 

このように体系の集団(ensemble)を想像し,その統計的平均が実測

の物理量に一致すると考えます。

 

そしてM個の同じ体系の全体をM×(2f)次元の相空間として一般

化座標を((1,(1),21,(2),..,(M),(M))とする代表点

全体の空間をΓ0空間と呼ぶことにします。

 

ただし,(j)=(q(j)1,..,q(j)),(j)=(p(j)1,..,p(j)), 

(j=1,2,..,M) とします。

 

Γ空間と空間とΓ0空間の関係は,前回の孤立系でのμ空間とΓ空間

の関係と同じです。

 

そこで,以下,このM個の系全体の作る孤立英について前回と同じ 

小正準集団の方法を適用します。

 

すなわち,Γ空間を2f次元体積がaの高々可算個の細胞に分割

し,k番目の細胞のエネルギーをEとします。

M個の系の集合の要素のうち,このk番目の状態を取るものの個数

をMとします。

 

そして,M個の系全体の不変なエネルギーをE0=Σ

します。

 

 M個をM,M,..,M,..個に分割する配置の数Wは, 

 W=M!/( M!M!..,M!..)で与えられます。

 1つの分け方:(M,M,..,M,..)がΓ0空間のaの体積

の細胞に相当するので,あらゆる(M,M,..,M,..)の組を

与えるΓ0空間内の体積はWaとなります。

 

「エルゴード仮設」,または「等重率の原理」によってΓ0空間内

のある体積を占める確率はその体積に比例します。

 

それ故,今の場合Σk=1=M,および,E0=Σの条件下

でlnWを最大にする分布が最も実現確率が高く,熱平衡時にはそうした(M,M,..,M,..)が実現されているはずです。

 

このMを求めると前回と同じ手順でM=M exp(-βE)/Z

なる分布が得られます。

 

 これをM=Σに代入すると, 

 Z=Σexp(-βE)

 =(1/a)∫exp(-βE)dq..dqdp..dp 

が得られます。

 このZを分配関数(partition function),または,

状態和(sum over states)と呼びます。

 前と同じくLagrangeの未定係数法で仮に与えた係数パラメータ

βは,熱平衡にあるM個の体系に共通なので.絶対温度Tと同じく

熱平衡の尺度と予期されます。

 いいかえると,βはTのみの関数です。

 そして,今.この想定された莫大な個数M個の系のうちで唯一つ

の系だけに着目すると,それがエネルギーEの状態となる確率が

exp(-βE)に比例することがわかります。

 

 実際の正しい確率はexp(-βE)/{Σexp(-βE)z}で

与えられます。

 

着目した体系以外の,それに接触した残りの巨大な(M-1)個の

系全体を温度を一定に保つ熱浴と見なせば,これは熱浴に漬かって

いる,または,恒温槽に接触している体系の確率分布を記述するもの

と考えられます。

 

 このような[分布を正準分布(canonical distribution)と呼び,

この方法を正準集団の方法といいます。

 

 前と同様,∂(lnZ)/∂β=-Σexp(-βε)/Z

 =-E0/M=<E>です。

 

 一方,Helmholtzの自由エネルギーをF=U-TSとすると,

 [∂{(F/T)}/∂T]V=-U/Tなので,内部エネルギー

 Uが系のエネルギーの平均値 <E>に一致する,

 つまり,U=<E>として比較することから,kをある定数

としてF=-kTlnZ,β=1/(kT)が得られます。

 

 さて,N粒子の理想気体ではf=3Nで

 Z=(1/a)∫exp(-βE)dq..dqdp..dp 

 =(1/a){∫d∫exp(-2/(2mkT))d}

 =V(2πんkT)3N/2/aなので,

 

 F=-kTlnZ=kTln{V(2πmkT)3N/2/a}

となります。

 そこで,(∂F/∂V)=-NkT/Vですが,F=U-TS

なので(∂F/∂V)=-(∂S/∂V)=-Pです。

 

 すなわち,PV=NkTを得ますから,やはり,k=k

であり,結局,F=-kTlnZ,β=1/(kT) です。

 

 しかし,実は上記のFの表式:

F=-kTlnZ=kTln{V(2πんkT)3N/2/a}は体積V 

への依存性が正しくないです。

つまり,N/V一定でV→2VならF→2Fとなるべきなのにそう

なっていません。

 

 熱力学量には温度T,圧力Pのように物質量とは無関係な示強性

の量と,エネルギーやエントロピーのように物質量に比例して1モル

が2モルになればそれらも倍の値になる,という示量性の量があり

ますが自由エネルギーFは1つの示量性の量です。

 

 N/Vが一定でV→ 2VというのはN→ 2NでVV→ 2Vで

あり,単純に同じ系が2つあるのと同じ意味ですから,常識で

考えてもF→ 2Fとなるはずです。

 

 つまり,F=-kln{V(2πmkT)3N/2/a}

 =-NkTln{V(2πmkT)3/2/a/(1/N)} 

→ -kln{(2V)2N(2πmkT)3N/a}

=2NkTln{2V(2πmkT)3/2/a1{1/(2N)}} 

2F-2NkTln2-kTlnaとなって常識と相容れません。

 

 仮に,Z=(1/a){∫d∫exp(-2/(2mkT))d}

 =V(2πんkT)3N/2/aの右辺をN!で割れば,

 -ln(N)=-NlnN+Nであり,N→ 2Nに対して, 

-NlnN+N→ -2NlnN+2N=2(-NlnN+N)-Nln2な

ので,ほぼこの矛盾が解消されます。

 

 これは古典統計の1つの欠陥であり,N個の粒子に1番,2番

と順序を付けて区別することはできない,という量子統計で解決

されることは後述しますが,量子統計に移行しなければならない

理由の1つです。

 

 前回の孤立系では,N個の粒子の個々のエネルギーが

ε,..,εE=ε+ε..+εと書けるとき,

Eが一定不変という条件つきでk番目の粒子がエネルギーε

を取る確率がexp{-ε/(kT)}/[Σjexp{-ε/(kT)}]

で与えられることを見ました。

 

 ここではN個の粒子系のエネルギーがE=ε+ε..+ε

を満たすエネルギー値Eを取る確率が

exp{-E/(kT)}/[Σexp{-E/(kT)}]

exp{-(Eε+ε+..+ε)/(kT)}

/[Σexp{-(Σε)/(kT)}]で与えられ,

, その際,E=ε+ε+..+ε=(一定)という制約条件

なしで各々の粒子のエネルギ^-εが独立に変わり得ると

してよいので,

 結局,k番目の粒子がエネルギーεを取る確率が

 exp{-ε/(kT)}/[Σjexp{-ε/(kT)}]で与えられる

という同じ結論が,孤立系で外界とエネルギーのやりとりがない

という条件抜きでも成立することがわかります。

 

したがって,今回の正準集団の方法による結果は前回の小正準集団

の方法を包摂したより一般的なものであることがわかりました。

 

 今日も次の大正準集団の項に入る予定でしたが,長くなり過ぎた

のでここで終わります。

 

参考文献: 中村 伝 著「統計力学」(岩波全書) 

阿部流蔵 著「統計力学(第2版)」(東京大学出版会), 

久保亮五 著「大学演習 熱学・統計力学」(裳華房)

| | コメント (0) | トラックバック (0)

2014年11月 5日 (水)

統計力学の基礎(1)(古典統計力学1)

 本ブログでは,統計力学についてこれまでは関連した断片的

トピックと,その応用例に言及する程度の記事しか書いてき

せんでした。


(↑※例えば,熱力学も含めると,一番最初のは2006年4/8の記事

基礎物理学講義①(温度と熱)」これ,昔池袋の専門学校での講義

用に作ったノートを書き写しただけです。

 次は2006年522の「エルゴード問題と次元」です。

また,2006年7/5の記事

可逆と不可逆のはざ間(エントロピー増大則)

や2006年8/6の「エントロピーの定義」など Pending)


 しかし,
今回,場理論による物性理論の扱いに関連した松原

温度Green関数についての記事を書くに当たって,急遽,統計力学

のエッセンスを紹介する気になりました。


 統計力学は,気体,液体,固体のように日常生活で我々が

いる物体を物理的対象と考え,そうした個々の物体を熱力学

的体系,または系と呼ぶとき,分子論に基づいて系を莫大な個数

の微視的粒子の集団と見るとき,それら粒子群の力学的運動の

総体から圧力や温度といった熱力学の巨視的状態量を説明する

ために生じた学問です。


 19世紀に分子論という思想が発生した当時,熱力学においては
,

熱とは何か?温度とは何か?などの根源的疑問がありました

,これらの解答が統計力学から得られることがと期待されて

いました。


当時,既に熱力学において熱はエネルギーの1形態であり温度

系の熱平衡を特徴付ける尺度であることなどは解明されて

いました。


普通に,大学で学ぶ統計力学というのは,熱平衡での系の温度

エネルギーがEのとき系の状態の出現確率はBoltzmann因子,

または,Gibbs因子と呼ばれるexp[-E/(kT)]なる量に比例

する,という法則を認めさえすれば,非常に便利な計算道具と

なります。

 
それ故,実験観測値を精密に計算で再現できる場の量子論で

くりこみ理論などと同じく.謂わゆる対症療法=有効理論

いう見地では,それほど難解であるという印象を受けないか

知れません。

 
しかし,私のように計算ができる云々より何故そうした方法

で計算可能なのか?という原理の方に興味を抱くような科学

的には変態?思想の持ち主にとっては,ずいぶん入り口の敷居

が高いものでした。

 
先達たちには,とりあえず,それを使って種々の必要な計算を

実施し年月と経験を経るうち自然に原理的なことも理解できる

ようになる,といった内容の忠告がしばしば成されましたが,

とにかく,理由がわからなければ一歩も先に進めないという,

厄介で意固地な性格でしたからね。。。

 
だからこそ,先端の研究課題までサクサクとは進めず,未だ

地べたを這いずりまわってるのでしょうネ。。>自分

 

というわけで,私のように大学生当時は学生運動に明け暮れたり,

むしろ数学の方に魅力を感じていて物理学では劣等生であり,

物理については,その頃講義で履修してもよく理解できず,後に

必要性を感じて独学で勉強を始めたものがほとんどなので,統計

力学はそうした頃から初学者であった私には理解に苦労しました。

 

どの本を読んでも常に,理論の導入部で躓いて迷路に入ること

の繰り返しで,10冊以上の専門書の精読と挫折の遍歴があり,

結局,中村伝著「統計力学」(岩波全書)に至って一応の納得を

得てゴールしたのはもう四十路半ばでした。。

 

熱力学については未だに遍歴中ですが。。

イヤ,馬鹿は死んでも。。かな?

 

専門書ではないですが,「磁力と重力の発見」(岩波書店)で最近

大仏次郎賞をとられたらしい山本義隆氏の恐らく最初の物理学史

のヒット著作である「熱学思想の史的展開」(現代思潮社)や

高林武彦氏の「熱学歴史」海鳴社)などは大いに理解の助けに

なりました。

 

また余談ですが,アチラは覚えてないでしょうが,山本義隆氏には

氏が東大全共闘議長だった頃,例えば1970年6月の安保粉砕闘争

のデモで当時封鎖占拠していた東大地震研の教室に,S大全共闘

一員としてゴロ寝一泊した際などにお会いしたこともありま

したネ。

 

さて,本文です。

系の熱平衡状態において観測される巨視的量の値が,微視的には

莫大な個数の運動する粒子から成る系の個々の粒子の物理量の

総和の長時間平均に一致するはずである,というのが統計力学の

概念であると考えられます。

 

そして,統計力学には「エルゴード仮設(仮説)」というものが

あります。

 

これは,古典的には粒子は非常に長時間には取り得る全ての軌道

を同じ回数だけ取るため,この全軌道平均(相空間平均)が長時間

平均に一致するだろう,というものです。

また,量子論においても

上記の軌道を量子状態に置き換えるとそのまま,同じことがいえ

るという仮説です。

 

これは,この多粒子系の一般座標と共役運動量の作る相空間の

領域に系を構成する粒子群が存在する確率が単純にその領域の

相空間の体積に比例するというような確率原理を仮定して相空間

平均を計算する方,長時間平均を計算するより,はるかに容易な

ため,長時間平均の代わりに便利な相空間平均を用いて計算して

よいことを保証するものです。

 

ここで上述の,「領域に存在する確率が,その領域の相空間の体積

に比例する。」という確率原理は,等確率の原理」あるいは

「等重率の原理と呼ばれています。

 

書物によっては相空間でなく位相空間という表現をとっている

ようですが,ここで相空間と呼んでいるのは,ここでの位相空間(phasespace)=相空間が,呼称は同じでも全く無関係な数学(幾何学)

での位相空間(topological space)と混同されるのを避ける

ためです。

 

系の粒子数NとエネルギーEが定まっていて外界と全くやりとり

のない孤立系では,系のエネルギーは全粒子のエネルギーの総和

で与えられます。


この全N-粒子系の自由度が
fの場合の一般化座標を

(q1,..,q),とし,その共役運動量=(p1,..,p)

します。

 

この多粒子系に対する系全体のHamiltonianをH(,)

とすると孤立系の条件はH(,)=E(一定)です。

 

そして,長時間では,あらゆる軌道を同じ回数だけ取ることから,

N粒子の自由度fの一般化座標,と共役運動量の組()

を2f次元Euclid空間の位置座標と見たときの空間を相空間,

特にΓ空間と呼び,()で表わされる状態をその代表点と

呼びます。

 

個々の粒子が構造を持たない質点と見なせるなら各々は3

つの空間座標という自由度しか持たないのでN個の粒子の

場合,f=3Nで,2f=6NなのでΓ空間6N次元空間です。

 

しかし,一般には個々の粒子(分子)は構造を持っているので,

必ずしも2f=6Nとは限りません。


いずれにしても,H(,)=Eは2f次元の空間における

等エネルギー面を示す曲面に過ぎないので,その次元は

(2f-1)であり,その上でのいかなる状態集合の領域も 

2f次元の体積という意味での測度はゼロです。

 

したがって,E≦H(,)≦E+ΔEなる微小な厚さが

ある領域を想定すれば,これは2f次元の体積を持つ領域です。

 

ところで,何故,粒子数の他の系に対する制約がエネルギーE

に対するものだけで運動量や角運動量などその他の力学変量

を問題にしないのか?というと,

 

統計力学で通常扱う系では系を構成する多粒子群は箱など

容器の有限の領域内にと閉じ込められているケースがほとんど

,もはや空間の一様性や等方性は失われ,それらに伴なう

運動量や角運動量は平衡状態では保存量でないからです。

残っている対称性はは時間の一様性だけで,これに基づく

エネルギーの保存があるだけです。

 

さて,代表点の位置ベクトルを=,()と置けば3次元

空間のベクトル解析における公式:ΔE=ΔH=∇H・Δ

が成立するはずです。

 

ただし,ここでの一般化された勾配 ∇H=gradHは, 

∇H

(∂H/∂q1,..∂H/∂q,∂H/∂p1,..∂H/∂p) 

で与えられます。

 

そして,ΔE=∇H・Δから,ΔE=|∇H||Δ|

より,|=ΔE/|ΔH|と書けます。

 

また,等エネルギー面H(,)=E上の代表点(,)に

おけるこの面の法線単位ベクトルは∇H/|∇H|であること

がわかります。

 

したがって,Γ空間の等エネルギー面H(,)=E上の面積

がSの領域を法線方向にΔだけ切り取った,

E≦H(,)≦E+ΔEの部分の体積は,

|SΔE/|ΔH|となります。

 

ここで,SもΔEも任意に取った量ですから,それらを共に単位量

1に取ることもできます。また,等重率の原理により代表点が

この領域に存在する重み確率は一様でSとΔEの双方に単純に

比例しますから,結局,存在確率は1/|ΔH|に比例することが

わかります。

 

また,=,()により,代表点の速度を

=d/dt=(,p)と定義すれば, 

(,p)

=(∂H/∂p1,..∂H/∂p,-∂H/∂q1,..-∂H/∂q) 

です。

 

故に,

=div

=Σj,kH/∂q∂p-∂H∂p∂pq=0

を得ます。つまり,代表点の集合を流体と見たときの流速

の発散がゼロです。

 

一方,この代表点から成る流体の運動においてどこにも

湧出し,吸込みがないとして連続の方程式:

∂ρ/∂t+∇(ρ)=0 が成立します。

 

よって代表点密度のLagrange微分(流体素片の運動に沿って

の微分)はDρ/Dt=∂ρ/∂t+∇ρ=∂ρ/∂t+∇(ρ)

=0 となります。

 

したがって,代表点の密度はその運動と共に時間tに

依らず不変であることが示されました。

これはリウヴィル(Liouville)の定理と呼ばれています。

この方向で話を進めていく方法もありますが,等重率の原理

を利用して別のカラメ手から熱平衡での粒子の状態分布を

求めてみます。

 

簡単に,自由粒子が体積Vの箱の中にN個閉じ込められている

ケースを考えます。(※実際には互いに衝突したり分子であれば

分子間力などの微弱であっても相互作用があるため,もはや本来

の自由粒子ではありませんが。)

 

また,簡単のため,系全体のHamiltonianが系を構成する各粒子

のそれの単純和であるとしてよい場合を論じます。

 

 すなわち,系全体のHamiltonianHは,H(,)

(k=12,..,N)を各粒子のそれとして,

(,)=Σ(,)で与えられ,H(,)=E

のときの全体のエネルギーEは,(,)=ε

満たすεによってE=Σε表わされるとします。

 

さらに,点(,)全体のつくる2f次元Γ空間の1粒子

(または1調和振動子)(kjkj)(k=1,.N)の各々が作る

Γ空間の部分相空間をμ空間と呼ぶことにします。

Γ空間はこれらN個のμ空間の直積です。

 

そして,例えばN個の1次元調和振動子の場合,μ空間は

2自由度(2次元)で(p,q)のHamiltoninは1次元調和振動子:

(p,q)=p/(2m)+mω/2で与えられるとします。

/(2m)+mω/2=εはμ空間の楕円です。

 

ここで,便宜上添字kを省略してμ空間の点の座標を

(q,p)で表わし等エネルギーの楕円を

/(2m)+mω/2=εと書きます。

 

楕円の長軸と短軸は(2mε)1/2,{2ε/(mω)}1/2です

から,この楕円の面積は2πε/ωです。

 

このμ空間の微小部分dqdpに代表点がn個あると

すると,n/Nが1個の振動子がこのdqdpに存在する

確率を表わします。

 

有限次元Euclid空間全体はパラ・コンパクト,または

第二可算的なので,2次元のμ空間を大きさaの細胞

に分けこの高々可算個の細胞群に番号を付けることが

できます。

各μ空間はほとんど独立でε,ほぼ一定ですが,μ空間

同士でほんの僅かなエネルギー交換があって軌道には,少し

のぼやけがありそれは楕円の面積程度ですから,それが細胞

の大きさaに相当します。

 

そこでaのオ-ダーは楕円の面積程度であり,a~ε

ですが,量子論ではこれは不確定性原理:ΔqΔp~hから,

aはPlanck定数h程度と古典論より明確に決まります。

逆に,aの値を固定すると,μ空間での2次元代表点の取る

エネルギーの拡がりはε~aωとkには依存せず共通に

なります。

 

さて,N個の振動子のうちk番目の細胞にn個の

代表点が存在するとします。

 

すると=Nですがこの(n,n,..)の分配の仕方

の総数Wは,W=N!/(n!n2!..)と書けます。

そしてこの分配のΓ空間における体積はWaです。

 

そこで,等重率の原理により,分配が(n,n,..)となる

確率はWに比例します。

 

それ故,平衡時の分布はWが最大になるような分布であると

考えられます。

 

Wが最大になることとlnWが最大になることは等価ですから

計算の簡単さの都合上lnWを最大にする分布:(n,n,..)

を求めることにします。

 

W=N!/(n!n2!..)の両辺の自然対数を取ると

lnW=ln(N!)-Σln(n!)ですが,N,および全ての

常温,常圧で1リットルの気体分子数のオーダー 

1022であるるように系の粒子数N,または自由度f

は莫大な数ですから,Mが大きいときのStirlingの公式:

ln(M!)~M(lnM-1)をln(N!),ln(n!)に適用します。

 

そうすると,近似式として, 

lnW=ln(N!)-Σln(n!)

~N(lnN-1))-Σ(lnn-1)

=N(lnN-Σlnnk を得ます。

 

この式を利用すると,N=Σ,および

E=Σεの制約条件付きlnWの最大値が

Lagrangeの未定係数法を用いて得られます。

 

すなわち,制約条件はδ(Σ)=Σδn=0, 

および,δ(Σε)=Σδnε=0 です。

 

一方,lnWが最大になるための必要条件は 

δlnW=-δ(Σlnn)

=-Σ(lnn+1)δn=0 です。

 

しかし,制約条件があるので,δnを全て任意の値と

することはできず,これらは独立ではありません。

 

したがって,全ての変分δnが独立になるように,

-δlnW+α(δ(Σ)+βδ(Σε)=0  

とします。


すると,
 

Σ(lnn+1)δn+αΣδn+βΣεδn=0

より,Σ(lnn+α+βε)δn=0 ですが,この式では

全ての変分δnは独立なので,係数をゼロとして

lnn+α+βε=0 を得ます。

つまり,n=exp(-α-βε)です。

 

N=Σexp(-α)Σexp(-βε)より,  

Z≡Σexp(-βε)と置けばexp(-α)=N/Z 

であり,=Nexp(-βε)/Z,あるいは,

/N=exp(-βε)/{Σexp(-βε)}

と書けます。


これはエネルギーεを取る粒子の確率が因子

exp(-βε)に比例することを意味します。

 

こうして,古典統計分布であるMaxwell-Boltzmann分布

(M-B分布)が得られました。

 

最初箱の中にN個の自由粒子があると仮定し,途中から

考察の都合上,自由粒子の変わりに1儀元調和振動子という

特殊例を用いたのですが,最終的な分布の導出方法には自由度

関係なくここで得た古典分布は一般的結論であり,自由粒子

しても同じです,上述の1次元調和振動子の場合,系の自由度

f=Nで2f=2Nです。

 

一方,最初に仮定した箱の中にN個の自由質点粒子がある

ケースならf=3Nで2f=6Nであることを断っておきます。

 

粒子数一定で系の外からエネルギーの流出入のない孤立系

においてΓ空間における代表点(微視状態)の集合を小正準

集合,または小正準集団(micro canonical ensemble)と呼び

上述の統計分布はミクロカノニカル分布とも呼びますが,

この方法を小正準集団(集合)の方法といいます。

 

後に紹介する正準集団,大正準集団の方法では,それぞれ,

孤立系ではなくエネルギーの流出入,粒子の流出入を許す

場合の,より一般的な系を想定した方法で同様な統計分布を

見出すことができます。

 

そして,孤立系はそれらの特別な場合に過ぎず,それらに

含まれています。

 

正準集団,大正準集団は,対象とする系と全く同じ条件の系が

大量に合計M個存在すると想定して,系が取る状態の確率は

その集団の中で同じ状態にある系の個数がnならn/Mで与え

られるというような統計集団であるのに対し,小正準集団は

統計集団というよりも状態の集合ですから意味はちょっと

違いますね。

 

今日はここで終ります。

 

悪いクセで脱線しだすと際限がないですが,昔と違って全て

通ってきた道で指針となる覚え書きノートも揃っているので,

すぐ収拾はつきます。

 


参考文献:中村伝 著「統計力学」(岩波全書),

阿部流蔵 著「統計力学(第2版)」(東京大学出版会), 

久保亮五 著「大学演習 熱学・統計力学」(裳華房)

| | コメント (0) | トラックバック (0)

2009年10月11日 (日)

空気分子の大きさ(アインシュタインとブラウン運動)

 地球大気の温室効果について段階的に論じていこうと思います。

 

 まずは,私自身が"自己満足"して納得するため,2006年11/21の記事「地球の平均気温とステファン・ボルツマンの法則」において,温室効果を無視した地球平均気温(=約-18℃)という定量的評価のための材料とした「太陽からの輻射に対するアルベド(albedo:反射率)が30~31%であること」の理論的根拠などから論じてゆきます。

  

 今日は,まず太陽から地球に放射される光が地球表面で散乱され減衰することと関連して,大気圏での主要な散乱体である"主に窒素と酸素から構成された空気分子"の大きさを評価することから始めます。

 

 そのための方法として植物学者ブラウン(Brown)の「花粉の水中運動の永久性=ブラウン運動"の発見(1829)」に対する「アインシュタイン(Einstein)の解明(1905)」の話から始めようと思います。

 媒質中のブラウン粒子の運動を位置座標のx成分の軌道で代表させると,その運動方程式はu≡dx/dtとしてm(du/dt)=X(t)-u/μで与えられると考えられます。

ここにmは対象粒子の質量,X(t)は媒質からその粒子にかかる力のx成分,μは易動度と呼ばれる量で,u/μが媒質の速度に比例する摩擦力となるようなケースでの比例係数の逆数です。

 

m(du/dt)=X(t)-u/μの左辺がゼロの定常状態に達した場合ならu=μXですが,これが易動度という言葉の意味ですね。

流体力学のストークス(Stokes)の抵抗法則によれば,ブラウン粒子を半径aの球と仮定しηを媒質の粘性率とすると,この程度のレイノルズ数では 1/μ=6πηaと書けるはずです。

ストークスの法則の詳細は,2007年7/27の「遅い粘性流(1)(Stokes近似)」,および2007年7/28の記事「遅い粘性流(2)(Stokes近似)

そして,それに続く2007年7/31の記事「遅い粘性流(5)(Stokes近似)」を参照してください。

ブラウン運動の方程式:m(du/dt)=X(t)-u/μはランジュバン(Langevin)方程式:du/dt=-γu-R(t)/m (R(t)は"ゆらぎ=揺動",または雑音の影響を表わす量)と同じものです。

ランジュバン方程式については非平衡統計力学の線型応答理論に関連した2007年7/20の記事「揺動散逸定理 」を参照してください。

運動方程式:m(du/dt)=X(t)-u/μ (u=dx/dt)の両辺にx=x(t)を掛けてtについて0からtまで積分すると,左辺はm∫0t(xdu/dt)dt=m∫0t(xd2x/dt2)dt=[mxdx/dt]0t-m∫0t(dx/dt)2dtです。

 

一方,右辺は∫0tXxdt-(1/μ)∫0t(xdx/dt)dt=∫0tXxdt-(1/μ)∫0t[d/dt{x2/2}]dt=∫0tXxdt-(1/μ)[x2/2]0tとなります。

したがって,[mxu]0t-m∫0t2dt=∫0tXxdt-(1/μ)[x2/2]0tを得ます。

 

ところが,tが十分長ければ右辺第1項∫0tXxdtはX(t)が正負の全くランダムな値を取るために消えるはずです。

そして,統計力学のエネルギー等分配の法則により長時間平均の意味でm∫0t2dt=(t<mu2>)t→∞=kBTtとなります。ここにkBはボルツマン(Boltzmann)定数,Tは絶対温度です。

そこで,t→ ∞では[mxu]0t-kBTt=-(1/μ)[x2/2]0tです。それ故,大きいtで[mxu]0tが省略できる場合には<x2AV≡([x2]0t)t→∞=2μkBTtとなります。

 

こちらの<x2AVは,前の<mu2>のような長時間平均ではなく時刻tにおける相空間平均(確率平均)です。

実際,<u21/2~ (kBT/m)1/2ですがエルゴード性により時刻tにおける空間平均という意味でもu~ (kBT/m)1/2すから,x ~(2μkBTt)1/2であれば,<mxu>AV=[mxu]0t~kBT(2mμt)1/2となります。

  

そこで,[mxu]0tの値は大きいtに対してt1/2のオーダーですからtと比較して省略できることがわかります。

一方,ブラウン運動を時間τごとに微小長さを進む酔歩の問題と考えると次のように考察されます。

これは私のブログでは既にずいぶん前に考察済みです。

 

2006年9/14のブログ記事「酔歩(ランダム・ウォーク) 」から該当部分を多少修正して再掲します。

(再掲開始)

 

1次元ではx軸の上で左右どちらにも1歩ずつ動くことができて1歩の長さが一定値λであるとします。そして左右どちらかに進む確率は両側共に1/2であるとします。

x軸の原点から出発しててN歩の後にx=nλ(-N≦n≦N)の位置にいる確率をP(n,N)とすると,正の向きにN+,負の向きにN-歩いてxに到達するとした場合の数がN!/(N+!N-!)ですから,P(n,N)={N!/(N+!N-!)}(1/2)Nとなるはずです。

ただし,N++N-=N,+-N-=nなので単純に計算すると+(N+n)/2,-(N-n)/2ですが,N+nとN-nは一方が奇数なら他方も奇数,一方が偶数なら他方も偶数です。これらが偶数でなければ負でない整数なることが必要な+もN-も存在しません。 

したがって,N-nが偶数のときには(n,N)={N!/(N+!N-!)}(1/2)N(N+(N+n)/2,-(N-n)/2)ですが,N-nが奇数のときには(n,N)=0 です。

ここで,nが非常に大きいときのスタ-リングの公式:n!~(2π)1/2exp(-n)n(n+1/2),あるいはlog(n!)~(1/2)log(2π)+(n+1/2)log(n)-nを使用します。

N-nが偶数でNが非常に大きいとすれば,+(N+n)/2,-(N-n)/2も非常に大きいことになってlog{(n,N)}~ -Nlog2+NlogN-N+logN+-N-logN-(1/2)log(2π)+(1/2)(logN-logN+logN-)=(1/2)log{2/(πN)}-(N/2)[{1+(n+1)/N}log{1+(n/N)}+{1+(1-n)/N}log{1-(n/N)}]です。

ここで,n<<Nと考えて上のα=n/Nの対数関数において,αの2次までの近似展開:log(1-α)~ -α-α2/2,log(1+α)~ α-α2/2を利用します。

すると,log{(n,N)}~(1/2)log{2/πN)}-(N/2)(n/N)2より(n,N) ~{2/(πN)}1/2exp{-n2/(2N)}です。

ここで,x=nλ(-N≦n≦N)の酔歩の1歩の長さλは非常に小さいとしてN-nが偶数と奇数の両方の場合を考慮すれば,xがxとx+dxの間にある確率は(x,N)dx=(1/2){2/(πN)}1/2exp{-x2/(2Nλ2)}(dx/λ)(2πNλ2)-1/2exp{-x2/(2Nλ2)}dxです。

これは,Nλ→ ∞,λ→ 0,かつNλ2→σ2(有限)の条件で,xについて積分すると1になるので,確かに確率密度の条件を満たしています。

 一方,2次元での確率密度は,モデルが等方的なので単純に上の1次元の式で2をr2≡x2+y2で置き換えるだけでいいと考えられるところですが,実は1歩の各方向への成分Δx,ΔyがΔx2+Δy2=λ2を満たす必要があるので修正が必要です。

x方向とy方向を対等に扱うと,Δx2=Δy2=λ2/2なのでN歩で位置=(x,y)に到達する確率密度は,全平面で1になるように規格化してP(x,y,N)dxdy=(πNλ2)-1exp{-x2/(Nλ2)}exp{-y2/(Nλ2)}dxdy=(πNλ2)-1exp{-r2/(Nλ2)}2になります。 

同様に,3次元ではr2=x2+y2+z2としてΔx2=Δy2=Δz2=λ2/3により,位置=(x,y,z)に到達する確率は(x,y,z,N)dxdydz={(2/3)πNλ2}-3/2exp[-r2/{(2/3)Nλ2}]3になると考えられます。 

  特に,3次元の一般式でt=Nτ,D=λ2/(6τ)と置けば,4Dt=(2/3)Nλ2となるため,時刻tに位置に存在する確率はP(,t)=(x,y,z,N)=(4πDt)-3/2exp{-r2/(4Dt)}と書けます。

これは,丁度拡散係数がDの拡散方程式∂P/∂t=D∇2Pにおいて初期時刻t=0 に確率密度が原点に集中しているときの解,すなわち,初期値がP(,0)=δ3()のときのP(,t)の一意解に一致しています。(再掲終わり)※

そこで,この確率密度(,t)に基づいて計算すれば,x方向の"ゆらぎ=揺動",つまりt=0 に確率1で原点=0 にあった場合の時刻tでの平均位置(=0)からのずれxの2乗平均値は,<x2AV=∫-∞[x2(,t)]d3(4πDt)-1/2-∞[2exp{-x2/(4Dt)}dx=2Dtとなります。

これを,先にブラウン運動の1次元方程式から求めたxのゆらぎの表現:<x2AV([x2]0t)t→∞=2μkBTtと比較すれば,D=μBT(アインシュタインの関係式)が得られます。

,媒質の粘性率をηとし拡散粒子の半径をaとすれば,先に求めた易動度μに対するストークスの式:1/μ=6πηaから,D=kBT/(6πηa)なる等式が得られます。これをアインシュタイン・ストークスの関係式と呼びます。

ここで理科年表によると,空気の粘性率は温度が常温25℃=298Kでη=18.2×10-3Ns/m2,またkB=R/NA=1.38×10-23J/Kです。R~8.31J/(Kmol)は気体定数,NA~6.02×10-23/molはアヴォガドロ数(ロシュミット数)です。 

また,質量がmの気体分子の速度をとすると,統計力学によって平衡状態での2乗平均の速度は(<2AV)1/2=(3B/m)1/2です。一方,速度の絶対値||の平均は<||>AV={8B/(πm)}1/2です。

 

空気をO2とN2の1:4の混合気体とみると分子質量はm=28.8/NA(g)~ 4.78×10-26(kg)ですから,常温T=298Kでの2乗平均速度は21/2 ~ (3B/m)1/2~ 約508(m/s)です。

 

一方,絶対値平均速度で見ると,<||>AV={8B/(πm)}1/2~ 約468(m/s),ですです。

さて,気体分子を直径がdの剛体球とモデル化すると,2つの気体分子の中心間距離がdになるときにこれらは衝突します。

 

そして,1つの分子が衝突するまでに分子が移動する平均の距離を平均自由行程と呼びます。これは先に述べた酔歩の1歩に相当するので同じ記号λで表わすことにします。

 

1つの気体分子から見るとその中心を底面中心として体積がπd2λの円筒内に他の分子が1個入るという勘定になります。

 

したがって,気体分子数密度をnとするとπd2λn=1ですから平均自由行程はλ=1/(πd2n)と評価されます。

  

さて,速度=(u,v,w)=(vx,vy,vz)の各成分がvx~vx+dvx,vy~vy+dvy,vz~vz+dvzの間にある分子数をf(vx,vy,vz)dvxdvydvz=f()d3とします。

 

全分子数をNとするとf()は全速度空間で積分して∫f()d3=Nとなるように規格化されています。このように定義されるf()=f(vx,vy,vz)を速度の分布関数と呼びます。

 

そして,絶対温度がTの熱平衡状態では,f()がマクスウェル(Maxwell)分布:f()=N{m/(2πkBT)}3/2exp{-m2/(2kBT)}で与えられることがわかっています。

 

分子数密度がn=N/Vの熱平衡状態を仮定します。

 

通常のxyz空間のz=0 の面の単位面積を通って単位時間にz>0 の側からz<0 の側に移動する分子数をI+とすると,I+=n{m/(2πkBT)}3/2-∞dvx-∞dvy0dvzzexp{-m2/(2kBT)}=(n/4){8kBT/(πm)}1/2と計算されます。

 

ところが,先述したように<||>AV={8kBT/(πm)}1/2です。また,対称性からz=0 の単位面積を通ってz<0 の側からz>0 の側に移動する単位時間当たりの分子数をI-にとすると,これはI+は等しいのでI+=I-=(n/4)<||>AVと書けます。

 

そして,巨視的な粘性率を微視的な分子から統計的平均量として見積もるために,z=0 の面を通して下方から上方へと輸送されるx方向の運動量を評価してみます。

 

これは,z=0 の面から平均自由行程程度の下方の運動量が上方に運ばれ,これから下方に運ばれる平均自由行程程度の上方の運動量を引いた差で与えられると考えられます。

 

その平均自由行程程度のz座標を±αλ(0<α<1)とします。

 

まず,z=-αλから入ってくる運動量のx成分はx方向の速度成分をzだけの関数としてu(z)と表わせばI+mu(-αλ)=(mn/4)<||>AV{u(0)-αλ(∂u/∂z)}と見積もられます。

 

同様にz=0 の面を通って上方から下方へと輸送されるx方向の運動量はI-mu(αλ)=(mn/4)<||>AV{u(0)+αλ(∂u/∂z)}と見積もられます。

 

結局,z=0 の面を通って下方から上方へと輸送される正味の運動量のx成分は,I+mu(-αλ)-I-mu(αλ)=-(αλρ/2)<||>AV(∂u/∂z)となります。ただし,ρ≡mnは媒質の気体の密度です。

 

得られた単位時間当たりの輸送量-(αλρ/2)<||>AV(∂u/∂z)が,流体力学における現象論的粘性応力:-η(∂u/∂z)に一致すると考えられるので,粘性率ηに対してη=(αλmn/2){8kBT/(πm)}1/2なる評価式が得られました。

 

 これにλ=1/(πd2n)を代入するとη=(α/d2)(2πmkBT)1/2となります。それ故,D=kBT/(6πηa)={kBTd2/(6πaα)}(2πmkBT)-1/2=(2π)-3/2(kBTd4/m)1/2/(3aα)が得られます。

  

この式によれば,αがわかれば半径aが既知のブラウン粒子の空気中での分子拡散係数Dを測定すれば,空気分子の平均直径dを計算により評価できることがわかります。

 

実験等から得られる空気分子の径の評価値は0.4~0.8μmで可視光線の波長と同程度だそうです。

  

今日はここで終わります。

 

次回は電離層などプラズマの影響,大気層における空気分子や雲(水滴)によるレイリー散乱,ミイ散乱なども考慮して,地球面頂上でのフレネル反射や吸収の寄与による"太陽輻射(太陽定数)の減衰=アルベド(albed)"の定量的評価を論じることを予定しています。

 

参考文献:中村 伝 著「統計力学」(岩波書店),北原和夫 著「非平衡統計力学」(岩波書店),クドリャフツェフ 著(豊田博慈 訳)「熱と分子の物理学」(東京図書)

 

 ← クリックして投票してください。(1クリック=1投票です。1人1日1投票しかできません。クリックすると「人気blogランキング」に跳びます。)

にほんブログ村 科学ブログへ にほんブログ村 科学ブログ 物理学へ クリックして投票してください。(ブログ村科学ブログランキング投票です。1クリック=1投票です。1人1日1投票しかできません。クリックするとブログ村の人気ランキング一覧のホ-ムページームに跳びます。)

http://www.mediator.co.jp/category/pages.php?id=115「中古パソコン!メディエーター巣鴨店」

iconヤーマン プラチナゲルマローラー 1日3分コロコロエステ!ローラー型プラチナ配合美顔器  

ブックオフオンライン 

お売りください。ブックオフオンラインのインターネット買取 展開へ! ▼コミック 尾田栄一郎 「ONE PIECE(52)」 icon ▼コミック 「ONE PIECE」をオトナ買い icon

三国志特集 ▼コミック 横山光輝 「三国志全巻セット」 icon 「三国志(文庫版)全巻セット」 icon  「三国志(ワイド版)全巻セット」 icon  ▼書籍 「三国志」/吉川英治 icon  「三国志」/北方謙三 icon  「三国志」/宮城谷昌光

iconオンライン書店 boople.com(ブープル)

| | コメント (0) | トラックバック (0)

より以前の記事一覧

その他のカテゴリー

001. 目次 002. 募金・ボランティア 003. 日記・回想 004 訃報 005. 心身・思想・哲学 006. 社会・経済・政治 007. 病気(診察・薬) 008. 恋愛・異性 009 宗教・神話 010 歴史(日本,世界) 011. 将棋 012. TV(ニュース・ドラマ) 013 スポーツ(ニュース・イベント) 014 ノン・フィクション 015 小説・詩・評論 016 漫画・劇画・アニメ 017 演劇・映画・舞踊 018 音楽(日本・西洋・他) 019 タレント(俳優・お笑い) 020 ミュージシャン 021 アイドル・ヒーロー 022 創作 023 シャレ・ギャグ等 024 競馬・toto・賭け事 025 ファッション・風俗 100. 物理学一般 101 教育・学校(物理) 102. 力学・解析力学 103. 電磁気学・光学 104. 熱力学・統計力学 105. 相対性理論 106. 星・ブラックホール・一般相対性 107. 重力・宇宙・一般相対性 108. 連続体・流体力学 109. 物性物理 110. 複雑系・確率過程・非線型・非平衡 111. 量子論 112. 原子・分子物理 113. 原子核物理 114 . 場理論・QED 115. 素粒子論 116. 弦理論 118. 観測問題・量子もつれ 119. 電気回路 200. 問題・解答 201. 自然科学一般 202. 気象・地学・環境 203. 生物学・生理学・生化学 204. 経済学(ミクロ・マクロ・マルクス) 300 数学一般・算数 301. 集合・位相 302. 論理学・数学基礎論 303. 代数学・数論 304. 解析学 305. 複素数・複素関数論 306. 線型代数学 307. 幾何学(トポロジー・他) 308. 微分方程式 309. 確率・統計 310. 関数解析・超関数 311 .数値計算・調和解析・離散数学 312. 公式・特殊関数 501. 商用宣伝・アフィリエイト