Avalanches and the renormalization group for pinned charge-density waves

Theoretical foundations of studying criticality in the brain

Criticality is hypothesized as a physical mechanism underlying efficient transitions between cortical states and remarkable information-processing capacities in the brain. While considerable evidence generally supports this hypothesis, nonnegligible controversies persist regarding the ubiquity of criticality in neural dynamics and its role in information processing. Validity issues frequently arise during identifying potential brain criticality from empirical data. Moreover, the functional benefits implied by brain criticality are frequently misconceived or unduly generalized. These problems stem from the nontriviality and immaturity of the physical theories that analytically derive brain criticality and the statistic techniques that estimate brain criticality from empirical data. To help solve these problems, we present a systematic review and reformulate the foundations of studying brain criticality, that is, ordinary criticality (OC), quasi-criticality (qC), self-organized criticality (SOC), and self-organized quasi-criticality (SOqC), using the terminology of neuroscience. We offer accessible explanations of the physical theories and statistical techniques of brain criticality, providing step-by-step derivations to characterize neural dynamics as a physical system with avalanches. We summarize error-prone details and existing limitations in brain criticality analysis and suggest possible solutions. Moreover, we present a forward-looking perspective on how optimizing the foundations of studying brain criticality can deepen our understanding of various neuroscience questions.

Theory and experiments for disordered elastic manifolds, depinning, avalanches, and sandpiles

Domain walls in magnets, vortex lattices in superconductors, contact lines at depinning, and many other systems can be modeled as an elastic system subject to quenched disorder. The ensuing field theory possesses a well-controlled perturbative expansion around its upper critical dimension. Contrary to standard field theory, the renormalization group (RG) flow involves a function, the disorder correlator Δ(w), and is therefore termed the functional RG. Δ(w) is a physical observable, the auto-correlation function of the center of mass of the elastic manifold. In this review, we give a pedagogical introduction into its phenomenology and techniques. This allows us to treat both equilibrium (statics), and depinning (dynamics). Building on these techniques, avalanche observables are accessible: distributions of size, duration, and velocity, as well as the spatial and temporal shape. Various equivalences between disordered elastic manifolds, and sandpile models exist: an elastic string driven at a point and the Oslo model; disordered elastic manifolds and Manna sandpiles; charge density waves and Abelian sandpiles or loop-erased random walks. Each of the mappings between these systems requires specific techniques, which we develop, including modeling of discrete stochastic systems via coarse-grained stochastic equations of motion, super-symmetry techniques, and cellular automata. Stronger than quadratic nearest-neighbor interactions lead to directed percolation, and non-linear surface growth with additional Kardar-Parisi-Zhang (KPZ) terms. On the other hand, KPZ without disorder can be mapped back to disordered elastic manifolds, either on the directed polymer for its steady state, or a single particle for its decay. Other topics covered are the relation between functional RG and replica symmetry breaking, and random-field magnets. Emphasis is given to numerical and experimental tests of the theory.

Depinning Transition of Charge-Density Waves: Mapping onto O(n) Symmetric ϕ^{4} Theory with n→-2 and Loop-Erased Random Walks

Driven periodic elastic systems such as charge-density waves (CDWs) pinned by impurities show a nontrivial, glassy dynamical critical behavior. Their proper theoretical description requires the functional renormalization group. We show that their critical behavior close to the depinning transition is related to a much simpler model, O(n) symmetric ϕ^{4} theory in the unusual limit of n→-2. We demonstrate that both theories yield identical results to four-loop order and give both a perturbative and a nonperturbative proof of their equivalence. As we show, both theories can be used to describe loop-erased random walks (LERWs), the trace of a random walk where loops are erased as soon as they are formed. Remarkably, two famous models of non-self-intersecting random walks, self-avoiding walks and LERWs, can both be mapped onto ϕ^{4} theory, taken with formally n=0 and n→-2 components. This mapping allows us to compute the dynamic critical exponent of CDWs at the depinning transition and the fractal dimension of LERWs in d=3 with unprecedented accuracy, z(d=3)=1.6243±0.001, in excellent agreement with the estimate z=1.62400±0.00005 of numerical simulations.

