움직이는 교과서 · Interactive Textbook

통계 열역학

Atkins, Physical Chemistry, Chapter 13. 분배함수에서 열역학 성질이 어떻게 나오는지 시각적으로 탐색합니다.

Partition Function Boltzmann Distribution Heat Capacity Einstein Model

분배함수와 Boltzmann 분포

분배함수(partition function) $q$는 통계역학의 중심 개념으로, 모든 가능한 에너지 상태에 대한 Boltzmann 인자의 합이다: $q = \sum_i g_i \exp(-\varepsilon_i / k_BT)$. 이 하나의 양으로부터 내부 에너지, 엔트로피, 자유 에너지 등 모든 열역학 함수를 유도할 수 있다. 각 에너지 준위 $\varepsilon_i$에 있을 확률은 $p_i = g_i \exp(-\varepsilon_i / k_BT) / q$로, 이것이 Boltzmann 분포이다.

$$q = \sum_i g_i \exp\!\left(-\frac{\varepsilon_i}{k_BT}\right)$$
$$p_i = \frac{g_i \exp(-\varepsilon_i / k_BT)}{q}$$
$$\langle E \rangle = -\frac{\partial \ln q}{\partial \beta}, \quad \beta = \frac{1}{k_BT}$$

진동 분배함수와 에너지 준위 점유

조화 진동자의 분배함수 $q_\text{vib} = 1/(1 - e^{-\Theta_\text{vib}/T})$에서 $\Theta_\text{vib} = hc\tilde{\nu}/k_B$는 특성 진동 온도이다. 에너지 간격 $\tilde{\nu}$가 클수록(강한 결합, 가벼운 원자) $\Theta_\text{vib}$가 높아져 실온에서 $q_\text{vib} \approx 1$이 되고, 거의 모든 분자가 바닥 상태에 남는다. 반대로 낮은 $\tilde{\nu}$에서는 여러 준위가 점유되어 $q_\text{vib}$가 커진다.

$$q_\text{vib} = \frac{1}{1 - e^{-\Theta_\text{vib}/T}}, \quad \Theta_\text{vib} = \frac{hc\tilde{\nu}}{k_B}$$
Interactive Figure 13.1 분배함수 $q$ vs $T$ 및 준위 점유 분포
$\tilde{\nu}$ (cm-1) 600 cm⁻¹
T (K) 300 K
Θvib 863 K
qvib 1.06
p(v=0) 94.3%

분배함수와 열용량

진동 모드의 열용량 기여를 Einstein 모델로 계산하면 $C_{V,\text{vib}}/R = (\Theta_\text{vib}/T)^2 \cdot e^{\Theta_\text{vib}/T}/(e^{\Theta_\text{vib}/T} - 1)^2$이다. 낮은 온도($T \ll \Theta_\text{vib}$)에서 진동 모드는 "동결"되어 열용량에 기여하지 못하고, 높은 온도($T \gg \Theta_\text{vib}$)에서 등분배 정리에 따라 모드당 $R$에 수렴한다. 이것이 기체의 열용량이 온도에 따라 계단식으로 증가하는 이유이다: 병진(3/2 R), 회전(R 또는 3/2 R), 진동(R per mode) 순서로 각 자유도가 "해동"된다.

$$\frac{C_{V,\text{vib}}}{R} = \left(\frac{\Theta_\text{vib}}{T}\right)^2 \frac{e^{\Theta_\text{vib}/T}}{(e^{\Theta_\text{vib}/T} - 1)^2}$$
$$C_V(\text{total}) = \tfrac{3}{2}R\;(\text{trans}) + R_\text{rot} + \sum_{\text{modes}} C_{V,\text{vib}}$$
Interactive Figure 13.2 Einstein 모델: 열용량 vs 온도
Θvib (K) 3340 K
단원자 ($\frac{3}{2}R$)
이원자 (병진+회전+진동)
진동 기여만 ($C_{V,\text{vib}}/R$)
Θvib 3340 K
예시 H₂

열역학 함수의 통계역학적 표현

분배함수로부터 모든 열역학 함수를 유도할 수 있다. 내부 에너지는 $U - U(0) = -(\partial \ln Q / \partial \beta)_V$이고, 엔트로피는 $S = k_B \ln Q + U/T$이며, Helmholtz 자유 에너지는 $A = -k_BT \ln Q$이다. 여기서 $Q$는 정준 분배함수(canonical partition function)로, 구별 불가능한 입자의 경우 $Q = q^N / N!$이다.

$$U - U(0) = -\left(\frac{\partial \ln Q}{\partial \beta}\right)_V$$
$$S = \frac{U - U(0)}{T} + k_B \ln Q$$
$$A = -k_BT \ln Q, \quad Q = \frac{q^N}{N!}$$
$\Theta_\text{vib} = 3340\;\text{K}$인 H$_2$에서 실온(300 K)의 진동 분배함수 $q_\text{vib}$는 약 얼마인가?
$\Theta_\text{vib}/T = 3340/300 \approx 11.1$이므로 $e^{-11.1} \approx 1.5 \times 10^{-5}$. 따라서 $q_\text{vib} = 1/(1 - e^{-11.1}) \approx 1.00$. 실온에서 H$_2$는 진동 바닥 상태에만 있다.
이원자 기체의 열용량이 $C_V = \frac{5}{2}R$인 온도 영역에서 어떤 운동 자유도가 활성인가?
$\frac{5}{2}R = \frac{3}{2}R\,(\text{trans}) + R\,(\text{rot})$. 이원자 분자의 2개 회전 자유도가 각각 $\frac{1}{2}R$을 기여하여 총 $R$을 더한다. 진동 모드는 아직 "동결"된 상태이다.