벨루소프-자보틴스키 반응

계산화학 시뮬레이션 물리 마지막 수정: 2026-07-27 05:08:52

1. 개요[편집]

벨루소프-자보틴스키 반응(Belousov-Zhabotinsky reaction, BZ 반응)은 브로민산 이온이 말론산을 산화시키는 과정에서 촉매의 산화 상태가 주기적으로 뒤집히며, 용액의 색이 몇십 초 주기로 반복해서 변하는 대표적인 화학 진동자(chemical oscillator)다. 잘 저어준 비커에서는 시계처럼 규칙적인 색 진동이, 얇게 깔아놓은 페트리 접시에서는 동심원 표적 패턴과 나선파가 저절로 생겨난다. 평형에서 멀리 떨어진 열린 계가 어떻게 시간적·공간적 질서를 스스로 만들어내는지를 눈으로 보여주는, 자기조직화 분야의 국민 실험이다.

핵심은 되먹임이다. 어떤 중간체가 자기 자신의 생성을 가속하고(활성인자), 다른 중간체가 그것을 억누르되 한 박자 늦게 작동한다(억제인자). 이 구조가 분기 이론에서 말하는 호프 분기를 거쳐 극한 순환을 만든다. 여기에 확산이 붙으면 진행파가 되고, 파면이 끊어지면 나선이 감긴다.¹

2. 벨루소프의 거절당한 논문[편집]

1951년경, 소련의 생물물리학자 보리스 벨루소프는 크렙스 회로(구연산 회로)의 무기화학적 모형을 만들려고 브로민산칼륨·구연산·세륨염을 황산에 섞었다. 그런데 용액이 노란색(Ce $^{4+}$ )과 무색(Ce $^{3+}$ )을 왔다 갔다 하기 시작했다. 그는 이걸 논문으로 썼고, 두 번 거절당했다.

거절 사유가 전설적이다. 심사자는 “가역적이지 않은 반응이 평형으로 단조롭게 가지 않고 진동하는 것은 열역학 제2법칙 위반”이라며 반려했다. 물론 이건 오해다. 제2법칙이 금지하는 것은 닫힌 계가 평형에 도달한 뒤 다시 벗어나는 것이지, 평형에서 멀리 떨어진 계가 자유에너지를 태우며 진동하는 것이 아니다. BZ 반응도 결국 브로민산과 말론산을 다 써버리면 진동이 멈추고 조용히 평형으로 간다. 진동은 “연료를 태우는 동안”만 지속된다.²

벨루소프는 결국 1959년 방사선의학 학술대회 초록집이라는 아무도 안 보는 지면에 짧게 실은 뒤 연구를 접었다. 몇 년 뒤 시몬 시놀의 대학원생이던 아나톨 자보틴스키가 이 실험을 재현하고, 구연산을 말론산으로, 세륨을 페로인(ferroin, $\mathrm{Fe(phen)}_3^{2+/3+}$ )으로 바꿔 붉은색↔파란색의 극적인 색 대비를 만들어냈다. 이 조합이 지금 유튜브에 널린 그 영상이다. 벨루소프는 1970년에 사망했고, 1980년에 자보틴스키 등과 함께 레닌상을 추서받았다.

3. 화학 — FKN 메커니즘[편집]

전체 반응은 브로민산 이온이 촉매의 도움으로 말론산을 산화시켜 이산화탄소와 브로민화 유도체를 만드는 것이다.

3\,\mathrm{BrO_3^-} + 5\,\mathrm{CH_2(COOH)_2} + 3\mathrm{H^+} \;\longrightarrow\; 3\,\mathrm{BrCH(COOH)_2} + 4\,\mathrm{CO_2} + 2\,\mathrm{HCOOH} + 5\,\mathrm{H_2O}

이게 왜 진동하는지는 1972년 필드·쾨뢰시·노이에스가 제시한 FKN 메커니즘이 설명한다. 열여덟 개 남짓한 소반응을 세 뭉치로 묶은 것이 요령이다.

과정 A (브로민화물 소모): $\mathrm{Br^-}$ 농도가 높을 때 지배적. 브로민산이 브로민화물과 만나 산소 원자 이동 경로로 천천히 진행한다. 이 과정은 브로민화물을 갉아먹는다.
과정 B (자기촉매): $\mathrm{Br^-}$ 가 어떤 문턱 아래로 떨어지는 순간 스위치가 넘어간다. $\mathrm{BrO_3^-} + \mathrm{HBrO_2} + \mathrm{H^+} \to 2\,\mathrm{BrO_2^\bullet} + \mathrm{H_2O}$ 로 라디칼이 만들어지고, 이것이 촉매를 산화시키며 $\mathrm{HBrO_2}$ 를 두 배로 되돌려준다. 아브로민산 $\mathrm{HBrO_2}$ 의 폭발적 자기증식. 색이 확 바뀌는 순간이 여기다.
과정 C (재생): 산화된 촉매( $\mathrm{Ce^{4+}}$ 또는 페리인)가 말론산과 그 브로민화물을 산화시키면서 $\mathrm{Br^-}$ 를 다시 뱉어낸다. 브로민화물 농도가 올라가면 과정 B가 꺼지고, 다시 과정 A로 돌아간다.

