Модель абелевой кучи - Abelian sandpile model

Идентификационный элемент песчаной группы прямоугольной сетки. Желтые пиксели соответствуют вершинам, несущим три частицы, сиреневые - двум частицам, зеленые - одной и черные - нулю.

В Модель абелевой кучи, также известный как Модель Бак – Танга – Визенфельда, был первым обнаруженным примером динамическая система отображение самоорганизованная критичность. Он был представлен Per Bak, Чао Тан и Курт Визенфельд в статье 1987 года.[1]

Модель представляет собой клеточный автомат. В своей первоначальной формулировке каждому участку на конечной сетке соответствует значение, соответствующее уклону сваи. Этот уклон нарастает по мере того, как «песчинки» (или «щепки») случайным образом кладут на сваю, пока уклон не превысит определенное пороговое значение, когда это место обрушится, перенося песок на соседние участки, увеличивая их уклон. Бак, Тан и Визенфельд рассмотрели процесс последовательного случайного размещения песчинок на сетке; каждое такое размещение песка на определенном участке может не иметь никакого эффекта или вызвать каскадную реакцию, которая затронет многие участки.

С тех пор модель изучалась на бесконечной решетке, на других (неквадратных) решетках и на произвольных графах (включая ориентированные мультиграфы).[2] Это тесно связано с долларовая игра, вариант игра-стрелялка представленный Биггсом.[3]

Определение (прямоугольные сетки)

Модель песчаной кучи клеточный автомат первоначально определено на прямоугольная сетка (шахматная доска) из стандартная квадратная решетка .К каждой вершине (сторона, поле) сетки ставим в соответствие значение (песчинки, склон, частицы) , с называется (исходной) конфигурацией песчаной кучи.

Динамика автомата на итерации затем определяются следующим образом:

  1. Выберите случайную вершину по некоторому распределению вероятностей (обычно равномерному).
  2. Добавьте одну песчинку в эту вершину, оставляя номера зерен для всех остальных вершин неизменными, т.е. установите
    и
    для всех .
  3. Если все вершины стабильный, т.е. для всех , а также конфигурация считается стабильным. В этом случае перейдите к следующей итерации.
  4. Если хотя бы одна вершина неустойчивый, т.е. для некоторых , вся конфигурация считается нестабильным. В этом случае выберите любую неустойчивую вершину наугад. Опрокинуть эту вершину, уменьшив количество зерен на четыре и увеличив количество зерен каждого из (максимум четырех) прямых соседей на единицу, т.е.
    , и
    если .
    Если вершина на границе домена опрокидывается, это приводит к чистой потере зерен (два зерна в углу сетки, одно зерно в противном случае).
  5. Из-за перераспределения зерен опрокидывание одной вершины может сделать другие вершины нестабильными. Таким образом, повторяем процедуру опрокидывания до тех пор, пока все вершины в конечном итоге станет стабильным и продолжит выполнение следующей итерации.

Падение нескольких вершин за одну итерацию называется лавина. Гарантируется, что каждая лавина в конечном итоге остановится, т.е. после конечного числа опрокидываний достигается некоторая устойчивая конфигурация, при которой автомат определен правильно. Более того, хотя часто будет много возможных вариантов порядка опрокидывания вершин, окончательная стабильная конфигурация не зависит от выбранного порядка; в этом смысле песочная куча абелевский. Точно так же количество раз, когда каждая вершина опрокидывается во время каждой итерации, также не зависит от выбора порядка опрокидывания.

Определение (неориентированные конечные мультиграфы)

Обобщить модель песочницы с прямоугольной сетки стандартной квадратной решетки на произвольный неориентированный конечный мультиграф , особая вершина называется раковина указано, что нельзя опрокидывать. А конфигурация (состояние) модели тогда является функцией подсчет неотрицательного числа зерен на каждой вершине, не являющейся стоком. Вершина, не являющаяся стоком с

нестабильно; он может быть опрокинут, что отправит одну из своих крупинок каждому из своих (не принимающих) соседей:

для всех , .

Затем клеточный автомат работает, как и раньше, то есть путем добавления на каждой итерации одной частицы к случайно выбранной вершине, не являющейся стоком, и опрокидыванию, пока все вершины не станут стабильными.

