We study the scaling properties of avalanche activity in the two-dimensional Abelian sandpile model. Instead of the conventional avalanche size distribution, we analyze the site activity distribution, which measures how often a site participates in avalanches when grains are added across the lattice. Using numerical simulations for system sizes up to \(L = 160\), averaged over \(10^4\) configurations, we determine the probability distribution \(P(A, L)\) of site activities. The results show that \(P(A, L)\) follows a finite-size scaling form \[ P(A, L) \sim L^{-2} F\Big(\frac{A}{L^2}\Big). \] For small values \(A \ll L^2\) the scaling function behaves as \[ F(u) \sim u^{-1/2}, \quad \text{corresponding to} \quad P(A) \sim \frac{1}{L}, \] while for large activities \(A \sim O(L^2)\) the distribution decays as \[ F(u) \sim \exp\big(-c_3 u - c_4 u^2\big). \] The crossover between these two regimes occurs at \[ A^* \sim 0.1 \, L^2, \] marking the threshold between typical and highly excitable sites. This characterization of local avalanche activity provides complementary information to the usual avalanche size statistics, highlighting how local regions serve as frequent conduits for critical dynamics. These results may help connect sandpile models to real-world self-organized critical systems where only partial local activity can be observed.
- ID de l'article: 2510.09631
- Titre: Scaling Properties of Avalanche Activity in the Two-Dimensional Abelian Sandpile Model
- Auteur: Anubhav Ganguly (IISER Mohali)
- Classification: cond-mat.stat-mech nlin.CG
- Date de publication: 28 septembre 2025
- Lien de l'article: https://arxiv.org/abs/2510.09631
Cet article étudie les propriétés d'échelle de l'activité d'avalanche dans le modèle de tas de sable abélien bidimensionnel. Contrairement aux études traditionnelles sur la distribution de la taille des avalanches, l'auteur analyse la distribution d'activité des sites, qui mesure la fréquence de participation des sites aux avalanches lors de l'ajout de particules sur l'ensemble du réseau. Par le biais de simulations numériques pour des tailles de système jusqu'à L=160, moyennées sur 10⁴ configurations, la distribution de probabilité P(A,L) de l'activité des sites a été déterminée. Les résultats montrent que P(A,L) suit la forme d'échelle de taille finie P(A,L)~L⁻²F(A/L²). Pour les petites valeurs A≪L², la fonction d'échelle se comporte comme F(u)~u⁻¹/², correspondant à P(A)1/L ; pour les grandes activités AO(L²), la distribution décroît selon F(u)~exp(-c₃u-c₄u²). La transition entre les deux mécanismes se produit à A*~0.1L², marquant un seuil entre les sites typiques et les sites hautement excitables.
Le problème fondamental que cette recherche vise à résoudre est la compréhension des propriétés statistiques de l'activité locale des sites dans le modèle de tas de sable abélien (ASM). Les études traditionnelles se concentrent principalement sur la distribution de la taille des avalanches s(R) (taille totale de l'avalanche initiée au site R), tandis que cet article propose d'étudier la fonction d'activité A(R) (nombre total de fois où le site R participe à toutes les avalanches possibles).
- Compréhension approfondie de la criticité auto-organisée (SOC): L'ASM est un modèle classique de SOC, et la compréhension de ses modèles d'activité locale est cruciale pour comprendre les phénomènes critiques
- Valeur d'application pratique: Dans les systèmes réels, on ne peut généralement observer que l'activité locale, et non les statistiques globales des avalanches
- Complément théorique: Fournit des informations complémentaires aux statistiques traditionnelles de taille d'avalanche, révélant quelles régions servent de canaux de propagation fréquents pour les avalanches
- Les études traditionnelles se concentrent principalement sur la distribution en loi de puissance Prob(s)~s⁻τ de la taille des avalanches s(R)
- Absence d'analyse systématique de la fréquence de participation des sites locaux aux avalanches
- La distribution d'activité Prob(A,L) n'a pas été systématiquement étudiée dans la littérature
- Première étude systématique de la distribution d'activité des sites: Définition et analyse de la distribution de probabilité de la fonction d'activité A(R)=∑_{R'}N(R,R')
- Découverte d'une nouvelle forme d'échelle: Établissement de la relation d'échelle de taille finie P(A,L)~L⁻²F(A/L²)
- Identification d'un double mécanisme d'échelle: Découverte de F(u)~u⁻¹/² pour la région de petite activité et de la décroissance exponentielle F(u)~exp(-c₃u-c₄u²) pour la région de grande activité
- Détermination de l'échelle de transition: Identification de A*~0.1L² comme point de démarcation entre les deux mécanismes
- Fourniture d'un cadre théorique: Dérivation des exposants d'échelle α=β=2 par le biais de contraintes de normalisation et de valeur attendue
Étant donné une configuration stable C du modèle de tas de sable abélien L×L bidimensionnel, on définit:
- Fonction d'excitation: s(R) = taille de l'avalanche déclenchée par l'ajout d'une particule au site R
- Fonction d'activité: A(R) = ∑_{R'} N(R,R'), où N(R,R') est le nombre d'effondrements du site R lors de l'ajout d'une particule au site R'
L'objectif de la recherche est d'analyser la distribution de probabilité P(A,L) de la fonction d'activité A(R) et ses propriétés d'échelle.
- Configuration du système: Réseau carré L×L avec conditions aux limites ouvertes
- Variables d'état: Chaque site (x,y) possède une hauteur Δ_{x,y}∈{0,1,2,3} (état stable), Δ_{x,y}≥4 étant instable
- Règles de dynamique:
- Ajout de particule: Δ_{x,y} → Δ_{x,y} + 1
- Règle d'effondrement: Δ_{x,y} → Δ_{x,y} - 4, chaque voisin reçoit une particule
- Perte de particules aux limites
- Fixation d'une configuration stable C
- Pour chaque site R':
- Ajout d'une particule à R'
- Enregistrement du nombre d'effondrements N(R,R') du site R
- Réinitialisation à la configuration C
- Calcul de A(R) = ∑_{R'} N(R,R')
- Statistique de la distribution A(R) pour tous les sites
En supposant une forme d'échelle générale P(A,L)~L⁻ᵅF(a/Lᵝ), en utilisant:
- Condition de normalisation: ∑_a P(A,L) = 1 → α = β
- Contrainte de valeur attendue: ⟨A⟩~L² → 2β - α = 2
On obtient α = β = 2, c'est-à-dire P(A,L)~L⁻²F(A/L²).
- Région de petite activité (u≪0.1): F(u)~u⁻¹/² → P(A)~1/L
- Région de grande activité (u≳0.1): F(u)~exp(-c₃u-c₄u²)
- Forme unifiée: F(u) = (c₁u⁻⁰·⁵ + c₂)e⁻ᶜ³ᵘ⁻ᶜ⁴ᵘ²
- Tailles de système: L = 20 à 160
- Conditions aux limites: Limites ouvertes
- Échantillons statistiques: 10 000 configurations par taille
- Ressources de calcul: Cluster MARIO
- Le système atteint l'état critique
- Pour chaque configuration C, calcul de A(R) pour tous les sites
- Statistique de la distribution de probabilité empirique P(A=a,L)
- Analyse d'échelle multi-taille
- Exposants d'échelle: Détermination des exposants en loi de puissance par ajustement sur graphique log-log
- Qualité de l'effondrement des données: R² = 0.9839
- Vérification de la valeur attendue: ⟨A⟩ ≈ 0.085L² en accord avec les prédictions théoriques
- Pour une petite valeur fixe a, P(A=a,L)~1/L, avec une pente de -1±0.015
- Fonction d'échelle g(a) = LP(A,L) ~ a⁻⁰·³⁸⁵±⁰·⁰³
- Peut être ajustée sous la forme g(x) ≈ c₁x⁻¹/² + c₂
- La distribution complète obéit à P(A,L)~L⁻²F(A/L²)
- Fonction d'échelle F(u) = (c₁u⁻⁰·⁵ + c₂)e⁻ᶜ³ᵘ⁻ᶜ⁴ᵘ²
- Estimation des paramètres: c₃~10, c₄=O(1)
- Échelle de transition A*~0.1L²
- En dessous de A*: comportement en loi de puissance dominant
- Au-dessus de A*: décroissance exponentielle dominante
Pour les configurations rares C_ε avec masse dépassant la configuration standard de εL²:
- ⟨A⟩_C ~ εL³ (terme dominant)
- Fraction de ces configurations ~exp(-ε²L²)
- Conduisant à F(x~εL) ≳ exp(-x²) pour grand x
Les figures 6 et 9 montrent un bon effondrement des données de plusieurs tailles de système sous la forme d'échelle, vérifiant l'exactitude des prédictions théoriques.
- Modèle Bak-Tang-Wiesenfeld 1,2: Travail fondateur de la SOC
- Correspondance Majumdar-Dhar 9: Équivalence entre l'ASM et le modèle de Potts q→0
- Description par théorie des champs conformes: Charge centrale c=-2, indice de dynamique z=5/4
- Recherches traditionnelles: Concentration sur la loi de puissance Prob(s)~s⁻τ de la distribution de taille d'avalanche s(R)
- Modèles déterministes vs aléatoires: Le modèle de Manna et le modèle BTW appartiennent à des classes d'universalité différentes 12-18
- Analyse multifractale 17,18: Calcul précis des exposants critiques
- Innovation de l'observable: Passage de l'étude de s(R) à celle de A(R)
- Théorie d'échelle: Établissement d'une théorie d'échelle complète pour la distribution d'activité
- Valeur pratique: Fourniture d'une base théorique pour les systèmes à observation partielle
- Établissement de la forme d'échelle: P(A,L)~L⁻²F(A/L²) comme forme universelle de la distribution d'activité
- Identification du double mécanisme: Effets de limite pour petite activité (~u⁻¹/²) et effets de volume pour grande activité (décroissance exponentielle)
- Échelle de transition: A*~0.1L² marquant la transition du régime dominé par les limites au régime dominé par le volume
- Cohérence théorique: Satisfaction des contraintes de normalisation et de valeur attendue
- Local vs global: A(R) révèle les modèles de participation locale, complétant les statistiques globales d'avalanche
- Identification des points chauds: Les régions de haute activité servent de canaux de propagation fréquents pour les avalanches
- Propriétés critiques: Les effets de limite produisent des contributions d'échelle libre, mais la moyenne du système est dominée par les grandes avalanches
- Effets de taille finie: La taille maximale du système L=160 peut être insuffisante pour éliminer complètement les effets de taille finie
- Dépendance des conditions aux limites: Le choix des conditions aux limites ouvertes peut affecter le comportement de petite activité
- Précision numérique: L'écart significatif à A=2 reflète les contraintes combinatoires des séquences d'effondrement
- Systèmes plus grands: Extension à des tailles plus grandes pour vérifier le comportement asymptotique
- Conditions aux limites différentes: Étude des différences de comportement sous conditions aux limites périodiques
- Vérification expérimentale: Application de la théorie aux observations locales dans les systèmes SOC réels
- Théorie analytique: Développement d'une théorie analytique rigoureuse pour la distribution d'activité
- Innovation conceptuelle: La perspective d'étude de l'activité des sites plutôt que de la taille traditionnelle des avalanches a une valeur importante
- Complétude théorique: Dérivation rigoureuse des relations d'échelle par le biais de contraintes de normalisation et de valeur attendue
- Vérification numérique suffisante: Bon effondrement des données multi-taille vérifiant les prédictions théoriques
- Valeur pratique élevée: Fourniture d'une base théorique pour les systèmes à observation partielle
- Clarté de la rédaction: Structure logique claire et dérivations mathématiques rigoureuses
- Limitation de la taille du système: La taille maximale L=160 peut être insuffisante pour vérifier complètement le comportement asymptotique
- Traitement des effets de limite: La séparation des effets de limite et de volume nécessite une analyse plus approfondie
- Manque de théorie analytique: Dépendance principalement sur les résultats numériques, absence de dérivation analytique rigoureuse
- Vérification d'universalité: Nécessité de vérifier les comportements similaires dans d'autres modèles SOC
- Contribution théorique: Fourniture d'une nouvelle observable et d'un cadre d'analyse pour la théorie SOC
- Valeur méthodologique: Les méthodes d'analyse d'échelle peuvent être généralisées à d'autres systèmes complexes
- Perspectives d'application: Valeur d'application potentielle en neurosciences, sismologie et autres domaines
- Reproductibilité: Description détaillée des méthodes facilitant la reproduction et l'extension
- Réseaux de neurones: Analyse des modèles d'activité locale des neurones
- Systèmes sismiques: Étude des statistiques d'activité des capteurs de déformation locale
- Marchés financiers: Analyse des modèles de participation aux fluctuations du marché local
- Écosystèmes: Étude de la distribution du degré de participation des espèces dans les réseaux écologiques
Cet article cite les travaux importants du domaine de la SOC, notamment:
- Les travaux fondateurs de Bak, Tang, Wiesenfeld 1,2
- Les contributions théoriques de Dhar 4,5
- Les résultats exacts de Majumdar-Dhar 9
- Les progrès expérimentaux récents 22,23
Évaluation Générale: Cet article constitue un travail innovant dans l'étude du modèle de tas de sable abélien, fournissant une nouvelle perspective à la théorie SOC en introduisant le concept de distribution d'activité des sites. Bien qu'il présente certaines limitations en termes de taille de système et de théorie analytique, son cadre théorique est complet et sa vérification numérique est suffisante, offrant une valeur importante pour la compréhension des phénomènes critiques locaux.