정리하면 $\mathrm{HBrO_2}$ 는 활성인자(자기촉매), $\mathrm{Br^-}$ 는 억제인자이며, 억제인자의 재생이 산화된 촉매를 경유하기 때문에 반드시 시간 지연이 생긴다. 자기촉매 + 지연된 음의 되먹임 = 진동. 생물의 대사·유전자 회로가 진동하는 원리와 정확히 같은 골격이다.

4. 오레고네이터와 이완 진동[편집]

오리건 대학에서 나왔다고 해서 이름이 오레고네이터(Oregonator)인 이 5단계 축약 모형이 BZ 수치 시뮬레이션의 표준이다. $X=\mathrm{HBrO_2}$ , $Y=\mathrm{Br^-}$ , $Z$ = 산화된 촉매로 두고 무차원화하면 다음 3변수계가 된다.

\varepsilon\frac{dx}{d\tau} = qy - xy + x(1-x)

\delta\frac{dy}{d\tau} = -qy - xy + f z

\frac{dz}{d\tau} = x - z

전형적인 파라미터가 $\varepsilon \approx 4\times10^{-2}$ , $\delta \approx 4\times10^{-4}$ , $q \approx 8\times10^{-4}$ 다. 화학량론 인자 $f$ 는 촉매 한 분자가 환원될 때 몇 개의 브로민화물이 나오는지를 뭉뚱그린 조절 손잡이인데, 대략 $0.5 < f < 1+\sqrt{2} \approx 2.414$ 구간에서 정상상태가 불안정해지고 극한 순환이 생긴다. 그 양쪽 경계가 호프 분기점이다.

$\varepsilon$ 과 $\delta$ 가 $1$ 보다 몇 자릿수 작다는 사실이 이 계의 성격을 결정한다. $x$ 는 아주 빠르게, $z$ 는 느리게 움직이는 느림-빠름(slow-fast) 계이고, 그 결과 궤도는 사인파가 아니라 톱니에 가까운 이완 진동(relaxation oscillation)이 된다. 오래 참았다가 갑자기 튀고, 다시 천천히 회복한다. 심장 박동, 뉴런 발화, 릴랙세이션 발진 회로가 전부 같은 모양이다.

$\delta$ 가 극히 작으므로 $y$ 를 준정상상태로 두고 소거해 2변수 모형으로 줄이는 것이 흔한 관행이며, 이렇게 하면 위상평면 위에서 널클라인 두 개로 흥분성/진동성 여부를 손으로 판정할 수 있다.

5. 확산이 붙으면 — 표적, 나선, 그리고 심장[편집]

교반을 멈추고 얇은 층으로 깔면 각 화학종이 확산으로 이웃과 대화한다. 반응 항에 확산 항을 더한 편미분방정식이 된다.

\frac{\partial \mathbf{c}}{\partial t} = \mathbf{R}(\mathbf{c}) + \mathbf{D}\nabla^2 \mathbf{c}

이때 나타나는 구조들.

표적 패턴(target pattern): 먼지 입자나 기포 같은 국소 불균일이 다른 곳보다 빨리 발화하는 박동원(pacemaker)이 되어, 동심원 파가 주기적으로 뿜어져 나온다. 주기가 짧은 박동원이 이웃을 잡아먹는다.
나선파(spiral wave): 파면을 유리막대로 살짝 끊어주면 끊긴 끝(자유단)이 안쪽으로 말려 들어가며 나선이 감긴다. 한번 감기면 외부 박동원 없이 스스로 회전하며 영원히 파를 뿜는다. 3차원으로 확장하면 나선의 회전 중심이 선(필라멘트)이 되고, 이를 스크롤 웨이브(scroll wave)라 한다.
흥분성 매질의 성질: BZ 파는 선형 파동이 아니다. 문턱 이하의 자극은 아무 일도 일으키지 않고(전부 아니면 전무), 지나간 자리는 불응기(refractory period) 동안 재발화하지 않으며, 두 파면이 부딪히면 간섭하는 게 아니라 서로 소멸한다. 반사도 회절도 없다. 대류-확산 방정식의 직관을 그대로 가져오면 틀린다.

