Statistics of dislocation avalanches in FCC and BCC metals: dislocation mechanisms and mean swept distances across microsample sizes and temperatures

Plastic deformation in crystalline materials consists of an ensemble of collective dislocation glide processes, which lead to strain burst emissions in micro-scale samples. To unravel the combined role of crystalline structure, sample size and temperature on these processes, we performed a comprehensive set of strict displacement-controlled micropillar compression experiments in conjunction with large-scale molecular dynamics and physics-based discrete dislocation dynamics simulations. The results indicate that plastic strain bursts consist of numerous individual dislocation glide events, which span over minuscule time intervals. The size distributions of these events exhibit a gradual transition from an incipient power-law slip regime (spanning \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\approx$$\end{document}≈ 2.5 decades of slip sizes) to a large avalanche domain (spanning \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\approx$$\end{document}≈ 4 decades of emission probability) at a cut-off slip magnitude \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${s}_{\mathrm{c}}$$\end{document}sc. This cut-off slip provides a statistical measure to the characteristic mean dislocation swept distance, which allows for the scaling of the avalanche distributions vis-à-vis the archetypal dislocation mechanisms in face-centered cubic (FCC) and body-centered cubic (BCC) metals. Our statistical findings provide a new pathway to characterizing metal plasticity and towards comprehension of the sample size effects that limit the mechanical reliability in small-scale structures.

Crystal plasticity is characterized by the onset of intermittent dislocation glide events or avalanches, which lead to prominent strain burst emissions when the sample size is reduced to a few micrometers. The emission probability, P , of a dislocation burst of size s is ascribed to the power-law relation P(s) ∝ s − , where the exponent, , is a universal parameter with scale-invariant value ≈ 1.5 [1][2][3][4][5][6][7][8][9][10][11] . This arguably led to the interpretation in that the size distribution of the dislocation bursts fulfils a mean-field depinning phase transition paradigm characterized by a critical exponent 3,7-10, [12][13][14] . Truncation of intermittent responses is a key feature in the scaling of physical phenomena ranging from earthquakes to neural activity 15 . In crystal plasticity, P(s) is exponentially truncated when the strain bursts approach a cut-off value. The notion of stress-tuned criticality (STC) then arises when this cut-off slip is governed by the applied stress, while volume-tuned criticality (VTC) implies that the cutoff depends on sample size. In both models, tuned criticality may arise at critical values of stress and sample volume 3,7,8,12,13 . Alternatively, in the advent of self-organization processes of the dislocation structure during straining, dislocation glide may occur at a fixed criticality level. The slip distribution then becomes independent of the applied stress, leading to the hypothesis of self-organized critically (SOC) 12,[15][16][17][18] . Along these lines, it may be argued that the self-similar scaling of the dislocation cells arising at large shear strains during single crystal deformation may tentatively result in SOC.

