Criticalidad autoorganizada: carga lenta y descarga colectiva

COMPLEXA · Ciencias de la Complejidad y Teoría del Caos
Módulo 18. Organización fuera del equilibrio y redes de reacción

179.1. Criticalidad ajustada externamente y autoorganización

La criticalidad autoorganizada describe una dinámica que puede aproximarse a un régimen crítico mediante la interacción entre alimentación, transporte y disipación, en la que el recurso se acumula lentamente y se redistribuye en episodios de actividad colectiva. La pila de arena permite formular con precisión qué reglas producen esa organización, qué límites justifican hablar de criticalidad y qué observaciones distinguen este mecanismo de otras fuentes de fluctuaciones amplias.

En una transición crítica ajustada externamente, un parámetro de control se fija cerca de un valor especial. En un sistema magnético, por ejemplo, aproximarse a la temperatura crítica modifica las correlaciones y la respuesta. En un sistema autoorganizado, una variable interna ligada a la disponibilidad del recurso evoluciona por el impulso y las pérdidas; la proximidad al régimen crítico depende de esa retroalimentación dinámica y de condiciones físicas específicas.

Cuando la actividad es escasa, la entrada aumenta el recurso almacenado; si una perturbación desencadena redistribuciones extensas, el transporte hacia regiones disipativas lo reduce. Este acoplamiento puede mantener el sistema cerca de una región de elevada susceptibilidad, pero debe especificarse junto con el tamaño, las fronteras, la tasa de entrada y cómo cambian al considerar sistemas mayores.

Bak, Tang y Wiesenfeld (1987) propusieron la criticalidad autoorganizada como explicación de comportamientos sin escala característica dominante y del ruido de baja frecuencia. Su propuesta convirtió la intuición de acumulación y descarga en un programa de modelos verificables. La denominación abarca desarrollos posteriores diversos, por lo que las propiedades demostradas para una pila concreta deben distinguirse de las hipótesis sobre sistemas naturales.

El caos se caracteriza por la sensibilidad a condiciones iniciales en una evolución determinista; la criticalidad, por la organización colectiva y el escalamiento. Pueden coexistir, pero ninguna implica la otra. La complejidad en general tampoco identifica un mecanismo crítico; distribuciones heterogéneas o señales intermitentes permiten plantear preguntas, pero discriminan poco entre modelos dinámicos.

De modo semejante, un espectro próximo al ruido $1/f$, donde $f$ es la frecuencia, admite mecanismos alternativos y depende del observable. El análisis separa reglas microscópicas, observables de avalancha e inferencias estadísticas para tratar la pila como objeto matemático y evaluar con cautela su pertinencia empírica. La interpretación debe explicitar escalas temporales, balance del recurso y convenciones para identificar eventos.

179.2. El modelo de la pila de arena BTW en dos dimensiones

El modelo BTW bidimensional se define sobre una red cuadrada de lado $L$, con $L^2$ sitios. Cada sitio $i$ contiene una altura entera no negativa $h_i$, que cuenta los granos almacenados. La configuración es estable si todas las alturas están entre cero y tres. El umbral de inestabilidad es cuatro en el interior, los bordes y las esquinas.

Una adición incrementa en una unidad la altura de un sitio elegido. Si alcanza o supera cuatro, el sitio vuelca (pierde cuatro granos y envía uno a cada vecino existente). Los sitios inestables se relajan con la misma regla hasta recuperar estabilidad; la secuencia de vuelcos desencadenada es una avalancha. Una adición también puede no provocar ningún vuelco.

Aquí $h_j'$ es la altura posterior al vuelco de $i$; $j\sim i$ indica vecindad entre sitios existentes y $\Delta$ es la matriz de vuelco. La diagonal vale cuatro incluso cuando faltan vecinos. Con fronteras abiertas, los granos dirigidos al exterior abandonan el sistema; así, un vuelco interior conserva la masa, mientras que uno en la frontera puede reducirla.

