Cooperative deformation in high-entropy alloys at ultralow temperatures

See allHide authors and affiliations

Science Advances  27 Mar 2020:
Vol. 6, no. 13, eaax4002
DOI: 10.1126/sciadv.aax4002


High-entropy alloys exhibit exceptional mechanical properties at cryogenic temperatures, due to the activation of twinning in addition to dislocation slip. The coexistence of multiple deformation pathways raises an important question regarding how individual deformation mechanisms compete or synergize during plastic deformation. Using in situ neutron diffraction, we demonstrate the interaction of a rich variety of deformation mechanisms in high-entropy alloys at 15 K, which began with dislocation slip, followed by stacking faults and twinning, before transitioning to inhomogeneous deformation by serrations. Quantitative analysis showed that the cooperation of these different deformation mechanisms led to extreme work hardening. The low stacking fault energy plus the stable face-centered cubic structure at ultralow temperatures, enabled by the high-entropy alloying, played a pivotal role bridging dislocation slip and serration. Insights from the in situ experiments point to the role of entropy in the design of structural materials with superior properties.


Materials’ mechanical behaviors under extreme conditions have been a subject of intensive research. At high temperatures, plastic deformation is dictated by atomic diffusion, which causes degradation of strength, elongation, phase transformation, and precipitation (1). At low temperatures, a ductile-to-brittle transition is expected because the atomic mobility vanishes, limiting plastic deformation through, e.g., movement of dislocations (2). However, in complex materials, other deformation mechanisms become competitive at low temperatures, providing alternative mechanisms to deform. High-entropy alloy (HEA) is a case in point.

HEAs, which consist of multi-principal elements, are an intriguing class of structural materials with excellent strength—ductility combination, high fracture toughness, resistance against corrosion and hydrogen embrittlement (39). They also have good stability against radiation damage, hence deemed potential materials for advanced reactor applications (10, 11). Refractory HEAs are promising candidates for high-temperature applications, e.g., for aero-engines, due to their thermal stability at elevated temperatures (3, 12, 13). Superconductivity has also been reported in Ta-Nb-Hf-Zr-Ti HEA, thus extending their possible applications beyond structural materials (14). Despite the complex chemistry, HEAs can form a single-phase solid solution with an incredibly simple lattice. For example, CrMnFeCoNi HEA, also known as Cantor alloy, has a face-centered cubic (fcc) structure (4, 6).

At room temperature, the deformation in the quinary CrMnFeCoNi and quaternary CrFeCoNi alloys is dominated by dislocation activities (4, 15, 16). At 1000 K, the major deformation mechanism is diffusion-controlled dislocation creep (17). At liquid-nitrogen temperature, CrMnFeCoNi shows both increased strength and ductility (4, 15). Although the activation of twinning, in addition to dislocation slip, is considered the main reason for the unusual ductility (4, 15), other mechanisms may also contribute, especially when the temperature is further lowered. In stainless steels, for example, martensitic transformation can occur, leading to the so-called TRIP (transformation-induced plasticity) effect (18). First-principle ab initio calculations by Huang et al. (19) suggested fcc to hcp (hexagonal close-packed) phase transformation in CrMnFeCoNi HEA below 77 K. At very low temperatures, such as liquid-helium temperature, serrated deformation has been observed, which usually deteriorates the ductility (20, 21). The different deformation mechanisms depend sensitively on the composition and microstructure. A key issue, in the presence of multiple deformation mechanisms, is at what stage each deformation mechanism emerges and how they contribute to the hardening behaviors. To this end, in situ observations are vitally important to determine the activation of different deformation mechanisms and understand their interactions during subsequent deformation. In particular, for deformation involving a phase transformation, an in situ study at the test temperature is the only way to observe the influence of phase transformation on deformation.

