Scaling exponents for barkhausen avalanches in polycrystalline and amorphous ferromagnets

Uncovering avalanche sources via acceleration measurements

Avalanche sources describe rapid and local events that govern deformation processes in various materials. The fundamental differences between an avalanche source and its associated measured acoustic emission (AE) signal are encoded in the acoustic transfer function, which undesirably modifies the properties of the source. Consequently, information about the physical characteristics of avalanche sources is scarce and its exposure poses a great challenge. We introduce a novel experimental method based on acceleration measurements, which eliminates the effect of the transfer function and distills the avalanche source. Applying this method to deformation twinning in magnesium shows that the amplitudes and characteristic times of avalanche sources are unrelated by a clear physical law. Conversely, the amplitudes and durations of AE signals are related by a power law, which is attributed to the transfer function. Using our method, we identify and compute a new feature of avalanche sources, which is directly linked to the growth rate of the twinned volume. This feature displays a power-law distribution, implying an unpredicted behavior at dynamic criticality. Simultaneously, the characteristic times of avalanche sources possess an intrinsic upper bound, indicating a predicted limit that relates to the underlying physical process of twinning.

Simple model for the prediction of seizure durations

A simple model is used to simulate seizures in a population of spiking excitatory neurons experiencing a uniform effect from inhibitory neurons. A key feature is introduced into the model, i.e., a mechanism that weakens the firing thresholds. This weakening mechanism adds memory to the dynamics. We find a seizure-prone state in a "mode-switching" phase. In this phase, the system can suddenly switch from a "healthy" state with small scale-free avalanches to a "seizure" state with almost periodic large avalanches ("seizures"). Simulations of the model predict statistics for the average time spent in the seizure state (the seizure "duration") that agree with experiments and theoretical examples of similar behavior in neuronal systems. Our study points to. different connections between seizures and fracture and also offers an alternative view on the type of critical point controlling neuronal avalanches.

Acoustic emissions from spin crossover complexes

Acoustic emission from the compounds [Fe(HB(tz)3)2] and [Fe(Htrz)(trz)2]BF4 was detected during the thermally induced spin transition and is correlated with simultaneously recorded calorimetric signals. We ascribe this phenomenon to elastic waves produced by microstructural and volume changes accompanying the spin transition. Despite the perfect reversibility of the spin state switching (seen by the calorimeter), the acoustic emission activity decreases for successive thermal cycles, revealing thus irreversible microstructural evolution of the samples. The acoustic emission signal amplitude and energy probability distribution functions followed power-law behavior and the characteristic exponents were found to be similar for the two samples both on heating and cooling, indicating the universal character, which is further substantiated by the well scaled average temporal shapes of the avalanches.

Barkhausen noise in disordered striplike ferromagnets: Experiment versus simulations

In this work, we present a systematic comparison of the results obtained from the low-frequency Barkhausen noise recordings in nanocrystalline samples with those from the numerical simulations of the random-field Ising model systems. We performed measurements at room temperature on a field-driven metallic glass stripe made of VITROPERM 800 R, a nanocrystalline iron-based material with an excellent combination of soft and magnetic properties, making it a cutting-edge material for a wide range of applications. Given that the Barkhausen noise emissions emerging along a hysteresis curve are stochastic and depend in general on a variety of factors (such as distribution of disorder due to impurities or defects, varied size of crystal grains, type of domain structure, driving rate of the external magnetic field, sample shape and temperature, etc.), adequate theoretical modeling is essential for their interpretation and prediction. Here the Random field Ising model, specifically its athermal nonequilibrium version with the finite driving rate, stands out as an appropriate choice due to the material's nanocrystalline structure and high Curie temperature. We performed a systematic analysis of the signal properties and magnetization avalanches comparing the outcomes of the numerical model and experiments carried out in a two-decade-wide range of the external magnetic field driving rates. Our results reveal that with a suitable choice of parameters, a considerable match with the experimental results is achieved, indicating that this model can accurately describe the Barkhausen noise features in nanocrystalline samples.

Mechanical excitation and marginal triggering during avalanches in sheared amorphous solids

We study plastic strain during individual avalanches in overdamped particle-scale molecular dynamics (MD) and mesoscale elastoplastic models (EPM) for amorphous solids sheared in the athermal quasistatic limit. We show that the spatial correlations in plastic activity exhibit a short length scale that grows as t^{3/4} in MD and ballistically in EPM, which is generated by mechanical excitation of nearby sites not necessarily close to their stability thresholds, and a longer lengthscale that grows diffusively for both models and is associated with remote marginally stable sites. These similarities in spatial correlations explain why simple EPMs accurately capture the size distribution of avalanches observed in MD, though the temporal profiles and dynamical critical exponents are quite different.

Force Correlations in Disordered Magnets

We present a proof of principle for the validity of the functional renormalization group, by measuring the force correlations in Barkhausen-noise experiments. Our samples are soft ferromagnets in two distinct universality classes, differing in the range of spin interactions, and the effects of eddy currents. We show that the force correlations have a universal form predicted by the functional renormalization group, distinct for short-range and long-range elasticity, and mostly independent of eddy currents. In all cases correlations grow linearly at small distances, as in mean-field models, but in contrast to the latter are bounded at large distances. As a consequence, avalanches are anti-correlated. We derive bounds for these anticorrelations, which are saturated in the experiments, showing that the multiple domain walls in our samples effectively behave as a single wall.

Denouement of the Energy-Amplitude and Size-Amplitude Enigma for Acoustic-Emission Investigations of Materials

