双原子分子的自由度：分子质心的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θ1​pϕ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​=25​NkB​，Cp​=27​NkB​，γ=57​

g = g l g r g v g=g_lg_rg_v g=gl​gr​gv​ Z = Z t Z r Z v Z=Z_tZ_rZ_v Z=Zt​Zr​Zv​ 为何原子中的电子、质子和种子对理想气体的热容没有贡献？ 答：原子中的电子能级差在eV级别，1eV对应的特征温度为 1 0 4 K 10^4K 104K,是一个很高的温度，常温下，电子几乎都处于基态，因此对热容没有贡献。

e α ⩾ 1 e^{\alpha}\geqslant1 eα⩾1

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} NKB​TpV​=1∓25/2g1​VNλ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} 23​NkB​TU​=1∓25/2g1​VNλT3​​ 其中 λ T = 2 π ℏ 2 π m k B T \lambda_T=\frac{2\pi\hbar}{\sqrt{2\pi mk_BT}} λT​=2πmkB​T ​2πℏ​

统计关联：有效的相互作用   对于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=31​u 列别捷夫的光压强实验可以导出，也可由热力学与统计物理理论导出。证明过程 将空窖内平衡辐射看作热力徐系统，选取温度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=34​aT3V 可逆绝热过程，熵不变 将以上代入 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=41​cu=41​caT4=σ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得到光子的能量动量关系。

  配分函数是热力学与统计物理当中最重要的一个物理量，计算配分函数的一种方法便是量子态密度法，这里给出量子态密度的解论公式以及计算方法。

注：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πV​p2dp 之后根据能量和动量关系，用能量替换动量即可得到量子态密度，不再赘述。

相关知识点 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=∑l​al​dϵl​ 热量的统计表达 d Q = ∑ l ϵ l d a l dQ=\sum_l {\epsilon_lda_l} dQ=∑l​ϵl​dal​ β \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} α=−kB​Tμ​,β=kB​T1​ 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πmkB​T ​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ϵ=23​NkB​T=23​nRT 定容热容 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​=23​NKB​=23​nR 定压热容 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​=25​NKB​=25​nR 热容比$ γ = C p C V \gamma=\frac{C_p}{C_V} γ=CV​Cp​​ 对于Bose系统， e α > = 1 e^\alpha>=1 eα>=1

计算微正则系综热力学函数的流程 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=kB​lnΩ(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)=kB​TCV​ 系统能量的相对涨落 Δ 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​​ ​∝N ​1​表明系统能量的相对涨落随这系统的粒子数增加而衰减。因此只有系统的粒子数足够多，各种物理量的涨落才可以忽略。

正则系综推导范德瓦尔斯方程

相关证明 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∑​gl​ln(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∑​gl​ln(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​也是热辐射系统压强与内能关系