흥분성 매질

1. 개요[편집]

건드리기 전엔 아무 일도 없다. 문제는 조금 세게 건드렸을 때다.

흥분성 매질(excitable medium)은 정지 상태가 안정하지만, 역치(threshold)를 넘는 자극을 받으면 자극 크기와 거의 무관하게 크고 정형화된 반응을 일으킨 뒤, 불응기(refractory period)를 거쳐 다시 정지 상태로 돌아오는 확산 결합계다. 신경 축삭, 심근, 벨루소프-자보틴스키 반응 용액, 세포 내 칼슘 방출계, 심지어 산불이 번지는 숲까지 전부 이 범주에 든다.

핵심 정체성은 역치와 불응기 두 개다. 역치 아래 자극은 지수적으로 스러지고(선형 감쇠), 역치 위 자극은 진폭이 정해진 펄스로 증폭되어 감쇠 없이 전파된다. 그리고 지나간 자리는 한동안 다시 흥분할 수 없다 — 이 두 성질이 결합하면 파동이 되돌아오지 못하고, 충돌하면 서로를 지워버리며, 파면이 끊기면 나선으로 감긴다.

2. 피츠휴-나구모 모형[편집]

가장 널리 쓰이는 최소 모형은 피츠휴-나구모(FitzHugh–Nagumo) 방정식이다. 빠른 활성화 변수 $u$(막전위 격)와 느린 회복 변수 $v$(불응 게이트 격) 두 개면 흥분성의 본질이 전부 재현된다.

$$\epsilon \dot{u} = u - \frac{u^3}{3} - v + I, \qquad \dot{v} = u + a - b v$$

여기서 $0 < \epsilon \ll 1$이 시간척도 분리 계수다. $u$는 $v$보다 $1/\epsilon$배 빠르게 움직이므로, 위상평면에서 궤도는 거의 수평으로 튀었다가 느린 곡선을 따라 기어가는 이완진동(relaxation oscillation) 형태가 된다.

기하는 널클라인으로 읽는다.

• $u$-널클라인: $v = u - u^3/3 + I$ — 뒤집힌 N자 3차 곡선. 왼쪽·오른쪽 가지는 끌개(각각 정지 상태, 흥분 상태), 가운데 가지는 척력이라 이 가운데 가지가 곧 역치선이다.
• $v$-널클라인: $v = (u+a)/b$ — 직선.

두 널클라인의 교점이 유일하고 그것이 왼쪽 가지 위에 있으면 정지 상태가 안정한 흥분성이다. 파라미터를 밀어 교점이 가운데 가지로 올라가면 호프 분기가 일어나 정지 상태가 불안정해지고 계는 스스로 진동하는 극한 순환 모드(진동성 매질)로 넘어간다. 즉 흥분성과 진동성은 분기 이론적으로 이웃 사이다.

여기에 확산을 붙이면($\epsilon \dot u = \dots + D\nabla^2 u$) 진행파가 생긴다. 전파 속도는 대략 $c \sim \sqrt{D/\tau}$로 스케일하며($\tau$는 활성화 시간척도), 확산계수를 4배로 올리면 속도는 2배가 된다.

3. 호지킨-헉슬리에서 바클리까지[편집]

계보를 짚으면 이렇다.

• 호지킨-헉슬리(1952). 오징어 거대축삭의 막전위를 $V, m, h, n$ 4변수 상미분방정식으로 기술. 1963년 노벨 생리의학상. 정량적으로 정확하지만 변수 넷에 지수함수 게이트가 잔뜩 붙어 있어 위상평면 직관이 안 나온다.
• 피츠휴(1961)·나구모(1962). 빠른 $m$을 준정상 근사로 소거하고 $h+n \approx$ 상수라는 관측을 써서 2변수로 축약했다. 나구모는 같은 방정식을 터널 다이오드 회로로 실제 구현해 보였다.
• 바클리 모형(1991). 수치 효율에 특화된 모형.

$$\frac{\partial u}{\partial t} = \frac{1}{\epsilon} u (1-u)\left(u - \frac{v+b}{a}\right) + \nabla^2 u, \qquad \frac{\partial v}{\partial t} = u - v$$

