케이블 방정식

1. 개요[편집]

케이블 방정식 Cable equation
형태	1차원 반응-확산형 선형 PDE
표준형	τm ∂V/∂t = λ² ∂²V/∂x² − V
뿌리	켈빈 경(1855), 대서양 해저전신 케이블
신경 적용	Hodgkin & Rushton (1946), Rall (1959~)
이산화	중심차분 + 크랭크-니콜슨법
선형 solver	토머스 알고리즘 / Hines 행렬

케이블 방정식(cable equation)은 축삭이나 수상돌기처럼 가늘고 긴 원통형 막에서 막전위 $V(x,t)$ 가 축 방향으로 퍼지고 막을 통해 새어 나가는 과정을 기술하는 1차원 선형 편미분방정식이다. 표준형은

$$\tau_m \frac{\partial V}{\partial t} = \lambda^2 \frac{\partial^2 V}{\partial x^2} - V$$

이며, 여기서 $\lambda$ 는 공간상수, $\tau_m$ 은 시간상수다. 우변의 확산항이 전류를 옆으로 퍼뜨리고 $-V$ 항이 막 누설로 그것을 계속 깎아 먹는다. 즉 대류-확산 방정식에서 대류를 빼고 1차 소멸항을 붙인 꼴이라, 수치적으로는 지극히 얌전한 문제에 속한다.

이 방정식의 기원은 신경이 아니라 전신선이다. 1855년 켈빈 경(윌리엄 톰슨)이 대서양 횡단 해저 케이블에서 신호가 왜 그렇게 뭉개지고 늦게 도착하는지를 설명하려고 세운 식이 사실상 같은 식이었다.¹ 인덕턴스를 되살린 완전한 전송선 이론의 전신 방정식(telegrapher’s equation)에서 유도항을 지우면 그대로 케이블 방정식이 된다. 신경생리에 본격 도입한 것은 호지킨과 러슈턴(1946)이고, 이것을 수상돌기 나무 전체로 확장해 하나의 이론 체계로 세운 사람이 윌프리드 랄(Wilfrid Rall)이다.

2. 등가회로에서의 유도[편집]

막을 단위길이당 축전기 $c_m$ 과 누설저항 $r_m$ 이 병렬로 붙은 사다리 회로로 놓고, 그 사다리를 단위길이당 축저항 $r_a$ 로 연결한다. 축방향 전류는 옴의 법칙에서

$$i_a = -\frac{1}{r_a}\frac{\partial V}{\partial x}$$

이고, 축전류가 줄어든 만큼이 막으로 빠져나가므로 $i_m = -\partial i_a/\partial x$ 다. 막전류는 축전 성분과 누설 성분의 합 $i_m = c_m \partial V/\partial t + V/r_m$ 이니, 둘을 합치면

$$\frac{1}{r_a}\frac{\partial^2 V}{\partial x^2} = c_m \frac{\partial V}{\partial t} + \frac{V}{r_m}$$

이 된다. 양변에 $r_m$ 을 곱하면 서두의 표준형이 나오고, 두 상수의 정체가 드러난다.

$$\lambda = \sqrt{\frac{r_m}{r_a}} = \sqrt{\frac{d\,R_m}{4 R_a}}, \qquad \tau_m = r_m c_m = R_m C_m$$

여기서 $R_m$ 은 비막저항($\Omega\cdot\mathrm{cm}^2$), $C_m$ 은 비막용량($\mu\mathrm{F/cm^2}$), $R_a$ 는 세포질 비저항($\Omega\cdot\mathrm{cm}$), $d$ 는 지름이다. 두 상수가 말해주는 것이 이 이론의 거의 전부다.

• $\tau_m = R_m C_m$ 은 지름과 무관하다. 전형적으로 $R_m \sim 10^4\ \Omega\,\mathrm{cm^2}$, $C_m \approx 1\ \mu\mathrm{F/cm^2}$ 이므로 $\tau_m \sim 10$ ms 급이다.
• $\lambda \propto \sqrt{d}$ 다. 굵을수록 신호가 멀리 간다. 오징어 거대축삭(지름 0.5 mm급)에서는 $\lambda$ 가 수 mm에 이르고, 가느다란 수상돌기에서는 수백 $\mu$m 수준으로 떨어진다.

무차원화 $X = x/\lambda$, $T = t/\tau_m$ 을 쓰면 방정식은 파라미터가 하나도 없는 $\partial V/\partial T = \partial^2 V/\partial X^2 - V$ 로 정리된다. 서로 다른 굵기·재질의 케이블이 전부 같은 방정식으로 붕괴한다는 뜻이고, 실제로 신경 실무에서는 케이블 길이를 cm가 아니라 전기긴장 길이(electrotonic length) $L = \ell/\lambda$ 로 말한다.

3. 정상상태 해와 경계조건[편집]

