神经场

神经场方程是组织层面的模型，描述了粗粒度变量的时空演变，例如神经元群体中的突触或放电率活动。

简介

即使是一小块皮层中的神经元和突触的数量也是巨大的。因此，一种流行的建模方法是采用连续极限并研究神经网络（神经网络中，空间是连续的，并且宏观状态变量等于平均发放率）。首次尝试开发神经活动的连续近似，或可追溯到 1950 年代的 Beurle (1956) 和 1960 年代的 Griffith (1963, 1965)。通过关注给定体积的模型大脑中每单位时间被激活的神经元比例组织，Beurle 能够分析大规模大脑活动的触发和传导。然而，这项工作只处理没有难治或恢复变量的兴奋性神经元网络。Wilson和 Cowan(1972,1973) 在 1970 年代将 Beurle 的工作扩展到包括抑制性和兴奋性神经元以及不应期。Amari (1975, 1977) 在对连通性和放电率函数的自然假设下，在对经活动的连续模型中进行了进一步工作，特别是关于模式形成。Amari 考虑了局部激发和远端抑制对于具有典型皮质连接（通常称为“墨西哥帽连接”）的相互作用的抑制性和兴奋性神经元的混合群体，这是一个有效的模型。

由于这些对动态神经场理论的开创性贡献，类似的模型已被用于研究EEG节律、视觉幻觉、短期记忆和运动感知机制。

生理学动机

神经场考虑嵌入在粗粒度空间区域中的神经元群体(Wilson & Cowan 1973)。这种场的颗粒反映了在初级感觉区域中观察到的微观或宏观柱状体，例如啮齿动物的桶状皮层 (Petersen 2007) 或哺乳动物的视觉皮层 (Hubel & Wiesel 1962, Saez et al. 1998)。此外，神经领域考虑瞬时种群放电率（种群编码），即在几毫秒的某个短时间间隔内放电神经元的数量。神经场模型中的时间变量是这个时间间隔的倍数。因此，神经场在时间和空间上是粗粒度的，并表示平均场模型。重要的是要提到与神经质量模型的相似性，后者忽略了空间扩展，但涉及到关于群体编码和时间粒度的相同基本假设。

Beurle(1956)和 Wilson & Cowan(1972, 1973)是最早从数学上推导出神经场模型方程的人之一。最近 Faugeras 等人 (2009) 和 Bressloff (2009) 考虑到单个神经元的统计特性和随机动力学，给出了不同的推导。有趣的是，他们（Faugeras et al. & Bressloff）提出了一个扩展模型，同时考虑了平均活动及其方差的动态，而以前的模型只考虑了平均活动。以下段落给出了神经场模型更具启发性的推导，其中考虑了神经作用的三个主要要素：输入脉冲在化学突触处引起的树突电流、受到树突电流影响的神经元的放电以及沿轴突分支的动作电位传输。诱发的树枝状电流 I(t) 服从

[math]\displaystyle{ I(t) = \int_{-\infty}^t {\rm d}\tau^\prime h(t-t^\prime) P(t^\prime)\label{eqn_1} }[/math] (1)

具有化学受体对单个传入脉冲和传入瞬时脉冲序列 P(t) 的突触响应函数 h(t)。 单个突触的实验研究表明某些函数 h(t) 可以很好地近似突触反应行为。 在神经领域，人们考虑一个有效的平均突触响应函数，它描述了许多突触对神经群体网络中许多传入脉冲序列的平均响应。 当然，一组突触的平均时间响应可能与单个突触受体的响应函数不同，但形状可能相同。因此，方程(1)也适用于群体树突电流 I(t)，其中涉及有效参数和与突触群体响应函数 h(t)相关的群体放电率 P(t)。

在数学上，从 P(t) 到 I(t) 的转换是具有核 h(t) 的 Volterra 积分方程。 考虑微分方程通常比考虑积分方程更方便，因此可以尝试计算积分算子的逆。 这对于某些积分核函数 h(t) 是可能的。 一个简单且应用广泛的模型是 [math]\displaystyle{ h(t)=\exp(-t/\tau)/\tau }[/math] 衰减时间常数为 [math]\displaystyle{ \tau }[/math]。然后