Up to now, deformation behavior at ultralow temperatures, where a myriad of deformation mechanisms is expected, has remained a largely unexplored territory, and no in situ study has been reported. Conventional experimental methods, such as transmission electron microscope (TEM), cannot meet the demand, because in situ experiments are difficult to implement at extremely low temperatures. In addition, the field of view by TEM is limited (22). In situ neutron diffraction is a powerful technique to identify deformation mechanisms (23, 24). In particular, because neutrons are highly penetrating, the neutron diffraction measurements are representative of the bulk of the sample rather than from the surface or local areas. While this method has been routinely used at room temperature, recent advances have made low-temperature measurements possible (25). Here, in situ neutron diffraction measurements were conducted on three representative fcc HEAs, i.e., CrMnFeCoNi, CrFeCoNi, and CrCoNi, which all showed a multistage deformation process. The prototypic deformation behavior is illustrated with the data of CrMnFeCoNi HEA, unless otherwise noted. The aim of this study was to establish the deformation pathway for HEAs at ultralow temperatures and to elucidate the interaction and competition between different deformation mechanisms.


Crystal structure and load curves at 15 K

The neutron detector banks at ±90° provided diffraction data along the loading and transverse directions, designated as LD and TD, respectively. The CrMnFeCoNi HEA has a single-phase fcc crystal structure at room temperature, which remained stable on cooling to 15 K and during deformation at 15 K to failure. The change in the lattice parameter during cooling is depicted in fig. S1. The evolution of the selected diffraction patterns is shown in Fig. 1A. No additional peaks were found, either during cooling or deformation, indicating the structural stability down to 15 K. The complete structure evolution during cooling and loading is presented in fig. S2. During the deformation, pronounced texture development was observed in addition to a shift in peak positions. Notably, (220)//LD disappeared almost entirely toward the end of deformation. Furthermore, the diffraction peaks also show significant broadening at large deformations.

Fig. 1 Crystal structure and deformation behavior of CrMnFeCoNi alloy at low temperatures.

(A) Selected diffraction patterns at room temperature and at 15 K during deformation, showing a clean single-phase fcc structure. (B) True stress-strain curves at room temperature, 140 K, and 15 K. The trend of temperature fluctuations (right-hand axis) due to serrations at 15 K is also superimposed. (C) Two serrations are shown with temperature (plotted in reverse scale for better comparison) to illustrate details of serrated deformation. (D) Plots of strain-hardening rate (SHR) at three temperatures. (E) Scanning electron microscopy image of the fractured sample at 15 K showing the 45° wedge.

The true stress-strain curves are presented in Fig. 1B. On cooling, the strength and ductility both increased. At 15 K, the alloy showed an extraordinarily high strength of ~2.5 GPa combined with ~62% in ductility. Serrated deformation, characterized by a sharp load drop and a temperature rise, can be seen beyond 21% strain. The amplitude of the load drop became increasingly larger with the deformation, approaching 570 MPa near the ultimate tensile strength. Note that this drop is even higher than the yield strength at room temperature and 140 K. The trend of load drop with deformation is shown in fig. S3. Figure 1C shows details of load drop and temperature rise in two consecutive serration events. The load drop was sharp, followed by a linear recovery over a time span of ~300 s. Continued loading triggers the next serration. It can be seen that the measured temperature started to increase after the stress dropped to the minimum value.

Figure 1D presents a plot of the strain-hardening rate at three temperatures. The initial rapid drop is attributed to dynamic recovery. As the temperature decreased, the HEA samples showed a higher and more stable strain hardening. The sample at 15 K exhibited the highest strain-hardening rate after the initial recovery, followed by the sample at 140 K, and the lowest hardening rate was at 295 K. At large strains, the strain-hardening rate for the 295 K sample continued to decrease with the deformation, while for the 15 K sample, the strain-hardening rate was stable until fracture.

Figure 1E shows the fractured sample at 15 K, with the fracture in the form of a wedge at nearly 45°. There was no significant necking even after 62% true strain, as commonly observed in conventional ductile alloys. By contrast, at 15 K, the ductile HEA sample was fractured by two shear bands at ~90° angles.

Lattice strain and texture evolution

The upper panels of Fig. 2 (A to C) show the evolution of lattice strains in LD at room temperature, 140 and 15 K, respectively. The lattice strain evolution for fcc HEAs at room temperature has been well established (16, 26) and has many of the characteristics of that due to dislocation slip (23, 24). For fcc HEAs, (220) is the stiffest orientation and (220)//LD grains carry the highest load, causing (220)//LD grains to reach the yield strength faster than other grains. As a result, the (220) elastic lattice strain ceased to increase after yield. The yielding of the (220)//LD grains transferred the load to the more compliant orientation, e.g., (200)//LD, where large elastic strains developed.

