First-principles study of the ternary effects on the plasticity of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\upgamma $$\end{document}γ-TiAl crystals

We studied the effects of important ternary elements, such as Cr, Nb, and V, on the plasticity of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\upgamma $$\end{document}γ-TiAl crystals by calculating the point defect formation energy and the change in the generalized stacking fault energy (GSFE) surface from first-principles calculations. For all three elements, the point defect formation energies of the substitutional defects are lower in the Ti site than in the Al site, which implies that substitution on the Ti site is energetically more stable. We computed the GSFE surfaces with and without a substitutional solute and obtained the ideal critical resolved shear stress (ICRSS) of each partial slip. The change in the GSFE surface indicates that the substitution of Ti with Cr, Nb, or V results in an increase in the yield strength because the ICRSS of the superlattice intrinsic stacking fault (SISF) partial slip increases. Interestingly, we find that Cr substitution on an Al site could occur owing to the small difference between the substitutional defect formation energies of the Ti and Al sites. In that case, the reduction of ICRSSs of the SISF partial slip and twinning would lead to improved twinnability. We discuss the implications of the computational predictions by comparing them with experimental results in the literature.

First-principles study of the ternary effects on the plasticity of γ-TiAl crystals Taegu Lee 1 , Seong-Woong Kim 2 , Ji Young Kim 2,3 , Won-Seok Ko 4* & Seunghwa Ryu 1* We studied the effects of important ternary elements, such as Cr, Nb, and V, on the plasticity of γ-TiAl crystals by calculating the point defect formation energy and the change in the generalized stacking fault energy (GSFE) surface from first-principles calculations. For all three elements, the point defect formation energies of the substitutional defects are lower in the Ti site than in the Al site, which implies that substitution on the Ti site is energetically more stable. We computed the GSFE surfaces with and without a substitutional solute and obtained the ideal critical resolved shear stress (ICRSS) of each partial slip. The change in the GSFE surface indicates that the substitution of Ti with Cr, Nb, or V results in an increase in the yield strength because the ICRSS of the superlattice intrinsic stacking fault (SISF) partial slip increases. Interestingly, we find that Cr substitution on an Al site could occur owing to the small difference between the substitutional defect formation energies of the Ti and Al sites. In that case, the reduction of ICRSSs of the SISF partial slip and twinning would lead to improved twinnability. We discuss the implications of the computational predictions by comparing them with experimental results in the literature.
TiAl alloy crystals are intermetallic compounds in which the alpha phase (Ti 3 Al) and the gamma phase (TiAl) form a laminated structure with special orientation relationships. The TiAl alloy is used primarily in the aircraft industry because of its excellent physical properties, such as low density, high strength ( σ YS = 300 − 800MPa ), and good heat resistance for wide range of temperature up to 800 °C 1,2 . However, wider application of the TiAl alloy has been limited because of a critical drawback: low ductility at room temperature that often leads to catastrophic failure 1,2 . To overcome these limitations, many experimental studies have been performed to reveal the deformation mechanism of the TiAl alloy [3][4][5][6] . Moreover, previous researchers tested the feasibility of tuning the mechanical properties by adding a third element, the so called the ternary effect. Because the alpha phase is relatively more brittle owing to its limited slip systems [7][8][9] , most studies related to ternary effects focused on the relatively ductile gamma phase or duplex structures. Previous studies found that the addition of Nb or V (while keeping the atomic fraction of Al constant) increases the yield strength of the gamma TiAl crystal [10][11][12][13][14][15] , and more deformation twins form with the addition of Mn 16 . Other studies reported that the addition of Cr or V atoms only improves the ductility of the duplex TiAl structure, while not affecting the ductility of single phase gamma 17,18 . Despite the extensive investigations on the ternary effect, a fundamental understanding is still lacking on how a ternary atom affects the mechanical properties of the TiAl gamma phase.
To complement the experimental observations, a couple of theoretical studies reported the first-principle density functional theory (DFT) calculation results of the stacking fault energy (SFE) of the gamma phase or the interface energy of the lamellar structures [19][20][21] . The SFE of each can be obtained from γ = E sf −E 0 A , where E sf and E 0 are the energies of crystals with and without the stacking fault, respectively, and A is the area of the slip plane in which a stacking fault lies. Three different stacking faults may form in the gamma TiAl crystal: superlattice intrinsic stacking faults (SISFs), antiphase boundaries (APBs), and complex stacking faults (CSFs). The generalized stacking fault energy (GSFE) surface refers to the entire two-dimensional map made by arbitrary fault vectors that are not necessarily at a local minimum, which offers a more complete picture on the slip system. Moreover, previous research revealed the deformation mechanism of the single nanowire 22 and the lamellar structure 23,24 of the binary TiAl alloy. In addition, other studies found that the occupation sites of ternary interstitial sites and the ternary effects on the changes of the lattice constant 25,26 . However, to deepen the understanding of the ternary www.nature.com/scientificreports/ effect on the mechanical properties, the changes in the stacking faults induced by the ternary atom and the resulting effect on the plastic deformation mechanism need to be predicted. Although the previous calculations show that the changes in the GSFE of the gamma-TiAl-based alloy are due to the ternary substitution 27,28 , the effect of the GSFE change on the plasticity has not been systematically analyzed yet. Furthermore, the relationship between the changes in the GSFE and the mechanical properties by the ternary effects remain to be revealed. In this study, we investigated the changes in the plastic behaviors of the gamma TiAl crystal upon the addition of Cr, Nb, or V by combining computationally obtained GSFEs and theoretical analyses on the plastic deformation mechanism. First, we calculated the point defect formation energy to determine the preferred site (Ti or Al) for a substitutional ternary atom defect. We then calculated the GSFE surface with and without the substitutional ternary defect. The ideal critical resolved shear stress (ICRSS) of each partial slip was obtained from the maximum slope of a portion of the GSFE surface along the partial Burgers vector. Based on these calculations, we discuss the relationship between the GSFE change and plastic behaviors by considering major plastic deformation mechanisms involving the SISF partial slip, the CSF partial slip, and twinning. We compare the predictions with experimental reports on the ternary effect, such as changes in the yield strength and ductility. This paper is organized as follows: In Section 2-A, we describe the slip system of the gamma phase in detail and explain the relationship between the GSFE and deformation mechanisms. In addition, we describe the DFT calculation method in Section 2-B. In Section 3, we show the calculation results of the point defect formation energy and the GSFE. We present the ICRSS of each partial dislocation with and without ternary atoms and discuss the ternary effect on the mechanical properties, such as yield strength and ductility. We summarize the results and discuss the outlook in the final section.

