Crackling dynamics in material failure as the signature of a self-organized dynamic phase transition

Crack Propagation by Activated Avalanches during Creep and Fatigue from Elastic Interface Theory

The growth of cracks combines materials science, fracture mechanics, and statistical physics. The importance of fluctuations in the crack velocity is fundamental since it signals that the crack overcomes local barriers such as tough spots by avalanches. In ductile materials the omnipresent plasticity close to the crack tip influences the growth by history effects, which we here study in polymethylmetacrylate by various fatigue and creep protocols. We show how the crack tip local history may be encompassed in a time- and protocol-dependent length scale, which allows us to apply a statistical fracture description to the time-dependent crack growth rate, resolving the well-known paradox why fatigue cracks grow faster if the stress during a cycle is let to relax more from the peak value. The results open up novel directions for understanding fracture by statistical physics.

Fracture process of composite materials in a spring network model

We analyze a two-dimensional spring network model comprising breakable and unbreakable springs. Computer simulations showed this system to exhibit intermittent stress drops in a larger strain regime, and these stress drops resulted in ductilelike behavior. The scaling analysis reveals that the avalanche size distribution demonstrates a cutoff, depending on its internal structure. This study also investigates the relationship between cluster growth and stress drop, and we show that the amount of stress drop increases in terms of power law, corresponding to crack growth. The crack length distribution also demonstrates a cutoff depending on its internal structure. The results show that both the cluster growth-stress drop relationship and the crack size distribution are scaled by the quantity related to the internal structure, and the relevance of the exponent that scales the cluster growth-stress drop relationship to the exponent that scales crack size distribution is verified.

Verification of scaling behavior near dynamic phase transitions for nonantisymmetric field sequences

We investigate the scaling behavior of the magnetic dynamic order parameter Q in the vicinity of the dynamic phase transition (DPT) in the presence of temporal field sequences H(t) that are periodic with period P but lack half-wave antisymmetry. We verify by means of mean-field calculations that the scaling of Q is preserved in the vicinity of the second-order phase transition if one defines a suitable generalized conjugate field H^{*} that reestablishes the proper time-reversal symmetry. For the purpose of our quantitative data analysis, we employ the dynamic equivalent of the Arrott-Noakes equation of state, which allows for a simultaneous scaling analysis of the period P and the conjugate-field H^{*} dependence of Q. By doing so, we demonstrate that both the scaling behavior and universality are preserved, even if the dynamics is driven by a more general applied field sequence that lacks antisymmetry.

Anomalous phase diagram of the elastic interface with nonlocal hydrodynamic interactions in the presence of quenched disorder

We investigate the influence of quenched disorder on the steady states of driven systems of the elastic interface with nonlocal hydrodynamic interactions. The generalized elastic model (GEM), which has been used to characterize numerous physical systems such as polymers, membranes, single-file systems, rough interfaces, and fluctuating surfaces, is a standard approach to studying the dynamics of elastic interfaces with nonlocal hydrodynamic interactions. The criticality and phase transition of the quenched generalized elastic model are investigated numerically and the results are presented in a phase diagram spanned by two tuning parameters. We demonstrate that in the one-dimensional disordered driven GEM, three qualitatively different behavior regimes are possible with a proper specification of the order parameter (mean velocity) for this system. In the vanishing order parameter regime, the steady-state order parameter approaches zero in the thermodynamic limit. A system with a nonzero mean velocity can be in either the continuous regime, which is characterized by a second-order phase transition, or the discontinuous regime, which is characterized by a first-order phase transition. The focus of this research is to investigate the critical scaling features near the pinning-depinning threshold. The behavior of the quenched generalized elastic model at the critical depinning force is explored. Near the depinning threshold, the critical exponent is obtained numerically.

Analysis of the heat transfer fluctuations in the Rayleigh-Bénard convection of concentrated emulsions with finite-size droplets

Employing numerical simulations, we provide an accurate insight into the heat transfer mechanism in the Rayleigh-Bénard convection of concentrated emulsions with finite-size droplets. We focus on the unsteady dynamics characterizing the thermal convection of these complex fluids close to the transition from conductive to convective states, where the heat transfer phenomenon, expressed in terms of the Nusselt number Nu, is characterized by pronounced fluctuations triggered by collective droplet motion [F. Pelusi et al., Soft Matter, 2021, 17(13), 3709-3721]. By systematically increasing the droplet concentration, we show how these fluctuations emerge along with the segregation of "extreme events" in the boundary layers, causing intermittent bursts in the heat flux fluctuations. Furthermore, we quantify the extension S and the duration of the coherent droplet motion accompanying these extreme events via a suitable statistical analysis involving the droplet displacements. We show how the increase in droplet concentration results in a power-law behaviour of the probability distribution function of S and and how this outcome is robust at changing the analysis protocol. Our work offers a comprehensive picture, linking macroscopic heat transfer fluctuations with the statistics of droplets at the mesoscale.

