반응속도론

계산화학 분자동역학 마지막 수정: 2026-07-24 04:39:28

1. 개요[편집]

반응속도론 Chemical Kinetics
다루는 것	반응의 속도 · 차수 · 메커니즘
기본량	속도 상수 $k$, 반응 차수, 활성화 에너지
온도 의존성	아레니우스 방정식
이론적 기반	전이상태 이론
수치적 성격	강성(stiff) 연립 ODE

열역학은 “그 반응이 일어날 수 있나”를 묻는다. 반응속도론은 “그래서 대체 언제 끝나나”를 묻는다.

반응속도론(chemical kinetics)은 화학 반응이 얼마나 빠르게 진행되는지를 정량적으로 기술하고, 그 속도를 결정하는 미시적 경로(메커니즘)를 밝히는 분야다. 열역학이 반응의 방향과 최종 평형 상태만 말해 준다면, 반응속도론은 거기까지 가는 경로와 시간표를 다룬다. 다이아몬드가 흑연보다 열역학적으로 불안정한데도 반지에 박힌 채 멀쩡한 이유는 전적으로 반응속도론의 관할이다.¹

이 문서는 분야 전체 — 속도식, 반응 차수, 율속 단계, 메커니즘, 그리고 그것을 컴퓨터로 적분하는 문제 — 를 다룬다. 반면 속도 상수의 온도 의존성이라는 하나의 경험식은 아레니우스 방정식 문서가 따로 맡는다. 여기서는 그 식이 반응속도론이라는 큰 지도에서 정확히 어느 칸에 놓이는지만 짚고 넘어간다.

2. 속도식과 반응 차수[편집]

출발점은 속도식(rate law)이다. 반응물 A, B가 참여하는 반응의 속도는 보통 농도의 거듭제곱 곱으로 쓴다.

r = -\frac{1}{a}\frac{d[\mathrm{A}]}{dt} = k\,[\mathrm{A}]^{m}[\mathrm{B}]^{n}

$k$ 는 속도 상수, 지수 $m$ 과 $n$ 은 각 성분에 대한 반응 차수(order)이며 그 합이 전체 차수다. 여기서 입문자가 가장 자주 밟는 지뢰가 있다. 반응 차수는 화학반응식의 계수와 원칙적으로 아무 관계가 없다. 차수는 실험으로 결정되는 양이고, 계수는 종이 위의 물질수지일 뿐이다. 둘이 일치하는 경우는 반응이 진짜로 한 번의 충돌로 끝나는 소반응(elementary reaction)일 때뿐이다.

차수별 적분형 속도식은 실험 데이터를 직선으로 펴서 차수를 판별하는 도구로 쓰인다.

차수	적분형	직선이 되는 축	반감기
0차	$[\mathrm{A}] = [\mathrm{A}]_0 - kt$	농도 대 시간	초기 농도에 비례
1차	$\ln[\mathrm{A}] = \ln[\mathrm{A}]_0 - kt$	로그 농도 대 시간	초기 농도와 무관
2차	$1/[\mathrm{A}] = 1/[\mathrm{A}]_0 + kt$	역농도 대 시간	초기 농도에 반비례

1차 반응의 반감기가 초기 농도와 무관하다는 성질은 방사성 붕괴나 약물 대사에서 “반감기”라는 단일 숫자가 통용될 수 있는 근거다. 반대로 2차 반응에서 반감기를 한 숫자로 보고하는 문서를 보면 일단 의심하고 봐야 한다.

회분식 반응기의 연속반응 A→B→C 를 RK4로 적분한 결과다. 위쪽은 농도-시간 곡선으로, k2/k1 을 키우면 중간체 B의 봉우리가 앞으로 당겨지면서 낮아진다. 아래쪽은 같은 [A](t) 데이터를 ln[A] 대 t 와 1/[A] 대 t 두 축에 올리고 최소제곱 직선을 적합한 차수 판정 플롯이며, 차수 n 을 1에서 2로 바꾸면 R2 가 1이 되는 쪽이 왼쪽에서 오른쪽으로 넘어간다. 머리글의 QSSA 오차는 준정상상태 근사 B≈k1[A]^n/k2 의 최대 편차로, k2/k1 이 커질수록 줄어든다. 등온·정용 조건이다.