There are many systems producing crackling noise (avalanches) in materials. Temporal shapes of avalanches, U(t) (U is the detected voltage signal, t is the time), have self-similar behaviour and the normalized U(t) function (e.g., dividing both the values of U and t by S^1/2, where S is the avalanche area), averaged for fixed S, should be the same, independently of the type of materials or avalanche mechanisms. However, there are experimental evidences that the temporal shapes of avalanches do not scale completely in a universal way. The self-similarity also leads to universal power-law-scaling relations, e.g., between the energy, E, and the peak amplitude, A_m, or between S and A_m. There are well-known enigmas, where the above exponents in acoustic emission measurements are rather close to 2 and 1, respectively, instead of E~Am3 and S~Am2, obtained from the mean field theory, MFT. We show, using a theoretically predicted averaged function for the fixed avalanche area, U(t)=atexp(−bt2) (where a and b are non-universal, material-dependent constants), that the scaling exponents can be different from the MFT values. Normalizing U by A_m and t by t_m (the time belonging to the A_m: rise time), we obtain tm~Am1−φ (the MFT values can be obtained only if φ would be zero). Here, φ is expected to be material-independent and to be the same for the same mechanism. Using experimental results on martensitic transformations in two different shape-memory single-crystals, φ = 0.8 ± 0.1 was obtained (φ is the same for both alloys). Thus, dividing U by A_m as well as t by Am1−φ (~t_m) leads to the same common, normalized temporal shape for different, fixed values of S. This normalization can also be used in general for other experimental results (not only for acoustic emission), which provide information about jerky noises in materials.

Temporal evolution of failure avalanches of the fiber bundle model on complex networks

We investigate how the interplay of the topology of the network of load transmitting connections and the amount of disorder of the strength of the connected elements determines the temporal evolution of failure cascades driven by the redistribution of load following local failure events. We use the fiber bundle model of materials' breakdown assigning fibers to the sites of a square lattice, which is then randomly rewired using the Watts-Strogatz technique. Gradually increasing the rewiring probability, we demonstrate that the bundle undergoes a transition from the localized to the mean field universality class of breakdown phenomena. Computer simulations revealed that both the size and the duration of failure cascades are power law distributed on all network topologies with a crossover between two regimes of different exponents. The temporal evolution of cascades is described by a parabolic profile with a right handed asymmetry, which implies that cascades start slowly, then accelerate, and eventually stop suddenly. The degree of asymmetry proved to be characteristic of the network topology gradually decreasing with increasing rewiring probability. Reducing the variance of fibers' strength, the exponents of the size and the duration distribution of cascades increase in the localized regime of the failure process, while the localized to mean field transition becomes more abrupt. The consistency of the results is supported by a scaling analysis relating the characteristic exponents of the statistics and dynamics of cascades.

Thermally activated intermittent dynamics of creeping crack fronts along disordered interfaces

We present a subcritical fracture growth model, coupled with the elastic redistribution of the acting mechanical stress along rugous rupture fronts. We show the ability of this model to quantitatively reproduce the intermittent dynamics of cracks propagating along weak disordered interfaces. To this end, we assume that the fracture energy of such interfaces (in the sense of a critical energy release rate) follows a spatially correlated normal distribution. We compare various statistical features from the obtained fracture dynamics to that from cracks propagating in sintered polymethylmethacrylate (PMMA) interfaces. In previous works, it has been demonstrated that such an approach could reproduce the mean advance of fractures and their local front velocity distribution. Here, we go further by showing that the proposed model also quantitatively accounts for the complex self-affine scaling morphology of crack fronts and their temporal evolution, for the spatial and temporal correlations of the local velocity fields and for the avalanches size distribution of the intermittent growth dynamics. We thus provide new evidence that an Arrhenius-like subcritical growth is particularly suitable for the description of creeping cracks.

Effects of external noise on threshold-induced correlations in ferromagnetic systems

In the present paper we investigate the impact of the external noise and detection threshold level on the simulation data for the systems that evolve through metastable states. As a representative model of such systems we chose the nonequilibrium athermal random-field Ising model with two types of the external noise, uniform white noise and Gaussian white noise with various different standard deviations, imposed on the original response signal obtained in model simulations. We applied a wide range of detection threshold levels in analysis of the signal and show how these quantities affect the values of exponent γ_{S/T} (describing the scaling of the average avalanche size with duration), the shift of waiting time between the avalanches, and finally the collapses of the waiting time distributions. The results are obtained via extensive numerical simulations on the equilateral three-dimensional cubic lattices of various sizes and disorders.

Mean-field theories for depinning and their experimental signatures

Mean-field theory is an approximation replacing an extended system by a few variables. For depinning of elastic manifolds, these are the position u of its center of mass and the statistics of the forces F(u). There are two proposals how to model the latter: as a random walk (ABBM model), or as uncorrelated forces at integer u (discretized particle model, DPM). While for many experiments the ABBM model (in the literature misleadingly equated with mean-field theory) makes quantitatively correct predictions for the distributions of velocities, or avalanche size and duration, the microscopic disorder force-force correlations cannot grow linearly, and thus unboundedly as a random walk, with distance. Even the effective (renormalized) disorder forces which do so at small distances are bounded at large distances. To describe both regimes, we model forces as an Ornstein-Uhlenbeck process. The latter has the statistics of a random walk at small scales, and is uncorrelated at large scales. By connecting to results in both limits, we solve the model largely analytically, allowing us to describe in all regimes the distributions of velocity, avalanche size, and duration. To establish experimental signatures of this transition, we study the response function, and the correlation function of position u, velocity u[over ̇], and forces F under slow driving with velocity v>0. While at v=0 force or position correlations have a cusp at the origin and then decay at least exponentially fast to zero, this cusp is rounded at a finite driving velocity. We give a detailed analytic analysis for this rounding by velocity, which allows us, given experimental data, to extract the timescale of the response function, and to reconstruct the force-force correlator at v=0. The latter is the central object of the field theory, and as such contains detailed information about the universality class in question. We test our predictions by careful numerical simulations extending over up to ten orders in magnitude.

