1
2.ミクロカノニカル分布とエントロピー 2-1理想気体
分子間力が全く働かない気体(実在はしない)
2-1.1 量子力学の簡単な考え方
古典 量子力学
位置ベクトル r 波動関数
r
確率密度
r 21次元粒子を考える
x ae
ikx
(2.1)k
p
(2.2)有限領域に閉じ込められた粒子
壁の効果 → 壁の外では波動関数が0に なる
壁の位置を
0 , L
とすると、(教科書は弦の場合、
右図は1次元量子量子井戸)
0 L 0
(2.3) 境界条件弦の振動の場合、固定端の条件下で起こる基準振動は定常波。 弦の変位
u x
は x u kx
u
osin
(2.4) 1 , 2 , 3 ,
n
L
k n
(2.5) 図2-1 参照波動関数も同様にして、
x a sin kx
(2.6) または e
ikxe
ikx
i
x a
2
(2.7)(2.7) の第1項、2項は各々 x 軸正方向、負方向に進む波を表す。 この重ね合わせ。
この時、粒子の取り得る運動量は
1 , 2 , 3 ,
n
L k n
p
(2.8)
と、飛び飛びの数となる。(離散化、量子化)
E
0 L x
2 2-1.2 周期的境界条件
結晶の様に周期的に同じ原子配列を取っている場合、ある距離 L を進むと元と同じ状態が実現
x x L
(2.9)これに波動関数の定義 (2.1) を考えると
x L
ik
ikx
e
e
∴e
ikL 1
0 , 1 , 2 , 3 ,
2
n n
k L
(2.10)よって運動量は量子化され
0 , 1 , 2 , 3 ,
2
n n
p L
(2.11)即ち、固体の中では運動量は
p 2 L
を単位として量子化され、離散的。粒子のエネルギー
は、m を粒子の質量としてm p 2
2
(2.12)である。 従って、エネルギーも量子化され、
2 2 2 22 2 2
2
1 n
n mL L m
(2.13)エネルギーが
より小さい量子状態の数を
とすれば、これは運動量がm
p 2
の量子状態の数に等しい。 これは、
(a) 境界条件 (2.9) では
2 m p 2 m
に間隔L h L
2
で分布(b) 境界条件 (2.3) では
0 p 2 m
に間隔L h L
2
で分布物理的には、図2-1の両端を固定した弦振動で n が奇数のものは周期的境界条件を適用すると、
連続になっていない事に注意。
0 0
|p|
p h/L
h/2L
3 いずれの場合も
L m
2
(2.14)2-1.3 3次元粒子の波動関数とエネルギー 1次元の場合と同様に周期的境界条件を考える
x , y , z x L , y , z x , y L , z x , y , z L
(2.15)すると、波動関数として、平面波を考える事ができ
x , y , z ae
i
kxxkyykzz
(2.16)が量子状態として存在。 波数ベクトル
k k
x, k
y, k
z
は境界条件 (2.15) より 1
ikL ikLL
ikx
e
ye
ze
∴
2 0 , 1 , 2 ,
2 ,
2 ,
x y y z z ix
n n
k L L n
k L n
k
(2.17)∴
, 2 , 1 , 2 0
2 ,
2 ,
x y y z z ix
n n
p L L n
p L n
p
(2.18)即ち、3次元の運動量空間
p
x, p
y, p
z
を考えると体積2
3
L
に1個の割合で存在する。
注:教科書の図2-3 は2次元の図だが、議論は3次元で行っている事に留意。
3次元の粒子のエネルギー
と、運動量p
の関係はm 2
p
2
, p2 p
x2 p
y2 p
z2 (2.19) 00
|p|
p h/L
h/2L
x x L
0 L 0
4 よって、エネルギーは量子化され、
22
2 2 2
2 2
z y
x
n n
mL n
(2.20)エネルギーが
より小さい量子状態の数を
は運動量空間における半径2 m
の球の中にある状態数で、体積
V L
3 の容器内を考えると
3 2
323 4
2
V m
(2.21)☆ 量子状態の数を数える事は、即ち、図2-3における格子点の数を数える事に帰着。
(参考)同じ事は不確定性関係
x p
x~ h 2
からも導ける事に留意。(補足)位相空間における離散的状態
物質の運動ないし状態は、一般にその物質の位置
q x , y , z
と運動量 p
px, py, pz
という6つのパラメーターで表す事が出来る。 とすると、粒子
N
個のアンサンブルでは、6 N
という自 由度を持った空間で表現が可能である。 が、量子力学の不確定性関係より、 x p
z~ h
であるので(h は 既出の Planck 定数:
6 . 626 10
34 J s
)、h
3N が位相空間の単位体積である。この事をもう少し言い方を変え、再度説明する。
図2A-1 3次元速度空間と1方向における許容される波。
0 L
0 L
0 L
0 L
n = 1/2
n = 1
n = 3/2
n = 2 vz
vx
vy
5 ここで、原子・分子や電子と云ったミクロな粒子を扱うには、「量子力学」の適用が必要である。 量
子力学に従うと、運動量は
mv h
p
(
は de Broglie(物質)波長) (2A.1) と書き表される。一辺
L
の立方体に閉じ込められた気体は、前に取り扱った様に、壁に衝突したり、また、分子間 の衝突により、個々の原子(分子)は速度空間の中での「点」を飛び移る。従って
mv
h
また、 m
v h
(2A.2)問題にしている空間を定めた(ここでは一辺
L
の立方体)時、L
が
の整数倍になっている事が 量子力学の要請である(教科書の (2.9) 式、図2-2(a) に相当)。 すると、mv n h n
L
∴n mL mL n h
v 2
(2A.3)この結果、上のピンクの点が実は連続的ではなく、とびとびの値しか取らなくなるのである。
z z
y y
x
x
n
v mL mL n
v mL n
v
2
2 ,
2 ,
(2A.4)この結果、速度は量子化されて、
mL
2
の整数倍の値を取る事になる。 図を分かり易くするため、y
x
の2次元速度空間を考えると、 v
x, v
y
は下図の様に、碁盤の目の値を取る事になる。注: 連続的な空間の体積でなく、「碁盤の目=格子点」を考える事が量子化の意味する処。 しかし、
碁盤の目が細かくなれば、両者は(単に比例定数を掛けるというだけで)次第に同じ結果に落ち着く。
例えば、水素原子を考え、一辺
10 cm
の空間に存在しているとすると、] [ 10 66 . 10 1 02 . 6
] [
10
2723 3
kg kg
m
(2A.5)速度は x(無論、
y
or z でも構わないが) 方向成分で考えるとT k v
m
x B2 1 2
1
2
∴m T
v
x2 k
B (2A.6)室温
T 300 K
であれば、∴
1 . 58 10 [ / ]
] [ 10 66 . 1
] [ 300 ] [
10 38 .
1
327 1 23
s kg m
K K
J m
T
v
xk
B
(2A.7)のオーダーである。 一方、考える空間
L n
であり、L 10
1 m
なので、6
] / [ 10 99 . ] 3 [ 10 ] [ 10 66 . 1
] [ 10 626 .
6
61 27
34
s m m
kg s J mL
v
ih
(2A.8)のオーダーで量子化されている。
図2A-2 2次元速度空間の量子化。 原子・分子の速度は格子点を占める。
さて、各点における運動エネルギーは、
2 2 2
2 1
z y x
k
mv mv mv
E
(2A.9)と書き表され、原点からの距離の自乗に運動エネルギーは比例する。
(因みに、
p
空間でも、v 空間でも、係数 m を考慮すれば同じ事になる)従って、あるエネルギー範囲
E
i E
j E
k に入る状態の数は同心球(3次元の場合。 図2-3A は2次元なので円)m v E v
v
x2
y2
z2 2
i とm v E v
v
x2
y2
z2 2
k の間の格子点の数 (2A.10)すると、運動量 p は
p mv 4 . 54 10
24[ kg m / s ]
、この時の de Broglie 波長は] [ 10 46 . ] 1 / [ 10 54 . 4
] [ 10 626 .
6
1024 34
m m s J
s J p
h
(2A.11)である。
vy
vx
h/mL h/mL
7
図2A-3 2次元速度空間の格子点と等エネルギー面
ところで、冒頭に述べた通り、量子力学においては、不確定性関係が
x p
x~ h 2
と、実空間位置(座標) vs. 運動量 で定義されているので、運動量空間で考えるのが通常である。
従って、以下、p
px, py, pz
を用いる。 すると、mE p
p
p
x2
y2
z2 2
(2A.12) である。図2A-4 2次元運動量空間の量子化と等エネルギー面
py
px
h/L h/L
vy
vx
h/mL h/mL
8 ここに十分碁盤の大きい時、面積・体積を求める事が碁盤の目の数を数える事に相当。
2次元:
2 mE p
x2 p
2y p
2 2 m ( E dE )
にある点の数を考えると、dE h m
dp L h p
L 2 2
2 2
(2A.13)3次元:
2 mE p
2x p
2y p
2z p
2 2 m ( E dE )
dE E h m
dp L h p
L
33 2
3
2 4
4
(2A.14)ところで、
L
は始めに適当に定めた空間だが、3次元でL
3 V
(体積)であり、結局、状態 の数は定めた体積に比例する事になる。 要は、マクロな空間 → 離散的運動量(速度)が連続的に(厳密に連続ではない)
という事で、運動量空間での状態の数を計算する事にある。
運動量の離散の具合は、
10 cm
空間では] / [ 10 626 . ] 6
[ 10
] [ 10 626 .
6
331 34
s m m kg
s J L
p
i h
(2A.15)という微小量となる。
(補足 終)
9
2-2 理想気体のエントロピー
多粒子系の量子状態とエントロピー
独立に運動する多粒子系で、構成粒子に番号をつけ、各々の運動量を定めると、
p
1, p
2, p
3, p
N
がN
粒子系の量子状態を定める事になる または、各粒子の運動量ベクトル p
x, p
y, p
z
を粒子1:
p
1, p
2, p
3
粒子2: p
4, p
5, p
6
... とすると、N
3
次元ベクトル p
1, p
2, p
3, , p
3N
と表す事ができる 周期的境界条件により
0 , 1 , 2 ,
2
i ii
n n
p L
(2.22)と量子化され、
3 N
次元の運動量空間に、体積 NN
V L 2
3
に1個の割合で分布。
N
粒子系全体のエネルギーは
Ni
p
iE m
3
1 2
2
1
(2.23)エネルギーが
E
より小さい量子状態の数 E
は、2 mE
の3 N
次元球の体積を量子化した 時の格子点の数。 一般に、半径R
の n 次元球の体積は p.277 の公式 (A.7)
nn
n
R
n R n
V 2
2
2
(2.24) (課題) (2.24) 式を各自導出せよ。
但し、
0
1
e dt t
z
z t (A.4) はガンマ関数で、
! 2
! 2 2 1
2
n n
nn
(A.6)eg. 例えば、3次元球では
n 3
で、
32 3 32
3 32 2 3 32 3 33
3
4 2
! 4 3
2
! 2
! 1 2 3
2 3 2
1 3
2 2
3 3
2 R R R R R
R
V
よって、
N
粒子系のエネルギーがE
より小さい量子状態の数は E
は
3
3 2 2
3 22 3 3
2 2
N N N
N
N mE N
E V
(2.25)10 エネルギーが
E
とE E
の間にある量子状態の数W E
は
dE E E E d
W
(2.26)なので、
E mE E
N E V
W
NN N
N
3
3 22
3 22 2 3
(2.27)従って
V k N k E E
N Nk mE
E k E
k N V N
Nk mE
E k E
k N Nk N
mE V Nk
E k E
k N N N k N
mE Nk V
E k E
k N k N
mE Nk V
E k E
k N V mE
Nk
E mE E
N k V
E W k
E S
B B
B
B B
B
B B
B B
B B
B B
B B
B B
B B
B
N N N N B
B
2 log log 3 2
log 3 2
3 log 4 2 3
2 log log 3 2
3 2 log 3 2 log 3 2
log 2 2 3
2 log log 3 2
3 2 log 3 2 log 3
2 log 2
2 log log 3 2
3 2 log 3 2 3 2
log 2
2 log log 3 2 !
log 3 2
log 2
log
! 2 1 3 log 1 2
2 log
2 2 2 3
log log
2 2
32 2
3 2 3 3
2 3
32 3
2 3
3 2 2
3 3
(2.27)’
1章の議論と同様に、上の (2.27)’ 式の
2 log 3 2
3 log 4 2 3
2
V
N Nk
BmE
はNk
B のオーダーE k E
k
BN
Blog 2
log 3
はk
Blog N
のオーダーよって、第2項は無視でき、
2
log 3 2
3 log 4 2 3
2
V
N Nk mE
E
S
B
(2.28)さて、
N 2 N , E 2 E , V 2 V
として、エントロピーを再考すると11
NN
NmE
NN N
E V
33 3
2
3 4 6
2 2
2 2
(2.25)’
E mE E N
E V dE
E E d
W
NN N
N
6 3 32
3 4 2
2
(2.26)なので、
2 2 log 2
2 log 3 2 2 log
3 log 4 2 2 3
2 log log 3 2
3 3 2 log 2 3
2 log log 4 2 2 3
log 3
log 3
3 log 3 2
log 2
log 4 2
log 3
log 3
3 log 2 3
4 log 2 2
log 3
log
! 3 2 log
4 log 2 2
! log 1 3 log 1 4
2 log 2 2
3 4 2
log 2 log
2 2
32 2
3 2 3 3
2 3
2 3 3
2 3
3 3
6 2
B B
B B
B
B B
B B
B B
B B
B B
B B
B B
B
N N
N N B
B
Nk E
S N V
Nk mE
E k E
k N N
mE V Nk
E k E
N k
N Nk
mE V Nk
E k E
N k
N N N mE k
Nk V
E k E
N k
N mE k
Nk V
E k E
k N V mE
Nk
E mE E N k V
E W k
E S
即ち、
S 2 E S E S E
であり、示量変数となっていない!!何がマズイか? → 気体系では、粒子が識別できないものと考えるのが相当
すると状態数を計算し直すと、
E mE E
N N E V
W
NN N
N
3
3 22
3 22 3 2 !
(2.29)(分母にある
N !
がポイント)上と同様にエントロピーの計算をすると
12
mE V N N k N k E E
Nk
E k E
k N N N
Nk mE V
Nk
E k E
k N N N N N
N k N
mE Nk V
E k E
N N k
N k mE k
Nk V
E mE E
N N k V
E W k
E S
B B
B
B B
B B
B B
B B
B B
B B
B
N N N
N B
B
2 log log 3 2
log 5 2 log 3 2 log 3 2
log 2 2 3
2 log log 3 2
log 5 2 log 3 2 log 3
2 log 2
2 log log 3 2 log
3 2 log 3 2 3 2
log 2
log
! 2 log
log 3 2 !
log 3 2
log 2
2 2 3 2 !
log log
2 32 2
3 2 3 3
2 3
3 2 2
3 3
2
log 5 log 3
2 log 5 2
3 log 4 2
3
322
2
N
V N
Nk mE N
V N
Nk
BmE
B
(2.30)さて、
N 2 N , E 2 E , V 2 V
として、エントロピーを再考すると
E mE E N
N E V
dE E E d
W
NN N
N
6 3 32
3 4
! 2 2
2
(2.26)’
S E
N V N
Nk mE
E k E
N k
N N
Nk mE V
Nk
E k E
N k
N N N
N N N
mE k Nk V
E k E
N k
N k
N mE k
Nk V
E mE E N N k V
E W k
E S
B
B B
B B
B B
B B
B B
B B
B
N N
N N B
B
2 2 log 5 2
3 log 4 2 2 3
log 3
2 log 2 5
log 3
2 log 2 3
2 log 2
log 4 2
log 3
log 2
2 log 2 3 3 log 3 2
4 log 2 2
log
! 2 log 3
log
! 3 log 2
4 log 2 2
3 4
! 2 2
log 2 log
2 32 2
3 2 3 3
2 3
3 3
6 2
よって、エントロピーがきちんと示量変数となっている。
但し、この議論も1点ごまかしがある。 2個以上の粒子が同じ運動量を持つ状態は重複して数 えていない。 しかし、一つの量子状態を占める確率が十分小さければ2個以上の粒子が同じ運動量 を持つ確率は(その自乗)更に小さいので、無視する。 → 次節で議論(粗視化)。
13 エントロピーの熱力学的定義 (1.44) 式
T dE dS 1
より
Nk E dE
E dS
T
B1 2
3
1
よって、理想気体のエネルギーとして、
T Nk
E
B2
3
(2.31)気体分子分配のエントロピー
前に議論した事を思い起こし、2つの領域に分けられた気体のエントロピーを考える。(図2-4)
領域1,2に エネルギーが
E
1, E
2
、分子数が N
1, N
2
に配分された部分平衡状態の エントロピーは E
1, E
2; N
1, N
2 S
1 E
1, N
1 S
2 E
2, N
2
S
(2.32)ここに、
E E
1 E
2 const
,N N
1 N
2 const
である。よって、(1.37), (1,38) を求めた時と同様、エントロピー最大となるのは
2 2 2 2 1
1 1 1 1
2 2 2 1
1 1 1 1
2 1 2
1
, ; , , , , ,
E N E S E
N E S E
N E S E
N E S E
N N E E S
(1.37)’これが 0 となるのは、
2 2 2 2 1
1 1
1
, ,
E N E S E
N E S
(1.38)’(1.23) と同様にして、
N E N E N
E
2 2 1
1 (2.33)
2
log 5 2
3 log 4 2 3 2
log 5 2
3 log 4 2 , 3
1 1 2
1 2
1 1 1
1
1
N
V N
k mE N N
V N
k mE N N E
S
B B
V 1 V 2
E
1, N
1E
2, N
2図2-4
14
2 2
1 2 1
2 2 2
1 1 2
2 2 2 1 1 1 2 1 2 1
log 2 log
5 2
3 log 4 2 3
2 log 5
2 3 log 4 2 3 2
log 5 2
3 log 4 2 3
, ,
,
; ,
N N V
N N V k N
Nk mE
N V N
k mE N N
V N
k mE N
N E S N E S N N E E S
B B
B B
(2.34)この内、分子数配分に依存する部分は
N V N
N N
N N
N N
Nk N N
N V N N V
k
Blog log
Blog log log
2 2
1 1
2 2
1 1
これは1-1節で2項定理を用いて求めた
N
n N N
n N N
n N N n
n
P
Nlog 2 log log
log
(1.5)と定数を除き一致している。
2-3 理想気体の速度分布則
理想気体分子1個のエネルギーは、(2.31) より、
k T N
E 2
B 3
(第1回の演習も参照)しかし、実際には個々の粒子のエネルギーは分布を持っている。 さて、どの様に?
粗視化された分子分布
ある瞬間、理想気体の個別の状態を見たとする → 分子が様々な量子状態に分布
→ 分子の状態を粗視化
1)1粒子状態をエネルギーによりグループ分けし、グループにエネルギーの低い方から
1, 2, 3, … の番号をつける(図2-5)
(1) エネルギーの幅は十分小さく、1つのグループに属する量子状態のエネルギーは 中心の値で代表させる
(2) 各グループに属する量子状態の数は十分に多い
グループ l に属する量子状態の数を
M
l, エネルギーをE
l, そこを占める分子の数をN
l と する。 N
1, N
2, N
3,
が粗視化された分子分布。 全粒子数N
は、
l
N
lN
(2.35)全エネルギー
E
は、
l l l
E N
E
(2.36)15 分子分布の関数としてのエントロピー
特定の分子分布
N
1, N
2, N
3,
の下での全系の量子状態の数を求める事を考えるエネルギー値
E
l を有する(で代表される)量子状態数をM
l として、ここにN
l 個の分子を分 配する事を考える。(上図は図2-5 と等価)この様な分配の場合の数は
M
lNl であり、分子の数が状態数に比べて十分少ない 1
l l
M
N
(2.37)とすると、各々の量子状態には分子は0個若しくは1個しか入っていないと考える事ができる。
これらの分子が区別できないとすると、場合の数は
l
!
N l
N M
l通り
全てのグループについて同じ議論が出来、分布
N
1, N
2, N
3,
の下での全系の量子状態の数は
l l
N l
N N M
N N W
l
, ! , ,
2 31 (2.38)
Stirling の公式を用い、
l
l l l l l
B
B
W N N N k N M N N N
k N
N N
S
1,
2,
3,
log
1,
2,
3,
log log
l l
l l
B l
l l l l
B
N
N M k N N
N M
k log log 1
(2.39)(参考)少し異なる考え方でも上記は導ける
l :粒子数
l :量子準位数
l
16 分子が区別できるとすると
1 1 1
!
!
l l l l ll l
l N
M
M M M M N
N M P M
l
l
通り十分状態数が大きく、占有数が少ないとする。 即ち、
1
l l
M
N
(2.37)とすると、 Ml Nl
M
l lN
l M
lNlP M
!
!
通りここで、分子が区別出来ない場合には、分配の仕方は、
l
!
N l
N
M
l
MlC
Nl
通り (参考終)Lagrange の未定係数法 (2変数の場合) p.49-50 参照
2変数 x, y の関数
x y f
u ,
(2.40)を考える。 この関数の極値を
x y C const
g ,
(2.41)という条件の下で求めるには、未定係数 を導入して関数
x y g x y f
u ~ , ,
(2.42)について(条件無しで)極値を求めて、係数 は極値を与える
x
0, y
0
が条件 (2.41) を満たす 様に定めれば良い。(証明)(2.42) が極値を取る条件、即ち
, 0
~ ,
, 0
~ ,
y y x g y
y x f y u
x y x g x
y x f x u
(2.63)
が (2.41) を満たし、かつ、(2.40) が極値を取る条件を満たしていればよい。
l
l
17 ところで、(2.41) より y は x の関数。 よって、 (2.40) が極値を取る条件は
, , , 0
dx
dy y
y x f x
y x x f
y x dx f
d dx
du
(2.61)一方、(2.41) を x で微分すると
, , , 0
dx dy y
y x g x
y x g dx
y x dg
これを (2.61) に代入して、
, , , , 0
x y x g y
y x f y
y x g x
y x
f
(2.62)(2.63) を (2.62) 左辺に代入すると
, , , , 0
x y x g y
y x g y
y x g x
y x
g
よって (2.62) は成立。
即ち
x y g x y f
u ~ , ,
(2.42)が極値を取る条件は、
g x , y C const
の条件下でu f x , y
の極値を求める事となる。(証明終わり)
熱平衡の分子分布
熱平衡状態の分子状態は、(2.39) のエントロピーを、式 (2.35), (2.36) の条件の下で最大にする 事で求まる。 → Lagrange の未定係数法を用いる。
☆ 極値を求めるべき物理量
l l
l l
B
N
N M k N
N N
S
1,
2,
3,
l o g 1
(2.39)★ 束縛条件
l
N
lN
(2.35)
l l l
E N
E
(2.36)この時、条件は二つ(粒子数、エネルギー)あるので、2つの未定係数 a, b を導入し
l l l l
l
l l
l B
l
a N b N E
N k M
N N
N N
S ~ , , , log 1
3 2
1 (2.43)
を最大にする事を考える事に帰着。