3. 메커니즘과 율속 단계[편집]

실제 반응은 대개 여러 소반응이 사슬처럼 이어진 메커니즘(mechanism)으로 진행된다. 관측된 전체 속도식은 이 사슬이 만들어낸 결과물이지, 처음부터 주어진 것이 아니다. 메커니즘에서 속도식을 뽑아내는 표준 도구는 둘이다.

율속 단계 근사(rate-determining step) — 사슬 중 압도적으로 느린 한 단계가 전체 속도를 결정한다고 본다. 왕복 8차선 고속도로라도 1차선 공사 구간 하나가 통행량을 정하는 것과 같은 논리다.
정상상태 근사(steady-state approximation) — 반응성이 큰 중간체는 생기는 즉시 소모되어 농도가 거의 변하지 않는다고 두고( $d[\mathrm{I}]/dt \approx 0$ ) 중간체를 대수적으로 소거한다.

정상상태 근사의 가장 유명한 산물이 효소 반응의 미하엘리스-멘텐 식이다.

r = \frac{V_{\max}\,[\mathrm{S}]}{K_M + [\mathrm{S}]}

기질 농도가 낮을 때는 1차처럼, 충분히 높을 때는 0차처럼 보인다. 겉보기 차수는 조건에 따라 변한다는 반응속도론의 교훈이 한 식에 압축돼 있는 셈이다. 촉매는 이 지도 위에서 활성화 장벽이 더 낮은 우회로를 새로 뚫는 존재이며, 평형 위치는 전혀 건드리지 못하고 도달 시간만 줄인다. 촉매가 열역학을 이길 수 없다는 말이 이 뜻이다.

4. 온도 의존성 — 어디에 아레니우스가 들어가는가[편집]

농도 의존성이 차수로 정리되고 나면 속도식에 남는 미지수는 속도 상수 $k$ 하나다. 이 $k$ 가 온도에 어떻게 매달리는지를 기술하는 경험식이 아레니우스 방정식이고, 같은 내용을 분자 수준에서 다시 세운 것이 전이상태 이론이다. 전이상태 이론은 반응계가 포텐셜 에너지 표면 위의 안장점을 넘어가는 빈도로 $k$ 를 계산하며, 그 안장점을 실제로 찾아내는 계산 절차가 넛지 탄성 밴드법이다.

즉 현대적 반응속도론은 2층 구조다. 위층에서는 실험으로 차수와 $k$ 를 재고, 아래층에서는 밀도범함수이론으로 장벽 높이를 계산해 $k$ 를 예측한다. 두 층의 숫자가 맞아떨어지면 메커니즘 가설이 살아남고, 어긋나면 둘 중 하나가 틀린 것이다.² 분자동역학은 그 사이를 잇는 도구로, 희귀 사건인 반응을 그냥 기다리는 대신 자유에너지 섭동이나 편향 표본추출로 장벽 통과 빈도를 뽑아낸다.

5. 수치적으로 푼다는 것[편집]

여러 반응이 동시에 도는 계의 시간 변화는 결국 연립 상미분방정식이다.

\frac{d\mathbf{c}}{dt} = \boldsymbol{\nu}^{\mathsf{T}}\,\mathbf{r}(\mathbf{c}, T)