La masa es la suma de las alturas. Un sitio de borde que no sea esquina tiene tres vecinos y pierde netamente un grano al volcar; una esquina tiene dos y pierde dos. Este drenaje permite estabilizar cualquier avalancha en una red finita conectada a la frontera. Con fronteras periódicas y sin otro mecanismo disipativo, esa garantía desaparece y el protocolo requiere otro análisis.

La metáfora granular no implica que la altura reproduzca una superficie real ni que el umbral represente directamente un ángulo de reposo. El modelo idealiza almacenamiento local, activación por umbral, transferencia a vecinos y salida del recurso para estudiar sus consecuencias colectivas.

La actualización puede elegir sucesivamente sitios inestables o proceder en generaciones paralelas. Ambas convenciones usan la misma regla local, pero asignan tiempos diferentes; la equivalencia del estado final requiere condiciones que se precisarán al discutir la propiedad abeliana. Conviene distinguir la configuración de almacenamiento de la historia de vuelcos y transporte, porque una misma altura final puede ocultar una redistribución extensa durante la relajación.

179.3. Impulso lento y separación de escalas

El protocolo de impulso lento añade un grano solo cuando la avalancha anterior ha terminado. Se elige un sitio aleatorio, se incrementa su altura y se suspende la alimentación durante toda la relajación. La dinámica alterna así configuraciones estables y episodios de actividad. Esta definición separa con precisión la perturbación externa de la respuesta interna; así, cada vuelco puede asignarse a una adición determinada y el catálogo de avalanchas resulta inequívoco.

Una configuración estable es absorbente para la dinámica de relajación, pues sin una nueva entrada permanece inmóvil indefinidamente. El sistema impulsado, en cambio, puede seguir recorriendo configuraciones porque las adiciones reactivan su evolución. Esta distinción evita confundir el reposo entre eventos con el final de toda la dinámica y permite describir el proceso como una sucesión de estados estables conectados por operadores de adición y estabilización.

Si se introduce un reloj físico, la separación exige que la probabilidad de recibir nuevas entradas durante una avalancha sea despreciable. Para una alimentación aleatoria con tasa total $r$ y una avalancha de duración $T$, el producto de ambas cantidades es el número esperado de entradas durante ese episodio. El criterio debe aplicarse al intervalo de duraciones relevante y, especialmente, a la escala temporal que crece con el tamaño.

En esta expresión, $T_c(L)$ es una duración característica de corte para una red de lado $L$. Si se especifica una tasa por sitio $\lambda$, la tasa total es $r=\lambda L^2$. Por tanto, una alimentación local aparentemente pequeña puede producir solapamiento al aumentar la red. La lentitud del impulso tiene un significado relativo, pues depende del reloj de actualización, del número de sitios y de las duraciones colectivas accesibles.

El protocolo condicionado al reposo impone la separación exactamente para cada tamaño finito. En cambio, mantener una tasa física constante y aumentar indefinidamente el sistema define otro problema. El resultado depende del orden de los límites, porque las avalanchas más largas pueden crecer y encontrar nuevas entradas antes de extinguirse. Comparar simulaciones exige indicar si se usa alimentación suspendida, intervalos deterministas entre adiciones o una fuente aleatoria que sigue actuando durante la actividad.

Tamaño, disipación e impulso son controles distintos, aunque sus efectos se acoplen. El tamaño limita las distancias de propagación; la disipación regula cuánto recurso desaparece; el impulso determina cuándo se introduce recurso nuevo. Durante una avalancha lenta, la configuración se descarga y reorganiza sin aportes adicionales; entre avalanchas, la acumulación prepara otra respuesta. La descripción crítica requiere analizar conjuntamente estos controles y declarar cuáles permanecen fijos y cuáles tienden a cero o a infinito.

179.4. Un ejemplo exacto en una cadena unidimensional

Considérese, como ejemplo didáctico, una cadena de tres sitios con umbral dos y extremos abiertos. Cada vuelco retira dos granos del sitio activo y entrega uno a cada vecino existente; el grano dirigido fuera de la cadena se pierde. Las alturas estables son cero y uno. Desde la configuración inicial estable (1,1,1), añadir un grano al centro produce