Experimental Observation of Critical Scaling in Magnetic Dynamic Phase Transitions

We explore the critical behavior of dynamic phase transitions in ultrathin uniaxial Co films. Our data demonstrate the occurrence of critical fluctuations, which define the critical regime, and in which we conduct a scaling analysis of the dynamic order parameter Q, utilizing a dynamic analog to the Arrott-Noakes equation of state. Our results show dynamic critical exponents that agree with the 2D Ising model as theoretically predicted. However, equilibrium critical exponents of our sample agree with the 3D Ising model. We argue that these differences between dynamic and thermodynamic behavior are due to fundamentally different length scales at which dimensional crossovers occur.

Asymmetric Roughness of Elastic Interfaces at the Depinning Threshold

Roughness of driven elastic interfaces in random media is typically understood to be characterized by a single roughness exponent ζ. We show that at the depinning threshold, due to symmetry breaking caused by the direction of the driving force, elastic interfaces with local, long-range, and mean-field elasticity exhibit asymmetric roughness. It is manifested as a skewed distribution of the local interface heights, and can be quantified by using detrended fluctuation analysis to compute a spectrum of local, segment-level scaling exponents. The asymmetry is observed as approximately linear dependence of the local scaling exponents on the difference of the segment height from the mean interface height.

Clusters in an Epidemic Model with Long-Range Dispersal

In the presence of long-range dispersal, epidemics spread in spatially disconnected regions known as clusters. Here, we characterize exactly their statistical properties in a solvable model, in both the supercritical (outbreak) and critical regimes. We identify two diverging length scales, corresponding to the bulk and the outskirt of the epidemic. We reveal a nontrivial critical exponent that governs the cluster number and the distribution of their sizes and of the distances between them. We also discuss applications to depinning avalanches with long-range elasticity.

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.

Preventing a global transition to thermoacoustic instability by targeting local dynamics

The burning of fossil fuels to generate power produces harmful emissions. Lowering such emissions in gas turbine engines is possible by operating them at fuel-lean conditions. However, such strategies often fail because, under fuel-lean conditions, the combustors are prone to catastrophic high-amplitude oscillations known as thermoacoustic instability. We reveal that, as an operating parameter is varied in time, the transition to thermoacoustic instability is initiated at specific spatial regions before it is observed in larger regions of the combustor. We use two indicators to discover such inceptive regions: the growth of variance of fluctuations in spatially resolved heat release rate and its spatiotemporal evolution. In this study, we report experimental evidence of suppression of the global transition to thermoacoustic instability through targeted modification of local dynamics at the inceptive regions. We strategically arrange slots on the flame anchor, which, in turn, reduce the local heat release rate fluctuations at the inceptive regions and thus suppress the global transition to thermoacoustic instability. Our results open new perspectives for combustors that are more environmental-friendly.

Effect of thresholding on avalanches and their clustering for interfaces with long-range elasticity

Avalanches are often defined as signals higher than some detection level in bursty systems. The choice of the detection threshold affects the number of avalanches, but it can also affect their temporal correlations. We simulated the depinning of a long-range elastic interface and applied different thresholds including a zero one on the data to see how the sizes and durations of events change and how this affects temporal avalanche clustering. Higher thresholds result in steeper size and duration distributions and cause the avalanches to cluster temporally. Using methods from seismology, the frequency of the events in the clusters was found to decrease as a power-law of time, and the size of an event in a cluster was found to help predict how many events it is followed by. The results bring closer theoretical studies of this class of models to real experiments, but also highlight how different phenomena can be obtained from the same set of data.

Energy exponents of avalanches and Hausdorff dimensions of collapse patterns