Fig. 2 Lattice strain and texture development during deformation.

Evolution of the lattice strain and the normalized integrated intensity for (hkl)//LD grains at (A) room temperature, (B) 140 K, and (C) 15 K for CrMnFeCoNi alloy. Below room temperature, the lattice strains for (111)//LD and (222)//LD split at large deformations. The same goes for (200)//LD and (400)//LD. The macroscopic yield strengths σy at different temperatures are indicated by vertical dashed lines.

Examination of the lattice strains measured with different orders of Miller indices yielded new information. At low temperatures, a split between ε111 and ε222 was observed (Fig. 2, B and C). In particular, ε111 < ε222, and the gap widened at larger deformation strains. Similarly, a gap was also seen between (200) and (400) reflections, but with the opposite sign, i.e., ε200 > ε400. The observed splits are signatures of deformation by stacking faults (SFs) (27).

The evolution of the integrated intensity for each reflection signifies the development of texture. The integrated intensity data, normalized by the initial stress-free values, are shown in the lower panels of Fig. 2 (A to C). In general, (111)//LD shows the largest increase in intensity, which means that the (111) grains are orientated toward the LD with increasing deformation. The intensity of (200)//LD also increased, but with a smaller magnitude, while the intensity of (220)//LD decreased. These experimental observations are signatures of the texture development due to dislocation slip in fcc metals and alloys (23, 24). At low temperatures, the texture developed slowly and plateaued due to the activation of additional deformation modes, which is discussed further below.

Deformation pathways at 15 K

Stacking fault probability (SFP), a measure of deformation faulting, is defined as the probability of finding a deformation fault between any two layers in the fcc stacking sequence (27). From the difference between ε222 and ε111, the SFP was calculated, following Warren (27). The evolution of SFP, together with the integrated intensity of (111)//LD in Fig. 3, provided insights into the transition of different deformation modes. As shown in Fig. 3, SFP is almost zero at room temperature. However, at 15 K, SFP developed early, when the stress > 1075 MPa, and increased to 27 × 10−3 before fracture.

Fig. 3 Deformation pathway of CrMnFeCoNi HEA at 15 K.

SFP and normalized integrated intensity of (111)//LD as a function of true stress are compared for room temperature and 15 K. Vertical dashed lines are drawn to pinpoint the changes in deformation behaviors in the 15 K dataset: (1) macroscopic yield point at σy,15K ~ 718 MPa; (2) start of SFs (σ ~ 1075 MPa); (3) first sign of serration (σ ~ 1270 MPa), followed by a rapid increase in SFP and change in the slope of texture development; and (4) massive serrations coincided with the saturation of texture (σ ~ 2000 MPa).

The lower panel of Fig. 3 compares the texture development for (111)//LD at room temperature and 15 K. As indicated earlier, the texture development is characteristic of dislocation slip, so the (111)//LD texture was used as a measure of the dislocation activity. After yielding, the (111) intensity exhibits a sharp increase at room temperature. At 15 K, the (111) intensity also increases but at a slower pace. At large deformations, three distinct regions can be identified from the slope of the intensity curve. The vertical dashed lines in Fig. 3 show the critical points where different deformation modes kick in. The sample yielded at σ ~ 718 MPa (σy,15K), which can be seen from the increase in the (111)//LD texture. At σ ~ 1075 MPa, SFs started; before this point, the deformation was mainly by dislocations. The transition of dislocation slip to SF activity suggests that the generation of full dislocations is replaced by Shockley partials (28). At σ ~ 1270 MPa, the texture development slows down, but the SFP started to increase rapidly; this is also where serrations began (after a true strain of ~21%; see Fig. 1B). Thus, it appears that the serrated deformation coincided with the rapid development of SFs and a slowing down in dislocation activity. The last point of interest is at σ ~ 2000 MPa, where the (111) intensity reached a plateau, a sign of saturation of dislocations and the deformation proceeding mainly by SFs and twinning. This is also the point where massive serrations started with a load drop exceeding 200 MPa (fig. S3).

