沙堆模型

来自集智百科 - 复杂系统|人工智能|复杂科学|复杂网络|自组织
陈关荣（讨论 | 贡献）2020年11月16日 (一) 15:22的版本
跳到导航 跳到搜索


沙堆在矩形网格上的标识。黄色像素对应三颗沙粒的顶点，淡紫色代表两颗沙粒，绿色表示一颗沙粒，黑色表示零颗沙粒。


阿贝尔沙堆模型 Abelian sandpile model，也被称为 Bak-Tang-Wiesenfeld 模型，是第一个发现的动力系统展现自组织临界性的例子。它是由 Per Bak, Chao Tang 和 Kurt Wiesenfeld 在1987年的一篇论文[1]中提出的。


这个模型是一种 元胞自动机模型 Cellular automaton。在最初的公式中，有限网格上的每个位置都有一个与沙堆的坡度相对应的关联值。当“沙粒”（或“碎片”）被随机放置在沙堆上时，放置位置的斜坡就会堆积起来，直到倾斜程度超过一个特定的阈值，这个位置倒塌，沙子会转移到邻近的位置，增加它们的斜坡。Bak, Tang和 Wiesenfeld考虑了在网格上连续随机放置沙粒的过程; 每次这样在特定位置放置沙粒有可能不会产生影响，也有可能会引起级联反应，影响到周围的其他位置。


该模型已经在无限栅格、其他（非方形）栅格和任意图（包括有向多重图）上进行了研究。[2]它与美元博弈游戏密切相关，美元游戏是Biggs引入的一种chip-firing游戏的变体。[3]


定义（矩形网格）

沙堆模型是一个元胞自动机cellular automaton ，它最初定义在[math]\displaystyle{ N\times M }[/math]矩形网格（棋盘格）上，其顶点在标准的正方形格子[math]\displaystyle{ \mathbb{Z}^2 }[/math]上[math]\displaystyle{ \Gamma\subset\mathbb{Z}^2 }[/math]。


对于栅格中的每个顶点[math]\displaystyle{ (x,y)\in\Gamma }[/math]，我们关联一个值（沙粒、坡度、颗粒）[math]\displaystyle{ z_0(x,y)\in\{0,1,2,3\} }[/math]，这样所有顶点的初始状态[math]\displaystyle{ z_0\in\{0,1,2,3\}^\Gamma }[/math]被称为沙堆的（初始）构型。


自动机的动力学过程在第[math]\displaystyle{ i\in\mathbb{N} }[/math]次迭代时的定义如下:

  1. 根据某种概率分布（通常为均匀分布）随机选择一个顶点。根据一些概率分布（通常是均匀的）选择一个随机顶点[math]\displaystyle{ (x_i,y_i)\in\Gamma }[/math]。
  2. 向这个顶点添加一粒沙子，同时让其他顶点的沙粒数保持不变，也就是对于所有的[math]\displaystyle{ (x,y)\neq(x_i,y_i) }[/math]，设定
    [math]\displaystyle{ z_i(x_i,y_i)=z_{i-1}(x_i,y_i)+1 }[/math] 和
    [math]\displaystyle{ z_i(x,y)=z_{i-1}(x,y) }[/math]。
  3. 如果所有的顶点都是稳定的，即如果对于[math]\displaystyle{ (x,y)\in\Gamma }[/math]，[math]\displaystyle{ z_i(x,y)\lt 4 }[/math]，那么构型[math]\displaystyle{ z_i }[/math]被认为是稳定的。在这种情况下，继续下一轮迭代。
  4. 如果至少有一个顶点是不稳定的，即对于一些[math]\displaystyle{ (x_u,y_u)\in\Gamma }[/math]，[math]\displaystyle{ z_i(x_u,y_u)\geq 4 }[/math]，[math]\displaystyle{ z_i }[/math]被认为是不稳定的。在这种情况下，随机选择任意不稳定顶点[math]\displaystyle{ (x_u,y_u)\in\Gamma }[/math]。将该顶点的沙粒数减少4个，清空这个顶点，并将其每个（最多4个）直接邻居的沙粒数增加1个。即：
    [math]\displaystyle{ z_i(x_u,y_u) \rightarrow z_i(x_u,y_u) - 4, }[/math]， 如果 [math]\displaystyle{ ( x_u \pm 1, y_u\pm 1)\in\Gamma }[/math].
    ，
    [math]\displaystyle{ z_i( x_u \pm 1, y_u \pm 1) \rightarrow z_i( x_u \pm 1, y_u\pm 1) + 1 }[/math]。如果一个在边界的顶点产生崩塌，这将导致沙粒的净损失（两粒在网格的角落，一粒在其他地方）。 (two grains at the corner of the grid, one grain otherwise)
  5. 由于沙粒的重新分布，一个顶点的崩塌会使其他顶点不稳定。这样，重复崩塌的过程，直到[math]\displaystyle{ z_i }[/math]状态下的所有顶点最终稳定下来，继续下一轮迭代。