El centro es el único sitio inestable. Su vuelco da (2,0,2); ambos extremos alcanzan el umbral y vuelcan en la generación siguiente. Cada uno envía un grano al centro y pierde otro por su frontera abierta. El centro vuelve a ser inestable y su segundo vuelco produce una configuración estable, según la secuencia

La historia contiene cuatro vuelcos (dos en el centro y uno en cada extremo). El área es de tres sitios, pues todos vuelcan al menos una vez; se pierden dos granos. El tamaño cuenta operaciones y el área, sitios distintos; ambos difieren cuando un sitio se activa repetidamente.

Con $M_{\mathrm{ini}}$ y $M_{\mathrm{fin}}$ como masas inicial y final, $I$ como entrada y $\ell$ como pérdida, el balance es

El tiempo depende de la convención. La actualización paralela dura tres generaciones activas (centro, extremos y centro); la secuencial, cuatro pasos si cada paso corresponde a un vuelco. Los extremos pueden procesarse en cualquier orden legal. En este modelo abeliano, el orden cambia las configuraciones intermedias, mientras que la estabilización final y el número de vuelcos de cada sitio permanecen iguales.

Estas cifras son resultados exactos de un ejemplo didáctico; no son observaciones empíricas ni evidencia de exponentes o de propiedades asintóticas de una pila bidimensional. El ejemplo fija el vocabulario operativo y muestra que la duración depende del reloj. Informar «duración cuatro» sin unidad temporal, regla de actualización y criterio de inicio y final deja incompleto el registro.

179.5. Tamaño, área, duración y masa perdida

Una avalancha admite varias medidas que describen aspectos diferentes de su historia. El tamaño $s$ cuenta todos los vuelcos (incluidas las reactivaciones de un mismo sitio) y el área $a$, los sitios distintos que vuelcan. La duración $T$ mide la extensión temporal según un reloj explícito. La pérdida $\ell$ cuenta los granos que abandonan el sistema. Ninguna de estas variables puede sustituirse automáticamente por otra al interpretar un catálogo de eventos.

Sea $n_i$ el número de vuelcos del sitio $i$ durante una avalancha. El tamaño se obtiene sumando esos números y el área sumando indicadores de actividad. Si $b_i$ es el número de vecinos ausentes del sitio en la red cuadrada abierta, cada vuelco suyo pierde $b_i$ granos. Esta descripción expresa tres observables en función de la historia microscópica sin introducir una convención temporal particular.

El símbolo $\mathbf{1}_{\{n_i>0\}}$ vale uno cuando el sitio ha volcado y cero en caso contrario. De las definiciones se deduce que el área nunca supera al tamaño ni al número de sitios de la red. La pérdida, en cambio, puede ser cero en una avalancha que redistribuye mucho recurso sin salir del interior y positiva en otra de menor tamaño que alcanza enseguida la frontera.

Con actualización paralela, la duración cuenta las generaciones con al menos un sitio activo; el tamaño agrega todos los vuelcos de esas generaciones. Con un reloj secuencial que asigna una unidad a cada vuelco, duración y tamaño coinciden por construcción; esa coincidencia no es una propiedad independiente de la dinámica. Los tiempos de espera entre adiciones son otra variable y deben distinguirse del tiempo interno de relajación de una avalancha.

Para cada observable pueden estimarse distribuciones marginales, así como distribuciones conjuntas que describan sus dependencias. Condicionar a una pérdida positiva, por ejemplo, selecciona los eventos que alcanzaron una región disipativa y modifica así el conjunto estudiado. También hay que decidir si se incluyen las adiciones que no provocan vuelcos, cuyo tamaño y área son cero, o si se trabaja solo con el catálogo condicionado a la aparición de actividad.

En sistemas finitos aparecen escalas de corte que limitan el rango de comportamiento aproximadamente invariante de escala. El corte del área está restringido por la geometría; los cortes de tamaño y duración incorporan además recurrencia local y propagación temporal. Una cola amplia observada en un solo observable da información parcial. La comparación entre tamaños de red y entre medidas compatibles permite investigar si los cortes responden a una organización colectiva o a restricciones de observación y definición.