Local Probe Comparison of Ferroelectric Switching Event Statistics in the Creep and Depinning Regimes in Pb(Zr_{0.2}Ti_{0.8})O_{3} Thin Films

Ferroelectric materials provide a useful model system to explore the jerky, highly nonlinear dynamics of elastic interfaces in disordered media. The distribution of nanoscale switching event sizes is studied in two Pb(Zr_{0.2}Ti_{0.8})O_{3} thin films with different disorder landscapes using piezoresponse force microscopy. While the switching event statistics show the expected power-law scaling, significant variations in the value of the scaling exponent τ are seen, possibly as a consequence of the different intrinsic disorder landscapes in the samples and of further alterations under high tip bias applied during domain writing. Importantly, higher exponent values (1.98-2.87) are observed when crackling statistics are acquired only for events occurring in the creep regime. The exponents are systematically lowered when all events across both creep and depinning regimes are considered-the first time such a distinction is made in studies of ferroelectric materials. These results show that distinguishing the two regimes is of crucial importance, significantly affecting the exponent value and potentially leading to incorrect assignment of universality class.

Propagating bands of plastic deformation in a metal alloy as critical avalanches

The plastic deformation of metal alloys localizes in the Portevin-Le Chatelier effect in bands of different types, including propagating, or type "A" bands, usually characterized by their width and a typical propagation velocity. This plastic instability arises from collective dynamics of dislocations interacting with mobile solute atoms, but the resulting sensitivity to the strain rate lacks fundamental understanding. Here, we show, by using high-resolution imaging in tensile deformation experiments of an aluminum alloy, that the band velocities exhibit large fluctuations. Each band produces a velocity signal reminiscent of crackling noise bursts observed in numerous driven avalanching systems from propagating cracks in fracture to the Barkhausen effect in ferromagnets. The statistical features of these velocity bursts including their average shapes and size distributions obey predictions of a simple mean-field model of critical avalanche dynamics. Our results thus reveal a previously unknown paradigm of criticality in the localization of deformation.

Avalanche properties in striplike ferromagnetic systems

We present numerical findings on the behavior of the athermal nonequilibrium random-field Ising model of spins at the thin striplike L_{1}×L_{2}×L_{3} cubic lattices with L_{1}<L_{2}<L_{3}. Changing of system sizes highly influences the evolution and shape of avalanches. The smallest avalanches [classified as three-dimension- (3D) like] are unaffected by the system boundaries, the larger are sandwiched between the top and bottom system faces so are 2D-like, while the largest are extended over the system lateral cross section and propagate along the length L_{3} like in 1D systems. Such a structure of avalanches causes double power-law distributions of their size, duration, and energy with larger effective critical exponent corresponding to 3D-like and smaller to 2D-like avalanches. The distributions scale with thickness L_{1} and are collapsible following the proposed scaling predictions which, together with the distributions' shape, might be important for analysis of the Barkhausen noise experimental data for striplike samples. Finally, the impact of system size on external field that triggers the largest avalanche for a given disorder is presented and discussed.

Correlations between avalanches in the depinning dynamics of elastic interfaces

We study the correlations between avalanches in the depinning dynamics of elastic interfaces driven on a random substrate. In the mean-field theory (the Brownian force model), it is known that the avalanches are uncorrelated. Here we obtain a simple field theory which describes the first deviations from this uncorrelated behavior in a ε=d_{c}-d expansion below the upper critical dimension d_{c} of the model. We apply it to calculate the correlations between (i) avalanche sizes (ii) avalanche dynamics in two successive avalanches, or more generally, in two avalanches separated by a uniform displacement W of the interface. For (i) we obtain the correlations of the total sizes, of the local sizes, and of the total sizes with given seeds (starting points). For (ii) we obtain the correlations of the velocities, of the durations, and of the avalanche shapes. In general we find that the avalanches are anticorrelated, the occurrence of a larger avalanche making more likely the occurrence of a smaller one, and vice versa. Examining the universality of our results leads us to conjecture several exact scaling relations for the critical exponents that characterize the different distributions of correlations. The avalanche size predictions are confronted to numerical simulations for a d=1 interface with short range elasticity. They are also compared to our recent related work on static avalanches (shocks). Finally we show that the naive extrapolation of our result into the thermally activated creep regime at finite temperature predicts strong positive correlations between the forward motion events, as recently observed in numerical simulations.

Geometry-induced nonequilibrium phase transition in sandpiles

We study the sandpile model on three-dimensional spanning Ising clusters with the temperature T treated as the control parameter. By analyzing the three-dimensional avalanches and their two-dimensional projections (which show scale-invariant behavior for all temperatures), we uncover two universality classes with different exponents (an ordinary BTW class, and SOC_{T=∞}), along with a tricritical point (at T_{c}, the critical temperature of the host) between them. The transition between these two criticalities is induced by the transition in the support. The SOC_{T=∞} universality class is characterized by the exponent of the avalanche size distribution τ^{T=∞}=1.27±0.03, consistent with the exponent of the size distribution of the Barkhausen avalanches in amorphous ferromagnets Durin and Zapperi [Phys. Rev. Lett. 84, 4705 (2000)PRLTAO0031-900710.1103/PhysRevLett.84.4705]. The tricritical point is characterized by its own critical exponents. In addition to the avalanche exponents, some other quantities like the average height, the spanning avalanche probability (SAP), and the average coordination number of the Ising clusters change significantly the behavior at this point, and also exhibit power-law behavior in terms of ε≡T-T_{c}/T_{c}, defining further critical exponents. Importantly, the finite-size analysis for the activity (number of topplings) per site shows the scaling behavior with exponents β=0.19±0.02 and ν=0.75±0.05. A similar behavior is also seen for the SAP and the average avalanche height. The fractal dimension of the external perimeter of the two-dimensional projections of avalanches is shown to be robust against T with the numerical value D_{f}=1.25±0.01.