A simple numerical model to simulate athermal avalanches is presented. The model is inspired by the "porous collapse" process where the compression of porous materials generates collapse cascades, leading to power law distributed avalanches. The energy (E), amplitude (A_{max}), and size (S) exponents are derived by computer simulation in two approximations. Time-dependent "jerk" spectra are calculated in a single avalanche model where each avalanche is simulated separately from other avalanches. The average avalanche profile is parabolic, the scaling between energy and amplitude follows E∼A_{max}^{2}, and the energy exponent is ε = 1.33. Adding a general noise term in a continuous event model generates infinite avalanche sequences which allow the evaluation of waiting time distributions and pattern formation. We find the validity of the Omori law and the same exponents as in the single avalanche model. We then add spatial correlations by stipulating the ratio G/N between growth processes G (linked to a previous event location) and nucleation processes N (with new, randomly chosen nucleation sites). We found, in good approximation, a power law correlation between the energy exponent ε and the Hausdorff dimension H_{D} of the resulting collapse pattern H_{D}-1∼ɛ^{-3}. The evolving patterns depend strongly on G/N with the distribution of collapse sites equally power law distributed. Its exponent ɛ_{topo} would be linked to the dynamical exponent ε if each collapse carried an energy equivalent to the size of the collapse. A complex correlation between ɛ,ɛ_{topo}, and H_{D} emerges, depending strongly on the relative occupancy of the collapse sites in the simulation box.

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.

General existence and determination of conjugate fields in dynamically ordered magnetic systems

We investigate experimentally as well as theoretically the dynamic magnetic phase diagram and its associated order parameter Q upon the application of a non-antisymmetric magnetic field sequence composed of a fundamental harmonic component H_{0}, a constant bias field H_{b}, and a second-harmonic component H_{2}. The broken time antisymmetry introduced by the second-harmonic field component H_{2} leads to an effective bias effect that is superimposed onto the influence of the static bias H_{b}. Despite this interference, we can demonstrate the existence of a generalized conjugate field H^{*} for the dynamic order parameter Q, to which both the static bias field and the second-harmonic Fourier amplitude of the field sequence contribute. Hereby, we observed that especially the conventional paramagnetic dynamic phase is very susceptible to the impact of the second-harmonic field component H_{2}, whereas this additional field component leads to only very minor phase-space modifications in the ferromagnetic and anomalous paramagnetic regions. In contrast to prior studies, we also observe that the critical point of the phase transition is shifted upon introducing a second-harmonic field component H_{2}, illustrating that the overall dynamic behavior of such magnetic systems is being driven by the total effective amplitude of the field sequence.

Controlling crackling dynamics by triggering low-intensity avalanches

We examine the effect of small, spatially localized excitations applied periodically in different manners, on the crackling dynamics of a brittle crack driven slowly in a heterogeneous solid. When properly adjusted, these excitations are observed to radically modify avalanche statistics and considerably limit the magnitude of the largest events. Surprisingly, this does not require information on the front loading state at the time of excitation; applying it either at a random location or at the most loaded point gives the same results. Subsequently, we unravel how the excitation amplitude, spatial extent, and frequency govern the effect. We find that the excitation efficiency is ruled by a single reduced parameter, namely the injected power per unit front length; the suppression of extreme avalanches is maximum at a well-defined optimal value of this control parameter. analysis opens another way to control the largest events in crackling dynamics. Beyond fracture problems, it may be relevant for crackling systems described by models of the same universality class, such as the wetting of heterogeneous substrates or magnetic walls in amorphous magnets.

Dynamic crack growth along heterogeneous planar interfaces: Interaction with unidimensional strips

We examine theoretically and numerically fast propagation of a tensile crack along unidimensional strips with periodically evolving toughness. In such dynamic fracture regimes, crack front waves form and transport front disturbances along the crack edge at speed less than the Rayleigh wave speed and depending on the crack speed. In this configuration, standing front waves dictate the spatiotemporal evolution of the local crack front speed, which takes a specific scaling form. Analytical examination of both the short-time and long-time limits of the problem reveals the parameter dependency with strip wavelength, toughness contrast and overall fracture speed. Implications and generalization to unidimensional strips of arbitrary shape are lastly discussed.

Spatial Clustering of Depinning Avalanches in Presence of Long-Range Interactions

Disordered elastic interfaces display avalanche dynamics at the depinning transition. For short-range interactions, avalanches correspond to compact reorganizations of the interface well described by the depinning theory. For long-range elasticity, an avalanche is a collection of spatially disconnected clusters. In this Letter we determine the scaling properties of the clusters and relate them to the roughness exponent of the interface. The key observation of our analysis is the identification of a Bienaymé-Galton-Watson process describing the statistics of the number of clusters. Our work has concrete importance for experimental applications where the cluster statistics is a key probe of avalanche dynamics.

Thermally rounded depinning of an elastic interface on a washboard potential

The thermal rounding of the depinning transition of an elastic interface sliding on a washboard potential is studied through analytic arguments and very accurate numerical simulations. We confirm the standard view that well below the depinning threshold the average velocity can be calculated considering thermally activated nucleation of defects. However, we find that the straightforward extension of this analysis to near or above the depinning threshold does not fully describe the physics of the thermally assisted motion. In particular, we find that exactly at the depinning point the average velocity does not follow a pure power law of the temperature as naively expected by the analogy with standard phase transitions but presents subtle logarithmic corrections. We explain the physical mechanisms behind these corrections and argue that they are nonpeculiar collective effects which may also apply to the case of interfaces sliding on uncorrelated disordered landscapes.