바클리 모형이 빠른 이유는 구조에 있다. $u$가 $[0,1]$ 안에 갇히고 넓은 영역에서 정확히 $u=0$이라, 그 영역에서는 $v$가 $v \leftarrow v e^{-\Delta t}$로 해석적으로 갱신되므로 반응항 계산을 통째로 건너뛸 수 있다. 나선파 시뮬레이션에서 활성 영역은 전체의 몇 %뿐이므로 이 최적화만으로 한 자릿수 이상 빨라진다.

4. 진행파, 충돌 소멸, 전도 블록[편집]

흥분성 매질의 파동은 선형 파동이 아니다. 이 차이가 실무에서 가장 중요하다.

• 감쇠하지 않는다. 에너지를 매질에서 뽑아 쓰므로 진폭이 거리에 따라 줄지 않는다. 축삭 1 m 끝까지 활동전위 크기가 유지되는 이유.
• 중첩되지 않고 충돌하면 소멸한다(annihilation). 마주 오는 두 펄스는 서로를 통과하지 못한다. 각 파면 뒤가 불응 상태라 상대의 진행 방향에 흥분 가능한 매질이 없기 때문이다. 선형 파동 방정식의 중첩·간섭과 정면으로 대비되는 지점.
• 반사되지 않는다. 경계에 부딪히면 그냥 사라진다.
• 되돌아오지 못한다. 불응 꼬리 덕분에 파면은 항상 앞으로만 간다.

불응기는 둘로 나뉜다. 절대 불응기 동안에는 자극 크기와 무관하게 아무 반응도 못 하고(회복 변수가 아직 높아 3차 곡선 왼쪽 가지에 갇힘), 상대 불응기에는 평소보다 큰 자극이면 반응하되 진폭과 전도 속도가 낮아진다.

여기서 전도 블록이 나온다. 자극 주기를 불응기보다 짧게 하면 매질이 모든 자극을 따라가지 못하고 2:1 블록(두 번에 한 번만 통과)으로 떨어진다. 주기를 더 줄이면 3:1, 웬케바흐형 부정 주기 등 사다리꼴로 이어진다. 심전도에서 보이는 방실 블록이 바로 이 현상이고, 수치 실험에서도 자극 주기를 쓸어보면 같은 계단이 그대로 재현된다.

5. 나선파, 스크롤파, 그리고 부정맥[편집]

파면이 어딘가에서 끊기면(불응 영역이나 장애물을 만나 파단이 생기면) 자유단이 감기면서 나선파(spiral wave)가 만들어진다. 나선은 외부 자극 없이 스스로 회전하며 매질을 계속 재점화하는 자립 파원이 되고, 대개 원래의 정상 박동보다 빠르다. 3차원에서는 나선의 회전 중심선이 필라멘트가 되고 파면은 스크롤파가 된다.

• 심장. 정상 심장은 동방결절이라는 하나의 페이스메이커가 지휘하는 흥분성 매질이다. 여기에 나선파가 자리 잡으면 심실빈맥, 나선이 불안정해져 여러 조각으로 깨지면(spiral breakup) 무질서한 다중 파원 상태 — 즉 심실세동이 된다. 세동은 “심장이 랜덤하게 뛰는” 게 아니라 결정론적 반응-확산계의 시공간 카오스에 가깝다는 것이 계산 심장학의 표준 관점이다.¹ 제세동은 매질 전체를 동시에 불응 상태로 만들어 모든 나선을 한 번에 지우는 조작이다.
• BZ 반응. 얇은 페트리 접시의 벨루소프-자보틴스키 반응 용액에서 동심원 표적 패턴과 나선이 육안으로 보인다. 흥분성 매질 연구의 실험용 표준 시료.
• 칼슘 파. 난모세포·심근세포 안에서 IP3 유도 칼슘 방출이 흥분성 동역학을 이루어 세포 내부를 가로지르는 파를 만든다.
• 점균·산불. 세포성 점균의 cAMP 신호, 그리고 셀룰러 오토마타로 흔히 모형화하는 산불 확산도 역치+불응 구조를 공유한다.

6. 튜링 패턴과 무엇이 다른가[편집]

튜링 패턴과 흥분성 매질은 둘 다 반응-확산 방정식이라 자주 뭉뚱그려지는데, 출발점이 정반대다.