Определение модели песчаной кучи, данное выше для конечных прямоугольных сеток стандартной квадратной решетки в таком случае можно рассматривать как частный случай этого определения: рассмотрим граф который получается из добавляя дополнительную вершину, сток, и рисуя дополнительные ребра от стока к каждой граничной вершине так что степень каждой не принимающей вершины четыре. Таким образом могут быть определены также модели песчаных куч на непрямоугольных сетках стандартной квадратной решетки (или любой другой решетки): пересечение некоторого ограниченного подмножества из с . Сжимайте каждую грань из чьи две конечные точки не находятся в . Единственная оставшаяся вершина вне затем составляет сток результирующего графа песчаной кучи.

Переходные и повторяющиеся конфигурации

В определенной выше динамике песочного автомата некоторые устойчивые конфигурации ( для всех ) появляются бесконечно часто, в то время как другие могут появляться только конечное число раз (если вообще появляются). Первые упоминаются как повторяющиеся конфигурации, а последние называются переходные конфигурации. Таким образом, повторяющиеся конфигурации состоят из всех стабильных неотрицательных конфигураций, которые могут быть достигнуты из любой другой стабильной конфигурации путем многократного добавления песчинок в вершины и опрокидывания. Нетрудно заметить, что минимально стабильная конфигурация , где каждая вершина несет песчинки, достижимо из любой другой стабильной конфигурации (добавить зерна в каждую вершину). Таким образом, эквивалентным образом повторяющиеся конфигурации - это именно те конфигурации, которые могут быть достигнуты из минимально стабильной конфигурации, только добавляя песчинки и стабилизируя.

Не всякая неотрицательная стабильная конфигурация повторяется. Например, в каждой модели кучи песка на графе, состоящем как минимум из двух связанных вершин, не являющихся раковиной, каждая устойчивая конфигурация, в которой обе вершины несут ноль песчинок, не повторяется. Чтобы доказать это, сначала обратите внимание, что добавление песчинок может только увеличить общее количество песчинок, переносимых двумя вершинами вместе. Таким образом, чтобы достичь конфигурации, в которой обе вершины несут нулевые частицы из конфигурации, где это не так, необходимо выполнить шаги, на которых по крайней мере одна из двух вершин опрокидывается. Рассмотрим последний из этих шагов. На этом этапе одна из двух вершин должна опрокинуться последней. Поскольку при опрокидывании песчинка переносится в каждую соседнюю вершину, это означает, что общее количество песчинок, переносимых обеими вершинами вместе, не может быть меньше единицы, что завершает доказательство.

Песчаная группа

Учитывая конфигурацию , для всех , опрокидывание нестабильных вершин, не являющихся стоком, на конечном связном графе до тех пор, пока не останется неустойчивых вершин, не являющихся стоком, приводит к единственному стабильный конфигурация , который называется стабилизация из . Учитывая две стабильные конфигурации и , мы можем определить операцию , что соответствует добавлению зерен по вершинам с последующей стабилизацией полученной кучи песка.

При произвольном, но фиксированном порядке вершин, не являющихся приемником, несколько операций опрокидывания, которые могут, например, возникают во время стабилизации нестабильной конфигурации, могут быть эффективно закодированы с помощью граф лапласиан , куда это матрица степеней и это матрица смежности графа. Удаление строки и столбца соответствующий сток дает редуцированный граф лапласиан . Затем, начиная с конфигурации и опрокидывая каждую вершину Всего раз дает конфигурацию , куда - продукт сокращения. Кроме того, если соответствует количеству раз, когда каждая вершина опрокидывается во время стабилизации данной конфигурации , тогда

В этом случае, называется падение или же функция одометра (стабилизации ).

Под операцией , набор повторяющихся конфигураций образует абелева группа изоморфна коядру лапласиана редуцированного графа , т.е. к , Посредством чего обозначает количество вершин (включая раковину). В более общем смысле, набор стабильных конфигураций (переходных и повторяющихся) образует коммутативный моноид под операцией . Минимальный идеальный этого моноида тогда изоморфна группе рекуррентных конфигураций.

Группа, образованная повторяющимися конфигурациями, а также группа которому первый изоморфен, чаще всего называют песчаная группа. Другие общие названия той же группы: критическая группа, Группа якобиана или (реже) Группа Пикард. Обратите внимание, однако, что некоторые авторы обозначают только группу, образованную рекуррентными конфигурациями, как группу песочницы, оставляя за собой название якобиева группа или критическая группа для (изоморфной) группы, определенной формулой (или для связанных изоморфных определений). Наконец, некоторые авторы используют название группа Пикара для обозначения прямого продукта группы песчаных куч и , который естественным образом проявляется в клеточном автомате, тесно связанном с моделью кучи песка, называемой «запуском чипов» или «долларовой игрой».