Impact-induced transition from damage to perforation

We investigate the impact-induced damage and fracture of a bar-shaped specimen of heterogeneous materials focusing on how the system approaches perforation as the impact energy is gradually increased. A simple model is constructed which represents the bar as two rigid blocks coupled by a breakable interface with disordered local strength. The bar is clamped at the two ends, and the fracture process is initiated by an impactor hitting the bar in the middle. Our calculations revealed that depending on the imparted energy, the system has two phases: at low impact energies the bar suffers damage but keeps its integrity, while at sufficiently high energies, complete perforation occurs. We demonstrate that the transition from damage to perforation occurs analogous to continuous phase transitions. Approaching the critical point from below, the intact fraction of the interface goes to zero, while the deformation rate of the bar diverges according to power laws as function of the distance from the critical energy. As the degree of disorder increases, farther from the transition point the critical exponents agree with their zero disorder counterparts; however, close to the critical point a crossover occurs to a higher exponent.

Condensation Induced Blistering as a Measurement Technique for the Adhesion Energy of Nanoscale Polymer Films

Polymeric coatings having micro-to-nanoscale thickness show immense promise for enhancing thermal transport, catalysis, energy conversion, and water collection. Characterizing the work of adhesion (G) between these coatings and their substrates is key to understanding transport physics as well as mechanical reliability. Here, we demonstrate that water vapor condensation blistering is capable of in situ measurement of work of adhesion at the interface of polymer thin films with micrometer spatial resolution. We use our method to characterize adhesion of interfaces with controlled chemistry such as fluorocarbon/fluorocarbon (CF_n/CF_m, n, m = 0-3), fluorocarbon/hydrocarbon (CF_n/CH_m), fluorocarbon/silica (CF_n/SiO₂), and hydrocarbon/silica (CH_n/SiO₂) interfaces showing excellent agreement with adhesion energy measured by the contact angle approach. We demonstrate the capability of our condensation blister test to achieve measurement spatial resolutions as low as 10 μm with uncertainties of ∼10%. The outcomes of this work establish a simple tool to study interfacial adhesion.

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.

Universal Scaling of the Velocity Field in Crack Front Propagation

The propagation of a crack front in disordered materials is jerky and characterized by bursts of activity, called avalanches. These phenomena are the manifestation of an out-of-equilibrium phase transition originated by the disorder. As a result avalanches display universal scalings which are, however, difficult to characterize in experiments at a finite drive. Here, we show that the correlation functions of the velocity field along the front allow us to extract the critical exponents of the transition and to identify the universality class of the system. We employ these correlations to characterize the universal behavior of the transition in simulations and in an experiment of crack propagation. This analysis is robust, efficient, and can be extended to all systems displaying avalanche dynamics.

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.

Seismicity in sheared granular matter

We report on experiments investigating the dynamics of a slider that is pulled by a spring across a granular medium consisting of a vertical layer of photoelastic disks. The motion proceeds through a sequence of discrete events, analogous to seismic shocks, in which elastic energy stored in the spring is rapidly released. We measure the statistics of several properties of the individual events: the energy loss in the spring, the duration of the movement, and the temporal profile of the slider motion. We also study certain conditional probabilities and the statistics of mainshock-aftershock sequences. At low driving rates, we observe crackling with Omori-Utsu, Båth, and waiting time laws similar to those observed in seismic dynamics. At higher driving rates, where the sequence of events shows strong periodicity, we observe scaling laws and asymmetrical event shapes that are clearly distinguishable from those in the crackling regime.

Crackling to periodic dynamics in granular media

We study the local and global dynamics of sheared granular materials in a stick-slip experiment, using a slider and a spring. The system crackles, with intermittent slip avalanches, or exhibits irregular or periodic dynamics, depending on the shear rate and loading stiffness. The global force on the slider during shearing captures the transitions from the crackling to the periodic regime. We deduce a dynamic phase diagram as a function of the shear rate and the loading stiffness and associated scaling laws. Using photoelastic particles, we also capture the grain-scale stress evolution, and investigate the microscopic behavior in the different regimes.

Compressive Failure as a Critical Transition: Experimental Evidence and Mapping onto the Universality Class of Depinning

