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: 2510.09631
- शीर्षक: द्वि-आयामी एबेलियन सैंडपाइल मॉडल में हिमस्खलन गतिविधि के स्केलिंग गुण
- लेखक: अनुभव गांगुली (IISER मोहाली)
- वर्गीकरण: cond-mat.stat-mech nlin.CG
- प्रकाशन तिथि: 28 सितंबर 2025
- पेपर लिंक: https://arxiv.org/abs/2510.09631
यह पेपर द्वि-आयामी एबेलियन सैंडपाइल मॉडल में हिमस्खलन गतिविधि के स्केलिंग गुणों का अध्ययन करता है। परंपरागत हिमस्खलन आकार वितरण अनुसंधान के विपरीत, लेखक स्थल गतिविधि वितरण का विश्लेषण करते हैं, जो मापता है कि जब कणों को पूरी जाली में जोड़ा जाता है तो स्थल कितनी बार हिमस्खलन में भाग लेते हैं। L=160 तक की प्रणाली के आकार के साथ संख्यात्मक सिमुलेशन के माध्यम से, 10⁴ कॉन्फ़िगरेशन पर औसत, स्थल गतिविधि की संभाव्यता वितरण P(A,L) निर्धारित की गई है। परिणाम दिखाते हैं कि P(A,L) परिमित आकार स्केलिंग रूप P(A,L)~L⁻²F(A/L²) का पालन करता है। छोटे मानों A≪L² के लिए, स्केलिंग फ़ंक्शन F(u)~u⁻¹/² के रूप में व्यवहार करता है, जो P(A)1/L के अनुरूप है; बड़ी गतिविधि AO(L²) के लिए, वितरण F(u)~exp(-c₃u-c₄u²) में क्षय होता है। दोनों तंत्रों के बीच क्रॉसओवर A*~0.1L² पर होता है, जो विशिष्ट स्थलों और उच्च-उत्तेजनीय स्थलों के बीच सीमा को चिह्नित करता है।
यह अनुसंधान एबेलियन सैंडपाइल मॉडल (ASM) में स्थानीय स्थल गतिविधि के सांख्यिकीय गुणों को समझने की मूल समस्या को हल करने का प्रयास करता है। परंपरागत अनुसंधान मुख्य रूप से हिमस्खलन आकार वितरण s(R) (स्थल R से शुरू होने वाले हिमस्खलन का कुल आकार) पर केंद्रित है, जबकि यह पेपर गतिविधि फ़ंक्शन A(R) (स्थल R सभी संभावित हिमस्खलन में भाग लेने की कुल संख्या) का अध्ययन करने का प्रस्ताव देता है।
- आत्म-संगठित आलोचनात्मकता (SOC) की गहन समझ: ASM, SOC का एक शास्त्रीय मॉडल है, इसकी स्थानीय गतिविधि पैटर्न को समझना आलोचनात्मक घटना को समझने के लिए महत्वपूर्ण है
- व्यावहारिक अनुप्रयोग मूल्य: वास्तविक प्रणालियों में अक्सर केवल स्थानीय गतिविधि का अवलोकन किया जा सकता है, न कि वैश्विक हिमस्खलन सांख्यिकी
- सैद्धांतिक पूरक: परंपरागत हिमस्खलन आकार सांख्यिकी के साथ पूरक जानकारी प्रदान करता है, यह प्रकट करता है कि कौन से क्षेत्र हिमस्खलन प्रसार के लिए आवृत्ति चैनल के रूप में कार्य करते हैं
- परंपरागत अनुसंधान मुख्य रूप से हिमस्खलन आकार s(R) के शक्ति-नियम वितरण Prob(s)~s⁻τ पर केंद्रित है
- स्थानीय स्थलों के हिमस्खलन में भाग लेने की आवृत्ति का व्यवस्थित विश्लेषण अभाव है
- गतिविधि वितरण Prob(A,L) साहित्य में व्यवस्थित रूप से अध्ययन नहीं किया गया है
- स्थल गतिविधि वितरण का पहला व्यवस्थित अध्ययन: गतिविधि फ़ंक्शन A(R)=∑_{R'}N(R,R') की संभाव्यता वितरण को परिभाषित और विश्लेषण करता है
- नया स्केलिंग रूप खोजना: P(A,L)~L⁻²F(A/L²) का परिमित आकार स्केलिंग संबंध स्थापित करता है
- दोहरी स्केलिंग तंत्र की पहचान: छोटी गतिविधि क्षेत्र F(u)~u⁻¹/² और बड़ी गतिविधि क्षेत्र के घातीय क्षय F(u)~exp(-c₃u-c₄u²) की खोज करता है
- क्रॉसओवर स्केल निर्धारण: A*~0.1L² को दोनों तंत्रों के विभाजन बिंदु के रूप में पहचानता है
- सैद्धांतिक ढांचा प्रदान करना: सामान्यीकरण और अपेक्षा मूल्य बाधाओं के माध्यम से स्केलिंग घातांक α=β=2 प्राप्त करता है
द्वि-आयामी L×L एबेलियन सैंडपाइल मॉडल के स्थिर कॉन्फ़िगरेशन C को देखते हुए, परिभाषित करें:
- उत्तेजन फ़ंक्शन: s(R) = स्थल R पर कण जोड़ने से ट्रिगर किए गए हिमस्खलन का आकार
- गतिविधि फ़ंक्शन: A(R) = ∑_{R'} N(R,R'), जहां N(R,R') स्थल R' पर कण जोड़ते समय स्थल R के पतन की संख्या है
अनुसंधान का लक्ष्य गतिविधि फ़ंक्शन A(R) की संभाव्यता वितरण P(A,L) और इसके स्केलिंग गुणों का विश्लेषण करना है।
- प्रणाली सेटअप: L×L वर्गाकार जाली, खुली सीमा शर्तें
- स्थिति चर: प्रत्येक स्थल (x,y) की ऊंचाई Δ_{x,y}∈{0,1,2,3} (स्थिर अवस्था), Δ_{x,y}≥4 अस्थिर
- गतिविधि नियम:
- कण जोड़ना: Δ_{x,y} → Δ_{x,y} + 1
- पतन नियम: Δ_{x,y} → Δ_{x,y} - 4, चार पड़ोसियों को प्रत्येक एक कण मिलता है
- सीमा कणों का प्रवाह
- स्थिर कॉन्फ़िगरेशन C को ठीक करें
- प्रत्येक स्थल R' के लिए:
- R' पर कण जोड़ें
- स्थल R के पतन की संख्या N(R,R') रिकॉर्ड करें
- कॉन्फ़िगरेशन C पर रीसेट करें
- A(R) = ∑_{R'} N(R,R') की गणना करें
- सभी स्थलों के A(R) वितरण को सांख्यिकीय करें
सामान्य स्केलिंग रूप P(A,L)~L⁻ᵅF(a/Lᵝ) मानते हुए, उपयोग करें:
- सामान्यीकरण शर्त: ∑_a P(A,L) = 1 → α = β
- अपेक्षा मूल्य बाधा: ⟨A⟩~L² → 2β - α = 2
प्राप्त करें α = β = 2, अर्थात् P(A,L)~L⁻²F(A/L²)।
- छोटी गतिविधि क्षेत्र (u≪0.1): F(u)~u⁻¹/² → P(A)~1/L
- बड़ी गतिविधि क्षेत्र (u≳0.1): F(u)~exp(-c₃u-c₄u²)
- एकीकृत रूप: F(u) = (c₁u⁻⁰·⁵ + c₂)e⁻ᶜ³ᵘ⁻ᶜ⁴ᵘ²
- प्रणाली आकार: L = 20 से 160
- सीमा शर्तें: खुली सीमा
- सांख्यिकीय नमूने: प्रत्येक आकार के लिए 10,000 कॉन्फ़िगरेशन
- कम्प्यूटेशनल संसाधन: MARIO क्लस्टर
- प्रणाली आलोचनात्मक अवस्था तक पहुंचती है
- प्रत्येक कॉन्फ़िगरेशन C के लिए सभी स्थलों के A(R) की गणना करें
- अनुभवजन्य संभाव्यता वितरण P(A=a,L) को सांख्यिकीय करें
- बहु-आकार स्केलिंग विश्लेषण करें
- स्केलिंग घातांक: log-log ग्राफ फिटिंग के माध्यम से शक्ति-नियम घातांक निर्धारित करें
- डेटा पतन गुणवत्ता: R² = 0.9839
- अपेक्षा मूल्य सत्यापन: ⟨A⟩ ≈ 0.085L² सैद्धांतिक अपेक्षा के साथ सुसंगत
- निश्चित छोटे मान a के लिए, P(A=a,L)~1/L, ढलान -1±0.015
- स्केलिंग फ़ंक्शन g(a) = LP(A,L) ~ a⁻⁰·³⁸⁵±⁰·⁰³
- g(x) ≈ c₁x⁻¹/² + c₂ रूप में फिट किया जा सकता है
- पूर्ण वितरण P(A,L)~L⁻²F(A/L²) का पालन करता है
- स्केलिंग फ़ंक्शन F(u) = (c₁u⁻⁰·⁵ + c₂)e⁻ᶜ³ᵘ⁻ᶜ⁴ᵘ²
- पैरामीटर अनुमान: c₃~10, c₄=O(1)
- क्रॉसओवर स्केल A*~0.1L²
- A* से कम: शक्ति-नियम व्यवहार प्रभावी
- A* से अधिक: घातीय क्षय प्रभावी
मानक कॉन्फ़िगरेशन से अधिक द्रव्यमान εL² वाले दुर्लभ कॉन्फ़िगरेशन C_ε के लिए:
- ⟨A⟩_C ~ εL³ (प्रमुख पद)
- इस प्रकार के कॉन्फ़िगरेशन का अंश ~exp(-ε²L²)
- F(x~εL) ≳ exp(-x²) बड़े x के लिए
चित्र 6 और 9 स्केलिंग रूप के तहत कई प्रणाली आकारों के डेटा के अच्छे पतन को दिखाते हैं, जो सैद्धांतिक भविष्यवाणी की सही्ता को सत्यापित करते हैं।
- Bak-Tang-Wiesenfeld मॉडल1,2: SOC का अग्रणी कार्य
- Majumdar-Dhar मानचित्रण9: ASM और q→0 Potts मॉडल की समानता
- अनुरूप क्षेत्र सिद्धांत विवरण: केंद्रीय आवेश c=-2, गतिविधि सूचकांक z=5/4
- परंपरागत अनुसंधान: हिमस्खलन आकार वितरण s(R) के शक्ति-नियम Prob(s)~s⁻τ पर केंद्रित
- निर्धारणीय बनाम यादृच्छिक मॉडल: Manna मॉडल और BTW मॉडल विभिन्न सार्वभौमिकता वर्गों से संबंधित हैं12-18
- बहु-भग्न विश्लेषण17,18: आलोचनात्मक घातांकों की सटीक गणना
- अवलोकन मात्रा नवाचार: s(R) से A(R) अनुसंधान की ओर स्थानांतरण
- स्केलिंग सिद्धांत: गतिविधि वितरण का पूर्ण स्केलिंग सिद्धांत स्थापित करता है
- व्यावहारिक मूल्य: आंशिक अवलोकन प्रणालियों के लिए सैद्धांतिक आधार प्रदान करता है
- स्केलिंग रूप स्थापना: P(A,L)~L⁻²F(A/L²) गतिविधि वितरण का सार्वभौमिक रूप है
- दोहरी तंत्र पहचान: छोटी गतिविधि के सीमा प्रभाव (~u⁻¹/²) और बड़ी गतिविधि के आयतन प्रभाव (घातीय क्षय)
- क्रॉसओवर स्केल: A*~0.1L² सीमा-प्रभावी से आयतन-प्रभावी में परिवर्तन को चिह्नित करता है
- सैद्धांतिक सुसंगतता: सामान्यीकरण और अपेक्षा मूल्य बाधाओं को संतुष्ट करता है
- स्थानीय बनाम वैश्विक: A(R) स्थानीय भाग लेने के पैटर्न को प्रकट करता है, वैश्विक हिमस्खलन सांख्यिकी को पूरक करता है
- हॉटस्पॉट पहचान: उच्च गतिविधि क्षेत्र हिमस्खलन प्रसार के लिए आवृत्ति चैनल के रूप में कार्य करते हैं
- आलोचनात्मक गुण: सीमा प्रभाव स्केलिंग-मुक्त योगदान उत्पन्न करते हैं, लेकिन प्रणाली औसत बड़े हिमस्खलन द्वारा प्रभावित होता है
- परिमित आकार प्रभाव: अधिकतम प्रणाली L=160 परिमित आकार प्रभाव को पूरी तरह से समाप्त करने के लिए अपर्याप्त हो सकता है
- सीमा शर्त निर्भरता: खुली सीमा शर्त की पसंद छोटी गतिविधि व्यवहार को प्रभावित कर सकती है
- संख्यात्मक सटीकता: A=2 पर उल्लेखनीय विचलन पतन अनुक्रम की संयोजक बाधाओं को प्रतिबिंबित करता है
- बड़ी प्रणालियां: अनंतस्पर्शी व्यवहार को सत्यापित करने के लिए बड़े आकार तक विस्तार करें
- विभिन्न सीमा शर्तें: आवधिक सीमा शर्तों के तहत व्यवहार अंतर का अध्ययन करें
- प्रायोगिक सत्यापन: सिद्धांत को वास्तविक SOC प्रणालियों के स्थानीय अवलोकन पर लागू करें
- विश्लेषणात्मक सिद्धांत: गतिविधि वितरण का कठोर विश्लेषणात्मक सिद्धांत विकसित करें
- अवधारणा नवाचार: परंपरागत हिमस्खलन आकार से स्थल गतिविधि अनुसंधान में स्थानांतरण महत्वपूर्ण मूल्य रखता है
- सैद्धांतिक पूर्णता: सामान्यीकरण और अपेक्षा मूल्य बाधाओं के माध्यम से स्केलिंग संबंध को कठोरता से प्राप्त करता है
- पर्याप्त संख्यात्मक सत्यापन: बहु-आकार डेटा का अच्छा पतन सैद्धांतिक भविष्यवाणी को सत्यापित करता है
- उच्च व्यावहारिक मूल्य: आंशिक अवलोकन प्रणालियों के लिए सैद्धांतिक आधार प्रदान करता है
- स्पष्ट लेखन: तार्किक संरचना स्पष्ट, गणितीय व्युत्पत्ति कठोर
- प्रणाली आकार सीमा: L=160 का अधिकतम आकार अनंतस्पर्शी व्यवहार को पूरी तरह से सत्यापित करने के लिए अपर्याप्त हो सकता है
- सीमा प्रभाव उपचार: सीमा प्रभाव और आयतन प्रभाव के पृथक्करण के लिए गहन विश्लेषण की आवश्यकता है
- विश्लेषणात्मक सिद्धांत की कमी: मुख्य रूप से संख्यात्मक परिणामों पर निर्भर, कठोर विश्लेषणात्मक व्युत्पत्ति की कमी
- सार्वभौमिकता सत्यापन: अन्य SOC मॉडलों में समान व्यवहार को सत्यापित करने की आवश्यकता है
- सैद्धांतिक योगदान: SOC सिद्धांत के लिए नई अवलोकन मात्रा और विश्लेषण ढांचा प्रदान करता है
- पद्धति मूल्य: स्केलिंग विश्लेषण विधि अन्य जटिल प्रणालियों में सामान्यीकृत की जा सकती है
- अनुप्रयोग संभावनाएं: तंत्रिका विज्ञान, भूकंप विज्ञान आदि क्षेत्रों में संभावित अनुप्रयोग मूल्य
- पुनरुत्पादनीयता: विधि विवरण विस्तृत, पुनरुत्पादन और विस्तार के लिए सुविधाजनक
- तंत्रिका नेटवर्क: तंत्रिका कोशिका निर्वहन के स्थानीय गतिविधि पैटर्न का विश्लेषण करें
- भूकंप प्रणाली: स्थानीय तनाव सेंसर की गतिविधि सांख्यिकी का अध्ययन करें
- वित्तीय बाजार: स्थानीय बाजार के उतार-चढ़ाव भाग लेने के पैटर्न का विश्लेषण करें
- पारिस्थितिक प्रणाली: पारिस्थितिक नेटवर्क में प्रजातियों की भाग लेने की डिग्री वितरण का अध्ययन करें
यह पेपर SOC क्षेत्र के महत्वपूर्ण साहित्य का हवाला देता है, जिसमें शामिल हैं:
- Bak, Tang, Wiesenfeld का अग्रणी कार्य1,2
- Dhar का सैद्धांतिक योगदान4,5
- Majumdar-Dhar के सटीक परिणाम9
- हाल के प्रायोगिक प्रगति22,23
समग्र मूल्यांकन: यह एबेलियन सैंडपाइल मॉडल अनुसंधान में एक नवीन कार्य है, जो स्थल गतिविधि वितरण की अवधारणा को प्रस्तुत करके SOC सिद्धांत के लिए एक नया दृष्टिकोण प्रदान करता है। यद्यपि प्रणाली आकार और विश्लेषणात्मक सिद्धांत के संदर्भ में कुछ सीमाएं हैं, लेकिन इसका सैद्धांतिक ढांचा पूर्ण है, संख्यात्मक सत्यापन पर्याप्त है, और स्थानीय आलोचनात्मक घटना को समझने के लिए महत्वपूर्ण मूल्य रखता है।