Seismiclike organization of avalanches in a driven long-range elastic string as a paradigm of brittle cracks

Crack growth in heterogeneous materials sometimes exhibits crackling dynamics, made of successive impulselike events with specific scale-invariant time and size organization reminiscent of earthquakes. Here, we examine this dynamics in a model which identifies the crack front with a long-range elastic line driven in a random potential. We demonstrate that, under some circumstances, fracture grows intermittently, via scale-free impulse organized into aftershock sequences obeying the fundamental laws of statistical seismology. We examine the effects of the driving rate and system overall stiffness (unloading factor) onto the scaling exponents and cutoffs associated with the time and size organization. We unravel the specific conditions required to observe a seismiclike organization in the crack propagation problem. Beyond failure problems, implications of these results to other crackling systems are finally discussed.

Avalanche statistics during coarsening dynamics

We study the coarsening dynamics of a two-dimensional system via numerical simulations. The system under consideration is a biphasic system consisting of domains of a dispersed phase closely packed together in a continuous phase and separated by thin interfaces. Such a system is elastic and typically out of equilibrium. The equilibrium state is attained via the coarsening dynamics, wherein the dispersed phase slowly diffuses through the interfaces, causing the domains to change in size and eventually rearrange abruptly. The effect of rearrangements is propagated throughout the system via the intrinsic elastic interactions and may cause rearrangements elsewhere, resulting in intermittent bursts of activity and avalanche behaviour. Here we aim at quantitatively characterizing the corresponding avalanche statistics (i.e. size, duration, and inter-avalanche time). Despite the coarsening dynamics is triggered by an internal driving mechanism, we find quantitative indications that such avalanche statistics displays scaling-laws very similar to those observed in the response of disordered materials to external loads.

Playing with universality classes of Barkhausen avalanches

Many systems crackle, from earthquakes and financial markets to Barkhausen effect in ferromagnetic materials. Despite the diversity in essence, the noise emitted in these dynamical systems consists of avalanche-like events with broad range of sizes and durations, characterized by power-law avalanche distributions and typical average avalanche shape that are fingerprints describing the universality class of the underlying avalanche dynamics. Here we focus on the crackling noise in ferromagnets and scrutinize the traditional statistics of Barkhausen avalanches in polycrystalline and amorphous ferromagnetic films having different thicknesses. We show how scaling exponents and average shape of the avalanches evolve with the structural character of the materials and film thickness. We find quantitative agreement between experiment and theoretical predictions of models for the magnetic domain wall dynamics, and then elucidate the universality classes of Barkhausen avalanches in ferromagnetic films. Thereby, we observe for the first time the dimensional crossover in the domain wall dynamics and the outcomes of the interplay between system dimensionality and range of interactions governing the domain wall dynamics on Barkhausen avalanches.

Gaps between avalanches in one-dimensional random-field Ising models

We analyze the statistics of gaps (ΔH) between successive avalanches in one-dimensional random-field Ising models (RFIMs) in an external field H at zero temperature. In the first part of the paper we study the nearest-neighbor ferromagnetic RFIM. We map the sequence of avalanches in this system to a nonhomogeneous Poisson process with an H-dependent rate ρ(H). We use this to analytically compute the distribution of gaps P(ΔH) between avalanches as the field is increased monotonically from -∞ to +∞. We show that P(ΔH) tends to a constant C(R) as ΔH→0^{+}, which displays a nontrivial behavior with the strength of disorder R. We verify our predictions with numerical simulations. In the second part of the paper, motivated by avalanche gap distributions in driven disordered amorphous solids, we study a long-range antiferromagnetic RFIM. This model displays a gapped behavior P(ΔH)=0 up to a system size dependent offset value ΔH_{off}, and P(ΔH)∼(ΔH-ΔH_{off})^{θ} as ΔH→H_{off}^{+}. We perform numerical simulations on this model and determine θ≈0.95(5). We also discuss mechanisms which would lead to a nonzero exponent θ for general spin models with quenched random fields.

Spatial shape of avalanches

In disordered elastic systems, driven by displacing a parabolic confining potential adiabatically slowly, all advance of the system is in bursts, termed avalanches. Avalanches have a finite extension in time, which is much smaller than the waiting time between them. Avalanches also have a finite extension ℓ in space, i.e., only a part of the interface of size ℓ moves during an avalanche. Here we study their spatial shape 〈S(x)〉_{ℓ} given ℓ, as well as its fluctuations encoded in the second cumulant 〈S^{2}(x)〉_{ℓ}^{c}. We establish scaling relations governing the behavior close to the boundary. We then give analytic results for the Brownian force model, in which the microscopic disorder for each degree of freedom is a random walk. Finally, we confirm these results with numerical simulations. To do this properly we elucidate the influence of discretization effects, which also confirms the assumptions entering into the scaling ansatz. This allows us to reach the scaling limit already for avalanches of moderate size. We find excellent agreement for the universal shape and its fluctuations, including all amplitudes.

Universal temporal characteristics and vanishing of multifractality in Barkhausen avalanches