Theoretical background and methodology
Slip system and deformation mechanism of the gamma TiAl crystals. The gamma phase has an L1 0 intermetallic structure that is similar to that of the face-centered cubic (FCC) crystal, except that the lattice constant along the c-axis (represented as the [001] direction in this study) is approximately two percent longer than the lattice constant along other axes, and two different atomic species are stacked alternately along the c-axis. In this study, we adopt the modified Miller indices, which are often used to express the directions and planes in the L1 0 crystal 29 . Similar to FCC crystals, slip in L1 0 crystals occurs in the {111) close packed plane. However, because of the different GSFE surface (e.g., three different stacking faults, including SISF, APB, and CSF, may form in the L1 0 intermetallic structure), the plastic deformation mechanism differs significantly from that of FCC crystals.
In the FCC structure, slip on the {111} plane would begin with the 1/6 < 112 > Burgers vector (or the three equivalent directions of 1/6 11  Fig. 1. Another slip on the same plane along the trailing partial direction would lead to perfect slip, or alternate consecutive slips along the identical direction on the adjacent slip planes would lead to twinning. In contrast, in the L1 0 structure, the first partial slip would result in either a SISF or CSF, and then subsequent slips would lead to five distinct plastic deformation behaviors. Specifically, the first slip on the {111) plane can occur along two distinguishable directions; slip along the 1/6 < 112] direction leads to the formation of a SISF, while slip along either the 1/6 < 211] or 1/6 < 121] direction would form a CSF, as shown in Fig. 2. When the CSF is initially formed, it can form < 110] ordinary slip by additional slip along 1/6 < 121] (or 1/6 < 211]), which is similar to perfect slip in the FCC lattice. Alternatively, a super slip can form if the additional slip following the CSF occurs to complete the < 011] or 1/2 < 112] super slip, which will be referred to as an inverse super slip in this paper. On the other hand, if a SISF is formed in the first step, the slip process may proceed to form either < 011] or 1/2 < 112] super slip, which will be referred as a forward super slip in this paper. Instead, if the initial SISF is followed by www.nature.com/scientificreports/ consecutive slips on the adjacent slip planes along the 1/6 < 112] direction, deformation twinning would occur. Alternatively, instead of the deformation twin, multiple SISFs can form on slip planes that are not adjacent to each other. While the perfect dislocation of the FCC structure consists of two partial dislocations, the L1 0 structure may contain super slip composed of four partial slips or ordinary slip composed of two partial slips 30 as follows: In summary, the plastic deformation mechanisms can be categorized into the five different types as follows (see Fig. 2): (i) inverse super slip and (ii) ordinary slip in Step 2-1, and (iii) forward super slip, (iv) deformation twin, and (v) SISFs in Step 2-2. In a previous study 22 , we suggested a framework for predicting the deformation mechanism of a single crystal subjected to uniaxial loading, depending on the loading type (compression or tension) and direction. To determine the type of plastic deformation, the following properties must be considered: the ICRSS, the minimum resolved shear stress required to form the partial slip in the absence of thermal www.nature.com/scientificreports/ activation, and the Schmid factor, the geometrical factor relating the slip direction and the loading direction. The ICRSS divided by the Schmid factor yields the ideal critical stress required to initiate plastic deformation in a crystal subjected to external loading in the absence of thermal activation. The ICRSS can be obtained from the maximum slope of the GSFE curve along the corresponding partial slip direction. It was shown that the ICRSS has a positive correlation with the critical resolved shear stress and also the critical stress at finite temperatures 31 . For a single crystal subjected to external loading beyond yield, the plastic deformation begins with either the CSF partial slip or the SISF partial slip, which has the lower critical stress of the two. The deformation following the initial partial slip can also be determined by comparing the critical stresses. For instance, if a CSF is formed first, then the slip plane moves in the direction with the lower critical stress between the inverse super slip and the ordinary slip, as shown in Step 2-1. Conversely, if SISF partial slip initially occurs, the plastic deformation proceeds with the mechanism having the lower critical stress, either the forward super slip or the twin in Step 2-2. Based on calculations that show that the CRSS of the twin is always lower than the CRSS of the SISF, we excluded case (v). In our previous study, where full details on the deformation mechanism prediction can be found, we showed that the deformation mechanism predicted based on this framework matches well with molecular dynamics simulations of single crystal nanowires under uniaxial tensile and compressive loading 22 .
Because the plastic deformation behavior is closely related to the mechanical properties of the material, the ICRSS can serve a qualitative indicator for the mechanical properties. The ICRSS of the SISF or the CSF, which occurs first at the onset of plastic deformation, is not only the criterion for dislocation nucleation but also highly related to the ideal shear strength in bulk materials 31 . In particular, the ICRSS of the SISF is lower than that of the CSF (regardless of the ternary atom substitution, as will be shown later), and the change in the yield strength of a polycrystal due to a ternary atom can be inferred from the change in the ICRSS of the SISF. Second, the propensity for deformation twinning is reported to increase the ductility because the twin boundary serves as a barrier for dislocation motion [32][33][34] . Twin boundaries limit the mean free path of mobile dislocations, and thus plastic deformation becomes more distributed, which enhances the overall ductility. Consequently, it is reasonable to suspect that if the ICRSSs of the SISF and twinning are reduced because of a ternary atom, more twin boundaries would form, and thus the ductility of the material would increase.
Calculation details of the point defect formation energy and the GSFE. The DFT calculation of the GSFE surface was performed by employing the Vienna Ab initio Simulation Package (VASP) software 35 with the projector augmented wave (PAW) potentials 36 . The generalized gradient approximation (GGA) as parametrized by Perdew-Burke-Ernzerhof (PBE) 37 was used to treat electron change and correlation. In the PAW potential for Ti, 3p electrons were treated as part of the valence. An energy cutoff of 450 eV and a 5 × 5 × 5 Monkhorst-Pack k-point mesh 38 were used for the formation energy calculations for all supercells considered in the study, and the same energy cutoff and a 9 × 7 × 3 Monkhorst-Pack k-point mesh were used for the SFE calculations.
Three different elements (Cr, Nb, and V), which have been widely used in experiments, were chosen as the ternary elements in the present study. In order to determine whether a Ti or an Al site is more stable for the ternary element substitution, we defined the point defect formation energy at 0 K which is the difference between the formation energies of a ternary alloy and the stoichiometric gamma TiAl as follows: Ti N Al N indicates stoichiometric TiAl consisting of N atoms for each element, and Al N−1 or Ti N−1 indicates that one atom of the corresponding sublattice site is replaced by an X (Cr, Nb, or V) atom. The energy of a single atom was calculated for a stable structure under the ambient condition (i.e., FCC for Al, HCP for Ti, and BCC for Cr, Nb, and V).
To check possible artifacts from the finite size effect, L1 0 supercells (x: [ 100 ], y: [ 010 ], z: [ 001 ]) consisting of 32 atoms (ternary element ~ 3 at.%) and 108 atoms (ternary element ~ 1 at.%) were used. Both the unrelaxed and relaxed point defect formation energies of all three species of atoms were lower when the ternary atom occupied a Ti site (The detailed explanation is in the Results and discussion section (Section 3)).
Next, the GSFE surface at 0 K was computed with and without the three different ternary elements. Figure 3a shows the L1 0 supercell (x: [ 112 ], y: [ 110 ], z: [ 111 ]) consisting of 24 atoms that was used to calculate the SFE of the ternary TiAl-based alloys. Owing to the computational cost, smaller cells were used for the SFE calculations than the one used for the formation energy calculations. In addition, the calculation of the formation energy was repeated for a smaller supercell. Despite the size and orientation of the supercell, the relative order between the formation energies for the Ti and Al sites did not change. When the ternary atom was substituted for a Ti atom or an Al atom in the supercell, the atomic percent of the ternary atom was approximately 4 at.%, as shown in Fig. 3b,c. Because the GSFE was obtained by tilting the periodic boundaries while maintaining the real positions of atoms, as depicted in Fig. 4, the topmost Ti or Al atom was substituted with the ternary atom to locate the substitutional defect adjacent to the designated slip plane. Given the three vectors of the initial supercell periodic boundaries (namely, − → a 1 , − → a 2 , and − → a 3 , where − → a 1 and − → a 2 are parallel to the slip plane), the total energy of the system was computed as − → a 3 was changed following − → a 3 ′ = n − → a 3 + − → b , where the relative displacement, − → b = l − → a 1 + m − → a 2 , was intended to scan the entire slip vector on the slip plane, and n − → a 3 was adjusted to obtain relaxed GSFEs at a given slip vector, − → b . In detail, − → a 1 , − → a 2 , and − → a 3 were 4.953, 5.636, and 14 Å, respectively. Values of l , m , and n were 0, 1 20 , 2 20 , . . . , 1 , 0, 1 12 , 2 12 , . . . , 1 , and 23 28 , 25 28 , . . . , 33 28 , respectively. In the other words, the www.nature.com/scientificreports/ relaxed GSFE surface was obtained by computing the energies for different values (21, 13, and 6) of l, m, andn , respectively (in total, 1638 cases). Afterward, starting from a simulation cell possessing a SISF, the GSFE curve was obtained by considering the slip along the slip plane adjacent to the SISF to calculate the ICRSS of twinning. For the twinning, two scenarios were considered, as visualized in Fig. 5a,b. The first scenario considered the case in which the ternary atom was located next to the initial slip plane but apart from the subsequent slip plane (Fig. 5a), while the second scenario assumed that the ternary atom was located between two consecutive slip planes (Fig. 5b). When predicting the deformation mechanism, the ICRSS of twinning was set to the smaller maximum slope between the two scenarios. Table 1 lists the point defect formation energies obtained from the first-principles DFT calculations. The point defect formation energies obtained from the supercells used for the GSFE are presented in the Supplementary  Information (Fig. S2 and Tables S1-S2). The results reveal that Cr, Nb, and V defects preferentially substitute on the Ti site over the Al site as demonstrated by the lower point defect formation energies of substitutional point defects for the Ti site compared with those for the Al site. The difference between the substitutional energies for Ti and Al sites are the largest for Nb and the smallest for Cr. For Nb and V, the difference between the point defect formation energies for the Ti and Al sites is significantly larger than the thermal energy at room temperature (~ 0.026 eV). Hence, unless a ternary gamma TiAl alloy with Nb or V is quenched at an extremely fast cooling rate, the ternary element is likely to occupy the Ti sites. In contrast, the point defect formation energy difference for Cr substitutional defects is comparable to the thermal energy, and a non-negligible fraction of Cr would  www.nature.com/scientificreports/ occupy some portion of the Al sites. The formation energy calculation is consistent with previous ALCHEMI (atom location by channeling-enhanced microanalysis) experiments on a single gamma phase indicating that V occupies only the Ti sites and Cr occupies both the Ti and Al sites 17,18 . The calculated GSFE surface is shown in Fig. 6a. The blue dashed line indicates the slip path for forward super slip, the red dashed line indicates the slip path for inverse super slip, and the black dashed line indicates the slip path for ordinary slip. By taking the GSFEs along the paths and augmenting the GSFE curve for twinning, a set of GSFE curves can be extracted, as shown in Fig. 6b,c. Table 2 summarizes the SFE values for the SISFs, APBs, and CSFs obtained for each ternary alloy. The GSFE curves for when a Cr atom occupies either a Ti or an Al site are shown in Figs. 6 (d)-(g). The GSFE curves for Nb and V substitution are depicted in Fig. S3 in the Supplementary Information. Here, the curve from the first twinning scenario appears as a solid line, while the curve from the second twin scenario is shown as a dotted line. The ICRSS was obtained for various deformation mechanisms by computing the maximum slopes for three parts (SISF formation, twinning, and forward super slip) of the GSFE curves starting from the SISF formation, as well as maximum slopes for three parts (CSF formation, ordinary slip, and inverse super slip) of the GSFE curves starting from the CSF formation. Table 3 summarizes the ICRSSs for all cases and the relative changes in the ICRSSs of ternary alloys (for both Ti and Al site substitution) with respect to the values of the stoichiometric gamma TiAl alloy. For all ternary elements, two remarkable common trends are present. First, the substitution of a ternary atom in a Ti site increases the ICRSS of the SISF partial slip. Second, the substitution of a ternary atom in an Al site significantly decreases the ICRSS of the SISF and twinning. Out of six total cases (three substitution atoms in the Ti or Al site), only four Figure 5. Two possible models for twinning: (a) the ternary alloying element X (Cr, Nb, or V) is located in a single slip plane and (b) the ternary alloying element X is located between both slip planes. Table 1. Relaxed (without the parentheses) and unrelaxed (in the parentheses) point defect formation energies of each ternary alloy considering the substitution of a Ti or an Al site by an alloying element X (Cr, Nb, or V). Results for supercells with 32 and 108 atoms are presented.