이 그림이 그대로 심장으로 넘어간다. 심근도 흥분성 매질이고, 정상 심장은 동방결절이라는 단일 박동원이 만드는 표적 패턴으로 동작한다. 그런데 국소 전도 지연이나 흉터 조직 때문에 파면이 끊기면 재진입 회로(reentry)가 생기고, 나선파가 자리를 잡으면 동방결절보다 빠르게 회전하며 심장을 장악한다 — 이것이 심실빈맥이다. 나선이 불안정해져 여러 개로 부서지면 심실세동이 되고, 제세동기는 강한 전기 충격으로 매질 전체를 동시에 불응 상태로 만들어 나선을 지워버리는 장치다. 아서 윈프리가 BZ 접시와 심장 부정맥을 같은 언어로 묶은 이 통찰이, 화학 실험을 임상 심장학의 계산 모델로 승격시켰다.³ 반대 방향으로, 파면의 충돌·소멸 자체를 계산 자원으로 쓰려는 시도도 있다.⁴

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

ODE 단계는 강성(stiff) 문제다. $\varepsilon$ 과 $q$ 가 $10^{-4}$ 규모라 시간 척도가 네 자릿수 이상 벌어진다. 명시적 룽게-쿠타법으로 풀면 안정성 조건 때문에 시간 간격이 말도 안 되게 작아진다. BDF나 로젠브록 계열의 암시적 적분기(CVODE/LSODE 부류)를 쓰는 게 정석이다.
PDE는 연산자 분할이 국룰이다. 한 시간 스텝을 확산 반 스텝 → 반응 한 스텝 → 확산 반 스텝으로 쪼개는 스트랑 분할을 쓰면 2차 정확도를 유지하면서, 확산은 크랭크-니콜슨법 같은 암시적 도식으로, 반응은 국소 강성 적분기로 각각 최적화할 수 있다. 반응 항이 격자점끼리 결합돼 있지 않다는 점이 이 분할을 아주 싸게 만든다.
파면을 해상하라. 흥분성 파면은 두께가 수십 마이크로미터 수준으로 얇다. 격자가 파면 두께를 못 담으면 파속이 격자 간격에 따라 변하고, 심하면 수치적으로 파가 아예 전파되지 못하는(propagation failure) 가짜 현상이 생긴다. 격자 수렴 확인이 선택이 아니라 필수다.
튜링 패턴과 혼동하지 말 것. 튜링 불안정성은 균질 정상상태가 안정한데 확산 때문에 불안정해져 정지한 주기 구조가 생기는 현상이고, 활성인자와 억제인자의 확산계수 차이가 결정적이다. BZ의 표적·나선파는 국소 동역학이 이미 진동성 또는 흥분성이고 확산은 그 위상을 전달하는 역할이며, 패턴이 움직인다. 둘 다 반응-확산계지만 원인이 다르다.

7. 관련 문서[편집]

BZ 용액은 지금도 시약만 사면 부엌에서 재현 가능하다. 다만 브로민산염과 진한 황산을 다루는 실험이므로 후드와 보안경 없이 “유튜브 보고 따라했다”를 시도하면 안 된다. 예쁜 색을 보려다 브로민 증기를 마시는 건 국룰이 아니라 사고다. ↩
사실 1910년대에 로트카가 진동 반응의 수학적 모형을 이미 내놓았고, 1920년대에는 브레이 반응(과산화수소-아이오딘산)이 보고돼 있었다. 그런데도 화학계는 40년 넘게 “균질 용액은 진동할 수 없다”는 통념을 고수했다. 벨루소프의 심사자만 유별났던 게 아니라 그게 그 시절의 표준 상식이었다는 점이 더 무섭다. ↩
윈프리는 나선파의 회전 중심(위상 특이점)을 “위상이 정의되지 않는 점”으로 다뤘는데, 이건 수학적으로 위상수학적 결함(topological defect)이고 초전도체의 소용돌이나 액정의 디스클리네이션과 같은 부류다. 화학 접시와 심장과 초전도체가 같은 방정식의 다른 얼굴이라는 게 이 분야의 매력이다. ↩
BZ 반응을 계산 소자로 쓰려는 시도도 진지하게 존재한다. 파면의 충돌·소멸을 논리 게이트로 삼는 “화학 컴퓨터” 연구인데, 클럭이 초 단위라 실용성은 없다. 대신 미로 최단경로 찾기 같은 문제는 파면이 알아서 병렬로 풀어준다 — 물리가 곧 알고리즘인 사례. ↩