TEM images were collected on fractured samples and are shown in Fig. 4 to corroborate the findings by neutron diffraction. From Fig. 4A, the primary deformation mechanism in CrMnFeCoNi alloy at room temperature is the dislocation slip. Sparse twins were also observed in some regions (Fig. 4B), which is a direct result of the relatively high stacking fault energy (SFE) at room temperature (29, 30). For the sample at 15 K, SFs were widely observed (Fig. 4C) and a very high density of twins is also evident from Fig. 4D. These observations, along with the high SFP obtained by neutron diffraction, indicate that the SFE decreased appreciably on cooling to 15 K. In addition, the formation of intersecting twins was observed, which resulted in a twinning network due to the activation of multiple twinning systems (Fig. 4E). The selected-area diffraction (SAD) pattern is presented as an inset of Fig. 4E, where the additional bright spots demonstrated clear evidence of twinning. The high-resolution TEM (HR-TEM) image in Fig. 4F also shows twins at 15 K.

Fig. 4 Electron microscopy results of fractured samples.

TEM images at room temperature showing the dislocation pileups (A), while the dark-field image also shows some sparse twins (B). The 15 K sample has a large number of SFs (C) and a very high density of closely spaced twins (D). An intersecting twinning network is formed at 15 K (E), while the SAD pattern in the inset shows the additional bright spots due to twinning. (F) HR-TEM taken from the 15 K sample shows deformation twins, with the twin boundaries marked.

Deformation diagram for fcc HEAs

While the low-temperature deformation behaviors are illustrated in detail with a quinary CrMnFeCoNi alloy, similar deformation behaviors were also observed in quaternary CrFeCoNi and ternary CrCoNi alloys. Taking these data together, a deformation diagram (Fig. 5) was constructed on the basis of the activation strains for SF and serration. The critical strain for SF activation in Fig. 5 shows a systematic decrease with decreasing temperature, demonstrating the lowering of SFE on cooling. When compared to the quinary alloy, the entire curve of SF activation for CrFeCoNi shifted downward, a manifestation of the composition dependence of the SFE. The trend continued for the ternary CrCoNi, where the SFs started at lower strain values. For the activation strain for serration, the values for the quinary and quaternary alloys are close to each other, while in the ternary alloy the serration starts early, at ~11% strain.

Fig. 5 Deformation diagram of fcc HEAs.

The critical strain at which the SFs and serrated deformation were activated. In quinary CrMnFeCoNi, the SFs are observed only at a strain value closer to fracture (36.8% in true strain) at 295 K, while at low temperatures, the SFs appear earlier, after ~22.2% at 140 K and ~15.8% at 15 K, respectively. The SFs form more easily in the quaternary CrFeCoNi alloy, where the SFs start at lower deformation strains for all three temperatures. The trend continues for the ternary CrCoNi alloy, where SFs begin to form even at an early stage of deformation at room temperature. The open symbols are the activation strain for serrations at 15 K (the data points are displaced to the left-hand side for clarity), which are substantially higher than the critical strains for SFs.


It is intriguing that HEAs exhibit extreme work hardening, while other alloys suffer from reduced ductility at ultralow temperatures. We believe that this is related to the chemistry and structure of HEAs. The fcc structure facilitates slip, so plastic deformation can easily start with dislocations. However, in contrast to conventional fcc alloys, the crystal lattice of HEAs is distorted due to the presence of multicomponents, which renders the magnitude of Burgers vector ill defined with segmented wiggled dislocation lines and offers resistance against dislocation motion (31). As a result, the propagation of dislocations is hindered, and this has led to the multiplication and hence high density of dislocations. Room temperature study shows that the dislocation density is already reaching ~1015 m−2 level (16), which could be even higher when the temperature is lowered.

Because of the favorable composition, the SFE in these HEAs is low (29, 30, 32), thereby allowing a transition from deformation by full dislocations to partials with SFs and twinning (28, 33). As the neutron diffraction data demonstrated, the critical strain for SF decreases nearly linearly with temperature. Thus, at ultralow temperatures, SF becomes a competitive deformation mode, with a high probability of SFs. Continuing deformation led to cross-slip and cross-twinning, ultimately causing plastic instability and hence the serration. Normally, serrations at low temperatures limit the ductility (20, 21). It is thus remarkable that the ductility of HEAs continued to increase in the presence of serrations. It is likely that in HEAs, the local disorder (e.g., stress field caused by size mismatch) (34) resulting from the multicomponent nature of HEAs successfully prevented the propagation of serration flow, thus leading to continued strengthening until major shear bands form fracturing the sample.