Point defect formation energy (eV)
Cr Nb V www.nature.com/scientificreports/ cases were consider by excluding two cases (Nb in Al sites, V in Al sites) that are not likely to occur in realistic experiments on gamma TiAl crystals 17,18 . As mentioned in the previous section (Section 2-B), the ICRSS of the SISF is notably smaller than the ICRSS of the CSF for all the gamma TiAl alloys considered in the study. Hence, the ICRSS of the SISF is likely related to the yield strength of polycrystalline gamma TiAl. Because ternary atoms would occupy the Ti sites and the ICRSS of SISF slip increases in the presence of ternary atoms (as depicted in Table 3 and Fig. 6), our combined computational and theoretical analysis indicate that the yield strengths of ternary alloys are likely to increase. Our prediction is consistent with previous experimental studies 10-14 that report strengthening due to the ternary atom addition. We note that the SISF energy of the ternary alloy is similar to the SISF energy of the stoichiometric gamma TiAl alloy ( Fig. 6 and Table 2); therefore, the results confirm that the energy barrier and the maximum slope of the SISF are more relevant criteria than the specific value of the SISF energy. www.nature.com/scientificreports/ In contrast, for substitutional defect formation in Al sites, the ICRSS of SISF slip decreased by 30-47%, and the ICRSS of twinning decreased by 32-60%, as presented in Table 3 and Fig. 7. In addition, the ICRSS of forward super slip also decreases. For the quantitative analysis of the propensity for deformation twinning, the twinnability factor (T.F.) is defined as follows: This is the product of the ratio of the ICRSS of the CSF ( σ CRSS CSF ) and the SISF ( σ CRSS SISF ) with the ratio of the ICRSS of forward super slip ( σ CRSS FSuper ) and twinning ( σ CRSS Twin ). If this factor is high, twinning is generally preferred   www.nature.com/scientificreports/ over dislocation slip. Figure 8 shows the effect of the substitution of a Ti or an Al site by an alloying element (Cr, Nb, or V) on the twinnability factor. The substitution of the Ti site exhibits only a marginal difference in the twinnability factor compared to the stoichiometry TiAl alloy, while the substitution of the Al site exhibits a more significant difference. It was found that ternary substitution in an Al site would result in an increase in the twinnability factor. However, because only a small portion of the Cr is expected to occupy the Al sites instead of the Ti sites, it is not straightforward to discuss the effect of the Cr on the plastic deformation mechanism of gamma TiAl crystals.
In addition, to visualize and confirm the ternary effect, the ideal plastic deformation mechanism was predicted from the ideal shear strength comparison method 22 for Cr addition, which occupies both the Ti and the Al sites. First, the [ 11 − 2 ] vector was set as x-axis on the ( 111 ) slip plane, and the shear stress applied on the slip plane from 0 to 180° counterclockwise to the x-axis was examined. The predicted deformation behaviors and the ideal yield shear strength of the stoichiometric L1 0 TiAl and the ternary TiAl-based alloys are presented in Fig. 9a,b, respectively. In Fig. 9b, the blue region represents conditions for dominant SISF formation, while the red region represents conditions for dominant CSF slip, and the intermediate white region represents transient conditions affected by the ternary elements. As shown in Fig. 9b, the substitution of the Ti sites by Cr atoms facilitates the forward super slip with an increased yield shear strength. Conversely, the substitution of the Al sites by Cr atoms facilitates twinning with a decreased yield shear strength. A schematic that summarizes the expected deformation behavior is shown in Fig. 9c.
To summarize the discussion, the GSFE calculations show that a ternary atom substituted for Ti increases the yield strength, while a ternary atom substituted for Al increases the ductility. In other words, the GSFE change pattern is depending on the Al content in TiAl and the kind of the ternary atom affects the magnitude of GSFE changes. These results are consistent with previous experimental papers. 13,14 Our prediction on the yield strength due to a ternary atom occupying a Ti site is consistent with existing experimental studies 10-14 that report a strengthening effect induced by ternary alloying elements. In contrast, from GSFE surface analysis, it was found that all three ternary atoms in the Ti site would have a negligible effect on the ductility, while a ternary atom in the Al site may improve the ductility. However, because the portion of ternary atoms in the Al sites is significantly smaller than the ternary atoms in the Ti sites, the ternary effect on ductility in gamma TiAl phase is expected to be insignificant. Our results are consistent with the existing experimental studies 17,18 , where the ductility of pure gamma TiAl alloy does not change significantly upon the addition of ternary alloying elements. Our study indicates that, to explain the ductility increase in the duplex structure by ternary atoms, it is necessary to consider not only the gamma phase but also the interfacial effects between the gamma phase and the alpha phase as well as microstructural effects, such as grain size and lamellar size and spacing.

