El modelo de pila de arena abeliano (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 criticidad autoorganizada . Fue introducido 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 pila de arena BTW sigue una dinámica abeliana y, por lo tanto, se refirió a este modelo como el modelo de pila de arena abeliano. [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 de la pila. Esta pendiente se acumula a medida que se colocan aleatoriamente "granos de arena" (o "astillas") 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 un proceso de colocación aleatoria sucesiva de granos de arena en la cuadrícula; cada colocación de arena en un sitio en particular puede no tener 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 derribos que se siguen durante la avalancha. Como consecuencia directa de este hecho, se muestra que si se agregan 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, la configuración estable final de los granos de arena resulta ser exactamente la misma. Cuando se agrega un grano de arena a una configuración estable de montón de arena, da como resultado 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 un operador, cuando actúa sobre una configuración estable, produce otra configuración estable. Dhar mostró que todos estos operadores de adición forman un grupo abeliano, de ahí el nombre de modelo de montón de arena abeliano. [3] [4] Desde entonces, el modelo se ha estudiado en la red infinita, en otras redes (no cuadradas) y en grafos arbitrarios (incluidos los multigrafos dirigidos ). [5] Está estrechamente relacionado con el juego del dólar , una variante del juego de fichas introducido por Biggs. [6]
El modelo de pila de arena es un autómata celular definido originalmente en una cuadrícula rectangular ( damero ) 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 ) , al que denominamos configuración (inicial) de la pila de arena.
La dinámica del autómata en la iteración se define entonces de la siguiente manera:
El derribo de varios vértices durante una iteración se denomina avalancha . Se garantiza que cada avalancha se detendrá eventualmente, es decir, después de un número finito de derribos se alcanza una 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 vértices, la configuración estable final no depende del orden elegido; este 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 se derrumba durante cada iteración también es independiente de la elección del orden de derribo.
Para generalizar el modelo de pila de arena de 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 que no puede volcarse. Una configuración (estado) del modelo es entonces una función que cuenta el número no negativo de granos en cada vértice que no es sumidero. Un vértice que no es sumidero con
es inestable; puede volcarse, lo que envía uno de sus granos a cada uno de sus vecinos (no hundidos):
El autómata celular luego progresa como antes, es decir, añadiendo, en cada iteración, una partícula a un vértice no hundido elegido aleatoriamente y cayendo hasta que todos los vértices sean 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 de añadiendo 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 sea sumidero de sea cuatro. De esta manera, también se pueden definir modelos de pila de arena en cuadrículas no rectangulares de la red cuadrada estándar (o de cualquier otra red): Intersecar algún subconjunto acotado de con . Contraer cada arista de cuyos dos puntos finales 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 definido anteriormente, algunas configuraciones estables ( para todos ) aparecen infinitamente a menudo, mientras que otras solo pueden aparecer un número finito de veces (si es que aparecen). Las primeras se denominan configuraciones recurrentes , mientras que las últimas se denominan configuraciones transitorias . Las configuraciones recurrentes consisten, por tanto, en todas las configuraciones estables no negativas a las que se puede llegar desde cualquier otra configuración estable añadiendo repetidamente granos de arena a los vértices y derribando. 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 (añadir granos a cada vértice). Por tanto, de forma equivalente, las configuraciones recurrentes son exactamente aquellas configuraciones a las que se puede llegar desde la configuración mínimamente estable añadiendo únicamente granos de arena y estabilizándose.
No todas las configuraciones estables no negativas son recurrentes. Por ejemplo, en cada modelo de pila de arena en un grafo que consta de al menos dos vértices no sumideros conectados, cada configuración estable donde ambos vértices llevan cero granos de arena es no recurrente. Para probar esto, primero note que la adición de granos de arena solo puede aumentar el número total de granos llevados por los dos vértices juntos. Para llegar a una configuración donde ambos vértices llevan cero partículas a partir de una configuración donde este no es el caso, entonces, necesariamente implica pasos donde al menos uno de los dos vértices es derribado. Considere el último de estos pasos. En este paso, uno de los dos vértices tiene que derribarse último. Dado que el derribo transfiere un grano de arena a cada vértice vecino, esto implica que el número total de granos llevados por ambos vértices juntos no puede ser menor que uno, lo que concluye la prueba.
Dada una configuración , para todos , el derribo de vértices inestables no sumideros en un grafo conexo finito hasta que no quede ningún vértice inestable no sumidero conduce a una única configuración estable , que se denomina estabilización de . Dadas dos configuraciones estables y , podemos definir la operación , correspondiente a la adición de granos por vértice 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, las operaciones de vuelco múltiples, que pueden ocurrir, por ejemplo, durante la estabilización de una configuración inestable, se pueden codificar de manera eficiente utilizando el laplaciano de grafos , donde es la matriz de grados y es la matriz de adyacencia del grafo. Al eliminar la fila y la columna de correspondientes con el sumidero se obtiene el laplaciano de grafos reducido . Luego, al comenzar con una configuración y derribar cada vértice un total de veces se obtiene la configuración , donde es el producto de contracción. Además, si corresponde al número de veces que se derriba cada vértice durante la estabilización de una configuración dada , entonces
En este caso se habla de función de vuelco o de odómetro (de estabilización de ).
Bajo la operación , el conjunto de configuraciones recurrentes forma un grupo abeliano isomorfo al conúcleo del grafo reducido Laplaciano , es decir a , por lo que 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 el primero es isomorfo, se conoce más comúnmente como el grupo de pila de arena . Otros nombres comunes para el mismo grupo son grupo crítico , grupo jacobiano o (menos a menudo) grupo de Picard . Nótese, sin embargo, que algunos autores solo denotan al grupo formado por las configuraciones recurrentes como el grupo de pila de arena, mientras que reservan el nombre de grupo jacobiano o grupo crítico para el grupo (isomorfo) definido por (o para definiciones isomorfas relacionadas). Finalmente, algunos autores usan el nombre de grupo de Picard para referirse al producto directo del grupo de pila de arena y , que aparece naturalmente en un autómata celular estrechamente relacionado con el modelo de pila de arena, conocido como el juego de disparo de chips o dólar.
Dados los isomorfismos establecidos 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 de expansión del gráfico.
El interés original detrás del modelo surgió del hecho de que en simulaciones en redes, es atraído a su estado crítico , en cuyo punto la longitud de correlación del sistema y el tiempo de correlación del sistema tienden 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 solo se puede alcanzar mediante un ajuste preciso (por ejemplo, de la temperatura). Por lo tanto, en el modelo de pila 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 la misma. En general, esto significa que dejar caer otro grano de arena sobre la pila puede no provocar que pase nada o puede provocar que toda la pila se derrumbe en un deslizamiento masivo. El modelo también muestra ruido 1/ ƒ , una característica común a muchos sistemas complejos de la naturaleza.
Este modelo solo 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 hacia 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 campo conforme asociada consiste en fermiones simplécticos con una carga central c = −2. [7]
La estabilización de las configuraciones de chips obedece a una forma de principio de mínima acción : cada vértice no se derrumba más de lo necesario en el curso de la estabilización. [8] Esto se puede formalizar de la siguiente manera. Llamemos a una secuencia de derrumbes legal si solo derrumba 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, derrumbando mientras sea posible. Dicha 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 derrumbamiento; es decir, para cualquier vértice , el número de veces que se derrumba es el mismo en todas las secuencias estabilizadoras legales. Según el principio de mínima acción, una secuencia estabilizadora mínima también es equivalente hasta la permutación del orden de derrumbamiento a una secuencia legal (y aún estabilizadora). En particular, la configuración resultante de una secuencia estabilizadora mínima es la misma que la 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 vuelca durante la estabilización (a través del vértice inestable) de una configuración de chip , y es un vector integral (no necesariamente no negativo) tal que es estable, entonces 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 , por lo que 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 "converger a una imagen continua". Matemáticamente, esto sugiere la existencia de límites de escala de la identidad de la pila 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, la existencia de límites de escala de configuraciones recurrentes de pilas de arena ha sido probada por Wesley Pegden y Charles Smart. [9] [10] En un trabajo conjunto posterior con Lionel Levine, utilizan el límite de escala para explicar la estructura fractal de la pila 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 los trabajos 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 para cuadrículas infinitas. Un desafío en tales generalizaciones es que, en general, ya no se garantiza que todas las avalanchas se detengan eventualmente. Por lo tanto, varias de las generalizaciones solo consideran la estabilización de configuraciones para las cuales esto se puede garantizar.
Un modelo bastante popular sobre 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 lleva una gran cantidad de granos de arena, la configuración después de la relajación forma patrones fractales (ver figura). Al permitir que la cantidad 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 arbitrarios dirigidos. Las reglas son que cualquier vértice con
es inestable; al derribarlo nuevamente se envían chips a cada uno de sus vecinos, uno a lo largo de cada borde saliente:
y, para cada uno :
¿Dónde está el número de aristas de a ?
En este caso, la matriz laplaciana no es simétrica. Si especificamos un sumidero tal que exista un camino desde cada otro vértice hasta , entonces la operación de estabilización en grafos 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 con raíz en el sumidero.
Para comprender mejor la estructura del grupo de pilas de arena para diferentes cuadrículas convexas finitas 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, los vértices en el interior de la cuadrícula todavía solo pueden llevar números enteros de granos. Las reglas de derribo permanecen sin cambios, es decir, se supone que tanto los vértices interiores como los de los límites se vuelven inestables y se derrumban si el número de granos alcanza o supera cuatro.
Además, las configuraciones recurrentes del modelo de pila de arena extendida forman un grupo abeliano, denominado grupo de pila de arena extendida , del cual el grupo de pila de arena usual es un subgrupo discreto . A diferencia del grupo de pila de arena usual, el grupo de pila de arena extendida es, sin embargo, un grupo de Lie continuo . Dado que se genera simplemente añadiendo granos de arena al límite de la cuadrícula, el grupo de pila de arena extendida tiene además la topología de un toro de dimensión y un volumen dados por el orden del grupo de pila de arena usual. [14]
De particular interés 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 por la identidad. Esta cuestión conduce a la definición de la dinámica de la pila de arena.
respectivamente
inducida 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 orden bajo, la dinámica de las pilas de arena se caracteriza por la transformación suave y la conservación aparente de los parches que constituyen la identidad de la pila 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 usuales, lo que significa un mapeo de configuraciones recurrentes en una cuadrícula dada a configuraciones recurrentes en una subcuadrícula. De manera informal, esta renormalización simplemente asigna configuraciones que aparecen en un momento dado en la dinámica de la pila 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 de la pila de arena inducida por la restricción de la subcuadrícula respectiva. [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 resulta en la situación de que, si un sitio está lleno en el momento , estará lleno para todos los momentos posteriores.
El montón de arena de Bak-Tang-Wiesenfeld fue mencionado en el episodio "Rampage" de Numb3rs , cuando 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 pila de arena abeliana sobre una cuadrícula hexagonal finita donde, en lugar de colocar los granos de forma aleatoria, los jugadores los colocan.