1. 引言
噪声在神经系统中是无处不在的 [1] 。它有很多不同的来源,也会产生很多不同的影响。神经元中噪声的来源主要包括通道噪声,即离子通道的随机开关,以及突触噪声,即突触随机释放神经递质 [2] 。通常,人们会认为噪声产生的影响是消极的。然而,大量的实验结果表明:在某些特定条件下,噪声对非线性系统的动力学有着积极影响 [2] 。其中两个积极影响就是诱导神经元产生随机共振和相干共振 [3] [4] 。随机共振是描述非线性系统中噪声的存在可以增加系统对外部刺激的输出响应 [5] [6] [7] 。相干共振与随机共振不同,即使在没有外部周期性刺激的情况下,噪声也可以使非线性系统呈现出类随机共振 [8] [9] [10] 。实验证明,生物神经系统中普遍存在着随机共振和相干共振。例如,噪声可以诱导小龙虾尾部神经元产生随机共振现象 [11] ,也可以诱导猫的中枢神经系统发生相干共振现象 [12] ,这表明噪声在神经系统的信息编码中起着重要作用。
一直以来,随机共振和相干共振现象都受到广泛学者的关注。例如,Benzi等人用随机共振现象来解释第四纪冰川问题 [13] 。Lee等人研究了无外部周期信号的Hodgkin-Huxley (HH)模型,发现相干共振发生在一定噪声强度处 [14] 。同时,神经系统中的不同因素对噪声诱导的随机共振和相干共振的影响也受到了广泛学者的关注。例如,随着突触时滞的增加,神经元网络会产生多次随机共振或相干共振现象 [15] ;对于较弱耦合强度的神经元网络,相位噪声能诱导产生两次相干共振现象 [16] ;对于较强耦合强度的神经元网络,相位噪声能诱导产生一次相干共振现象 [16] ;对于中等耦合强度的神经元网络,小世界网络结构可以增强随机共振现象 [17] ;在HH神经元网络中,抑制性自突触增强相干共振 [18] 。在HH神经元中,抑制性自突触增强随机共振 [3] 。
最近,神经系统中电磁感应效应对随机共振和相干共振的影响受到了广泛关注。神经系统中的电磁感应效应主要来源于膜电位的波动 [19] [20] 。有研究利用忆阻器介导电磁感应在神经系统仿真随机共振和相干共振现象 [21] [22] [23] [24] 。在这些研究中,关注的是电磁感应效应对单神经元的相干共振的影响。然而,电磁感应效应对神经元网络的相干共振的影响以及神经系统中多种因素对相干共振的影响受到的关注相对较少。
本文基于FitzHugh-Nagumo (FHN)神经元建立了一个静息态的神经元网络模型,并仿真相干共振现象,进一步研究了电磁感应效应的强度对神经元网络中相干共振的影响。
2. 模型与方法
2.1. FHN神经元模型
FHN神经元模型可以看作是HH神经元模型的简化 [25] ,其方程描述如下:
(1)
其中v表示膜电位,w表示膜恢复变量。参数ε是FHN模型的系统参数,可以保证神经元的膜电位随着时间快速变化,本文将其固定为0.01。兴奋性参数a决定了FHN神经元的放电行为,当a > 1时,模型有一个稳定的焦点,其对应于模型的静息态,如图1(a)所示。当a < 1时,模型有一个稳定的极限环,其对应于模型的周期放电态,如图1(b)所示。当a = 1时,模型发生Hopf分岔。本文将参数a固定为1.1,模型处于静息态。
Figure 1. Plots of membrane potential over time. (a) a = 0.9; (b) a = 1.1
图1. 膜电位随时间的变化图。(a) a = 0.9;(b) a = 1.1
2.2. 神经元网络模型
构建含高斯白噪声和电磁感应的神经元网络模型如下:
(2)
模型中,下标i (
)表示神经元网络中的第i个神经元。本文的研究考虑WS小世界神经元网络 [26] ,神经元的总数固定为N = 100。k0代表电磁感应强度,由磁通量φ控制的忆阻器的方程为:
,参数g表示耦合强度,
表示为第i个神经元接受到的耦合电流。
是神经元网络的连接矩阵,当
时,Cij = 0,否则,Cij = 1。
表示高斯白噪声,其统计特征为:
和
。其中D表示噪声强度,
表示狄拉克函数。方程(1)是一个含高斯白噪声的高维的随机微分方程。
本文中,将FHN模型的系统参数ε固定为0.01,a固定为1.1,与电磁感应相关的参数固定为:α = 0.1,β = 0.1,k1 = 0.1,k2 = 1.0。WS小世界网络的重连概率p固定为0.2,不同神经元间的耦合强度g固定为0.05,重连概率p和耦合强度g的其它取值不会定性影响本文的结果 [27] [28] 。
2.3. 方法
本文使用欧拉方法数值求解随机微分方程(2),积分步长为0.001。所有神经元的峰峰间期(interspike interval, ISI)的变差系数的倒数λ用于定量描述神经元网络的放电规则性。λ定义为ISI的均值与ISI的标准差的比值,即
。由定义可知:较大的λ值代表神经元网络的放电规则性强。如果λ的值随噪声强度D的增大而先增大后减小,那么神经元网络产生相干共振。
3. 结果
3.1. 当k0 = 0时,神经元网络模型中的相干共振现象
图2描刻画了当神经元网络中的电磁感应强度为0时,在不同噪声强度下神经元网络放电的时空图,每个黑点对应于神经元的一次放电。当噪声强度D很小时(例如D = 0.001),网络中神经元的放电不规则且放电率较低(图2(a))。当噪声强度D增大至D = 0.0063时,网络中神经元的放电规则性增强且放电率也增大(图2(b))。当噪声强度D进一步增大至D = 1时,网络中神经元的放电规则性变弱,但放电率继续增大(图2(c))。图2表明,神经元网络的放电规则性随着噪声强度的增大而呈现出先增大后减小的趋势,即产生相干共振现象。
为了进行定量验证噪声可以诱导神经元网络产生相干共振现象,图3描绘了在k0 = 0的条件下,λ随着噪声强度D的变化趋势(黑色圆形标注)。这条曲线显示,λ随着噪声强度D的增大而呈现出先增大再减小的趋势,并且在D约为0.0063时达到最大值,神经元网络产生了相干共振。这与图2中的分析结果相一致。
3.2. 当k0 < 0时,神经元网络模型中的相干共振现象
当神经元网络中的电磁感应强度k0 = −1时,图3描绘了k0 = −1时λ随噪声强度D的变化趋势(蓝色三角形标注)。这条曲线显示,λ随噪声强度D的增大而先增大再减小,并在D约为0.0063时达到最大
Figure 2. When the electromagnetic induction intensity k0 = 0, space-time plots of the neuronal networks for different D values. (a) D = 0.001; (b) D = 0.0063; (c) D = 1
图2. 当电磁感应强度k0 = 0时,不同噪声强度D下神经元网络的时空图。(a) D = 0.001;(b) D = 0.0063;(c) D = 1
Figure 3. Under different k0 values λ as the noise intensity D changes. k0 = −1 (blue triangle), k0 = 0 (black circle), k0 = 1 (red square)
图3. 不同k0值下λ随着噪声强度D的变化规律。k0 = −1 (蓝色三角形),k0 = 0 (黑色圆圈),k0 = 1 (红色正方形)
值,神经元网络产生了相干共振。另外,图4刻画了不同噪声强度下神经元网络放电的时空图。当噪声强度很小时(例如,D = 0.0016时),网络中神经元的放电不规则且神经元的放电率较低(图4(a))。当噪声强度增大至D = 0.0063时,神经元放电的规则性变强且神经元的放电率增大(图4(b))。当噪声强度进一步增大至D = 1时,神经元放电的规则性减弱,但神经元的放电率继续增大(图4(c))。当D = 0.0063时,神经元的放电状态最规则,神经元网络产生相干共振现象。这与图3中蓝色三角形标注的曲线一致。
Figure 4. When the electromagnetic induction intensity k0 = −1, space-time plots of the neuronal networks for different D values. (a) D = 0.0016; (b) D = 0.0063; (c) D = 1
图4. 当电磁感应强度k0= −1时,不同噪声强度D下神经元网络的时空图。(a) D = 0.0016;(b) D = 0.0063;(c) D = 1
3.3. 当k0 > 0时,神经元网络模型中的相干共振现象
图3也描绘了当神经元网络中的电磁感应强度k0 = 1时,λ随噪声强度D的变化曲线(红色正方形标注)。这条曲线显示,λ随噪声强度D增大而呈现出先增大再减小的趋势,并在D约为0.0025时达到最大值,神经元网络产生了相干共振现象。图5刻画了在不同的噪声强度下神经元网络放电的时空图。D很小时(例如,D = 0.0006时),网络中神经元的放电不规律且神经元放电率较低(图5(a))。当D增大至0.0025时,神经元放电的规则性变强且放电率增大(图5(b))。当噪声强度D进一步增大至1时,神经元放电的规则性变弱,但放电率继续增大(图5(c))。噪声强度D = 0.0025时,神经元的放电状态达到最规则,神经元网络产生了相干共振。这与图3中红色正方形标注的曲线一致。
Figure 5. When the electromagnetic induction intensity k0 = 1, space-time plots of the neuronal networks for different D values. (a) D = 0.0006; (b) D = 0.0025; (c) D = 1
图5. 当电磁感应强度k0 = 1时,不同噪声强度D下神经元网络的时空图。(a) D = 0.0006;(b) D = 0.0025;(c) D = 1
3.4. 不同电磁感应强度下神经元网络的相干共振
上述研究结果表明:在三个电磁感应强度下(即k0 = 0、k0 = −1、k0 = 1),神经元网络均会产生相干共振现象。此外,针对这三个电磁感应强度k0,λ均有先增大后减小的趋势,这进一步说明,神经元网络在三个不同的电磁感应强度下均会产生相干共振现象。同时,电磁感应强度k0从−1增大至1,发生相干共振现象对应的λ值也会随之增大。这说明相干共振的强度会随着电磁感应强度的增大而增大。
为了进一步描述电磁感应强度k0对神经元网络的相干共振现象的影响,图6刻画了λ随电磁感应强度k0和噪声强度D的变化趋势。图6中的红色区域表示较大的λ值,蓝色区域表示较小的λ值。图6表明,在电磁感应强度k0保持不变的情况下,λ随D的增大而呈现出先增大后减小的趋势。因此,神经元网络在不同电磁感应强度下均会产生相干共振现象。同时,随着电磁感应强度的增加,产生相干共振现象对应的D值会减小,λ值会增大。这表示随着网络中电磁感应强度的增大,产生相干共振现象对应的噪声强度会减小,同时相干共振的强度会随之增强。根据之前的研究,诱导相干共振的噪声强度减小以及相干共振的强度增大是因为随着电磁感应强度的增大,单个FHN神经元的Hopf分岔点向静息态移动 [28] 。
Figure 6. The indicator λ varies with the changes in noise intensity D and electromagnetic induction intensity k0
图6. 指标λ随着噪声强度D和电磁感应强度k0的变化规律
4. 结论
本文研究了噪声诱导静息态神经元网络产生相干共振现象,研究结果发现在不同的电磁感应强度下,中等强度的高斯白噪声均可以诱导静息态神经元网络发生相干共振。尽管前人也研究了噪声能诱导神经元发生相干共振,但是前人的研究中发现噪声可以诱导单神经元发生相干共振 [29] 。本研究揭示了噪声可以诱导神经元网络发生相干共振,本研究进一步丰富了神经元网络的随机动力学。
此外,本文还揭示了电磁感应的强度对神经元网络的相干共振有重要影响。电磁感应强度的增大会使发生相干共振对应的噪声强度的临界值减小,同时λ的最大值会增大,即随着电磁感应的增强,相干共振的强度也会随之增强。本文的结果揭示调节电磁感应强度是调控神经元网络相干共振现象的重要手段。