Quantitative analysis of the experimental data suggests that SFs played a critical role bridging dislocation slip and serration. Without SFs, the alloys would have gone from dislocation slip straight to serrations and would fail early with a reduced ductility, like other alloys at low temperature. Deformation by SF was made possible by two factors.

(i) The alloys remained in the fcc phase at 15 K, as the neutron diffraction data demonstrated for all samples. This was a surprising finding because the first-principle calculations predicted a stable hcp phase at low temperatures for a solid solution with a random distribution of constituent elements (19, 35). One possible explanation is that, as a recent theoretical study suggested, the SFE varies widely depending on the short-range chemical order and can even become positive depending on the local configurations (36). Another possibility, as Ma et al. (35) argued, is that the fcc phase may be stabilized by kinetics at low temperatures.

(ii) The SFE is sufficiently low so that SFs can start at an immediate step before serrations. The experimental values of SFE at room temperature for CrMnFeCoNi, CrFeCoNi, and CrCoNi are 30 ± 5 mJ/m2 (29), 27 ± 4 mJ/m2 (32), and 22 ± 4 mJ/m2 (30), respectively. These values can be compared with the SFE of Ni (the only fcc element in these alloys), which is 125 ± 5 mJ/m2 (37). At low temperatures, the SFEs were further suppressed, as is evident from the early activation of SFs at 140 and 15 K. Recent studies by first-principle calculations suggested that the temperature dependence of SFE is strongly linked to the entropy terms, including configurational entropy, vibrational entropy, and magnetic entropy (35, 38). Among these entropy contributions, the vibrational entropy due to high-frequency phonon modes has been highlighted as a term dominating the temperature dependence (38).

The observed SFs, twinning, and serrations at 15 K are the source of the increased strain-hardening rate at low temperatures. At room temperature, the deformation is dominated by dislocation activities, which gives rise to a characteristic strain-hardening rate. At low temperature, dislocation activities are still important, as the deformation starts with dislocation slip, but the role of additional deformation modes is also evident from the strain-hardening rate curves in Fig. 1D, where a higher and stable hardening rate was observed, especially for the sample at 15 K. The additional deformation mechanisms resulted in minimum necking (Fig. 1E) and hence more significant uniform elongation and higher strength.

As the in situ data in Fig. 3 demonstrated, deformation at 15 K started with dislocation slip. SFs and twinning became active later. This observation provided important clues regarding the microscopic origin of the increased work hardening at 15 K. It has been predicted by first-principle calculations that, for the HEAs studied here, the stable phase is hcp (19, 38, 39). An fcc to hcp phase transition has been reported in CrMnFeCoNi under pressure, demonstrating the metastable nature of the fcc phase at room temperature (40). The fcc phase, a stable phase at high temperature, becomes increasingly unstable as the temperature approaches absolute zero. It is important to note that the pathway from fcc to hcp is through SFs, as the hcp phase is formed by changing the stacking sequence from A-B-C-A-B-C in fcc to A-B-A-B in hcp (28). While the fcc phase starts to deform by dislocation slip, SFs are activated at large strains when the fcc becomes unstable. Deformation by SFs, and subsequently twinning, creates an additional boundary. It has been shown by Zhu et al. (41) that there are significant energy barriers for a dislocation to cross the SF or twin boundary. The additional stress to overcome the barrier, due to the interaction between dislocations and SFs and twin boundaries, impedes the propagation of dislocations and gives rise to the extreme work hardening, leading to the unusually large ductility at 15 K. To this end, the in situ neutron diffraction experiments underscored the underlying link between the preferred deformation mechanism and lattice instability.