Micropillar compression experiments. Our microcrystal compression experiments, performed under
strict displacement control, enabled the identification of individual dislocation glide events. These are the constituent plastic intermittencies which occur within the strain bursts observed in experiments driven with conventional dynamically-reactive testing systems (see Supplementary Material, Section 1, and the current Results and discussion section for a description on the influence of the straining mode and testing system on the plastic intermittencies). The cylindrical pillar-like microcrystals used in our compression experiments were focused ion beam (FIB) milled from [111]-oriented Cu and Ta single crystals, a [100]-oriented Al single crystal, and a [100]-oriented W single crystal. Pillar-like rectangular microcrystals with square cross-sections were also FIBmilled from a [100]-oriented Cu single crystal and a [123]-oriented Al surface. Beam current in the Ga source was reduced to less than 200 pA in all final machining steps to reduce damage. The cylindrical Al microcrystals had a diameter D ≈ 8 µm whereas the Cu, Ta and W counterparts had D ranging from 0.5 to 20 µm. The effective diameters of the rectangular Cu microcrystals ranged from D ≈ 0.5 to 4 µm. In all cases, the sample height, l , was ≈ 3 D . The rectangular samples were heat treated at 600 °C for 1 h to remove any FIB-induced dislocation loops 31 .
The compression experiments were performed using a modified Alemnis nanoindenter (Alemnis GmbH, Thun, Switzerland) fitted inside the high vacuum chamber of a Zeiss DMS 962 scanning electron microscope (SEM). This system is capable of operating under strict, intrinsic displacement control, thus maintaining indenter tip position during sudden avalanche emissions (Section 1 in Supplementary Material). The nanoindenter was custom modified to enable testing at high temperatures with the sample and indenter temperature matched to minimize any influence of thermal expansion drift on applied displacements 32 . The applied strain rates, ε , ranged from 10 −4 to 10 −3 , which is sufficiently small so as to neglect additional straining produced during avalanche propagation (Section 2 in Supplementary Material). A detailed account of some of the stress-strain curves is published elsewhere 18,31,33 . www.nature.com/scientificreports/ plasticity at extreme dislocation densities. The simulation cell size was 70 × 70 × 40 nm comprising ≈ 11 million atoms for Al and ≈ 10 million atoms for Ta, modelled through the embedded-atom method (EAM) potentials in Refs. 34,35 . All simulations were carried out with the LAMMPS code 36 under an NPT ensemble, where the uniaxial stress σ is obtained from the pressure tensor in the direction of the applied displacement (with a frequency of 200 fs) while the barostat acts in the other main directions of the simulation cell. The Al cells were deformed at RT and the Ta counterparts at RT and 630 °C with ε of 1.5 ×10 6 s −1 . The networks were introduced into the computational cells by bringing one of the free surfaces into contact against a spherical indenter until a dense defect network attained (see 37 and Sects. 6 and 7 in the Supplementary Material). Upon application of periodic boundary conditions to the simulation cells, the stress-strain curves were taken to represent bulk material responses. The freestanding MD simulations were then exclusively employed to assess slip trace formation at the surface. Details of the strengthening provided by the incepted networks (see Ref. 38 ) are given in Sect. 6 of the Supplementary Material. The 3D DDD simulations were conducted with an in-house version of the open source code ParaDis 30,39 to capture confined plasticity responses where dislocation annihilation occurs at the free surfaces. All dislocation reactions are planar in this in-house version of the open source code, further incorporating a set of atomisticallyinformed, physics-based cross-slip mechanisms, the details of which are described in 30 . The assumed elastic material properties matched those of FCC Ni (Young's modulus E = 210 GPa, and Poisson's ratio ν = 0.31). The stress-strain curves were obtained by uniaxially compressing rectangular computational cells with a square cross-section under a nominal strain rate control of ε = 200 s −1 (where the stress σ along the straining direction is evaluated with a frequency of 5 ns). Free-standing [100]-oriented microcrystals were compressed with D = 2 and 5 μm with aspect ratio l/D = 3 . The initial dislocation microstructure was introduced as Frank-Reed (FR) sources having random line directions randomly assigned to one of the twelve FCC slip systems. The dislocation density, ρ , was varied in the range of 10 14 m −2 to reproduce the attainment of bulk-like plasticity and the slip distributions from the experiments with Cu microcrystals, while the applied ε remained fixed at 2 ×10 2 s −1 . DDD simulations under severe dislocation starvation were also performed in smaller microcrystals with D = 1 μm ( ρ = 10 12 m −2 ). Since the long-range image field resulting from the traction-free boundary conditions give rise to relatively small corrections in the current simulation cells 40,41 , this image field is ignored here.   31 . Notice the development of stress drops in the stress-strain curves, a feature that is associated with recurrent single-ended dislocation source operation and the occurrence of sudden dislocation network destabilization processes, which lead to large values of the Pearson correlation coefficient, r (see Fig. 5). (c) BCC plasticity across sample sizes 33 . Increasingly fluctuating stress-strain curves, characterized by the largest r values, are observed with decreasing micropillar diameter and increasing temperature. The inset illustrates the emission of an individual slip event under progressive straining (with the current, intrinsic displacement controlled testing system) characterized by Young's modulus, E, stress drop, Δσ, and plastic strain increment, ε. The schematic representation of a cylindrical micropillar illustrates the alternating screw/edge dislocation characters of an active singleended source revolving at the micropillar center. This leads to the onset of wavy slip as shown for D = 0.5 μm at 400 °C. Heterogeneous slip trace distributions along the <111> {112 } slip system family is also illustrated for D = 5 μm at 400 °C. www.nature.com/scientificreports/ tions, an individual slip event involves sudden instability of the dislocation arrangement producing stress drop, �σ , at fixed displacement (see inset to Fig. 1c). The onset of a slip magnitude, s , thus reduces the sample height, l , so as to induce stress drop, �σ , as the elastic energy is released 37,42 . This process is described through

Measurement
where SF is the Schmid factor and E is the effective elastic modulus of the microcrystal (see Sect. 2 in Supplementary Material). Equation (1) indicates that the stress drop magnitude �σ increases when the same slip event size occurs in samples of a smaller height. This supports the notion in that the plastic intermittencies may only become noticeable in small (microcrystalline) sample sizes. According to coarse-grained crystal plasticity 43 , the following correlation between the slip magnitudes and the shear strain γ can then be expressed where b is the magnitude of the Burgers vector, N is the total number of mobilized dislocations in a given avalanche, L is the mean distance swept by each dislocation, A s = LD is the swept area, and V is the sample volume. Equation