항목	튜링 패턴	흥분성 매질
균일 정지해	불안정(확산 유도 불안정)	안정
발동 조건	조건 없음, 미소 요동이 저절로 자람	역치를 넘는 유한 자극 필요
확산비	억제자 확산이 훨씬 빨라야 함	조건 없음(단일 확산도 가능)
결과	정지한 정상 패턴(줄무늬·점)	진행파, 나선, 스크롤

한 줄로 요약하면 튜링은 가만히 둬도 저절로 무늬가 생기는 계이고, 흥분성 매질은 가만히 두면 아무 일도 없지만 한 번 건드리면 크게 반응하는 계다.

7. 수치해법에서의 함정[편집]

반응-확산계를 명시적으로 푸는 초보자가 반드시 밟는 지뢰가 있다.

• 두 개의 상한이 따로 논다. 명시적 오일러에는 제약이 둘 붙는다. 확산 제한 $\Delta t \le h^2/(4D)$는 격자 간격에 따라 움직이고, 반응 야코비안의 고윳값이 $O(1/\epsilon)$이라는 데서 오는 반응 제한 $\Delta t \lesssim 2\epsilon$은 격자와 무관하게 고정된 천장이다. 어느 쪽이 이기는지는 해상도가 정한다 — 그리고 전선을 제대로 분해하려면 전선 두께 $\sqrt{D\epsilon}$보다 훨씬 잘게, 즉 $h \ll \sqrt{D\epsilon}$로 잡아야 하므로 $h^2/(4D) \ll \epsilon/4$가 되어 실제 계산에서는 확산 쪽이 먼저 걸린다. 반응 천장이 이기게 하려면 전선을 한두 셀에 뭉개야 하고, 그러면 전도 속도가 눈에 띄게 틀어진다.²
• 연산자 분리. 확산과 반응을 나눠 각자 최적의 방법으로 푼다. 확산은 ADI나 암시적 방법으로 안정하게, 반응은 국소 상미분방정식이므로 격자점마다 독립적으로 작은 스텝으로 밀거나 준해석적으로 푼다. 1차 리(Lie) 분리는 스텝당 $O(\Delta t)$ 오차, 대칭 스트랑 분리는 $O(\Delta t^2)$.
• 게이트 변수 준해석 적분. 심장 모형의 표준인 러시-라슨(Rush–Larsen) 법은 게이트 방정식 $\dot g = (g_\infty - g)/\tau_g$를 계수 동결 후 지수 적분으로 정확히 풀어, 반응항 뻣뻣함의 상당 부분을 제거한다.
• 파면 해상도. 파면 두께는 $\sqrt{D\epsilon}$ 규모다. 격자 간격이 이보다 크면 파속이 격자에 의존해 틀리고, 심하면 파면이 아예 전파를 멈추는 수치적 전도 블록이 생긴다. 메시 독립성 연구를 반드시 하되, 판정 지표를 파형이 아니라 전파 속도로 잡는 것이 이 바닥의 국룰이다.

8. 관련 문서[편집]

• 벨루소프-자보틴스키 반응 · 튜링 패턴
• 극한 순환 · 분기 이론 · 카오스 이론
• 셀룰러 오토마타 · 자기조직화
• 대류-확산 방정식 · 편미분방정식
• 유한차분법 · 크랭크-니콜슨법 · CFL 조건
• 반응속도론 · 메시 독립성 연구

9. Footnotes[편집]

1. 그래서 심장 전기생리 시뮬레이션은 격자 수천만 점 규모의 본격 HPC 문제다. “심장 하나 뛰는 것 계산하는 게 뭐가 어렵냐”는 말은 파면 두께 0.5 mm를 심실 전체 해상도로 잡아야 한다는 계산을 해보기 전까지만 유효하다. ↩

2. 그래서 “격자를 성기게 하면 스텝을 키울 수 있다”는 반쯤만 맞는 말이다. 확산 상한은 $h^2$로 실제로 풀리지만 $2\epsilon$ 천장에 닿는 순간 더는 안 올라가고, 그쯤이면 이미 전선이 뭉개져 전도 속도가 틀린 답을 내고 있다. 아래 시뮬의 파라미터($\epsilon=0.2$, $h=0.12$, $D\le1.6$)에서도 확산 상한 $4.5\times10^{-3}$이 반응 상한 $1.2\times10^{-1}$보다 26배 빡빡하다. ↩