LBan-V
Avalanches at the Depinning Transition
In the previous lesson, we studied the dynamics of an interface in a disordered medium under a uniform external force . In a fully connected model, we derived the force–velocity characteristic and identified the critical depinning force 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} .
In this lesson, we focus on the avalanches that occur precisely at the depinning transition. To do so, we introduce a new driving protocol: instead of controlling the external force 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} , we control the position of the interface by coupling it to a parabolic potential. Each block is attracted toward a prescribed position through a spring of stiffness 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 k_0} .
For simplicity, we restrict to the fully connected model, where the local force acting on block 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 i} is
Here 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 h_{CM}} is the center-of-mass position of the interface. Because the sum of all internal elastic forces vanishes, the total external force is balanced solely by the pinning forces. The effective external force can thus be written as a function 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} :
As 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} is increased quasistatically, the force 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} would increase if were fixed. When an avalanche takes place, 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 h_{CM}} jumps forward and 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(w)} suddenly decreases. However, in the steady state and in the thermodynamic limit , the force recovers a well-defined value. 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 k_0 \to 0} , this force tends to the critical depinning force 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} ; at finite it lies slightly below 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} .
As in the previous lesson, it is convenient to introduce the variables 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} that measure the distance of block from its local instability threshold:
These variables are the natural starting point for describing avalanches and their statistics.
Quasi-Static Protocol and Avalanche Definition
To study avalanches, the position 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} is increased quasi-statically: it is shifted by an infinitesimal amount so that the block closest to its instability threshold reaches it, i.e.
This block is the epicenter of the avalanche: it becomes unstable and jumps to the next well. Its stabilization induces a redistribution of the stress over all other blocks, which may in turn become unstable:
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 \begin{cases} x_i = 0 \;\longrightarrow\; x_i = \Delta (1 + k_0), \\[6pt] x_j \;\longrightarrow\; x_j - \dfrac{\Delta}{L} \quad \text{for } j \neq i. \end{cases} }
The key feature of the quasi-static protocol is that does not evolve during the avalanche: all subsequent destabilizations are triggered exclusively by previously unstable blocks.
It is convenient to organize the avalanche into generations of unstable sites:
- The first generation consists of the single epicenter block.
- The second generation is formed by the blocks destabilized by the stabilization of the epicenter.
- The third generation consists of the blocks destabilized by the stabilization of the second-generation blocks, and so on.
This hierarchical picture allows us to characterize the size and temporal structure of avalanches.
Avalanche Size
Once the avalanche is over, its size 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 S} is a random variable defined as:
where is the total number of instabilities (itself a random variable) and 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 \Delta_i} is the jump of block 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 i} .
Approximating the sum as makes explicit the proportionality between the number of instabilities and the total avalanche size:
This relationship highlights that larger avalanches correspond to a larger number of destabilized blocks.
Derivation of the Evolution Equation
Our goal is to determine the 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_w(x)} of the distances to threshold of all blocks, given their initial distribution 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} . We perform an infinitesimal change in the position of the parabolic potential 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 + \mathrm{d}w} .
Recall the expression of the distance to threshold of block just before the change:
After the change, a complex reorganization of the 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} variables occurs. We organize this reorganization generation by generation, indexing them by a discrete time :
- At time 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 t = 1} (first generation):
The center of mass is still 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 h_{CM}(w)} and two things can happen:
1. **Stable blocks:** if , the block approaches its threshold:
2. **Unstable blocks:** if 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 0 < x_i(w) < k_0 \, \mathrm{d}w}
, the block becomes unstable and is stabilized.
Since is infinitesimal and 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(w) \approx 0}
, we simply have:
This event occurs with probability 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(0) \, k_0 \, \mathrm{d}w} .
- At time 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 t = 2} (second generation):
The position of the parabolic potential remains fixed 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 w + \mathrm{d}w} , but the center of mass of the interface advances by
Again, two cases are possible:
1. **Stable blocks:** if 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^{t=0}(w+\mathrm{d}w) > k_0 \, \mathrm{d}w} , the block approaches its threshold:
2. **Unstable blocks:** if 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 0 < x_i^{t=0}(w+\mathrm{d}w) < k_0 \, \mathrm{d}w} , the block becomes unstable and is stabilized. As before, since 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 \mathrm{d}w} is infinitesimal and 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^{t=0}(w+\mathrm{d}w) \approx 0} :
This event occurs with probability 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 \overline{\Delta} \, P_w(0) \, k_0 \, \mathrm{d}w} .
This procedure can be iterated to higher generations 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 t = 3, 4, \dots} until the avalanche stops.
Dynamics
Our goal is thus to determine the 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_w(x)} of all blocks, 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 } . Let's decompose in steps the dynamics
- 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 block decreases its distance to threshold
.
As a consequence
- Stabilization : 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. The stabilization induces the change 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) \to m^2 d w P_w(0) g(x) } . Hence, one writes
The stabilization of the unstable blocks induce a drop of the force per unit length
\
- Redistribution This drop is (partially) compensated by the redistribution. The force acting on all points is increased:
Again, most of the distribution will be driven to instability while a fraction of the blocks become unstable... we can write
and finally:
Stationary solution
Increasing the drive the distribution converge to the fixed point:
- Determine 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?
We consider an avalanche starting from a single unstable site 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_0=0 } and the sequence of sites more close to instabitity 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_1< x_2<x_3\ldots } . For each unstable block, all the blocks receive a random kick:
with 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 \Delta_1,\Delta_2,\Delta_3, \ldots } drwan from 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 g(\Delta) } Are these kick able to destabilize other blocks?
Given the initial condition and 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 }
, 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) }
. From the extreme values theory we know the equation setting the average position of the most unstable block is
Hence, for large systems we have
Hence we need to compare the mean value of the kick with the mean gap between nearest unstable sites:
Note that 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 } simplifies. We expect three possibilities:
- if the mean kick is smaller than the mean gap the system is subcritical and avalanches quickly stops.
- if the mean kick is equal to the mean gap the system is critical and avalanches are power law distributed
- if the mean kick is larger of the mean gap the system is super-critical and avalanches are unstable.
Note that in the stationary regime the ratio between mean kick and mean gap is 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/(1+m^2) } . Hence, 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 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 }
Mapping to the Brownian motion
Let's define the random jumps and the associated random walk
An avalanche is active until 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_n } is positive. Hence, the size of the avalanche identifies with first passage time of the random walk.
- Critical case : In this case the jump distribution is symmetric and we can set 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_0=0} . Under these hypothesis the Sparre-Andersen theorem state that the probability that the random walk remains positive for 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 n} steps is independent on the jump disribution and for a large number of steps becomes 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 Q(n) \sim \frac{1}{\sqrt{\pi n}}} . Hence, the distribution avalanche size is
This power law is of Gutenberg–Richter type. The universal exponent is 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 \tau=3/2}
- Stationary regime: Replacing 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 \frac{1}{LP_w(0)}} with 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 \frac{1}{LP_{\text{stat}}(0)} = \frac{\overline{\Delta}}{L} } we get . For small m, the random walk is only sliglty tilted. The avalanche distribution will be power law distributed with 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 \tau=3/2} until a cut-off