시간 미분을 0으로 놓으면 $d^2V/dX^2 = V$ 라는 2계 상미분방정식만 남고, 해는 $e^{\pm X}$ 의 조합이다. 원점에 $V_0$ 를 물려놓았을 때 답이 경계조건에 따라 이렇게 갈린다.

• 반무한 케이블: $V(X) = V_0 e^{-X}$. 공간상수 한 칸마다 $1/e \approx 37\%$ 로 줄어든다. $\lambda$ 의 정의가 바로 이것이다.
• 밀봉 말단(sealed end, $dV/dX = 0$): $V(X) = V_0 \cosh(L-X)/\cosh(L)$.
• 개방·손상 말단(killed/open end, $V(L) = 0$): $V(X) = V_0 \sinh(L-X)/\sinh(L)$.

같은 그림을 임피던스 쪽에서 보면 입력저항이 나온다. 반무한 케이블의 입력저항은 $R_\infty = \sqrt{r_m r_a}$ 이고, 유한 케이블에서는

$$R_{\rm in} = R_\infty \coth(L) \ (\text{밀봉}), \qquad R_{\rm in} = R_\infty \tanh(L) \ (\text{개방})$$

이다. $L \to \infty$ 면 둘 다 $R_\infty$ 로 수렴하고, $L \to 0$ 이면 밀봉 말단은 발산(전류가 나갈 데가 없다), 개방 말단은 0으로 간다. 실험에서 재는 값은 늘 이 둘 사이 어딘가라, 말단 처리를 잘못 가정하면 $R_m$ 추정이 통째로 틀어진다.

시간에 의존하는 해도 닫힌 형태가 있다. 무한 케이블에 $t=0$ 에서 전하를 순간 주입하면

$$V(X,T) \propto \frac{1}{\sqrt{T}}\exp\!\left(-\frac{X^2}{4T} - T\right)$$

로, 확산 커널에 지수 감쇠를 곱한 꼴이다. 여기서 중요한 결론이 하나 나온다. 수동 케이블에는 진짜 의미의 전도 속도가 없다 — 어떤 거리에서도 신호는 원리적으로 즉시 시작하고, 다만 정점이 늦게·낮게 온다. 속도라는 개념은 능동 채널이 붙어 진행파가 생겨야 비로소 정의된다.

4. Rall의 3/2 제곱 법칙[편집]

수상돌기는 원통 하나가 아니라 계속 갈라지는 나무다. 랄은 여기서 놀라운 축약을 찾아냈다. 반무한 케이블의 입력 컨덕턴스가

$$\frac{1}{R_\infty} = \frac{\pi}{2\sqrt{R_m R_a}}\, d^{3/2}$$

즉 지름의 3/2제곱에 비례하므로, 어떤 분지점에서 자식 가지들이

$$\sum_i d_i^{3/2} = d_{\rm parent}^{3/2}$$

를 만족하면 부모 쪽에서 본 임피던스가 분지 전후로 정확히 정합된다. 여기에 (i) 모든 말단 가지가 같은 전기긴장 거리 $L$ 에서 끝나고 (ii) $R_m, R_a, C_m$ 이 균일하다는 조건을 더하면, 가지가 수백 개인 수상돌기 나무 전체가 단 하나의 등가 원통과 전기적으로 구별되지 않는다. 이것이 3/2 제곱 법칙과 등가 원통(equivalent cylinder) 축약이며, 컴퓨터가 귀하던 시절 수상돌기 이론을 손으로 풀 수 있게 만든 물건이다.

물론 실제 뉴런이 이 조건을 정확히 지키는 경우는 드물다.² 오늘날의 의미는 축약 도구라기보다 형태가 전기적 기능을 어떻게 결정하는가를 보여주는 기준선이다.

5. 수치해법과 Hines 행렬[편집]

케이블 방정식을 푸는 표준 절차는 유한차분법 교과서 그대로다. 공간을 등전위 구획으로 자르고 2차 중심차분

$$\left.\frac{\partial^2 V}{\partial x^2}\right|_j \approx \frac{V_{j-1} - 2V_j + V_{j+1}}{\Delta x^2}$$

을 쓴 뒤, 시간은 크랭크-니콜슨법이나 후진 오일러로 암시적으로 전진시킨다. 명시적 도식은 확산항 때문에 $\Delta t \lesssim \Delta x^2/(2D)$ 라는 잔인한 제약을 받으므로 실무에서는 거의 쓰지 않는다. 한 줄짜리 케이블이면 결과 행렬은 삼중대각이라 토머스 알고리즘으로 $O(N)$ 에 끝난다.

문제는 분지가 있을 때다. 나무 구조에서는 한 구획이 자식을 여럿 가질 수 있어 행렬이 더 이상 삼중대각이 아니다. 마이클 하인스(1984)는 구획을 부모가 자식보다 항상 앞서도록 번호 매기면(뿌리부터 깊이 우선 전위 순회) 가우스 소거에서 채움(fill-in)이 전혀 발생하지 않아 여전히 $O(N)$ 에 정확히 풀린다는 것을 보였다. 이 정렬로 만든 희소행렬을 Hines 행렬이라 부르고, NEURON을 비롯한 구획 모형 시뮬레이터의 심장이 정확히 이 루틴이다.

