## Abstract

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.

## INTRODUCTION

Complexity in materials deformation is important for engineering and involves fundamental nonequilibrium physics. Such phenomena are encountered when samples are loaded beyond the regime of linear, elastic response. Then, metals yield and the plastic deformation before failure is now known to exhibit very complex properties on various scales in time and space (*1*–*4*). The challenges this brings up range from avalanches of plastic deformation to the statistical fluctuations of the yield stress in finite samples to deformation localization. A typical manifestation of localization is the appearance of shear bands and, here, we study the Portevin–Le Chatelier (PLC) effect (*5*, *6*).

The PLC effect implies the creation of deformation bands in a sample (Fig. 1) when it is loaded beyond the yield point: Such bands nucleate and may or may not propagate depending on the class of PLC instability present (*7*, *8*) (in the common classification, type A denotes propagating and types B and C denote nonpropagating bands). The deformation bands are accompanied by material instabilities; in the case of tensile tests, stress drops, which then produce serrated stress-strain curves (Fig. 1B). This kind of strain rate sensitivity (*9*, *10*) arises as a strain rate–dependent phenomenon; moreover, its character and presence are dependent on the temperature. The PLC effect is attributed to dynamic strain aging (DSA) (*11*–*13*), and the crucial physics is in the interaction of the dislocations as the fundamental carriers of plastic deformation with the solute atoms in the alloy (*14*–*16*). On the mesoscopic level, theories of increasing complexity have been proposed such that they would account for the necessary dislocation physics: elementary classes of immobile and “aging,” solute bound dislocations, and mobile ones producing plastic deformation. Such models and a multitude of experiments have been recently introduced to explore the physics of the PLC effect: phases in the band nucleation (*17*–*20*) and dynamics including serrations in the stress-strain curves (*8*, *10*, *21*–*25*), acoustic emission (*25*–*28*) from the effect, and so forth.

## RESULTS

Here, we take a fundamentally different approach of coarse graining, where the bands are reduced to zero-dimensional “particles.” This amounts to studying the propagation velocity signals *v _{b}*(

*t*) of each individual propagating (“type A”) band during a deformation experiment. Our high-resolution experiments based on speckle imaging of the deforming sample (see Fig. 1 and Materials and Methods for details) reveal that the

*v*(

_{b}*t*) signals are reminiscent of crackling noise bursts found in numerous driven systems ranging from propagating cracks (

*29*) and fluid fronts invading porous media (

*30*) to the jerky field-driven motion of domain walls in ferromagnets (see Fig. 1C) (

*31*–

*33*). This is in contrast to the traditional viewpoint where one would characterize the movement of the bands only via their average velocity

*v*(

_{b}*t*) signals, we recover the known phenomenology in that

*34*–

*37*). In this case, we found a power-law increase with the strain rate and an exponential decrease with strain so that they can be summarized as

*p*= 0.6 and ϵ

_{0}= 0.16.

To characterize the properties of the *v _{b}*(

*t*) signals/velocity bursts corresponding to individual bands, we start by considering their average shapes

*T*(where

*t*

_{0}is the start of the band propagation); this is one of the standard quantities used to characterize crackling noise bursts. We find that short-lived bands exhibit an approximately parabolic shape, while considering bands with a longer

*T*results in

How can one theoretically understand the origin and properties of the crackling noise-like *v _{b}*(

*t*) band propagation velocity signals, exhibiting such average temporal velocity profiles? The starting point of our analysis is the empirical observation that the bands tend to propagate essentially as “rigid bodies,” and, hence, a description based on a single degree of freedom, the band position

*x*, is appropriate. This rigid body then moves via overdamped dynamics due to the forces acting on it. As the sample is strained with a constant strain rate

_{b}*k*, which includes the hardening of the sample, which can be incorporated in the simplest form as a linear dependence to the strain

*k*∝ ϵ. As the band propagates along the long axis of the specimen, it samples the random dislocation microstructure it encounters during motion, resulting in a position-dependent random force

*W*(

*x*), with Brownian correlations, 〈

_{b}*W*(

*x*)

_{b}*W*(

*x*

_{b}^{′})〉 = ∣

*x*−

_{b}*x*

_{b}^{′}∣. Collecting these terms, one arrives at an equation of motion for

*x*, which has the same form as the Alessandro-Beatrice-Bertotti-Montorsi (ABBM) model (

_{b}*38*) used as the mean-field description of domain wall depinning in disordered ferromagnets, i.e.

*D*is the disorder strength. The ABBM model (Eq. 2) is known to produce crackling noise or avalanches with power law–distributed sizes and durations (

*31*,

*32*), characterized by

*c*-dependent exponents; for instance, the size distribution scales as

*39*), Eq. 2 can be transformed to a form including a time-dependent noise term

*t*)ξ(

*t*

^{′})〉 = δ(

*t*−

*t*

^{′}). This Eq. 3 has the advantage of allowing one to analytically solve quantities like the average burst shape in the

*k*= 0,

*c*= 0 limit, resulting in an inverted parabola, while a finite

*k*gives rise to a flattening of the shape for long avalanches (

*39*).

### Comparing the model and the experiments