[math]\displaystyle{ \begin{array}{lcl} \frac{dI(t)}{dt}&=&-\frac{1}{\tau}\int_{-\infty}^t {\rm d}t^\prime \exp\left(-(t-t^\prime)/\tau\right)/\tau P(t^\prime) + \int_{-\infty}^t \left({\rm d}t^\prime/dt\right)|_{\infty}^t \exp\left(-(t-t^\prime)/\tau\right)/\tau P(t^\prime)\\ &=&-\frac{1}{\tau}I(t) + \frac{1}{\tau} P(t) \end{array} }[/math]

应用链式法则和方程(1) 可以写成微分方程

[math]\displaystyle{ \hat{L} I(t)=P(t)\quad,\quad \hat{L}=\tau\frac{d}{dt}+1 . }[/math]

更详细的响应函数及其对应的微分算子是

[math]\displaystyle{ h(t)=\frac{1}{\tau_1-\tau_2}\left(e^{-t/\tau_1}-e^{t/\tau_2}\right) \quad , \quad \hat{L}=\tau_1\tau_2\frac{d^2}{dt^2}+\left(\tau_1+\tau_2\right)\frac{d}{dt}+1 }[/math]

具有较短的上升时间 [math]\displaystyle{ \tau_2 }[/math]和较长的衰减时间 [math]\displaystyle{ \tau_1 }[/math]。 对于 [math]\displaystyle{ \tau_1=\tau_2=\tau }[/math], 响应函数就是所谓的阿尔法函数

[math]\displaystyle{ h(t)=te^{-t/\tau}/\tau\quad , \quad \hat{L}=\tau^2\frac{d^2}{dt^2}+2\tau\frac{d}{dt}+1~. }[/math]

数学框架

在神经组织中电活动传播的许多连续模型中，假设突触输入电流是突触前放电率函数的函数(Wilson & Cowan 1973)。 这些无限维动力系统通常是形式上的变体

[math]\displaystyle{ \frac{1}{\alpha} \frac{\partial u(x,t)}{\partial t} = -u + \int_{-\infty}^\infty {\rm d} y w(y) f(u(x-y,t - |y|/v)) . }[/math]

图1 神经场示意图。 这样的模型不仅可以包括轴突延迟的影响，还可以包括突触和树突对放电率信号的过滤。

在这里，u(x,t) 被解释为一个神经场，表示在位置 x 和时间 t 的一组神经元的局部活动。右边的第二项代表突触输入，f 被解释为单个神经元的放电率函数。间隔 y 的神经元之间的连接强度表示为 w(y)，函数 w 通常称为突触足迹。该公式假设系统在空间上是均匀的和各向同性的。通常 w 反映全局激发 (w>0)、全局抑制 (w<0)、局部激发 - 侧抑制 (墨西哥帽) 描述，例如视觉皮层的方向调整（Somers et al. 1995, Ben-Yishai et al. 1995），或反映抑制性中间神经元的短程相互作用和兴奋性锥体细胞的长程相互作用的局部抑制 - 横向激发（逆墨西哥帽）。周期性 w 也引起了一些关注（Ben-Yishai et al. 1995, Laing and Troy 2003）。参数 [math]\displaystyle{ \alpha }[/math] 是突触的时间衰减率。空间积分下 u 的延迟参数表示由于信号在距离 y 上传播的有限速度引起的轴突传导延迟（Wilson & Cowan 1972；Nunez 1974；Jirsa & Haken 1997）；即 |y|/v 其中 v 是沿轴突纤维的动作电位的速度。最近的扩展涉及轴突传输速度 v 的分布（Atay & Hutt 2006）。

发射率函数有多种自然选择，最简单的是 Heaviside 阶跃函数。 在这种情况下，神经元最大程度地激发（以绝对不应期设定的速率）或根本不激发，这取决于突触活动是否高于或低于某个阈值。 在制定平均场神经方程的统计力学方法中，这种全有或全无的响应被平滑的 S 形形式所取代（Wilson & Cowan 1972；Amari 1975）。

