El modelo abeliano de pilas de arena (ASM) es el nombre más popular del modelo original de Bak-Tang-Wiesenfeld (BTW). El modelo BTW fue el primer ejemplo descubierto de un sistema dinámico que muestra una criticidad autoorganizada . Fue presentado por Per Bak , Chao Tang y Kurt Wiesenfeld en un artículo de 1987. [1]
Tres años más tarde, Deepak Dhar descubrió que el modelo de pilas de arena, por cierto, sigue la dinámica abeliana y, por lo tanto, se refirió a este modelo como modelo de pilas de arena abelianas. [2]
El modelo es un autómata celular . En su formulación original, cada sitio en una cuadrícula finita tiene un valor asociado que corresponde a la pendiente del pilote. Esta pendiente se acumula a medida que "granos de arena" (o "astillas") se colocan aleatoriamente sobre la pila, hasta que la pendiente excede un valor umbral específico, momento en el que ese sitio colapsa transfiriendo arena a los sitios adyacentes, aumentando su pendiente. Bak, Tang y Wiesenfeld consideraron el proceso de colocación aleatoria sucesiva de granos de arena en la rejilla; cada una de estas colocaciones de arena en un sitio en particular puede no tener ningún efecto o puede causar una reacción en cascada que afectará a muchos sitios.
Dhar ha demostrado que la configuración estable final del montón de arena después de que termina la avalancha es independiente de la secuencia precisa de derrumbes que se sigue durante la avalancha. Como consecuencia directa de este hecho, se muestra que si se añaden dos granos de arena a la configuración estable en dos órdenes diferentes, por ejemplo, primero en el sitio A y luego en el sitio B, y primero en B y luego en A, el resultado final. La configuración estable de los granos de arena resulta ser exactamente la misma. Cuando se agrega un grano de arena a una configuración estable de un montón de arena, se produce una avalancha que finalmente se detiene y conduce a otra configuración estable. Dhar propuso que la adición de un grano de arena puede considerarse como un operador: cuando actúa sobre una configuración estable, produce otra configuración estable. Dhar demostró que todos estos operadores de suma forman un grupo abeliano, de ahí el nombre de modelo de pila de arena abeliano. [3] [4] Desde entonces, el modelo se ha estudiado en la red infinita, en otras redes (no cuadradas) y en gráficos arbitrarios (incluidos los multigrafos dirigidos ). [5] Está estrechamente relacionado con el juego del dólar , una variante del juego de disparo de fichas introducido por Biggs. [6]
El modelo de pila de arena es un autómata celular originalmente definido en una cuadrícula rectangular ( tablero de ajedrez ) de la red cuadrada estándar . A cada vértice ( lado , campo ) de la cuadrícula, asociamos un valor ( granos de arena , pendiente , partículas ) , a lo que nos referimos como configuración (inicial) del montón de arena.
La dinámica del autómata en la iteración se define entonces de la siguiente manera:
El derrumbe de varios vértices durante una iteración se denomina avalancha . Se garantiza que cada avalancha finalmente se detendrá, es decir, después de un número finito de derribos se alcanza alguna configuración estable tal que el autómata esté bien definido. Además, aunque a menudo habrá muchas opciones posibles para el orden en el que derribar los vértices, la configuración estable final no depende del orden elegido; éste es un sentido en el que el montón de arena es abeliano . De manera similar, el número de veces que cada vértice cae durante cada iteración también es independiente de la elección del orden de caída.
Para generalizar el modelo de pila de arena desde la cuadrícula rectangular de la red cuadrada estándar a un multigrafo finito arbitrario no dirigido , se especifica un vértice especial llamado sumidero al que no se le permite derribar. Una configuración (estado) del modelo es entonces una función que cuenta el número no negativo de granos en cada vértice no sumidero. Un vértice no sumidero con
es inestable; se puede derribar, lo que envía uno de sus granos a cada uno de sus vecinos (no sumidero):
Luego, el autómata celular progresa como antes, es decir, agregando, en cada iteración, una partícula a un vértice no sumidero elegido al azar y derribándose hasta que todos los vértices estén estables.
La definición del modelo de pila de arena dada anteriormente para cuadrículas rectangulares finitas de la red cuadrada estándar puede verse entonces como un caso especial de esta definición: considere el gráfico que se obtiene agregando un vértice adicional, el sumidero, y dibujando aristas adicionales. desde el sumidero hasta cada vértice límite de tal que el grado de cada vértice que no es sumidero es cuatro. De esta manera, también se pueden definir modelos de pilas de arena en cuadrículas no rectangulares de la red cuadrada estándar (o de cualquier otra red): Intersecta algún subconjunto acotado de con . Contrae cada arista cuyos dos extremos no estén en . El único vértice restante fuera de entonces constituye el sumidero del gráfico de pila de arena resultante.
En la dinámica del autómata de pila de arena definida anteriormente, algunas configuraciones estables ( para todos ) aparecen con una frecuencia infinita, mientras que otras sólo pueden aparecer un número finito de veces (si es que aparecen). Las primeras se denominan configuraciones recurrentes , mientras que las segundas se denominan configuraciones transitorias . Por lo tanto, las configuraciones recurrentes consisten en todas las configuraciones estables no negativas que se pueden alcanzar desde cualquier otra configuración estable agregando repetidamente granos de arena a los vértices y derribándolos. Es fácil ver que la configuración mínimamente estable , donde cada vértice lleva granos de arena, es alcanzable desde cualquier otra configuración estable (agregue granos a cada vértice). Por lo tanto, de manera equivalente, las configuraciones recurrentes son exactamente aquellas configuraciones que se pueden alcanzar a partir de la configuración mínimamente estable agregando solo granos de arena y estabilizando.
No todas las configuraciones estables no negativas son recurrentes. Por ejemplo, en cada modelo de montón de arena en un gráfico que consta de al menos dos vértices no sumideros conectados, cada configuración estable donde ambos vértices llevan cero granos de arena no es recurrente. Para probar esto, primero observe que la adición de granos de arena solo puede aumentar el número total de granos transportados por los dos vértices juntos. Llegar a una configuración en la que ambos vértices llevan cero partículas desde una configuración en la que este no es el caso implica necesariamente pasos en los que al menos uno de los dos vértices se derriba. Considere el último de estos pasos. En este paso, uno de los dos vértices debe caer en último lugar. Dado que el derrumbe transfiere un grano de arena a cada vértice vecino, esto implica que el número total de granos transportados por ambos vértices juntos no puede ser inferior a uno, lo que concluye la prueba.
Dada una configuración , para todos , derribar los vértices inestables no sumidero en un gráfico finito conectado hasta que no quede ningún vértice inestable no sumidero conduce a una configuración estable única , que se llama estabilización de . Dadas dos configuraciones estables y , podemos definir la operación , correspondiente a la adición de granos en los vértices seguida de la estabilización del montón de arena resultante.
Dado un orden arbitrario pero fijo de los vértices que no son sumidero, múltiples operaciones de derrumbe, que pueden ocurrir, por ejemplo, durante la estabilización de una configuración inestable, se pueden codificar eficientemente usando el gráfico laplaciano , donde es la matriz de grados y es la matriz de adyacencia de la gráfica. Al eliminar la fila y la columna correspondientes al sumidero se obtiene el gráfico reducido Laplaciano . Luego, al comenzar con una configuración y derribar cada vértice un total de veces se obtiene la configuración , donde está el producto de contracción. Además, si corresponde al número de veces que cada vértice se derriba durante la estabilización de una configuración determinada , entonces
En este caso, se denomina función de derribo o de odómetro (de estabilización de ).
Bajo la operación , el conjunto de configuraciones recurrentes forma un grupo abeliano isomorfo al cokernel del grafo reducido laplaciano , es decir, a , donde denota el número de vértices (incluido el sumidero). De manera más general, el conjunto de configuraciones estables (transitorias y recurrentes) forma un monoide conmutativo bajo la operación . El ideal mínimo de este monoide es entonces isomorfo al grupo de configuraciones recurrentes.
El grupo formado por las configuraciones recurrentes, así como el grupo al que la primera es isomorfa, se denomina más comúnmente grupo de pilas de arena . Otros nombres comunes para un mismo grupo son grupo crítico , grupo jacobiano o (con menor frecuencia) grupo Picard . Tenga en cuenta, sin embargo, que algunos autores sólo denotan el grupo formado por las configuraciones recurrentes como grupo de pilas de arena, reservando el nombre de grupo jacobiano o grupo crítico para el grupo (isomorfo) definido por (o para definiciones isomórficas relacionadas). Finalmente, algunos autores utilizan el nombre grupo Picard para referirse al producto directo del grupo sandpile y , que naturalmente aparece en un autómata celular muy relacionado con el modelo sandpile, denominado disparo de fichas o juego del dólar.
Dados los isomorfismos indicados anteriormente, el orden del grupo de pilas de arena es el determinante de , que según el teorema del árbol matricial es el número de árboles generadores del gráfico.
El interés original detrás del modelo surgió del hecho de que en simulaciones en redes, se siente atraído por su estado crítico , en cuyo punto la longitud de correlación del sistema y el tiempo de correlación del sistema llegan al infinito, sin ningún ajuste fino de un parámetro del sistema. Esto contrasta con ejemplos anteriores de fenómenos críticos, como las transiciones de fase entre sólido y líquido, o líquido y gas, donde el punto crítico sólo puede alcanzarse mediante un ajuste preciso (por ejemplo, de la temperatura). Por tanto, en el modelo de montón de arena podemos decir que la criticidad es autoorganizada .
Una vez que el modelo de pila de arena alcanza su estado crítico, no existe correlación entre la respuesta del sistema a una perturbación y los detalles de una perturbación. Generalmente, esto significa que dejar caer otro grano de arena sobre la pila puede causar que no pase nada, o puede provocar que toda la pila colapse en un deslizamiento masivo. El modelo también muestra ruido 1/ ƒ , una característica común a muchos sistemas complejos en la naturaleza.
Este modelo sólo muestra un comportamiento crítico en dos o más dimensiones. El modelo de pila de arena se puede expresar en 1D; sin embargo, en lugar de evolucionar a su estado crítico, el modelo de pila de arena 1D alcanza un estado mínimamente estable donde cada sitio de la red se dirige hacia la pendiente crítica.
Para dos dimensiones, se ha planteado la hipótesis de que la teoría de campos conforme asociada consta de fermiones simplécticos con una carga central c = −2. [7]
La estabilización de las configuraciones de chips obedece a una especie de principio de mínima acción : cada vértice no se cae más de lo necesario durante la estabilización. [8] Esto puede formalizarse de la siguiente manera. Llame legal a una secuencia de derribos si solo derriba vértices inestables y estabilizadora si da como resultado una configuración estable. La forma estándar de estabilizar el montón de arena es encontrar una secuencia legal máxima; es decir, derribándose tanto como sea posible. Esta secuencia es obviamente estabilizadora, y la propiedad abeliana del montón de arena es que todas esas secuencias son equivalentes hasta la permutación del orden de derrumbe; es decir, para cualquier vértice , el número de caídas es el mismo en todas las secuencias estabilizadoras legales. Según el principio de acción mínima, una secuencia estabilizadora mínima también equivale, hasta la permutación de la orden de derrocamiento, a una secuencia legal (y aún estabilizadora). En particular, la configuración resultante de una secuencia estabilizadora mínima es la misma que resulta de una secuencia legal máxima.
Más formalmente, si es un vector tal que es el número de veces que el vértice se cae durante la estabilización (mediante el derrumbe de vértices inestables) de una configuración de chip , y es un vector integral (no necesariamente no negativo) tal que es estable, luego para todos los vértices .
La animación muestra la configuración recurrente correspondiente a la identidad del grupo de pilas de arena en diferentes cuadrículas cuadradas de tamaños crecientes , mediante las cuales las configuraciones se reescalan para tener siempre la misma dimensión física. Visualmente, las identidades en cuadrículas más grandes parecen volverse cada vez más detalladas y "convergir en una imagen continua". Matemáticamente, esto sugiere la existencia de límites de escala de la identidad del montón de arena en cuadrículas cuadradas basadas en la noción de convergencia débil* (o alguna otra noción generalizada de convergencia). De hecho, Wesley Pegden y Charles Smart han demostrado la existencia de límites de escala de configuraciones recurrentes de montones de arena. [9] [10] En un trabajo conjunto posterior con Lionel Levine, utilizan el límite de escala para explicar la estructura fractal del montón de arena en cuadrículas cuadradas. [11] Otro límite de escala, cuando las relajaciones de una perturbación del estado estable máximo convergen a una imagen definida por curvas tropicales , se establece en las obras de Nikita Kalinin y Mikhail Shkolnikov. [12]
Los montones de arena abelianos en tres o más dimensiones se pueden utilizar para simular una máquina de Turing y, por lo tanto, son Turing completos . [13]
Existen varias generalizaciones del modelo de pila de arena a cuadrículas infinitas. Un desafío en tales generalizaciones es que, en general, ya no está garantizado que cada avalancha eventualmente se detenga. Por lo tanto, varias generalizaciones sólo consideran la estabilización de configuraciones para las que esto puede estar garantizado.
Un modelo bastante popular en la red cuadrada (infinita) con sitios se define de la siguiente manera:
Comience con alguna configuración no negativa de valores que sea finita, es decir
Cualquier sitio con
es inestable y puede caerse (o dispararse ), enviando uno de sus chips a cada uno de sus cuatro vecinos:
Dado que la configuración inicial es finita, se garantiza que el proceso terminará y los granos se dispersarán hacia afuera.
Un caso especial popular de este modelo se da cuando la configuración inicial es cero para todos los vértices excepto el origen. Si el origen contiene una gran cantidad de granos de arena, la configuración después de la relajación forma patrones fractales (ver figura). Al dejar que el número inicial de granos en el origen llegue al infinito, se demostró que las configuraciones estabilizadas reescaladas convergen a un límite único. [10] [11]
El modelo de pila de arena se puede generalizar a multigrafos dirigidos arbitrarios. Las reglas son que cualquier vértice con
es inestable; al caer nuevamente envía chips a cada uno de sus vecinos, uno a lo largo de cada borde de salida:
y, para cada uno :
¿Dónde está el número de aristas desde hasta ?
En este caso la matriz laplaciana no es simétrica. Si especificamos un sumidero tal que haya un camino desde cada dos vértices hasta , entonces la operación de estabilización en gráficos finitos está bien definida y el grupo de pilas de arena se puede escribir
como antes.
El orden del grupo de pilas de arena es nuevamente el determinante de , que según la versión general del teorema del árbol matricial es el número de árboles de expansión orientados enraizados en el sumidero.
Para comprender mejor la estructura del grupo de pilas de arena para diferentes cuadrículas finitas convexas de la red cuadrada estándar , Lang y Shkolnikov introdujeron el modelo de pila de arena extendido en 2019. [14] El modelo de pila de arena extendido se define casi exactamente igual que el modelo de pila de arena habitual ( es decir, el modelo original de Bak-Tang-Wiesenfeld [1] ), excepto que ahora se permite que los vértices en el límite de la cuadrícula lleven un número real no negativo de granos. Por el contrario, a los vértices en el interior de la cuadrícula solo se les permite transportar números enteros de granos. Las reglas de caída permanecen sin cambios, es decir, se supone que tanto el vértice interior como el límite se vuelven inestables y se caen si el número de granos alcanza o excede cuatro.
Además, las configuraciones recurrentes del modelo de pilotes de arena extendidos forman un grupo abeliano, denominado grupo de pilotes de arena extendido , del cual el grupo de pilotes de arena habitual es un subgrupo discreto . Sin embargo, a diferencia del grupo de montones de arena habitual, el grupo de montones de arena extendido es un grupo de Lie continuo . Dado que se genera añadiendo únicamente granos de arena al límite de la cuadrícula, el grupo de pilas de arena extendido tiene además la topología de un toro de dimensión y un volumen dado por el orden del grupo de pilas de arena habitual. [14]
De interés específico es la cuestión de cómo las configuraciones recurrentes cambian dinámicamente a lo largo de las geodésicas continuas de este toro que pasa a través de la identidad. Esta pregunta lleva a la definición de la dinámica de los montones de arena.
respectivamente
inducido por la función armónica de valor entero en el tiempo , con la identidad del grupo de pilas de arena y la función de piso. [14] Para funciones armónicas polinómicas de bajo orden, la dinámica del montón de arena se caracteriza por la suave transformación y la aparente conservación de los parches que constituyen la identidad del montón de arena. Por ejemplo, la dinámica armónica inducida por se asemeja al "estiramiento suave" de la identidad a lo largo de las diagonales principales visualizadas en la animación. Además, se conjeturó que las configuraciones que aparecen en la dinámica inducida por la misma función armónica en cuadrículas cuadradas de diferentes tamaños convergen débilmente, lo que significa que supuestamente existen límites de escala para ellas. [14] Esto propone una renormalización natural para los grupos de pilas de arena extendidos y habituales, es decir, un mapeo de configuraciones recurrentes en una red determinada a configuraciones recurrentes en una subred. De manera informal, esta renormalización simplemente mapea las configuraciones que aparecen en un momento dado en la dinámica del montón de arena inducida por alguna función armónica en la cuadrícula más grande a las configuraciones correspondientes que aparecen al mismo tiempo en la dinámica del montón de arena inducida por la restricción de a la respectiva subred. . [14]
Un modelo fuertemente relacionado es el llamado modelo de pila de arena divisible , introducido por Levine y Peres en 2008, [15] en el que, en lugar de un número discreto de partículas en cada sitio , hay un número real que representa la cantidad de masa en el sitio. En caso de que dicha masa sea negativa, se puede entender como un agujero. El derrumbe ocurre siempre que un sitio tiene una masa mayor que 1; derriba el exceso de manera uniforme entre sus vecinos, lo que da como resultado la situación de que, si un sitio está lleno en un momento , estará lleno en todos los momentos posteriores.
El montón de arena Bak-Tang-Wiesenfeld se mencionó en el episodio "Rampage" de Numb3rs , mientras el matemático Charlie Eppes explica a sus colegas una solución a una investigación criminal.
El juego de computadora Hexplode se basa en el modelo de montón de arena abeliano en una cuadrícula hexagonal finita donde, en lugar de colocar granos al azar, los jugadores colocan los granos.