To compare this model with experimental data, we simulate it by “nucleating bands” at random initial positions *L*, and let them propagate in a random direction according to Eq. 3. To mimic effects due to the finite length of the sample, we consider only bands that stop before the end of the sample. This leads to an *L*-dependent cutoff to the “avalanche” distribution; in addition, a cutoff could, in principle, be due to *k* in Eq. 3, but, here, *k* is sufficiently small such that the *L*-dependent cutoff dominates. This then results in a scaling form for the avalanche size distribution (see Materials and Methods for details)

Fixing the model parameters to *c*_{0} = 50 m/s), *k* = *k*_{0}ϵ (where *k*_{0} = 650 s^{−1}), and *D* = 700 mm/s^{3} reproduces the experimentally seen band velocities well (Fig. 2), both in terms of the behavior with increasing strain rate and strain. The number of bands generated with a given strain and strain rate is proportional to the statistics of the experiments. All this is achieved with a very simple linear proportionality of the model parameters to the experimental parameters ϵ and

The model also leads to an excellent agreement in terms of the velocity profiles (Fig. 3): inverted parabola for small *T* and increasingly flattened shape for larger *T*. The slight asymmetry (*29*, *40*) seen in the experimental profiles, where the velocities are slightly larger towards the start of the band lifetime, is naturally not reproduced by the model.

We then proceed to a more extensive statistical characterization of the band propagation velocity bursts, again exploiting the remarkable analogy with standard crackling noise systems. First, we look at the scaling of the avalanche size *S* with the duration *T* in the experiments and in the simulations (where the parameters were fitted to reproduce the band velocity behavior with increasing strain and strain rate). Typically, one would expect 〈*S*〉 ∼ *T*^{2} (*31*, *32*) from the ABBM model but, here, the finite sample of length *L* = 28 mm (and the exclusion of bands with *S* < 2 mm: see Materials and Methods for details) restricts this scaling. In addition, with different strains and strain rates, one would have different prefactors, and averaging over these would not result in the typical scaling form. After performing this average, we see very similar behavior of 〈*S*〉 increasing with duration (Fig. 4A) for both the experiments and the ABBM model, and the increase is slower than 〈*S*〉 ∼ *T*^{2}.

To connect the predictions of the ABBM model to the stress-strain curve, one can also study the scaling of the average stress-drop size 〈Δσ〉 from the stress-strain curves with duration. The one-to-one correspondence between deformation bands and stress drops is broken by the observed multiple simultaneously propagating bands (Fig. 1), but the average scaling seems to be similar for the avalanche sizes and stress-drop sizes (Fig. 4A). The shorter stress-drop durations are due to the simultaneous bands and the short *S* < 2 mm bands that are otherwise neglected from the analysis.

The prediction for the avalanche size distribution in a finite-size sample (Eq. 4) has a driving rate–dependent exponent *c* and *k* in

The instantaneous (band) velocity distribution in the ABBM model is known to be of the form (*33*, *38*, *39*, *41*)

## DISCUSSION

We have studied the statistics of the PLC deformation bands by using a fast imaging technique and by a simple mean-field model of avalanche dynamics. Empirically, we found that the average band velocities scale on average as *p* = 0.6 and ϵ_{0} = 0.16 and exhibit a flattening of the average velocity profile

The finite size of the sample and the dependence of the model parameters on both strain and strain rate hide the known scaling form of the avalanche size. However, our simulations show that both the average PLC band propagation distances and ABBM avalanche sizes scale similarly with duration. The same scaling can also be seen for the stress-drop sizes in the stress-strain curve. Although the one-to-one correspondence between bands and stress drops is lost with multiple simultaneous bands, the average scaling remains the same.

We analytically show that the finite size of the sample introduces a *c* → 0; however, we have shown here that it can be also used to explain the behavior of fronts under strong drive, here the deformation bands. Looking at these results from the viewpoint of classical theories of the PLC effect, how to modify and adapt such models of DSA (*42*) so that they correctly reproduce the kind of stochasticity seen in the band dynamics is an important question. This may be restated so that the “correct” model should be able to reduce to the ABBM used here.

What our results show is that interacting, mobile dislocations create avalanches of deformation in metal alloys. Here, the necessary conditions for this are temperature and strain rate values within a specific window such that propagating or type A PLC bands are observed. Given this, the avalanches follow the paradigm of the mean-field–like ABBM model. The eventual stopping of the band is a random fluctuation and depends on the local, heterogeneous material properties. Thus, the physics of these bands arises from a mixture of external drive, local randomness, and the coarse-grained, collective response of many dislocations. More work is needed in understanding the implications to other PLC band types, and what the practical predictions or consequences are for alloys with different composition (“disorder”) and for samples of different sizes. It is likely that the ABBM exponent *43*). In the same vein, propagating bands of deformation with serrations of the stress-strain curve are also seen in the plastic deformation of amorphous materials (*44*, *45*). An obvious question would be if these also can be shown to follow ABBM-like dynamics with a careful study, but then again, if such bands do not follow this simplest paradigm, that is also of profound interest.

## MATERIALS AND METHODS

### Experimental methods

