L-6

Avalanches and Bienaymé-Galton-Watson process

Goal: We solve the mean field version of the cellular automaton, derive its avalanche statistics and make a connection with the Bienaymé-Galton-Watson process used to describe an epidemic outbreak.

Fully connected (mean field) model for the cellular automaton

Let's study the mean field version of the cellular automata introduced in the previous lecture. We introduce two approximations:

• Replace the Laplacian, which is short range, with a mean field fully connected interction

$\sigma _{i}=h_{CM}-h_{i}+m^{2}(w-h_{i}),\quad$

.

• The local threshold are all equal. In particular we set

$\sigma _{i}^{th}=1,\quad \forall i$

.

As a consequence, in the limit $L\to \infty$, the statistical properties of the system are described by the distribution of the local stresses $\sigma _{i}$. For simplicity, instead of the stresses, we study the distance from threshold

$x_{i}=1-\sigma _{i}$

Our goal is thus to determine their distribution $P_{w}(x)$, given their intial distribution, $P_{0}(x)$, and a value of $w$.

Dynamics

Let's rewrite the dynamics with the new variables

• Drive: Increasing Failed to parse (unknown function "\dd"): {\displaystyle w \to w + \dd w} each point decreases its distance to threshold

Failed to parse (unknown function "\dd"): {\displaystyle x_i \to x_i - m^2 \dd w }

.

As a consequence

$P_{w+dw}(x)=P_{w}(x+dw)\sim P_{w}(x)+m^{2}dw\partial _{x}P_{w}(x)$

• Instability 1: Stress drop The instability occurs when a point is at $x_{i}=0$. Then, the point is stabilized (stress drop):

$x_{i}=0\to x_{i}=\Delta$

Increasing Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle w \to w + \dd w} , a fraction $m^{2}dwP_{w}(0)$ of the points of the interface is unstable. Due to the stress drop, their distance to instability will be $m^{2}dwP_{w}(0)g(x)$. Hence, one writes

$\partial _{w}P_{w}(x)\sim m^{2}\left[\partial _{x}P_{w}(x)+P_{w}(0)g(x)\right]$ ${\begin{cases}\sigma _{i}=\sigma _{i}-\Delta \quad {\text{stress drop}}\\\\\sigma _{i\pm 1}=\sigma _{i\pm 1}+{\frac {1}{2}}{\frac {\Delta }{1+m^{2}}}\quad {\text{stress redistribution}}\\\end{cases}}$

Note that $\Delta$ is a positive random variable drwan from $g(\Delta )$.

Hence, an unstable point, $x_{i}<0$, is stabilized to a value $x>0$ drawn from $g(x)$. The stress redistribution induced on each bloch is ${\frac {1}{L}}{\frac {x}{1+m^{2}}}$

we define  and write its evolution equation. 

• Drive: Changing $w\to w+dw$ gives
• Instability: This shift is stable far from the origin, however for a fraction $m^{2}dwP_{w}(0)$ of the points of the interface is unstable. Due to the stress drop, their distance to instability will be $m^{2}dwP_{w}(0)g(x)$. Hence, one writes

$\partial _{w}P_{w}(x)\sim m^{2}\left[\partial _{x}P_{w}(x)+P_{w}(0)g(x)\right]$

• Stress redistribution: as a consequence all points move to the origin of

$m^{2}dwP_{w}(0){\frac {\overline {x}}{1+m^{2}}}$

where ${\overline {x}}=\int dxxg(x)$

• Avalanche: Let us call we can write

$\partial _{w}P_{w}(x)=m^{2}\left[\partial _{x}P_{w}(x)+P_{w}(0)g(x)\right]\left[1+P_{w}(0){\frac {\overline {x}}{1+m^{2}}}+(P_{w}(0){\frac {\overline {x}}{1+m^{2}}})^{2}+\ldots \right]$

and finally:

$\partial _{w}P_{w}(x)={\frac {m^{2}}{1-P_{w}(0){\frac {\overline {x}}{1+m^{2}}}}}\left[\partial _{x}P_{w}(x)+P_{w}(0)g(x)\right]$

Stationary solution

Increasing the drive the distribution converge to the fixed point:

$0=\partial _{x}P_{\text{stat}}(x)+P_{\text{stat}}(0)g(x)$

• Determne $P_{\text{stat}}(0)={\frac {1}{\overline {x}}}$ using

$1=\int _{0}^{\infty }dx\,P_{\text{stat}}(x)=-\int _{0}^{\infty }dx\,x\partial _{x}P_{\text{stat}}(x)$

• Show

$P_{\text{stat}}(x)={\frac {1}{\overline {x}}}\int _{x}^{\infty }g(z)dz$