Crossover behavior in interface depinning

We study the crossover scaling behavior of the height-height correlation function in interface depinning in random media. We analyze experimental data from a fracture experiment and simulate an elastic line model with nonlinear couplings and disorder. Both exhibit a crossover between two different universality classes. For the experiment, we fit a functional form to the universal crossover scaling function. For the model, we vary the system size and the strength of the nonlinear term and describe the crossover between the two universality classes with a multiparameter scaling function. Our method provides a general strategy to extract scaling properties in depinning systems exhibiting crossover phenomena.

Exact mapping of the stochastic field theory for Manna sandpiles to interfaces in random media

We show that the stochastic field theory for directed percolation in the presence of an additional conservation law [the conserved directed-percolation (C-DP) class] can be mapped exactly to the continuum theory for the depinning of an elastic interface in short-range correlated quenched disorder. Along one line of the parameters commonly studied, this mapping leads to the simplest overdamped dynamics. Away from this line, an additional memory term arises in the interface dynamics; we argue that this does not change the universality class. Since C-DP is believed to describe the Manna class of self-organized criticality, this shows that Manna stochastic sandpiles and disordered elastic interfaces (i.e., the quenched Edwards-Wilkinson model) share the same universal large-scale behavior.

Scaling description of the yielding transition in soft amorphous solids at zero temperature

Yield stress materials flow if a sufficiently large shear stress is applied. Although such materials are ubiquitous and relevant for industry, there is no accepted microscopic description of how they yield, even in the simplest situations in which temperature is negligible and in which flow inhomogeneities such as shear bands or fractures are absent. Here we propose a scaling description of the yielding transition in amorphous solids made of soft particles at zero temperature. Our description makes a connection between the Herschel-Bulkley exponent characterizing the singularity of the flow curve near the yield stress Σc, the extension and duration of the avalanches of plasticity observed at threshold, and the density P(x) of soft spots, or shear transformation zones, as a function of the stress increment x beyond which they yield. We argue that the critical exponents of the yielding transition may be expressed in terms of three independent exponents, θ, df, and z, characterizing, respectively, the density of soft spots, the fractal dimension of the avalanches, and their duration. Our description shares some similarity with the depinning transition that occurs when an elastic manifold is driven through a random potential, but also presents some striking differences. We test our arguments in an elasto-plastic model, an automaton model similar to those used in depinning, but with a different interaction kernel, and find satisfying agreement with our predictions in both two and three dimensions.

Avalanche dynamics of elastic interfaces

