Effects of Alloying on the Interface Energy of the γ′′\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma ''$$\end{document}-Phase in Nickel-Based Superalloys

Nickel-based superalloys of the 718 type have been extensively used in the production of engine parts which operate at high temperatures. The alloy is strengthened by the precipitation of small γ′′\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma ''$$\end{document}-phase particles with a large lattice distortion in the γ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma $$\end{document} matrix. These γ′′\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma ''$$\end{document}-phase particles, however, may transform to the more thermodynamically stable δ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\delta $$\end{document}-phase. Such transformation might be mitigated by adding other alloying elements. In this work, we used density functional theory calculations at 0 K to study the influence of relevant alloying elements on the interface energy of the γ′′/γ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma ''/\gamma $$\end{document} coherent interfaces. By replacing an atom either in the γ′′\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma ''$$\end{document}-phase or at the interface and calculating the corresponding interface segregation energy, we assessed the favorable replacement position of each alloying element. Implications of our results on the γ′′\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma ''$$\end{document}-phase growth are also discussed.


I. INTRODUCTION
DUE to their high strength at elevated temperatures, wrought nickel-based alloy-718-like superalloys have been widely used in the manufacture of highly loaded parts of aircraft engines. [1][2][3][4] The excellent strength is due to the precipitation of the strengthening c 00 -Ni 3 Nb-phase [3] which usually takes the form of disk-shaped particles in the c matrix. Because of the orientation relationships ð001Þ c 00 jjf001g c and ½100 c 00 jjh100i c , c 00 -phase particles usually have (001) and (100) interfaces with the c matrix. [5] At temperatures above 650 C, the c 00 -Ni 3 Nb phase (lattice structure D0 22 ), however, transforms to the d-Ni 3 Nb phase (lattice structure D0 a ) which is more thermodynamically stable. Since the d-phase does not confer the high strength provided by the c 00 -phase, other alloying elements are usually added to stabilize the c 00 -phase. [6] Previous studies using density functional theory (DFT) have reported the effect of aluminum and transition group elements on the stability of the c 00 -phase in the c-matrix. [7,8] Of the elements relevant to c 00 -strengthened alloys (Al, Co, Cr, Hf, Mo, Ta, Ti, W, and Zr), only Co has a lower energy on a Ni site than on a Nb site in the c 00 -phase, but its energy is still lower when it is in the c matrix. [7] Other alloying elements prefer to replace a Nb atom. Moreover, while Ti, Al (c 0 -phase strengtheners), and Cr can be be used for stabilizing c 00 -phase, both Mo and W, which are usually used as solid solution strengtheners in superalloys, [9] prefer the d-phase over the c 00 -phase.
Recently, we have studied the interface energy of c=c 00 using DFT. [10] With the magnetic effects included in the calculations, the interface energies of the c=c 00 001 and 100 interfaces were put at $ 48 and $ 24 mJ=m 2 , respectively, which are within the range of the dynamic interface energy reported by Moore et al. [11] and are, however, smaller than the experimental value given in Reference 12 Moreover, switching off the magnetic effects (by using non-polarized calculations) led to negative values of the interface energy. We note that any comparisons between these results are to be taken with caution since only fully coherent pure c=c 00 interfaces were investigated in the mentioned DFT study, which are very different from experimental semi-coherent interfaces where other alloying elements exist in the c matrix as well as in the c 00 -phase precipitates.
The effects of alloying elements on the interface energy and phase stability of Ni-based superalloys have been reported extensively for variations of interfaces. For example, Xue et al. [13] found that solute segregation can lead to a reduction in grain boundary energy and, thus, stabilize grain boundaries in nickel-based superalloys. Using DFT calculations to investigate cÀNi=c 0 ÀNi 3 Al interface, Ahmed et al. [14] also showed that many alloying elements improve the stability of the interface and increase its fracture strength. For the c=d and c=g interfaces, the effects of alloying on the interface BAO energy have also been reported [15] recently, together with a model to relate the relevant interface quantities to the solution energy and elastic energy.
In this study, we carried out DFT calculations at 0 K to study the influence of the relevant alloying elements on the c=c 00 interface energy. These calculations allow us to see whether these alloying elements prefer to stay in the c 00 bulk phase or at the interface. We also discuss the implication of our results, in connection with previous studies, [7,16] on the design of alloys where c 00 -strengthening is of interest.