Results and discussion
Individual avalanche emissions versus dynamically-driven bursts. An important finding from the current MD and DDD simulations is that microscale plasticity is essentially characterized by the onset of individual avalanche events associated with recurrent stress drops �σ . Our simulations show that these individual avalanches are emitted regardless of the externally applied deformation rate, which may become several decades greater than those in conventional dynamically-driven experiments. In these experiments, the sample is rapidly strained after the onset a plastic instability, the magnitude of which is then controlled by the activation the feedback loop and mechanical dynamics of the testing system (see the discussion given in Sect. 1 of the Supplementary Material).
Since the MD simulations indicate that the propagation time frame is ≈ 20 ns for a medium-sized avalanche event sweeping across a microcrystal with D = 2 µm, it is argued that a larger, dynamically-driven, strain burst ( �ε ≈ 0.1) occurring over a time frame several decades greater would be comprised of more than a hundred of individual avalanches. Within the range of time frames under consideration, each constituent individual avalanche would essentially propagate at fixed displacement, as reproduced with the current experiments, which results in the C(s) distributions being much narrower than those of the dynamically-driven strain bursts. A correspondence is therefore anticipated between the C(s) distributions from experiments and simulations performed under strict displacement control and those of the constituent individual avalanche events concealed in the dynamically-driven experiments. Much smaller ν(s) are however reached in the current displacement controlled experiments, where low strain-rates of ε = 10 −3 to 10 −4 s −1 are externally applied, than in the computational simulations or in the dynamically-driven strain bursts where ε is several decades greater.
Finally, the MD simulations illustrate the energetics of the current individual slip events. It is found that within the large avalanche domain, the elastic energy release rate is 20 times greater than the elastic energy storage rate prior to the onset of the intermittency. The ratio between these elastic energy rates is reduced to 4 within the small slip domain, while it approaches unity for s < 0.1b where plastic deformation proceeds through smooth dislocation bowing processes ( �σ → 0). The cut-off slip size s c then sets the exponential truncation for s > s c , which results in the onset of a large avalanche regime with reduced emission probability. It is noted that while a similar formulation as in Eq. (3) was proposed in Ref. 4 , the distribution of the individual dislocation glide events is in better accord with the currently employed relation   Fig. 2 using a different set of [100]-oriented Cu microcrystals tested as a function of temperature and size. These microcrystals were heat treated to remove the FIB-induced dislocations, and exhibited confined plasticity when subjected to uniaxial compression (see Fig. 1b). As compared to in the Cu micropillars in Fig. 2, the current incipient slip regime (κ ≈ 0.4) is clearly prolonged to larger cut-off slip sizes in the range of s c ≈ 10b-20b. This feature becomes more evident as the testing temperature is increased (D = 3.5 μm at T = 400 °C). Notice that when the microcrystal diameter is reduced towards sub-micrometer sizes, the cut-off slip s c is drastically reduced. This is supported by the DDD simulations performed with a reduced number of dislocation segments to capture such source starvation responses (grey circles, where D = 1 μm, ρ = 1 × 10 12 m −2 and N = 1 along with the mechanisms illustrated in Fig. 4e,f). The inset provides avalanche emission frequency function C(s) from the experiments in the main figure. The scattering of the slip events as a function of the applied shear stress τ in the active slip systems in shown in (b), in support of stress-tuned criticality (STC). Larger slip events thus occur when τ > τ c .  (2), where it follows that s c = b(NL) c /D . The substitution of this relation into Eq. (3) yields the explicit dependency of the slip distributions upon the microsample diameter, D , and the total distance travelled by the mobile dislocations in an intermittency with s = s c , as given by the product (NL) c . Along the seminal dislocation storage-recovery model 45 , detailed investigations in FCC crystals have revealed the intricate connection between a mean, deterministic value of the dislocation swept distance L and the strain hardening response produced by the expansion of pinned dislocation segments across an arrangement of forest dislocations, leading to the formation of orientation-dependent junctions 19,46,47 . In the light of these investigations, we find L = 0.7 μm and 2.3 μm for the Cu and Al microsamples with D = 2 μm and 8 μm, respectively. This is consonant with the estimated ρ ≈ 7 × 10 14 m −2 and ρ ≈ 8 × 10 13 m −2 which develop under the applied shear stresses of τ • ≈ 65 MPa and 27 MPa in the microcrystal compression experiments of Cu and Al, respectively (see Sect. 6 in the Supplementary Material). It is noted that the strain hardening parameter θ is then reduced as compared to that attaining at smaller dislocation densities 18 . Mechanisms and statistics of bulk-like plasticity in FCC microcrystals. Figure 2 shows the C(s) distributions of the individual slip events from experiments and simulations of bulk-like microscale plasticity, where a rather uniform arrangement of non-localized, fainting slip traces develops along the microcrystal height in conjunction with plastic barreling (Fig. 1a). A master slip distribution is constructed for microcrystals that exhibit intense interactions between the mobile dislocations and the forest dislocations. This distribution is fitted with s c ≈ 1 b , and incorporates the experimental results from Cu microcrystals ( D ≈ 2 μm) compressed along the [111] orientation, the DDD simulations for freestanding Ni microcrystals containing a large dislocation density, and the MD simulations for an extremely entangled Al bulk modeled through periodic boundary conditions.
The dislocation mechanisms underlying the above master slip distribution are shown in Fig. 4. The MD simulations illustrate that the incipient slip domain is distinguished by correlated dislocation glide events involving the expansion of pinned loops over L < 30 nm. These processes occur throughout the entire cell volume (Fig. 4a-d). In the MD simulations, stress drops �σ are no longer detected below the minimum slip size s min ≈ 0.07b where the stress-strain curves acquire a continuous appearance and deformation is accumulated through reversible bowing of pinned dislocation segments along with the onset of gradual, smooth dislocation glide events over L ≈ 5 nm. As the sudden release of elastic energy vanishes in these processes, they are not regarded as true avalanches.
While the above observations concern bulk-like responses, the DDD simulations incorporate the role of the free surfaces in avalanche propagation at dislocation density levels rendering similar C(s) distributions as those measured experimentally with the Cu microcrystals. Overall, the DDD simulations indicate that when s → s c , the individual slip events are comprised of a total of N = 2 to 3 mobile dislocations, where L 0.5 μm according to Eq. (2). At the cut-off slip size s c ≈ 1 b , we then find that N c ≈ 3 so that L c ≈ 0.7 μm (Eq. (2)), which is in excellent accord with the deterministic value of the mean swept distance estimated in the above section. With increasing slip size ( s > s c ), single-ended source operation accounts for an increasing fraction of the total slip magnitude in the freestanding microcrystal. All of the observed slip events producing s > s c then consist of N ≥ 3 mobile dislocations, reaching the maximum value of N = 17 with L ≈ 1.2 μm for the largest measured avalanche event. The DDD simulations finally show that while dislocation storage and dislocation annihilation processes develop within the incipient slip regime, large avalanche emissions ( s > s c ) are characterized by net dislocation storage (Fig. 4j).
The general slip distribution function in Eq. (3) is maintained for [100]-oriented Al microcrystals with D = 8 µm, where s c increases from ≈ 1 b to ≈ 3.5b indicating the onset of moderate dislocation interactions along with a further increase in L c from ≈ 0.7 μm to ≈ 2.3 μm (as described above). These microcrystals deform through bulk-like plasticity, which is evidenced by the occurrence of non-localized slip traces and sample barreling (Fig. 1a). Since the strain hardening response is further influenced by the dislocation stacking fault width, d , it is concluded that Al microcrystals with d ≈ 1 nm, thus exhibiting enhanced cross-slip and reduced interactions between mobile and forest dislocations, will tend to produce greater slip events than Cu microcrystals with d ≈ 3 nm. According to the dislocation storage-recovery model, cross-slip would increase the estimated L c by ≈ 15% in the Al crystals.
The following discussion concerns the scaling of the slip distributions as a function of the microsample size. In this sense, bulk-like microcrystal plasticity is known to fulfill a classic continuum mechanics description where the deformed state is prescribed by the shear strain γ in the active slip systems. This implies that for a given slip event occurring in a hypothetical, representative material volume element, the number of mobile dislocations, N , increase in proportion with the element's diameter, D , so that the same slip magnitude, s , is produced (Eq. (2)). Hence, the shear strain γ = s/l is preserved irrespective of sample diameter D . Since a rather homogeneous slip distribution develops along the sample height l (as evidenced from the diffuse slip traces in Fig. 1a), increasing microcrystal height gives rise to a proportional increase in the avalanche emission rate. In this sense, it is noted that in microcrystals deformed with the same ε , a three-fold increase in l raises ν(s) by the same three-fold factor (inset to Fig. 2). The shear strain γ is consequently maintained irrespective of l . The aforementioned scaling between N and D obviously breaks down in sufficiently small samples where an incipient slip event ( s < s c ) is www.nature.com/scientificreports/ www.nature.com/scientificreports/ already comprised of a single dislocation. This is inconsequential to our statistical analyses since the incipient slip domain adheres to a scale invariant power-law distribution.
Equations (2) and (3) prescribe that the avalanche distributions scale as where the cut-off, mean dislocation glide distance L c can be estimated through the measured stress-strain curves, as described above. As illustrated through the scaling of the slip distributions of Cu and Al microcrystals with different diameter D via Eq. (4), see Fig. 2, it is found that the linear density, N c /D , of mobile dislocations carrying a slip magnitude s c remains constant within the current range of ρ ≈10 14 to 10 15 m -2 where bulk-like microcrystal plasticity develops.
Statistical transition towards confined plasticity in FCC microcrystals. Figure 2 contains the slip distribution from Cu microcrystals with the [100] orientation exhibiting a heterogeneous arrangement of localized slip traces along the sample height l . The stress-strain curves under such confined plasticity responses where dislocation network development is counterbalanced by surface annihilation are then characterized by the onset of large stress drops, �σ , and θ = 0 . Under these conditions, our results show that the cut-off slip increases from the value of s c ≈ 1 b , measured for the above bulk-like strain hardening microcrystals, to the current s c ≈ 12 b in confined plasticity. As discussed next, this is consistent with the conception in that when strain hardening is hindered, plasticity is increasingly sustained by the activation of a limited number of single-ended sources in a steady, forest dislocation arrangement. The measured function ν(s) then shifts towards larger levels while dislocation glide becomes highly localized along the sample height (inset to Fig. 2). The slip distributions from a second set of experiments in Cu microcrystals, tested as a function of the characteristic sample size D and temperature (Fig. 3), are key to the statistical analysis of confined plasticity. These microcrystals also exhibit large �σ and the heterogeneous slip trace arrangement in Fig. 1b. The measured slip distributions indicate that when D is reduced from 3.5 µ m to 1 µ m, a two-fold decrease in s c is attained (see Fig. 3), which may only occur if the product (NL) c exhibits a seven-fold decrease, as prescribed by Eq. (2). Even under the assumption that the revolution of single-ended sources about their pinning points in smaller microcrystals eventually leads to a smaller critical distance L c , the very significant seven-fold decrease in (NL) c suggests that the number of mobile sources, N c , would also have to decrease when D is reduced to 1µ m. Along the lines of the dislocation mechanisms, the DDD simulations indicate that within the large avalanche domain, the total dislocation swept distance that is sustained through single-ended source operation becomes several times greater than the sample diameter (implying that NL ≫ D ). This is illustrated in Fig. 4g-i, where the pivoting of single-ended sources about multiple junctions with immobile dislocations results in the activation of dislocation segments at the intercepting cross-slip systems. Finally, following the thermally-activated nature of cross-slip, an increase in the total dislocation swept distance (NL) c is anticipated to occur at elevated temperatures, as dislocation annihilation precludes stable junction formations and the mobile segments are disseminated across the active slip systems throughout the sample height. This is supported by the increase of the measured s c and the attendant rise in ν(s) (see Fig. 3 for D = 3.5 µm) along the onset of more noticeable localized slip (Fig. 1b). Increases in temperature thus facilitate the attainment of confined plasticity.
As the sample diameter decreases further towards submicrometer sizes (i.e. D ≈ 0.6 µm in Fig. 3), the slip distributions exhibit drastic reductions in s c and ν(s) which are attributed to severe dislocation source starvation. Therefore, the marked decrease in the number of mobile dislocations, N c , in the transition from micrometer to submicrometer sizes significantly reduces the product (NL) c , which results in the currently measured smaller s c values. This is supported by the DDD simulations performed with D = 1 µm, where only one single-ended source recursively operates ( ρ = 10 12 m -2 ; N = 1 ) under a fixed value of the shear stress, resulting in markedly serrated stress-strain curves and a similar s c ≈ 4b as in the above experiments with D ≈ 0.6 µm (see Fig. 4e,f). These simulations eventually show the onset of cross-slip events of the active source and the occurrence of collinear interactions due to the presence of a few secondary dislocation segments. The latter leads to the segmentation of the slip events, which reduces s c .
Assessment of the correlation between the magnitude of the slip events and the applied shear stresses in the active slip systems is finally illustrated in Fig. 2b,c. We choose to evaluate this feature through Pearson correlation coefficient, r , rising towards unity when an increasingly linear relationship is reached between the scattered slip population and the stress, and the p-value prescribing the likelihood in that this relationship is merely coincidental. In the FCC microcrystals in Fig. 2c, it is found that the same scattered slip distribution develops irrespective of the applied stress level, a feature prescribed by small r < 0.2 or by large p > 20% values. This characterizes the development of bulk-like plasticity both under intense and moderate dislocation interactions, where the size of the individual avalanche events is statistically unaffected by the applied stress, along the lines of SOC. It is then noted that the greatest r values from our analysis pertain to the smaller microcrystals deforming well within the confined plasticity regime, where r lies in the 0.3 to 0.4 range while p < 5%. These values suggest the onset of stress dependent slip distributions, where the likelihood for large avalanche emissions ( s > s c ) correlates with the application of large shear stresses that are greater than a critical τ c value (see Fig. 3b). Therefore, the C(s) distributions measured at larger stress levels would tentatively exhibit a higher s c , in support of STC.
Along the above analyses, it is noted that a single displacement-controlled experiment does not provide a significant number of slip events within a window of stress so as to enable assertion of the possible stress-binning of function C(s) . Therefore, it is not statistically sound to evaluate whether the avalanches adhere to STC or SOC by comparing between the averaged slip magnitudes attained within different windows of stress. The variability in  Fig. 6, it is argued in the following discussion that the measurement of different s c values in BCC microcrystals is the signature of the distinct mobility of the screw dislocations as a function of temperature, sample size and the applied stress level. In submicrometer-sized pillars deformed at RT, plastic deformation proceeds through single-ended source operation 26,27 at applied stresses reaching the distinctive large Peierls barrier [48][49][50] . Plasticity thus sets-in at σ 2000 MPa, as shown in Fig. 1c for D ≤ 0.5 µm, where the character of the mobilized single ended-sources alternates from pure edge to pure screw (inset to Fig. 1c). It is then argued that the recurrent cross-kinking of screw segments triggered under such large applied stresses 51 leads to the wavy slip trace emissions observed in the W microcrystals (see Fig. 1c where similar slip patterns arise in the RT tests).
When D 3 μm, full dislocation loops can be effectively accommodated within the microsample. At RT, the MD simulations suggest that the plastic intermittencies are produced under the smaller applied stresses yielding anisotropic loop expansions along the highly mobile edge segments of the loops (Fig. 6d,e). Although the slip traces remain globally straight (Fig. 1c), atomistically wavy patterns arise through cross-kinking mechanisms of the surface-intercepted screw segments (Fig. 6a). It is noted that while the enhancement of cross-kinking at large applied stresses 51 is consistent with the observed wavy glide in the W microcrystals (3 μm < D < 5 μm) deformed at RT, where σ ≈ 1500 MPa, planar glide was found to prevail in the Ta microcrystals of similar size that exhibited plastic deformation at the smaller stresses, σ ≈ 300 MPa, where cross-kinking is hindered.
A highly dense dislocation network may also develop in the larger microcrystals ( D ≈ 5 μm), as modelled through the MD simulations with periodic boundaries that mimic bulk-like responses. This RT dislocation network is characterized by long screw segments (Fig. 6c) in accord with experimental findings in macroscale samples 48,52 . The MD simulations with free-standing cells that contain similar entangled dislocation networks then predict the predominance of nanoscale, wavy slip features as the meandering mobile dislocation segments www.nature.com/scientificreports/ released from the dislocation structure intercept the surface (Fig. 6b). Although these wavy glide features facilitate slip trace dissemination throughout the sample length along with the activation of cross-kinking in the W www.nature.com/scientificreports/ microcrystals deforming at large stresses ( σ ≈ 1500 MPa and D ≈ 5 μm), the localization of these wavy deformation features within individual slip traces results in the planar deformation mode that was experimentally observed in the Ta microcrystals of similar size, where plasticity was triggered at the smaller σ ≈ 300 MPa in the absence of cross-kinking. The development of planar, localized glide has been observed in bulk BCC crystals 53 which also exhibit the smooth stress-strain curves in Fig. 1c for D ≈ 5 μm. Therefore, slip localization is not necessarily associated with confined plasticity.
Concerning the statistical analysis of the slip distributions, Fig. 5 shows the onset of higher s c and ν(s) in small BCC microcrystals (0.5 μm < D < 2 μm) where plasticity is sustained at the large applied stresses in Fig. 1c. This implies the attainment of a larger, characteristic total dislocation swept distance (which is statistically measured through the product (NL) c ) as compared to in microcrystals with D ≈ 5 μm, deforming under smaller stresses. It is noted that the larger s c in the smaller microcrystals is not a direct outcome of the reduction in the diameter D (prescribed by the relation s c = [b(NL) c /D] ∝ 1/D ), as the measured three-fold increase in s c resulting under the maximum ten-fold decrease in D can only occur under a significant decrease of the product (NL) c ). Finally, smaller s c values are measured within the range of D < 3 μm when the activation of cross-kinking and the onset of wavy glide become more evident (cf. RT experiments where D ranges from 2 to 3 μm in Fig. 5a for Ta-planar glide at larger s c -with Fig. 5b for W-wavy glide at smaller s c ). At fixed microcrystal size, cross kinking is thus associated with a reduction in (NL) c .
The MD simulations reveal that irrespective of the testing temperature, the incipient slip domain is governed by the gliding of individual screw and edge dislocation segments over a mean free path length L of more than 5 atomic lattices, yielding s 0.2b . The distinctive features of the high temperature MD simulations are the increased waviness of the screw segments (cf. Fig. 6c,f) owing to the thermally and stress-activated nature of kink-formation and cross-kinking phenomena 48,51 , as well as the marked rise in the velocity of the mobile dislocations. The large avalanche domain is then associated with the gliding of several dislocation segments along predominantly edge sides, a feature that is facilitated by the enhanced mobility of the lateral, partnering screws. The edge segments thus sweep across L values approaching a hundred of atomic lattices, where the activation of the first dislocation segment triggers gliding of secondary segments. Moreover, the dislocation entanglement is reduced during high temperature avalanche emissions as the mobile segments become mutually annihilated (cf. Fig. 6f,g along with the Supporting Movies S2 and S3). Large stress drops are consequently produced (Fig. 6h) where avalanche arrest occurs in less entangled arrangements that withstand smaller applied stresses. The above results are consistent with the experimentally measured increase in the cut-off slip magnitude s c in the high temperature avalanche distributions, along with a surge of large avalanche emissions as evidenced through the measured larger ν(s) (see the Fig. 5a,b). It is thus argued that the high temperature avalanche distribution is characterized by greater N c and L c values which markedly increase parameter (NL) c .
The influence of sample size on the slip distribution is less clear at elevated temperatures (Fig. 5a), where even a reversed trend is found indicating the occurrence of larger slip events (characterized by higher s c values) as the microcrystal size is increased (Fig. 5b). To rationalize these results, it is important to consider that a thermallyassisted transition from homogeneous (wavy glide) to heterogeneous (planar glide) slip trace distributions is found in BCC W with increasing microcrystal size (Fig. 1c). The wavy glide features observed in the smaller W microcrystals ( D < 1 μm) are then consistent with the onset of cross-kinking at the large applied stresses ( σ > 2000 MPa) needed to sustain plastic deformation for T = 400 ºC ( T/T m ≈ 0.11 ). On the other hand, planar glide prevails in the larger W microcrystals ( D ≈ 5 μm) where cross-kinking is hindered under smaller applied stresses ( σ < 1000 MPa). In accord with the above discussion for the RT experiments where planar glide is associated with a higher s c , the larger W microcrystals tested at elevated temperatures then exhibit higher s c and (NL) c values than the smaller microcrystals (Fig. 5b). Finally, the role of temperature in facilitating large avalanche emissions also becomes evident when comparing the slip distributions from Ta and W microcrystals that exhibited similar planar glide responses (cf. Fig. 5a,b for D > 3 μm). Higher s c and (NL) c are thus measured in the W microsamples tested at a greater homologous temperature ( T/T m ≈ 0.11) than the Ta counterparts ( T/T m ≈ 0.07).
Our RT experiments finally show that function C(s) is increasingly binned by stress with decreasing sample size (cf. Fig. 5d,e), so that larger avalanches are emitted when the applied stress is greater than some critical level ( τ > τ c ). Since this feature gradually vanishes in the larger BCC microcrystals with small s c values, the slip distributions tentatively exhibit a transition from STC to SOC. Increasing temperature then results in statistically larger avalanches for τ > τ c , supporting the prevalence of STC (Fig. 5c).