==Zcy（讨论） (two grains at the corner of the grid, one grain otherwise)翻译存疑 ==Zcy（讨论）


在一次迭代中多个顶点的崩塌被称为”雪崩“。每一次雪崩最终都会停止，也就是说，经过有限数量的顶点崩塌，会达到某种稳定的构型，这样自动机就得到了很好的定义。此外，尽管顶点崩塌的顺序常常有许多可能的选择，但最终的稳定状态并不依赖于所选择的顺序; 这是沙堆模型具有的可交换性质。类似地，在每次迭代过程中，每个顶点的崩塌次数也与崩塌顺序的选择是无关的。

定义（无向有限多重图）

为了将沙堆模型从标准方格的矩形网格推广到任意无向有限多重图 [math]\displaystyle{ G=(V,E) }[/math] ，在 [math]\displaystyle{ V }[/math] 中指定了一个不允许崩塌的特殊”沉没 sink“ 顶点[math]\displaystyle{ s }[/math]。模型的构型（状态）服从函数[math]\displaystyle{ z:V\setminus\{s\}\rightarrow\mathbb{N}_0 }[/math]，计算每个非沉没顶点上的非负沙粒数。非沉没顶点[math]\displaystyle{ v\in V\setminus\{s\} }[/math]当满足[math]\displaystyle{ z(v)\geq \deg(v) }[/math]时是不稳定的，它会产生崩塌，向给它的每个（非沉没）邻居分发一颗沙粒：

[math]\displaystyle{ z(v) \to z(v) - \deg(v) }[/math]
[math]\displaystyle{ z(u) \to z(u) + 1 }[/math]对于所有的[math]\displaystyle{ u\sim v }[/math], [math]\displaystyle{ u\neq s }[/math]


元胞自动机像之前一样进行，即在每次迭代中，向随机选择的非沉没顶点添加一个沙粒，不断进行崩塌过程，直到所有顶点都稳定。


上面给出的沙堆模型的定义，是在标准正方形网格[math]\displaystyle{ \mathbb{Z}^2 }[/math]上的有限矩形网格[math]\displaystyle{ \Gamma\subset\mathbb{Z}^2 }[/math]上，它可以看作是下面定义的一个特例：考虑图[math]\displaystyle{ G=(V,E) }[/math]，从[math]\displaystyle{ \Gamma }[/math]添加一个沉没顶点，并添加从沉没顶点到每个边界顶点的边，使得[math]\displaystyle{ G }[/math]的每个非沉没顶点的度数为4。以这种方式，也可以定义标准正方形网格（或任何其他类型网格）的非矩形格上的沙堆模型: 将[math]\displaystyle{ \mathbb{R}^2 }[/math]的一些有界子集[math]\displaystyle{ S }[/math]与[math]\displaystyle{ \mathbb{R}^2 }[/math]相交。收缩[math]\displaystyle{ \mathbb{Z}^2 }[/math]的每条边，其两个端点不在[math]\displaystyle{ S\cap\mathbb{Z}^2 }[/math]中。[math]\displaystyle{ S\cap\mathbb{Z}^2 }[/math]之外的一个单独剩余顶点构成了最终沙堆图的沉没顶点。

瞬态和循环构型