Conclusions
We calculated the point defect formation energy and the GSFE of TiAl-based ternary alloys to analyze the changes in the mechanical properties by adding a third element, such as Cr, Nb, or V, to the gamma TiAl. First, from the point defect formation energy calculation, we found that all ternary substitutional atoms are more stable in the Ti site than in the Al site. In the case of Cr, the point defect formation energy difference between the Ti and Al sites is comparable to the thermal energy at ambient conditions. Our calculations are consistent with experimental reports that Nb and V substitutional defects occupy Ti sites, while Cr substitution is found in both Ti and Al sites. Next, we calculated the ternary-atom-induced change in the GSFE and the ICRSS for each partial dislocation via the maximum slope of the GSFE curve. The results suggest that ternary atoms occupying Ti sites increases the yield strength of the alloy because of the increase in the ICRSS of the SISF. In contrast, Figure 8. Calculated twinnability factor of each ternary alloy considering the substitution of a Ti (X Ti ) or Al (X Al ) site by an alloying element X (Cr, Nb, or V).

Scientific Reports
| (2020) 10:21614 | https://doi.org/10.1038/s41598-020-77891-5 www.nature.com/scientificreports/ Cr substitution in Al sites decreases the ICRSS of the SISF and twin, increasing the propensity for twinning. However, the effect on the ductility of the gamma TiAl alloy is limited because a significantly larger portion of the Cr atoms would occupy Ti sites. In summary, we theoretically explain the previous experimental studies reporting that the yield strength increases due to substitutional defects, while the effect on the ductility of the single phase gamma is limited. Our study indicates that further study is necessary to deepen our understanding on the ductility improvement of the duplex structure. We note that the framework used in the present study can be applied to investigate the ternary effect on a variety of intermetallic systems. www.nature.com/scientificreports/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creat iveco mmons .org/licen ses/by/4.0/.