1. 개요 및 연구 배경

2022 개정 생명과학 II 신경계 단원(뉴런의 구조, 흥분의 발생과 전도, 휴지 전위 및 활동 전위)과 물리학 II 전자기학 단원(RC 회로 시상수 및 전신 방정식 케이블 이론)을 융합하는 뇌신경생물리학적 시뮬레이션입니다. 척추동물의 축삭돌기(Axon)는 슈반 세포(Schwann Cell) 또는 희소돌기아교세포(Oligodendrocyte)에 의해 지질 성분의 말이집(Myelin Sheath)으로 둘러싸여 전기적 절연층을 형성합니다.

말이집은 세포막 전기용량(Capacitance, \(C_m\))을 극적으로 낮추고 막저항(\(R_m\))을 높여 활동전위가 감쇠 없이 수동적으로 축삭 내부를 따라 빛의 속도에 가깝게 전파되도록 합니다. 말이집 사이의 노출 부위인 **랑비에 결절(Node of Ranvier)**에는 전압개폐성 나트륨 통로(\(\text{Na}_v 1.6\))가 초고밀도(\(\sim 1000/\mu\text{m}^2\))로 밀집되어 있어, 탈분극 신호를 폭발적으로 재생성하여 활동전위가 결절과 결절 사이를 뛰어넘는 **도약 전도(Saltatory Conduction)**를 구현합니다.

본 연구에서는 뇌-컴퓨터 인터페이스(BCI) 및 탈수초 신경병증(다발성 경화증 등) 신경 보철 치료를 지원하기 위해, 0-Alloc C# 10,000Hz 케이블 방정식 및 랑비에 결절 이온 통로 릴레이 서보 엔진을 수립합니다.

2. 수학적 및 생체전기 케이블 이론 모델링

수초화 축삭의 막전위 \(V(x, t)\)는 1차원 비균질 케이블 편미분방정식으로 기술됩니다:

\[ \frac{a}{2 \rho_i} \frac{\partial^2 V}{\partial x^2} = C_m(x) \frac{\partial V}{\partial t} + I_{\text{ion}}(V, x) \]

여기서 \(a\)는 축삭 반경, \(\rho_i\)는 세포내액 비저항입니다. 말이집 절간(Internode) 영역에서는 \(C_m \to C_m / N_{\text{lamellae}}\), \(I_{\text{ion}} \approx G_{\text{leak}} (V - V_L)\)로 수동적 감쇠를 따르며, 결절(Node) 영역에서는 초고밀도 \(\text{Na}_v\) 이온전류가 작동합니다:

\[ I_{\text{ion}}^{\text{node}} = \bar{g}_{\text{Na}} m^3 h (V - E_{\text{Na}}) + \bar{g}_K n^4 (V - E_K) + g_L (V - E_L) \]

결절 간 거리 \(L\)과 축삭 외경 \(D\) 간의 비례 관계 \(L \approx 100 D\) 및 내경/외경 최적 비 \(d/D \approx 0.6\)에서 전도 속도(Conduction Velocity) \(\theta\)가 극대화됩니다:

\[ \theta \propto \frac{d}{\sqrt{\rho_i C_m}} \ln\left(\frac{D}{d}\right) \]

3. 0-Alloc C# 초고속 10,000Hz 랑비에 결절 도약전도 서보

다중 랑비에 결절 격자에서 활동전위의 전파 지연과 자극 펄스를 힙 할당 없이 100µs 단위로 추적하는 0-Alloc 스택 서보 엔진입니다: 

using System;
using System.Runtime.CompilerServices;

namespace ZeliDesk.NeuralSaltatoryConduction
{
    public readonly ref struct NodeState
    {
        public readonly double MembraneVoltageMv;
        public readonly double MGate;
        public readonly double HGate;

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public NodeState(double v, double m, double h)
        {
            MembraneVoltageMv = v;
            MGate = m;
            HGate = h;
        }
    }

    public static class RanvierNodeServoEngine
    {
        private const double ENa = 55.0; // mV
        private const double EL = -70.0; // mV
        private const double GNaBar = 1200.0; // mS/cm^2 (고밀도 결절)
        private const double GL = 10.0;

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public static NodeState StepNode(
            in NodeState current,
            double electrotonicCurrentFromAdjacent,
            double dt)
        {
            double v = current.MembraneVoltageMv;
            double m = current.MGate;
            double h = current.HGate;

            // Na 게이트 전이율
            double alphaM = 0.1 * (v + 40.0) / (1.0 - Math.Exp(-(v + 40.0) / 10.0));
            double betaM = 4.0 * Math.Exp(-(v + 65.0) / 18.0);
            double alphaH = 0.07 * Math.Exp(-(v + 65.0) / 20.0);
            double betaH = 1.0 / (1.0 + Math.Exp(-(v + 35.0) / 10.0));

            m += dt * (alphaM * (1.0 - m) - betaM * m);
            h += dt * (alphaH * (1.0 - h) - betaH * h);

            double iNa = GNaBar * (m * m * m) * h * (v - ENa);
            double iLeak = GL * (v - EL);

            // 막전위 갱신 (Cm = 1 uF/cm^2)
            double dV_dt = - (iNa + iLeak) + electrotonicCurrentFromAdjacent;
            double nextV = v + dV_dt * dt;

            return new NodeState(nextV, m, h);
        }
    }
}

4. 2022 개정 교육과정 연계 및 신경의공학 응용

  • 생명과학 II: 활동 전위의 발생 메커니즘과 도약 전도의 속도학적 우위(무수신경 vs 유수신경) 정량 분석.
  • 물리학 II: 절연체와 도체 결합 구조에서의 RC 시상수 단축 및 신호 감쇠율 이론적 검증.
  • 국방 인간공학 & BCI: 고성능 군사용 엑소스켈레톤 착용자 신경근 신호 초저지연 판독 및 신경 자극 인터페이스.