Avalanche dynamics, surface roughening, and self-organized criticality: Experiments on a three-dimensional pile of rice

Retinal neovascularization as self-organized criticality on ultra-widefield fluorescein angiography imaging of diabetic retinopathy

BACKGROUND/OBJECTIVES

Do the distributions of surface area of non-perfusion (NP) and neovascularization (NV) on ultra-widefield fluorescein angiography (UWF FA) in patients with diabetic retinopathy (DR) differ significantly?

SUBJECTS/METHODS

Inclusion criteria were patients who had a UWF FA taken for DR at the Kellogg Eye Center from January 2009 to May 2018. Exclusion criteria included previous panretinal photocoagulation and significant media opacity (e.g., vitreous haemorrhage or significant cataract). UWF FAs were manually segmented for surface areas of NP and NV. The total areas per patient were organized in a histogram, and logarithmically binned to test against power law and exponential distributions. Then, a computational model was constructed in Python 3.7 to suggest a mechanistic explanation for the observed distributions.

RESULTS

Analysis of areas of NV across 189 images demonstrated a superior fit by the least squares method to a power law distribution (p = 0.014) with an R2 fit of 0.9672. Areas of NP over 794 images demonstrated a superior fit with an exponential distribution instead (p = 0.011). When the far periphery was excluded, the R2 fit for the exponential distribution was 0.9618. A computational model following the principles of self-organized criticality (SOC), akin to earthquake and forest fire models, matched these datasets.

CONCLUSIONS

These distributions inform what useful statistics may be applied to study of these imaging characteristics. Further, the difference in event distribution between NV and NP suggests that the two phenomena are mechanistically distinct. NV may follow SOC, propagating as a catastrophic event in an unpredictable manner.

Anomalous Self-Organization in Active Piles

Inspired by recent observations on active self-organized critical (SOC) systems, we designed an active pile (or ant pile) model with two ingredients: beyond-threshold toppling and under-threshold active motions. By including the latter component, we were able to replace the typical power-law distribution for geometric observables with a stretched exponential fat-tailed distribution, where the exponent and decay rate are dependent on the activity's strength (ζ). This observation helped us to uncover a hidden connection between active SOC systems and α-stable Levy systems. We demonstrate that one can partially sweep α-stable Levy distributions by changing ζ. The system undergoes a crossover towards Bak-Tang-Weisenfeld (BTW) sandpiles with a power-law behavior (SOC fixed point) below a crossover point ζ<ζ*≈0.1.

Scaling Behavior of Quasi-One-Dimensional Vortex Avalanches in Superconducting Films

Scaling behaviour of dynamically driven vortex avalanches in superconducting YBa2Cu3O7-δ films deposited on tilted crystalline substrates has been observed using quantitative magneto-optical imaging. Two films with different tilt angles are characterized by the probability distributions of avalanche size in terms of the number of moving vortices. It is found in both samples that these distributions follow power-laws over up to three decades, and have exponents ranging between 1.0 and 1.4. The distributions also show clear finite-size scaling, when the system size is defined by the depth of the flux penetration front - a signature of self-organized criticality. A scaling relation between the avalanche size exponent and the fractal dimension, previously derived theoretically from conservation of the number of magnetic vortices in the stationary state and shown in numerical simulations, is here shown to be satisfied also experimentally.

Observation of Two-Dimensional Localized Jones-Roberts Solitons in Bose-Einstein Condensates

Jones-Roberts solitons are the only known class of stable dark solitonic solutions of the nonlinear Schrödinger equation in two and three dimensions. They feature a distinctive elongated elliptical shape that allows them to travel without change of form. By imprinting a triangular phase pattern, we experimentally generate two-dimensional Jones-Roberts solitons in a three-dimensional atomic Bose-Einstein condensate. We monitor their dynamics, observing that this kind of soliton is indeed not affected by dynamic (snaking) or thermodynamic instabilities, that instead make other classes of dark solitons unstable in dimensions higher than one. Our results confirm the prediction that Jones-Roberts solitons are stable solutions of the nonlinear Schrödinger equation and promote them for applications beyond matter wave physics, like energy and information transport in noisy and inhomogeneous environments.