The laser speckle technique (*35*) was used to observe the bands in a commercial aluminum alloy AW-5754 sample. The samples were laser cut to a flat dog bone shape with the dimensions 28 mm by 4 mm by 0.5 mm for the gauge volume. The samples have a polycrystalline structure with an average grain size of 38 ± 14 μm. The experimental setup is illustrated in Fig. 5.

The samples were tensile-loaded with Instron ElectroPuls E1000 using an Instron Dynacell load cell with a constant displacement rate. The stress and strain were calculated from the displacement and force data provided by the machine. These were recorded with an acquisition rate of 500 Hz, and the samples were held using an initial force of 4 N.

The speckle pattern was recorded with ProtoRhino FlexRHINO DynaMat system, which includes a high-speed camera, a laser, and an field-programmable gate array chip–based unit for data acquisition and storage. The camera had an electronic freeze-frame shutter and a Navitar MVL7000 objective with a macro zoom lens, an aperture of f/2.5, and a spatial resolution of 54 μm. The laser used was a collimated laser diode with a wavelength of 638 nm and a power of 200 mW. The acquisition rates varied around 0.5 to 2.0 kHz.

The speckle images were analyzed using the equal interval subtracting method [similar to (*35*)] where the subtraction was done for consecutive images or with the highest acquisition rate for every other image. A one-dimensional (1D) projection was taken from these subtracted images in a direction perpendicular to the band with the two different band inclinations. This provides two different effective strain rate maps, where the measured quantity

As the band angles and widths were observed to remain very close to constant (the band widths are 0.9 ± 0.1 mm, which is of the order of the sample thickness 0.5 mm), these effective strain rate maps were used to track the band movement as a 1D rigid body. The maximum value of *x _{b}*, therefore also determining the propagation distance

*S*and band duration

*T*. As we are considering type A band dynamics, a propagation distance cutoff of 2 mm was imposed. Bands that propagate less than 2 mm correspond more to the type B regime of nonpropagating bands and were excluded from the analysis. The band velocity signal

*v*was then obtained by numerically differentiating the band position signal, and the average band velocity was calculated simply as

_{b}There are sometimes multiple simultaneous bands present in the sample (see Fig. 1, A and C) that can collide with each other. These collisions usually lead to the disappearance of both bands (except for a few cases right before the sample failure where bands with different inclinations can interact; for more details, see movie S1), and therefore, we have chosen to consider each band as an independent event.

### Fitting the model parameters

The model parameters *c*, *k*, and *D* were fitted to the experimental data in the following way. First, a value for *c* was chosen high enough (with some arbitrary *k* and *D*) so that the average band velocities *k* was chosen so that the normalized average band velocities *c* was increased and the value for *k* was chosen again. Last, the value for *D* was chosen so that that the actual values (not just the slopes) of

In the end, all three of the values were perturbed around the chosen values to make sure that the values represent a local minimum of the difference between the simulation results and experimental data.

### Derivation of the avalanche size distribution

One can calculate the effect of the finite sample length on the avalanche size distribution or, in other words, the band propagation distance by a simple construction. Nucleating bands at a random position on the 1D sample, letting them propagate, and constraining the propagation to the size of the sample directly give the cutoff induced by the finite size sample.

Let *Y* be the event of a band starting at a certain position *X* be the event of the band stopping at a certain position *S* by calculating two convolutions

The conditional probability is handled most simply by splitting it into two portions. After starting, the band goes in either direction with equal probability and, as is known for the ABBM model, travels a distance that is power law–distributed. The joint distribution is then

Normalizing this distribution (from a minimum value *S*_{0} to *L*) gives the full functional form

For the special case of *k* = 0, one can obtain an analytic solution for the exponent in the ABBM model as *31*, *32*), which would here give the avalanche size distribution

## SUPPLEMENTARY MATERIALS

Supplementary material for this article is available at http://advances.sciencemag.org/cgi/content/full/6/41/eabc7350/DC1

This is an open-access article distributed under the terms of the Creative Commons Attribution license, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## REFERENCES AND NOTES

**Acknowledgments:**We thank I. Lomakin for the determination of the sample grain structure.

**Funding:**We acknowledge the financial support from the Academy of Finland through the Centers of Excellence program (project no. 251748), an Academy Research Fellowship (L.L., project no. 268302), and the Academy Projects COPLAST (L.L., project no. 322405) and FLUFRA (M.A., project no. 317464). T.M. acknowledges the support of The Finnish Foundation for Technology Promotion. M.A. acknowledges support from the European Union Horizon 2020 research and innovation program under grant agreement no. 857470 and from the European Regional Development Fund via Foundation for Polish Science International Research Agenda PLUS program grant no. MAB PLUS/2018/8. We acknowledge the computational resources provided by the Aalto University School of Science “Science-IT” project.

**Author contributions:**T.M., P.K., and M.O. performed the experiments, and T.M. analyzed the data and performed the simulations. T.M., L.L., and M.J.A. wrote the manuscript.

**Competing interests:**The authors declare that they have no competing interests.

**Data availability:**All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper may be requested from the corresponding authors.

- Copyright © 2020 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution License 4.0 (CC BY).