L-6: Difference between revisions
Line 88: | Line 88: | ||
Let's define the random jumps and the associated random walk | Let's define the random jumps and the associated random walk | ||
<center><math> \eta_1 = \frac{\Delta_1}{(1+m^2)L}- x_1, \; \eta_2=\frac{\Delta_2}{(1+m^2)L}- (x_2-x_1) \; \eta_3=\frac{\Delta_3}{(1+m^2)L}- (x_3-x_2) \ldots | <center><math> \eta_1 = \frac{\Delta_1}{(1+m^2)L}- x_1, \; \eta_2=\frac{\Delta_2}{(1+m^2)L}- (x_2-x_1), \; \eta_3=\frac{\Delta_3}{(1+m^2)L}- (x_3-x_2) \ldots | ||
</math></center> | |||
<center><math> | |||
X_t= \sum_{i=1}^t \eta_i \; \test{with} \; \langle \eta_i\rangle = \left(\frac{\overline{\Delta}}{(1+m^2)} -\frac{1}{P_w(0)}\right)/L | |||
</math></center> | </math></center> |
Revision as of 19:59, 2 March 2024
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
.
- The local threshold are all equal. In particular we set
.
As a consequence, in the limit 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 L\to \infty}
, the statistical properties of the system are described by the distribution of the local stresses . For simplicity, instead of the stresses, we study the distance from threshold
Our goal is thus to determine their distribution , given their intial distribution, 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 P_0(x)} , and a value of 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 } .
Dynamics
Let's rewrite the dynamics with the new variables
- Drive: 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 + dw} each point decreases its distance to threshold
.
As a consequence
- Instability 1: Stress drop The instability occurs when a point is at 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 x_i =0 } . Then, the point is stabilized (stress drop):
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 + dw} , a fraction 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 m^2 d w P_w(0) } of the blocks is unstable. Due to the stress drop, their distance to threshold becomes . Hence, one writes
- Instability 2: Stress redistribution The stress drop of a single block induces a stress redistribution where all blocks approach threshold.
The total stress drop is hence all points move to the origin of
part of them shifts, part of them become unstable... we can write
and finally:
Stationary solution
Increasing the drive the distribution converge to the fixed point:
- Determne 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 P_{\text{stat}}(0) =\frac{1}{\overline{\Delta}} } using
- Show
which is well normalized.
Critical Force
The average distance from the threshold gives a simple relation for the critical force, namely 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 1-f_c= \overline{x} } . Hence for the automata model we obtain:
Exercise:
Let's assume an exponential distribution of the thresholds and show
- 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 P_{\text{stat}}(x)= e^{-x/\overline{\Delta}}/\overline{\Delta} }
- 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 f_c= 1- \overline{\Delta}}
Avalanches or instability?
Given the initial condition and , the state of the system is described by 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 P_w(x) } . For each unstable block, all the blocks receive a kick. The mean value of the kick is
Is this kick able to destabilize another block? The equation setting the average position of the most unstable block is
Hence, for large systems we have
We expect three possibilities:
- if the mean kick, 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 \sim \overline{\Delta}/(1+m^2) } is smaller than the mean gap , the system is subcritical and avalanches quickly stops.
- if the mean kick, 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 \sim \overline{\Delta}/(1+m^2) } is equal to the mean gap 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 \sim 1 /P_w(0)} , the system is critical and avalanches are power law distributed
- if the mean kick, is larger of the mean gap 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 \sim 1 /P_w(0)} , the system is super-critical and avalanches are unstable.
Noe that in the stationary regime the system is subcritical when 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 m>0 } and critical for
Mapping to the Brownian motion
Let's define the random jumps and the associated random walk