공간 이산화 간격은 “구획 길이를 국소 $\lambda$ 의 0.1배 이하”로 잡는 것이 관례다.³ 시간 적분에서는 크랭크-니콜슨이 2차 정확도를 주지만 급한 전위 변화에서 진동(ringing)이 생길 수 있어, NEURON의 기본값은 안전한 후진 오일러이고 2차 정확도는 옵션으로 켜게 돼 있다.

6. 능동 케이블로 넘어가기[편집]

지금까지의 모든 결과는 막이 수동 저항이라는 가정 위에 서 있다. 실제 막에 전압 의존 채널이 붙으면 $-V/r_m$ 자리에 비선형 이온전류 $I_{\rm ion}(V, m, h, n)$ 이 들어앉고, 방정식은

$$\frac{a}{2R_a}\frac{\partial^2 V}{\partial x^2} = C_m \frac{\partial V}{\partial t} + I_{\rm ion}$$

가 된다. 이 순간 계는 선형 확산에서 흥분성 매질로 바뀌고, 감쇠 없이 일정 속도로 전파하는 진행파 — 즉 활동전위 — 가 해로 등장한다. 정량적 형태는 호지킨-헉슬리 모형이, 위상평면 직관은 피츠휴-나구모 모형이 담당한다. 게이팅 변수의 시상수가 막 시상수보다 두 자릿수 작아서 계는 전형적인 강성 방정식이 되고, 그래서 게이트는 지수 적분기로, 전압은 암시적 도식으로 푸는 연산자 분리가 국룰이다.

그래도 수동 케이블 이론은 여전히 살아 있다. 시냅스 입력이 세포체까지 얼마나 살아서 도달하는지, 수상돌기 필터링이 어느 주파수를 깎는지, 전극으로 잰 $R_{\rm in}$ 에서 $R_m$ 을 어떻게 역산할지 — 이런 질문은 전부 선형 케이블 해로 답한다.⁴ 능동 채널을 다 켜놓고 시뮬레이션을 돌린 뒤 “왜 이런 결과가 나왔는지” 설명해야 할 때, 결국 꺼내 드는 것은 $\lambda$ 와 $L$ 두 숫자다.

7. 관련 문서[편집]

• 활동전위 · 호지킨-헉슬리 모형 · 피츠휴-나구모 모형
• 흥분성 매질 · 계산 신경과학 · 구획 모형
• 전송선 이론 · 대류-확산 방정식 · 편미분방정식
• 크랭크-니콜슨법 · 토머스 알고리즘 · 희소행렬
• 강성 방정식 · 지수 적분기 · NEURON

8. Footnotes[편집]

1. 켈빈의 결론은 “신호 지연이 케이블 길이의 제곱에 비례한다”는 것이었고, 이는 확산 방정식의 $x^2 \sim 4Dt$ 스케일링 그 자체다. 당시 케이블 회사 기술고문이던 휘트하우스는 이 계산을 무시하고 전압을 2000 V까지 올려 밀어붙였다가 1858년 케이블의 절연을 태워 먹었다. 수치해석 무시의 역사적 청구서 중에서도 액수가 큰 편이다. ↩

2. 실측 형태를 3D 재구성해 놓고 각 분지점에서 $\sum d_i^{3/2} / d_{\rm parent}^{3/2}$ 를 찍어 보면 세포 종류에 따라 1 근처에서 꽤 넓게 흩어진다. 말단이 같은 $L$ 에서 끝나야 한다는 조건은 더 잘 깨진다. “등가 원통으로 줄여도 되나요?”라는 질문의 정직한 답은 “재구성 데이터로 직접 확인해 보라”다. ↩

3. NEURON의 d_lambda 규칙은 정상상태 $\lambda$ 가 아니라 100 Hz 교류 신호의 감쇠 길이를 기준으로 구획 수를 정한다. 시냅스 전위처럼 빠른 입력은 DC보다 훨씬 짧은 거리에서 깎이기 때문이다. 격자를 두 배로 늘려도 답이 안 바뀌는지 보는 메시 독립성 연구를 신경 모형에서도 똑같이 해야 한다. ↩

4. 수상돌기는 통과대역이 낮은 저역통과 필터다. 원위부 시냅스가 만든 날카로운 전위는 세포체에 도달할 즈음이면 진폭이 깎이고 봉우리가 뭉툭해진다. “가까운 시냅스가 먼 시냅스보다 세다”는 흔한 요약은 반은 맞고 반은 틀렸는데, 실제 뉴런은 원위부 채널 밀도를 올려 이 감쇠를 부분적으로 보상하기 때문이다. ↩