Extreme events scaling in self-organized critical models

Basic Information

• Paper ID: 2510.08733
• Title: Extreme events scaling in self-organized critical models
• Authors: Abdul Quadir, Haider Hasan Jafri (Department of Physics, Aligarh Muslim University)
• Classification: cond-mat.stat-mech (Statistical Mechanics)
• Publication Date: October 9, 2025 (arXiv preprint)
• Paper Link: https://arxiv.org/abs/2510.08733

Abstract

This study investigates extreme value events in avalanche activity within finite-size two-dimensional self-organized critical (SOC) models, with particular focus on the stochastic Manna model (SMM) and the Bak-Tang-Wiesenfeld (BTW) sandpile model. Employing the block maxima method, the numerical study reveals that the distributions of extreme avalanche sizes and areas follow the generalized extreme value (GEV) distribution. Extreme avalanche sizes follow a Gumbel distribution with shape parameter ξ=0, while extreme avalanche areas exhibit ξ>0. The study proposes scaling functions for extreme avalanche activity connecting different length scales, and estimates precise values of critical exponents through data collapse, providing a framework for understanding the extreme value behavior of SOC models.

Research Background and Motivation

Research Questions

The core questions addressed in this study are: How do extreme events scale with system size in finite-size self-organized critical systems, and how can the statistical properties of these extreme events be described and predicted using extreme value theory?

Significance

1. Universality in Complex Systems: Self-organized criticality phenomena are ubiquitous in natural and social systems, such as earthquakes, forest fires, and financial market crashes
2. Risk Assessment: Understanding extreme events is crucial for risk management, disaster prevention, and system resilience evaluation
3. Theoretical Refinement: Existing SOC theory primarily focuses on infinite systems, with relatively limited research on finite-size effects

Limitations of Existing Approaches

• Traditional SOC theory assumes infinite system size, neglecting finite-size effects in realistic systems
• Lack of systematic statistical description framework for extreme events
• Insufficient understanding of universality in extreme value behavior across different SOC models

Core Contributions

1. Established scaling theory for extreme events in SOC models: Proposed scaling functions connecting different system sizes
2. Discovered universality in extreme value distributions: Demonstrated that extreme avalanche sizes and areas follow different GEV distribution types
3. Precisely measured critical exponents: Obtained accurate scaling exponent values through data collapse techniques
4. Provided a unified theoretical framework: Combined extreme value theory with SOC theory, offering new perspectives for understanding extreme value behavior in complex systems

Methodology Details

Task Definition

Investigation of extreme value statistical properties of avalanche activity in finite-size two-dimensional SOC models, including:

• Input: Avalanche time series for different system sizes N=L²
• Output: Extreme value distribution parameters, scaling functions, and critical exponents
• Constraints: Two-dimensional periodic boundary conditions, open boundary dissipation

Model Architecture

1. SOC Models

Stochastic Manna Model (SMM):

• Toppling occurs when particle number hi ≥ 2
• Toppling rule: hi → hi-2, with one particle randomly given to each of two nearest neighbors
• Threshold hc = 2

BTW Sandpile Model:

• Toppling occurs when particle number hi ≥ 4
• Toppling rule: hi → hi-4, with one particle given to each of four nearest neighbors
• Threshold hc = 4, boundary dissipation

2. Extreme Value Analysis Methods

Block Maxima Technique:

• Divide time series of length T into M blocks
• Extract maximum value from each block to construct extreme value sequence
• Analyze statistical properties of extreme value sequence

Generalized Extreme Value Distribution: Cumulative distribution function:

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

where μ, β, ξ are location, scale, and shape parameters, respectively.

Technical Innovations

1. Scaling variable construction: Introduce scaling variable u = (x-M)/σ, where M is the mode and σ is the standard deviation
2. Unified scaling relations: Discovered that all statistical quantities follow the same scaling exponent α₁=α₂=α₃=α₄=α
3. Data collapse technique: Achieved universal description of data from different system sizes through rescaling

