Modeling of Stacking Fault Energy in Hexagonal-Close-Packed Metals

The deformation of metals is known to be largely affected by their stacking fault energies (SFEs). In the review, we examine the theoretical background of three normally used models, supercell model, Ising model, and bond orientation model, for the calculation of SFE of hexagonal-close-packed (hcp) metals and their alloys. To predict the nature of slip in nanocrystalline metals, we further review the generalized stacking fault (GSF) energy curves in hcp metals and alloys. We conclude by discussing the outstanding challenges in the modeling of SFE and GSF energy for studying the mechanical properties of metals.


Introduction
Due to the high fracture toughness, uniform property in all directions, conductivity, and recyclability, metals are dominant in engineering applications.Strength and ductility are two important mechanical properties for metal structure materials; however, it is difficult for a material to have both high strength and high ductility [1].For face-centered cubic (fcc) metals, such as Cu alloys, deformation twinning is a useful mechanism that can simultaneously enhance the strength and ductility [2,3].The factors affecting the deformation twinning include strain rate [4], deformation temperature [5], grain size [6], and stacking fault energy (SFE) [7,8].In contrast, the hexagonal-close-packed (hcp) metals usually have low plasticity due to limited independent slip systems.In hcp systems, Kwaśniak et al. [9] revealed that the interplanar distance of prismatic and pyramidal slip planes is different.They listed the possible slip system of hcp metals as ⟨⟩ basal, ⟨⟩ prismatic I, ⟨⟩ prismatic II, ⟨⟩ pyramidal I, ⟨⟩ pyramidal II, ⟨ + ⟩ pyramidal I, and ⟨ + ⟩ pyramidal II, where I and II denote large or small interplanar distance.This dual nature of pyramidal and prismatic planes has not been reported in experimental or theoretical studies.
SFE is commonly accepted as a useful quantity for understanding the plastic deformation mechanisms of metals.Deformation twinning is favored at low SFE and dislocation slip at high SFE.In high-manganese austenitic steels, for instance, the mechanical twinning is for low values, while dislocation slip is for high values of the SFE [10].The introduction of additional alloy elements is also an effective approach to improve the mechanical properties of metals.By adding alloy atoms, the SFE of Cu-Zn alloy has been readily adjusted, producing high strength and high ductility.Suzuki et al. [11] claimed that, with the addition of Zn, the creep strength of Mg-Y alloys improved at 550-650 K and decreased the SFE values.In addition, Kawamura et al. [12] synthesized Mg 97 Y 2 Zn 1 alloys by a rapidly solidified powder metallurgy technique and achieved a high tensile strength of 600 MPa and elongation of 5%, along with a decrease of the SFE value.
However, the effect of solute atoms on SFE remains unclear, which is mainly because the magnitudes of SFE are small, being difficult for experimental measurement.Accordingly, density-functional theory (DFT) becomes an effective and popular method to calculate the SFEs for metal surfaces.
In this review, we will focus on the modeling of the SFE in hcp metals.In Section 2, we will briefly introduce the definition of stacking fault and generalized stacking fault (GSF) energy.The models used for first principle calculations will be elucidated in Section 3. In Section 4, we will compare 2 Advances in Materials Science and Engineering the GSF energies determined by different DFT functional.The challenge and outlook for the modeling methods will be shown in Section 5.

Stacking Fault Energy and Generalized Stacking Fault Energy
On the microscale, atoms were described as hard spheres which change size from one element to another.Perfect crystals can be described by the hard spheres stacking at three-dimensional space in a regular sequence.The atoms in one layer are identical; the sequence of atom layers in the crystal is described by labeling one layer and other layers with atoms in identical positions as an A layer; layers of atoms in other positions in the stack are referred to as B or C layers.
Stacking faults, which appear at the close-packed planes in the close-packed structures (fcc and hcp), describe the change of local region in crystals by removal or introducing external layers.Stacking faults destroy the perfection of the ideal structure, and the accompanied energy per unit area of faults is known as SFE  [13].
With the model described above, the energy of an atom in the system is determined by the local environment (layers close to it).In perfect crystals, hcp structures have a local hcp-like environment and fcc structures have a local fcclike environment.The main contribution to  arises from the changes of the local environment of atoms.For example, in the ideal hcp structure, the stacking sequence is . . .ABABAB . .., but, in the faulted sequences, some planes have a local fcc-like environment.There are one fcc-like atom in  1 , two in  2 and  2 , and three in ; then the SFE satisfies [13].Experimental estimates of  often show large uncertainties and are different from the theoretical estimate.For instance, the SFEs measured for magnesium are reported to range from 50 mJ/m 2 to 280 mJ/m 2 [14].Fortunately, the DFT calculations could give reliable estimates of the SFEs at reasonable computer cost.For hcp metals, the basal close-packed plane is (0001) plane and the close-packed direction is ⟨1120⟩.Along this direction, the dislocation glide is frequently observed.Mg has a / ratio of 1.632, which is close to ideal close-packed atomic spheres.Its dislocation can also easily slip along the (1/3)⟨1120⟩ direction on the first-order prism planes {1010} and pyramidal plane {1011} [13].
Since the stacking fault described above only locates at stable sites, it is thus named "stable SFE".In contrast, the SFE related to the slip process is referred to as the "unstable SFE".When studying the deformation properties or dislocation core, it is important to export the evolution of SFEs from one stable position to another.The generalized stacking fault (GSF) energy surface (-surface) can provide a description of stacking faults evolution.The GSF energy was first introduced by Vítek, for the description of the dislocation core structure and mobility [15].The GSF energy is described by the difference of two supercells, which can be used to model stacking faults with the burger vectors  = 0 (stable position) and  ̸ = 0 (unstable position).

