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.
- معرّف الورقة: 2510.09631
- العنوان: Scaling Properties of Avalanche Activity in the Two-Dimensional Abelian Sandpile Model
- المؤلف: Anubhav Ganguly (IISER Mohali)
- التصنيف: 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، وفهم أنماط نشاطه المحلي أمر حاسم لفهم الظواهر الحرجة
- القيمة التطبيقية العملية: في الأنظمة الحقيقية، غالباً ما يكون من الممكن ملاحظة النشاط المحلي فقط، وليس إحصائيات الانهيار العالمية
- التكامل النظري: يوفر معلومات تكميلية لإحصائيات حجم الانهيار التقليدية، ويكشف المناطق التي تعمل كقنوات متكررة لانتشار الانهيار
- ركزت الدراسات التقليدية بشكل أساسي على توزيع قانون القوة Prob(s)~s⁻τ لحجم الانهيار s(R)
- تفتقر إلى التحليل المنهجي لتكرار مشاركة الموقع في الانهيارات
- لم يتم دراسة توزيع النشاط 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'
الهدف من البحث هو تحليل التوزيع الاحتمالي لدالة النشاط 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
- وصول النظام إلى الحالة الحرجة
- حساب A(R) لجميع المواقع لكل تكوين C
- إحصاء التوزيع الاحتمالي التجريبي 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*: يهيمن التناقص الأسي
بالنسبة للتكوينات النادرة C_ε ذات الكتلة الزائدة εL²:
- ⟨A⟩_C ~ εL³ (الحد الأساسي)
- كسر هذه التكوينات ~exp(-ε²L²)
- يؤدي إلى F(x~εL) ≳ exp(-x²) للـ x الكبير
تُظهر الأشكال 6 و9 انهياراً جيداً للبيانات من أحجام نظام متعددة تحت شكل التحجيم، مما يتحقق من صحة التنبؤات النظرية.
- نموذج Bak-Tang-Wiesenfeld1,2: العمل الرائد في SOC
- تعيين Majumdar-Dhar9: التكافؤ بين ASM ونموذج Potts عند q→0
- وصف نظرية المجال المطابق: الشحنة المركزية c=-2، الأس الديناميكي z=5/4
- الدراسات التقليدية: تركزت على قانون القوة Prob(s)~s⁻τ لتوزيع حجم الانهيار s(R)
- النماذج الحتمية مقابل العشوائية: نموذج 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 و Wiesenfeld1,2
- المساهمات النظرية لـ Dhar4,5
- النتائج الدقيقة لـ Majumdar-Dhar9
- التطورات التجريبية الحديثة22,23
التقييم الشامل: هذا عمل مبتكر في دراسة نموذج كومة الرمل الأبيلي، حيث يوفر منظوراً جديداً لنظرية SOC من خلال إدخال مفهوم توزيع نشاط الموقع. على الرغم من وجود بعض القيود في حجم النظام والنظرية التحليلية، فإن إطاره النظري كامل، والتحقق العددي كافٍ، وله قيمة مهمة لفهم الظواهر الحرجة المحلية.