Barkhausen effect in ferromagnetic materials provides an excellent area for investigating scaling phenomena found in disordered systems exhibiting crackling noise. The critical dynamics is characterized by random pulses or avalanches with scale-invariant properties, power-law distributions, and universal features. However, the traditional Barkhausen avalanches statistics may not be sufficient to fully characterize the complex temporal correlation of the magnetic domain walls dynamics. Here we focus on the multifractal scenario to quantify the temporal scaling characteristics of Barkhausen avalanches in polycrystalline and amorphous ferromagnetic films with thicknesses from 50 to 1000 nm. We show that the multifractal properties are dependent on film thickness, although they seem to be insensitive to the structural character of the materials. Moreover, we observe for the first time the vanishing of the multifractality in the domain walls dynamics. As the thickness is reduced, the multifractal behavior gives place to a monofractal one over the entire range of time scales. This reorganization in the temporal scaling characteristics of Barkhausen avalanches is understood as a universal restructuring associated to the dimensional crossover, from three- to two-dimensional magnetization dynamics.

Spatiotemporal Patterns in Ultraslow Domain Wall Creep Dynamics

In the presence of impurities, ferromagnetic and ferroelectric domain walls slide only above a finite external field. Close to this depinning threshold, they proceed by large and abrupt jumps called avalanches, while, at much smaller fields, these interfaces creep by thermal activation. In this Letter, we develop a novel numerical technique that captures the ultraslow creep regime over huge time scales. We point out the existence of activated events that involve collective reorganizations similar to avalanches, but, at variance with them, display correlated spatiotemporal patterns that resemble the complex sequence of aftershocks observed after a large earthquake. Remarkably, we show that events assemble in independent clusters that display at large scales the same statistics as critical depinning avalanches. We foresee these correlated dynamics being experimentally accessible by magnetooptical imaging of ferromagnetic films.

Interevent Correlations from Avalanches Hiding Below the Detection Threshold

Numerous systems ranging from deformation of materials to earthquakes exhibit bursty dynamics, which consist of a sequence of events with a broad event size distribution. Very often these events are observed to be temporally correlated or clustered, evidenced by power-law-distributed waiting times separating two consecutive activity bursts. We show how such interevent correlations arise simply because of a finite detection threshold, created by the limited sensitivity of the measurement apparatus, or used to subtract background activity or noise from the activity signal. Data from crack-propagation experiments and numerical simulations of a nonequilibrium crack-line model demonstrate how thresholding leads to correlated bursts of activity by separating the avalanche events into subavalanches. The resulting temporal subavalanche correlations are well described by our general scaling description of thresholding-induced correlations in crackling noise.

Bursts of activity in collective cell migration

Dense monolayers of living cells display intriguing relaxation dynamics, reminiscent of soft and glassy materials close to the jamming transition, and migrate collectively when space is available, as in wound healing or in cancer invasion. Here we show that collective cell migration occurs in bursts that are similar to those recorded in the propagation of cracks, fluid fronts in porous media, and ferromagnetic domain walls. In analogy with these systems, the distribution of activity bursts displays scaling laws that are universal in different cell types and for cells moving on different substrates. The main features of the invasion dynamics are quantitatively captured by a model of interacting active particles moving in a disordered landscape. Our results illustrate that collective motion of living cells is analogous to the corresponding dynamics in driven, but inanimate, systems.

Universal correlations between shocks in the ground state of elastic interfaces in disordered media

The ground state of an elastic interface in a disordered medium undergoes collective jumps upon variation of external parameters. These mesoscopic jumps are called shocks, or static avalanches. Submitting the interface to a parabolic potential centered at w, we study the avalanches which occur as w is varied. We are interested in the correlations between the avalanche sizes S_{1} and S_{2} occurring at positions w_{1} and w_{2}. Using the functional renormalization group (FRG), we show that correlations exist for realistic interface models below their upper critical dimension. Notably, the connected moment 〈S_{1}S_{2}〉^{c} is up to a prefactor exactly the renormalized disorder correlator, itself a function of |w_{2}-w_{1}|. The latter is the universal function at the center of the FRG; hence, correlations between shocks are universal as well. All moments and the full joint probability distribution are computed to first nontrivial order in an ε expansion below the upper critical dimension. To quantify the local nature of the coupling between avalanches, we calculate the correlations of their local jumps. We finally test our predictions against simulations of a particle in random-bond and random-force disorder, with surprisingly good agreement.

Quantitative Scaling of Magnetic Avalanches

We provide the first quantitative comparison between Barkhausen noise experiments and recent predictions from the theory of avalanches for pinned interfaces, both in and beyond mean field. We study different classes of soft magnetic materials, including polycrystals and amorphous samples-which are characterized by long-range and short-range elasticity, respectively-both for thick and thin samples, i.e., with and without eddy currents. The temporal avalanche shape at fixed size as well as observables related to the joint distribution of sizes and durations are analyzed in detail. Both long-range and short-range samples with no eddy currents are fitted extremely well by the theoretical predictions. In particular, the short-range samples provide the first reliable test of the theory beyond mean field. The thick samples show systematic deviations from the scaling theory, providing unambiguous signatures for the presence of eddy currents.

Driving Rate Dependence of Avalanche Statistics and Shapes at the Yielding Transition

We study stress time series caused by plastic avalanches in athermally sheared disordered materials. Using particle-based simulations and a mesoscopic elastoplastic model, we analyze system size and shear-rate dependence of the stress-drop duration and size distributions together with their average temporal shape. We find critical exponents different from mean-field predictions, and a clear asymmetry for individual avalanches. We probe scaling relations for the rate dependency of the dynamics and we report a crossover towards mean-field results for strong driving.

Crossover from Localized to Cascade Relaxations in Metallic Glasses

Thermally activated deformation is investigated in two metallic glass systems with different cooling histories. By probing the atomic displacements and stress changes on the potential energy landscape, two deformation modes, a localized process and cascade process, have observed. The localized deformation involves fewer than 30 atoms and appears in both systems, and its size is invariant with cooling history. However, the cascade deformation is more frequently observed in the fast quenched system than in the slowly quenched system. The origin of the cascade process in the fast quenched system is attributed to the higher density of local minima on the underlying potential energy landscape.