Recently, Sun et al. (42) conducted an ex situ TEM study on ultrafine-grained CrMnFeCoNi HEA at 77 K, where HR-TEM images showed deformation twins and some evidence of SFs. Our in situ experiments at 15 K, a lower temperature, demonstrated that SFs became abundant and an additional deformation mechanism, serration, emerged. Unlike dislocations, SFs, and twinning that occur across the entire test specimen, serration is a form of inhomogeneous deformation localized in particular regions. Serrations are commonly observed in metallic glasses, where other deformation mechanisms associated with crystalline materials are not available (43). Serration behavior has also been reported in Al-containing HEAs (3, 44). However, the serration events in Al-containing HEAs occur at high temperatures, which is qualitatively different from our observations at 15 K. High-temperature serrations are associated with dynamic strain aging, arising from the interaction between diffusing solutes (in this case, Al atoms) and dislocations (44). The magnitude of the load drop for serration at high temperatures, also known as the Portevin–Le Chatelier (PLC) effect, is extremely small compared to that observed for HEAs at 15 K. At ultralow temperatures, where the diffusion is frozen, the serration occurs through different mechanisms (20, 21).

There are presently two schools of thought regarding the mechanism of serration at low temperatures: (i) mechanical origin due to the instability of plastic deformation (e.g., dislocation pileups) or (ii) a manifestation of thermodynamic instability where the specific heat and thermal conductivity vanish (20, 21, 45). In our experiments, as can be seen from Fig. 1C, the temperature rise follows the load drop, suggesting that the temperature rise is a consequence of the load drop. However, it is cautioned that, in this experiment, the thermocouples were located at the sample grips, not on the sample. Thus, there could be a delay in the temperature sensor. For the same reason, the local temperature rise in the deformed region could also be much higher. It is possible then that locally, the HEAs were heated to a very high temperature state by the adiabatic heat release, where some of the deformation was annealed during the recovery (~50 s from the temperature data), before the region of serration was cooled back to ~15 K.

In summary, in situ neutron diffraction measurements revealed that the extreme work hardening observed in HEAs at 15 K is the result of the cooperation of multiple deformation mechanisms at low temperatures. Like at room temperature, the deformation at 15 K started with dislocation slip, but SFs and twinning were activated at larger strains before transitioning to serrations. The activation of these additional deformation mechanisms is the source of the extreme work hardening and the exceedingly large ductility observed at 15 K. The intermediate step, SF, is critically important and highlights the influence of entropy and lattice instability on deformation and mechanical behaviors. As our study shows, the unique combination of complex chemistry and simple crystal structure may hold the key to the design of new materials for applications at low temperatures, where the competitive mechanisms can be deployed in synergy to produce materials with high strength and superior ductility.


Sample preparation

Ingots of Cr20Mn20Fe20Co20Ni20, Cr25Fe25Co25Ni25, and Cr33.3Co33.3Ni33.3 (atomic %) alloys were fabricated by arc melting a mixture of pure elements (purity > 99.9%) in a Ti-gettered high-purity argon atmosphere. To ensure composition homogeneity, the ingots were flipped over and remelted eight times. Last, they were drop-cast into rectangular Cu-mold of dimensions 15 mm × 15 mm × 65 mm. The ingots were homogenized by annealing at 1200°C for 24 hours under inert argon atmosphere followed by quenching in iced water. Cold rolling was conducted to ~80% reduction in thickness (15 to ~3.1 mm) to break down the dendritic structure, followed by recrystallization treatment at 800°C for 1 hour and iced-water quenching. The equiaxed grains were obtained with an average grain size of ~5 μm. Rectangular dog-bone shaped tensile specimens with gauge dimensions of 3 mm × 4 mm × 25 mm (thickness × width × length) were prepared by electrical discharge machining.

In situ neutron diffraction measurements

The in situ neutron diffraction experiments under tensile loading were conducted at the Engineering Materials Diffractometer TAKUMI (beamline 19) at the Materials and Life Science Experimental Facility (MLF) of Japan Proton Accelerator Research Complex (J-PARC) using a cryogenic load frame. The sample, 45° to the incident neutron beam, was mounted inside the cryogenic loading cell chamber, which was sealed and evacuated to <10−3 kPa before the cooling was started. A Gifford-McMahon (GM) cooler installed in the cryogenic load frame was used to cool the grips and the specimen. The GM cooler had two stages of 50 and 4.2 K, respectively. The 4.2 K stage wound with a heater wire was connected to both grips using copper wires. The calibrated Cernox thermosensors were mounted on the Cu-Be grips to monitor and control the temperature. The details of this cryogenic load frame for the neutron diffractometer have been reported in (46). To monitor the sample temperature, we attached two thermocouples on both sides of sample mount. Note that the thermocouple could not be fixed on the sample gauge section to avoid the incident neutron beam. The sample was cooled down and held for ~1 hour at a stable temperature of 15 K (vacuum ~ 10−5 kPa) before the loading was commenced. Samples were loaded in 10 equal steps in the elastic region and held for 10 min at each step to obtain good statistics. Continuous loading was carried out in the plastic regime at a strain rate of 0.04 mm/min (2.67 × 10−5 s−1). The time-of-flight (TOF) neutron diffraction data were acquired simultaneously for the loading and transverse directions by two detector banks, at ±90° to the incident neutron beam. The data at 140 K and room temperature were also obtained in a similar manner.