在上面定义的沙堆自动机的动力学过程中，一些稳定状态的构型（对于所有[math]\displaystyle{ v\in G\setminus\{s\} }[/math]，[math]\displaystyle{ 0\leq z(v)\lt 4 }[/math]）经常无限次出现，而另一些则只能出现有限次（如果真的发生的话）。前者被称为“循环构型”，而后者被称为“瞬态构型”。因此，周期性构型由所有稳定的非负构型组成，这些构型可以从任何其他稳定构型中，通过反复向顶点添加沙粒，产生崩塌而得到。很容易看出，“最小稳定构型”[math]\displaystyle{ zum }[/math]，其中每个顶点放置[math]\displaystyle{ z_m(v)=deg(v)-1 }[/math]颗沙粒，可从任何其他稳定构型得到（通过向每个顶点添加[math]\displaystyle{ deg(v)-z(v)-1\geq 0 }[/math]颗沙粒）。因此，也就是说，周期性构型可以从最小稳定构型开始，通过添加沙粒，再稳定化得到。


并非所有非负稳定构型都是循环的。例如，在一个至少包含由两个连通的非沉没顶点的图结构的沙堆模型中，如果这两个顶点没有放置沙粒，那么这个稳定结构是非循环的。为了证明这一点，首先要注意的是，沙粒的增加只能增加两个顶点放置的沙粒的总数。为了达到两个顶点都不放置沙粒的构型，从一个不是这种情况的构型出发，必然涉及到两个顶点中至少有一个崩塌的步骤。考虑这些步骤中的最后一个，在这个步骤中，两个顶点中的一个必须最后崩塌。由于崩塌会将对每个相邻的顶点转移一颗沙粒，这意味着两个顶点共同放置的沙粒总数不能低于一颗，因而得证。

沙堆群

给定一个构型[math]\displaystyle{ z }[/math]，[math]\displaystyle{ z(v)\in\mathbb{N}_0 }[/math]对于所有[math]\displaystyle{ v\in G\setminus\{s\} }[/math]，在有限连通图上使不稳定的非沉没顶点崩塌，直到没有不稳定的非汇顶点保留，这将导致唯一的“稳定”构型[math]\displaystyle{ z^\circ }[/math]，这就是[math]\displaystyle{ z }[/math]的“稳定化”。给定两个稳定构型[math]\displaystyle{ z }[/math]和[math]\displaystyle{ w }[/math]，我们可以定义运算[math]\displaystyle{ z*w \to (z+w)^\circ }[/math]，对应于沙粒的顶点方向相加，然后稳定得到的沙堆。corresponding to the vertex-wise addition of grains followed by the stabilization of the resulting sandpile. 

==Zcy（讨论）corresponding to the vertex-wise addition of grains翻译存疑 ==Zcy（讨论）


给定一个任意但固定的非沉没顶点的顺序，进行多个崩塌操作，可以通过使用拉普拉斯矩阵[math]\displaystyle{ \Delta=D-A }[/math]高效地编码多个崩塌操作（例如，在不稳定构型的稳定过程中可能发生），其中[math]\displaystyle{ D }[/math]是图的度矩阵，[math]\displaystyle{ A }[/math]是图的邻接矩阵。


删除与沉没顶点相对应的 [math]\displaystyle{ \Delta }[/math]的行和列，得到简化图拉普拉斯矩阵 [math]\displaystyle{ \Delta' }[/math]。然后，当以一个构型[math]\displaystyle{ z }[/math] 开始，并将每个顶点[math]\displaystyle{ v }[/math]进行总共[math]\displaystyle{ \mathbf{x}(v)\in\mathbb{N}_0 }[/math]次的崩塌操作时，产生[math]\displaystyle{ z-\Delta'\boldsymbol{\cdot}~\mathbf{x} }[/math]构型，其中[math]\displaystyle{ \boldsymbol{\cdot} }[/math]是收缩积。此外，如果 [math]\displaystyle{ \mathbf{x} }[/math] 对应于在给定构型[math]\displaystyle{ z }[/math]的稳定过程中每个顶点产生崩塌的次数，则

