统计热力学，最不化学的一章。
（ 本节涉及大量数学运算，计算细节均省略，仅给出最终结果。 ）
一切粒子都在不停地运动，每个粒子的运动都遵守力学规律。但是想用力学中的微分方程去描述系统的运动状态是不可能的，必须用统计学的方法，推求极大数目粒子运动的统计平均。人们最早是用统计的方法来研究气体的行为，现在通称为气体分子动理论，以后才发展为统计力学（statistical mechanics）。
统计力学把系统的微观性质和宏观性质联系起来，从这个角度来说，又可以称为统计热力学（statistical thermodynamics）。
根据对物质结构的某些基本假定，以及实验获得的光谱数据，可以求出物质的一些基本常数，利用这些数据可以算出配分函数（partition function），然后再求出物质的热力学性质。这就是统计热力学的基本任务。
对于简单分子，使用统计热力学得到的结果相当令人满意，但是对于大的游离分子或凝聚系统，应用统计热力学的结果还存在很大困难。
本书主要介绍Boltzmann统计，但是已经采用能量量子化的概念加以修正，再根据情况过渡到经典统计所能适用的公式。
先快速引入一些基本概念：
按照统计单位（粒子）是否可以分辨，把系统分为定位系统（localized system）和非定位系统（non-localized system）。前者的粒子可以彼此分辨，如晶体中相对固定的粒子；后者的粒子不能彼此分辨，如气体分子始终在混乱运动。
按照统计单位之间有无相互作用，把系统分为近独立粒子系统（assembly of independent particles）和非独立粒子系统（assembly of interacting particles）。前者粒子之间的相互作用可以忽略，系统总能量等于各粒子能量之和；后者粒子之间的相互作用不能忽略，是各粒子坐标的函数。
系统的热力学概率Ω \Omega Ω S = k ln Ω S=k\ln\Omega S = k ln Ω
（ 注意Ω \Omega Ω ）
考虑一个系统有N N N N N N
不过无论哪一种分布方式都需要满足如下两个条件：
∑ i N i = N ∑ i N i ϵ i = U \sum_iN_i=N \\[3mm] \sum_iN_i\epsilon_i=U
i ∑ N i = N i ∑ N i ϵ i = U
此时需要考虑如何求解Ω \Omega Ω t m t_m t m Ω \Omega Ω
结合上述条件，通过一系列数学计算，我们就能得到微观状态数最多的一种分布，称为最概然分布（most probable distribution）：
N i ∗ = e α + β ϵ i N_i^*=e^{\alpha+\beta\epsilon_i}
N i ∗ = e α + β ϵ i
上式中还包含两个待定因子α , β \alpha,\beta α , β
N i ∗ = N e − ϵ i / k B T ∑ e − ϵ i / k B T N_i^*=N\frac {e^{-\epsilon_i/k_BT}}{\sum e^{-\epsilon_i/k_BT}}
N i ∗ = N ∑ e − ϵ i / k B T e − ϵ i / k B T
我们还能获得定位系统的熵和Helmholtz自由能的表达式：
S = k B N ln ∑ e − ϵ i / k B T + U T S=k_BN\ln\sum e^{-\epsilon_i/k_BT}+\frac UT
S = k B N ln ∑ e − ϵ i / k B T + T U
A = − N k B T ln ∑ i e − ϵ i / k B T A=-Nk_BT\ln\sum_ie^{-\epsilon_i/k_BT}
A = − N k B T ln i ∑ e − ϵ i / k B T
注意到ϵ i \epsilon_i ϵ i A A A ( T , V , N ) (T,V,N) ( T , V , N ) A A A
在量子力学中，我们把该能级可能有的微观状态数称为该能级的简并度（degeneracy），用g i g_i g i
N i ∗ = N g i e − ϵ i / k B T ∑ i g i e − ϵ i / k B T S 定 位 = k B N ln ∑ i g i e − ϵ i / k B T + U T A 定 位 = − N k B T ln ∑ i g i e − ϵ i / k B T N_i^* = N\frac {g_ie^{-\epsilon_i/k_BT}}{\sum_ig_ie^{-\epsilon_i/k_BT}} \\[3mm] S_{定位}=k_BN\ln\sum_ig_ie^{-\epsilon_i/k_BT}+\frac UT \\[3mm]
A_{定位}=-Nk_BT\ln\sum_ig_ie^{-\epsilon_i/k_BT}
N i ∗ = N ∑ i g i e − ϵ i / k B T g i e − ϵ i / k B T S 定 位 = k B N ln i ∑ g i e − ϵ i / k B T + T U A 定 位 = − N k B T ln i ∑ g i e − ϵ i / k B T
对于非定位系统，我们需要对上述结果做一些修正：
N i ∗ = N g i e − ϵ i / k B T ∑ i g i e − ϵ i / k B T S 非 定 位 = k B ln ( ∑ i g i e − ϵ i / k B T ) N N ! + U T A 非 定 位 = − k B T ln ( ∑ i g i e − ϵ i / k B T ) N N ! N_i^* = N\frac {g_ie^{-\epsilon_i/k_BT}}{\sum_ig_ie^{-\epsilon_i/k_BT}} \\[3mm] S_{非定位}=k_B\ln \frac {\left(\sum_ig_ie^{-\epsilon_i/k_BT} \right)^N}{N!}+\frac UT \\[3mm]
A_{非定位}=-k_BT \ln\frac {\left(\sum_ig_ie^{-\epsilon_i/k_BT} \right)^N}{N!}
N i ∗ = N ∑ i g i e − ϵ i / k B T g i e − ϵ i / k B T S 非 定 位 = k B ln N ! ( ∑ i g i e − ϵ i / k B T ) N + T U A 非 定 位 = − k B T ln N ! ( ∑ i g i e − ϵ i / k B T ) N
无论定位或非定位系统，最概然分布的公式是一样的，但在S S S A A A
Boltzmann公式常转化为不同形式：
N i ∗ N j ∗ = g i e − ϵ i / k B T g j e − ϵ j / k B T \frac {N_i^*}{N_j^*}=\frac {g_ie^{-\epsilon_i/k_BT}}{g_je^{-\epsilon_j/k_BT}}
N j ∗ N i ∗ = g j e − ϵ j / k B T g i e − ϵ i / k B T
假定最低能级为ϵ 0 \epsilon_0 ϵ 0 N 0 N_0 N 0
N i = N 0 e − Δ ϵ i / k B T N_i=N_0e^{-\Delta\epsilon_i/k_BT}
N i = N 0 e − Δ ϵ i / k B T
在推导Boltzmann公式时，曾作以下假定：
在所有的分布方式中，有一种分布方式热力学概率最大，称为最概然分布。
最概然分布的状态数最多，基本上可以用它代替总的微观状态，也就是说最概然分布实质上可以代表一切分布，最概然分布实际上也就是平衡分布。
我们在这里加以说明。
数学推导
根据计算可得：
ln ( t m + δ t t m ) = − 1 2 ∑ t ( δ N i ) 2 N i ∗ \ln \left(\frac {t_m+\delta t}{t_m} \right)=-\frac 12\sum_t\frac {(\delta N_i)^2}{N_i^*}
ln ( t m t m + δ t ) = − 2 1 t ∑ N i ∗ ( δ N i ) 2
上式表明，不论偏差δ N i \delta N_i δ N i t m t_m t m ( t m + δ t ) (t_m+\delta t) ( t m + δ t ) N i ∗ = 3 × 1 0 19 N_i^*=3\times10^{19} N i ∗ = 3 × 1 0 1 9 δ N i / N i ∗ = 0.01 \delta N_i/N_i^*=0.01 δ N i / N i ∗ = 0 . 0 1
t m + δ t t m = exp ( − 3 × 1 0 15 ) \frac {t_m+\delta t}{t_m}=\exp(-3\times10^{15})
t m t m + δ t = exp ( − 3 × 1 0 1 5 )
据此可发现偏移量是极其小的，反过来就证明了最概然分布的概率最大。
进行计算和化简近似后还可以得到：
ln t m ≤ ln Ω ≤ ln t m + ln N \ln t_m \le \ln\Omega \le \ln t_m+\ln N
ln t m ≤ ln Ω ≤ ln t m + ln N
假设N ≈ 1 0 24 N \approx 10^{24} N ≈ 1 0 2 4 ln t m \ln t_m ln t m 1 0 23 10^{23} 1 0 2 3 ln N = 55.262 \ln N=55.262 ln N = 5 5 . 2 6 2 ln N \ln N ln N ln t m = ln Ω \ln t_m = \ln \Omega ln t m = ln Ω
可以发现当N N N
根据量子力学原理，已知基本粒子和由奇数个基本粒子组成的原子和分子必须遵循Pauli不相容原理，即每一个量子状态最多只能容纳一个粒子。但对光子以及由偶数个基本粒子组成的原子和分子则不受Pauli不相容原理的制约，即每个量子状态所能容纳的粒子数没有限制。
这两类粒子适用的统计方法是不同的，对于前者组成的等同粒子系统服从Fermi-Dirac统计，对于后者组成的等同粒子系统则服从Bose-Einstein统计。
B o s e − E i n s t e i n : N i = g i e − α − β ϵ i − 1 或 g i N i = e − α − β ϵ i − 1 Bose-Einstein: \\[3mm] N_i=\frac {g_i}{e^{-\alpha-\beta\epsilon_i}-1} \quad 或 \quad \frac {g_i}{N_i}=e^{-\alpha-\beta\epsilon_i}-1
B o s e − E i n s t e i n : N i = e − α − β ϵ i − 1 g i 或 N i g i = e − α − β ϵ i − 1
F e r m i − D i r a c : N i = g i e − α − β ϵ i + 1 或 g i N i = e − α − β ϵ i + 1 Fermi-Dirac: \\[3mm] N_i=\frac {g_i}{e^{-\alpha-\beta\epsilon_i}+1} \quad 或 \quad \frac {g_i}{N_i}=e^{-\alpha-\beta\epsilon_i}+1
F e r m i − D i r a c : N i = e − α − β ϵ i + 1 g i 或 N i g i = e − α − β ϵ i + 1
这两种统计方法都是建立在量子力学的基础之上，因此通称量子统计。与Boltzmann统计的公式比较一下不难发现，它们只差了一个± 1 \pm1 ± 1 g i ≫ N i g_i \gg N_i g i ≫ N i
本书只讨论Boltzmann统计。
已知最概然分布公式：
N i = N g i e − ϵ i / k B T ∑ i g i e − ϵ i / k B T N_i = N\frac {g_ie^{-\epsilon_i/k_BT}}{\sum_ig_ie^{-\epsilon_i/k_BT}}
N i = N ∑ i g i e − ϵ i / k B T g i e − ϵ i / k B T
令分母：
∑ i g i e − ϵ i / k B T ≡ q \sum_ig_ie^{-\epsilon_i/k_BT} \equiv q
i ∑ g i e − ϵ i / k B T ≡ q
q q q
统计热力学的最重要的任务之一就是通过配分函数来计算系统的热力学函数。因此配分函数可以看作是联系微观和宏观的桥梁。
A 定 位 = − k B T ln q N S 定 位 = N k B [ ∂ ∂ T ( T ln q ) ] V , N = N k B ln q + N k B T ( ∂ ln q ∂ T ) V , N U 定 位 = N k B T 2 ( ∂ ln q ∂ T ) V , N G 定 位 = A + p V = − k B T ln q N + N k B T V ( ∂ ln q ∂ V ) T , N H 定 位 = U + p V = N k B T 2 ( ∂ ln q ∂ T ) V , N + N k B T V ( ∂ ln q ∂ T ) T , N C V , 定 位 = ∂ ∂ T [ N k B T 2 ( ∂ ln q ∂ T ) V , N ] V A_{定位}=-k_BT\ln q^N \\[3mm]
S_{定位}=Nk_B \left[\frac {\partial}{\partial T}(T\ln q) \right]_{V,N}=Nk_B\ln q+Nk_BT \left(\frac {\partial \ln q}{\partial T} \right)_{V,N} \\[3mm]
U_{定位}=Nk_BT^2 \left(\frac {\partial \ln q}{\partial T} \right)_{V,N} \\[3mm]
G_{定位}=A+pV=-k_BT\ln q^N+Nk_BTV \left(\frac {\partial \ln q}{\partial V} \right)_{T,N} \\[3mm]
H_{定位}=U+pV=Nk_BT^2 \left(\frac {\partial \ln q}{\partial T} \right)_{V,N}+Nk_BTV \left(\frac {\partial \ln q}{\partial T} \right)_{T,N} \\[3mm]
C_{V,定位}=\frac {\partial}{\partial T} \left[Nk_BT^2 \left(\frac {\partial \ln q}{\partial T} \right)_{V,N} \right]_{V}
A 定 位 = − k B T ln q N S 定 位 = N k B [ ∂ T ∂ ( T ln q ) ] V , N = N k B ln q + N k B T ( ∂ T ∂ ln q ) V , N U 定 位 = N k B T 2 ( ∂ T ∂ ln q ) V , N G 定 位 = A + p V = − k B T ln q N + N k B T V ( ∂ V ∂ ln q ) T , N H 定 位 = U + p V = N k B T 2 ( ∂ T ∂ ln q ) V , N + N k B T V ( ∂ T ∂ ln q ) T , N C V , 定 位 = ∂ T ∂ [ N k B T 2 ( ∂ T ∂ ln q ) V , N ] V
接着我们来进行配分函数的分离。一个分子的能量可以认为是分子整体运动的能量即平动能（ϵ t \epsilon_t ϵ t ϵ r \epsilon_r ϵ r ϵ v \epsilon_v ϵ v ϵ e \epsilon_e ϵ e ϵ n \epsilon_n ϵ n
这几个能级的大小次序是：
ϵ n > ϵ e > ϵ v > ϵ r > ϵ t \epsilon_n>\epsilon_e>\epsilon_v>\epsilon_r>\epsilon_t
ϵ n > ϵ e > ϵ v > ϵ r > ϵ t
可以推得：
q = q t ⋅ q r ⋅ q v ⋅ q e ⋅ q n q=q_t\cdot q_r \cdot q_v \cdot q_e \cdot q_n
q = q t ⋅ q r ⋅ q v ⋅ q e ⋅ q n
原子核配分函数可以写作（近似后结果）：
q n = g n , 0 exp ( − ϵ n , 0 k B T ) q_n=g_{n,0}\,\exp \left(-\frac {\epsilon_{n,0}}{k_BT} \right)
q n = g n , 0 exp ( − k B T ϵ n , 0 )
除了核反应外，通常的化学反应过程中，原子核总是处于基态，若将核基态能量选作零，则：
q n = g n , 0 q_n=g_{n,0}
q n = g n , 0
核能级的简并度来源自原子核的自旋，但是自旋方向不同的各态之间没有显著能量差别，只有在超精细结构中才能看出。
电子配分函数可以写作（近似后结果，且规定最低能态能量为零）：
q e = g e , 0 q_e=g_{e,0}
q e = g e , 0
电子绕核运动的简并度来源自电子的总动量矩，与核能级简并度类似。
核配分函数与电子配分函数对于U , H , C V U,H,C_V U , H , C V A , G , S A,G,S A , G , S
注意到有些原子中，电子的基态和第一激发态之间的间隔不是很大，此时需要考虑q e q_e q e
平动配分函数可以写作：
q t = ∑ i g i , t exp ( − ϵ i , t k B T ) q_t=\sum_ig_{i,t}\,\exp \left(-\frac {\epsilon_{i,t}}{k_BT} \right)
q t = i ∑ g i , t exp ( − k B T ϵ i , t )
根据势箱模型，可以从波动方程求出平动能的能量ϵ i , t \epsilon_{i,t} ϵ i , t
q t = ( 2 π m k B T h 2 ) 3 / 2 ⋅ V q_t=\left(\frac {2\pi mk_BT}{h^2} \right)^{3/2} \cdot V
q t = ( h 2 2 π m k B T ) 3 / 2 ⋅ V
平动配分函数对热力学函数的贡献如下：
A t = − N k B T ln [ ( 2 π m k B T h 2 ) 3 / 2 V ] + N k B T ln N − N k B T S t = N k B ( ln q t N + 5 2 ) A_t=-Nk_BT\ln \left[\left(\frac {2\pi mk_BT}{h^2} \right)^{3/2} V\right]+Nk_BT\ln N-Nk_BT \\[3mm]
S_t=Nk_B\left(\ln\frac {q_t}{N}+\frac 52 \right)
A t = − N k B T ln [ ( h 2 2 π m k B T ) 3 / 2 V ] + N k B T ln N − N k B T S t = N k B ( ln N q t + 2 5 )
这个公式称为Sackur-Tetrode公式，可用来计算理想气体的平动熵。
对于双原子分子，转动配分函数可以写作（忽略转动和振动的相互作用）：
q r = 8 π 2 I k B T σ h 2 = T Θ r σ q_r=\frac {8\pi^2Ik_BT}{\sigma h^2}=\frac {T}{\Theta_r\sigma}
q r = σ h 2 8 π 2 I k B T = Θ r σ T
式中σ \sigma σ σ = 2 \sigma=2 σ = 2 Θ r \Theta_r Θ r K K K
（ 对称数也就是对称轴的最高轴次，或者说是主轴轴次。 ）
对于转动特征温度较高的分子，经常需考虑激发态的量子态。
对于非线性多原子分子，可以证明其转动配分函数为：
q r = 8 π 2 ( 2 π k B T ) 3 / 2 σ h 3 ⋅ ( I x ⋅ I y ⋅ I z ) 1 / 2 q_r=\frac {8\pi^2(2\pi k_BT)^{3/2}}{\sigma h^3}\cdot (I_x\cdot I_y \cdot I_z)^{1/2}
q r = σ h 3 8 π 2 ( 2 π k B T ) 3 / 2 ⋅ ( I x ⋅ I y ⋅ I z ) 1 / 2
对于双原子分子，它的振动配分函数可以写作（近似后结果）：
q v ′ = 1 1 − e − h ν / k B T q_v'=\frac 1{1-e^{-h\nu/k_BT}}
q v ′ = 1 − e − h ν / k B T 1
q v ′ q_v' q v ′ ϵ v , 0 = 1 2 h ν \epsilon_{v,0}=\frac 12h\nu ϵ v , 0 = 2 1 h ν
对于多原子分子，需要考虑分子的自由度，描述一个分子质心在空间的位置需要三个坐标参数，描述分子骨架的空间取向需要三个角度参数（线性分子只需两个），考虑分子中的n n n 3 n 3n 3 n 3 n − 3 − 3 3n-3-3 3 n − 3 − 3 3 n − 3 − 2 3n-3-2 3 n − 3 − 2
有了自由度，就可以给出振动配分函数，对于非线形多原子分子：
q v = ∏ i = 1 3 n − 6 e − h ν i / 2 k B T 1 − e − h ν i / k B T q_v=\prod_{i=1}^{3n-6}\frac {e^{-h\nu_i/2k_BT}}{1-e^{-h\nu_i/k_BT}}
q v = i = 1 ∏ 3 n − 6 1 − e − h ν i / k B T e − h ν i / 2 k B T
经典理论认为，晶体中的原子在晶格上做简谐振动，按照能量均分定理，推出其摩尔热容为：
C V , m = ( ∂ U ∂ T ) V = 3 R = 24.94 J ⋅ m o l − 1 ⋅ K − 1 C_{V,m}=\left(\frac {\partial U}{\partial T} \right)_V=3R=24.94J\cdot mol^{-1}\cdot K^{-1}
C V , m = ( ∂ T ∂ U ) V = 3 R = 2 4 . 9 4 J ⋅ m o l − 1 ⋅ K − 1
这是一个经验定律，称为Dulong-Petit定律，即固体物质的摩尔热容大致相等。
但是实际情况不是这样，当温度降低，摩尔热容逐渐减小，且不同单质的摩尔热容随温度变化率也不相同。经典理论无法解释。
1907年，Einstein用量子论来解释晶体的热容，对这个问题给出了一个初步的解释。但是Einstein的结果与实验曲线不完全相合，这说明他的假设中存在些许问题。
Debye对Einstein理论提出了修正，并得到了令人满意的结果。当温度趋于零时：
C V , m = 1943 ( T Θ D ) 3 J ⋅ m o l − 1 ⋅ K − 1 = c o n s t a n t × T 3 \begin{aligned} C_{V,m} &=1943 \left(\frac {T}{\Theta_D} \right)^3 J\cdot mol^{-1} \cdot K^{-1} \\[3mm] &=constant \times T^3 \end{aligned}
C V , m = 1 9 4 3 ( Θ D T ) 3 J ⋅ m o l − 1 ⋅ K − 1 = c o n s t a n t × T 3
即T 3 T^3 T 3 Θ D \Theta_D Θ D
将上述讨论的几种配分函数合并，就得到分子的全配分函数：
q 总 = q n ⋅ q e ⋅ q t ⋅ q r ⋅ q v q_总=q_n \cdot q_e \cdot q_t \cdot q_r \cdot q_v
q 总 = q n ⋅ q e ⋅ q t ⋅ q r ⋅ q v
对于单原子分子：
q 总 = [ g n , 0 exp ( − ϵ n , 0 k B T ) ] ⋅ [ g e , 0 exp ( − ϵ e , 0 k B T ) ] ⋅ [ ( 2 π m k B T ) 3 / 2 h 3 V ] q_总= \left[g_{n,0}\exp \left(-\frac {\epsilon_{n,0}}{k_BT} \right)\right]\cdot \left[g_{e,0}\exp \left(-\frac {\epsilon_{e,0}}{k_BT}\right)\right] \cdot \left[\frac {(2\pi mk_BT)^{3/2}}{h^3}V \right]
q 总 = [ g n , 0 exp ( − k B T ϵ n , 0 ) ] ⋅ [ g e , 0 exp ( − k B T ϵ e , 0 ) ] ⋅ [ h 3 ( 2 π m k B T ) 3 / 2 V ]
对于双原子分子：
q 总 = [ g n , 0 exp ( − ϵ n , 0 k B T ) ] ⋅ [ g e , 0 exp ( − ϵ e , 0 k B T ) ] ⋅ [ ( 2 π m k B T ) 3 / 2 h 3 V ] ⋅ 8 π 2 I k B T σ h 2 ⋅ exp ( − 1 2 h ν k B T ) 1 − exp ( − h ν k B T ) q_总= \left[g_{n,0}\exp \left(-\frac {\epsilon_{n,0}}{k_BT} \right)\right]\cdot \left[g_{e,0}\exp \left(-\frac {\epsilon_{e,0}}{k_BT}\right)\right] \cdot \left[\frac {(2\pi mk_BT)^{3/2}}{h^3}V \right] \cdot \frac {8\pi^2Ik_BT}{\sigma h^2} \cdot \frac {\exp \left(-\frac 12\frac {h\nu}{k_BT} \right)}{1-\exp \left(-\frac {h\nu}{k_BT} \right)}
q 总 = [ g n , 0 exp ( − k B T ϵ n , 0 ) ] ⋅ [ g e , 0 exp ( − k B T ϵ e , 0 ) ] ⋅ [ h 3 ( 2 π m k B T ) 3 / 2 V ] ⋅ σ h 2 8 π 2 I k B T ⋅ 1 − exp ( − k B T h ν ) exp ( − 2 1 k B T h ν )
对于线形多原子分子：
q 总 = [ g n , 0 exp ( − ϵ n , 0 k B T ) ] ⋅ [ g e , 0 exp ( − ϵ e , 0 k B T ) ] ⋅ [ ( 2 π m k B T ) 3 / 2 h 3 V ] ⋅ 8 π 2 I k B T σ h 2 ⋅ ∏ i = 1 3 n − 5 exp ( − 1 2 h ν i k B T ) 1 − exp ( − h ν i k B T ) q_总= \left[g_{n,0}\exp \left(-\frac {\epsilon_{n,0}}{k_BT} \right)\right]\cdot \left[g_{e,0}\exp \left(-\frac {\epsilon_{e,0}}{k_BT}\right)\right] \cdot \left[\frac {(2\pi mk_BT)^{3/2}}{h^3}V \right] \cdot \frac {8\pi^2Ik_BT}{\sigma h^2} \cdot \prod_{i=1}^{3n-5} \frac {\exp \left(-\frac 12\frac {h\nu_i}{k_BT} \right)}{1-\exp \left(-\frac {h\nu_i}{k_BT} \right)}
q 总 = [ g n , 0 exp ( − k B T ϵ n , 0 ) ] ⋅ [ g e , 0 exp ( − k B T ϵ e , 0 ) ] ⋅ [ h 3 ( 2 π m k B T ) 3 / 2 V ] ⋅ σ h 2 8 π 2 I k B T ⋅ i = 1 ∏ 3 n − 5 1 − exp ( − k B T h ν i ) exp ( − 2 1 k B T h ν i )
对于非线形多原子分子：
q 总 = [ g n , 0 exp ( − ϵ n , 0 k B T ) ] ⋅ [ g e , 0 exp ( − ϵ e , 0 k B T ) ] ⋅ [ ( 2 π m k B T ) 3 / 2 h 3 V ] ⋅ [ 8 π 2 ( 2 π k B T ) 3 / 2 σ h 3 ( I x ⋅ I y ⋅ I z ) 1 / 2 ] ⋅ ∏ i = 1 3 n − 6 exp ( − 1 2 h ν i k B T ) 1 − exp ( − h ν i k B T ) q_总= \left[g_{n,0}\exp \left(-\frac {\epsilon_{n,0}}{k_BT} \right)\right]\cdot \left[g_{e,0}\exp \left(-\frac {\epsilon_{e,0}}{k_BT}\right)\right] \cdot \left[\frac {(2\pi mk_BT)^{3/2}}{h^3}V \right] \cdot \left[\frac {8\pi^2(2\pi k_BT)^{3/2}}{\sigma h^3}(I_x \cdot I_y \cdot I_z)^{1/2} \right] \cdot \prod_{i=1}^{3n-6} \frac {\exp \left(-\frac 12\frac {h\nu_i}{k_BT} \right)}{1-\exp \left(-\frac {h\nu_i}{k_BT} \right)}
q 总 = [ g n , 0 exp ( − k B T ϵ n , 0 ) ] ⋅ [ g e , 0 exp ( − k B T ϵ e , 0 ) ] ⋅ [ h 3 ( 2 π m k B T ) 3 / 2 V ] ⋅ [ σ h 3 8 π 2 ( 2 π k B T ) 3 / 2 ( I x ⋅ I y ⋅ I z ) 1 / 2 ] ⋅ i = 1 ∏ 3 n − 6 1 − exp ( − k B T h ν i ) exp ( − 2 1 k B T h ν i )
Δ r G m ⊖ \Delta_rG_m^{\ominus} Δ r G m ⊖ 在化学反应中，当几种物质共同存在时，不得不选取一个公共的能量零点才能衡量不同物质的能量。于是，按公共能量零点计算的各分子的能量：
ϵ = ϵ j + ϵ 0 \epsilon = \epsilon_j + \epsilon_0
ϵ = ϵ j + ϵ 0
式中ϵ \epsilon ϵ ϵ j \epsilon_j ϵ j ϵ 0 \epsilon_0 ϵ 0
有了公共能量标度，计算出的配分函数应为：
q ′ = ∑ j g j e − ( ϵ 0 + ϵ j ) / k B T = e − ϵ 0 / k B T ⋅ q q'=\sum_j g_j e^{-(\epsilon_0+\epsilon_j)/k_BT} = e^{-\epsilon_0/k_BT}\cdot q
q ′ = j ∑ g j e − ( ϵ 0 + ϵ j ) / k B T = e − ϵ 0 / k B T ⋅ q
然后就可以计算出相应的热力学函数（结果此处略去），发现S , C V , p S,C_V,p S , C V , p A , G , H , U A,G,H,U A , G , H , U U 0 U_0 U 0 U 0 = N ϵ 0 U_0=N\epsilon_0 U 0 = N ϵ 0 N N N
单独拿出G G G
G = A + p V = − N k B T ln q N + U 0 G=A+pV=-Nk_BT\ln\frac qN + U_0
G = A + p V = − N k B T ln N q + U 0
变形可得：
G ( T ) − U 0 T = − N k B ln q N \frac {G(T)-U_0}{T} = -Nk_B\ln\frac qN
T G ( T ) − U 0 = − N k B ln N q
上式等号左侧称为自由能函数（free energy function）。据此可见，根据配分函数就可以求出自由能函数，然后就可以进一步计算反应的平衡常数。
本章可以说是最难的一章，因为全是物理的内容，学校授课也是处于“老师讲不好，学生学不懂”的状态，所以凑合凑合看吧。
这篇拖了一个月…