Учитывая изоморфизмы, указанные выше, порядок группы песочницы является определителем , что по теорема о матричном дереве - количество остовных деревьев графа.

Самоорганизованная критичность

Первоначальный интерес к модели вызван тем, что при моделировании на решетках ее привлекают критическое состояние, в этот момент длина корреляции системы и время корреляции системы стремятся к бесконечности без какой-либо точной настройки параметра системы. Это контрастирует с более ранними примерами критических явлений, такими как фазовые переходы между твердым телом и жидкостью или жидкостью и газом, где критическая точка может быть достигнута только путем точной настройки (например, температуры). Следовательно, в модели песчаной кучи можно сказать, что критичность самоорганизованный.

Когда модель песчаной кучи достигает своего критического состояния, корреляция между реакцией системы на возмущение и детали возмущения. Обычно это означает, что падение еще одной песчинки на сваю может ничего не вызвать или может привести к обрушению всей сваи в результате массивного оползня. Модель также отображает 1/ƒ шум, характерная для многих сложных систем в природе.

Эта модель отображает критическое поведение только в двух или более измерениях. Модель песчаной кучи может быть выражена в 1D; однако, вместо того, чтобы развиваться до своего критического состояния, одномерная модель песчаной кучи достигает минимально устойчивого состояния, когда каждый узел решетки приближается к критическому уклону.

Для двух измерений была выдвинута гипотеза, что связанные конформная теория поля состоит из симплектические фермионы с центральный заряд c = −2.[4]

Характеристики

Принцип наименьшего действия

Стабилизация конфигураций микросхем подчиняется форме принцип наименьшего действия: каждая вершина опрокидывается не больше, чем необходимо в процессе стабилизации.[5] Это можно формализовать следующим образом. Назовите последовательность падений законный если он только опрокидывает неустойчивые вершины, и стабилизирующий если это приведет к стабильной конфигурации. Стандартный способ стабилизации песчаной кучи - найти максимальную законную последовательность; то есть, опрокидывая до тех пор, пока это возможно. Такая последовательность, очевидно, является стабилизирующей, и абелевым свойством песочной кучи является то, что все такие последовательности эквивалентны с точностью до перестановки порядка опрокидывания; то есть для любой вершины , количество раз topples одинаково во всех законных стабилизирующих последовательностях. Согласно принципу наименьшего действия минимальная стабилизация Последовательность также эквивалентна перестановке порядка отмены в допустимую (и все еще стабилизирующую) последовательность. В частности, конфигурация, полученная в результате минимальной стабилизирующей последовательности, такая же, как и в результате максимальной правовой последовательности.

Более формально, если вектор такой, что - количество раз, когда вершина опрокидывается при стабилизации (за счет опрокидывания нестабильных вершин) конфигурации микросхемы , и является целым вектором (не обязательно неотрицательным) такой, что стабильно, то для всех вершин .

Пределы масштабирования

Анимация идентичности песчаной кучи на квадратных сетках увеличивающегося размера. Черным цветом обозначены вершины с 0 зернами, зеленым - 1, фиолетовым - 2, а золотым - 3.

Анимация показывает повторяющуюся конфигурацию, соответствующую идентичности группы кучи на разных квадратные сетки увеличивающихся размеров , в результате чего конфигурации масштабируются, чтобы всегда иметь один и тот же физический размер. Визуально идентичности на более крупных сетках, кажется, становятся все более детализированными и «сходятся в непрерывное изображение». Математически это предполагает существование масштабных пределов идентичности песчаной кучи на квадратных сетках, основанных на понятии слабой * сходимости (или некотором другом обобщенном понятии сходимости). Действительно, существование пределов масштабирования повторяющихся конфигураций песчаных куч было доказано Уэсли Пегденом и Чарльзом Смартом.[6].[7] В дальнейшей совместной работе с Лайонелом Левином они использовали предел масштабирования для объяснения фрактальной структуры кучи песка на квадратных сетках.[8]

Обобщения и родственные модели

Модели песчаных куч на бесконечных сетках

30 миллионов зерен упали на участок бесконечной квадратной сетки, а затем опрокинулись в соответствии с правилами модели кучи песка. Белым цветом обозначены участки с 0 зернами, зеленым - 1, фиолетовым - 2, золотым - 3. Ограничивающая рамка - 3967 × 3967.

