Модель купи піску Бака-Тана-Візенфельда (BTW, 1987) є канонічним прикладом самоорганізованої критичності (СОК): система «накопичення-дисипація», що самостійно налаштовується до критичного стану без зовнішнього підбору параметрів. Зерна піску додаються по одному до випадкових клітинок; коли клітинка накопичує чотири або більше зерен, вона «перекидається» — розподіляє по одному зерну кожному з чотирьох сусідів (зерна на межі губляться). Розподіл розмірів лавин підпорядковується степеневому закону P(s) ∝ s−3/2 — без характерного масштабу.
Додавайте зерна вручну або безперервно і спостерігайте каскади перекидань. Лог-лог-графік розмірів лавин накопичується в реальному часі, виявляючи сигнатуру критичності — прямолінійний ступеневий закон.
Що таке самоорганізована критичність (СОК)?
СОК описує клас динамічних систем, що природно еволюціонують до критичного стану — межі між порядком і хаосом — без точного підбору зовнішніх параметрів. На критичності флуктуації всіх розмірів відбуваються зі степеневою частотою. СОК запропонована Пером Баком, Чао Таном і Куртом Візенфельдом у 1987 р. як єдине пояснення степеневих розподілів у природі: від магнітуд землетрусів (закон Гутенберга-Ріхтера) до сонячних спалахів і нейронних лавин у мозку.
Які правила перекидання в моделі BTW?
Кожна клітинка (i, j) сітки містить ціле число зерен z(i,j). Коли z(i,j) ≥ 4, клітинка перекидається: z(i,j) → z(i,j) − 4, і кожен з чотирьох сусідів отримує по 1 зерну. Граничні клітинки втрачають зерна, що виходять за межі (ця дисипація необхідна для досягнення стаціонарного стану). Перекидання повторюються — одне може спровокувати сусідів — доти, доки всі клітинки не матимуть менше 4 зерен. Весь каскад рахується як одна лавина.
Чому розподіл розмірів лавин підпорядковується степеневому закону?
У стані СОК система є інваріантною за масштабом: характерного розміру лавини не існує. Малі збурення можуть необмежено посилюватися, бо купа всюди водночас перебуває на межі стабільності. Ця масштабна інваріантність — ознака фазового переходу другого роду, але тут система досягає критичності автоматично, без підбору параметра керування. Показник τ ≈ 1,5 (P(s) ∝ s−τ) є універсальним для 2D моделі BTW.
Абелева купа піску на скінченній сітці з відкритими межами аналізується за допомогою теорії абелевих груп: правила перекидання комутують, що робить модель «абелевою купою піску». Алгоритм горіння Дхара (1990) ефективно визначає, чи є конфігурація рекурентною. У термодинамічній границі критичні показники та кореляційні функції отримані методами конформної теорії поля (відповідність логарифмічній КТП c=−2).
СОК запропонована для: землетрусів (закон Гутенберга-Ріхтера: log N ∝ −bM з b ≈ 1); сонячних спалахів (степеневий розподіл енергії протягом 8 порядків); нейронних лавин у корі (каскади P(s) ∝ s−3/2); лісових пожеж; волатильності фондового ринку; і подій вимирання в палеонтології. В кожному разі степеневий закон охоплює кілька десятків — що свідчить про спільний механізм, хоча для деяких систем СОК як пояснення залишається дискусійною.
Реальні купи піску мають характерний розмір лавин, визначений формою зерен, тертям та геометрією купи — степеневих розподілів у простих експериментах не спостерігається. Однак досліди Фретте з видовженими зернами рису (1996) дійсно виявили лавини за степеневим законом — найближча до BTW фізична реалізація СОК.
СОК була запропонована частково для пояснення рожевого шуму (1/f): спектральна потужність масштабується як 1/f^α (α ≈ 1) в широкому діапазоні частот. Рожевий шум зустрічається у флуктуаціях напруги в резисторах, інтервалах серцебиття, музиці та витраті річок. Модель BTW породжує 1/f-подібні часові флуктуації як наслідок безмасштабної динаміки лавин, пов'язуючи просторовий степеневий закон з часовим.
Абелева група куп піску має унікальний тотожний елемент: конфігурацію zid таку, що додавання zid до будь-якої рекурентної конфігурації та стабілізація повертає ту саму конфігурацію. У 2D тотожній конфігурації властиві складні фрактальні візерунки з точною тристоронньою симетрією, що виникають суто з алгебраїчної структури — без жодного «ручного» задавання.
Нейронні лавини в культурах кіркових зрізів підпорядковуються степеневим законам, узгодженим з СОК (P(s) ∝ s−3/2), — вперше показано Беггсом і Пленцем у 2003 р. Вважається, що кора самоорганізується до критичності, бо критичний стан максимізує передачу інформації, динамічний діапазон і чутливість до входів. Порушення степеневого закону пов'язане з епілептичними нападами (надкритичний стан) та анестезією (підкритичний).
При звичайних рівноважних фазових переходах (наприклад, феромагнетик Ізінга) критичність виникає при єдиній точно налаштованій температурі Tc. Відхилення T від Tc руйнує критичність. У системах із СОК повільне накопичення (додавання зерен) врівноважується дисипацією (втрати на межі), утворюючи автоматичну зворотну петлю: при підкритичному стані зерна накопичуються до досягнення критичності; при надкритичному — великі лавини швидко знижують критичність.
Так — багата сімейство: модель Манни (випадковий напрямок перекидання) має інші критичні показники; моделі лісових пожеж демонструють СОК через ріст дерев та запалення від блискавки; модель Оламі-Федера-Крістенсена (OFC) використовує безперервні дійсні напруження і є стандартною моделлю СОК для землетрусів. Кожна модель демонструє безмасштабну динаміку з різними показниками, геометрією та класами універсальності.
Ця симуляція відтворює модель купи піску Бака-Тана-Візенфельда — класичний приклад самоорганізованої критичності. Зерна накопичуються на сітці, і щойно клітинка досягає критичної висоти K, вона перекидається, віддаючи по одному зерну кожному з чотирьох сусідів — ця подія може розростися в лавину будь-якого розміру, від однієї клітинки до каскаду, що охоплює всю сітку. Розподіл розмірів лавин підпорядковується степеневому закону без жодного характерного масштабу — ознака системи, що сама вивела себе на межу стійкості.
Кожна клітинка забарвлена відповідно до кількості зерен у ній: від майже чорного (порожньо) через зелений і жовтогарячий до червоного біля критичної висоти. Каскади перекидань розходяться сплесками кольору, а графік під сіткою в реальному часі фіксує розмір кожної лавини — його майже пряма лінія і є видимою ознакою степеневого закону.
Оберіть режим скидання — Центр, Випадково або Клік, — щоб задати, куди падають зерна. Повзунок «Критична висота K» (2–8) визначає поріг перекидання, «Швидкість» (1–50 зерен за кадр) — темп додавання піску, а «Сітка» (від 30 на 30 до 120 на 120) змінює розмір поля. Кнопки Пауза/Старт і Скинути — нижче, а в режимі Клік достатньо натиснути на полотно, щоб висипати невелику купку зерен саме там.
Коли Пер Бак, Чао Тан і Курт Візенфельд представили цю модель 1987 року, вони запропонували самоорганізовану критичність як пояснення степеневих закономірностей, що трапляються в природі повсюдно — від магнітуди землетрусів до сонячних спалахів — і виникають без жодного зовнішнього налаштування.
K визначає, скільки зерен може втримати клітинка, перш ніж перекинутися. Збільшення K відтягує момент перекидання, дозволяючи зернам накопичуватися довше перед лавиною, а зменшення (аж до 2) змушує перекидання — і каскади — починатися набагато швидше.
Симуляція використовує відкриту межу: коли перекидання відправляє зерно за край поля, воно просто покидає систему, а не з'являється з протилежного боку. Саме ця втрата дозволяє купі встановитися у стійкий критичний стан, а не рости безмежно.
Центр завжди додає зерна в одну клітинку посередині сітки — класичний спосіб вирощувати купу з точкового джерела. Випадково розсіює зерна по всій сітці щокадру, а Клік дозволяє висипати невелику купку з двадцяти зерен саме там, де ви натиснете, щоб викликати лавину на власний вибір.
Графік показує, як часто трапляються лавини кожного розміру, причому обидві осі логарифмічні. Оскільки розміри лавин підпорядковані степеневому закону, залежність кількості від розміру в таких осях утворює майже пряму лінію — це доказ того, що лавини найрізноманітніших розмірів підпорядковуються тому самому простому правилу перекидання.
Більша сітка дає лавинам більше простору для розростання, тож найбільші каскади на полі 120 на 120 можуть у рази перевищувати ті, що трапляються на полі 30 на 30, хоча сам показник степеневого закону лишається тим самим. Зміна розміру сітки скидає симуляцію, адже конфігурація зерен має сенс лише за початкового розміру поля.