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

自己組織化臨界モデルにおける極値事象のスケーリング

基本情報

  • 論文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、ランダムに選択された2つの最近接格子点が各々1つの粒子を獲得
  • 閾値hc = 2

BTW砂山モデル：

  • 格子点上の粒子数hi ≥ 4の場合に転倒が発生
  • 転倒ルール：hi → hi-4、4つの最近接格子点が各々1つの粒子を獲得
  • 閾値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. モデル範囲: 2つの古典的SOCモデルのみを研究しました
  3. 理論的解釈: 異なる極値挙動の深層的な物理メカニズムに関する説明が不足しています

今後の方向性

  1. 三次元および他のSOCモデルへの拡張
  2. 駆動速度が極値挙動に与える影響の研究
  3. 実際の地球物理および金融システムへの応用

深度評価

利点

  1. 方法論の革新: 極値理論とSOC理論を巧妙に結合しました
  2. 数値精度: 高品質の数値シミュレーションと正確なパラメータ推定
  3. 理論的完全性: 微視的動力学から巨視的統計への完全な枠組みを確立しました
  4. 普遍性の発見: 異なるモデル間の共通統計則を明らかにしました

不足点

  1. 物理的解釈の不足: 雪崩サイズと面積が異なる極値挙動を示す理由に関する深い物理的解釈が欠けています
  2. 応用検証の欠如: 実際のシステムにおける検証が不足しています
  3. 理論的導出: 主に数値観察に基づいており、理論的導出は相対的に弱いです

影響力

  1. 理論的貢献: SOC理論に極値統計の次元を追加しました
  2. 応用の見通し: リスク評価、災害予測などの分野に応用可能です
  3. 学際的価値: 統計物理と極値統計の2つの分野を結びつけました

適用シーン

  • 地震などの地球物理現象の極値分析
  • 金融市場の極端な変動のリスク評価
  • 生態系の崩壊など複雑系の極値事象予測
  • 工学システムの信頼性分析

参考文献

本論文は58篇の重要な文献を引用しており、SOC理論、極値統計、複雑ネットワークなど複数の分野の古典的および最先端の研究をカバーし、研究に堅実な理論的基礎を提供しています。

総合評価: これは統計物理学の高品質な理論論文であり、SOC理論と極値統計の交差領域で重要な貢献をしています。研究方法は科学的で厳密であり、数値結果は信頼でき、理論枠組みは完全であり、重要な学術的価値と応用の見通しを有しています。