上面的简单数学模型可以自然地扩展到描述多个群体、皮质层、尖峰频率适应、神经调节、慢离子电流以及更复杂的突触和树突处理形式，如下面的综述文章中所述。

动力学行为

在神经场模型中通常观察到的动态行为类型包括超出图灵不稳定性的空间和时间周期性模式 (Ermentrout 1979; Tass 1995)、局部活动区域，例如颠簸 (Kishimoto 1979) 和行波 (Ermentrout 1993; Pinto & Ermentrout 2001)。在后一种情况下，可以使用多电极记录和成像方法通过实验观察到相应的现象。特别是可以电刺激取自皮层的经过药物处理的组织切片 (Chervin et al. 1988; Golomb & Amitai 1997, Wu et al. 1999)、海马体 (Miles et al. 1995) 和丘脑 (Kim et al. 1995)。 在脑切片中，这些波可以采取癫痫发作期间的同步放电形式 (Connors & Amitai 1993) 和与感觉处理相关的传播兴奋 (Ermentrout & Kleinfeld 2001)。前额叶皮层的工作记忆（大脑中的信息的临时存储）(Colby et al. 1995, Goldman-Rakic 1995)，头部方向系统的表征(Zhang 1996)，以及视觉皮层的特征选择性，其中凸块的形成与特定神经元反应的调节有关（Ben-Yishai et al. 1995）。

模式形成

图2 具有短程激发和长程抑制的二维神经场模型中出现超出图灵不稳定性的六边形图案的活动曲线。

神经场模型是非线性空间扩展系统，因此具有支持模式形成的所有必要成分。 这种行为的分析通常结合线性图灵不稳定性理论、弱非线性微扰分析和数值模拟进行。 在一维中，具有墨西哥帽连通性的单一种群模型可以支持全球周期性静止模式。 对于一个以上的人口，非固定（行波）模式也是可能的。 在二维中可能会出现许多其他有趣的模式，例如螺旋波（Laing 2005）、目标波和双周期模式。 后面的这些图案采用条纹和棋盘状图案的形式，并且已被 Ermentrout & Cowan (1979) 和 Bressloff et al. (2001) 与药物引起的视觉幻觉联系起来。

图3 二维神经场模型中的空间局部化 3-bump 解。

具有短程激发和长程抑制的神经场模型也能够支持空间局部解决方案，通常称为颠簸或多颠簸。对于触发率函数是具有阈值 h 的 Heaviside 阶跃函数的情况，Amari (1977) 能够构造如下形式的显式单凸点解。

[math]\displaystyle{ u(x) = \int_0^\Delta w(x-y) {\rm d y}, \qquad u(0)=h=u(\Delta) , }[/math]

这样，在某个临界阈值以下，宽和窄的解决方案并存。在这两者中，更广泛的解决方案是稳定的。对于平滑的 sigmoidal 发射率，没有已知的封闭形式的空间局部解决方案，尽管使用最初为研究四阶图案形成系统而开发的技术已经获得了对多凹凸解决方案形式的深入了解 (Laing & Troy 2003)。静止的活动颠簸可以表现出各种动态不稳定性，包括 Hopf 分岔到空间局部振荡解或呼吸器(Folias & Bressloff 2004，Coombes & Owen 2004)。

活动凸起的一个可能的计算角色是根据空间结构化网络中凸起的峰值位置对一组刺激特征进行编码。在同构网络的情况下，允许的特征集将形成反映网络底层拓扑的连续流形（吸引子）。因此，周期性刺激特征可以由同质环网络中的活动凸起编码（Ben-Yishai et al. 1995，Zhang 1995）。然而，同质网络的后果之一是，相对于与连续流形相切的空间平移，凹凸将略微稳定。这意味着在存在任意小的噪声水平的情况下，活动颠簸会随着时间缓慢漂移（Laing & Chow 2001）。构建对噪声具有鲁棒性的神经场模型的一种方法是引入某种形式的细胞双稳态（Camperi & Wang 1998，Fall et al. 2004）。

行波