Slowly driven elastic interfaces, such as domain walls in dirty magnets, contact lines wetting a nonhomogeneous substrate, or cracks in brittle disordered material proceed via intermittent motion, called avalanches. Here we develop a field-theoretic treatment to calculate, from first principles, the space-time statistics of instantaneous velocities within an avalanche. For elastic interfaces at (or above) their (internal) upper critical dimension d≥d(uc) (d(uc)=2,4 respectively for long-ranged and short-ranged elasticity) we show that the field theory for the center of mass reduces to the motion of a point particle in a random-force landscape, which is itself a random walk [Alessandro, Beatrice, Bertotti, and Montorsi (ABBM) model]. Furthermore, the full spatial dependence of the velocity correlations is described by the Brownian-force model (BFM) where each point of the interface sees an independent Brownian-force landscape. Both ABBM and BFM can be solved exactly in any dimension d (for monotonous driving) by summing tree graphs, equivalent to solving a (nonlinear) instanton equation. We focus on the limit of slow uniform driving. This tree approximation is the mean-field theory (MFT) for realistic interfaces in short-ranged disorder, up to the renormalization of two parameters at d=d(uc). We calculate a number of observables of direct experimental interest: Both for the center of mass, and for a given Fourier mode q, we obtain various correlations and probability distribution functions (PDF's) of the velocity inside an avalanche, as well as the avalanche shape and its fluctuations (second shape). Within MFT we find that velocity correlations at nonzero q are asymmetric under time reversal. Next we calculate, beyond MFT, i.e., including loop corrections, the one-time PDF of the center-of-mass velocity u[over ·] for dimension d<d(uc). The singularity at small velocity P(u[over ·])~1/u[over ·](a) is substantially reduced from a=1 (MFT) to a=1-2/9(4-d)+... (short-ranged elasticity) and a=1-4/9(2-d)+... (long-ranged elasticity). We show how the dynamical theory recovers the avalanche-size distribution, and how the instanton relates to the response to an infinitesimal step in the force.

Quantum mechanisms of density wave transport

We report on new developments in the quantum picture of correlated electron transport in charge and spin density waves. The model treats the condensate as a quantum fluid in which charge soliton domain wall pairs nucleate above a Coulomb blockade threshold field. We employ a time-correlated soliton tunneling model, analogous to the theory of time-correlated single electron tunneling, to interpret the voltage oscillations and nonlinear current-voltage characteristics above threshold. An inverse scaling relationship between threshold field and dielectric response, originally proposed by Grüner, emerges naturally from the model. Flat dielectric and other ac responses below threshold in NbSe(3) and TaS(3), as well as small density wave phase displacements, indicate that the measured threshold is often much smaller than the classical depinning field. In some materials, the existence of two distinct threshold fields suggests that both soliton nucleation and classical depinning may occur. In our model, the ratio of electrostatic charging to pinning energy helps determine whether soliton nucleation or classical depinning dominates.

Evidence for a Peierls phase-transition in a three-dimensional multiple charge-density waves solid

The effect of dimensionality on materials properties has become strikingly evident with the recent discovery of graphene. Charge ordering phenomena can be induced in one dimension by periodic distortions of a material's crystal structure, termed Peierls ordering transition. Charge-density waves can also be induced in solids by strong coulomb repulsion between carriers, and at the extreme limit, Wigner predicted that crystallization itself can be induced in an electrons gas in free space close to the absolute zero of temperature. Similar phenomena are observed also in higher dimensions, but the microscopic description of the corresponding phase transition is often controversial, and remains an open field of research for fundamental physics. Here, we photoinduce the melting of the charge ordering in a complex three-dimensional solid and monitor the consequent charge redistribution by probing the optical response over a broad spectral range with ultrashort laser pulses. Although the photoinduced electronic temperature far exceeds the critical value, the charge-density wave is preserved until the lattice is sufficiently distorted to induce the phase transition. Combining this result with ab initio electronic structure calculations, we identified the Peierls origin of multiple charge-density waves in a three-dimensional system for the first time.

Seismic cycles, size of the largest events, and the avalanche size distribution in a model of seismicity

We address several questions on the behavior of a numerical model recently introduced to study seismic phenomena, which includes relaxation in the plates as a key ingredient. First, we make an analysis of the scaling of the largest events with system size and show that, when parameters are appropriately interpreted, the typical size of the largest events scale as the system size, without the necessity to tune any parameter. Second, we show that the temporal activity in the model is inherently nonstationary and obtain from here justification and support for the concept of a "seismic cycle" in the temporal evolution of seismic activity. Finally, we ask for the reasons that make the model display a realistic value of the decaying exponent b in the Gutenberg-Richter law for the avalanche size distribution. We explain why relaxation induces a systematic increase in b from its value b≃0.4 observed in the absence of relaxation. However, we have not been able to justify the actual robustness of the model in displaying a consistent b value around the experimentally observed value b≃1.

Minimal integer automaton behind crystal plasticity

Power law fluctuations and scale-free spatial patterns are known to characterize steady state plastic flow in crystalline materials. In this Letter we study the emergence of correlations in a simple Frenkel-Kontorova-type model of 2D plasticity which is largely free of arbitrariness, amenable to analytical study, and is capable of generating critical exponents matching experiments. Our main observation concerns the possibility to reduce continuum plasticity to an integer-valued automaton revealing inherent discreteness of the plastic flow.

Statistics of static avalanches in a random pinning landscape

We study the minimum-energy configuration of a d -dimensional elastic interface in a random potential tied to a harmonic spring. As a function of the spring position, the center of mass of the interface changes in discrete jumps, also called shocks or "static avalanches." We obtain analytically the distribution of avalanche sizes and its cumulants within an =4-d expansion from a tree and one-loop resummation using functional renormalization. This is compared with exact numerical minimizations of interface energies for random-field disorder in d=2,3 . Connections to dynamic avalanches are mentioned.

Size distributions of shocks and static avalanches from the functional renormalization group

Interfaces pinned by quenched disorder are often used to model jerky self-organized critical motion. We study static avalanches, or shocks, defined here as jumps between distinct global minima upon changing an external field. We show how the full statistics of these jumps is encoded in the functional-renormalization-group fixed-point functions. This allows us to obtain the size distribution P(S) of static avalanches in an expansion in the internal dimension d of the interface. Near and above d=4 this yields the mean-field distribution P(S) approximately S;{-3/2}e;{-S4S_{m}} , where S_{m} is a large-scale cutoff, in some cases calculable. Resumming all one-loop contributions, we find P(S) approximately S;{-tau}exp(C(SS_{m});{1/2}-B/4(S/S_{m});{delta}) , where B , C , delta , and tau are obtained to first order in =4-d . Our result is consistent to O() with the relation tau=tau_{zeta}:=2-2/d+zeta , where zeta is the static roughness exponent, often conjectured to hold at depinning. Our calculation applies to all static universality classes, including random-bond, random-field, and random-periodic disorders. Extended to long-range elastic systems, it yields a different size distribution for the case of contact-line elasticity, with an exponent compatible with tau=2-1/d+zeta to O(=2-d) . We discuss consequences for avalanches at depinning and for sandpile models, relations to Burgers turbulence and the possibility that the relation tau=tau_{zeta} be violated to higher loop order. Finally, we show that the avalanche-size distribution on a hyperplane of codimension one is in mean field (valid close to and above d=4 ) given by P(S) approximately K_{13}(S)S , where K is the Bessel- K function, thus tau_{hyperplane}=4/3 .

Cusps, self-organization, and absorbing states

Elastic interfaces embedded in (quenched) random media exhibit metastability and stick-slip dynamics. These nontrivial dynamical features have been shown to be associated with cusp singularities of the coarse-grained disorder correlator. Here we show that annealed systems with many absorbing states and a conservation law but no quenched disorder exhibit identical cusps. On the other hand, similar nonconserved systems in the directed percolation class are also shown to exhibit cusps but of a different type. These results are obtained both by a recent method to explicitly measure disorder correlators and by defining an alternative new protocol inspired by self-organized criticality, which opens the door to easily accessible experimental realizations.

Confirming and extending the hypothesis of universality in sandpiles

Stochastic sandpiles self-organize to an absorbing-state critical point with scaling behavior different from directed percolation (DP) and characterized by the presence of an additional conservation law. This is usually called the C-DP or Manna universality class. There remains, however, an exception to this universality principle: a sandpile automaton introduced by Maslov and Zhang, which was claimed to be in the DP class despite the existence of a conservation law. We show, by means of careful numerical simulations as well as by constructing and analyzing a field theory, that (contrarily to what was previously thought) this sandpile is also in the C-DP or Manna class. This confirms the hypothesis of universality for stochastic sandpiles and gives rise to a fully coherent picture of self-organized criticality in systems with conservation. In passing, we obtain a number of results for the C-DP class and introduce a strategy to easily discriminate between DP and C-DP scaling.

Absorbing states and elastic interfaces in random media: two equivalent descriptions of self-organized criticality

We elucidate a long-standing puzzle about the nonequilibrium universality classes describing self-organized criticality in sandpile models. We show that depinning transitions of linear interfaces in random media and absorbing phase transitions (with a conserved nondiffusive field) are two equivalent languages to describe sandpile criticality. This is so despite the fact that local roughening properties can be radically different in the two pictures, as explained here. Experimental implications of our work as well as promising paths for future theoretical investigations are also discussed.

Irrational mode locking in quasiperiodic systems

A model for ac-driven systems, based on the Tang-Wiesenfeld-Bak-Coppersmith-Littlewood automaton for an elastic medium, exhibits mode-locked steps with frequencies that are irrational multiples of the drive frequency, when the pinning is spatially quasiperiodic. Detailed numerical evidence is presented for the large-system-size convergence of such a mode-locked step. The irrational mode locking is stable to small thermal noise and weak disorder. Continuous-time models with irrational mode locking and possible experimental realizations are discussed.

Dynamics below the depinning threshold in disordered elastic systems

We study the steady-state low-temperature dynamics of an elastic line in a disordered medium below the depinning threshold. Analogously to the equilibrium dynamics, in the limit T-->0, the steady state is dominated by a single configuration which is occupied with probability 1. We develop an exact algorithm to target this dominant configuration and to analyze its geometrical properties as a function of the driving force. The roughness exponent of the line at large scales is identical to the one at depinning. No length scale diverges in the steady-state regime as the depinning threshold is approached from below. We do find a divergent length, but it is associated only with the transient relaxation between metastable states.

Exact solution of the one-dimensional deterministic fixed-energy sandpile

By reason of the strongly nonergodic dynamical behavior, universality properties of deterministic fixed-energy sandpiles are still an open and debated issue. We investigate the one-dimensional model, whose microscopical dynamics can be solved exactly, and provide a deeper understanding of the origin of the nonergodicity. By means of exact arguments, we prove the occurrence of orbits of well-defined periods and their dependence on the conserved energy density. Further statistical estimates of the size of the attraction's basins of the different periodic orbits lead to a complete characterization of the activity vs energy density phase diagram in the limit of large system's size.

Interface depinning versus absorbing-state phase transitions

According to recent numerical results from lattice models, the critical exponents of systems with many absorbing states and order parameter coupled to a nondiffusive conserved field coincide with those of the linear interface depinning model within computational accuracy. In this paper the connection between absorbing-state phase transitions and interface pinning in quenched disordered media is investigated. For that, we present an heuristic mapping of the interface dynamics in a disordered medium into a Langevin equation for the active-site density and show that a Reggeon-field-theory-like description, in which the order parameter appears coupled to an additional nondiffusive conserved field, emerges rather naturally. Reciprocally, we construct a mapping from a discrete model belonging in the absorbing state with a conserved-field class to a discrete interface equation, and show how a quenched disorder, typical of the interface representation is originated. We discuss the character of the possible noise terms in both representations, and overview the critical exponent relations. Evidence is provided that, at least for dimensions larger that one, both universality classes are just two different representations of the same underlying physics.

Depinning with dynamic stress overshoots: mean field theory

An infinite-range model of an elastic manifold pulled through a random potential by a force F is analyzed focusing on inertial effects. When the inertial parameter M is small, there is a continuous depinning transition from a small- F static phase to a large- F moving phase. When M is increased to M(c), a novel tricritical point occurs. For M>M(c), the depinning transition becomes discontinuous with hysteresis. Yet, the distribution of discrete "avalanche"-like events as the force is increased in the static phase for M>M(c) has an unusual mixture of first-order-like and critical features. The results may be relevant for the onset of crack propagation and for dynamics of geological faults.

Aging and self-organized criticality in driven dissipative systems

We study the noisy dynamics of a close relative to the original sandpile model. Depending on the type of noise and the time scale of observation, we find stationary fluctuations (similar to self-organized criticality) or an aging dynamics with punctuated equilibria, a decreasing rate of events and reset properties qualitatively similar to those of glassy systems, evolution models, and vibrated granular media.

First-passage-time exponent for higher-order random walks: using Lévy flights

We present a heuristic derivation of the first-passage-time exponent for the integral of a random walk [Y. G. Sinai, Theor. Math. Phys. 90, 219 (1992)]. Building on this derivation, we construct an estimation scheme to understand the first-passage-time exponent for the integral of the integral of a random walk, which is numerically observed to be 0.220+/-0.001. We discuss the implications of this estimation scheme for the nth integral of a random walk. For completeness, we also address the n=infinity case. Finally, we explore an application of these processes to an extended, elastic object being pulled through a random potential by a uniform applied force. In so doing, we demonstrate a time reparametrization freedom in the Langevin equation that maps nonlinear stochastic processes into linear ones.

Anomalous scaling in depinning transitions

It is demonstrated that the renormalization group (RG) flows of depinning transitions do not depend on whether the driving force or the system velocity is kept constant. This allows for a comparison between RG results and corresponding self-organized critical models. However, close to the critical point, scaling functions cross over to forms that can have singular behavior not seen in equilibrium thermal phase transitions. These can be different for the constant force and constant velocity driving modes, leading to different apparent critical exponents. This is illustrated by comparing extremal dynamics for interface depinning with RG results, deriving the change in apparent exponents. Thus, care has to be exercised in such comparisons.

Interface scaling in the contact process

A correspondence between lattice models with absorbing states and models of pinned interfaces in random media can be established by defining local height variables h(x,t) as integrals of the activity at point x up to time t. Within this context we study the interface representation of a prototypical model with absorbing states, the contact process, in dimensions 1-3. Simulations confirm the scaling relation beta(W)=1-straight theta between the interface-width growth exponent beta(W) and the exponent straight theta governing the decay of the order parameter. A scaling property of the height distribution, which serves as the basis for this relation, is also verified. The height-height correlation function shows clear signs of anomalous scaling, in accord with Lopez' analysis [Phys. Rev. Lett. 83, 4594 (1999)], but no evidence of multiscaling.

Absorbing-state phase transitions in fixed-energy sandpiles

We study sandpile models as closed systems, with the conserved energy density zeta playing the role of an external parameter. The critical energy density zeta(c) marks a nonequilibrium phase transition between active and absorbing states. Several fixed-energy sandpiles are studied in extensive simulations of stationary and transient properties, as well as the dynamics of roughening in an interface-height representation. Our primary goal is to identify the universality classes of such models, in hopes of assessing the validity of two recently proposed approaches to sandpiles: a phenomenological continuum Langevin description with absorbing states, and a mapping to driven interface dynamics in random media.

Time-correlated soliton tunneling in charge and spin density waves

We consider a model in which an electric field induces quantum nucleation of kink-antikink pairs in a pinned charge or spin density wave. Pair nucleation events, prevented by Coulomb blockade below a pair creation threshold, become correlated in time above threshold. The model provides a natural explanation for the observed (i) small density wave polarization below threshold in NbSe (3), (ii) narrow band noise, (iii) coherent oscillations, and (iv) mode-locking at high drift frequencies.

Highly optimized tolerance: a mechanism for power laws in designed systems

We introduce a mechanism for generating power law distributions, referred to as highly optimized tolerance (HOT), which is motivated by biological organisms and advanced engineering technologies. Our focus is on systems which are optimized, either through natural selection or engineering design, to provide robust performance despite uncertain environments. We suggest that power laws in these systems are due to tradeoffs between yield, cost of resources, and tolerance to risks. These tradeoffs lead to highly optimized designs that allow for occasional large events. We investigate the mechanism in the context of percolation and sand pile models in order to emphasize the sharp contrasts between HOT and self-organized criticality (SOC), which has been widely suggested as the origin for power laws in complex systems. Like SOC, HOT produces power laws. However, compared to SOC, HOT states exist for densities which are higher than the critical density, and the power laws are not restricted to special values of the density. The characteristic features of HOT systems include: (1) high efficiency, performance, and robustness to designed-for uncertainties; (2) hypersensitivity to design flaws and unanticipated perturbations; (3) nongeneric, specialized, structured configurations; and (4) power laws. The first three of these are in contrast to the traditional hallmarks of criticality, and are obtained by simply adding the element of design to percolation and sand pile models, which completely changes their characteristics.

Self-organized criticality: sandpiles, singularities, and scaling

We present an overview of the statistical mechanics of self-organized criticality. We focus on the successes and failures of hydrodynamic description of transport, which consists of singular diffusion equations. When this description applies, it can predict the scaling features associated with these systems. We also identify a hard driving regime where singular diffusion hydrodynamics fails due to fluctuations and give an explicit criterion for when this failure occurs.