2025-11-16T12:46:12.439899

Extreme events scaling in self-organized critical models

Quadir, Jafri
We study extreme events of avalanche activities in finite-size two-dimensional self-organized critical (SOC) models, specifically the stochastic Manna model (SMM) and the Bak-Tang-Weisenfeld (BTW) sandpile model. Employing the approach of block maxima, the study numerically reveals that the distributions for extreme avalanche size and area follow the generalized extreme value (GEV) distribution. The extreme avalanche size follows the Gumbel distribution with shape parameter $ξ=0$ while in the case of the extreme avalanche area, we report $ξ>0$. We propose scaling functions for extreme avalanche activities that connect the activities on different length scales. With the help of data collapse, we estimate the precise values of these critical exponents. The scaling functions provide an understanding of the intricate dynamics for different variants of the sandpile model, shedding light on the relationship between system size and extreme event characteristics. Our findings give insight into the extreme behavior of SOC models and offer a framework to understand the statistical properties of extreme events.
academic

Extreme events scaling in self-organized critical models

基本信息

  • 论文ID: 2510.08733
  • 标题: Extreme events scaling in self-organized critical models
  • 作者: Abdul Quadir, Haider Hasan Jafri (阿里格尔穆斯林大学物理系)
  • 分类: cond-mat.stat-mech (统计力学)
  • 发表时间: 2025年10月9日 (arXiv预印本)
  • 论文链接: https://arxiv.org/abs/2510.08733

摘要

本研究探讨了有限尺寸二维自组织临界(SOC)模型中雪崩活动的极值事件，特别关注随机Manna模型(SMM)和Bak-Tang-Weisenfeld(BTW)沙堆模型。采用块最大值方法，研究数值揭示了极值雪崩尺寸和面积的分布遵循广义极值(GEV)分布。极值雪崩尺寸遵循形状参数ξ=0的Gumbel分布，而极值雪崩面积显示ξ>0。研究提出了连接不同长度尺度活动的极值雪崩活动标度函数，通过数据坍缩估算了临界指数的精确值，为理解SOC模型的极值行为提供了框架。

研究背景与动机

研究问题

本研究要解决的核心问题是：在有限尺寸的自组织临界系统中，极值事件如何随系统尺寸变化，以及这些极值事件的统计性质如何用极值理论来描述和预测。

重要性

  1. 复杂系统普遍性: 自组织临界现象广泛存在于自然和社会系统中，如地震、森林火灾、金融市场崩盘等
  2. 风险评估: 理解极值事件对风险管理、灾难预防和系统韧性评估具有重要意义
  3. 理论完善: 现有SOC理论主要关注无限系统，有限尺寸效应的研究相对不足

现有方法局限性

  • 传统SOC理论假设无限系统尺寸，忽略了实际系统的有限尺寸效应
  • 缺乏系统性的极值事件统计描述框架
  • 对不同SOC模型间极值行为的普适性认识不足

核心贡献

  1. 建立了SOC模型中极值事件的标度理论：提出了连接不同系统尺寸的标度函数
  2. 发现了极值分布的普适性：证明了极值雪崩尺寸和面积分别遵循不同的GEV分布类型
  3. 精确测定了临界指数：通过数据坍缩技术获得了精确的标度指数值
  4. 提供了统一的理论框架：将极值理论与SOC理论相结合，为理解复杂系统极值行为提供新视角

方法详解

任务定义

研究有限尺寸二维SOC模型中雪崩活动的极值统计性质，包括：

  • 输入：不同系统尺寸N=L²的雪崩时间序列
  • 输出：极值分布参数、标度函数和临界指数
  • 约束：二维周期边界条件，开放边界耗散

模型架构

1. SOC模型

随机Manna模型(SMM)：

  • 格点上粒子数hi ≥ 2时发生倾倒
  • 倾倒规则：hi → hi-2，随机选择两个最近邻格点各获得一个粒子
  • 阈值hc = 2

BTW沙堆模型：

  • 格点上粒子数hi ≥ 4时发生倾倒
  • 倾倒规则：hi → hi-4，四个最近邻格点各获得一个粒子
  • 阈值hc = 4，边界耗散

2. 极值分析方法

块最大值技术：

  • 将长度为T的时间序列分成M个块
  • 提取每个块的最大值构成极值序列
  • 分析极值序列的统计性质

广义极值分布： 累积分布函数：

F(x,μ,β,ξ) = exp{-[1 + ξ(x-μ)/β]^(-1/ξ)}