Slip statistics of dislocation avalanches under different loading modes

Slowly compressed microcrystals deform via intermittent slip events, observed as displacement jumps or stress drops. Experiments often use one of two loading modes: an increasing applied stress (stress driven, soft), or a constant strain rate (strain driven, hard). In this work we experimentally test the influence of the deformation loading conditions on the scaling behavior of slip events. It is found that these common deformation modes strongly affect time series properties, but not the scaling behavior of the slip statistics when analyzed with a mean-field model. With increasing plastic strain, the slip events are found to be smaller and more frequent when strain driven, and the slip-size distributions obtained for both drives collapse onto the same scaling function with the same exponents. The experimental results agree with the predictions of the used mean-field model, linking the slip behavior under different loading modes.

Stochastic formation of magnetic vortex structures in asymmetric disks triggered by chaotic dynamics

The non-trivial spin configuration in a magnetic vortex is a prototype for fundamental studies of nanoscale spin behaviour with potential applications in magnetic information technologies. Arrays of magnetic vortices interfacing with perpendicular thin films have recently been proposed as enabler for skyrmionic structures at room temperature, which has opened exciting perspectives on practical applications of skyrmions. An important milestone for achieving not only such skyrmion materials but also general applications of magnetic vortices is a reliable control of vortex structures. However, controlling magnetic processes is hampered by stochastic behaviour, which is associated with thermal fluctuations in general. Here we show that the dynamics in the initial stages of vortex formation on an ultrafast timescale plays a dominating role for the stochastic behaviour observed at steady state. Our results show that the intrinsic stochastic nature of vortex creation can be controlled by adjusting the interdisk distance in asymmetric disk arrays.

Statistical properties of Barkhausen noise in amorphous ferromagnetic films

We investigate the statistical properties of the Barkhausen noise in amorphous ferromagnetic films with thicknesses in the range between 100 and 1000 nm. From Barkhausen noise time series measured with the traditional inductive technique, we perform a wide statistical analysis and establish the scaling exponents τ,α,1/σνz, and ϑ. We also focus on the average shape of the avalanches, which gives further indications on the domain-wall dynamics. Based on experimental results, we group the amorphous films in a single universality class, characterized by scaling exponents τ=1.28±0.02,α=1.52±0.3, and 1/σνz=ϑ=1.83±0.03, values compatible with that obtained for several bulk amorphous magnetic materials. Besides, we verify that the avalanche shape depends on the universality class. By considering the theoretical models for the dynamics of a ferromagnetic domain wall driven by an external magnetic field through a disordered medium found in literature, we interpret the results and identify an experimental evidence that these amorphous films, within this thickness range, present a typical three-dimensional magnetic behavior with predominant short-range elastic interactions governing the domain-wall dynamics. Moreover, we provide experimental support for the validity of a general scaling form for the average avalanche shape for non-mean-field systems.

Unzipping bird feathers

The bird feather vane can be separated into two parts by pulling the barbs apart. The original state can be re-established easily by lightly stroking through the feather. Hooklets responsible for holding vane barbs together are not damaged by multiple zipping and unzipping cycles. Because numerous microhooks keep the integrity of the feather, their properties are of great interest for understanding mechanics of the entire feather structure. This study was undertaken to estimate the separation force of single hooklets and their arrays using force measurement of an unzipping feather vane. The hooklets usually separate in some number synchronously (20 on average) with the highest observed separation force of 1.74 mN (average force 0.27 mN), whereas the single hooklet separation force was 14 μN. A simple numerical model was suggested for a better understanding of zipping and unzipping behaviour in feathers. The model demonstrates features similar to those observed in experiments.

Role of superposition of dislocation avalanches in the statistics of acoustic emission during plastic deformation

Various dynamical systems with many degrees of freedom display avalanche dynamics, which is characterized by scale invariance reflected in power-law statistics. The superposition of avalanche processes in real systems driven at a finite velocity may influence the experimental determination of the underlying power law. The present paper reports results of an investigation of this effect using the example of acoustic emission (AE) accompanying plastic deformation of crystals. Indeed, recent studies of AE did not only prove that the dynamics of crystal defects obeys power-law statistics, but also led to a hypothesis of universality of the scaling law. We examine the sensitivity of the apparent statistics of AE to the parameters applied to individualize AE events. Two different alloys, MgZr and AlMg, both displaying strong AE but characterized by different plasticity mechanisms, are investigated. It is shown that the power-law indices display a good robustness in wide ranges of parameters even in the conditions leading to very strong superposition of AE events, although some deviations from the persistent values are also detected. The totality of the results confirms the scale-invariant character of deformation processes on the scale relevant to AE, but uncovers essential differences between the power-law exponents found for two kinds of alloys.

Universal properties of magnetization dynamics in polycrystalline ferromagnetic films

We investigate the scaling behavior in the statistical properties of Barkhausen noise in ferromagnetic films. We apply the statistical treatment usually employed for bulk materials in experimental Barkhausen noise time series measured with the traditional inductive technique in polycrystalline ferromagnetic films having different thickness from 100 to 1000 nm and determine the scaling exponents. Based on this procedure, we group the samples in a single universality class, since the scaling behavior of Barkhausen avalanches is characterized by exponents τ∼1.5, α∼2.0, and 1/σνz∼ϑ∼2.0 for all the films. We interpret these results in terms of theoretical models and provide experimental evidence that a well-known mean-field model for the dynamics of a ferromagnetic domain wall in three-dimensional ferromagnets can be extended for films. We identify that the films present an universal three-dimensional magnetization dynamics, governed by long-range dipolar interactions, even at the smallest thicknesses, indicating that the two-dimensional magnetic behavior commonly verified for films cannot be generalized for all thickness ranges.