179.6. Hipótesis de escalamiento y exponentes

Una hipótesis habitual describe las avalanchas mediante una ley de potencias aproximada en un régimen intermedio, con un corte que depende del tamaño de la red y suprime los eventos mayores. Su utilidad reside en las restricciones conjuntas que impone a distribuciones medidas para distintos valores de $L$.

Aquí $P(s;L)$ es la probabilidad del tamaño $s$ en una red de lado $L$; $\tau$ caracteriza el régimen algebraico, $D$ controla el crecimiento del corte y $f$ decae para argumentos grandes. El ansatz describe un régimen asintótico, con escala microscópica inferior y posibles correcciones. Deben especificarse la normalización y el tratamiento de eventos nulos antes de ajustar datos.

Los momentos son cantidades integradas que permiten examinar la hipótesis. Para un orden real $q$, se pondera la distribución con $s^q$. Si domina el extremo superior y converge la integral de la función de corte, cambiar la variable al cociente entre tamaño y corte produce una dependencia algebraica con $L$; esta condición es esencial para interpretar la relación.

La función $\sigma(q)$ expresa el escalamiento de cada momento. La segunda igualdad corresponde al ansatz simple bajo esas condiciones; en el límite pueden aparecer términos logarítmicos y, para otros órdenes, pueden dominar los tamaños pequeños. Una dependencia no lineal de los exponentes de los momentos con el orden puede indicar multiescalamiento, aunque tamaños insuficientes, transitorios o correcciones también pueden explicar desviaciones aparentes.

De Menech, Stella y Tebaldi (1998) documentaron la ruptura del escalamiento simple de los tamaños de avalancha en el modelo abeliano bidimensional, con un papel destacado de los eventos raros que alcanzan la frontera. Según este resultado, describir todas esas estadísticas mediante un único par de exponentes puede fallar, aunque algunas relaciones de momentos y otros observables admitan leyes más sencillas. Una recta en coordenadas logarítmicas tampoco demuestra el ansatz; para contrastarlo se examinan el movimiento del corte, los momentos y la estabilidad de las estimaciones al variar el intervalo. Aquí los exponentes se dejan en forma simbólica, pues asignarles valores requiere precisar observable, modelo y procedimiento de estimación.

Las pilas estocásticas modifican las reglas de redistribución y pueden pertenecer a clases de universalidad distintas. La comparación requiere identificar simetrías, conservación, dimensionalidad y estructura del ruido, ya que colas semejantes no prueban identidad de clase. Las relaciones de escalamiento se evalúan primero dentro de una formulación concreta; compartir umbrales, carga y descarga no justifica transferir exponentes entre dinámicas.

179.7. Conservación, estacionariedad y propiedad abeliana

En el BTW abierto, la conservación es local en el interior, donde cada vuelco redistribuye cuatro granos sin cambiar la masa total; las pérdidas se concentran en la frontera, donde parte del transporte sale de la red. Si se añade disipación en el interior, cambia la dinámica y debe incorporarse explícitamente una contribución adicional a las pérdidas.

Sea $M_n$ la masa de la configuración estable anterior a la adición $n$, $I_n$ la entrada y $\ell_n$ la pérdida durante su relajación. El balance vale para cada evento, con independencia de que la distribución ya sea estacionaria. Al promediar en un régimen estacionario, la masa media deja de variar y la entrada media debe igualar la pérdida media. El promedio incluye las adiciones que no desencadenan actividad.

Para una entrada de un grano por adición, la pérdida media por adición es uno. Esto permite que muchas avalanchas no pierdan masa y que otras compensen mediante descargas mayores. El balance es una restricción exacta, pero deja sin determinar la distribución completa de pérdidas y la forma de las colas. Tampoco implica independencia estadística entre eventos; de hecho, la configuración heredada de una avalancha condiciona la respuesta a las adiciones posteriores.

