뇌파 동기화 억제를 위한 가변 지연 피드백과 힌드마시 로즈 진동기 네트워크

초록

본 연구는 평균장(mean‑field)으로 상호 연결된 다수의 동일한 힌드마시‑로즈(HR) 뉴런 모델에 피라가스(Pyragas)형 가변 지연 피드백을 도입해, 기존 고정 지연 방식보다 넓은 파라미터 구역에서 효과적인 탈동기화를 달성함을 보인다. 가변 지연은 피드백 신호가 과거 평균장의 값과 현재 평균장의 차이에 비례하도록 하며, 이는 불안정한 고정점(stable steady state) 억제 기법과 유사한 원리다. 수치 시뮬레이션 결과, 동기화 지표인 전역 순서 매개변수(order parameter)가 크게 감소하고, 개별 진동기의 위상 분포가 넓어지는 것이 확인되었다.

상세 분석

본 논문은 신경계 질환에서 관찰되는 과도한 동기화 현상을 억제하기 위한 수학적 모델링 접근법을 제시한다. 기본 모델은 N개의 동일 HR 진동기를 평균장 coupling 형태로 연결한 것으로, 각 진동기의 3차원 상태 변수(x, y, z)가 전형적인 스파이킹-버스트(bursting) 동작을 보인다. 기존 연구(Rosenblum·Pikovsky, 2004)는 현재 평균장과 일정 시간 τ 전의 평균장 차이에 비례하는 피라가스 피드백을 도입해 탈동기화를 유도했으며, 파라미터 K(피드백 강도)와 τ의 정확한 조합이 필요했다. 저자들은 여기서 피드백 지연을 고정값이 아니라 시간에 따라 변하는 함수 τ(t)로 확장한다. 구체적으로 τ(t)=τ0+ε·sin(ωt)와 같은 주기적 변조 혹은 상태 의존적 변조를 고려했으며, 이는 시스템이 특정 지연값에 국한되지 않고 넓은 지연 구간을 탐색하도록 만든다. 가변 지연은 비선형 제어 이론에서 알려진 ‘delay‑modulated feedback’과 유사하게, 시스템의 고유 고정점(stable steady state)을 넓은 영역에서 안정화한다. 수치 실험에서는 K를 0.11.0, τ0를 0.52.0(시간 단위) 범위로 스캔했으며, ε와 ω를 적절히 선택하면 전역 순서 매개변수 R이 0.9 수준에서 0.2 이하로 급격히 감소한다. 특히, 고정 지연 방식에서는 R이 0.4 이하로 떨어지기 위해 K와 τ가 매우 정밀히 맞춰져야 했지만, 가변 지연에서는 K와 τ의 조합이 넓어도 유사한 탈동기화 효과가 관찰된다. 또한, Lyapunov 지수 분석을 통해 가변 지연 피드백이 시스템의 최대 리아프노프 지수를 음수 영역으로 이동시켜, 동기화된 궤도 자체가 불안정해짐을 확인했다. 이러한 결과는 신경계 치료용 전기 자극(예: 깊은 뇌 자극, DBS)에서 고정된 파라미터 설정에 의존하지 않고, 환자별 혹은 시간에 따라 변하는 최적 자극을 제공할 수 있는 이론적 근거를 제공한다.