Acoustic emission (AE) measurements performed during the compressive loading of concrete samples with three different microstructures (aggregate sizes and porosity) and four sample sizes revealed that failure is preceded by an acceleration of the rate of fracturing events, power law distributions of AE energies and durations near failure, and a divergence of the fracturing correlation length and time towards failure. This argues for an interpretation of compressive failure of disordered materials as a critical transition between an intact and a failed state. The associated critical exponents were found to be independent of sample size and microstructural disorder and close to mean-field depinning values. Although compressive failure differs from classical depinning in several respects, including the nature of the elastic redistribution kernel, an analogy between the two processes allows deriving (finite-) sizing effects on strength that match our extensive data set. This critical interpretation of failure may have also important consequences in terms of natural hazards forecasting, such as volcanic eruptions, landslides, or cliff collapses.

Depinning Dynamics of Crack Fronts

We investigate experimentally and theoretically the dynamics of a crack front during the microinstabilities taking place in heterogeneous materials between two successive equilibrium positions. We focus specifically on the spatiotemporal evolution of the front, as it relaxes to a straight configuration, after depinning from a single obstacle of controlled strength and size. We show that this depinning dynamics is not controlled by inertia, but instead by the rate dependency of the dissipative mechanisms taking place within the fracture process zone. This implies that the crack speed fluctuations around its average value v_{m} can be predicted from an overdamped equation of motion (v-v_{m})/v_{0}=[G-G_{c}(v_{m})]/G_{c}(v_{m}) involving the characteristic material speed v_{0}=G_{c}(v_{m})/G_{c}^{'}(v_{m}) that emerges from the variation of fracture energy with crack speed. Our findings pave the way to a quantitative description of the critical depinning dynamics of cracks in disordered solids and open up new perspectives for the prediction of the effective failure properties of heterogeneous materials.

Crack growth and energy dissipation in paper

Here, we follow the stable propagation of a roughening crack using simultaneously Digital Image Correlation and Infra-Red imaging. In a quasi-two-dimensional paper sample, the crack tip and ahead of that the fracture process zone follow the slowly, diffusively moving “hot spot” ahead of the tip. This also holds when the crack starts to roughen during propagation. The well-established intermittency of the crack advancement and the roughening of the crack in paper are thus subject to the dissipation and decohesion in the hot spot zone. They are therefore not only a result of the depinning of the crack in a heterogeneous material.

Avalanches and extreme value statistics in interfacial crackling dynamics

We study the avalanche and extreme statistics of the global velocity of a crack front, propagating slowly along a weak heterogeneous interface of a transparent polymethyl methacrylate block. The different loading conditions used (imposed constant velocity or creep relaxation) lead to a broad range of average crack front velocities. Our high-resolution and large dataset allows one to characterize in detail the observed intermittent crackling dynamics. We specifically measure the size S, the duration D, as well as the maximum amplitude [Formula: see text] of the global avalanches, defined as bursts in the interfacial crack global velocity time series. Those quantities characterizing the crackling dynamics follow robust power-law distributions, with scaling exponents in agreement with the values predicted and obtained in numerical simulations of the critical depinning of a long-range elastic string, slowly driven in a random medium. Nevertheless, our experimental results also set the limit of such model which cannot reproduce the power-law distribution of the maximum amplitudes of avalanches of a given duration reminiscent of the underlying fat-tail statistics of the local crack front velocities.This article is part of the theme issue 'Statistical physics of fracture and earthquakes'.

Crack growth in heterogeneous brittle solids: intermittency, crackling and induced seismicity

Crack growth is the basic mechanism leading to the failure of brittle materials. Engineering addresses this problem within the framework of continuum mechanics, which links deterministically the crack motion to the applied loading. Such an idealization, however, fails in several situations and in particular cannot capture the highly erratic (earthquake-like) dynamics sometimes observed in slowly fracturing heterogeneous solids. Here, we examine this problem by means of innovative experiments of crack growth in artificial rocks of controlled microstructure. The dynamical events are analysed at both global and local scales, from the time fluctuation of the spatially averaged crack speed and the induced acoustic emission, respectively. Their statistics are characterized and compared with the predictions of a recent approach mapping fracture onset to the depinning of an elastic interface. Finally, the overall time-size organization of the events is characterized to shed light on the mechanisms underlying the scaling laws observed in seismology.This article is part of the theme issue 'Statistical physics of fracture and earthquakes'.

Thermally activated crack fronts propagating in pinning disorder: simultaneous brittle/creep behaviour depending on scale

We study theoretically the propagation of a crack front in mode I along an interface in a disordered elastic medium, with a numerical model considering a thermally activated rheology, toughness disorder and long-range elastic interactions. This model reproduces not only the large-scale dynamics of the crack front position in fast or creep loading regimes, but also the small-scale self-affine behaviour of the front. Two different scaling laws are predicted for the front morphology, with a Hurst exponent of 0.5 at small scales and a logarithmic scaling law at large scales, consistently with experiments. The prefactor of these scaling laws is expressed as a function of the temperature, and of the quenched disorder characteristics. The cross-over between these regimes is expressed as a function of the quenched disorder amplitude, and is proportional to the average energy release rate, and to the inverse of temperature. This model captures as well the experimentally observed local velocity fluctuation probability distribution, with a high-velocity tail P(v)∼v -2.6 This feature is shown to arise when the quenched disorder is sufficiently large, whereas smaller toughness fluctuations lead to a lognormal-like velocity distribution. Overall, the system is shown to obey a scaling determined by two distinct mechanisms as a function of scale: namely, the large scales display fluctuations similar to an elastic line in an annealed noise excited as the average front travels through the pinning landscape, while small scales display a balance between thresholds in possible elastic forces and quenched disorder.This article is part of the theme issue 'Statistical physics of fracture and earthquakes'.

Crackling Dynamics in the Mechanical Response of Knitted Fabrics

Crackling noise, which occurs in a wide range of situations, is characterized by discrete events of various sizes, often correlated in the form of avalanches. We report experimental evidence that the mechanical response of a knitted fabric displays such broadly distributed events both in the force signal and in the deformation field, with statistics analogous to that of earthquakes or soft amorphous materials. A knit consists of a regular network of frictional contacts, linked by the elasticity of the yarn. When deformed, the fabric displays spatially extended avalanchelike yielding events resulting from collective interyarn contact slips. We measure the size distribution of these avalanches, at the stitch level from the analysis of nonelastic displacement fields and externally from force fluctuations. The two measurements yield consistent power law distributions reminiscent of those found in other avalanching systems. Our study shows that a knitted fabric is not only a thread-based metamaterial with highly sought after mechanical properties, but also an original, model system, with topologically protected structural order, where an intermittent, scale-invariant response emerges from minimal ingredients, and thus a significant landmark in the study of out-of-equilibrium universality.

Aftershock sequences and seismic-like organization of acoustic events produced by a single propagating crack

Brittle fractures of inhomogeneous materials like rocks, concrete, or ceramics are of two types: Nominally brittle and driven by the propagation of a single dominant crack or quasi-brittle and resulting from the accumulation of many microcracks. The latter goes along with acoustic noise, whose analysis has revealed that events form aftershock sequences obeying characteristic laws reminiscent of those in seismology. Yet, their origin lacks explanation. Here we show that such a statistical organization is not only specific to the multi-cracking situations of quasi-brittle failure and seismology, but also rules the acoustic events produced by a propagating crack. This simpler situation has permitted us to relate these laws to the overall scale-free distribution of inter-event time and energy and to uncover their selection by the crack speed. These results provide a comprehensive picture of how acoustic events are organized upon material failure in the most fundamental of fracture states: single propagating cracks.

The multiple microcracking events underlying damage in inhomogeneous brittle materials form characteristic aftershocks sequences obeying similar laws to those in seismology. Here, Barés et al. evidence and explain the same organization in the acoustic noise produced by a single propagating crack.

Heat treatment effect on the mechanical properties, roughness and bone ingrowth capacity of 3D printing porous titanium alloy

The weak mechanical strength and biological inertia of Ti-6Al-4V porous titanium alloy limit its clinical application in the field of orthopedics. The present study investigated the influence of different solution temperatures (e.g. 800 °C, 950 °C and 1000 °C) on the mechanical properties, roughness and bone ingrowth capacity of Ti-6Al-4V porous titanium alloy prepared by Electron Beam Melting. It was found that the compressive and shear strength were promoted with the increase of solution temperature because of the transformed crystallinity of Ti-6Al-4V titanium alloy and phase changes from TiAl to TiAl + TiV. In addition, the topological morphology, surface roughness and wettability of the porous titanium alloy scaffolds were improved after heat treatment, and in turn, the adhesion rate and cell proliferation of bone marrow mesenchymal stem cells were enhanced. Compared with the scaffolds before and after heat treatment at 800 °C, the scaffolds heat-treated at 950 °C and 1000 °C achieved better bone ingrowth, extracellular matrix deposition and osseointegration. These findings indicate the great potential of heat treatment in possessing Ti-6Al-4V porous titanium alloy for orthopedic implant.

Localization of soft modes at the depinning transition

We characterize the soft modes of the dynamical matrix at the depinning transition, and compare the matrix with the properties of the Anderson model (and long-range generalizations). The density of states at the edge of the spectrum displays a universal linear tail, different from the Lifshitz tails. The eigenvectors are instead very similar in the two matrix ensembles. We focus on the ground state (soft mode), which represents the epicenter of avalanche instabilities. We expect it to be localized in all finite dimensions, and make a clear connection between its localization length and the Larkin length of the depinning model. In the fully connected model, we show that the weak-strong pinning transition coincides with a peculiar localization transition of the ground state.

Nonlinear Focusing in Dynamic Crack Fronts and the Microbranching Transition

Cracks in brittle materials produce two types of generic surface structures: facets at low velocities and microbranches at higher ones. Here we observe a transition from faceting to microbranching in polyacrylamide gels that is characterized by nonlinear dynamic localization of crack fronts. To better understand this process we derive a first-principles nonlinear equation of motion for crack fronts in the context of scalar elasticity. Its solution shows that nonlinear focusing coupled to rate dependence of dissipation governs the transition to microbranching.

Local and global avalanches in a two-dimensional sheared granular medium

We present the experimental and numerical studies of a two-dimensional sheared amorphous material composed of bidisperse photoelastic disks. We analyze the statistics of avalanches during shear including the local and global fluctuations in energy and changes in particle positions and orientations. We find scale-free distributions for these global and local avalanches denoted by power laws whose cutoffs vary with interparticle friction and packing fraction. Different exponents are found for these power laws depending on the quantity from which variations are extracted. An asymmetry in time of the avalanche shapes is evidenced along with the fact that avalanches are mainly triggered by the shear bands. A simple relation independent of the intensity is found between the number of local avalanches and the global avalanches they form. We also compare these experimental and numerical results for both local and global fluctuations to predictions from mean-field and depinning theories.

Scaling characteristics of one-dimensional fractional diffusion processes in the presence of power-law distributed random noise

Here, we present results of numerical simulations and the scaling characteristics of one-dimensional random fluctuations with heavy-tailed probability distribution functions. Assuming that the distribution function of the random fluctuations obeys Lévy statistics with a power-law scaling exponent, we investigate the fractional diffusion equation in the presence of μ-stable Lévy noise. We study the scaling properties of the global width and two-point correlation functions and then compare the analytical and numerical results for the growth exponent β and the roughness exponent α. We also investigate the fractional Fokker-Planck equation for heavy-tailed random fluctuations. We show that the fractional diffusion processes in the presence of μ-stable Lévy noise display special scaling properties in the probability distribution function (PDF). Finally, we numerically study the scaling properties of the heavy-tailed random fluctuations by using the diffusion entropy analysis. This method is based on the evaluation of the Shannon entropy of the PDF generated by the random fluctuations, rather than on the measurement of the global width of the process. We apply the diffusion entropy analysis to extract the growth exponent β and to confirm the validity of our numerical analysis.

Universal features of amorphous plasticity

Plastic yielding of amorphous solids occurs by power-law distributed deformation avalanches whose universality is still debated. Experiments and molecular dynamics simulations are hampered by limited statistical samples, and although existing stochastic models give precise exponents, they require strong assumptions about fixed deformation directions, at odds with the statistical isotropy of amorphous materials. Here, we introduce a fully tensorial, stochastic mesoscale model for amorphous plasticity that links the statistical physics of plastic yielding to engineering mechanics. It captures the complex shear patterning observed for a wide variety of deformation modes, as well as the avalanche dynamics of plastic flow. Avalanches are described by universal size exponents and scaling functions, avalanche shapes, and local stability distributions, independent of system dimensionality, boundary and loading conditions, and stress state. Our predictions consistently differ from those of mean-field depinning models, providing evidence that plastic yielding is a distinct type of critical phenomenon.

The universality class for plastic yield in amorphous materials remains controversial. Here authors present a tensorial mesoscale model that captures both complex shear patterns and avalanche scaling behaviour, which differs from mean-field models and suggests a distinct type of critical phenomenon.

Crack propagation through disordered materials as a depinning transition: A critical test of the theory

The dynamics of a planar crack propagating within a brittle disordered material is investigated numerically. The fracture front evolution is described as the depinning of an elastic line in a random field of toughness. The relevance of this approach is critically tested through the comparison of the roughness front properties, the statistics of avalanches, and the local crack velocity distribution with experimental results. Our simulations capture the main features of the fracture front evolution as measured experimentally. However, some experimental observations such as the velocity distribution are not consistent with the behavior of an elastic line close to the depinning transition. This discrepancy suggests the presence of another failure mechanism not included in our model of brittle failure.

Avalanches dynamics in reaction fronts in disordered flows

We report on numerical studies of avalanches of an autocatalytic reaction front in a porous medium. The front propagation is controlled by an adverse flow resulting in upstream, static, or downstream regimes. In an earlier study focusing on front shape, we identified three different universality classes associated with this system by following the front dynamics experimentally and numerically. Here, using numerical simulations in the vicinity of the second-order transition, we identify an avalanche dynamics characterized by power-law distributions of avalanche sizes, durations, and lateral extensions. The related exponents agree well with the quenched-Kardar-Parisi-Zhang theory, which describes the front dynamics. However, the geometry of the propagating front differs slightly from that of the theoretical one. We show that this discrepancy can be understood in terms of the nonquasistatic correction induced by the finite front velocity.

Metamagnetic Anomalies near Dynamic Phase Transitions

We report the existence of anomalous metamagnetic fluctuations in the vicinity of the dynamic phase transition (DPT) that do not occur for the corresponding thermodynamic behavior of simple ferromagnets. Our results demonstrate that key characteristics associated with the DPT are qualitatively different from conventional thermodynamic phase transitions. We also provide evidence that these differences are tunable by showing that the presence of metamagnetic fluctuations and the size of the critical scaling regime depend strongly on the amplitude of the oscillating field that is driving the DPT in the first place.

Inertia and universality of avalanche statistics: The case of slowly deformed amorphous solids

By means of a finite elements technique we solve numerically the dynamics of an amorphous solid under deformation in the quasistatic driving limit. We study the noise statistics of the stress-strain signal in the steady-state plastic flow, focusing on systems with low internal dissipation. We analyze the distributions of avalanche sizes and durations and the density of shear transformations when varying the damping strength. In contrast to avalanches in the overdamped case, dominated by the yielding point universal exponents, inertial avalanches are controlled by a nonuniversal damping-dependent feedback mechanism, eventually turning negligible the role of correlations. Still, some general properties of avalanches persist and new scaling relations can be proposed.

Linking scales in sea ice mechanics

Mechanics plays a key role in the evolution of the sea ice cover through its control on drift, on momentum and thermal energy exchanges between the polar oceans and the atmosphere along cracks and faults, and on ice thickness distribution through opening and ridging processes. At the local scale, a significant variability of the mechanical strength is associated with the microstructural heterogeneity of saline ice, however characterized by a small correlation length, below the ice thickness scale. Conversely, the sea ice mechanical fields (velocity, strain and stress) are characterized by long-ranged (more than 1000 km) and long-lasting (approx. few months) correlations. The associated space and time scaling laws are the signature of the brittle character of sea ice mechanics, with deformation resulting from a multi-scale accumulation of episodic fracturing and faulting events. To translate the short-range-correlated disorder on strength into long-range-correlated mechanical fields, several key ingredients are identified: long-ranged elastic interactions, slow driving conditions, a slow viscous-like relaxation of elastic stresses and a restoring/healing mechanism. These ingredients constrained the development of a new continuum mechanics modelling framework for the sea ice cover, called Maxwell-elasto-brittle. Idealized simulations without advection demonstrate that this rheological framework reproduces the main characteristics of sea ice mechanics, including anisotropy, spatial localization and intermittency, as well as the associated scaling laws.This article is part of the themed issue 'Microdynamics of ice'.

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.

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.

Universal Pinning Energy Barrier for Driven Domain Walls in Thin Ferromagnetic Films

We report a comparative study of magnetic field driven domain wall motion in thin films made of different magnetic materials for a wide range of field and temperature. The full thermally activated creep motion, observed below the depinning threshold, is shown to be described by a unique universal energy barrier function. Our findings should be relevant for other systems whose dynamics can be modeled by elastic interfaces moving on disordered energy landscapes.

Distribution of joint local and total size and of extension for avalanches in the Brownian force model

The Brownian force model is a mean-field model for local velocities during avalanches in elastic interfaces of internal space dimension d, driven in a random medium. It is exactly solvable via a nonlinear differential equation. We study avalanches following a kick, i.e., a step in the driving force. We first recall the calculation of the distributions of the global size (total swept area) and of the local jump size for an arbitrary kick amplitude. We extend this calculation to the joint density of local and global sizes within a single avalanche in the limit of an infinitesimal kick. When the interface is driven by a single point, we find new exponents τ_{0}=5/3 and τ=7/4, depending on whether the force or the displacement is imposed. We show that the extension of a "single avalanche" along one internal direction (i.e., the total length in d=1) is finite, and we calculate its distribution following either a local or a global kick. In all cases, it exhibits a divergence P(ℓ)∼ℓ^{-3} at small ℓ. Most of our results are tested in a numerical simulation in dimension d=1.