Существует несколько обобщений модели кучи песка на бесконечные сетки. Проблема в таких обобщениях состоит в том, что, как правило, больше не гарантируется, что каждая лавина в конечном итоге остановится. Таким образом, некоторые обобщения рассматривают только стабилизацию конфигураций, для которых это может быть гарантировано.

Достаточно популярная модель на (бесконечной) квадратной решетке с узлами определяется следующим образом:

Начнем с некоторой неотрицательной конфигурации значений что конечно, то есть

Любой сайт с

является неустойчивый и может опрокидывать (или же Огонь), отправив по одной из своих фишек каждому из четырех своих соседей:

Поскольку исходная конфигурация конечна, процесс гарантированно прекратится, и зерна разлетятся наружу.

Популярный частный случай этой модели дается, когда начальная конфигурация равна нулю для всех вершин, кроме начала координат. Если исходная точка несет огромное количество песчинок, конфигурация после релаксации образует фрактальные узоры (см. Рисунок). Когда исходное количество зерен в начале координат стремится к бесконечности, было показано, что масштабированные стабилизированные конфигурации сходятся к уникальному пределу.[7][8]

Модели песочницы на ориентированных графах

Модель песчаной кучи может быть обобщена на мультиграфы произвольной направленности. Правила таковы, что любая вершина с

нестабильно; опрокидывание снова отправляет фишки каждому из своих соседей, по одной на каждом исходящем ребре:

и для каждого :

куда количество ребер из к .

В этом случае матрица лапласиана не симметрична. Если указать раковину такой, что есть путь от любой другой вершины до , то операция стабилизации на конечных графах корректно определена и группу песочницы можно записать

как прежде.

Порядок группы кучи снова является определяющим фактором , что по общей версии теорема о матричном дереве количество ориентированных остовные деревья укоренился в раковине.

Расширенная модель песчаной кучи

Динамика песчаной насыпи, индуцированная гармонической функцией H = x * y на квадратной сетке 255x255.

Чтобы лучше понять структуру группы песчаных куч для различных конечных выпуклых сеток стандартной квадратной решетки , Ланг и Школьников представили расширенная модель песчаной кучи в 2019 году.[9] Расширенная модель песчаной кучи определяется почти так же, как и модель обычная модель песчаника (т.е. исходная модель Бака – Танга – Визенфельда [1]), за исключением того, что вершины на границе сетки теперь разрешено нести неотрицательное действительное количество зерен. В отличие от этого, вершины внутри сетки по-прежнему могут нести только целое число зерен. Правила опрокидывания остаются неизменными, т.е. предполагается, что как внутренние, так и граничные вершины становятся нестабильными и опрокидываются, если число зерен достигает или превышает четыре.

Также повторяющиеся конфигурации расширенной модели песчаной кучи образуют абелеву группу, называемую расширенная песчаная группа, из которых обычная группа песчаных дискретная подгруппа. В отличие от обычной песчаной группы, расширенная песчаная группа, однако, является непрерывной. Группа Ли. Поскольку он создается только добавлением песчинок к границе сетки расширенная группа песчаных куч, кроме того, имеет топология из тор измерения и объем дан по порядку обычной песчаной группы.[9]

Особый интерес представляет вопрос, как рекуррентные конфигурации динамически изменяются вдоль непрерывного геодезические этого тора, проходящего через единицу. Этот вопрос приводит к определению динамики песчаной насыпи.

(расширенная модель песчаной кучи)

соответственно

(обычная песочная модель)

индуцированный целочисленным гармоническая функция вовремя , с идентичность песчаной группы и функция пола.[9] Для полиномиальных гармонических функций низкого порядка динамика песчаной кучи характеризуется плавным преобразованием и очевидным сохранением пятен, составляющих идентичность песчаной кучи. Например, гармоническая динамика, индуцированная напоминают «плавное растяжение» идентичности по основным диагоналям, визуализированное в анимации. Конфигурации, возникающие в динамике, индуцированной одной и той же гармонической функцией на квадратных сетках разного размера, кроме того, предположительно имеют слабую * сходимость, что означает, что для них предположительно существуют пределы масштабирования.[9] Это предлагает естественный перенормировка для расширенных и обычных групп песочницы, что означает отображение повторяющихся конфигураций на заданной сетке на повторяющиеся конфигурации на подсетке. Неформально, эта перенормировка просто отображает конфигурации, появляющиеся в данный момент времени. в динамике песчаной насыпи, индуцированной некоторой гармонической функцией на большей сетке к соответствующим конфигурациям, которые одновременно проявляются в динамике песчаной насыпи, вызванной ограничением в соответствующую подсетку.[9]

