Norskov, Fundamental Concepts in Heterogeneous Catalysis, Chapter 4. 슬라이더를 움직이면 그래프가 실시간으로 반응합니다.
Activation EnergyTransition State TheoryArrhenius EquationPre-exponential FactorSurface Kinetics
Section 4.1
Potential Energy Surfaces
화학반응은 반응좌표(reaction coordinate)를 따라 Initial State(IS)에서 Transition State(TS)를 거쳐 Final State(FS)로 진행합니다. 활성화 에너지 $E_a$는 IS와 TS의 에너지 차이이며, 반응이 진행되기 위해 넘어야 하는 에너지 장벽입니다.
반응 에너지 $\Delta E$는 IS와 FS의 에너지 차이로, 반응의 열역학적 구동력을 결정합니다. 역반응의 활성화 에너지는 $E_{a,\text{reverse}} = E_a - \Delta E$로 주어집니다.
$$E_a = E_{\text{TS}} - E_{\text{IS}}$$
$$\Delta E = E_{\text{FS}} - E_{\text{IS}}$$
$$E_{a,\text{reverse}} = E_a - \Delta E$$
$E_a$는 항상 양수입니다. 음수가 되면 barrier가 없다는 의미이므로 barrierless 반응이 됩니다.
$\Delta E < 0$이면 발열 반응(exothermic), $\Delta E > 0$이면 흡열 반응(endothermic)입니다.
촉매는 $E_a$를 낮추지만 $\Delta E$는 변하지 않습니다. 반응의 열역학은 촉매로 바꿀 수 없습니다.
$\exp(-E_a/k_BT)$ = 성공 확률. $E_a$가 클수록, 온도가 낮을수록 이 값은 급격히 감소합니다.
$E_a$가 0.2 eV 증가하면 300 K에서 $k$가 약 1000배 감소합니다. 촉매의 핵심 역할입니다.
Figure 4.2Rate vs Temperature
Ea0.80 eV
k(300K)--
k(500K)--
k(800K)--
Ea가 0.2 eV 변하면 rate이 ~1000배 변합니다 (300 K)
이것이 촉매가 반응속도를 극적으로 바꾸는 이유입니다. 활성화 에너지의 작은 변화가 지수함수를 통해 rate의 거대한 변화로 증폭됩니다. 온도가 높아질수록 이 효과는 완화됩니다.
Ea = 0.3 eV로 설정하고 k(300K) 값을 확인하세요.
Ea = 1.5 eV로 변경하고 k(300K) 값을 확인하세요.
두 값의 차이가 몇 자릿수(orders of magnitude)인지 계산하세요.
같은 비교를 800 K에서 해보세요. 온도가 높아지면 Ea의 영향이 줄어드는 것을 확인할 수 있습니다.
Key Insight: 높은 온도에서는 더 많은 분자가 barrier를 넘을 수 있는 에너지를 갖고 있으므로, Ea의 차이가 rate에 미치는 영향이 줄어듭니다. 이것이 "고온에서는 촉매의 역할이 줄어든다"는 의미입니다.
Section 4.3
Arrhenius Analysis
실험에서 rate constant를 여러 온도에서 측정한 뒤, $\ln k$ vs $1/T$ 그래프(Arrhenius plot)를 그리면 직선이 됩니다. 이 직선의 기울기에서 활성화 에너지를, y절편에서 pre-exponential factor $A$를 추출할 수 있습니다.
$$\ln k = \ln A - \frac{E_a}{k_B}\cdot\frac{1}{T}$$
유도 과정 보기
1
TST rate: $k = (k_BT/h)\exp(-E_a/k_BT)$
출발점은 Transition State Theory의 rate expression입니다.
2
양변에 $\ln$: $\ln k = \ln(k_BT/h) - E_a/(k_BT)$
로그를 취하면 지수 항이 선형 항으로 바뀝니다.
3
$\ln(k_BT/h)$의 $T$ 의존성은 $\ln T$로 매우 약합니다. 이를 상수 $\ln A \approx 29.5$로 근사합니다.
Pre-exponential factor $A \approx k_BT/h$는 온도에 선형적으로만 의존하므로, 지수 항에 비해 변화가 작습니다.
4
최종: $\ln k = \ln A - (E_a/k_B)\cdot(1/T)$. 기울기 = $-E_a/k_B$
$1/T$에 대한 직선 관계이며, 기울기의 절대값에 $k_B$를 곱하면 $E_a$를 얻습니다.