II. METHODS
Here, we study the influence alloying elements on the interface energy of both coherent c=c 00 interfaces parallel to ð001Þ c 00 and ð100Þ c 00 surfaces (referred to as 001 and 100 interfaces, respectively). For each orientation, an initial interface supercell is created by stacking two individual c and c 00 phases along the direction normal to the interface (Figure 1), using lattice constants from [7]. The in-plane lattice constants of the two constituents of the initial supercell were averaged. These supercells contain 128 atoms.
For each interface, we replaced an atom in the supercell by an alloying element, either at the center of the c 00 phase or at the interface. As mentioned above, we studied the following relevant alloying elements: Al, Co, Cr, Hf, Mo, Ta, Ti, W, and Zr. We used results reported in [7] to choose the replacement site. Consequently, except for Co which replaced a Ni atom (marked b in Figure 1), all alloying elements replaced a Nb atom (marked a in Figure 1).
All DFT calculations were carried out using the VASP program [17][18][19] with the projector augmented wave method potentials. [20,21] Similar to a previous study, [7] we used the pseudo-potential with maximum number of electrons for each alloying element. The energy cutoff was set to 520 eV which is at least 1.5 times the largest energy cut off of the used potentials. The precision parameter was set to 'accurate.' The real-space projection operators were calculated to the precision of 10 À4 eV=atom. The electronic and ionic convergence criteria were set to EDIFF ¼ 4 Â 10 À7 eV and EDIFFG ¼ 4 Â 10 À5 eV, respectively.
For each alloying element in a particular supercell, we used the same procedure as in 7. First, we relaxed the ionic positions and cell shape at fixed supercell volume for different volumes and use a Birch-Murnaghan equation of state [22] to get the optimal scale factor. A final calculation was then carried out for the optimal scale factor with a k-point spacing of 0.1 Å . All calculations were spin polarized where c-phase atoms were initialized with a magnetic moment of 0:8 l B . The magnetic moments of c 00 -phase atoms were initialized to 0:05 l B , except for Co and Cr whose magnetic moments were initialized to a very large value of 2 l B . In all simulations, c and c 00 phases converged to ferro-magnetic and non-magnetic states, respectively.

