Scholarly article on topic 'Dislocation mechanisms and 3D twin architectures generate exceptional strength-ductility-toughness combination in CrCoNi medium-entropy alloy'

Dislocation mechanisms and 3D twin architectures generate exceptional strength-ductility-toughness combination in CrCoNi medium-entropy alloy Academic research paper on "Materials engineering"

Share paper
Academic journal
Nature Communications
OECD Field of science

Academic research paper on topic "Dislocation mechanisms and 3D twin architectures generate exceptional strength-ductility-toughness combination in CrCoNi medium-entropy alloy"



Received 14 Jun 2016 | Accepted 21 Dec 2016 | Published 20 Feb 2017^^HDOI:iai038/ncomms14390l OPEN

Dislocation mechanisms and 3D twin architectures generate exceptional strength-ductility-toughness combination in CrCoNi medium-entropy alloy

Zijiao Zhang1, Hongwei Sheng2, Zhangjie Wang3, Bernd Gludovatz4, Ze Zhang1, Easo P. George5, Qian Yu1, Scott X. Mao1'6 & Robert O. Ritchie4'7

Combinations of high strength and ductility are hard to attain in metals. Exceptions include materials exhibiting twinning-induced plasticity. To understand how the strength-ductility trade-off can be defeated, we apply in situ, and aberration-corrected scanning, transmission electron microscopy to examine deformation mechanisms in the medium-entropy alloy CrCoNi that exhibits one of the highest combinations of strength, ductility and toughness on record. Ab initio modelling suggests that it has negative stacking-fault energy at 0K and high propensity for twinning. With deformation we find that a three-dimensional (3D) hierarchical twin network forms from the activation of three twinning systems. This serves a dual function: conventional twin-boundary (TB) strengthening from blockage of dislocations impinging on TBs, coupled with the 3D twin network which offers pathways for dislocation glide along, and cross-slip between, intersecting TB-matrix interfaces. The stable twin architecture is not disrupted by interfacial dislocation glide, serving as a continuous source of strength, ductility and toughness.

1 Center of Electron Microscopy & State Key Laboratory of Silicon Materials, Department of Materials Science & Engineering, Zhejiang University,

Hangzhou 310027, China. 2 Department of Physics and Astronomy, George Mason University, Fairfax, Virginia 22030, USA. 3 Department of Materials

Science and Engineering, Xi'an Jiaotong University, Xi'an 710049, China. 4 Materials Sciences Division, Lawrence Berkeley National Laboratory, Berkeley,

California 94720, USA. 5 Materials Sciences & Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA. 6 Department of

Mechanical Engineering & Materials Science, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, USA. 7 Department of Materials Science & Engineering,