Crossover from growing to stationary interfaces in the Kardar-Parisi-Zhang class

This Letter reports on how the interfaces in the (1+1)-dimensional Kardar-Parisi-Zhang (KPZ) class undergo, in the course of time, a transition from the flat, growing regime to the stationary one. Simulations of the polynuclear growth model and experiments on turbulent liquid crystal reveal universal functions of the KPZ class governing this transition, which connect the distribution and correlation functions for the growing and stationary regimes. This in particular shows how interfaces realized in experiments and simulations actually approach the stationary regime, which is never attained unless a stationary interface is artificially given as an initial condition.

Stratification, segregation, and mixing of granular materials in quasi-two-dimensional bounded heaps

Segregation and mixing of granular mixtures during heap formation has important consequences in industry and agriculture. This research investigates three different final particle configurations of bidisperse granular mixtures--stratified, segregated and mixed--during filling of quasi-two-dimensional silos. We consider a large number and wide range of control parameters, including particle size ratio, flow rate, system size, and heap rise velocity. The boundary between stratified and unstratified states is primarily controlled by the two-dimensional flow rate, with the critical flow rate for the transition depending weakly on particle size ratio and flowing layer length. In contrast, the transition from segregated to mixed states is controlled by the rise velocity of the heap, a control parameter not previously considered. The critical rise velocity for the transition depends strongly on the particle size ratio.

Predictability of the coherent-noise model and its applications

We study the threshold distribution function of the coherent-noise model for the case of infinite number of agents. This function is piecewise constant with a finite number of steps n. The latter exhibits a 1/f behavior as a function of the order of occurrence of an avalanche and hence versus natural time. An analytic expression of the expectation value E(S) for the size S of the next avalanche is obtained and used for the prediction of the next avalanche. Apart from E(S), the number of steps n can also serve for this purpose. This enables the construction of a similar prediction scheme which can be applied to real earthquake aftershock data.

Relation between self-organized criticality and grain aspect ratio in granular piles

We investigate experimentally whether self-organized criticality (SOC) occurs in granular piles composed of different grains, namely, rice, lentils, quinoa, and mung beans. These four grains were selected to have different aspect ratios, from oblong to oblate. As a function of aspect ratio, we determined the growth (β) and roughness (α) exponents, the avalanche fractal dimension (D), the avalanche size distribution exponent (τ), the critical angle (γ), and its fluctuation. At superficial inspection, three types of grains seem to have power-law-distributed avalanches with a well-defined τ. However, only rice is truly SOC if we take three criteria into account: a power-law-shaped avalanche size distribution, finite size scaling, and a universal scaling relation relating characteristic exponents. We study SOC as a spatiotemporal fractal; in particular, we study the spatial structure of criticality from local observation of the slope angle. From the fluctuation of the slope angle we conclude that greater fluctuation (and thus bigger avalanches) happen in piles consisting of grains with larger aspect ratio.

Avalanche spatial structure and multivariable scaling functions: sizes, heights, widths, and views through windows

We introduce a systematic method for extracting multivariable universal scaling functions and critical exponents from data. We exemplify our insights by analyzing simulations of avalanches in an interface using simulations from a driven quenched Kardar-Parisi-Zhang (qKPZ) equation. We fully characterize the spatial structure of these avalanches--we report universal scaling functions for size, height, and width distributions, and also local front heights. Furthermore, we resolve a problem that arises in many imaging experiments of crackling noise and avalanche dynamics, where the observed distributions are strongly distorted by a limited field of view. Through artificially windowed data, we show these distributions and their multivariable scaling functions may be written in terms of two control parameters: the window size and the characteristic length scale of the dynamics. For the entire system and the windowed distributions we develop accurate parametrizations for the universal scaling functions, including corrections to scaling and systematic error bars, facilitated by a novel software environment SloppyScaling.

