시냅스 전류 동역학을 고려한 잡음 통합‑발화 뉴런의 발화율

초록

본 논문은 시냅스 전류의 동적 특성을 포함한 잡음 통합‑발화 모델의 발화율을 정확히 계산하는 해석식을 제시한다. 다중 시냅스 입력과 절대 불응기 효과까지 고려하여, 기존의 정적 시냅스 가정보다 현실적인 뉴런 발화 메커니즘을 수학적으로 기술한다.

상세 분석

논문은 먼저 전형적인 레키티컬 적분‑발화(Leaky Integrate‑and‑Fire, LIF) 모델에 시냅스 전류의 일차적 저역통과(First‑order low‑pass) 동역학을 결합한다. 시냅스 전류 I(t)는 τ_s · dI/dt = ‑I + σ ξ(t) 형태의 오리엔테이션을 갖으며, 여기서 ξ(t)는 백색 가우시안 잡음, σ는 잡음 강도, τ_s는 시냅스 시정수이다. 이러한 연속적 전류는 전압 V(t)의 미분 방정식 τ_m · dV/dt = ‑V + RI(t)와 결합되어 두 차원 마코프 과정으로 기술된다. 저자들은 이 시스템을 Fokker‑Planck 방정식으로 전이하고, 반사 경계(전압 리셋)와 흡수 경계(스파이크 임계값) 조건을 적용한다. 핵심은 전류와 전압이 결합된 확률 밀도 함수 p(V,I,t)의 정상 상태 해를 구하는 것이며, 이를 위해 라플라스 변환과 특수 함수(오일러‑맥클로린, 에어리 함수)를 활용한다.

다중 시냅스 입력을 다루기 위해 각 시냅스 채널을 독립적인 OU 과정으로 모델링하고, 총 전류는 이들의 선형 결합으로 표현한다. 이 경우 전체 시냅스 시정수들의 가중 평균 τ_eff가 새로운 효과적 시정수로 등장하며, 잡음의 분산은 각 채널의 σ_i^2의 합으로 나타난다. 절대 불응기 τ_ref는 스파이크 후 전압을 고정된 리셋값 V_r로 유지하는 기간으로, 이는 재시작 확률을 감소시키는 곱셈 인자로서 발화율 λ에 λ = 1/(τ_ref + ⟨T⟩) 형태로 포함된다. 여기서 ⟨T⟩는 시냅스 전류와 전압이 임계값에 도달하는 평균 첫 통과 시간이다.

저자들은 ⟨T⟩을 정확히 계산하기 위해 고전적 첫 통과 시간 이론을 확장한다. 특히, 전류와 전압이 결합된 2차원 확산 과정에서 경계 조건을 만족하는 Green’s function를 구축하고, 이를 시간 적분하여 평균 통과 시간을 도출한다. 결과적으로 발화율 λ는 τ_m, τ_s, σ, μ(평균 전류), V_th, V_r, τ_ref 등 모든 모델 파라미터를 명시적으로 포함하는 복합 함수 형태가 된다. 이 식은 수치 시뮬레이션과 일치함을 보이며, 특히 τ_s ≫ τ_m인 경우 전류가 느리게 변해 발화율이 크게 감소하고, τ_s ≪ τ_m인 경우 기존의 백색 잡음 근사와 동일한 결과를 재현한다는 중요한 한계를 확인한다.

이러한 분석은 신경 회로 모델링에서 시냅스 전류의 시간적 필터링 효과를 정량화하는 데 필수적이며, 다중 입력과 불응기까지 포함한 포괄적 프레임워크를 제공한다.