图4 具有尖峰频率自适应的一维神经场模型中的行进脉冲。

对于具有 sigmoid 发射率函数和兴奋耦合的一维模型，可以找到将激发态连接到静止态的波前 (Ermentrout & McLeod 1993)。 此外，在具有混合（兴奋性和抑制性）耦合的系统或具有自适应电流的兴奋性系统中，单独的行进脉冲也是可能的。 对于具有阈值 h 的 Heaviside 发射率函数，已经获得了许多关于行波的精确结果。 例如，在 w(x)=exp(-|x|)/2 的纯兴奋性网络中，稳定行进前沿的速度采用显式形式

[math]\displaystyle{ c = \frac{v(2h-1)}{2h-1 - 2 h v/\alpha} }[/math]

图5 切片实验中传播的电活动波，来自 Pinto Lab (经过允许)。 左：来自 16 个电极的线性阵列的局部场电位记录。 右：该字段的平滑颜色编码版本。

现在，Richardson 等人在真实的神经组织（浸泡在 GABA_A 阻滞剂印防己毒素中的大鼠皮质切片）中间接确定了波速对阈值 h 的强烈依赖性。 2005. 这些实验利用了这样一个事实，即皮质神经元具有长的顶端树突，并且很容易被电场极化，并且癫痫样爆发可以由刺激电极引发。跨切片施加的正（负）电场增加（降低）波传播速度，这与假设正（负）电场降低（增加）阈值 h 的神经场理论的理论预测一致。

行波的分岔结构可以使用所谓的埃文斯函数来分析。这最初是由 Evans (1975) 在关于霍奇金-赫胥黎类型的可兴奋神经轴突方程的稳定性定理的背景下制定的。这个复解析函数的零点决定了算子的法线谱，该算子是通过将系统关于其行波解进行线性化而获得的。神经场模型的扩展是最近才出现的，对于 Heaviside 放电率函数的特殊情况，现在已经研究了几个模型（Coombes & Owen 2004, 2005）。

大多数神经场模型中的一个常见假设是网络是同质的和各向同性的，也就是说，权重分布取决于网络内交互种群之间的距离。然而，由于在皮层表层中发现的长程水平连接的斑块性质，真正的皮层更现实地被建模为各向异性和不均匀的二维介质（Bosking et al. 1997）。重量分布的各向异性可能导致不同方向的波速变化，而不均匀性可能导致随时间变化的波剖面，并可能导致波传播失败（Bressloff 2001）。

参考文献

Internal references

• Yuri A. Kuznetsov (2006) Andronov-Hopf bifurcation. Scholarpedia, 1(10):1858.
• John W. Milnor (2006) Attractor. Scholarpedia, 1(11):1815.
• Eugene M. Izhikevich (2006) Bursting. Scholarpedia, 1(3):1300.
• Philip Holmes and Eric T. Shea-Brown (2006) Stability. Scholarpedia, 1(10):1838.
• S. Murray Sherman (2006) Thalamus. Scholarpedia, 1(9):1583.

参考 http://www.scholarpedia.org/article/Neural_fields

外部链接

• Stephen Coombes' website

近期综述文章

• P C Bressloff and S Coombes. Physics of the extended neuron. Int. J. Mod. Phys. B, 11:2343-2392, 1997.

Preprint available.

• B Ermentrout. Neural networks as spatiotemporal pattern-forming systems. Reports on Progress in Physics, 61:353-430, 1998.

Preprint available.

• V K Jirsa. Connectivity and Dynamics of Neural Information Processing. Neuroinformatics, 2:183-204, 2004.

Preprint available.

• T P Vogels TP, K Rajan, and L F Abbott LF. Neural network dynamics. Annual Reviews Neuroscience, 28:357-76, 2005.

Preprint available.

• S Coombes. Waves and bumps in neural field theories. Biological Cybernetics, 93:91-108, 2005.

Preprint available.

推荐阅读

• Electric Fields of the Brain: The Neurophysics of EEG, by Paul L. Nunez and Ramesh Srinivasan, Oxford University Press, 2006.

更多

Propagating Waves Pattern formation