神经复杂性经典综述:涌现的复杂神经动力学
导语
大脑活动呈现的丰富时空活动模式,构成了其适应性行为的基础。理解大脑的千亿数量级的神经元和万亿数量级的突触是如何以灵活的方式产生如此多种皮层连接机制,仍然是神经科学中的一个基本问题。一个可能的解决方案指向了动态系统中产生涌现的复杂现象的普遍机制,尤其指向了复杂系统二阶相变的临界点。这篇综述文章回顾了支持“大脑自然地处于临界状态”这一观点的理论和实证研究结果,以及该观点怎样帮助我们更好地理解大脑。文章于2010年发表于Nature Physics,进一步拓宽了使用物理方法对神经复杂性展开研究的思路,目前引用量已经超过千次,启发了近年来诸多工作。
研究领域:神经复杂性,大脑临界性,神经动力学,涌现
论文题目:Emergent complex neural dynamics论文地址:https://www.nature.com/articles/nphys1803
目录
一、概要
二、定义问题
三、涌现现象
四、自发的大脑活动是复杂的
五、涌现的复杂性总是处于临界态
六、小尺度:大脑皮层的“地震”
七、大尺度:涌现的复杂性
八、总结与展望
一、概要
一、概要
二、定义问题
二、定义问题
本文讨论的主要问题属于该领域中的“搁置问题”类别,即大量相互连接的神经元如何以灵活的、自组织的方式产生一系列特定行为。尽管有很多口语化解释,但当构建详细模型来解释这一问题时,问题中三个关键的组成部分都被系统性地违背了。要么是问题(1)——模型是所关注的神经结构的低维版本,要么是问题(2)——它只能产生单一行为(在系统中预设),因此造成问题(3)——它无法灵活地做其他任何事情。只有通过任意地改变神经元连接,大多数当前的模型才能执行一系列行为。这需要一种类似“辅助大脑”的机制,用于掌握和决定哪些神经元连接需要被重新布线。当然,目前还没有提供这种辅助大脑的证据,这就是一个迫切需要回答却鲜少被提及的问题。
出于各种历史和概念上的原因,我们对这个问题的方法仍然来自于连接主义范式,这种范式将神经动力学限制在只能通过连接产生,从而将我们的研究引导向了一个无效的方向。尽管集体属性已经被提及了很久,但对大多数人来说,它们的相关性仍然是次要的。即使是美国生物物理学家霍普菲尔德在1982年的开创性论文《神经网络和具有涌现集体计算能力的物理系统》中的呼吁,似乎也已被遗忘,这篇论文的观点可能被计算思想的吸引力和早期热潮所取代。“高等动物的大脑区域的许多结构,必须由大量具有明确功能定义的简单局部连接的扩展所构成。简单神经连接与更高级神经系统的复杂计算属性之间的桥梁,可能是由大量简单处理单元的集体行为中自发涌现的新计算能力。”
三、涌现现象
三、涌现现象
让我们回顾一下涌现现象是什么。涌现指的是大型复杂系统所展现出来的意料之外的集体时空模式。在这里,“意料之外”表明我们无法(不论数学还是其他方式)从描述系统各个部分动力学的方程中推导出这种涌现模式。正如相关研究详细讨论的那样,复杂系统通常是由相互作用的元素组成的大型集合体,每个元素都表现出某种非线性动力学。不深入细节,我们也知道相互作用也可以是间接的,例如通过某种平均场。通常,能量会进入系统,因此存在某种驱动力。这三个重点特征(即1.系统存在大量相互作用的元素;2.每一个元素的动力学是非线性的;3.系统中有能量输入,存在某种驱动力)是系统在某一点上展现出涌现行为的必要非充分条件。只要每个单独元素的动力学是非线性的,其他细节就不那么重要;例如,这些元素可以是人类,受食物和其他能源驱动,最终形成某种集体政治或社会结构。无论涌现出何种类型的结构,如果上述三个重点特征之一缺失,这种涌现现象很可能不会出现。例如,已经确定少数孤立的线性元素无法产生意料之外的行为(实际上,这种情况下,一切都可以在数学上预测)。
四、自发的大脑活动是复杂的
四、自发的大脑活动是复杂的
近期对大脑小规模和大规模节律的研究表明,健康大脑的自发动力学并非完全由随机活动模式或周期性振荡构成。通过对无明确输入情况下神经动力学的统计特性进行深入分析,揭示了之前忽视的复杂活动模式,这些模式曾被视为背景噪声动力学。实际上,不论是通过头皮电活动(脑电图)、功能性磁共振成像(fMRI)技术、振荡活动同步化,还是局部场电位峰值的统计特征来观察,大脑活动始终在本质上是无节律的。
已经有多项研究指出,在健康状态下,大脑活动的各个时间尺度(例如神经元活动的频率、节律等)都是重要的,没有某种特定时间尺度的特征可以像平均数那样具有代表性。这就是所谓“1/f噪声”现象,即各种频率的神经活动强度与其频率的倒数成正比。这是一种典型的幂律分布,或称为无标度特性(标度即有代表性的特征尺度),在人类认知、人类活动、动物活动中均有体现。然而直到最近,大脑无标度动力学的起源才得到了充分的研究,这可能是由于无标度特性(也称标度不变性)在自然界中的普遍存在。目前,人们对此越来越感兴趣,并认识到对大脑自发模式进行重新解读至少有两方面的重要意义。自发模式的存在为大脑电路组织提供了重要线索,因为过去的理论很难解释这些新发现。此外,我们观察到的复杂动力学似乎赋予了大脑以前未被认识到的鲁棒性质。这些方面将在下一节中以两个不同的尺度进行回顾。
五、涌现的复杂性总是处于临界态
五、涌现的复杂性总是处于临界态
图-框 处于临界点的复杂现象。二维伊辛模型的数值模拟结果显示了三个温度下(亚临界、临界和超临界),某一时刻的自旋构型。只有在临界温度下,系统才会呈现出二阶相变特征,显示出高度异质的相关域,而亚临界和超临界条件下则呈现均匀状态。 |
六、小尺度:大脑皮层的“地震”
六、小尺度:大脑皮层的“地震”
Beggs和Plenz首次提供了令人信服的证据,表明神经元群体可以展现出临界动力学。他们首次描述了一种被称为“神经元雪崩”的大脑皮层电活动。这些集体神经元模式位于两种先前已知的皮层模式之间:一种是振荡或波状的高度同步活动,另一种是异步和不连贯的尖峰放电。通常,每个雪崩涉及可变数量的神经元。奇特的是,这些雪崩遵循的统计模式破局特色。平均而言,我们观察到的小型雪崩远多于大型雪崩(例如,每个神经元雪崩有很大的概率只涉及几个神经元,并且有极低的概率扩散并激活整个皮层组织(见图1))。在这些实验中,估计了一些暗示临界性的属性。其中包括雪崩大小的无标度分布,遵循一个指数近似为3/2的逆幂律,这与之前由Zapperi和同事们提出的临界分支(critical branching)过程理论的预期完全一致。雪崩的持续时间统计也遵循了一个指数近似为2的指数幂律,这与级联活动理论的预期相符。
图1 神经元雪崩是复杂的。在成熟的皮层培养网络中,神经元雪崩的大小分布遵循接近3/2的幂律(虚线表示),并显示出有限尺度缩放。图示为覆盖指定数量电极的雪崩观察到的相对概率,针对三组网格尺寸(插图中n=15、30或60个感应电极,等间隔200微米)。统计数据源自总共70小时的数据记录,每小时累积约58,000(±55,000)次雪崩(均值±标准差)。
雪崩的初始标度不变性(即无论观察的尺度如何,该现象的基本特征始终保持一致)这一观察结果,已经在其他研究实验中得到了重复验证。此外,在各种不同的环境中也报道了类似的发现,包括活体的猴子皮层和成年猫。此外,雪崩的功能意义在于它们在皮层表层发育的最早阶段就被观察到,这需要神经调节物质(多巴胺)的存在以及兴奋性和抑制性传递之间保持一定平衡。
尽管模型构建工作正在进行(见框3),但导致观察到的无标度雪崩现象的精确神经机制尚未确定。因为除了临界动力学之外,还可以通过其他几种机制产生类似的统计结果。然而,目前尚未出现任何令人信服的替代实验分析或证据。最近报道的数值证据表明,非临界的其他替代理论方案给出的逆幂律指数比实验观察到的 3/2大一个数量级(>20)。
关于雪崩起源的讨论中,一直以来的难题在于,只能复制出其大小或持续时间的概率密度。然而,如果分析其他不变量,这场讨论可以置于更严格的背景下。在这方面,最近的研究结果提供了新证据,表明神经元雪崩具有五个与临界性相符的基本特性。这些特性是:(1)触发事件和雪崩本身的动力学之间的时间尺度分离(这通过雪崩大小和时间长度的逆幂律在慢化(slow driving)中保持不变所证实);(2)尽管雪崩频率波动巨大,但雪崩大小的统计保持稳定,排除了非齐次泊松过程;(3)大型雪崩前后的雪崩概率遵循地震学中的大森定律(主震发生后余震频率随时间衰减的规律。根据大森定律,余震发生的频率与时间的倒数成正比);(4)大型雪崩后的平均雪崩大小以逆幂律衰减;(5)雪崩的空间扩散表现出分形特征。
总体而言,这些结果支持了“神经元雪崩是大脑临界性表现”的观点,并在某些情况下,明确排除了文献中讨论的大多数非临界的替代理论机制。
六、大尺度:涌现的复杂性
六、大尺度:涌现的复杂性
为了能够快进到当前进展,我们必须需要省略大量研究工作。现在学界已经认识到,大脑正在自发地创建和重塑复杂的功能网络,以应对各脑区之间的神经交流活动。近期,人们使用功能性磁共振成像对这些网络进行了研究。这一领域的研究活动相当繁忙,“对解剖和功能皮层网络组织的大量现有数据进行研究讨论,就像‘往纽卡斯尔运煤’一般。”许多详尽的综述涵盖了越来越多用于分析这堆庞大复杂数据的统计物理方法。
为了本文的目标,我们将注意力限定在自发的“静息状态” 功能性磁共振成像动力学研究。大脑的“休息”被定义为无明显大脑输入或输出,或者无显著外部刺激的状态。受试者闭眼躺着进行扫描,被告知避免入睡。从这些实验中获取到的成千上万个信号(称为BOLD,代表血氧水平依赖),每一个都反映了一个典型小区域(约十几立方毫米)的神经活动量,使我们能够描绘出整个大脑的活动。图2中展示了一个序列作为例子,由于图形限制,只记录了众多切片中的一部分。细致观察数据后,似乎已经存在重要的时空相关性,类似于云层变化过程中捕捉到的图像。有趣的是,通过相对简单线性交叉相关分析BOLD信号时,会出现一些集体群落。图2底部所示聚类就证明了这一点,它们与响应各种不同激活条件下相同区域高度匹配。因此,在静息状态下,“流动的云”(即集体的时空动力学)会激活所有在给定的主动行为期间激活的脑区。这些发现的相关性进一步强调了一个事实,即这些网络在不同受试者之间具有显著一致性,甚至在睡眠或麻醉期间也是如此。
图2 大尺度涌现的脑网络。通过对自发活动的人脑时空模式中的相互作用进行分析,揭示了涌现网络。顶部图像序列以红色/蓝色显示了4分钟连续大脑休息期间相对于平均功能性磁共振成像 BOLD的增加/减少(单个受试者的连续数据,从左上角开始,每行为1分钟,图像间隔为2.5秒)。底部图像是计算网络内小区域活动与整个大脑之间线性相关性的结果(最亮的颜色表示最强的相关性)。这些网络对应于大脑的六个主要系统:视觉、听觉、感觉运动、默认模式、执行控制和背侧注意力。所表示切片来自蒙特利尔神经研究所,坐标为顶部序列的z=18,底部图像从左到右为z=0,8,44,24,26,44。
一个自然产生的问题是,这些大脑休息模式对应的是何种已知的动力学场景。这个问题最初在三篇近期的报告中得到了探讨。在第一篇报告中,Kitzbichler等人分析了正常志愿者在休息状态下使用功能性磁共振成像和脑磁图记录的数据,利用不同空间位置之间的相位同步。他们发现两个大脑区域平均保持锁定状态的时间长度遵循一种无标度分布。这种分布也在Ising模型和Kuramoto模型的临界状态下出现。
在第二篇报告中,Fraiman等人将典型的二维Ising模型在各种温度下与休息状态大脑数据进行了比较。他们通过计算所有格点上的Ising状态间的互相关,构建了相关网络,并在相关系数超过一定阈值时,在这些点(节点)之间建立连接。对于大脑功能性磁共振成像数据,也进行了相同的计算。在比较最具描述性的网络特性后,作者们得出结论:虽然亚临界和超临界温度下的Ising网络与大脑网络有很大不同,但在临界温度下产生的网络与大脑网络“无法区分”。图3显示了其中一个进行对比的特性——这些网络连接数量(即度)的分布。可以发现,在大脑数据中存在厚重尾部(幂律分布的一个特点),只有当Ising模型处于临界温度时才能复现这种现象。此外,计算具有正负相关性位置所占比例(一个与磁化相关的变量)表明数值接近一,这与此前研究一致,该结果将此平衡视为静息状态下健康大脑功能的指标。总体来说,这些结果显示,从人类大脑的功能性磁共振成像信号中得出的相关网络,与从临界温度下Ising模型中提取的网络无法区分。
图3 从脑部功能性磁共振成像(fMRI)数据中得出的复杂网络仅在临界温度下与Ising模型相似。图示为来自fMRI大脑休息状态数据(底部图例)和在温度T=2、T=2.3和T=3时的Ising模型(顶部三个图例)在不同平均度值下的网络度(k)分布。
第三篇报告旨在揭示静息状态下功能性磁共振成像动力学机制,作者是Expert等人。他们研究了在空间粗粒化——一种临界现象中广泛使用的重整化技术——连续步骤后的两点相关函数。研究结果表明了空间自相似性,再加上功率谱在时序上频率呈1/f行为,这都意味着临界动力学的存在。
从最初的功能性磁共振成像研究以来,我们已经取得了诸多方法进展,利用这些方法可以评估正常与病理条件下(包括阿尔茨海默病、精神分裂症和癫痫)的脑功能完整性。甚至长期慢性疼痛的影响似乎超出了疼痛本身的感觉,这进一步激发我们去更深入理解大脑静息状态下,大尺度自组织结构背后的基本机制。
我们已经回顾了近年来关于神经动力学中集体复杂性涌现的关键成果。值得注意的是,至今还没有一种理论(以文章开头段落的标准)能在不引入临界性的情况下全面解释这些结果。神经科学寻求理解控制其他复杂系统的物理定律,其动机之一是期望这种普适性能给予该领域优势。与“假装生物学有其特殊性,为大脑寻找特设规律”的观点相对,对普遍原理的深入理解可能会带来突破,因为大脑必然遵循一些基本的自然规律。过去十年与当前的一个主要区别在于,我们现在有了用以检验理论的时空大脑数据:这是一个等待物理学家接受挑战、揭示集体行为底层机制的未开垦之地。
八、总结与展望
八、总结与展望
计算神经科学读书会
详情请见:计算神经科学读书会启动:从复杂神经动力学到类脑人工智能
推荐阅读