En una red finita con acceso disipativo desde todos los sitios, cada avalancha termina. La terminación puede justificarse mediante un potencial positivo construido con la matriz de vuelco, pues existe una ponderación espacial para la cual cada vuelco disminuye el potencial total. Como las alturas permanecen no negativas, la relajación no admite una sucesión infinita de vuelcos. Esta finitud depende del drenaje; en una red cerrada que conserve toda la masa debe establecerse por separado.

La propiedad abeliana estándar afirma que, para una configuración y unas adiciones dadas, la estabilización final y el número de vuelcos de cada sitio son independientes del orden legal de actualización. Un orden es legal si solo vuelcan sitios que han alcanzado el umbral. La conmutatividad algebraica de las operaciones locales, junto con la terminación y las condiciones de estabilización, permite establecer esa independencia. Las configuraciones intermedias y el número de generaciones paralelas requieren información temporal adicional.

El balance local permite obtener un resultado más informativo que la igualdad global entre entrada y salida. Volvamos a la cadena abierta con umbral dos, ahora de $L$ sitios y con adiciones uniformes de un grano. Sea $u_i$ el número medio de vuelcos del sitio $i$ por adición, incluidas las adiciones que no provocan actividad. En un régimen estacionario, el almacenamiento medio de cada sitio deja de cambiar. Por tanto,

La entrada media local es $1/L$; los dos términos vecinos son el transporte recibido. La solución es una parábola discreta, de la que se obtiene el tamaño medio

Con tres sitios se obtiene $u=(1/2,2/3,1/2)$ y tamaño medio $5/3$. Este promedio no tiene que coincidir con el tamaño cuatro de la avalancha particular de 179.4. La pérdida media es $u_1+u_L=1$, de acuerdo con el balance global. La solución requiere estacionariedad y el protocolo de adición especificado, sin suponer independencia entre avalanchas ni una forma para su distribución.

El crecimiento cuadrático del tamaño medio en esta cadena tiene una explicación de transporte, ya que los granos añadidos deben recorrer el sistema antes de salir. Probar esa ley para el primer momento no demuestra que todos los momentos sigan una sola función de escalamiento. Si además se propone un ansatz simple con primer momento dominado por el corte y amplitudes de normalización compatibles, habría que satisfacer $D(2-\tau)=2$. La ecuación restringe exponentes candidatos; no los determina por separado ni valida el ansatz del que procede. Condicionar a avalanchas activas introduce además su probabilidad, que puede depender de $L$.

En una red general, el mismo razonamiento adopta la forma matricial $\Delta^{\mathsf T}u=p$, donde $p_i$ es la probabilidad de añadir en $i$ y la convención de $\Delta$ es la de 179.2. La inversa existe para el grafo finito disipativo considerado. Así, cambiar la localización de las entradas modifica los promedios aun sin cambiar la regla de vuelco. El enfoque algebraico de Dhar (1990) permite relacionar estas propiedades exactas con la estructura de las configuraciones recurrentes del modelo abeliano.

La distribución estacionaria expresa invariancia estadística bajo adición y relajación; no implica equilibrio termodinámico ni balance detallado. El sistema sigue alimentado y disipativo, con transporte dirigido desde la fuente hacia los sumideros. La abelianidad no garantiza duraciones iguales bajo relojes distintos. Conservación, estacionariedad y conmutatividad son propiedades distintas; separarlas permite formular resultados exactos sin atribuirles consecuencias dinámicas o estadísticas que no se derivan de ellas.

179.8. Impulso finito, disipación y cuasicriticalidad

Cuando la alimentación continúa durante la actividad, dos perturbaciones pueden contribuir a una misma región activa. La identificación de avalanchas deja de estar determinada solo por el retorno al reposo, pues un episodio puede reunir cascadas iniciadas por entradas distintas; sus estadísticas dependen entonces de la segmentación. La tasa finita introduce una escala temporal externa que compite con la propagación interna, modifica los cortes y sostiene actividad durante intervalos prolongados.

La disipación interior introduce otra escala. Si cada transferencia puede perder recurso antes de llegar a la frontera, la perturbación se atenúa durante su propagación. Mantener disipación interior no nula al aumentar el sistema puede impedir que las correlaciones crezcan indefinidamente. El resultado depende de las reglas concretas; debe distinguirse esta pérdida distribuida de la disipación fronteriza del BTW estándar y declararse cómo cambia cada mecanismo con el tamaño.