Data analysis

Rietveld refinement for multiple peak fitting of the neutron diffraction data was carried out using Z-Rietveld software (47). Lattice strain εhkl was calculated from the TOF data by the following methodεhkl=(dhkldhkl°)dhkl°=(TOFhklTOFhkl°)TOFhkl°where dhkl is the lattice spacing of hkl planes, dhkl° is the value for the stress-free sample, and TOFhkl is the Bragg position of the hkl reflection in the TOF data.

SFP was calculated from the lattice strain by following equation 13.67 provided by Warren (27)εhkl=3 α4πho2(u+b)b(±)Lowhere εhkl is the experimentally measured lattice strain for the specific hkl reflection, α is the SFP, ho2=h2+k2+l2 for the cubic system, (u + b) represents the total number of unbroadened and broadened components, Lo = 3N ± 1 for components broadened by faulting, and Lo = 3N for components not broadened by faulting. To calculate SFP from the lattice strain of (111) and (222), a simplified equation was obtained by inputting the 1ho2(u+b)b(±)Lo value of +14 for (111) reflection and 18 for (222) reflection (27). It followsSFP=32π33(ε222ε111)where ε222 and ε111 are the experimentally measured lattice strain values with applied stress for (222) and (111) reflections, respectively. Adjacent reflections were used from the lattice strain plots for the SFP calculation, as their difference due to faulting can be observed for the in situ neutron diffraction deformation data (Fig. 2). The measured lattice strain for (111) and (222) reflections should be the same as they are crystallographically equivalent. However, SFs led to the difference in their lattice strain due to deformation faulting (27). The SFP behavior was also verified from the (200) and (400) pair [SFP=16π33(ε400ε200)]; the effect was opposite as illustrated in Fig. 2 (B and C).

Transmission electron microscopy

For microstructural examination of fractured samples, a Tecnai F30 TEM operated at 300 kV acceleration voltage was used. TEM specimens with a diameter of 3 mm were cut from the deformed samples and then ground to around 50 μm by 2000-grit SiC paper under the press of a rubber plug, followed by twin-jet electropolishing using a mixed solution of HClO4:C2H6O (volume ratio, 1:9) with a direct voltage of 30 V and a current of 50 mA at a temperature of 243 K.


Supplementary material for this article is available at

Fig. S1. Lattice parameter change during cooling.

Fig. S2. Structure evolution of CrMnFeCoNi HEA with time.

Fig. S3. The magnitude of load drops with deformation.

This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.


Acknowledgments: Funding: This work was supported by grants from the Croucher Foundation (project no. CityU 9500034), the Research Grants Council of the Hong Kong Special Administrative Region (project no. CityU 11215917), and the National Natural Science Foundation of China (project no. 51571170). X.-L.W. acknowledges the support by the Shenzhen Science and Technology Innovation Committee (grant no. JCYJ20170413140446951) and the Ministry of Science and Technology of China (grant no. 2016YFA0401501). The in situ neutron diffraction experiments were carried out at the TAKUMI beamline of the MLF, J-PARC, under proposal numbers 2016B0254 and 2017B0142. Author contributions: M.N., H.He., S.H., T.K., and X.-L.W. designed and performed the neutron diffraction experiments using TAKUMI at J-PARC, with samples prepared by M.N., F.Z., and S.L. M.N. and H.Hu. carried out TEM examination. M.N. and H.He. analyzed the data. M.N. and X.-L.W. drafted the manuscript. All authors reviewed the results and provided input to this paper. X.-L.W. conceived and supervised the overall project. Competing interests: The authors declare that they have no competing interests. Data and materials 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 authors.
View Abstract

Stay Connected to Science Advances

Navigate This Article