Self-organized criticality induced by diversity

We have studied the collective behavior of a population of integrate-and-fire oscillators. We show that diversity, introduced in terms of a random distribution of natural periods, is the mechanism that permits one to observe self-organized criticality (SOC) in the long time regime. As diversity increases the system undergoes several transitions from a supercritical regime to a subcritical one, crossing the SOC region. Although there are resemblances with percolation, we give proofs that criticality takes place for a wide range of values of the control parameter instead of a single value.

Self-organized criticality and synchronization in a coupled map lattice model of integrate-and-fire oscillators

We introduce two coupled map lattice models with nonconservative interactions and a continuous nonlinear driving. Depending on both the degree of conservation and the convexity of the driving we find different behaviors, ranging from self-organized criticality, in the sense that the distribution of events (avalanches) obeys a power law, to a macroscopic synchronization of the population of oscillators, with avalanches of the size of the system.

Symmetries and fixed point stability of stochastic differential equations modeling self-organized criticality

A stochastic nonlinear partial differential equation is constructed for two different models exhibiting self-organized criticality: the Bak-Tang-Wiesenfeld (BTW) sandpile model [Phys. Rev. Lett. 59, 381 (1987); Phys. Rev. A 38, 364 (1988)] and the Zhang model [Phys. Rev. Lett. 63, 470 (1989)]. The dynamic renormalization group (DRG) enables one to compute the critical exponents. However, the nontrivial stable fixed point of the DRG transformation is unreachable for the original parameters of the models. We introduce an alternative regularization of the step function involved in the threshold condition, which breaks the symmetry of the BTW model. Although the symmetry properties of the two models are different, it is shown that they both belong to the same universality class. In this case the DRG procedure leads to a symmetric behavior for both models, restoring the broken symmetry, and makes accessible the nontrivial fixed point. This technique could also be applied to other problems with threshold dynamics.

Dynamical properties of the Zhang model of self-organized criticality

Critical exponents of the infinitely slowly driven Zhang model of self-organized criticality are computed for d=2 and 3, with particular emphasis devoted to the various roughening exponents. Besides confirming recent estimates of some exponents, new quantities are monitored, and their critical exponents computed. Among other results, it is shown that the three-dimensional exponents do not coincide with the Bak-Tang-Wiesenfeld [Phys. Rev. Lett. 59, 381 (1987); Phys. Rev. A 38, 364 (1988)] (Abelian) model, and that the dynamical exponent as computed from the correlation length and from the roughness of the energy profile do not necessarily coincide, as is usually implicitly assumed. An explanation for this is provided. The possibility of comparing these results with those obtained from renormalization group arguments is also briefly addressed.

Noise and dynamics of self-organized critical phenomena

Different microscopic models exhibiting self-organized criticality are studied numerically and analytically. Numerical simulations are performed to compute critical exponents, mainly the dynamical exponent, and to check universality classes. We find that various models lead to the same exponent, but this universality class is sensitive to disorder. From the dynamic microscopic rules we obtain continuum equations with different sources of noise, which we call internal and external. Different correlations of the noise give rise to different critical behavior. A model for external noise is proposed that makes the upper critical dimensionality equal to 4 and leads to the possible existence of a phase transition above d=4. Limitations of the approach of these models by a simple nonlinear equation are discussed.

Genesis of self-organized zebra textures in burial dolomites: Displacive veins, induced stress, and dolomitization

The dolomite veins making up rhythmites common in burial dolomites are not cement infillings of supposed cavities, as in the prevailing view, but are instead displacive veins, veins that pushed aside the host dolostone as they grew. Evidence that the veins are displacive includes a) small transform-fault-like displacements that could not have taken place if the veins were passive cements, and b) stylolites in host rock that formed as the veins grew in order to compensate for the volume added by the veins. Each zebra vein consists of crystals that grow inward from both sides, and displaces its walls via the local induced stress generated by the crystal growth itself. The petrographic criterion used in recent literature to interpret zebra veins in dolomites as cements - namely, that euhedral crystals can grow only in a prior void - disregards evidence to the contrary. The idea that flat voids did form in dolostones is incompatible with the observed optical continuity between the saddle dolomite euhedra of a vein and the replacive dolomite crystals of the host. The induced stress is also the key to the self-organization of zebra veins: In a set of many incipient, randomly-spaced, parallel veins just starting to grow in a host dolostone, each vein¿s induced stress prevents too-close neighbor veins from nucleating, or redissolves them by pressure-solution. The veins that survive this triage are those just outside their neighbors¿s induced stress haloes, now forming a set of equidistant veins, as observed.

Giant step bunching from self-organized coalescence of SrRuO3 islands

