《Chaos, Solitons & Fractals》：Self-organized criticality and structural robustness: Node-failure analysis of the sandpile model on rewired networks
编辑推荐：
自组织临界性（Self-Organized Criticality, SOC）描述了由多相互作用元素组成的驱动耗散系统中临界行为的自发涌现。本研究探讨了当支撑SOC动力学的基底经历渐进式结构退化时，这种临界行为能在多大程度上持续存在。研究人员采用了一种保度的重
自组织临界性（Self-Organized Criticality, SOC）描述了由多相互作用元素组成的驱动耗散系统中临界行为的自发涌现。本研究探讨了当支撑SOC动力学的基底经历渐进式结构退化时，这种临界行为能在多大程度上持续存在。研究人员采用了一种保度的重连（rewiring）程序，将规则晶格转变为从弱随机化到高度混合的网络拓扑，随后通过随机移除节点并消除封闭簇（即不含边界节点的连通分量）来模拟结构损伤。在Bak-Tang-Wiesenfeld（BTW）沙堆模型框架下，研究人员分析了不同拓扑结构对节点失效的结构与动力学鲁棒性的影响。结果表明，拓扑随机化显著延缓了结构碎裂：规则晶格在大约40%节点被移除后崩溃，而适度重连的网络能将功能性活动维持至约60%节点被移除。这种增强的鲁棒性源于到耗散边界的有效路径缩短，以及开放簇结构的更渐进退化。此外，平均颗粒密度和颗粒运动等状态与输运观测量在整个退化过程中提供了稳健的性能指标，即使在雪崩分布因有限尺寸效应偏离理想幂律形式时依然有效。这些发现表明，即使适度偏离严格正则性也能增强耗散雪崩动力学在结构损伤下的恢复力，这对理解神经、电力和金融系统中的级联现象具有启示意义。
论文解读文章
一、研究背景与问题提出
自组织临界性（Self-Organized Criticality, SOC）是复杂性科学的核心理论框架，描述了由大量相互作用单元构成的驱动耗散系统如何自发演化至临界态，产生跨越多个尺度的间歇性爆发活动。Bak-Tang-Wiesenfeld（BTW）沙堆模型作为SOC的经典范式，以其简洁的局域规则——缓慢添加颗粒、超过阈值时发生倾倒（toppling）、级联传播直至稳定——揭示了丰富的集体行为。然而，现实世界中展现雪崩动力学（avalanche dynamics）的系统，如神经网络、电力基础设施和金融网络，其底层交互拓扑既非规则也非结构不变。这些系统可能遭受故障而改变连接性，进而影响其动力学行为。
现有研究虽已证实网络拓扑显著影响SOC的雪崩统计特征，但对于SOC动力学如何在渐进式结构退化下保持鲁棒性仍缺乏系统性认识。特别是，当节点损失导致图碎裂、耗散路径改变和有限尺寸截断时，传统的幂律雪崩指数变得难以解释。因此，评估鲁棒性不仅需要追踪雪崩统计量，还需结合更直接的状态与输运观测量。本研究旨在填补这一空白，在一个统一可控的框架内比较不同拓扑对相同失效过程的响应。
二、研究方法概述
研究人员建立了一套结合拓扑重连与结构退化的分析框架。首先，采用基于Croes置换原理的保度重连（degree-preserving rewiring）算法，将二维方形晶格逐步转化为弱随机化乃至高度混合的类随机网络，同时保持所有节点度数不变及图的简单性。其次，实施两步退化过程：随机移除指定比例的节点，再删除所有不含原始晶格边界节点的封闭簇（closed clusters），仅保留允许颗粒耗散的开放簇（open clusters）。最后，根据幸存节点的实际邻居数更新其临界阈值（threshold），确保沙堆动力学与新结构兼容。所有模拟均在128×128大小的晶格上进行，并对多种重连率和退化率组合进行统计分析。
三、研究结果
3.1 重连对预退化沙堆动力学的影响
在引入退化前，研究人员验证了重连方法能否复现已有文献中的定性趋势。结果显示，雪崩持续时间分布的幂律指数在低重连水平下迅速增大，随后随拓扑无序程度增加趋于稳定，表明短程捷径的引入使小规模雪崩相对大规模雪崩的频率增加。与此同时，颗粒密度（grain density）仅随重连轻微变化，说明重连的主要效果并非改变系统的平均占据率，而是重塑雪崩活动的分布模式。
3.2 结构退化的拓扑依赖性
研究人员对比了规则晶格与不同程度重连网络在节点移除下的结构碎裂过程。规则晶格在大约40%节点被移除后，最大开放簇的相对规模急剧下降，系统结构性崩溃。相比之下，经过适度重连的网络能够将大型连通结构维持至约60%节点被移除，碎裂过程更为渐进。这一差异源于重连引入的长程捷径（long-range shortcuts）使得更多节点与边界相连，从而减少了封闭簇的形成。
3.3 动力学鲁棒性的增强机制
通过分析颗粒运动热力图，研究人员发现退化导致结构碎裂后，颗粒输运集中在最大的幸存开放簇内，远离边界的节点被频繁穿越。重连通过缩短节点到耗散边界的有效路径长度，使颗粒能够更快离开系统，从而延缓了动力学功能的丧失。即使在雪崩分布因有限尺寸效应偏离理想幂律时，平均颗粒密度和颗粒运动次数仍可作为稳健的退化指示器。
3.4 雪崩分布变形与有限尺寸效应
随着退化加剧，雪崩分布逐渐偏离纯幂律形式，出现明显的有限尺寸截断。研究人员采用截断分布描述符而非单纯依赖幂律拟合来表征这一变形，发现重连网络的雪崩分布退化更为平缓，进一步佐证了其动力学鲁棒性。
四、总结与讨论
本研究表明，拓扑随机化显著增强了BTW沙堆模型在节点失效下的结构与动力学鲁棒性。具体而言，适度重连使系统能够在高达60%节点被移除时仍维持功能性活动，而规则晶格在40%节点损失后即告崩溃。这种增强归因于重连带来的到耗散边界的更短路径以及更渐进的结构退化过程。
研究结论翻译：研究人员的结论表明，拓扑对鲁棒性具有显著影响。适度的拓扑随机化大幅延迟了碎裂，并使SOC相关的动力学特征在比规则网格高得多的节点移除水平下得以持续。这种增强的鲁棒性与到耗散边界的有效路径缩短以及主要动力学指标的更渐进退化相关联。更广泛地说，研究人员的发现表明，即使对严格正则性的有限偏离也能改善耗散雪崩动力学在结构损伤下的恢复力。这些结果对理解神经、电力和金融系统中的级联失效现象具有启示意义，并为设计更具韧性的复杂网络提供了理论依据。