Casimir-like effect on a granular pile

We investigate experimentally a Casimir-like effect in a three-dimensional pile of rice, which has a power-law avalanche size distribution. We observe the change in distance between two Plexiglas sheets placed on the pile parallel to each other and parallel to the mean avalanche flow direction, while rice grains are continuously and uniformly falling on top of the pile. The resulting avalanches are fluctuations, confinement of which is found to drive the two plates together. During 25-h experimental runs, for initial intersheet distances ranging from 20.0 to 90.0 mm we observe changes in the range from 6.0 mm to less than 1.0 mm. A similar distance dependence is obtained from a simple analytical model.

Avalanche prediction in a self-organized pile of beads

It is a common belief that power-law distributed avalanches are inherently unpredictable. This idea affects phenomena as diverse as evolution, earthquakes, superconducting vortices, stock markets, etc., from atomic to social scales. It mainly comes from the concept of "self-organized criticality" (SOC), where criticality is interpreted in the way that, at any moment, any small avalanche can eventually cascade into a large event. Nevertheless, this work demonstrates experimentally the possibility of avalanche prediction in the classical paradigm of SOC: a pile of grains. By knowing the position of every grain in a two-dimensional pile, avalanches of moving grains follow a distinct power-law distribution. Large avalanches, although uncorrelated, are on average preceded by continuous, detectable variations in the internal structure of the pile that are monitored in order to achieve prediction.

Edge effect on the power law distribution of granular avalanches

Many punctuated phenomena in nature are claimed [e.g., by the theory of self-organized criticality (SOC)] to be power-law distributed. In our experiments on a three-dimensional pile of long-grained rice, we find that by only changing the boundary condition of the system, we switch from such power-law-distributed avalanche sizes to quasiperiodic system-spanning avalanches. Conversely, by removing ledges the incidence of system-spanning avalanches is significantly reduced. This may offer a perspective on new avalanche prevention schemes. In addition, our findings may help to explain why the archetype of SOC, the sandpile, was found to have power-law-distributed avalanches in some experiments, while in other experiments quasiperiodic system-spanning avalanches were found.

Probability distribution of residence times of grains in models of rice piles

We study the probability distribution of residence time of a grain at a site, and its total residence time inside a pile, in different rice pile models. The tails of these distributions are dominated by the grains that get deeply buried in the pile. We show that, for a pile of size L, the probabilities that the residence time at a site or the total residence time is greater than t, both decay as 1/t(ln t)x for L(omega) << t << exp(L(gamma)) where gamma is an exponent > or = 1, and values of x and omega in the two cases are different. In the Oslo rice pile model we find that the probability of the residence time T(i) at a site i being greater than or equal to t is a nonmonotonic function of L for a fixed t and does not obey simple scaling. For model in d dimensions, we show that the probability of minimum slope configuration in the steady state, for large L, varies as exp(-kappaL(d+2)) where kappa is a constant, and hence gamma=d+2.

Extremal dynamics and the approach to the critical state: experiments on a three dimensional pile of rice

The evolution of the growth of a ricepile is studied in three dimensions. With time, the pile approaches a critical state with a certain slope. Assuming extremal dynamics in the evolution of the pile, the way the critical state is approached is dictated by the scaling properties of the critical state itself. Experimentally, we determine the envelope of the maximal slope, which is a measure for the distance from the critical state, as well as the growth of the average avalanche size with time. These quantities obey power-law scaling, where the experimental exponents are in good agreement with those obtained from an earlier determination of the critical state properties and extremal dynamics. Furthermore, we discuss the influence of the transient state on the avalanche size distribution, which may have applications in the prevention of large avalanches in natural systems.