El recurso disponible puede ser finito y recuperarse con una escala temporal propia. Tras una descarga extensa, las nuevas activaciones dependen del agotamiento y de la velocidad de recuperación. Estas dinámicas admiten oscilaciones, excursiones alrededor de una transición y actividad intermitente. La retroalimentación estabilizadora restringe la evolución, pero que la distancia a un punto crítico desaparezca en el límite apropiado requiere una demostración aparte.

Un sistema es cuasicrítico si permanece cerca de una región crítica sin concentrarse exactamente en ella. Las fluctuaciones de la variable de control pueden conservar una anchura apreciable y mezclar respuestas de distinta susceptibilidad. Un intervalo amplio de comportamiento aproximadamente algebraico puede coexistir con cortes persistentes. Para distinguir la cuasicriticalidad de la criticalidad asintótica, se analizan las distribuciones de avalanchas y la evolución y fluctuaciones del recurso que regula la actividad.

El orden de los límites termodinámicos es decisivo. Llevar primero el impulso a cero para cada red finita preserva la separación entre eventos; aumentar primero la red con alimentación puede destruirla. Del mismo modo, hacer desaparecer la disipación interior requiere especificar su relación con el crecimiento espacial. Decir «sistema grande e impulso pequeño» deja indeterminada la trayectoria en el espacio de parámetros; trayectorias distintas pueden conducir a regímenes estadísticos e inferencias incompatibles.

Bonachela y Muñoz (2009) examinan modelos en los que la autoorganización sin conservación produce una aproximación fluctuante a la región crítica y puede requerir ajuste adicional para obtener criticalidad asintótica. Su resultado identifica condiciones y fallos concretos de una explicación demasiado general, sin prohibir toda forma de criticalidad no conservativa. Ese trabajo obliga a precisar qué ajuste permanece implícito y cuál es la amplitud de las fluctuaciones alrededor de la transición. Evaluar la SOC exige especificar las condiciones del modelo y analizar los límites; en aplicaciones, fuente, disipación y recuperación deben medirse o justificarse independientemente, pues una cola amplia es compatible con varios mecanismos de regulación.

179.9. Contraste estadístico y evidencia del mecanismo

La evidencia de criticalidad autoorganizada combina estadísticas de eventos con información sobre el mecanismo que los produce. Una cola amplia es solo un punto de partida, cuya interpretación depende del rango, la definición de evento y las alternativas consideradas. En datos naturales, el tamaño efectivo del sistema y la intensidad del impulso pueden estar parcialmente ocultos, lo que limita la interpretación física de un ajuste.

Un contraste puede comparar una ley de potencias, una distribución lognormal y una exponencial truncada con las mismas observaciones y un soporte compatible. La razón de verosimilitud resume la preferencia relativa entre modelos ajustados y requiere estimar su incertidumbre y controlar la selección. Que favorezca una alternativa no prueba la bondad absoluta del ajuste; identificar el proceso dinámico exige además información sobre el mecanismo.

Aquí $x_k$ es la observación $k$ y $N$ el número de eventos; $p$ y $q$ son modelos con parámetros estimados $\widehat{\theta}$ y $\widehat{\phi}$. Una suma positiva favorece relativamente a $p$. La dependencia temporal puede invalidar errores estándar que suponen eventos independientes; la estimación de la incertidumbre debe conservar o modelar las correlaciones y la estructura del muestreo.

Comparar varios tamaños $L$ permite estudiar el desplazamiento de los cortes y las relaciones de escalamiento, aunque cambiar el campo observado puede significar solo muestrear una fracción del mismo sistema. El submuestreo altera la actividad registrada y fragmenta episodios conectados; Levina y Priesemann (2017) muestran cómo afecta a la interpretación de distribuciones parciales. El binning temporal también influye, porque los intervalos pequeños fragmentan cascadas y los grandes fusionan respuestas distintas. Umbrales de detección y límites instrumentales sesgan los extremos, por lo que un análisis reproducible documenta las decisiones, prueba la estabilidad ante elecciones razonables y conserva información sobre eventos descartados y cambios en la capacidad de observación.