$\boldsymbol{\nu}$ 는 화학량론 행렬, $\mathbf{r}$ 은 각 반응의 속도 벡터다. 식은 얌전해 보이지만 실제로 적분하면 악명 높은 강성(stiff) 문제가 된다. 나노초 만에 사라지는 라디칼과 몇 초에 걸쳐 진행되는 총괄 반응이 같은 벡터 안에 섞여 있어, 시간 척도가 열 자릿수 넘게 벌어지기 때문이다. 명시적 룽게-쿠타법으로 밀어붙이면 안정성 조건이 가장 빠른 성분에 묶여 시간 간격이 실질적으로 0으로 내려간다. 그래서 실무 표준은 암시적 BDF 계열 솔버이고, 매 스텝 자코비안 행렬을 만들어 뉴턴-랩슨법으로 푼다.³

유동과 결합하면 규모가 한 단계 더 커진다. 연소 시뮬레이션은 격자 하나하나마다 수십~수백 화학종의 강성 ODE를 풀어야 해서, 화학 항과 수송 항을 번갈아 푸는 연산자 분리, 메커니즘 축약, 화염면 테이블화(flamelet) 같은 우회로가 총동원된다. 반대로 분자 개수가 적어 통계 요동을 무시할 수 없는 계에서는 결정론적 ODE 자체가 부적절해지고, 반응 사건 하나하나를 몬테카를로 방법식으로 추첨하는 동역학 몬테카를로(깁스피 알고리즘)로 넘어간다.

6. 반응과 확산이 겨룰 때[편집]

반응이 아무리 빨라도 반응물이 서로 만나지 못하면 소용없다. 반응 속도와 물질 전달 속도가 비슷해지는 순간, 계의 거동은 국소 화학이 아니라 반응-확산 방정식이 지배한다.

\frac{\partial c_i}{\partial t} = D_i \nabla^2 c_i + R_i(\mathbf{c})

둘의 비를 재는 무차원수가 담쾰러 수(Damköhler number)다.⁴ $\mathrm{Da} \ll 1$ 이면 반응이 율속이라 계 전체가 균일한 농도로 천천히 반응하고, $\mathrm{Da} \gg 1$ 이면 확산이 율속이라 반응이 얇은 층에서만 일어난다. 촉매 입자의 유효 인자, 화염면 두께, 다공성 담체 설계가 전부 이 경쟁의 결과다.

이 경쟁이 극단으로 가면 비선형 반응 항 하나 때문에 균일한 초기 상태가 스스로 무늬로 갈라지는 튜링 불안정까지 나타난다. 반응 속도 상수를 조금 바꿨을 뿐인데 정성적으로 다른 정상 상태가 나온다는 점이 핵심이다. 실제 계에서 그 상수의 손잡이를 쥐고 있는 것이 온도이며, 그 예민함을 정량화한 것이 다시 아레니우스 방정식이다. 반응속도론이 한 바퀴 도는 지점이다.

7. 관련 문서[편집]

흑연으로의 전환은 열역학적으로 확실히 유리하지만 활성화 장벽이 워낙 높아 상온에서의 속도가 사실상 0이다. “다이아몬드는 영원하다”는 광고 문구는 열역학이 아니라 반응속도론이 보증하는 셈이다. ↩
실제로는 “둘 다 틀렸다”인 경우가 가장 흔하다. 계산 쪽은 범함수 선택으로 장벽이 수 kcal/mol씩 흔들리고, 실험 쪽은 미량 불순물이 촉매 노릇을 하고 있다. 두 오차가 우연히 상쇄돼 논문이 나오는 일도 없지는 않다. ↩
화학 반응 ODE를 처음 명시적 적분기로 돌려 본 사람은 대개 “발산했다”고 생각한다. 사실은 발산이 아니라 강성이며, 솔버를 바꾸면 같은 문제가 몇 초 만에 풀린다. 강성은 문제의 성질이지 코드의 버그가 아니다. ↩
담쾰러 수는 이름이 여러 개다. 어떤 교재는 반응/대류 비를 Da로, 반응/확산 비를 Da II로 부른다. 논문에서 Da라는 기호를 봤다면 정의부터 확인하는 것이 국룰이다. ↩