MB分布 | BE分布 | FD分布 | |
---|---|---|---|
微观状态数 | Ω M B = N ! ∏ l a l ! ∏ l g l a l \Omega_{MB}=\frac{N!}{\prod_l a_l!}\prod_lg_l^{a_l} ΩMB=∏lal!N!∏lglal | Ω B E = ∏ l ( g l + a l − 1 ) ! ( g l − 1 ) ! a l ! \Omega_{BE}=\prod_l\frac{(g_l+a_l-1)!}{(g_l-1)!a_l!} ΩBE=∏l(gl−1)!al!(gl+al−1)! | Ω F D = ∏ l g l ! ( g l − a l ) ! a l ! \Omega_{FD}=\prod_l\frac{g_l!}{(g_l-a_l)!a_l!} ΩFD=∏l(gl−al)!al!gl! |
分布 | a l = g l e α + β ϵ l a_l=\frac{g_l}{e^{\alpha+\beta\epsilon_l}} al=eα+βϵlgl |
a
l
=
g
l
e
α
+
β
ϵ
l
−
1
a_l=\frac{g_l}{e^{\alpha+\beta\epsilon_l }-1}
al=eα+βϵl−1gl = g l e ( ϵ l − μ ) / ( k B T ) − 1 =\frac{g_l}{e^{(\epsilon_l-\mu)/(k_BT) }-1} =e(ϵl−μ)/(kBT)−1gl | a l = g l e α + β ϵ l + 1 a_l=\frac{g_l}{e^{\alpha+\beta\epsilon_l }+1} al=eα+βϵl+1gl |
定域系统与满足经典极限条件的Bose及Fermi系统 | |||
特性函数 |
Z
=
∑
l
g
l
e
−
β
ϵ
l
=
∑
s
e
−
β
ϵ
s
Z=\sum_l {g_le^{-\beta\epsilon_l}}=\sum_s {e^{-\beta\epsilon_s}}
Z=∑lgle−βϵl=∑se−βϵs Z = V / λ T 3 Z=V/\lambda_T^3 Z=V/λT3 |
Ξ
:
=
∏
l
Ξ
l
=
∏
l
(
1
−
e
−
α
−
β
ϵ
l
)
−
g
l
\Xi:=\prod_l {\Xi_l}=\prod_l {(1-e^{-\alpha-\beta\epsilon_l})^{-g_l}}
Ξ:=∏lΞl=∏l(1−e−α−βϵl)−gl 常用 l n Ξ = − ∑ l g l l n ( 1 − e − α − β ϵ ) d ϵ ln\Xi=-\sum_lg_lln(1-e^{-\alpha-\beta\epsilon})d\epsilon lnΞ=−∑lglln(1−e−α−βϵ)dϵ l n Ξ = − ∫ D ( ϵ ) l n ( 1 − e − α − β ϵ ) d ϵ ln\Xi=-\int D(\epsilon)ln(1-e^{-\alpha-\beta\epsilon})d\epsilon lnΞ=−∫D(ϵ)ln(1−e−α−βϵ)dϵ |
Ξ
:
=
∏
l
Ξ
l
=
∏
l
(
1
+
e
−
α
−
β
ϵ
l
)
g
l
\Xi:=\prod_l {\Xi_l}=\prod_l {(1+e^{-\alpha-\beta\epsilon_l})^{g_l}}
Ξ:=∏lΞl=∏l(1+e−α−βϵl)gl 常用 l n Ξ = ∑ l g l l n ( 1 + e − α − β ϵ ) d ϵ ln\Xi=\sum_lg_lln(1+e^{-\alpha-\beta\epsilon})d\epsilon lnΞ=∑lglln(1+e−α−βϵ)dϵ l n Ξ = ∫ D ( ϵ ) l n ( 1 + e − α − β ϵ ) d ϵ ln\Xi=\int D(\epsilon)ln(1+e^{-\alpha-\beta\epsilon})d\epsilon lnΞ=∫D(ϵ)ln(1+e−α−βϵ)dϵ |
总粒子数 N N N | e − α Z e^{-\alpha}Z e−αZ | − ∂ ∂ α l n Ξ -\frac{\partial}{\partial \alpha}ln\Xi −∂α∂lnΞ | |
内能 U U U | − N ∂ ∂ β l n Z -N \frac {\partial }{\partial \beta}lnZ −N∂β∂lnZ | − ∂ ∂ β l n Ξ -\frac{\partial}{\partial \beta}ln\Xi −∂β∂lnΞ | |
广义力 Y λ Y_\lambda Yλ | − N ∂ β ∂ y λ l n Z , p = N ∂ β ∂ V l n Z -\frac{N\partial }{\beta \partial y_\lambda}lnZ,p=\frac{N\partial }{\beta \partial V}lnZ −β∂yλN∂lnZ,p=β∂VN∂lnZ | − 1 β ∂ ∂ y λ l n Ξ , p = 1 β ∂ ∂ V l n Ξ -\frac{1}{\beta}\frac{\partial }{\partial y_\lambda}ln\Xi,p=\frac{1}{\beta}\frac{\partial }{\partial V}ln\Xi −β1∂yλ∂lnΞ,p=β1∂V∂lnΞ | |
焓 H = U + p V H=U+pV H=U+pV | N ( V ∂ β ∂ V l n Z − ∂ ∂ β l n Z ) N(\frac{V \partial }{\beta\partial V}lnZ-\frac{\partial }{\partial \beta}lnZ) N(β∂VV∂lnZ−∂β∂lnZ) | k B ( l n Ξ − α ∂ ∂ α l n Ξ − β ∂ ∂ β l n Ξ ) k_B(ln\Xi-\alpha\frac{\partial }{\partial \alpha}ln\Xi-\beta\frac{\partial }{\partial \beta}ln\Xi) kB(lnΞ−α∂α∂lnΞ−β∂β∂lnΞ) | |
熵 S , d S = d Q T S,dS=\frac{dQ}{T} S,dS=TdQ | S = N k B ( l n Z − β ∂ ∂ β l n Z ) S=Nk_B(lnZ-\beta\frac{\partial }{\partial \beta}lnZ ) S=NkB(lnZ−β∂β∂lnZ) | ||
Helmholtz自由能 F = U − T S F=U-TS F=U−TS | − N β l n Z = − N k B T l n Z -\frac{N}{\beta}lnZ=-Nk_BTlnZ −βNlnZ=−NkBTlnZ | − 1 β ( l n Ξ − α ∂ ∂ α ) l n Ξ -\frac{1}{\beta}(ln\Xi-\alpha\frac{\partial }{\partial \alpha})ln\Xi −β1(lnΞ−α∂α∂)lnΞ | |
Gibbs自由能 G = F + p V G=F+pV G=F+pV G = N μ G=N\mu G=Nμ | − N β ( l n Z − V ∂ ∂ V l n Z ) -\frac{N}{\beta}(lnZ-V\frac{\partial }{\partial V}lnZ) −βN(lnZ−V∂V∂lnZ) | − 1 β ( l n Ξ − α ∂ ∂ α l n Ξ − V ∂ ∂ V l n Ξ ) -\frac{1}{\beta}(ln\Xi-\alpha\frac{\partial }{\partial \alpha}ln\Xi-V\frac{\partial }{\partial V}ln\Xi) −β1(lnΞ−α∂α∂lnΞ−V∂V∂lnΞ) | |
巨势 Ψ = F − G \Psi=F-G Ψ=F−G | Ψ = − N V β ∂ ∂ V l n Z \Psi=-\frac{NV}{\beta}\frac{\partial }{\partial V}lnZ Ψ=−βNV∂V∂lnZ | − 1 β l n Ξ = − V ∂ β ∂ V l n Ξ -\frac{1}{\beta}ln\Xi=-\frac{V\partial}{\beta\partial V}ln\Xi −β1lnΞ=−β∂VV∂lnΞ | |
Boltzmann公式 |
S
=
k
B
l
n
Ω
M
B
S=k_Bln\Omega_{MB}
S=kBlnΩMB 满足经典极限条件的Bose及Fermi系统 S = k B l n Ω M B N ! S=k_Bln\frac{\Omega_{MB}}{N!} S=kBlnN!ΩMB = N k B ( l n Z − β ∂ ∂ β l n Z ) − K B l n N ! =Nk_B(lnZ-\beta\frac{\partial }{\partial \beta}lnZ )-K_BlnN! =NkB(lnZ−β∂β∂lnZ)−KBlnN! | S = k B l n Ω B E S=k_Bln\Omega_{BE} S=kBlnΩBE | S = k B l n Ω F D S=k_Bln\Omega_{FD} S=kBlnΩFD |
配分函数计算 |
Z
=
∫
d
ω
(
2
π
ℏ
)
r
e
−
β
ϵ
Z=\int {\frac{d\omega}{(2\pi\hbar)^r}e^{-\beta\epsilon}}
Z=∫(2πℏ)rdωe−βϵ d ω = d q 1 ⋅ ⋅ ⋅ d q r d p 1 ⋅ ⋅ ⋅ d p r d\omega=dq_1···dq_rdp_1···dp_r dω=dq1⋅⋅⋅dqrdp1⋅⋅⋅dpr 或 Z = ∫ D ( ϵ ) e − β ϵ d ϵ Z=\int {D(\epsilon)e^{-\beta\epsilon}d\epsilon} Z=∫D(ϵ)e−βϵdϵ 配分函数Z为态密度的Laplace变换 |
单原子分子XX理想气体
非简并性条件:
e
α
≫
1
e^{\alpha}\gg1
eα≫1
对于单原子分子理想气体,用能量均分定理及理想气体状态方程即可解得其所有物理量。
**能量均分定理:**处于平衡态的温度为T的经典系统,其中非相对论性粒子的能量中每一个平方项的平均值均为
k
B
T
/
2
k_BT/2
kBT/2。
内能
U
=
N
ϵ
‾
=
3
2
N
k
B
T
=
3
2
n
R
T
U=N\overline{\epsilon}=\frac{3}{2}Nk_BT=\frac{3}{2}nRT
U=Nϵ=23NkBT=23nRT
焓
H
=
U
+
p
V
=
3
2
N
k
B
T
+
N
K
B
T
=
5
2
N
k
B
T
=
5
2
n
R
T
H=U+pV=\frac{3}{2}Nk_BT+NK_BT=\frac{5}{2}Nk_BT=\frac{5}{2}nRT
H=U+pV=23NkBT+NKBT=25NkBT=25nRT
定容热容
C
V
=
(
∂
U
∂
T
)
V
=
3
2
N
k
B
=
3
2
n
R
C_V=(\frac{\partial U}{\partial T})_V=\frac{3}{2}Nk_B=\frac{3}{2}nR
CV=(∂T∂U)V=23NkB=23nR
定压热容
C
p
=
(
∂
H
∂
T
)
p
=
5
2
N
k
B
=
5
2
n
R
C_p=(\frac{\partial H}{\partial T})_p=\frac{5}{2}Nk_B=\frac{5}{2}nR
Cp=(∂T∂H)p=25NkB=25nR
热容比
γ
=
C
p
C
V
=
5
3
\gamma=\frac{C_p}{C_V}=\frac{5}{3}
γ=CVCp=35
点击蓝色,对比单原子经典理想气体与理想Bose与Fermi气体的异同
理想气体熵的Sackur-Tetrode方程
S
=
N
k
B
(
l
n
V
(
2
π
m
k
B
T
)
3
/
2
N
(
2
π
ℏ
)
3
+
5
2
)
=
N
k
B
(
3
l
n
d
λ
T
+
5
2
)
S=Nk_B(ln\frac{V(2\pi mk_BT)^{3/2}}{N(2\pi \hbar)^3}+\frac{5}{2})=Nk_B(3ln\frac{d}{\lambda_T}+\frac{5}{2})
S=NkB(lnN(2πℏ)3V(2πmkBT)3/2+25)=NkB(3lnλTd+25)
表明理想气体熵随温度的增长趋势为
S
S
S~
l
n
T
lnT
lnT
双原子分子理想气体
双原子分子的自由度:分子质心的3个平动自由度;两原子绕过质心的轴旋转的2个自由度;沿两个原子连线振动的1个自由度。总共6个自由度。
双原子分子的能量
ϵ
=
1
2
M
(
p
x
2
+
p
y
2
+
p
z
2
)
=
1
2
I
(
p
θ
2
+
1
s
i
n
2
θ
p
ϕ
2
)
+
p
r
2
2
μ
+
V
(
r
)
\epsilon=\frac{1}{2M}(p_x^2+p_y^2+p_z^2)=\frac{1}{2I}(p_{\theta}^2+\frac{1}{sin^2\theta}p_{\phi}^2)+\frac{p_r^2}{2\mu}+V(r)
ϵ=2M1(px2+py2+pz2)=2I1(pθ2+sin2θ1pϕ2)+2μpr2+V(r)
C
V
=
5
2
N
k
B
,
C
p
=
7
2
N
k
B
,
γ
=
7
5
C_V=\frac{5}{2}Nk_B,C_p=\frac{7}{2}Nk_B,\gamma=\frac{7}{5}
CV=25NkB,Cp=27NkB,γ=57
平动 translation | 转动 retate | 振动 vibration | |
---|---|---|---|
能量 | 1 2 M ( p x 2 + p y 2 + p z 2 ) \frac{1}{2M}(p_x^2+p_y^2+p_z^2) 2M1(px2+py2+pz2) |
l
(
l
+
1
)
ℏ
2
2
I
(
l
=
0
,
1
,
2
,
.
.
.
)
\frac{l(l+1)\hbar^2}{2I}(l=0,1,2,...)
2Il(l+1)ℏ2(l=0,1,2,...) g l = 2 l + 1 g_l=2l+1 gl=2l+1 |
(
n
+
1
2
)
ℏ
ω
(
n
=
0
,
1
,
2
,
.
.
.
)
(n+\frac{1}{2})\hbar\omega (n=0,1,2,...)
(n+21)ℏω(n=0,1,2,...) g v = 1 g_v=1 gv=1 |
特征温度 |
θ
r
:
=
ℏ
2
2
I
k
B
\theta_r:=\frac{\hbar^2}{2Ik_B}
θr:=2IkBℏ2 10 10 10~ 100 K 100K 100K |
θ
v
:
=
ℏ
ω
k
B
\theta_v:=\frac{\hbar\omega}{k_B}
θv:=kBℏω 1 0 3 10^3 103~ 1 0 4 K 10^4K 104K | |
配分函数 | |||
热容 | |||
低温极限 | |||
高温极限 |
g
=
g
l
g
r
g
v
g=g_lg_rg_v
g=glgrgv
Z
=
Z
t
Z
r
Z
v
Z=Z_tZ_rZ_v
Z=ZtZrZv
为何原子中的电子、质子和种子对理想气体的热容没有贡献?
答:原子中的电子能级差在eV级别,1eV对应的特征温度为
1
0
4
K
10^4K
104K,是一个很高的温度,常温下,电子几乎都处于基态,因此对热容没有贡献。
e α ⩾ 1 e^{\alpha}\geqslant1 eα⩾1
单原子分子 | 理想Bose气体 | 理想Fermi气体 |
---|---|---|
态密度 | 2 π g V ( 2 π ℏ ) 3 ( 2 m ) 3 / 2 ϵ 1 / 2 \frac{2\pi gV}{(2\pi \hbar )^3}(2m)^{3/2}\epsilon^{1/2} (2πℏ)32πgV(2m)3/2ϵ1/2 | 2 π g V ( 2 π ℏ ) 3 ( 2 m ) 3 / 2 ϵ 1 / 2 \frac{2\pi gV}{(2\pi \hbar )^3}(2m)^{3/2}\epsilon^{1/2} (2πℏ)32πgV(2m)3/2ϵ1/2 |
巨配分函数严格解 | l n Ξ = g V λ T 3 ∑ n = 1 ∞ z n n 5 / 2 = g V λ T 3 g 5 / 2 ( z ) ln\Xi=\frac{gV}{\lambda_T^3}\sum_{n=1}^{\infty} {\frac{z^n}{n^{5/2}}}=\frac{gV}{\lambda_T^3}g_{5/2}(z) lnΞ=λT3gV∑n=1∞n5/2zn=λT3gVg5/2(z) | l n Ξ = g V λ T 3 ∑ n = 1 ∞ ( − 1 ) n − 1 z n n 5 / 2 = g V λ T 3 f 5 / 2 ( z ) ln\Xi=\frac{gV}{\lambda_T^3}\sum_{n=1}^{\infty} {\frac{(-1)^{n-1}z^n}{n^{5/2}}}=\frac{gV}{\lambda_T^3}f_{5/2}(z) lnΞ=λT3gV∑n=1∞n5/2(−1)n−1zn=λT3gVf5/2(z) |
BE(FD)函数 | g v ( z ) : = ∑ n = 1 ∞ z n n v g_v(z):=\sum_{n=1}^{\infty}\frac{z^n}{n^v} gv(z):=∑n=1∞nvzn | f v ( z ) : = ∑ n = 1 ∞ ( − 1 ) n − 1 z n n v f_v(z):=\sum_{n=1}^{\infty}\frac{(-1)^{n-1}z^n}{n^v} fv(z):=∑n=1∞nv(−1)n−1zn |
理想xx气体状态方程 | p V N k B T = g 5 / 2 ( z ) g 3 / 2 ( z ) \frac{pV}{Nk_BT}=\frac{g_{5/2}(z)}{g_{3/2}(z)} NkBTpV=g3/2(z)g5/2(z) | p V N k B T = f 5 / 2 ( z ) f 3 / 2 ( z ) \frac{pV}{Nk_BT}=\frac{f_{5/2}(z)}{f_{3/2}(z)} NkBTpV=f3/2(z)f5/2(z) |
内能 | U 3 2 N K B T = g 5 / 2 ( z ) g 3 / 2 ( z ) \frac{U}{\frac{3}{2}NK_BT}=\frac{g_{5/2}(z)}{g_{3/2}(z)} 23NKBTU=g3/2(z)g5/2(z) | U 3 2 N K B T = f 5 / 2 ( z ) f 3 / 2 ( z ) \frac{U}{\frac{3}{2}NK_BT}=\frac{f_{5/2}(z)}{f_{3/2}(z)} 23NKBTU=f3/2(z)f5/2(z) |
Riemam
ζ
\zeta
ζ函数
ζ
(
v
)
=
g
v
(
1
)
=
∑
n
=
1
∞
1
n
v
\zeta(v)=g_v(1)=\sum_{n=1}^{\infty}{\frac{1}{n^v}}
ζ(v)=gv(1)=n=1∑∞nv1
经典理想气体,单原子分子理想Bose与Fermi气体都有
p
=
2
U
3
V
p=\frac{2U}{3V}
p=3V2U
弱简并性条件:
e
α
>
1
e^{\alpha}>1
eα>1
计算过程:
z<1
单原子分子弱简并理想量子气体状态方程与内能:
p
V
N
K
B
T
=
1
∓
1
2
5
/
2
g
N
λ
T
3
V
\frac{pV}{NK_BT}=1\mp\frac{1}{2^{5/2}g}\frac{N\lambda_T^3}{V}
NKBTpV=1∓25/2g1VNλT3
U
3
2
N
k
B
T
=
1
∓
1
2
5
/
2
g
N
λ
T
3
V
\frac{U}{\frac{3}{2}Nk_BT}=1\mp\frac{1}{2^{5/2}g}\frac{N\lambda_T^3}{V}
23NkBTU=1∓25/2g1VNλT3
其中
λ
T
=
2
π
ℏ
2
π
m
k
B
T
\lambda_T=\frac{2\pi\hbar}{\sqrt{2\pi mk_BT}}
λT=2πmkBT2πℏ
统计关联:有效的相互作用
对于Bose系统,其压强小,内能小,这表明玻色子之间存在有效吸引作用。相对于经典粒子,玻色子更倾向于分布在较低能级上。
对于Fermi系统,其压强大,内能大,这表明费米子之间存在有效排斥作用。相对于经典粒子,费米子更倾向于分布在较高能级上。
这种有效的相互作用称为量子系统的统计关联。
e
α
∼
>
1
e^{\alpha}\overset{>}{\sim}1
eα∼>1
Bose-Einstein凝聚
光子气体
首先探讨热辐射系统的相关知识。
受热固体会辐射电磁波,称为热辐射。一般情形下热辐射的强度和强度按频率的分布与辐射体的温度和性质都有关系。如果辐射体对电磁波的吸收和辐射达到平衡,热辐射特性将只取决于温度,与辐射体的其他特性无关,称为平衡辐射,也称黑体辐射。
如何证明空窑辐射的内能和内能密度按频率的分布只与温度有关?
p65
辐射压强
p
p
p与辐射能量密度
u
u
u之间的关系
p
=
1
3
u
p=\frac{1}{3}u
p=31u
列别捷夫的光压强实验可以导出,也可由热力学与统计物理理论导出。证明过程
将空窖内平衡辐射看作热力徐系统,选取温度T与体积V为状态参量。由于空窖是均匀的,内能密度只能是温度的函数。则
U
(
T
,
V
)
=
u
(
T
)
V
U(T,V)=u(T)V
U(T,V)=u(T)V
利用热力学公式证明及记忆过程稍后补充
(
∂
U
∂
V
)
T
=
T
(
∂
p
∂
T
)
V
−
p
(\frac{\partial U}{\partial V})_T=T(\frac{\partial p}{\partial T})_V-p
(∂V∂U)T=T(∂T∂p)V−p
代入以上式子及利用p对T的偏导等于p对u,u再对T的偏导,积分即可得到
u
=
a
T
4
,
U
=
a
T
4
V
u=aT^4,U=aT^4V
u=aT4,U=aT4V
将以上所得代入
d
S
=
d
U
+
p
d
V
T
dS=\frac{dU+pdV}{T}
dS=TdU+pdV
可得到
S
=
4
3
a
T
3
V
S=\frac{4}{3}aT^3V
S=34aT3V
可逆绝热过程,熵不变
将以上代入
G
=
U
−
T
S
+
p
V
G=U-TS+pV
G=U−TS+pV可得
G
=
0
G=0
G=0(原因是光子数不守恒)
斯特藩-玻尔兹曼定律证明过程暂缺p66
J
=
1
4
c
u
=
1
4
c
a
T
4
=
σ
T
4
J=\frac{1}{4}cu=\frac{1}{4}caT^4=\sigma T^4
J=41cu=41caT4=σT4
斯特藩常量
σ
=
5.669
×
1
0
−
8
W
m
−
2
K
−
4
\sigma=5.669\times10^{-8}Wm^{-2}K^{-4}
σ=5.669×10−8Wm−2K−4
根据粒子的观点,可以将空窖内的辐射场看作光子气体。由于空窖内的辐射场可以分解为无穷多个单色平面波的叠加对于单色平面波点这里了解,待补充。
具有一定圆频率
ω
\omega
ω和波矢
k
k
k的单色平面波与具有一定能量
ϵ
\epsilon
ϵ和动量
p
p
p的光子相对应。即德布罗意关系
p
=
ℏ
k
,
ϵ
=
ℏ
ω
p=\hbar k,\epsilon=\hbar\omega
p=ℏk,ϵ=ℏω
由
ω
=
c
k
\omega=ck
ω=ck可以得到光子的能量动量关系
ϵ
=
c
p
\epsilon=cp
ϵ=cp
或者由相对论质能关系
ϵ
=
c
2
p
2
+
m
2
c
4
\epsilon=\sqrt{c^2p^2+m^2c^4}
ϵ=c2p2+m2c4及光子静质量
m
=
0
m=0
m=0得到光子的能量动量关系。
量子态密度 D ( ϵ ) D(\epsilon) D(ϵ)
配分函数是热力学与统计物理当中最重要的一个物理量,计算配分函数的一种方法便是量子态密度法,这里给出量子态密度的解论公式以及计算方法。
D ( ϵ ) D(\epsilon) D(ϵ) | 总体积 | 体积元 | ϵ = c p \epsilon = cp ϵ=cp | ϵ = p 2 2 m \epsilon = \frac{p^2}{2m} ϵ=2mp2 | ϵ = ( p c ) 2 + ( m c 2 ) 2 \epsilon =\sqrt{(pc)^2 +(mc^2)^2} ϵ=(pc)2+(mc2)2 |
---|---|---|---|---|---|
三维 | 4 π 3 p 3 \frac{4\pi}{3}p^3 34πp3 | ( 2 π ℏ / L ) 3 (2\pi \hbar /L)^3 (2πℏ/L)3 | 4 π g V ( 2 π ℏ c ) 3 ϵ 2 \frac{4\pi gV}{(2\pi \hbar c)^3}\epsilon^2 (2πℏc)34πgVϵ2 | 2 π g V ( 2 π ℏ ) 3 ( 2 m ) 3 / 2 ϵ 1 / 2 \frac{2\pi gV}{(2\pi \hbar )^3}(2m)^{3/2}\epsilon^{1/2} (2πℏ)32πgV(2m)3/2ϵ1/2 | |
二维 | π p 2 \pi p^2 πp2 | ( 2 π ℏ / L ) 2 (2\pi \hbar /L)^2 (2πℏ/L)2 | 2 π g A ( 2 π ℏ c ) 2 ϵ \frac{2\pi gA}{(2\pi \hbar c)^2}\epsilon (2πℏc)22πgAϵ | ||
一维 | 2p | ( 2 π ℏ / L ) 1 (2\pi \hbar /L)^1 (2πℏ/L)1 | 2 g L 2 π ℏ c \frac{2gL}{2\pi \hbar c} 2πℏc2gL |
注:g=2s+1 ,需要考虑自旋。
计算方法
1、根据总量子态数
N
(
ϵ
)
N(\epsilon)
N(ϵ)及
D
(
ϵ
)
=
d
N
(
ϵ
)
/
d
ϵ
D(\epsilon)=dN(\epsilon)/d\epsilon
D(ϵ)=dN(ϵ)/dϵ求出,总体积中的
p
p
p需根据能量动量关系代入能量。
N
(
ϵ
)
=
动
量
空
间
总
体
积
/
每
个
量
子
态
体
积
元
N(\epsilon)=动量空间总体积/每个量子态体积元
N(ϵ)=动量空间总体积/每个量子态体积元
2、从动量空间开始,这里以三维为例子。
动量空间中的球坐标系
p
x
=
p
s
i
n
θ
c
o
s
ϕ
p_x=psin\theta cos\phi
px=psinθcosϕ
p
y
=
p
s
i
n
θ
s
i
n
ϕ
p_y=psin\theta sin\phi
py=psinθsinϕ
p
z
=
p
c
o
s
θ
p_z=pcos\theta
pz=pcosθ
体积元为
p
2
s
i
n
2
θ
d
p
d
θ
d
ϕ
p^2sin^2\theta dpd\theta d\phi
p2sin2θdpdθdϕ
体积V内即
p
到
p
+
Δ
p
,
θ
到
θ
+
Δ
θ
,
ϕ
到
ϕ
+
Δ
ϕ
p到p+\Delta p,\theta到\theta+\Delta \theta,\phi到\phi+\Delta\phi
p到p+Δp,θ到θ+Δθ,ϕ到ϕ+Δϕ内的总微观状态数为
V
p
2
s
i
n
2
θ
d
p
d
θ
d
ϕ
/
h
3
Vp^2sin^2\theta dpd\theta d\phi/h^3
Vp2sin2θdpdθdϕ/h3
对
θ
和
ϕ
\theta和\phi
θ和ϕ积分即可得到
p
到
p
+
Δ
p
p到p+\Delta p
p到p+Δp内的微观状态数
4
π
V
h
3
p
2
d
p
\frac{4\pi V}{h^3}p^2dp
h34πVp2dp
之后根据能量和动量关系,用能量替换动量即可得到量子态密度,不再赘述。
相关知识点
1、相对论性自由粒子能量与动量关系
ϵ
=
c
p
\epsilon=cp
ϵ=cp证明
根据德布罗意关系
p
=
ℏ
k
,
ϵ
=
ℏ
ω
p=\hbar k,\epsilon=\hbar\omega
p=ℏk,ϵ=ℏω
以及
ω
=
c
k
\omega =ck
ω=ck得出
做功的统计表达
d
W
=
∑
l
a
l
d
ϵ
l
dW=\sum_l {a_ld\epsilon_l}
dW=∑laldϵl
热量的统计表达
d
Q
=
∑
l
ϵ
l
d
a
l
dQ=\sum_l {\epsilon_lda_l}
dQ=∑lϵldal
β
\beta
β与
1
/
T
1/T
1/T是
d
Q
dQ
dQ的积分因子
α
=
−
μ
k
B
T
,
β
=
1
k
B
T
\alpha=-\frac{\mu}{k_BT},\beta=\frac{1}{k_BT}
α=−kBTμ,β=kBT1
Boltzmann常数:
k
B
:
=
1.380
,
649
×
1
0
−
23
J
K
−
1
k_B:=1.380,649\times10^{-23}JK^{-1}
kB:=1.380,649×10−23JK−1
热波长
λ
T
=
2
π
ℏ
2
π
m
k
B
T
\lambda_T=\frac{2\pi\hbar}{\sqrt{2\pi mk_BT}}
λT=2πmkBT2πℏ
单原子分子理想气体内能
U
=
N
ϵ
‾
=
3
2
N
k
B
T
=
3
2
n
R
T
U=N\overline{\epsilon}=\frac{3}{2}Nk_BT=\frac{3}{2}nRT
U=Nϵ=23NkBT=23nRT
定容热容
C
V
=
(
∂
U
∂
T
)
V
=
3
2
N
K
B
=
3
2
n
R
C_V=(\frac{\partial U}{\partial T})_V=\frac{3}{2}NK_B=\frac{3}{2}nR
CV=(∂T∂U)V=23NKB=23nR
定压热容
C
p
=
(
∂
H
∂
T
)
p
=
5
2
N
K
B
=
5
2
n
R
C_p=(\frac{\partial H}{\partial T})_p=\frac{5}{2}NK_B=\frac{5}{2}nR
Cp=(∂T∂H)p=25NKB=25nR
热容比$
γ
=
C
p
C
V
\gamma=\frac{C_p}{C_V}
γ=CVCp
对于Bose系统,
e
α
>
=
1
e^\alpha>=1
eα>=1
微正则系综 | 正则系综 | 巨正则系综 | |
---|---|---|---|
定义 | 系统为孤立系,N、V、E不变 | N,V,T不变 | |
特性函数 | 熵或总微观状态数 | Z : = ∑ s e − β ϵ s Z:=\sum_s e^{-\beta\epsilon_s} Z:=∑se−βϵs | |
概率分布 | ρ = 1 Ω \rho=\frac{1}{\Omega} ρ=Ω1 | P s = e − β ϵ s Z P_s=\frac{e^{-\beta\epsilon_s}}{Z} Ps=Ze−βϵs | |
内能 U U U | − ∂ ∂ β l n Z -\frac{\partial }{\partial \beta}lnZ −∂β∂lnZ | ||
广义力 Y λ Y_{\lambda} Yλ | − 1 β ∂ ∂ y λ l n Z -\frac{1}{\beta}\frac{\partial }{\partial y_{\lambda}}lnZ −β1∂yλ∂lnZ | ||
熵 S S S | k B ( l n Z − β ∂ ∂ β ) l n Z k_B(lnZ-\beta\frac{\partial }{\partial \beta})lnZ kB(lnZ−β∂β∂)lnZ | ||
Helmholtz自由能 F F F | − k B T l n Z -k_BTlnZ −kBTlnZ |
计算微正则系综热力学函数的流程
1、先计算系统总微观状态数
Ω
(
N
,
V
,
E
)
\Omega(N,V,E)
Ω(N,V,E),并得到系统的熵
S
=
k
B
l
n
Ω
(
N
,
V
,
E
)
S=k_Bln\Omega(N,V,E)
S=kBlnΩ(N,V,E)反解出
E
=
E
(
N
,
V
,
S
)
E=E(N,V,S)
E=E(N,V,S);
2、由
T
=
(
∂
E
/
∂
S
)
N
,
V
T=(\partial E/\partial S)_{N,V}
T=(∂E/∂S)N,V,求出
T
=
T
(
N
,
V
,
S
)
T=T(N,V,S)
T=T(N,V,S)反解出
S
=
S
(
N
,
V
,
T
)
S=S(N,V,T)
S=S(N,V,T);代入E的得出
E
=
E
(
N
,
V
,
T
)
E=E(N,V,T)
E=E(N,V,T);
3、由
p
=
−
(
∂
E
/
∂
V
)
N
,
S
p=-(\partial E/\partial V)_{N,S}
p=−(∂E/∂V)N,S,求出
p
=
p
(
N
,
V
,
S
)
=
p
(
N
,
V
,
T
)
p=p(N,V,S)=p(N,V,T)
p=p(N,V,S)=p(N,V,T)。
至此,p、E、S均化为了N,V,T的函数,然后就可求出一切热力学量。
计算正则系综热力学函数流程
在热力学中,证明过以N,V,T为独立变量时,Helmholtz自由能F(N,V,T)是系统的特性函数,因此求得Z或F即可求得系统的一切热力学量。
一个物理量涨落的定义
Δ
O
=
O
−
O
‾
\Delta O=O-\overline{O}
ΔO=O−O
方差:
v
a
r
(
O
)
=
O
2
‾
−
O
‾
2
var(O)=\overline{O^2}-\overline{O}^2
var(O)=O2−O2
正则系综中能量的方差
v
a
r
(
E
)
=
k
B
T
C
V
var(E)=k_BTC_V
var(E)=kBTCV
系统能量的相对涨落
Δ
E
E
∼
(
E
−
E
‾
)
2
‾
E
‾
2
∝
1
N
\frac{\Delta E}{E}\sim\sqrt{\frac{\overline{(E-\overline{E})^2}}{\overline{E}^2}}\propto\frac{1}{\sqrt{N}}
EΔE∼E2(E−E)2∝N1表明系统能量的相对涨落随这系统的粒子数增加而衰减。因此只有系统的粒子数足够多,各种物理量的涨落才可以忽略。
正则系综推导范德瓦尔斯方程
相关证明
1、箱中的自由粒子(一维情况)
确定粒子的运动状态,首先确定德布罗意波在器壁的边界条件(有驻波条件或周期性边界条件)
周期性边界条件
L
=
∣
n
x
∣
λ
,
∣
n
x
∣
=
0
,
1
,
2
,
.
.
.
L=|n_x|\lambda,|n_x|=0,1,2,...
L=∣nx∣λ,∣nx∣=0,1,2,...
考虑波矢与波长的关系,及波在空间中两个传播方向
k
x
=
2
π
L
n
x
,
n
x
=
0
,
±
1
,
±
2
,
.
.
.
k_x=\frac{2\pi}{L}n_x,n_x=0,\pm1,\pm2,...
kx=L2πnx,nx=0,±1,±2,...代入
p
=
ℏ
k
p=\hbar k
p=ℏk得到
p
x
=
2
π
ℏ
L
n
x
,
n
x
=
0
,
±
1
,
±
2
,
.
.
.
p_x=\frac{2\pi\hbar}{L}n_x,n_x=0,\pm1,\pm2,...
px=L2πℏnx,nx=0,±1,±2,...
对于非相对论的自由粒子,
ϵ
=
p
2
/
2
m
\epsilon=p^2/2m
ϵ=p2/2m
ϵ
=
1
2
m
(
2
π
ℏ
L
)
2
(
n
x
2
+
n
y
2
+
n
z
2
)
\epsilon=\frac{1}{2m}(\frac{2\pi\hbar}{L})^2(n_x^2+n_y^2+n_z^2)
ϵ=2m1(L2πℏ)2(nx2+ny2+nz2)
对于相对论性粒子,
ϵ
=
c
p
\epsilon=cp
ϵ=cp
ϵ
=
2
π
ℏ
c
L
(
n
x
2
+
n
y
2
+
n
z
2
)
1
/
2
\epsilon=\frac{2\pi\hbar c}{L}(n_x^2+n_y^2+n_z^2)^{1/2}
ϵ=L2πℏc(nx2+ny2+nz2)1/2
讨论:若粒子局域在微观大小的空间中运动,则能量会显著分立,能级由
(
n
x
2
+
n
y
2
+
n
z
2
)
(n_x^2+n_y^2+n_z^2)
(nx2+ny2+nz2)给出,有简并度。
宏观状态运动 【p172】
相对论性理想Bose或Fermi气体压强与内能的关系。
相对论性气体的态密度为
D
(
ϵ
)
=
4
π
g
V
(
2
π
ℏ
c
)
2
ϵ
2
D(\epsilon)=\frac{4\pi gV}{(2\pi\hbar c)^2}\epsilon^2
D(ϵ)=(2πℏc)24πgVϵ2
Bose系统巨配分函数
l
n
Ξ
=
−
∑
l
g
l
l
n
(
1
−
e
−
α
−
β
ϵ
)
=
−
∫
D
(
ϵ
)
l
n
(
1
−
e
−
α
−
β
ϵ
)
d
ϵ
ln\Xi=-\sum_l g_lln(1-e^{-\alpha-\beta\epsilon})=-\int D(\epsilon)ln(1-e^{-\alpha-\beta\epsilon})d\epsilon
lnΞ=−l∑glln(1−e−α−βϵ)=−∫D(ϵ)ln(1−e−α−βϵ)dϵ
Fermi系统巨配分函数
l
n
Ξ
=
∑
l
g
l
l
n
(
1
+
e
−
α
−
β
ϵ
)
=
∫
D
(
ϵ
)
l
n
(
1
+
e
−
α
−
β
ϵ
)
d
ϵ
ln\Xi=\sum_l g_lln(1+e^{-\alpha-\beta\epsilon})=\int D(\epsilon)ln(1+e^{-\alpha-\beta\epsilon})d\epsilon
lnΞ=l∑glln(1+e−α−βϵ)=∫D(ϵ)ln(1+e−α−βϵ)dϵ
以及
p
=
1
β
∂
∂
V
l
n
Ξ
p=\frac{1}{\beta}\frac{\partial }{\partial V}ln\Xi
p=β1∂V∂lnΞ
U
=
−
∂
∂
β
l
n
Ξ
U=-\frac{\partial }{\partial \beta}ln\Xi
U=−∂β∂lnΞ
即可得到
p
=
U
3
V
p=\frac{U}{3V}
p=3VU也是热辐射系统压强与内能关系