[math]\displaystyle{ z^\circ=z-\Delta'\boldsymbol{\cdot}~\mathbf{x} }[/math]


在这种情况下，[math]\displaystyle{ \mathbf{x} }[/math]被称为崩塌或[math]\displaystyle{ z }[/math]的稳定过程的里程计函数 odometer function。 

==Zcy（讨论）odometer function翻译存疑 ==Zcy（讨论）


在[math]\displaystyle{ * }[/math]运算下，递归构型的集合构成一个与约化图拉普拉斯矩阵[math]\displaystyle{ \Delta' }[/math]的核同构的阿贝尔群。对于 [math]\displaystyle{ \mathbf{Z}^{n-1}/\mathbf{Z}^{n-1}\Delta' }[/math]，其中[math]\displaystyle{ n }[/math] 表示顶点数（包括沉没顶点）。更一般地说，稳定构型集（瞬态和循环）在[math]\displaystyle{ * }[/math].运算下形成 交换幺半群Commutative monoid。这个幺半群的最小理想同构于循环构型群。The minimal ideal of this monoid is then isomorphic to the group of recurrent configurations. 


==Zcy（讨论）ideal翻译存疑==Zcy（讨论）


由循环构型形成的群，以及与之同构的群[math]\displaystyle{ \mathbf{Z}^{n-1}/\mathbf{Z}^{n-1}\Delta' }[/math]，通常称为 沙堆群Sandpile group。相同群的其它常用名称是“临界群”、“Jacobian群”或（不常见的）“Picard群”。然而，要注意的是，有些作者只把循环构型形成的组称为沙堆组，而为由[math]\displaystyle{ \mathbf{Z}^{n-1}/\mathbf{Z}^{n-1}\Delta' }[/math]（或相关的同构定义）定义的同构群保留Jacobian群或临界群的称呼。最后，一些作者使用Picard群来指代沙堆群和[math]\displaystyle{ \mathbb{Z} }[/math]的直接产物，后者出现在与沙堆模型密切相关的元胞自动机中，被称为 碎片点火或美元博弈游戏Chip firing or Dollar game。

Given the isomorphisms stated above, the order of the sandpile group is the determinant of {\displaystyle \Delta '}\Delta ', which by the matrix tree theorem is the number of spanning trees of the graph.

给定上述同构，沙堆群的顺序是[math]\displaystyle{ \Delta' }[/math]的行列式，根据矩阵树定理，它是图的生成树数目。 

==Zcy（讨论）the order of the sandpile group is the determinant of [math]\displaystyle{  \Delta' }[/math] 的翻译存疑，order是否翻译成顺序

自组织临界性

模型最初源于这样一个事实，即在晶格上的模拟中，它被吸引到了它的临界状态，此时系统的关联长度和关联时间趋于无穷大，不需要对系统参数进行任何微调。这与早期临界现象的例子形成了对比，例如固体和液体之间，或液体和气体之间的相变 phase transition，其中临界点只能通过精确调节（例如，温度）来达到。因此，在沙堆模型中，我们可以说临界性是自组织的。

一旦沙堆模型达到其临界状态，系统对扰动的响应和扰动细节之间就没有关联。一般来说，这意味着再往沙堆形成的斜坡上撒一粒沙子可能不会导致任何事情发生，或者可能导致整个沙堆形成的斜坡在大规模滑坡中崩塌。该模型还显示了1/f噪声 [math]\displaystyle{ 1/f noise }[/math]，这是自然界中许多复杂系统的共同特征。


此模型仅在两个或多个维度中显示临界现象。沙堆模型可以用一维来表示; 然而，一维沙堆模型不是演化到临界状态，而是达到最小稳定状态，其中每个格点都趋向临界坡度。


对于二维，相关共形场理论被认为是中心电荷为c = −2的辛费米子 symplectic fermion组成的。[4]

属性/性质

最小作用原理

The stabilization of chip configurations obeys a form of least action principle: each vertex topples no more than necessary in the course of the stabilization.[5]

芯片构型的稳定遵循一种“最小作用原理”的形式：每个顶点在稳定过程中不超过必要的崩塌量。 

==Zcy（讨论）chip configurations的翻译存疑==Zcy（讨论）


这可以形式化如下。如果一个崩塌序列只推倒不稳定的顶点，则称其为“合法的”，使用一组合法的顶点序列，如果它的结果是一个稳定的构型，则称其为“稳定的”。稳定沙堆的标准方法是找到一个最大的合法崩塌序列，也就是说，让崩塌序列尽可能地长。这种序列具有明显的稳定性，沙堆的可交换性质是所有这些序列都等价于倾斜序列的置换，也就是说，对于任何顶点[math]\displaystyle{ v }[/math]，在所有合法的稳定序列中[math]\displaystyle{ v }[/math]的崩塌次数都是相同的。根据最小作用原理，最小稳定序列等价于合法的（且稳定的）崩塌序列的置换。特别地，由最小稳定序列产生的构型与由最大合法序列产生的构型是相同的。


More formally, if [math]\displaystyle{ \mathbf{u} }[/math] is a vector such that [math]\displaystyle{ \mathbf{u}(v) }[/math] is the number of times the vertex [math]\displaystyle{ v }[/math] topples during the stabilization (via the toppling of unstable vertices) of a chip configuration [math]\displaystyle{ z }[/math], and [math]\displaystyle{ \mathbf{n} }[/math] is an integral vector (not necessarily non-negative) such that [math]\displaystyle{ z-\mathbf{n}\Delta' }[/math] is stable, then [math]\displaystyle{ \mathbf{u}(v) \leq \mathbf{n}(v) }[/math] for all vertices [math]\displaystyle{ v }[/math].

更正式地说，如果[math]\displaystyle{ \mathbf{u} }[/math]是一个向量， [math]\displaystyle{ \mathbf{u}（v） }[/math]是碎片构型[math]\displaystyle{ z }[/math]在稳定过程中（通过不稳定顶点的崩塌）顶点[math]\displaystyle{ v }[/math]崩塌的次数，并且[math]\displaystyle{ \mathbf{n} }[/math]是一个积分向量（不一定是非负的），使得[math]\displaystyle{ z-\mathbf{n}\Delta' }[/math]是稳定的，那么对于所有顶点[math]\displaystyle{ v }[/math]，[math]\displaystyle{ \mathbf{u}(v) \leq \mathbf{n}(v) }[/math]。 

==Zcy（讨论）[math]\displaystyle{ \mathbf{n} }[/math] is an integral vector (not necessarily non-negative)这一句话的翻译存疑。==Zcy（讨论）

缩放极限

沙堆标识在方形网格不断增加的动画。 黑色表示沙粒数为0的顶点，绿色表示沙粒数为1，紫色表示沙粒数为2，金色表示沙粒数为3。

动画显示了对应网格尺寸[math]\displaystyle{ N\geq 1 }[/math]不断增大，不同大小的[math]\displaystyle{ N\times N }[/math]正方形网格上的沙堆群标识的重复构型，从而重新缩放构型，使其最终具有相同的物理尺度。从视觉上看，更大网格上的标识似乎变得越来越详细，并且“收敛到一个连续的图像”。从数学上讲，这表明基于弱收敛的概念（或其他一些广义的收敛概念），正方形网格上沙堆恒等式存在缩放极限。事实上，Wesley-Pegden和Charles-Smart已经证明了循环沙堆结构缩放极限的存在性。[6] .[7]在与Lionel Levine的进一步合作中，他们使用缩放极限来解释方形网格上沙堆的分形结构。[8] 

==Zcy（讨论）whereby the configurations are rescaled to always have the same physical dimension翻译存疑 ==Zcy（讨论）

归纳与相关模型

有向图上的沙堆模型

3000万粒沙粒落在无限方形网格的一个位置上，然后按照沙堆模型的规则产生崩塌。白色表示0颗沙粒位置，绿色表示1颗，紫色表示2颗，金色表示3颗。图中框的大小是3967×3967。


沙堆模型可以推广到无限网格中。这种归纳法的一个挑战是，一般来说，不能保证每次雪崩最终都会停止。因此，一些一般化方法只考虑了构型的稳定性，因为这一点是能保证的。


在（无限）方格上有一个相当流行的模型，其位置[math]\displaystyle{ (x,y)\in\mathbb{Z}^2 }[/math]定义如下：


从\mathbf{z}</math>中有限值[math]\displaystyle{ z(x,y)\in \mathbf{Z} }[/math]的一些非负构型开始，这意味着

[math]\displaystyle{ \sum_{x,y}z(x,y)\lt \infty. }[/math]


任何位置[math]\displaystyle{ (x,y) }[/math]有

[math]\displaystyle{ z(x,y)\geq 4 }[/math]

就是“不稳定的”，并且会产生崩塌，将它位置上的沙粒分发给它的四个邻居：

[math]\displaystyle{ z(x,y) \rightarrow z(x,y) - 4, }[/math]
[math]\displaystyle{ z( x \pm 1, y) \rightarrow z( x \pm 1, y) + 1, }[/math]
[math]\displaystyle{ z(x, y \pm 1) \rightarrow z( x, y \pm 1 ) + 1. }[/math]


由于初始构型是有限的，这一过程必然会随着沙粒向外散布终止。


给出了该模型的一个常见特例，即除原点外的所有顶点的初始构型都为零。如果原点放置大量沙粒，松弛后的构型形成分形图案(见图)。当初始颗粒数趋于无穷时，重缩放的稳定构型收敛到唯一极限。[7][8]

有向图上的沙堆模型

沙堆模型可以推广到任意有向多重图。规则是任何顶点[math]\displaystyle{ v }[/math]有

[math]\displaystyle{ z(v)\geq \deg^{+}(v) }[/math]

则是不稳定的；崩塌将它的碎片沿着边的输出方向，分发给它的邻居：

[math]\displaystyle{ z(v) \rightarrow z(v) - \deg^{+}(v) + \deg(v,v) }[/math]

并且，对于每个[math]\displaystyle{ u\neq v }[/math]：

[math]\displaystyle{ z(u) \rightarrow z(u) + \deg(v,u) }[/math]

其中，[math]\displaystyle{ \deg(v,u) }[/math]是从[math]\displaystyle{ v }[/math]到[math]\displaystyle{ u }[/math]的边数。


在这种情况下，拉普拉斯矩阵是不对称的。如果我们指定一个沉没顶点[math]\displaystyle{ s }[/math]，使得每一个顶点都有一条到[math]\displaystyle{ s }[/math]的路径，那么有限图上的稳定操作是定义良好的，并且沙堆群可以像之前一样被写出来

[math]\displaystyle{ \mathbf{Z}^{n-1}/\mathbf{Z}^{n-1}\Delta' }[/math]

The order of the sandpile group is again the determinant of [math]\displaystyle{ \Delta' }[/math], which by the general version of the matrix tree theorem is the number of oriented spanning trees rooted at the sink.

沙堆群的顺序又是[math]\displaystyle{ \Delta' }[/math]的行列式，根据矩阵树定理的一般版本，它是根在沉没顶点的有向spanning tree生成树的数目。 

==Zcy（讨论）The order of the sandpile group翻译存疑。==Zcy（讨论）

扩展沙堆模型

由谐波函数H=x*y在255x255方格网格上引起的沙堆动力学。


为了更好地理解不同有限凸网格[math]\displaystyle{ \Gamma\subset\mathbb{Z}^2 }[/math]的沙堆群的结构，Lang和Shkolnikov在2019年提出了“扩展沙堆模型”。[9]扩展沙堆模型的定义与“通常的沙堆模型”几乎完全相同（即原始的Bak–Tang–Wiesenfeld模型[1]），除了网格边界[math]\displaystyle{ \partial\Gamma }[/math]的顶点现在允许放置非负实数的沙粒。相比之下，网格内部的顶点仍然只允许放置整数个粒子。崩塌规则保持不变，即假设当沙粒数达到或超过4时，内部顶点和边界顶点都变得不稳定并发生崩塌。


Also the recurrent configurations of the extended sandpile model form an abelian group, referred to as the extended sandpile group, of which the usual sandpile group is a discrete subgroup. Different to the usual sandpile group, the extended sandpile group is however a continuous Lie group. Since it is generated by only adding grains of sand to the boundary [math]\displaystyle{ \partial\Gamma }[/math] of the grid, the extended sandpile group furthermore has the topology of a torus of dimension [math]\displaystyle{ |\partial\Gamma| }[/math] and a volume given by the order of the usual sandpile group.[9]

扩展沙堆模型的递归构型也形成了一个阿贝尔群，称为“扩展沙堆群”，通常的扩展沙堆群是一个离散子群。与通常的沙堆群不同，扩展沙堆群是一个连续的李群。只因为它是由添加沙粒到网格的边界[math]\displaystyle{ \partial\Gamma }[/math]上形成的，扩展后的沙堆群还具有维度[math]\displaystyle{ |\partial\Gamma| }[/math]的环面拓扑结构，并且按通常沙堆组的顺序给出的体积。[9] 

==Zcy（讨论）Since it is generated by only adding grains of sand to the boundary [math]\displaystyle{ \partial\Gamma }[/math] of the grid翻译存疑。a volume given by the order of the usual sandpile group.翻译存疑。==Zcy（讨论）

Of specific interest is the question how the recurrent configurations dynamically change along the continuous geodesics of this torus passing through the identity. This question leads to the definition of the sandpile dynamics

特别感兴趣的问题是循环构型如何通过恒等式，沿着这个环面的连续测地线动态变化的问题。这个问题引出了沙堆动力学的定义 

==Zcy（讨论）identity在整篇文章中的翻译需进行统一，如何翻译？？同一性，恒等式？？==Zcy（讨论）
[math]\displaystyle{ D_H(t)=(I-t\Delta H)^\circ }[/math] （扩展沙堆模型）

respectively

[math]\displaystyle{ \tilde{D}_H(t)=(I+\lfloor-t\Delta H\rfloor)^\circ }[/math] （普通沙堆模型）

This proposes a natural renormalization for the extended and usual sandpile groups, meaning a mapping of recurrent configurations on a given grid to recurrent configurations on a sub-grid. Informaly, this renormalization simply maps configurations appearing at a given time [math]\displaystyle{ t }[/math] in the sandpile dynamics induced by some harmonic function [math]\displaystyle{ H }[/math] on the larger grid to the corresponding configurations which appear at the same time in the sandpile dynamics induced by the restriction of [math]\displaystyle{ H }[/math] to the respective sub-grid.[9]

由整值调和函数[math]\displaystyle{ H }[/math]在时间[math]\displaystyle{ t\in\mathbb{R}\setminus\mathbb{Z} }[/math]，沙堆群的同一性[math]\displaystyle{ I }[/math]和底函数[math]\displaystyle{ \lfloor.\rfloor }[/math]导出的。[9]对于低阶多项式调和函数，沙堆动力学的特征是组成沙堆恒等式的斑块的光滑变换和明显守恒。例如，由[math]\displaystyle{ H=xy }[/math] 诱导的谐波动力学类似于动画中可视化的主对角线上恒等式的“平滑拉伸”。进一步推测了由相同的谐函数在不同尺寸的正方形网格上引起的动力学构型的弱收敛，这意味着可能存在缩放极限。[9]这为扩展的和普通的沙堆组提出了一个自然的重归一化，这意味着在给定网格上的重复构型映射到子网格上的重复构型。非正式地，重归一化简单地映射了沙堆动力学中给定时间[math]\displaystyle{ t }[/math]时的构型，动力学由大型网格上的谐波函数[math]\displaystyle{ H }[/math]导出到相应的构型，这种构型在[math]\displaystyle{ H }[/math]限制到各自子网格的沙堆动力学中时同时出现。[9] 

==Zcy（讨论）上面的整个长句翻译需要重新审校Informaly, this renormalization simply maps==Zcy（讨论）

可分割的沙堆

Levine和Peres在2008年提出了一个与之密切相关的模型，即所谓的“可分割的沙堆模型”。[10]与每个位置[math]\displaystyle{ x }[/math]上的沙粒数量为离散数不同，有一个实数[math]\displaystyle{ s(x) }[/math]代表位置的总质量。如果这个质量是负的，我们就可以把它理解为一个空洞。当一个位置上的质量大于1时，就会发生崩塌; 它将多余的部分均匀地分发给它的邻居，这就导致了如果一个位置在[math]\displaystyle{ t }[/math]的时刻质量是1，它在以后的所有时间质量都是1。

参考文献

  1. ↑ 1.0 1.1 Bak, P.; Tang, C.; Wiesenfeld, K. (1987). "Self-organized criticality: an explanation of 1/ƒ noise". Physical Review Letters. 59 (4): 381–384. Bibcode:1987PhRvL..59..381B. doi:10.1103/PhysRevLett.59.381. PMID 10035754.
  2. ↑ Holroyd, A.; Levine, L.; Mészáros, K.; Peres, Y.; Propp, J.; Wilson, B. (2008). Chip-Firing and Rotor-Routing on Directed Graphs. 60. pp. 331–364. arXiv:0801.3306. Bibcode 1987PhRvL..59..381B. doi:10.1007/978-3-7643-8786-0_17. ISBN 978-3-7643-8785-3. 
  3. ↑ Biggs, Norman L. (25 June 1997). "Chip-Firing and the Critical Group of a Graph" (PDF). Journal of Algebraic Combinatorics: 25–45. Retrieved 10 May 2014.
  4. ↑ S. Moghimi-Araghi; M. A. Rajabpour; S. Rouhani (2004). "Abelian Sandpile Model: a Conformal Field Theory Point of View". Nuclear Physics B. 718 (3): 362–370. arXiv:cond-mat/0410434. Bibcode:2005NuPhB.718..362M. doi:10.1016/j.nuclphysb.2005.04.002. S2CID 16233977.
  5. ↑ Fey, A.; Levine, L.; Peres, Y. (2010). "Growth Rates and Explosions in Sandpiles". Journal of Statistical Physics. 138 (1–3): 143–159. arXiv:0901.3805. Bibcode:2010JSP...138..143F. doi:10.1007/s10955-009-9899-6. ISSN 0022-4715. S2CID 7180488.
  6. ↑ Pegden, Wesley; Smart, Charles (2017). "Stability of patterns in the Abelian sandpile". arXiv:1708.09432 [math.AP].
  7. ↑ 7.0 7.1 Pegden, Wesley; Smart, Charles (2013). "Convergence of the Abelian sandpile". Duke Mathematical Journal. 162 (4): 627–642. arXiv:1105.0111. doi:10.1215/00127094-2079677. S2CID 13027232.
  8. ↑ 8.0 8.1 Levine, Lionel; Pegden, Wesley (2016). "Apollonian structure in the Abelian sandpile". Geometric and Functional Analysis. 26 (1): 306–336. doi:10.1007/s00039-016-0358-7. hdl:1721.1/106972. S2CID 119626417.
  9. ↑ 9.0 9.1 9.2 9.3 9.4 9.5 9.6 Lang, Moritz; Shkolnikov, Mikhail (2019-02-19). "Harmonic dynamics of the abelian sandpile". Proceedings of the National Academy of Sciences (in English). 116 (8): 2821–2830. doi:10.1073/pnas.1812015116. ISSN 0027-8424. PMC 6386721. PMID 30728300.
  10. ↑ Levine, Lionel; Peres, Yuval (2008-10-29). "Strong Spherical Asymptotics for Rotor-Router Aggregation and the Divisible Sandpile". Potential Analysis (in English). 30 (1): 1–27. arXiv:0704.0688. doi:10.1007/s11118-008-9104-6. ISSN 0926-2601. S2CID 2227479.

进一步阅读

拓展链接

编者推荐

复杂性与临界现象课程链接


复杂性与临界现象

本课程中，北京师范大学系统科学院长陈晓松介绍了复杂系统中的平衡态临界现象以及非平衡与复杂系统临界现象。


NetLogo NetLogo多主体建模课程

北京师范大学系统科学学院张江教授通过数个案例教会大家如何去动手搭建一个多主体仿真模型，以及如何利用NetLogo去实现。可以教你如何用Netlogo程序来做人工社会的仿真和模拟。


包括以下模型：
  • 用“生命游戏”认识Patch
  • 从Langton的蚂蚁看Turtle与Patch的交互
  • 从羊-草生态系统深入Turtle与Plot画图
  • 玩具经济模型与 Turtles间的互动
  • 从玩具经济模型学文件导出与复杂曲线绘制
  • 从玩具经济模型学习如何使用行为空间做实验
  • 透过人工鸟群Boid模型学习List的使用
  • 网络上病毒传播的SIR模型 用Links建模网络动力学
  • 重访羊-草模型与系统动力学建模


本中文词条由水流心不竞翻译，由Zcy审校，思无涯咿呀咿呀编辑欢迎在讨论页面留言。

本词条内容翻译自 wikipedia.org，遵守 CC3.0协议。

取自“https://wiki.swarma.org/index.php?title=沙堆模型&oldid=18457”