University of California, Berkeley, California 94720, USA. Correspondence and requests for materials should be addressed to Q.Y. (email:

or to S.X.M. (email: or to R.O.R. (email:

Strength and ductility (and hence toughness) are key mechanical properties of structural materials1-14, although these properties are often mutually exclusive15. In crystalline metallic materials, strength and ductility depend on the presence of crystal defects and how they move under mechanical loading. Dislocations serve as the prime carrier of plasticity in metals and their motion to create plastic deformation is strongly influenced by interactions with other dislocations as well as defects such as solute atoms14,16-20, stacking faults21,22, grain boundaries (GBs)23-25 and reinforcement particles or precipitates26,27. In general, strong interactions (obstacles) strengthen a material although this often limits its ability to plastically deform. Twinning, however, appears to be a deformation mechanism capable of defeating this 'conflict' between strength and ductility1-6,8,9,16,28,29. The presence of twins usually serves to impede dislocation motion and induce strengthening, but multiple twinning systems can also enhance ductility. This is apparent in twinning-induced plasticity (TWIP) steels, which are known to form multiple types of twins that result in high strength with substantial uniform ductility6,10-12. Another notable example includes certain equiatomic multielement alloys, termed high-entropy alloys30 (HEAs), based on the CrMnFeCoNi composition that have a single-phase face-centered cubic structure (fcc) structure31,32 with relatively low-stacking-fault energies (SFEs), comparable with 304 stainless steels33 and relatively large separations between the Shockley partials34. These alloys can display excellent strength9,35, ductility9,35 and toughness1 at room temperature, properties which are further enhanced or maintained at cryogenic temperatures where deformation nano-twinning becomes a more prominent mode of deformation. Despite our recent observations of various toughening mechanisms occurring on the nanometre length-scale in the vicinity of a crack tip in such an alloy8, the quantitative strengthening effect from deformation twinning and the nature of twinning-induced plasticity remain somewhat uncertain.

In this work, we show that the mechanical behaviour of the nominally equiatomic alloy CrCoNi, a medium-entropy alloy (MEA) with a single-phase, fcc structure16,36, and even better mechanical properties than the CrMnFeCoNi high-entropy alloy—specifically, a high strength (~1GPa) and ductility (~60% uniform elongation) combined with exceptional fracture toughness (KIc > 200 MPa^m)8,16—can be related to its capacity for twinning-controlled deformation. Our ab initio calculations based on the density-functional theory (DFT) suggest that this alloy has a negative SFE at 0 K, similar to what has been reported for some other systems37-41 and high propensity for twinning when compared with pure fcc metals42. We seek to demonstrate that a 3D twin network is formed within the grains of this fcc CrCoNi MEA and that the interactions of dislocations with these twins can lead to a simultaneous increase of strength and ductility, that is, to identify the origin of its excellent mechanical properties.


Numerical simulations. DFT-based ab initio calculations were performed to compute the formation energy of the equiatomic CrCoNi alloy with different crystal structures as well as their planar fault energies. In previous theoretical work33,43,44, the formation energies of medium- and high-entropy alloys were treated with the coherent potential approximation (CPA)44,45 or the special quasi-random structure (SQS)46 method in conjunction with ab initio treatments. In this study, the formation energies of CrCoNi crystals were explicitly assessed with the 'multiple randomly populated supercell' approach. The

configurational dependence of the formation energy of CrCoNi MEAs is shown in Supplementary Fig. 1. In our calculations, atoms of the different elements comprising the MEA were randomly assigned to the lattice points of the corresponding crystal structure, followed by geometrical optimization to obtain the ground state configuration of each structure. Although computationally demanding, this approach gives us an advantage by allowing direct evaluation of the stacking-fault energies and their formation energy barriers.

Supplementary Table 1 lists the formation energies of hypothetical face-centered cubic (fcc), hexagonal close-packed (hcp), double hcp (dhcp) and 9R structures for the CrCoNi alloy. Note that their respective stacking sequences are: fcc— ...ABC ABC..., hcp— ...AB AB..., dhcp— „.ABAC ... and 9R—... ABA BCB CAC...

It was found that the formation energy of hcp CrCoNi at 0 K is slightly lower than that of fcc CrCoNi, which implies that the hcp structure can be a competing phase during the synthesis of CrCoNi (Supplementary Fig. 1). The energy difference at 0 K between fcc and hcp CrCoNi was found to be ~7meV per atom (that is, 0.68kJmol_ 1). According to our free-energy calculations based on the quasi-harmonic approximation (Supplementary Fig. 2), the fcc structure gains a free-energy advantage at higher temperatures (for example, 500 K), due mainly to the contribution from lattice vibration. Similar results (namely that hcp is stable at low temperatures and fcc at higher temperatures) were obtained recently for the 5-element HEA, CrMnFeCoNi, when, in addition to configurational entropy, other entropic contributions including electronic, vibrational and magnetic were included in the total energy calculations44. Lattice distortions arising from size mismatches of the constituent atoms in the CrCoNi MEA with the fcc or hcp structure were also examined from the atomic pair-distribution functions (PDFs) of CrCoNi crystals; these PDFs of fcc and hcp CrCoNi at 0 K are shown in Supplementary Fig. 3. Further description of the analyses is provided in Supplementary Note 1.

The stacking faults in fcc CrCoNi were analysed with ab initio modelling; results are shown in Fig. 1. A slab of CrCoNi with the fcc structure (360 atoms) was first created by sequentially stacking 12 close-packed (111) atomic planes. The slabs have two free surfaces. In the simulation cell, in the direction perpendicular to the slab (z direction), a vacuum of thickness 1.6 nm is inserted between slabs to separate each slab from the adjoining ones in the next periodic cell (to minimize the interactions of the slabs). The atomic configurations of the slab were geometrically relaxed in the three directions. Next, an intrinsic stacking fault, of energy yisf, was created by shifting the top six layers along the (112) direction by the Burgers vector of the Shockley partial bs=a/g(112), where a is the lattice parameter; this resulted in a stacking sequence: ABCABC|BCABCA, shown in Fig. 1b. The stacking-fault energy is obtained by the energy difference of the two optimized structures normalized by the area of the stacking fault A.

Subsequent displacement of the top five layers in Fig. 1b along the (112) direction results in a two-layer stacking fault (Fig. 1c), which is called an extrinsic stacking fault (ABCABC|B|ABCAB). Further shifting of the fault in Fig. 1c leads to a three-layer fault, shown in Fig. 1d and so on. As such, continued shifting of the crystal will create a twinned structure with a twin boundary. The calculated energies of the intrinsic stacking fault, gisf, and twin boundary, gTB, are given in Table 1. The energy gus is that of the unstable stacking fault, and gut is the nucleation energy barrier for the formation of twin boundaries (Fig. 1f).

Alternatively, stacking-fault energies can be calculated through the use of a parametrized method, that is, the axial interaction model (AIM)47. Using AIM, the energy of an intrinsic stacking

• t «99

P • ^ ^

I • • • •

• I • • • • •

-9- ■•-•- -,

• • • • •

a V3 V2

A •• • C

t • w • • •

• • • • • i

: • • • « • • • • •


Intrinsic stacking fault

B • • • •• C

A • r • B •

C • • • • • A >

A A »1

C •>• C «

B B •

A A »

C / • • • • C i

B • • • • B ;

A • • • • • A !

Extrinsic stacking i

f • • # • • • • •

• • # V

• • •

Three-layer twin

Cr Co Ni

T 400 _ 350 rE 300 | 250 ^ 200 CÛ 150 § 100 = 50

Tus A T'ut A

A \ / \

y ' \ Tisf *. Tesf „ j \ V

0.0 0.5 1.0 1.5 2.0

Ux/\ 1/6 <112>l

2.5 3.0

Figure 1 | Stacking faults and their ab initio calculated energies for the CrCoNi alloy. (a-d) show the atomic configurations of the original fcc structure, intrinsic stacking fault, two-layer extrinsic fault and three-layer twin, respectively. (e) Top view image showing the close-packed (111) plane. (f) Energy barriers along the displacement pathway (112) direction in e that result in the two types of stacking fault and a twin boundary. The smallest displacement along the pathway 6 (112) is given by the magnitude of the Burgers vector of a Shockley partial dislocation, bs. Here, gisf and gesf represent the stacking-fault energies of intrinsic (b) and extrinsic (c) stacking faults, respectively. gus is the unstable stacking-fault energy, and gut denotes the energy barrier for the formation of the initial twin boundary. After gut is overcome, the subsequent energy barriers (for example, g'ut) for creating multi-layer twins are smaller than gut. The energy barriers along the displacement pathway were obtained with the nudged-elastic band method in our ab initio calculations.

Table 1 | Ab initio calculations of various fault energies in CrCoNi MEA at 0 K.


Fault energy, E (mJ m - 2) - 24 264 -17 310

fault can be expressed in terms of the energies (E) of the various structures and the fault areas (A) as:

gisf = (Ehcp + 2Edhcp - 3Efcc)/A,

gesf = 4(Edhcp - Efcc)/A, (1)

gTB = 2(Edhcp - Efcc)/A.

Substituting the values of the formation energies of the hcp, fcc and dhcp structures from Supplementary Table 1, we obtain gisf =- 23mJm- 2, gesf =- 20mJm- 2 and gTB = - 10mJm - 2, respectively, which are generally consistent with the results of our explicit calculations provided in Table 1.

An energy barrier must be overcome when an intrinsic stacking fault is created by displacing the upper half of the crystal, as demonstrated in Fig. 1f. This energy barrier corresponds to the unstable stacking-fault energy of fcc CrCoNi, gus. On further displacement of the crystal with a pre-existing fault (as shown in Fig. 1b), another energy barrier is crossed, namely gut, the energy barrier for twin-boundary formation. These energy barriers lie on the minimum energy paths between the stacking faults, and have profound implications for the deformation mechanism of the crystal.

In this work, we used the nudged-elastic band (NEB) method48 in the DFT calculation to obtain the unstable stacking-fault energies along the deformation paths. During the NEB simulation, the atoms are relaxed to the minimum energy paths. It should be mentioned that, without the atomic relaxation process, the unstable fault energy calculated for the rigidly shifted crystal is about 1.5-2 times higher than the correct values shown in Fig. 1f.

Having obtained the stacking-fault energy values in Fig. 1f, we are able to evaluate the propensity of twinning ('twinnability') in the CrCoNi crystal. On the basis of Tadmor and Bernstein's work on fcc metals42,49, an empirical expression for twinnability, tta, is:

tta = [1.136 - 0.151 f J y/^us/gut. (2)

Using this equation, the dimensionless twinnability for the present CrCoNi fcc alloy is tta — 1.07 ± 0.02, which is higher than that of all pure fcc metals42 (for example, the numerically derived tta for Ag is ~ 1.05), and comparable to that of Fe-Cr-Ni TWIP steels50 (1.08-1.12, based on equation (2) and the ab initio data from ref. 50). The high propensity for twinning of fcc CrCoNi appears to be a factor in the excellent strength, ductility and toughness of this MEA.

Nanostructural characterization. This calculated twinnability is consistent with the twinning observed in the CrCoNi alloy after plastic deformation8. To examine its specific atomic-scale mechanisms of deformation, we performed advanced structure characterization utilizing quantitative in situ transmission

Figure 2 | TEM of twin structures in the CrCoNi alloy. (a) Bright-field TEM image showing the hierarchical twinning architecture in a grain of the CrCoNi alloy. A grain boundary is marked by the yellow line near the top-left corner, and multiple twinning systems are labelled. Scale bar, 1 mm. (b) Low-magnification bright-field TEM image showing dislocation arrays on the twin boundary. Scale bar, 500 nm. (c) SAED pattern along <110> beam direction from the region on the CTB circled in blue in a showing extra spots which belong to the twin. (d) HAADF STEM image showing the structure of a CTB and an ITB which contains a 9R structure. This image was taken from an intersection of CTB and ITB of twin 1 in a. Scale bar, 500 pm.

electron microscopy (TEM) and aberration-corrected scanning transmission electron microscopy (STEM). We also performed quantitative in situ TEM compression tests to investigate the strengthening effect of a twin boundary in the CrCoNi alloy. By compressing micro-pillars with and without a twin boundary, we were able to quantify the strengthening effect from blocked dislocations at a twin boundary (TB). The dislocation-twin and twin-twin interactions during plastic deformation were systematically studied in situ in the TEM to reveal the influence of the 3D twin architecture on the plastic deformation of the material. Remarkably, while the twin boundaries contributed to strengthening by acting as barriers to dislocation motion, the interconnected twin boundaries in the 3D twin network can also generate significant ductility by offering multiple pathways for dislocation motion along the twin boundaries.

Post-mortem TEM investigations of the CrCoNi samples revealed that large numbers of deformation twinning events occurred in samples deformed in uniaxial tension; three twinning systems were activated within individual coarse grains, forming a 3D twin architecture in which multiple twins with different thicknesses form and intersect each other. The twins typically had thicknesses ranging from B100 nm to several micrometres (examples are shown in Fig. 2a), and were present within most grains that we examined. Two twin variants, comprising a compound twin structure that is common in fcc crystals, were observed to intersect each other; these are labelled as 'Twin 1' and 'Twin 2' in Fig. 2a. Selected area electron diffraction (SAED) patterns taken along the [110] zone axis (Fig. 2c) reveal separate sets of spots, belonging to the matrix and the twin, respectively. Two types of twin boundaries were observed: {111} coherent twin boundaries (CTBs) and {112} incoherent twin boundaries (ITBs). Dislocations interact with these TBs in various ways. For example, dislocations pile-up at CTBs when they glide on a {111} plane that is not parallel to the CTB, as shown in Fig. 2a. Moreover, dislocation

arrays can also form near the TBs, either at a CTB or the intersections of CTBs and ITBs (Fig. 2b). Corresponding atomic structures of a dislocation at the intersection of a CTB {111} and ITB {112} are shown in Fig. 2d, which is a typical high-angle annular dark field (HAADF) STEM image obtained in the aberration-corrected TEAM0.5 TEM microscope operating at 200 kV. A 9R structure composed of a set of Shockley partial dislocations with a repetitive sequence b2:b1:b3 (refs 51-54) on every (111) plane can be seen in the {112} ITB. The twin architecture is presumably formed due to twin-twin interactions. A smaller twin (for example, Twin 3 in Fig. 2) can be generated within large parent twins by the ITB of other twins (Twin 1) advancing and hitting the CTB of a thick twin (Twin 2) during plastic deformation, as shown in Fig. 2.

In situ mechanical testing. A high density of twins is expected to have a marked strengthening effect as twin boundaries can act as effective barriers to dislocation motion. Using micro-pillar compression experiments, we attempted to quantify the contribution of an individual twin boundary to the strength of this CrCoNi alloy. We used B300-nm diameter, focused-ion beam (FIB) machined micro-pillars, which were tested in a JEOL 3010 TEM at 300 kV, either with or without a {111} twin boundary (and no other boundary), as shown respectively in the dark field images in Fig. 3a,b and in the Supplementary Movie 1. The SAED pattern in the inset of Fig. 3a shows two sets of spots, confirming the existence of the twin. Engineering stress-displacement curves from the compression of these pillars (Fig. 3 c) clearly show the much higher strength of the pillar containing an initial twin. Specifically, the maximum strength of the pillar with a pre-existing twin exceeded 2.5 GPa, that is, about 1 GPa higher than that without a twin. As the sizes of the pillars were kept the same to minimize possible size effects, such an elevation in strength can be plausibly attributed to the

Displacement (nm)

Figure 3 | In situ compression of CrCoNi micro-pillars with and without a twin boundary. (a) Low-magnification dark field TEM image showing the structure of the pillar containing a twin; the inset in a shows the SAED pattern of the pillar. The spots from the matrix and twin can be distinguished. The twin boundary is marked by a white dashed line. The zone axis is [110]. Scale bar, 200 nm. (b) Low-magnification dark field TEM image showing the pillar with no twin (g = [220]). Scale bar, 200 nm. (c) Engineering stress-displacement curves of compression tests on the two pillars shown in a and b. The pillar containing a twin boundary is 67% stronger than the pillar without the boundary.

presence of the twin boundary. However, as these nano-pillar tests only investigate limited volumes, in situ TEM straining tests were also performed to systematically investigate the dislocation-twin and twin-twin interactions, with particular focus on the formation of the 3D twin architecture and how this architecture influences the dislocation activity.

The in situ straining experiments enabled direct observations of the formation of the twinning architecture and the interactions between dislocations and twin boundaries, as shown in Fig. 4 and Supplementary Movies 2-5. During straining, multiple types of deformation twins were quickly activated, forming a hierarchical architecture consisting of twins with different thicknesses. Figure 4a shows the dynamic formation of two twins, labelled as 'Twin 1' and 'Twin 2', and their interaction to form a twin junction. It was observed that the twins propagated very fast until they hit the boundary of other twins. During further straining, the hierarchical twin architecture remains stable with no noticeable thickening or de-twinning. Activity of screw dislocations was also observed at the TB. A dislocation, originally gliding on the boundary between 'Twin 2' and the matrix, can be seen to impinge at the intersection of 'Twin 1' and 'Twin 2', before subsequently cross-slipping onto the boundary between 'Twin 1' and the matrix, on another {111} plane. From the real-time observations, the twin boundaries clearly act as strong barriers to dislocation motion, preventing dislocations from moving across them and resulting in the pile-up of

dislocations (Supplementary Movie 3), as has also been reported previously55-59. Although impeding dislocation motion usually compromises ductility, the generation of multiple types of deformation twins contributes to the homogeneous plastic deformation. More importantly, the intersected boundaries of the hierarchical twins in the CrCoNi MEA were able to serve as fast channels for dislocation movement. Specifically, dislocations can cross-slip from one twin boundary to another using the 'overpass' constructed by the intersected twin boundaries, resulting in a more homogeneous distribution of dislocation activity in three dimensions. Although we did not capture the detailed process, sessile dislocation(s) should be involved in the cross-slip phenomenon of those imperfect dislocations at the intersection of the twin planes. We observed that, during further straining, the majority of the deformation inside individual grains occurred through gliding of dislocations along the twin boundaries in the 3D architecture. The dislocations moved fast along the boundaries, which is consistent with previous reports that demonstrated that the dislocation motion on a twin boundary would be faster than that in the matrix since the impediment from the image force is null if the Burgers vector of dislocations is parallel to the twin boundary41'60'61. Figure 4b shows a side view of a twin boundary in which the TBs are inclined and intersect both the top and bottom surfaces of the TEM film. Although cross-slip of partial dislocations62,63 at twin-twin junctions could not be fully observed due to the limited time-resolution of in situ TEM, in Fig. 4b we show the glissile boundary dislocations with Burgers vector parallel to <112> gliding on the twin boundary; however, the twin boundaries remained stable as these dislocations moved along them, that is, there was no evidence of either twin growth or de-twinning occurring during such deformation (Fig. 4b and Supplementary Movies 4 and 5). As such, a high ductility was achieved.


On the basis of the in situ straining of the CrCoNi, we have found that a three-dimensional hierarchical twin network is established within individual grains in this MEA, associated with its very low (negative) stacking-fault energy calculated at 0K; this network presents substantial barriers for dislocation motion, and contributes to its high strength and significant strain hardening. However, at the same time, the network provides multiple pathways for the easy motion of dislocations, which permits the simultaneous generation of significant plastic deformation. Importantly, the formation of the 3D twin network is achieved within relatively large grains (~5-50 mm in size8), and is found to be stable against de-twinning as deformation continues with the twin-boundary interfaces serving as pathways for dislocation slip. This observation is significantly different from that previously reported in twinning-deformation dominated materials3'4'6'7'64. As a result, high strength coupled with high ductility and continuous strain hardening is achieved in this alloy, which presents the perfect ingredients for its exceptionally high fracture toughness8. As lower (cryogenic) temperatures serve to elevate the strength and hence are expected to further promote twinning activity, the damage tolerance (strength with toughness) of this alloy is destined to be enhanced at such low temperatures, which is exactly what has been observed experimentally8. This dual effect on dislocation activity is a consequence of the three-dimensional hierarchical twinning architecture that is generated with plastic strain, which in turn results from the high twinnability of this particular alloy.

We note that 3D twinning networks resulting from the activation of three twinning systems may also form within

t = 0s b^r. * Twin 1 t = 2s Twin 2 Dislocations glide from twin2's boundary to twin1 Twin 1 I — tuaffl

t=0s •f ; . Side view of a CTB t=4s \ / /■ Glissile boundary dislocations with Burgers vector parallel to <112> B'

Figure 4 | In situ imaging shows the dislocation and twin sequences during deformation. (a) TEM image sequence showing the formation of the twinning architecture and the dynamic process of dislocations gliding from one twin boundary to another. Scale bar, 200 nm. (b) TEM images showing the glide of paired partial dislocations on a CTB. The inset on the right shows the change of stacking at the twin boundary from fcc to hcp due to the glide of a partial dislocation. Scale bar, 200 nm.

individual grains in TWIP steels, specifically under complex loading conditions such as torsion. The formation of a high density of twin lamellae serves to refine the grains, thereby hindering dislocation motion and generating strain hardening; ductility is considered to originate from the formation of the deformation twins. In comparison, the density of twins in CrCoNi is much lower, but the intersected twin boundaries of different twins act to construct an overpass, which enables fast motion of dislocations on the boundaries and cross-slip from one boundary to another, leading to more homogeneous plastic deformation even within individual grains.

We believe that documentation of the enhanced twinning-induced toughness of this MEA may serve not only to improve our understanding of the extraordinary structural behaviour of specific high-/MEAs but also to present new directions for the future design of metallic alloys with unprecedented combinations of mechanical properties.

the corresponding crystal structure, followed by geometrical optimization to obtain the ground state configuration of each structure. The formation energy of each structure corresponds to a statistical mean of 80 such calculations. The Purdew-Wang-type generalized gradient approximation (GGA)70 was employed for the exchange-correlation functional. Our calculations were mainly performed on non-magnetic structures (that is, non-spin-polarized), but the effect of magnetization was also analysed with the results presented in Supplementary Fig. 4. The unstable stacking-fault energy and the energy barrier for twinning were analysed with the nudged-elastic band method48 implemented in the VASP package. To obtain the formation energy, the initially created lattice structure was subjected to geometrical optimization with high-precision calculations (the energy cutoff used in this work was 337.0 eV) on the r point of the Brillouin zone. Geometrical optimization was fully achieved when the total energy converged to within 10 _ 4eV between successive steps. Having optimized the volume while allowing the atoms to relax to their ground state, the formation energy (cohesive energy) was derived by subtracting the sum of single-atom energies from the total energy.

Data availability. The data that support the findings of this study are available from the corresponding authors upon request.


Materials preparation and microstructural characterization. Samples for this study were extracted from a previously produced nominally equiatomic CrCoNi MEA whose microstructure and mechanical properties were reported in a recent paper8 where details of its processing and mechanical characterization can be found. Atomic structures were investigated using the aberration-corrected TEAM0.5 transmission electron microscope (operating at 200 kV), housed at the National Center for Electron Microscopy at the Lawrence Berkeley National Laboratory (LBNL), and the in situ compression tests were performed using a Hysitron PI95 nanoindenter in a JEOL 3010 microscope at 300 KV. The nanopillars for the in situ compression tests were produced using focused-ion beam techniques; details of sample preparation and in situ compression have been described in previous studies65,66. The in situ TEM tensile tests were conducted at room temperature using a Gatan model 654 single-tilt straining holder in an FEI Tecnai G2 F20 TEM operating at 200 kV. Roughly 12 samples, thinned by jet polishing and well attached to the substrate, were selected for in situ tensile straining and detailed TEM investigation as described in a previous paper22. Time-resolved TEM and HRTEM images of the regions of interest were recorded with a Gatan CCD camera at 10 frames per second.

Computational methods. Ab initio calculations were performed with the plane-wave based Vienna Ab initio Simulation Package (VASP)67,68. The valence electrons of Cr, Co and Ni are, respectively, 6, 9, 10 in the calculation and the projector augmented wave (PAW) method69 was used for the pseudo-potential treatment of core and semi-core electrons. In our calculation, the crystal structures, with periodic boundary conditions, typically contain 300-500 atoms, depending on the specific crystal structure. Atoms are randomly assigned to the lattice points of


1. Gludovatz, B. et al. A fracture-resistant high-entropy alloy for cryogenic applications. Science 345, 1153-1158 (2014).

2. Zhao, Y. H. et al. Simultaneously increasing the ductility and strength of ultra-fine-grained pure copper. Adv. Mater. 18, 2949-2953 (2006).

3. Lu, K., Lu, L. & Suresh, S. Strengthening materials by engineering coherent internal boundaries at the nanoscale. Science 324, 349-352 (2009).

4. Lu, L., Chen, X., Huang, X. & Lu, K. Revealing the maximum strength in nanotwinned copper. Science 323, 607-610 (2009).

5. Zhu, T., Li, J., Samanta, A., Kim, H. G. & Suresh, S. Interfacial plasticity governs strain rate sensitivity and ductility in nanostructured metals. Proc. Natl Acad. Sci. USA 104, 3031-3036 (2007).

6. Wei, Y. et al. Evading the strength-ductility trade-off dilemma in steel through gradient hierarchical nanotwins. Nat. Commun. 5, 3580 (2014).

7. Kou, H., Lu, J. & Li, Y. High-strength and high-ductility nanostructured and amorphous metallic materials. Adv. Mater. 26, 5518-5524 (2014).

8. Gludovatz, B. et al. Exceptional damage-tolerance of a medium-entropy alloy CrCoNi at cryogenic temperatures. Nat. Commun. 7, 10602 (2016).

9. Otto, F. et al. The influences of temperature and microstructure on the tensile properties of a CoCrFeMnNi high-entropy alloy. Acta Mater. 61, 5743-5755 (2013).

10. Frommeyer, G., Briix, U. & Neumann, P. Supra-ductile and high-strength manganese-TRIP/TWIP steels for high energy absorption purposes. ISIJInt. 43, 438-446 (2003).

11. Karaman, I., Sehitoglu, H., Gall, K., Chumlyakov, Y. & Maier, H. Deformation of single crystal Hadfield steel by twinning and slip. Acta Mater. 48, 1345-1359 (2000).

12. Karaman, I. et al. Modeling the deformation behavior of Hadfield steel single and polycrystals due to twinning and slip. Acta Mater. 48, 2031-2047 (2000).

13. Deng, Y. et al. Design of a twinning-induced plasticity high entropy alloy. Acta Mater. 94, 124-133 (2015).

14. Li, Z., Pradeep, K. G., Deng, Y., Raabe, D. & Tasan, C. C. Metastable high-entropy dual-phase alloys overcome the strength-ductility trade-off. Nature 534, 227-230 (2016).

15. Ritchie, R. O. The conflicts between strength and toughness. Nat. Mater. 10, 817-822 (2011).

16. Wu, Z., Bei, H., Pharr, G. M. & George, E. P. Temperature dependence of the mechanical properties of equiatomic solid solution alloys with face-centered cubic crystal structures. Acta Mater. 81, 428-441 (2014).

17. Romaner, L., Ambrosch-Draxl, C. & Pippan, R. Effect of rhenium on the dislocation core structure in tungsten. Phys. Rev. Lett. 104, 195503 (2010).

18. Leyson, G. P., Curtin, W. A., Hector, Jr. L. G. & Woodward, C. F. Quantitative prediction of solute strengthening in aluminium alloys. Nat. Mater. 9, 750-755 (2010).

19. Liddicoat, P. V. et al. Nanostructural hierarchy increases the strength of aluminium alloys. Nat. Commun. 1, 63 (2010).

20. Xie, K. Y. et al. The effect of pre-existing defects on the strength and deformation behavior of a-Fe nanopillars. Acta Mater. 61, 439-452 (2013).

21. Yu, Q., Li, S., Minor, A. M., Sun, J. & Ma, E. High-strength titanium alloy nanopillars with stacking faults and enhanced plastic flow. Appl. Phys. Lett. 100, 063109 (2012).

22. Zhang, Z. et al. Nanoscale origins of the damage tolerance of the high-entropy alloy CrMnFeCoNi. Nat. Commun. 6, 10143 (2015).

23. Kumar, K., Van Swygenhoven, H. & Suresh, S. Mechanical behavior of nanocrystalline metals and alloys. Acta Mater. 51, 5743-5774 (2003).

24. Chen, J., Lu, L. & Lu, K. Hardness and strain rate sensitivity of nanocrystalline Cu. Scripta Mater. 54, 1913-1918 (2006).

25. Li, H. & Ebrahimi, F. Tensile behavior of a nanocrystalline Ni-Fe alloy. Acta Mater. 54, 2877-2886 (2006).

26. Gladman, T. Precipitation hardening in metals. Mater. Sci. Technol. 15, 30-36 (1999).

27. Zimmerman, A., Palumbo, G., Aust, K. & Erb, U. Mechanical properties

of nickel silicon carbide nanocomposites. Mate. Sci. Eng. A 328, 137-146 (2002).

28. Dao, M., Lu, L., Shen, Y. F. & Suresh, S. Strength, strain-rate sensitivity and ductility of copper with nanoscale twins. Acta Mater. 54, 5421-5432 (2006).

29. Yan, F. K., Liu, G. Z., Tao, N. R. & Lu, K. Strength and ductility of 316L austenitic stainless steel strengthened by nano-scale twin bundles. Acta Mater.

60, 1059-1071 (2012).

30. Yeh, J. W. et al. Nanostructured High-entropy alloys with multiple principal elements: novel alloy design concepts and outcomes. Adv. Eng. Mater 6, 299-303 (2004).

31. Cantor, B., Chang, I. T. H., Knight, P. & Vincent, A. J. B. Microstructural development in equiatomic multicomponent alloys. Mater. Sci. Eng. A 375-377, 213-218 (2004).

32. Otto, F., Yang, Y., Bei, H. & George, E. P. Relative effects of enthalpy and entropy on the phase stability of equiatomic high-entropy alloys. Acta Mater.

61, 2628-2638 (2013).

33. Zaddach, A. J., Niu, C., Koch, C. C. & Irving, D. L. Mechanical properties and stacking fault energies of NiFeCrCoMn high-entropy alloy. JOM 65, 1780-1789 (2013).

34. Smith, T. et al. Atomic-scale characterization and modeling of 60° dislocations in a high-entropy alloy. Acta Mater. 110, 352-363 (2016).

35. Gali, A. & George, E. P. Tensile properties of high- and medium-entropy alloys. Intermetallics 39, 74-78 (2013).

36. Wu, Z., Bei, H., Otto, F., Pharr, G. M. & George, E. P. Recovery, recrysta-llization, grain growth and phase stability of a family of FCC-structured multi-component equiatomic solid solution alloys. Intermetallics 46, 131-140 (2014).

37. Ferreira, P. & Mullner, P. A thermodynamic model for the stacking-fault energy. Acta Mater. 46, 4479-4484 (1998).

38. Kibey, S., Liu, J., Curtis, M., Johnson, D. & Sehitoglu, H. Effect of nitrogen on generalized stacking fault energy and stacking fault widths in high nitrogen steels. Acta Mater. 54, 2991-3001 (2006).

39. Suzuki, H. Segregation of solute atoms to stacking faults. J. Phys. Soc. Jpn 17, 322-325 (1962).

40. Schramm, R. & Reed, R. Stacking fault energies of seven commercial austenitic stainless steels. Metall. Trans. A 6, 1345-1351 (1975).

41. Koizumi, Y. et al. Strain-induced martensitic transformation near twin boundaries in a biomedical Co-Cr-Mo alloy with negative stacking fault energy. Acta Mater. 61, 1648-1661 (2013).

42. Tadmor, E. & Bernstein, N. A first-principles measure for the twinnability of FCC metals. J. Mech. Phys. Solids 52, 2507-2519 (2004).

43. Niu, C. et al. Spin-driven ordering of Cr in the equiatomic high entropy alloy NiFeCrCo. Appl. Phys. Lett. 106, 161906 (2015).

44. Ma, D., Grabowski, B., Kormann, F., Neugebauer, J. & Raabe, D. Ab initio thermodynamics of the CoCrFeMnNi high entropy alloy: importance of entropy contributions beyond the configurational one. Acta Mater. 100, 90-97 (2015).

45. Vitos, L., Skriver, H. L., Johansson, B. & Kollar, J. Application of the exact muffin-tin orbitals theory: the spherical cell approximation. Comput. Mater. Sci. 18, 24-38 (2000).

46. Zunger, A., Wei, S.-H., Ferreira, L. & Bernard, J. E. Special quasirandom structures. Phys. Rev. Lett. 65, 353 (1990).

47. Vitos, L., Korzhavyi, P. A., Nilsson, J. & Johansson, B. Stacking fault energy and magnetism in austenitic stainless steels. Phys. Scr. 77, 065703 (2008).

48. Joansson, H., Mills, G. & Jacobsen, K. W. in Classical and Quantum Dynamics in Condensed Phase Simulations (eds Berne, B. J. et al.) 385 (World Scientific Publ. Co, 1998).

49. Bernstein, N. & Tadmor, E. Tight-binding calculations of stacking energies and twinnability in fcc metals. Phys. Rev. B 69, 094116 (2004).

50. Li, W. et al. First-principles prediction of the deformation modes in austenitic Fe-Cr-Ni alloys. Appl. Phys. Lett. 108, 081903 (2016).

51. Liu, L., Wang, J., Gong, S. K. & Mao, S. X. High resolution transmission electron microscope observation of zero-strain deformation twinning mechanisms in Ag. Phys. Rev. Lett. 106, 175504 (2011).

52. Ernst, F. et al. Theoretical prediction and direct observation of the 9R structure in Ag. Phys. Rev. Lett. 69, 620 (1992).

53. Wolf, U., Ernst, F., Muschik, T., Finnis, M. & Fischmeister, H. The influence of grain boundary inclination on the structure and energy of s = 3 grain boundaries in copper. Philos. Mag. A 66, 991-1016 (1992).

54. Schmidt, C., Finnis, M., Ernst, F. & Vitek, V. Theoretical and experimental investigations of structures and energies of S = 3,[112] tilt grain boundaries in copper. Philos. Mag. A 77, 1161-1184 (1998).

55. Rohatgi, A., Vecchio, K. S. & Gray, III G. T. The influence ofstacking fault energy on the mechanical behavior of Cu and Cu-Al alloys: Deformation twinning, work hardening, and dynamic recovery. Metall. Mater. Trans. A 32, 135-145 (2001).

56. Afanasyev, K. A. & Sansoz, F. Strengthening in gold nanopillars with nanoscale twins. Nano Lett. 7, 2056-2062 (2007).

57. Zheng, Y., Lu, J., Zhang, H. & Chen, Z. Strengthening and toughening by interface-mediated slip transfer reaction in nanotwinned copper. Scripta Mater. 60, 508-511 (2009).

58. Zhu, L., Qu, S., Guo, X. & Lu, J. Analysis of the twin spacing and grain size effects on mechanical properties in hierarchically nanotwinned face-centered cubic metals based on a mechanism-based plasticity model. J. Mech. Phys. Solids 76, 162-179 (2015).

59. Coujou, A. Deformation in situ d'un alliage a basse energie de faute d'empilement. Acta Metall. 31, 1505-1515 (1983).

60. Priester, L. Grain Boundaries: From Theory to Engineering Vol. 172 (Springer Science & Business Media, 2012).

61. Priester, L. & Khalfallah, O. Image force on a lattice dislocation due to a grain boundary in anisotropic fcc materials. Philos. Mag. A 69, 471-484 (1994).

62. Itakura, M., Kaburaki, H., Yamaguchi, M. & Tsuru, T. Novel Cross-slip mechanism of pyramidal screw dislocations in magnesium. Phys. Rev. Lett. 116, 225501 (2016).

63. Vanderschaeve, G. Cross-slip of partial dislocations via the stair rod mode. The zigzag propagation of deformation microtwins in ordering alloys. Phys. Status Sol. (A) 100, 59-68 (1987).

64. Shen, Y. F., Lu, L., Lu, Q. H., Jin, Z. H. & Lu, K. Tensile properties of copper with nano-scale twins. Scripta Mater. 52, 989-994 (2005).

65. Yu, Q. et al. Origin of dramatic oxygen solute strengthening effect in titanium. Science 347, 635-639 (2015).

66. Shan, Z., Mishra, R. K., Asif, S. S., Warren, O. L. & Minor, A. M. Mechanical annealing and source-limited deformation in submicrometre-diameter Ni crystals. Nat. Mater. 7, 115-119 (2008).

67. Kresse, G. & Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. Rev. B 47, 558 (1993).

68. Kresse, G. & Hafner, J. Ab initio molecular-dynamics simulation of the liquid-metal-amorphous-semiconductor transition in germanium. Phys. Rev. B 49, 14251 (1994).

69. Kresse, G. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. Rev. B 59, 1758 (1999).

70. Perdew, J. P. et al. Atoms, molecules, solids, and surfaces: applications of the generalized gradient approximation for exchange and correlation. Phys. Rev. B 46, 6671 (1992).


This work was supported by the Chinese 1000-Youth-Talent Plan (for Q.Y.) and the State Key Program for Basic Research in China under Grant No. 2015CB65930 (for Z.Z., Z.W., Z.Z. & S.X.M.), and by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, Materials Sciences and Engineering Division through the Materials Science and Technology Division at the Oak Ridge National Laboratory (for E.P.G.) and through the Mechanical Behaviour of Materials Program (KC13) under contract no. DE-AC02-05CH11231 at the Lawrence Berkeley National Laboratory

(LBNL) (for B.G. and R.O.R.). Work at GMU was supported by the NSF under Grant No. DMR-1611064. We acknowledge the use of the aberration-corrected TEAM0.5 transmission electron microscope at the National Center for Electron Microscopy at LBNL, which is supported by the Office of Science of the U.S. Department of Energy also under contract no. DE-AC02-05CH11231 and the Center for Computational Materials Science at the Institute for Materials Research, Tohoku University, for use of the SR16000 supercomputing facilities.

Author contributions

Q.Y. and R.O.R. conceived the study; H.S. conducted the numerical simulations; Z.Z. and Z.W. performed the experiments (directly supervised by Z.Z., Q.Y. and S.X.M.); H.S., B.G., Z.Z., E.P.G., Q.Y., S.X.M. and R.O.R. analysed the results, and H.S., Q.Y., S.X.M. and R.O.R. wrote the manuscript with comments and input from B.G. and E.P.G.

Additional information

Supplementary Information accompanies this paper at naturecommunications

Competing financial interests: The authors declare no competing financial interests.

Reprints and permission information is available online at reprintsandpermissions/

How to cite this article: Zhang, Z. et al. Dislocation mechanisms and 3D twin architectures generate exceptional strength-ductility-toughness combination in CrCoNi medium-entropy alloy. Nat. Commun. 8, 14390 doi: 10.1038/ncomms14390 (2017).

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit

© The Author(s) 2017