Concluding remarks
Prior investigations of plastic intermittencies involved the large strain bursts triggered with dynamically-reactive testing systems. Under these driving conditions, the microcrystal is rapidly compressed as a function of the system's dynamic response, which ultimately governs the termination of the instability and its overall size. We argue that a dynamically-induced plastic instability or strain burst is comprised of many of the currently measured individual dislocation avalanche events, propagating over minuscule time scales, which are ubiquitous to crystal plasticity. Our investigation provides in-depth statistical and mechanistic understandings to these individual ultrafast dislocation avalanche events, extending our knowledge of intermittent plasticity.
The following are the salient conclusions from our work: www.nature.com/scientificreports/ multiple dislocation segments that simultaneously glide at numerous locations in the microcrystal without significant dislocation density build-up. 2. A large avalanche domain with reduced emission probability sets-in upon the truncation of the incipient slip regime at the cut-off slip size s c , as prescribed though the universal avalanche distribution function in Eq. (3). It is shown that when the individual dislocation avalanche events are recorded, the cut-off slip s c becomes the essential parameter prescribing the influence of crystalline structure, testing temperature and microsample size-vis-à-vis the onset of collective dislocation glide, dislocation interaction and surface dislocation annihilation mechanisms-in intermittent crystal plasticity. This provides insight into the role of the underlying slip phenomenology upon the dislocation avalanche distribution 14,54-59 . 3. Bulk-like plasticity in FCC Al and Cu microcrystals is characterized by a homogeneous slip trace distribution along with sample barreling. Large avalanche propagation then results in the storage of dislocation segments through junction forming processes. In continuum microcrystal plasticity analyses, the cut-off slip s c attains through the emission of a slip event comprised of N c mobile dislocations which sweep across a mean distance L c in a sample of diameter D . Since s c = b(NL) c /D , a statistically sound measure to the mean dislocation swept distance is obtained through the characteristic slip size s c of the distribution. This distance is the key continuum crystal plasticity parameter to the assessment of dislocation junction formation and strain hardening within the seminal dislocation storage-recovery model 45 . In the absence of a marked increase of the dislocation density ρ with plastic straining, as in the case of the current microcrystals with moderate θ , the avalanche distributions exhibit fixed s c and L c . Since the cut-off slip s c is then independent of the effective shear stress τ in the active slip systems, intermittent plasticity becomes a spontaneous self-organized process, in support of SOC. It is emphasized that while L c and s c would gradually decrease if plastic straining involved pronounced dislocation storage at increasing applied stresses ( θ ≫ 0), this behavior clearly challenges the distinctive increase in s c triggered at high stresses under STC and confined plasticity, as summarized next. 4. Pure FCC microcrystals that deform under confined plasticity exhibit greater values of the cut-off slip s c along with a transition from uniformly distributed to highly localized slip trace patterns. Confined plasticity attains under insufficient dislocation interactions, so that network development is hindered by surface dislocation annihilation when L c > D/2. A marked increase in the total distance travelled by the mobile dislocations is observed due to recurrent single-ended source operation, resulting in higher values of the product (NL) c . Our experiments and DDD simulations further show that when the sample size is decreased towards the submicrometer scale, severe dislocation source starvation leads to a drastic reduction in s c . The natural scaling of the number of mobile dislocations with microcrystal diameter, characterizing continuum bulk-like crystal plasticity thus breaks down. Our experimental results suggest that the onset of the large avalanche domain under confined microcrystal plasticity is governed by the shear stress in the active slip systems. Since greater slip events result from single-ended source operation in conjunction with dislocation network destabilization processes at increasing applied stresses, a higher s c is measured when τ surpasses a critical τ c level, in support of STC. The parameter L c then becomes a function of the microsample size. It is finally noted that while the homogenization of the slip trace distribution is favored under wavy glide, the suppression of cross-slip in alloyed FCC metals that exhibit planar glide 60 and moderate θ facilitates slip localization, arguably triggering the onset of confined microcrystal plasticity. 5. The slip distributions in BCC microcrystals are governed by the sluggish glide of the screw dislocation segments. In spite of the markedly different BCC dislocation glide phenomenology, it is found that the range of s c values measured under bulk-like BCC plasticity ( D 5 μm) are similar to those from FCC counterparts exhibiting marked dislocation interactions. Larger values of the cut-off slip s c are then observed in the smaller BCC microcrystals, where plasticity is triggered through single-ended source operation under the large shear stress level surpassing the Peierls barrier of the screw segments. This indicates an increase of the product (NL) c . Moreover, a transition from planar to wavy slip is eventually observed in the smaller microsamples, denoting the onset of cross-kinking mechanisms. This is characterized by the decrease of the cut-off slip s c and the product (NL) c . Finally, a pronounced raise of s c is observed at elevated temperatures, evidencing an increase in the number of mobile dislocation segments, N c , in conjunction with a higher L c . This is in accord with the enhanced mobility of the screw segments observed in the high temperature MD simulations. In the line of the above slip phenomenology, the occurrence of planar glide at high temperatures facilitates the onset of stress-dependent avalanche distributions, where higher s c and L c are promoted for τ > τ c , in support of STC.