En neuronas y terremotos, el vocabulario de avalanchas debe apoyarse en datos sobre activación, transferencia, recuperación y disipación; la semejanza con una pila discrimina poco entre explicaciones causales. Son más informativas las predicciones conjuntas sobre cortes, correlaciones y respuesta a cambios del impulso. Si no pueden variarse el tamaño ni intervenirse el sistema, las conclusiones quedan limitadas por los datos disponibles y por la capacidad explicativa de modelos alternativos.

Ajuste de colas y controles de muestreo

Para tamaños discretos, una ley de potencias por encima de un umbral inferior $s_{\min}$ requiere una normalización discreta. Si no se introduce un corte superior, el modelo candidato es

Su logaritmo de verosimilitud, para $m$ observaciones incluidas en la cola, es $-\tau\sum_{k=1}^m\ln s_k-m\ln\zeta(\tau,s_{\min})$. Omitir la constante de normalización cambia la estimación. En una pila finita, este candidato sin corte solo puede describir un intervalo; si los datos alcanzan la escala máxima, hay que incorporar su supresión o restringir explícitamente el alcance del ajuste. La selección del umbral también añade incertidumbre y debe repetirse dentro del procedimiento de evaluación cuando se estime con los mismos datos.

Los histogramas con intervalos de anchura variable requieren dividir las frecuencias relativas por esa anchura para representar una densidad. En una cola ideal continua, una densidad proporcional a $s^{-\tau}$ tiene supervivencia proporcional a $s^{1-\tau}$; confundir ambas pendientes desplaza el exponente en una unidad. Esta relación necesita $\tau>1$ y un intervalo alejado del corte. En datos discretos se usa la suma acumulada correspondiente.

Las series deben conservar su orden temporal. Un remuestreo por bloques permite retener parte de la dependencia entre avalanchas, siempre que la longitud de bloque sea pertinente para las correlaciones observadas. Conviene contrastar varias longitudes, separar el transitorio de preparación y comprobar el balance de masa en el registro completo. El ajuste de una cola, su incertidumbre y el balance son controles complementarios, porque cada uno detecta fallos que los otros pueden dejar pasar.

179.10. Límites de la predicción

Las colas amplias no predicen el próximo gran evento. La evidencia puede imponer restricciones sobre el recurso, la recuperación y la disipación; la localización exacta de la siguiente descarga queda fuera de su alcance. La criticalidad tampoco implica una función óptima; esa valoración requiere definir objetivos, restricciones y alternativas comparables.

La evaluación de eventos raros exige estadística reproducible, criterios de detección explícitos y tratamiento de la incertidumbre. Hay que contrastar los mecanismos de alimentación y pérdida, comparar tamaños y examinar el submuestreo. El ejemplo BTW distingue resultados exactos (regla de vuelco, balance de masa e independencia del orden legal) de hipótesis de escalamiento; la duración depende además de la convención temporal. Establecer universalidad o capacidad predictiva en otros sistemas exige argumentos que van más allá de estas propiedades.

Referencias

Bak, P., Tang, C. y Wiesenfeld, K. (1987). Self-organized criticality: An explanation of the 1/f noise. Physical Review Letters, 59, 381–384. DOI.

Bonachela, J. A. y Muñoz, M. A. (2009). Self-organization without conservation: true or just apparent scale-invariance? Journal of Statistical Mechanics: Theory and Experiment, P09009. DOI. Texto en arXiv.

Levina, A. y Priesemann, V. (2017). Subsampling scaling. Nature Communications, 8, 15140. DOI.

De Menech, M., Stella, A. L. y Tebaldi, C. (1998). Rare events and breakdown of simple scaling in the Abelian sandpile model. Physical Review E, 58, R2677–R2680. DOI.

Dhar, D. (1990). Self-organized critical state of sandpile automaton models. Physical Review Letters, 64, 1613–1616. DOI.