其中μ、β、ξ分别为位置、尺度和形状参数。

技术创新点

  1. 标度变量构造：引入标度变量u = (x-M)/σ，其中M为众数，σ为标准差
  2. 统一标度关系：发现所有统计量遵循相同的标度指数α₁=α₂=α₃=α₄=α
  3. 数据坍缩技术：通过重标度实现不同系统尺寸数据的普适性描述

实验设置

数据集

  • SMM模型：系统尺寸L = 2³, 2⁴, ..., 2⁷
  • BTW模型：相同的系统尺寸范围
  • 时间序列长度：足够长以确保统计收敛
  • 观测量：雪崩尺寸s（倾倒总次数）和雪崩面积a（受影响格点数）

评价指标

  • 拟合优度：R² > 0.99
  • 标度指数：α₁(均值)、α₂(方差)、α₃(众数)、α₄(概率)
  • GEV参数：形状参数ξ、位置参数μ、尺度参数β

实现细节

  • 数值方法：蒙特卡罗模拟
  • 拟合算法：Levenberg-Marquardt算法
  • 计算资源：PARAM Rudra超级计算设施

实验结果

主要结果

标度指数测定

模型极值活动α₁α₂α₃α₄
SMM尺寸1.301.341.331.16
SMM面积1.001.101.011.01
BTW尺寸1.191.361.131.25
BTW面积1.01.100.991.10

GEV分布参数

模型极值活动ξμβ
SMM尺寸0.040.001.00
SMM面积0.210.001.00
BTW尺寸0.003-0.091.00
BTW面积0.153-0.110.99

关键发现

  1. 极值雪崩尺寸：两个模型都显示ξ≈0，遵循Gumbel分布
  2. 极值雪崩面积：显示ξ>0，属于Fréchet类
  3. 标度关系：验证了理论预测α₁=α₂=α₃=α₄
  4. 数据坍缩：成功实现了不同系统尺寸数据的普适性描述

标度函数

概率密度函数的标度形式：

P(x) ~ (1/N^α) F((x-M)/N^α)

其中F为普适标度函数。

相关工作

主要研究方向

  1. SOC理论：Bak-Tang-Wiesenfeld开创性工作建立了SOC概念
  2. 极值理论：Gumbel等建立的经典极值统计理论
  3. 有限尺寸标度：重整化群理论在临界现象中的应用

本文优势

  • 首次系统研究SOC模型中的极值事件标度
  • 结合极值理论与SOC理论的跨学科方法
  • 提供了精确的数值验证和理论预测

结论与讨论

主要结论

  1. 普适性：不同SOC模型的极值事件显示相似的标度行为
  2. 分类：极值雪崩尺寸和面积属于不同的极值分布类别
  3. 标度理论：建立了完整的有限尺寸标度理论框架

局限性

  1. 维度限制：仅考虑二维系统
  2. 模型范围：只研究了两个经典SOC模型
  3. 理论解释：缺乏对不同极值行为的深层物理机制解释

未来方向

  1. 扩展到三维和其他SOC模型
  2. 研究驱动速率对极值行为的影响
  3. 应用于实际地球物理和金融系统

深度评价

优点

  1. 方法创新：巧妙结合极值理论与SOC理论
  2. 数值精度：高质量的数值模拟和精确的参数估计
  3. 理论完整性：建立了从微观动力学到宏观统计的完整框架
  4. 普适性发现：揭示了不同模型间的共同统计规律

不足

  1. 物理解释不足：对为什么雪崩尺寸和面积显示不同极值行为缺乏深入的物理解释
  2. 应用验证缺失：缺乏在实际系统中的验证
  3. 理论推导：主要基于数值观察，理论推导相对薄弱

影响力

  1. 理论贡献：为SOC理论增加了极值统计维度
  2. 应用前景：可应用于风险评估、灾害预测等领域
  3. 跨学科价值：连接了统计物理与极值统计两个领域

适用场景

  • 地震等地球物理现象的极值分析
  • 金融市场极端波动的风险评估
  • 生态系统崩溃等复杂系统极值事件预测
  • 工程系统的可靠性分析

参考文献

本文引用了58篇重要文献，涵盖了SOC理论、极值统计、复杂网络等多个领域的经典和前沿工作，为研究提供了坚实的理论基础。

总体评价：这是一篇高质量的理论物理论文，在SOC理论与极值统计的交叉领域做出了重要贡献。研究方法科学严谨，数值结果可靠，理论框架完整，具有重要的学术价值和应用前景。