Experimental Setup

Datasets

• SMM model: System sizes L = 2³, 2⁴, ..., 2⁷
• BTW model: Same system size range
• Time series length: Sufficiently long to ensure statistical convergence
• Observables: Avalanche size s (total number of topplings) and avalanche area a (number of affected sites)

Evaluation Metrics

• Goodness of fit: R² > 0.99
• Scaling exponents: α₁ (mean), α₂ (variance), α₃ (mode), α₄ (probability)
• GEV parameters: Shape parameter ξ, location parameter μ, scale parameter β

Implementation Details

• Numerical method: Monte Carlo simulation
• Fitting algorithm: Levenberg-Marquardt algorithm
• Computational resources: PARAM Rudra supercomputing facility

Experimental Results

Main Results

Scaling Exponent Measurements

GEV Distribution Parameters

Key Findings

1. Extreme avalanche sizes: Both models show ξ≈0, following Gumbel distribution
2. Extreme avalanche areas: Show ξ>0, belonging to Fréchet class
3. Scaling relations: Verified theoretical prediction α₁=α₂=α₃=α₄
4. Data collapse: Successfully achieved universal description of data from different system sizes

Scaling Functions

Scaling form of probability density function:

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

where F is the universal scaling function.

Related Work

Main Research Directions

1. SOC Theory: Foundational work by Bak-Tang-Wiesenfeld establishing the SOC concept
2. Extreme Value Theory: Classical extreme value statistics theory established by Gumbel and others
3. Finite-Size Scaling: Application of renormalization group theory in critical phenomena

Advantages of This Work

• First systematic study of extreme event scaling in SOC models
• Interdisciplinary approach combining extreme value theory with SOC theory
• Provides precise numerical verification and theoretical predictions

Conclusions and Discussion

Main Conclusions

1. Universality: Extreme events in different SOC models display similar scaling behavior
2. Classification: Extreme avalanche sizes and areas belong to different extreme value distribution categories
3. Scaling Theory: Established a complete finite-size scaling theory framework

Limitations

1. Dimensional restriction: Only considers two-dimensional systems
2. Model scope: Only studies two classical SOC models
3. Theoretical explanation: Lacks deep physical mechanism explanation for different extreme value behaviors

Future Directions

1. Extension to three-dimensional and other SOC models
2. Investigation of driving rate effects on extreme value behavior
3. Application to real geophysical and financial systems

In-Depth Evaluation

Strengths

1. Methodological innovation: Skillfully combines extreme value theory with SOC theory
2. Numerical precision: High-quality numerical simulations and accurate parameter estimation
3. Theoretical completeness: Establishes a complete framework from microscopic dynamics to macroscopic statistics
4. Universality discovery: Reveals common statistical laws across different models

Weaknesses

1. Insufficient physical interpretation: Lacks deep physical explanation for why avalanche sizes and areas display different extreme value behaviors
2. Missing experimental validation: Lacks verification in real systems
3. Weak theoretical derivation: Primarily based on numerical observations with relatively limited theoretical derivation

Impact

1. Theoretical contribution: Adds extreme value statistics dimension to SOC theory
2. Application prospects: Applicable to risk assessment, disaster prediction, and other fields
3. Interdisciplinary value: Bridges statistical physics and extreme value statistics

Applicable Scenarios

• Extreme value analysis of geophysical phenomena such as earthquakes
• Risk assessment of extreme financial market fluctuations
• Prediction of extreme events in complex systems such as ecosystem collapse
• Reliability analysis of engineering systems

References

This paper cites 58 important references covering classical and cutting-edge works in SOC theory, extreme value statistics, complex networks, and other fields, providing a solid theoretical foundation for the research.

Overall Assessment: This is a high-quality theoretical physics paper that makes important contributions at the intersection of SOC theory and extreme value statistics. The research methodology is scientifically rigorous, numerical results are reliable, the theoretical framework is complete, and it possesses significant academic value and application prospects.