Theoretical Models for Stacking Fault Energy
For SFE calculations, there are three typical models: 3.1.Supercell Model.In 1980s, Chou et al. [16] placed stacking faults in a periodic repeated supercell, where two large supercells are required for SFE calculations.One is the structure of perfect crystal; the other is the structure containing stacking fault; calculate the energy of the two structures by first principle calculations directly.Then, the definition of SFE  can be expressed as where  SF denotes the energy of stacking fault structure,  0 the energy of perfect structure, and  the surface area.
Based on this model, four basal plane SFEs were calculated by Chetty and Weinert using LDA [17] and by Smith using GGA [18] (cf.Table 1).The  1 and  2 SFEs are found to be the most stable and to have a lower energy than the ideal fcc Mg.However, the results of SFE are closely related to the size of the supercell.For example, Wang et al. [19] revealed that when the supercell size increases from 8 to 32, the SFE reduced by 50%.Recently, several long-period stacking ordered (LPSO) structures have been found in Mg-Y alloys; these alloys show superplastic behavior at high temperatures [20,21].First principle calculations were used to investigate the LPSO structures and showed that the SFEs of these structures only contain  1 and  2 type stacking fault [22].The SFEs of other hcp metals (Be, Sc, Y, Ti, Zr, Tc, Re, Zn, Cd, etc. [23,24]) and fcc metals (Al, Ni, Cu, Ag, Au, etc. [25]) have also been studied by using DFT.

Ising Model.
Axial next-nearest-neighbor Ising model was firstly applied to describe the polytypic structure of Chou et al. [16].The analogue spin   has value 1 or −1 for layer ; the sign is determined by the connecting subsequent layer from the close-packed layers.For example, if  layer is A and  + 1 layer is B, then   = 1; if  layer is A and  + 1 layer is C, then   = −1.Based on this model, the energy of arbitrary structure can be described by a series expressed as where   is the interaction energy between two layers ( = 1 is the first-nearest-neighbor,  = 2 is the second-nearestneighbor, etc.). 0 is the energy where all the interlayer interactions are neglected.Denteneer and van Haeringen [26] provided the energies of three different structures.For the perfect hcp structure (. . .ABABAB . ..), where ( 3 ) represents the interaction energies between layers with distance over the second-nearest-neighbor.For the perfect fcc structure (. . .ABCABC . ..), For the double hcp (dhcp) structure (. . .ABAC . ..), Under this definition and neglecting the series over the second-nearest-neighbor, the energies of the systems of hcp, fcc, and dhcp structures can be calculated by the supercell model.Then the parameters  0 ,  1 , and  2 can be determined by solving (7) to (9): Thus, the energy of arbitrary stacking sequence can be expressed by   .
According to the Ising model, four types of basal plane SFEs of 14 hcp metals were calculated; the results divided hcp metals into three types for the different interlayer interaction features (the results of   ) [23].Sandlöbes et al. [20] calculated the SFE of Mg and Mg-Y alloys and found that SFE was decreased by the addition of yttrium.Lu et al. [27] investigated the composition dependence of the SFEs by adding different alloy atoms (Mn, Co, and Nb) to austenitic stainless steels.They found that Co decreases the SFE, whilst Nb increases the SFE value.Moreover, Mn decreases the SFE in alloys with less than 16 at.% Ni; beyond that the SFE rises slightly.Notably, the SFE of Ni was calculated by both the supercell model (127 mJ/m 2 ) and Ising model (129 mJ/m 2 ), and both methods agree nicely with experiments (125 mJ/m 2 ) [25].

Bond Orientation Model. Chetty and Weinert further
analyzed the local environment model and found that the energy also depends on the local environment of the nearestneighbors [17].Specifically, when considering the local environment of the nearest-neighbors, there are six different energies of a site shown in Table 2.With energies of hcp, fcc,  1 ,  2 , , and  2 structures calculated by the supercell model, the results of  2 to  6 (listed in Table 2) relative to hcp can be determined by The parameters  1 to  6 can be calculated by solving (11) and the results are listed in Table 2.
This model is effectively a second-neighbor model, and the energies of arbitrary fault structures can be described by these parameters.To test the ability of this model in

Methods for GSFE Calculation.
There are two commonly used methods for GSFE calculation: One is computed by direct crystal slip path, in which atoms are only relaxed along the directions perpendicular to the slip planes [28].The other one is named as climbing image nudged elastic band (CINEB) method, which describes the -surface with precise values of saddle points without constrained way of crystal slips [9].

