《Biological Cybernetics》:Adaptive conduction delays and phase locking in spiking Haken Lighthouse networks
编辑推荐:
研究人员利用Haken Lighthouse模型作为神经动力学的一种解析可追踪的事件驱动描述,发展了延迟脉冲网络中相位锁定活动的理论。对于具有固定延迟的网络,研究人员推导了相位锁定状态的自洽条件,以及直接根据脉冲时间扰动表述的线性稳定性理论。该框架通过一个延迟
研究人员利用Haken Lighthouse模型作为神经动力学的一种解析可追踪的事件驱动描述,发展了延迟脉冲网络中相位锁定活动的理论。对于具有固定延迟的网络,研究人员推导了相位锁定状态的自洽条件,以及直接根据脉冲时间扰动表述的线性稳定性理论。该框架通过一个延迟自突触、一个相互耦合的双细胞网络以及具有距离依赖性耦合和传导延迟的空间结构化环进行了展示,其中循环对称性使得稳定性可分解为傅里叶模式。随后,研究人员引入了一种活动依赖的白质可塑性规则,其中髓鞘调节轴突传导速度,进而调节通信延迟。这自然地产生了一个具有状态依赖性延迟的慢-快系统,其中冻结的相位锁定分支组织了自适应动力学。该可塑性规则选择相称的延迟-周期关系,为同步、其他频率锁定状态、竞争性相位锁定模式之间的慢切换以及异质延迟组织成分立的延迟-周期类别提供了机制。对事件驱动网络的直接模拟支持了分析预测,并展示了适应性传导如何重塑延迟脉冲网络的吸引子结构并产生长时间尺度转变。这些结果为研究活动依赖性髓鞘形成如何调节脉冲神经系统中的时间协调、同步性和通过相干性进行通信提供了易于处理的数学框架。
理论神经科学的核心目标之一是构建可处理的脉冲神经元网络模型。简单的脉冲网络模型如积分-发放单元虽然能展现丰富的涌现行为,但数学分析非常困难;而发放率模型和耦合相位振子模型虽易处理,却会平均化或压抑神经通信中基于事件的结构。Haken Lighthouse模型由Hermann Haken提出,它通过受突触输入驱动的相位类变量旋转决定脉冲发放,因此既保留脉冲事件结构,又利于解析研究。尽管近期研究在Lighthouse模型中发现了同步、行波、局部凸块等自组织行为,但缺乏对传导延迟的充分处理。然而,延迟是神经系统通信的基本成分:轴突传导、突触加工、树突滤波和长程解剖连接都会引入时间约束,影响同步、相位锁定和多稳定性。同时,白质可塑性研究表明,髓鞘形成具有活动依赖性,可缓慢调节轴突传导速度。因此,研究人员将Lighthouse模型拓展到包含传导延迟的情形,并引入活动依赖白质可塑性规则,以探讨适应性延迟如何影响脉冲网络的时序协调。
论文有两个主要贡献:第一,建立了固定延迟下相位锁定态的存在性与线性稳定性理论,直接基于脉冲时间扰动表述;第二,将该固定延迟理论作为骨干,提出一个现象学白质可塑性规则,使传导速度(和延迟)随时间缓慢演化。通过慢-快系统分析,研究人员证明了该规则可选择相称的延迟-周期关系,进而促进同步、频率锁定、慢切换,并将异质延迟组织为离散类别。这些结果为理解活动依赖髓鞘形成如何调节神经同步和长程相位通信提供了数学框架。
主要关键的技术方法包括:(1)事件时间映射:将Lighthouse模型从一次脉冲到下一次脉冲积分,导出精确的脉冲时间映射方程;(2)相位锁定方程:令脉冲时间满足\(T_i^m=mT+\phi_i T\),将存在性化为N维方程组;(3)线性稳定性分析:对脉冲时间扰动线性化,得到特征方程\(\mathcal{E}(\lambda)=0\),并利用傅里叶展开计算\(G_{ij}(\lambda)\);(4)循环矩阵对角化:对环网络使用空间傅里叶模式分解稳定性矩阵;(5)慢-快奇异扰动框架:将白质可塑性视为慢变量,以固定延迟分支为临界流形;(6)事件驱动数值模拟。样本队列不适用。
4 单节点自突触(A single node with an autapse):研究人员考虑具有瞬时和延迟自反馈的单神经元。正则发放分支由\(\mathcal{H}(T,\tau)=0\)隐式定义。线性稳定性分析发现,分支折叠点对应双零条件\(\mathcal{E}(0)=0=\mathcal{E}'(0)\),动态不稳定性对应\(\lambda=\mathrm{i}\omega\)处\(\mathcal{E}(\mathrm{i}\omega)=0\),其表现为发放间隔的振荡调制,即“抖动”。数值模拟确认了正则发放分支的折叠、多稳定性以及高周期脉冲间隔模式。
5 双向双节点网络中的相位锁定态(Phase-locked states in a reciprocal two-node network):对于对称耦合、公共延迟的双节点网络,同步解(\(\phi=0\))和反同步解(\(\phi=1/2\))满足单方程\(\mathcal{H}(T,\psi)=