Система, що самоорганізовується для досягнення критичної точки
Опустіть зерна піску поодинці, один за одним, на плоску сітку. Більшість падає тихо і нічого не відбувається. Іноді додавання одного зерна запускає лавину, яка розповсюджується по всій купі. Ніхто не налаштовує систему так, щоб вона поводилася таким чином — вона сама себе налаштовує, рухаючись до делікатної критичної точки, де можуть виникати лавини будь-якого розміру. Це самоорганізована критикальність, представлена у 1987 році фізиками Пером Баком, Чяо Тангом та Курт Візенфельдом у статті, яка стала однією з найпопулярніших у фізиці. Однакова статистична відбиток зустрічається в природі: під час землетрусів, лісових пожеж, фінансових криз і спалахів нейронної активності, тому скромна купа піску стала єдиним дзеркалом для розуміння того, як часто складні системи живуть на межі нестабільності, а не заспокоюються у тихій рівновазі.
Зворотній процес і критична висота K
Кожна клітинка сітки містить ціле число зернин, її висоту. Зерна додаються по черзі до обраної клітинки. Коли висота клітинки досягає критичного порогу K, вона стає нестабільною і падає: втрачає K зерен та додає 1 до кожного з чотирьох сусідніх ортогональних клітинок. Один падіж може скинути сусідню клітинку на її власний поріг, яка, в свою чергу, падає також, і ланцюгова реакція — лавина — поширюється назовні до тих пір, поки всі клітинки не стануть стабільними знову. Зерна, що падають за межі сітки, просто зникають, розсіюючи енергію, що вводиться кожним новим зерном, і зупиняючи нагромаження зернин у нескінченність.
Якщо z(x,y) >= K: z(x,y) -= K z(x±1,y) += 1 (кожна з 4 ортогональних сусідніх клітинок) z(x,y±1) += 1 P(s) ~ s^(-tau) розподіл розмірів лавин, tau ≈ 1–1.5 у 2D Зміна K змінює те, скільки зернин може накопичити клітинка перед тим, як випустити енергію до своїх сусідів, що переформує статистику лавин без зміни якісного зображення: незалежно від встановленого значення K, підтримка випадкового падіння все ще змушує сітку рухатися до свого власного критичного пункту.
if z(x,y) >= K:
z(x,y) -= K
z(x±1,y) += 1 (each of the 4 orthogonal neighbours)
z(x,y±1) += 1
P(s) ~ s^(-tau) avalanche-size distribution, tau ≈ 1–1.5 in 2D
Закон влади: чому немає «типового» обвалу
Коли пісок усередині поступово стабілізується, розмір наступного обвалу — виміряний за кількістю обвалів або площею, що охоплюється — фундаментально непередбачуваний. Якщо на графіку показати частоту зустрічання кожної величини, розподіл буде слідувати закону влади: P(s) ~ s^(-tau). На логарифмічній діаграмі це пряма лінія, що охоплює багато порядків величини, і наслідок є глибоким: немає характерного масштабу. Малі події, пов'язані з одним обвалом, є поширеними, тоді як катастрофи, що охоплюють всю систему, рідкісні, але ніколи неможливі. Цей інверсійний розподіл є ознакою системи, яка знаходиться точно у критичній точці — такому стані, який зазвичай вимагає ретельного налаштування параметра, такого як температура, в рівноважній фізиці. Повільний процес додавання піску у пісочницю, збалансований швидким розслабленням і дисперсією на межах, автоматично підтримує систему в критичному стані, без необхідності обертання будь-якої ручки.
Зде́сь ті самі статистичні дані повторюються
Самоорганізована критicità запропоновано як пояснювальний механізм у надзвичайно широкому діапазоні дисциплін. Землетруси дотримуються закону Гутенберга-Ріхтера, що є законом степеня щодо величини та частоти, який відображає цикл «накопичення та зсуву» піскової гори як накопичення тектонічного напруження та його вивільнення у вигляді зсувів різного розміру. Лісові пожежі часто демонструють законом степеня щодо частоти та площі, що спалено, оскільки рослинність поступово накопичується, а пожежа швидко поширюється. Дослідження показують, що каскади активності в кортиковальних мережах — нейронні зсуви — можуть слідувати подібним розподілам за розмірами степеня, що спонукає гіпотезу про те, що мозок працює поблизу критичної точки для оптимізації обробки інформації. Великі рухи цін на фінансових ринках і утворення заторів також вивчалися через цей же механізм малих місцевих взаємодій, які іноді запускають системні каскади.
Що не стверджує модель
Часто виникає непорозуміння, що модель пісковину точно описує реальний піщаний потік – це не так. Справжні зернисті матеріали включають інерцію, тертя та вологість і схильні до великих квазіпериодичних зсувів замість чіткої статистики із законами степеней. Самоорганізована криттєвість також не є єдиним універсальним законом, який пояснює всі складні системи; це один механізм серед багатьох, і відповідність за законами степенів само по собі не доводить криттєвості, оскільки інші процеси можуть імітувати такі розподіли. Нарешті, криттєвість не слід плутати з хаосом: самоорганізована критична система – це статистично стабільний стан із масштабованими флуктуаціями, а не детермінований хаос у сенсі логістичної карти або аттрактора Лоренца.
Frequently asked questions
Що таке самоорганізована критичність?
Самоорганізована критичність – це схильність певних повільно-приводів енергії, розсіяних систем до еволюції, без зовнішнього налаштування, до критичного стану між порядком і хаосом, де події будь-якого розміру відбуваються відповідно до закону потужності.
Чому зсуви у піщаній купі слідують за законом потужності?
Коли купа досягає свого критичного стабільного стану, більше немає типового розміру зсуву. Ймовірність зсуву розміром s масштабується як P(s) ~ s^(-tau), тому невеликі падіння є поширеними, а системно-масштабні каскади рідкісні, але ніколи неможливі – відбиток системи, яка знаходиться точно на критичній точці.
Чи описує модель піщаної купи справжній пісок?
Не зовсім. Реальні гранульові купки включають інерцію, тертя та вологість і схильні виробляти великі квазіперіодичні зсуви замість чистого закону потужності. Модель є абстрактною метафорою критичності, а не фізичною симуляцією гранул.
Спробуйте наживо
Усе, що вище, працює прямо у вашому браузері — відкрийте the simulation і змінюйте параметри під час роботи. Нічого не встановлюється, нічого не завантажується на сервер, уся модель живе в одній вкладці.
▶ Відкрити симуляцію the simulation