Делимая куча песка

Сильно связанная модель - это так называемая делимая модель песчаной кучи, представленный Левином и Пересом в 2008 году,[10] в котором вместо дискретного количества частиц в каждом узле , есть реальное число представляющий количество массы на сайте. Если такая масса отрицательна, то это можно понимать как дыру. Падение происходит всякий раз, когда сайт имеет массу больше 1; он распределяет излишки равномерно между своими соседями, что приводит к тому, что, если сайт заполнен , он будет полон на все более поздние времена.

Культурные ссылки

Песчаная куча Бак – Тан – Визенфельд упоминается на Numb3rs Эпизод «Буйство», как математик Чарли Эппес объясняет своим коллегам решение уголовного расследования.

В компьютерная игра Hexplode основан на абелевой модели песчаной кучи на конечной гексагональной сетке, где вместо случайного размещения зерен игроки размещают зерна.

Рекомендации

  1. ^ а б Бак, П.; Тан, К.; Визенфельд, К. (1987). «Самоорганизованная критичность: объяснение 1 /ƒ шум". Письма с физическими проверками. 59 (4): 381–384. Bibcode:1987ПхРвЛ..59..381Б. Дои:10.1103 / PhysRevLett.59.381. PMID  10035754.
  2. ^ Holroyd, A .; Levine, L .; Mészáros, K .; Перес, Ю .; Propp, J .; Уилсон, Б. (2008). Chip-Firing и Rotor-Routing на ориентированных графах. В равновесии и вне его 2. 60. С. 331–364. arXiv:0801.3306. Bibcode:1987ПхРвЛ..59..381Б. Дои:10.1007/978-3-7643-8786-0_17. ISBN  978-3-7643-8785-3. S2CID  7313023.
  3. ^ Биггс, Норман Л. (25 июня 1997 г.). "Запуск микросхем и критическая группа графа" (PDF). Журнал алгебраической комбинаторики: 25–45. Получено 10 мая 2014.
  4. ^ С. Могими-Араги; М. А. Раджабпур; С. Рухани (2004). "Абелева модель песчаной кучи: точка зрения теории конформного поля". Ядерная физика B. 718 (3): 362–370. arXiv:cond-mat / 0410434. Bibcode:2005НуФБ.718..362М. Дои:10.1016 / j.nuclphysb.2005.04.002. S2CID  16233977.
  5. ^ Fey, A .; Levine, L .; Перес Ю. (2010). «Темпы роста и взрывы в отвалах». Журнал статистической физики. 138 (1–3): 143–159. arXiv:0901.3805. Bibcode:2010JSP ... 138..143F. Дои:10.1007 / s10955-009-9899-6. ISSN  0022-4715. S2CID  7180488.
  6. ^ Пегден, Уэсли; Смарт, Чарльз (2017). «Устойчивость узоров в абелевой куче». arXiv:1708.09432 [math.AP ].
  7. ^ а б Пегден, Уэсли; Смарт, Чарльз (2013). «Схождение абелевой кучи». Математический журнал герцога. 162 (4): 627–642. arXiv:1105.0111. Дои:10.1215/00127094-2079677. S2CID  13027232.
  8. ^ а б Левин, Лайонел; Пегден, Уэсли (2016). «Аполлоническая структура в абелевой куче песка». Геометрический и функциональный анализ. 26 (1): 306–336. Дои:10.1007 / s00039-016-0358-7. HDL:1721.1/106972. S2CID  119626417.
  9. ^ а б c d е Ланг, Мориц; Школьников, Михаил (19.02.2019). «Гармоническая динамика абелевой кучи». Труды Национальной академии наук. 116 (8): 2821–2830. Дои:10.1073 / pnas.1812015116. ISSN  0027-8424. ЧВК  6386721. PMID  30728300.
  10. ^ Левин, Лайонел; Перес, Юваль (2008-10-29). "Сильная сферическая асимптотика для агрегации ротор-маршрутизатор и делимая куча песка". Возможный анализ. 30 (1): 1–27. arXiv:0704.0688. Дои:10.1007 / s11118-008-9104-6. ISSN  0926-2601. S2CID  2227479.

дальнейшее чтение

внешняя ссылка