Statistics of avalanches with relaxation and Barkhausen noise: a solvable model

We study a generalization of the Alessandro-Beatrice-Bertotti-Montorsi (ABBM) model of a particle in a Brownian force landscape, including retardation effects. We show that under monotonous driving the particle moves forward at all times, as it does in absence of retardation (Middleton's theorem). This remarkable property allows us to develop an analytical treatment. The model with an exponentially decaying memory kernel is realized in Barkhausen experiments with eddy-current relaxation and has previously been shown numerically to account for the experimentally observed asymmetry of Barkhausen pulse shapes. We elucidate another qualitatively new feature: the breakup of each avalanche of the standard ABBM model into a cluster of subavalanches, sharply delimited for slow relaxation under quasistatic driving. These conditions are typical for earthquake dynamics. With relaxation and aftershock clustering, the present model includes important ingredients for an effective description of earthquakes. We analyze quantitatively the limits of slow and fast relaxation for stationary driving with velocity v>0. The v-dependent power-law exponent for small velocities, and the critical driving velocity at which the particle velocity never vanishes, are modified. We also analyze nonstationary avalanches following a step in the driving magnetic field. Analytically, we obtain the mean avalanche shape at fixed size, the duration distribution of the first subavalanche, and the time dependence of the mean velocity. We propose to study these observables in experiments, allowing a direct measurement of the shape of the memory kernel and tracing eddy current relaxation in Barkhausen noise.

Crackling versus continuumlike dynamics in brittle failure

We study how the loading rate, specimen geometry, and microstructural texture select the dynamics of a crack moving through an heterogeneous elastic material in the quasistatic approximation. We find a transition, fully controlled by two dimensionless variables, between dynamics ruled by continuum fracture mechanics and crackling dynamics. Selection of the latter by the loading, microstructure, and specimen parameters is formulated in terms of scaling laws on the power spectrum of crack velocity. This analysis defines the experimental conditions required to observe crackling in fracture. Beyond failure problems, the results extend to a variety of situations described by models of the same universality class, e.g., the dynamics in wetting or of domain walls in amorphous ferromagnets.

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.

Mechanisms of evolution of avalanches in regular graphs

A mapping of avalanches occurring in the zero-temperature random-field Ising model to life periods of a population experiencing immigration is established. Such a mapping allows the microscopic criteria for the occurrence of an infinite avalanche in a q-regular graph to be determined. A key factor for an avalanche of spin flips to become infinite is that it interacts in an optimal way with previously flipped spins. Based on these criteria, we explain why an infinite avalanche can occur in q-regular graphs only for q>3 and suggest that this criterion might be relevant for other systems. The generating function techniques developed for branching processes are applied to obtain analytical expressions for the durations, pulse shapes, and power spectra of the avalanches. The results show that only very long avalanches exhibit a significant degree of universality.

Universal fluctuations and extreme statistics of avalanches near the depinning transition

We derive exact predictions for universal scaling exponents and scaling functions associated with the statistics of maximum velocities v(m) during avalanches described by the mean-field theory of the interface depinning transition. In particular, we find a robust power-law regime in the statistics of maximum events that can explain the observed distribution of the peak amplitudes in acoustic emission experiments of crystal plasticity. Our results are expected to be broadly applicable to a broad range of systems in the mean-field interface depinning universality class, ranging from magnets to earthquakes.

Evolution of the average avalanche shape with the universality class

A multitude of systems ranging from the Barkhausen effect in ferromagnetic materials to plastic deformation and earthquakes respond to slow external driving by exhibiting intermittent, scale-free avalanche dynamics or crackling noise. The avalanches are power-law distributed in size, and have a typical average shape: these are the two most important signatures of avalanching systems. Here we show how the average avalanche shape evolves with the universality class of the avalanche dynamics by employing a combination of scaling theory, extensive numerical simulations and data from crack propagation experiments. It follows a simple scaling form parameterized by two numbers, the scaling exponent relating the average avalanche size to its duration and a parameter characterizing the temporal asymmetry of the avalanches. The latter reflects a broken time-reversal symmetry in the avalanche dynamics, emerging from the local nature of the interaction kernel mediating the avalanche dynamics.

Multifractality in domain wall dynamics of a ferromagnetic film

We investigate the multifractal properties in the dynamics of domain walls of a ferromagnetic film. We apply the Multifractal Detrended Fluctuation Analysis method in experimental Barkhausen noise time series measured in a 1000-nm-thick Permalloy film under different driving magnetic field frequencies, and calculate the fluctuation function F_{q}(s), generalized Hurst exponent h(q), multifractal scaling exponent τ(q), and the multifractal spectrum f(α). Based on this procedure, we provide experimental evidence of multifractality in the dynamics of domain walls in ferromagnetic films and identify a rich and strong multifractal behavior, revealed by the changes of the scaling properties of over the entire Barkhausen noise signal, independently of the driving magnetic field rate employed in the experiment.

Distribution of velocities and acceleration for a particle in Brownian correlated disorder: inertial case

We study the motion of an elastic object driven in a disordered environment in presence of both dissipation and inertia. We consider random forces with the statistics of random walks and reduce the problem to a single degree of freedom. It is the extension of the mean-field Alessandro-Beatrice- Bertotti-Montorsi (ABBM) model in presence of an inertial mass m. While the ABBM model can be solved exactly, its extension to inertia exhibits complicated history dependence due to oscillations and backward motion. The characteristic scales for avalanche motion are studied from numerics and qualitative arguments. To make analytical progress, we consider two variants which coincide with the original model whenever the particle moves only forward. Using a combination of analytical and numerical methods together with simulations, we characterize the distributions of instantaneous acceleration and velocity, and compare them in these three models. We show that for large driving velocity, all three models share the same large-deviation function for positive velocities, which is obtained analytically for small and large m, as well as for m=6/25. The effect of small additional thermal and quantum fluctuations can be treated within an approximate method.

Nonstationary dynamics of the Alessandro-Beatrice-Bertotti-Montorsi model

We obtain an exact solution for the motion of a particle driven by a spring in a Brownian random-force landscape, the Alessandro-Beatrice-Bertotti-Montorsi (ABBM) model. Many experiments on quasistatic driving of elastic interfaces (Barkhausen noise in magnets, earthquake statistics, shear dynamics of granular matter) exhibit the same universal behavior as this model. It also appears as a limit in the field theory of elastic manifolds. Here we discuss predictions of the ABBM model for monotonous, but otherwise arbitrary, time-dependent driving. Our main result is an explicit formula for the generating functional of particle velocities and positions. We apply this to derive the particle-velocity distribution following a quench in the driving velocity. We also obtain the joint avalanche size and duration distribution and the mean avalanche shape following a jump in the position of the confining spring. Such nonstationary driving is easy to realize in experiments, and provides a way to test the ABBM model beyond the stationary, quasistatic regime. We study extensions to two elastically coupled layers, and to an elastic interface of internal dimension d, in the Brownian force landscape. The effective action of the field theory is equal to the action, up to one-loop corrections obtained exactly from a functional determinant. This provides a connection to renormalization-group methods.

Emergent rotational symmetries in disordered magnetic domain patterns

Uniaxial systems often form labyrinthine domains that exhibit short-range order but are macroscopically isotropic and would not be expected to exhibit precise symmetries. However, their underlying frustration results in a multitude of metastable configurations of comparable energy, and driving such a system externally might lead to pattern formation. We find that soft x-ray speckle diffraction patterns of the labyrinthine domains in CoPd/IrMn heterostructures reveal a diverse array of hidden rotational symmetries about the magnetization axis, thereby suggesting an unusual form of emergent order in an otherwise disordered system. These symmetries depend on applied magnetic field, magnetization history, and scattering wave vector. Maps of rotational symmetry exhibit intriguing structures that can be controlled by manipulating the applied magnetic field in concert with the exchange bias condition.

Tuning coupling: discrete changes in runaway avalanche sizes in disordered media

Hysteretic systems may exhibit a runaway avalanche in which a large fraction of the constituents of the system collectively change state. It would be very valuable to understand the role that interaction strength between constituents plays in the size of such catastrophic runaway avalanches. We use a simple model, the random field Ising model, to study how the size of the runaway avalanche changes as the coupling between spins, J, is tuned. In particular, we calculate P(S), the distribution of size changes S in the runaway avalanche size as J comes close to a critical value J(c), and find that the distribution scales as P(S)∼S(-τ)D(S(σ)(J/J(c)-1)), with τ and σ critical exponents and D(x) a universal scaling function. In mean field theory we find τ=3/2, σ=1/2, and D(x)=exp[-(3x)(1/σ)/2]. On the basis of these results and previous studies, we also predict that for three dimensions τ=1.6 and σ=0.24.

Universal hysteresis scaling for incommensurate magnetic order in dysprosium

We study the scaling properties of magnetic minor hysteresis loops in a polycrystalline dysprosium metal, varying temperature and magnetic-field amplitudes. We observe irreversibility-related hysteresis loss in the helical antiferromagnetic phase, which is related with remanent flux density as a power law with the same scaling exponent of 1.25±0.05 as that in ferromagnetic materials. In contrast to hysteresis scalings in ferromagnets associated with 180° Bloch walls, the observed law is governed by spiral walls which separate helical domains with oppositely rotating spins.

Tracking random walk of individual domain walls in cylindrical nanomagnets with resistance noise

The stochasticity of domain-wall (DW) motion in magnetic nanowires has been probed by measuring slow fluctuations, or noise, in electrical resistance at small magnetic fields. By controlled injection of DWs into isolated cylindrical nanowires of nickel, we have been able to track the motion of the DWs between the electrical leads by discrete steps in the resistance. Closer inspection of the time dependence of noise reveals a diffusive random walk of the DWs with a universal kinetic exponent. Our experiments outline a method with which electrical resistance is able to detect the kinetic state of the DWs inside the nanowires, which can be useful in DW-based memory designs.

Elasticity of a contact-line and avalanche-size distribution at depinning

Motivated by recent experiments, we extend the Joanny-deGennes calculation of the elasticity of a contact line to an arbitrary contact angle and an arbitrary plate inclination in presence of gravity. This requires a diagonalization of the elastic modes around the nonlinear equilibrium profile, which is carried out exactly. We then make detailed predictions for the avalanche-size distribution at quasistatic depinning: we study how the universal (i.e., short-scale independent) rescaled size distribution and the ratio of moments of local to global avalanches depend on the precise form of the elastic kernel.

Phenomenological approach to mechanical damage growth analysis

The problem of characterizing damage evolution in a generic material is addressed with the aim of tracing it back to existing growth models in other fields of research. Based on energetic considerations, a system evolution equation is derived for a generic damage indicator describing a material system subjected to an increasing external stress. The latter is found to fit into the framework of a recently developed phenomenological universality (PUN) approach and, more specifically, the so-called U2 class. Analytical results are confirmed by numerical simulations based on a fiber-bundle model and statistically assigned local strengths at the microscale. The fits with numerical data prove, with an excellent degree of reliability, that the typical evolution of the damage indicator belongs to the aforementioned PUN class. Applications of this result are briefly discussed and suggested.