Step bunching develops in the epitaxy of SrRuO3 on vicinal SrTiO3(001) substrates. We have investigated the formation mechanisms and we show here that step bunching forms by lateral coalescence of wedgelike three-dimensional islands that are nucleated at substrate steps. After coalescence, wedgelike islands become wider and straighter with growth, forming a self-organized network of parallel step bunches with altitudes exceeding 30 unit cells, separated by atomically flat terraces. The formation mechanism of step bunching in SrRuO3, from nucleated islands, radically differs from one-dimensional models used to describe bunching in semiconducting materials. These results illustrate that growth phenomena of complex oxides can be dramatically different to those in semiconducting or metallic systems.

Long-tailed trapping times and Lévy flights in a self-organized critical granular system

We present a continuous time random walk model for the scale-invariant transport found in a self-organized critical rice pile [K. Christensen et al., Phys. Rev. Lett. 77, 107 (1996)]. From our analytical results it is shown that the dynamics of the experiment can be explained in terms of Lvy flights for the grains and a long-tailed distribution of trapping times. Scaling relations for the exponents of these distributions are obtained. The predicted microscopic behavior is confirmed by means of a cellular automaton model.

Self-organized evolution in a socio-economic environment

We propose a general scenario to analyze technological changes in socio-economic environments. We illustrate the ideas with a model that incorporating the main trends is simple enough to extract analytical results and, at the same time, sufficiently complex to display a rich dynamic behavior. Our study shows that there exists a macroscopic observable that is maximized in a regime where the system is critical, in the sense that the distribution of events follow power laws. Computer simulations show that, in addition, the system always self-organizes to achieve the optimal performance in the stationary state.

What determines the wavelength of self-organized shoreline sand waves?

Shoreline undulations extending into the bathymetric contours with a length scale larger than that of the rhythmic surf zone bars are referred to as shoreline sand waves. Many observed undulations along sandy coasts display a wavelength in the order 1-7 km. Several models that are based on the hypothesis that sand waves emerge from a morphodynamic instability in case of very oblique wave incidence predict this range of wavelengths. Here we investigate the physical reasons for the wavelength selection and the main parametric trends of the wavelength in case of sand waves arising from such instability. It is shown that the existence of a minimum wavelength depends on an interplay between three factors affecting littoral drift: (A) the angle of wave fronts relative to local shoreline, which tends to cause maximum transport at the downdrift flank of the sand wave, (B) the refractive energy spreading which tends to cause maximum transport at the updrift flank and (C) wave focusing (de-focusing) by the capes (bays), which tends to cause maximum transport at the crest or slightly downdrift of it. Processes A and C cause decay of the sand waves while process B causes their growth. For low incidence angles, B is very weak so that a rectilinear shoreline is stable. For large angles and long sand waves, B is dominant and causes the growth of sand waves. For large angles and short sand waves C is dominant and the sand waves decay. Thus, wavelength selection depends on process C, which essentially depends on shoreline curvature. The growth rate of very long sand waves is weak because the alongshore gradients in sediment transport decrease with the wavelength. This is why there is an optimum or dominant wavelength. It is found that sand wave wavelength scales with λ0/β where λ0 is the water wave wavelength in deep water and β is the mean bed slope from shore to the wave base.

Self-accelerating dolomite-for-calcite replacement: Self-organized dynamics of burial dolomitization and associated mineralization

A new dynamic model of dolomitization predicts a multitude of textural, paragenetic, geochemical and other properties of burial dolomites. The model is based on two postulates, (1) that the dolomitizing brine is Mg-rich but under saturated with both calcite and dolomite, and (2) that the dolomite-for-calcite replacement happens not by dissolution-precipitation as usually assumed, but by dolomite-growth-driven pressure solution of the calcite host. Crucially, the dolomite-for-calcite replacement turns out to be self-accelerating via Ca2 : the Ca2 released by each replacement increment accelerates the rate of the next, and so on. As a result, both pore-fluid Ca2 and replacement rate grow exponentially.

Hyperbolic reaction-diffusion equations for a forest fire model

Forest fire models have been widely studied from the context of self-organized criticality and from the ecological properties of the forest and combustion. On the other hand, reaction-diffusion equations have interesting applications in biology and physics. We propose here a model for fire propagation in a forest by using hyperbolic reaction-diffusion equations. The dynamical and thermodynamical aspects of the model are analyzed in detail

Synchronization in a lattice model of pulse-coupled oscillators

We analyze the collective behavior of a lattice model of pulse-coupled oscillators. By means of computer simulations we find the relation between the intrinsic dynamics of each member of the population and their mutual interactions that ensures, in a general context, the existence of a fully synchronized regime. This condition turns out to be the same as that obtained for the globally coupled population. When the condition is not completely satisfied we find different spatial structures. This also gives some hints about self-organized criticality.

Invaded cluster dynamics for frustrated models

The invaded cluster (IC) dynamics introduced by Machta et al. [Phys. Rev. Lett. 75, 2792 (1995)] is extended to the fully frustrated Ising model on a square lattice. The properties of the dynamics that exhibits numerical evidence of self-organized criticality are studied. The fluctuations in the IC dynamics are shown to be intrinsic of the algorithm and the fluctuation-dissipation theorem is no longer valid. The relaxation time is found to be very short and does not present a critical size dependence.