Calculation of Generalized Stacking Fault Energy
4.1.Pure Metal.For Mg, the -surface of basal plane along the (1/3)⟨1010⟩ direction was listed by Smith [18].The respective stable SFE from the curve is 36 mJ/m 2 .Furthermore, the surface of four basal plane stacking faults ( 1 ,  2 , , and  2 ) and others of nonbasal plane stacking faults in the prismatic and pyramidal planes were gained using supercell model.The results show that basal plane stacking fault is the most likely slip direction; the slipping along ⟨1010⟩ could also form a stable stacking fault structure, but its SFE is higher than that of basal plane slip system.Thus, it is difficult to form nonbasal plane stacking faults [28].The GSF energies of 2layer (AB), 4-layer (dhcp), and 6-layer (ABACAB) structures on basal plane slip were calculated in [22]; the calculated results of unstable SFE for 4-layer and 6-layer structures are higher than that of 2-layer structure.The GSF energies of other hcp metals (Be, Ti, Zn, and Zr) were calculated and the -surface of {0001}⟨1120⟩ and {0001}⟨1010⟩ slip system was shown in Figure 1.The curves illustrate that for different hcp metals the unstable SFE presents at different positions and the values of unstable SFE along ⟨1010⟩ are always larger than that along ⟨1120⟩ direction.Then, compared by employing DFT within GGA and spin-polarized GGA, results show that the primitive cell volumes influence the GSF energy of hcp metals.For transition metals a spin-polarized calculation should be considered [24].

Metal Alloys.
The GSF energy surface of Mg alloys has been extensively studied by DFT.Muzyk et al. [29] calculated GSF energy curves of basal plane involving 13 types of metallic alloying elements (Figure 2).They found that the SFE of Mg reduced obviously when adding Pb and Sn atoms.Han et al. [30] considered the effects of alloy atoms on the GSF energy of the basal plane.They concluded that adding Al atoms tends to decrease the SFEs along the faulting pathways, but adding Li atoms tends to increase the SFEs along the faulting pathways.Kwasniak et al. [33] investigated the effect of C, H, N, and O on the GSFE of -Ti.They found that solute H can reduce the SFE considerably, solute C can improve partial dislocation emission and twin content, and solutes N and O can reduce the SFE in considered slip mode [33].Ghazisaeidi and Trinkle [34] considered the energetics of O interstitial interactions with a twin boundary and a prismatic stacking fault of Ti, proving that sites located away from the twin boundary are repulsive, while the sites at twin boundary are more attractive to O [34].

Challenge and Outlook
In this review, we summarized the theoretical models used for SFE calculations and their applications in metal systems.The available data clearly show that different exchangecorrelation functionals in the framework of DFT can give different results, with large error bars.We thus conclude that it remains difficult to accurately determine the SFE by using the currently used DFT functional.In fact, to accurately describe such systems, some important factors, such as van der Waals forces, many body effects, and dielectric screening effects, should also be carefully considered in DFT calculations.Therefore, the development and application of new DFT that properly include all the above effects will be highly demanded.

( 1 )
supercell model, (2) Ising model, and (3) bond orientation model.The most commonly used supercell model based on the geometry structure will be shown in Section 3.1.The other two models based on different parameters will be shown in Sections 3.2 and 3.3.

Table 1 :
Calculated stable SFE (mJ/m 2 ) and unstable SFE (mJ/m 2 ) for Mg and Mg alloys by DFT with different exchange-correlation functionals.

Table 2 :
The local environment of atoms  and  ± 1; hcp and fcc are atoms of local environment.The parameters  1 to  6 denote the energy of atom .+ 2 5 + 4 6 = 142 mJ/m 2 ) and supercell model (148 mJ/m 2 ).The energy of a randomly chosen structure, ABACBACBCBAC, is  = 2 2 + 2 3 +  4 + 4 5 + 3 8 = 186 mJ/m 2 , while the result calculated by supercell model is 197 mJ/m 2 .The results show that the bond orientation model and the supercell model predict very close SFE for Mg.
[32] etal.[31]considered the GSF energy of a new Mg-Al-Sn alloy of different slip systems.They found that the unstable SFE is reduced due to the doping of Al and Sn atoms.Zhang et al.[32]investigated the surface of  1 and  2 stacking fault when solute atoms Y are added.Compared with the -surface of pure Mg, the  1 and  2 -surfaces of Mg-Y alloy are very similar to that of Mg, as shown in Figure3.When adding Y and Zn simultaneously, the stable and unstable SFEs drastically decrease.
[22]experiment observation of Mg-Y alloys shows the same results.For the LPSO structure, when adding 2% Zn, the energy of 6-layer structure is compared with that of 2-layer one.The results show that addition of Zn further stabilizes the 6-layer structure[22].Another typical hcp metal is Ti.
from Chinese Ministry of Education, National Natural Science Foundation of China (nos.21403113, 51225102, and 2012CB932203), the Fundamental Research Funds for the Central Universities (no.30915011330), a project funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions, and the 8th Liuda Rencai Gaofeng (no.B932203) from Jiangsu Province, China.