III. RESULTS AND DISCUSSION
The main results of this study are plotted in Figure 2 and given in Table I (with the alloying elements sorted by their atomic number) which show the interface segregation energy, E seg , defined as the difference in the supercell energy when the alloying element replaced an atom at the interface and in the bulk, respectively. A positive E seg means that the alloying element prefers to stay in the c 00 -phase bulk. Conversely, that element prefers to segregate towards the interface if E seg <0.  seg and E 001 seg , defined as the difference in the supercell energy when the alloying element replaced an atom at the interface and in the bulk, respectively. A negative value means that the alloying element prefers to segregate towards the interface. Vice versa, the alloying element prefers to stay in the c 00 -phase bulk if its E seg is positive.
Note that E seg of a particular alloying element is also the difference in the interface energy (density) when the element moves from the c 00 -phase bulk to the interface. Figure 2 does not hint to any clear trend of the E seg . Take for example the 100 interface: while E 100 seg <0 for Mo and W and E 100 seg is almost 0 meV for Ta, the remaining E 100 seg values are positive and of the same order. Such an observation, however, does not occur in the E 001 seg of the 001 interface. It has been found that for many dilute transition-metal alloys, the surface segregation energies of single-transition-metal impurities can be related to their d-band filling [23,24] -though such a relation has not been confirmed in Ni(111) surface [25] or in bimetallic alloy surfaces. [26] Here, we also plotted the segregation energy versus d-band filling (Figure 3), with the idealized d-band filling taken from Reference 27 (Figure 3 therein). Our current dataset, however, does not reveal any relation between the interface segregation energies E seg and d-band fillings. For example, while the E 100 seg values of the investigated 3d elements are more or less constant, E 100 seg of 4d and 5d elements seems to decrease with respective to d-band occupancy.
Our results show that the elements Mo and W have the tendency to move to the 100 interface: their E seg are negative for both interfaces with jE 100 seg j>jE 001 seg j. Note that their diffusion in nickel is slow. At 900 C, their interdiffusion coefficients with nickel are only around 10 À17 m 2 =s (as a comparison, the diffusion coefficients of Ta and Hf in nickel, for example, are more than one order of magnitude larger). [28,29] Together with our current results, this means that adding W would hinder the growth of c 00 -phase particles in the 100 direction: Any growth of the c 00 -phase particles would require the interface to carry these alloying atoms along. Moreover previous studies [7,15,16] have already shown that, in terms of the stabilization of the c 00 -phase, it is not advisable to add Mo and W since they possibly destabilize the c 00 -phase in favor of the d-phase.
Ta is another element whose values of E seg are negative for both interfaces. However, these values are small (less than 10 meV) and thus Ta basically does not have any preferable positions and adding it will not affect the growth of c 00 -phase particles.
The c 0 -forming elements Al and Ti prefer to be at the 001 interface (E 100 seg >0 and E 001 seg <0). Thus, having a 001-interface bordering on a region rich of (Al, Ti) might be energetically favorable. This is in agreement with the experimental results which show composite c 00 =c 0 precipitates in the c matrix at low temperatures. [30] Co also prefers to be at the 001 interface (E 100 seg >0 and E 001 seg <0). However, Co has lower energy when it is in the c matrix [7] and, thus, it will not stay in the c 00 -phase but rather move to the c phase and diffuse away from the interface.
The remaining alloying elements (Cr, Zr, Hf) all prefer to stay in the bulk c 00 -phase (E 100 seg >0 and E 001 seg >0). Here, we note that while the Cr diffusion in nickel is slow (its diffusion activation energy is about 300kJ=mol which is on par with the diffusion activation energy of W in nickel), the diffusion of Zr and Hf in nickel is much faster. [29] Adding Cr will, thus, make the c 00 -phase growth in the 100 direction less energetically  favorable than in the 001 direction (its E 100 seg is much larger than its E 001 seg ). Conversely, adding Zr and Hf will probably not suppress the growth of the c 00 -phase in any directions. However, the contents of Zr and Hf must be limited since increasing their contents might hinder the grain coalescence at grain boundaries [31] or lead to the formation of chains of misoriented grains which causes the deterioration of creep resistance. [32] The present study was based on spin-polarized calculations with which the magnetic effects were included. As mentioned above, switching off the magnetic effects led to a negative c=c 00 interface energy. [10] Similar influence of the magnetic effects has been also reported previously by Woodward et al. [33] where the authors found that non-polarized calculations led to a negative interface energy of the cÀNi=c 0 ÀNi 3 Al interface. One, thus, expects that the results of the current study might change when the influence of magnetic effects is altered either by switching off all these effects altogether or by introducing other elements with strong magnetic interactions, such as Fe, Co, and Cr, into the c matrix.
Finally, we also emphasize that obvious simplifications of the c=c 00 interface have been made to ease the computing difficulties: Here, the interface is assumed to be fully coherent and pure-e.g., the c phase contains Ni only. The existence of other elements in the alloys, as in experiments, might lead to changes in, for example, elastic energy (due to difference in sizes) and magnetic interactions (see above) and, thus, might lead to changes in the segregation energy.

IV. SUMMARY
Using density functional theory at 0 K , we have calculated the interface segregation energy for c=c 00 interface for some alloying elements which are relevant to c 00 -phase strengthening in nickel-based superalloys. The implications of these results were then discussed: For example, while Mo and W should be avoided since these elements might hinder the growth of c 00 precipitates, addition of Al, Ti, and Cr might change the shape of c 00 precipitates. Further studies are necessary to assess the influence of magnetic effects on the current results.

ACKNOWLEDGMENTS
This work was supported by the German Research Foundation (DFG) through the project BA-1795/13-1. Computing resource was made available by HLRN through the compute project nii00170.

FUNDING
Open Access funding enabled and organized by Projekt DEAL.

CONFLICT OF INTEREST
On behalf of all authors, the corresponding author states that there is no conflict of interest.

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 ivecommons.org/licenses/by/4.0/.