arXiv:1306.1903v1 [cond-mat.mtrl-sci] 8 Jun 2013 nteceia nuty hl,teeaeconsiderable are there While, as well industry. air- chemical as in engineering, the applications aerospace in and important construction many craft has titanium gen, ex- [9–14]. the reach where cannot quantitative processes, measurements specific atomistic perimental microscopic the to of in the rise involved study give track values the and to in processes help diffusion used can widely diffusion, been solid-state have which sim- dynamic ulations, applications. molecular practical and and calculations First-principles understanding impor- quite basic in are for which involved tant steps, diffusion processes and de- and physical formation bubbles cannot the microscopic helium they the of generally termine experiments coefficients, energy the diffusion formation Although overall the metals give [1–8]. in can diffusion presented helium been and investigations have growth analytical bubble and helium experimental of Many 3]. swelling [2, and bubbles cave, embrittlement, helium Particu- intergranular temperatures, cause can homologous de- metals. high substantially of at may properties larly, which mechanical site, the migrate trapping teriorate atoms same helium the several to when agglomerate bubbles to implanta- helium tend they into direct transmutation, by As neutron or site. increases tion trapping concentration a helium reaches in it rapidly the until diffuse extremely region atoms has interstitial Helium the and metals. metals in to config- solubility electron inactive low closed-shell is the de- helium to radioactive Due uration, the [1]. via as tritium or of such ions cay particles light and energetic neutrons of reactor reactions nuclear through zhang ∗ uhrt hmcrepnec hudb drse.E-mail: addressed. be should correspondence whom to Author sa motn aeilt tr n ereehydro- retrieve and store to material important an As either metals in produced be could atoms Helium is-rnilssuyo eprtr-eedn diffusio temperature-dependent of study First-principles [email protected] 1 C,Isiueo ple hsc n opttoa Mathem Computational and Physics Applied of Institute LCP, ASnmes 32.k 32.- 66.30.-h 63.20.D-, diff 63.20.dk, the numbers: confines PACS and energy fo total the the that reduces found site is it octahedral Also, anisotropy. remarkable showing h otltie h lmigiaende lsi ad(CI str the point band saddle within elastic associated coefficients nudged and diffusion image pathways climbing vibra energy first- minimum The the the from between lattice. coupling obtained host full are the the coefficients including diffusion calculations impurity paramet the microscopic The of theory. state transition the using 2 h eprtr-eedn iuincecet finterst of coefficients diffusion temperature-dependent The ejn opttoa cec eerhCne,Biig10 Beijing Center, Research Science Computational Beijing .INTRODUCTION I. ogLu, Yong xy 1 ln ( plane ae Zheng, Fawei D xy α sawy ihrta htaogthe along that than higher always is ) -Ti 1 where where ee 1] h iuincecetcnb rte as written be can coefficient diffusion the [15], Zener k as expressed ntedrcino diffusion, of direction the in atom, interstitial diffusing the aresaddffso offiinsfo rtprinciples. first energy from corresponding coefficients the diffusion of and values barriers specific the obtain to ege hs rniinfrom transition phase dergoes ntemcaia rpris tro eprtr,ti- temperature, hexagonal room the At exhibits tanium properties. mechanical deterioration severe the a in in resulting bubbles into glomerate rO odt,teeprmna iuincoefficients diffusion experimental C, the in H, helium date, as To of such elements O. light or are other metals N to in contrast in behavior few diffusion very its the about within knowledge out the carried solely diffusion α be makes to Ti difficult rela- for This experiments temperature K. transition 1155 low around tively temperature the at structure) nweg,systematic knowledge, n rtu eoe oe with contain- doped titanium becomes then tritium decay, ing tritium into introduced through usually titanium are atoms Helium hydrogen-containing environments. under applications its in restrictions uinin fusion h tmcdffso ehns fhlu in helium of mechanism diffusion investigate comprehensively atomic to the is work present the of pose n igZhang Ping and pae utemr,det h nouiiyo helium, of insolubility the to due Furthermore, -phase. steBlzancntn.FloigteWr and Wert the Following constant. Boltzmann the is is coefficient diffusion the form, general most the In tc,Biig108,Pol’ eulco China of Republic People’s 100088, Beijing atics, r ntepefco n ciainenergy activation and pre-factor the in ers n Q inlmdso h iuigao and atom diffusing the of modes tional mto fhlu ie etrda the at centered dimer helium of rmation ta eimao in atom helium itial so fhlu atoms. helium of usion stenme fnaetniho tbestsfor sites stable nearest-neighbor of number the is steatvto energy, activation the is α cue.I sdmntae htthe that demonstrated is It uctures. rnilsttleeg n phonon and energy total principles I HOYADMETHODS AND THEORY II. E)mto sue osac for search to used is method NEB) 04 epesRpbi fChina of Republic People’s 0084, ii tl akn nteltrtr.Tepur- The literature. the in lacking still is Ti α iaimaerte cre[] n oour to and [4], scarce rather are titanium ,2, 1, offiinsfrhlu in helium for coefficients n D ∗ D = binitio ab = D α 0 nβd e T r predicted are -Ti d − Q/kT β stelnt ftejump the of length the is 2 α α Γ z tde fhdoe dif- hydrogen of studies 3 paeto -phase stejm probability jump the is , D e hc ed oag- to tends which He, hs,wees tun- it whereas, phase, xs( axis , 0 sapefco,and pre-factor, a is D z ), β pae(bcc -phase α T,and -Ti, (1) (2) 2 projected onto the direction of diffusion, and Γ is the jump rate between adjacent sites of the diffusing atom. According to the transition state theory (TST) [16, 17], the jump rate is written as

3N−6 0 kT 1 − exp(−hν /kT ) − Γ= i=1 i e ∆Hm/kT , (3) h Q3N−7  ∗  i=1 [[1 − exp(−hνi /kT )]] ∗ Q 0 where νi and νi are the vibrational frequencies at the transition state and the ground state, respectively. For the high temperature range (hνi/kT ≪ 1), all vibrational degrees of freedom approach a classical behavior and the jump rate reduces to, FIG. 1: (Color online) Visual structure of α-Ti matrix marked 3N−6 0 i=1 νi −∆Hm/kT with tetrahedral (T) site, octahedral (O) site, the center of Γ= − e . (4) Q3N 7 ν∗ equilateral trigonal face of octahedron (F) site, and the cen- i=1 i ter of equilateral trigonal face of tetrahedral (M) site. Large Q For the low temperature range (hνi/kT ≫ 1), the jump green and small red balls represent for titanium atoms and rate can be expressed as [10] interstitial sites for helium atom, respectively.

kT − Γ= e (∆Hm+∆Ezp)/kT , (5) h where Ezp is the difference in zero-point energies between ground state and transition state. Using the phonon free energy at different states ∞ hν Fvib = kT g(ν)ln 2sinh , (6) Z0  2kT  the jump rate can be simply expressed as

kT − − Γ= e ∆Fvib/kT e ∆Hm/kT , (7) h where the zero-point energy is included in the Fvib term. The density functional theory (DFT) calculations are carried out using the Vienna ab-initio simulation pack- age (VASP) [18, 19] with the projector-augmented-wave (PAW) potential method [20]. The exchange and correla- FIG. 2: (Color online) Contour plots of charge density dis- tion effects are described by the local density approxima- tributions of pure Ti matrix inside the (110) plane crossing T-site (a), O-site (b), F-site (c), and M-site (d), respectively. tion (LDA). The cutoff energy for the plane-wave basis set is 450 eV. We employ a 3 × 3 × 2 α-Ti supercell con- taining 36 host atoms to simulate helium migration in the III. SOLUTION OF HELIUM ATOM α-Ti matrix. To check the convergence of the formation energies, we have also considered a 4×4×3 α-Ti supercell containing 96 host atoms and one helium atom. The in- The formation energy of interstitial helium atom in tegration over the Brillouin zone is carried out on 5×5×5 titanium matrix is calculated by k-point meshes generated using the Monkhorst-Pack [21] method, which proves to be sufficient for energy conver- Ef (He) = E(Ti,nHe) − nµ − E(Ti), (8) gence of less than 1.0×10−4 eV per atom. During the where E(Ti,nHe) is the total energy of the supercell con- supercell calculations, the shape and size of the super- taining n interstitial helium atoms, µ is the chemical po- cell are fixed while all the ions are free to relax until the forces on them are less than 0.01 eV A˚−1. In order to cal- culate the migration energies of helium in titanium, each saddle-point structure and associated minimum energy TABLE I: Calculated formation energy (Ef ) for helium atom pathway (MEP) were calculated by employing the image in various solution sites in hcp-type Ti matrix. nudged elastic band (CINEB) method [22]. The vibra- T O F M Ef Ef Ef Ef tional frequencies were computed using the density func- 36-atom cell 2.501 2.746 2.385 2.677 tional perturbation theory with the help of PHONOPY 96-atom cell 2.499 2.747 2.382 2.675 package [23]. 3

(a) 0.4

O-site (TS)

(b) 0.3 TABLE II: Interaction energy (Eint) and deformation energy

0.2 (Edef ) for helium atom in different solution sites.

H = 0.358 eV

0.1

F-site F-site

0.0 F-site O-site T-site M-site

Migration enthalpy (eV)

0 1 2 3 4

Reaction coordinate Eint(He-Ti) 1.952 2.169 2.074 2.130

Edef (Ti) 0.432 0.578 0.427 0.547

TS (c) 0.6

0.4 M-site

1

H = 0.637 eV

0.2

3

F-site

0.0

tion with the lowest charge density is the F-site, and the Migration enthalpy (eV) c 0 1 2 3 4

Reaction coordinate

2 M-site occupies the highest charge density. The charge

0.3

TS Ti (d) density in O-site is somewhat higher than that in F-site,

0.2

F-site which can be further confirmed by comparing the refer- H = 0.262 eV

1

0.1

H = 0.145 eV T-site T-site 2 ence isosurface (not shown here). From the viewpoint of F-site

0.0

O-site

Migration enthalpy (eV) electron distribution, F-site is the most favorable site for 0 1 2 3 4 5 6

a Reaction coordinate helium atom, since the interstitial light elements tend to locate at the site with lowest electron density [28]. Fur- FIG. 3: (Color online) Visual structure of three different jump thermore, we calculated the interaction energy between paths characterized by their jump rates Γ1,Γ2, and Γ3 for the He (a), and the reference energy profiles for the He to diffuse interstitial helium atom and its surrounding host atoms along (b) F→O→F, (c) F→M, and (d) F→T paths. according to the following expression ∗ Eint(He − Ti) = E(Ti, He) − E(Ti) − E(He), (9) tential of helium atom, E(Ti) is the energy of the Ti where E(Ti, He) is the total energy of optimized Ti super- supercell without interstitial atoms. cell containing one He atom, E(Ti)∗ is the total energy The favorable solution sites for helium atom in met- of the supercell containing only Ti atoms with the same als are those positions of high-symmetry or low charge atom position and cell parameters as those of E(Ti, He), density. From an intuitional consensus, the octahedral and E(He) is the energy of one isolate helium atom. The site (O-site) with six nearest Ti atoms and the tetrahe- reference results are organized in Table II. Clearly, the dral site (T-site) with four nearest Ti atoms are strong interaction energy Eint in F-site is the lowest, while the candidates, as marked in Fig. 1. In addition to O-site highest energy appears at the O-site. By comparing the and T-site, there are two other interstitial sites in the energy, the F-site is still the most favorable site for he- hcp-type Ti matrix, i.e., the center of equilateral trigo- lium atom. The interaction between interstitial helium nal face shared by two adjacent O sites (F-site) and the atom and host atoms will cause the deformation of Ti center of equilateral trigonal face shared by two adjacent lattice, further increasing the formation energy. The de- T sites along z axis (M-site), as depicted in Fig. 1. We formation energy can be evaluated by calculated the formation energies of one helium atom in ∗ these four types of solution site using the 36-atom model Edef (Ti) = E(Ti) − E(Ti), (10) and 96-atom model respectively, and the results are or- ganized in Table I. Clearly, these two models show the and the results are also listed in Table II. Due to the similar convergence in energy, implying that the 36-atom strong interaction between He and Ti atoms, the O-site supercell is sufficient to assume the simulation task. has the largest deformation energy. The F-site and T-site The positive formation energies indicate that the dis- have low deformation energy, coinciding with the rela- solution of He in Ti is endothermic. The formation en- tively weak interaction between He and host atoms in ergy of helium in titanium matrix (ranging from 2.385 these two sites. eV to 2.746 eV) is much lower than that in bcc-type Fe (ranging from 4.22 eV to 5.36 eV) and fcc-type Ni metals (around 4.50 eV) [24–26]. The favorable site in energy IV. DIFFUSION IN TITANIUM for helium atom in titanium matrix is the F-site. Unex- pectedly, the O-site, supplying the maximum space for From the calculation of formation energy above, we helium atom among the four interstitial sites, has the have confirmed that the most favorable site for helium highest formation energy. This behavior of helium is dif- atom in α-Ti matrix is the F-site. In the following, we ferent from that of H atom in titanium, where O-site is will investigate the diffusion paths for interstitial helium more favorable in energy [27]. To a great extent, the elec- atom. tronic distribution and the interaction between He atom As illustrated in Fig. 3 (a), we set the F-site as initial and its surrounding Ti atoms make main contribution in site and choose three candidate jump types characterized formation energy. Thus, we first plot the charge density by their corresponding jump rates Γ1, Γ2, and Γ3. A dif- of (001) planes crossing the four different sites for pure fusive jump of Γ1 connects two nearest-neighbor F sites Ti matrix, as shown in Fig. 2. By comparing the four in z-axis direction. Γ2 and Γ3 are two jump types of he- contour plots of charge density [Figs. 2(a)-(d)], the posi- lium atom within the basal xy plane, where Γ2 connects 4

TABLE III: Calculated electronic energy Eel and zero-point energy Ezp for interstitial helium atom in 36-Ti matrix. All energies are in kJ/(mol of Ti36(He)). The electronic energy of the elements in their standard states, i.e., Ti36, is taken as reference.

Ti36 Ti36He(F) Ti36He(T) Ti36He(O) Ti36He(M) Ti36He(TSFT) Ti36He(TSFM) Eel 0 229.63 240.80 264.49 257.78 254.87 290.99 Ezp 119.39 124.98 123.11 113.59 121.82 122.73 121.86

xy z TABLE IV: Calculated diffusion pre-factors (D0 and D0 ) and activation energy (Qxy and Qz) for helium atom in α-Ti. The temperatures represent the ranges over which diffusion coefficients were fit to extract corresponding pre-factors and activation energies. For comparison, experimental and molecular simulation results are also listed.

xy 2 z 2 Method D0 (m /s) D0 (m /s) Qxy (kJ/mol) Qz (kJ/mol) T (K) − − This study 5.06×10 6 6.63×10 7 24.48 23.16 280-850 − − MDa 2.1×10 8 2.6×10 7 17.34 17.34 298-667 a Reference [12]

Temperature (K)

20 19.08 1000 800 600 400 300

20.45 (double)

17.53 -6

16.62 (double) 10

15 Ttr = 1155 K

-7 D xy 10

10

D z

5 ) -1

-8 s Frequency (THz)

2 10

0

M A M A M A M A M A

-1.44

-5

-6.18

-9

-6.83

10

-10

M-site F-site T-site O-site TS-FT Diffusivity (m

-10

FIG. 4: Calculated phonon dispersions of a 36-atom Ti su- 10 percell with helium atom in the F-site, T-site, O-site, M-site and the transition states along the F→T (TS-FT) and F→M

(TS-FM) paths, respectively. The modes with imaginary fre- -11

10

quencies at the transition states correspond to the motion of 1.0 1.5 2.0 2.5 3.0 3.5 the helium atom across the barriers. 1000/T (1/K)

FIG. 5: (Color online) Calculated diffusion coefficients of atomic helium in α-Ti. The Ttr=1155 K stands for the phase these two nearest-neighbor F sites directly and Γ3 via transition temperature of metal Ti from α-phase to β-phase. the T site. We shall determine the more favorable jump type between Γ2 and Γ3 by comparing their migration enthalpy. To estimate the migration enthalpy for helium ones of F→M for Γ2 path and F→T for Γ3 path respec- atom diffusion between interstitial sites, the three paths tively, as shown in Figs. 3(c) and 3(d). Evidently, the are evenly divided into several segments. We placed He migration enthalpy of the Γ3 path (0.262 eV) is much atom at the appropriate position every equal distance. lower than that of Γ2 (0.637 eV). That is, the helium By employing the CINEB method, the activation energy atom tends to diffuse between F sites through the Γ3 profiles vs diffusion distance are obtained and shown in path within the basal xy plane. Figs. 3(b)-(d). For the diffusion path Γ1, the energy is The lattice vibrational energy also plays a great role in found to display a single maximum, corresponding to a the migration energy, especially in the high temperature saddle point at the high-symmetry position located half range. In the following, we will calculate the phonon dis- way between neighboring F sites, i.e., the O site acts as persions of these different states using the density func- a transition state role. The energy profile can be well de- tional perturbational theory, and obtain the reference lat- scribed by a sinusoidal curve with an migration enthalpy tice vibrational energy. The 3×3×2 Ti supercell matrix of 0.358 eV. For jump Γ2, we find that the middle point of is employed and one helium impurity is placed at the diffusion path coincides with M site when the relaxation F-site, T-site, O-site, M-site, and the transition states process completed. Actually, the jump Γ2 thus connects of F→T (TS-FT), respectively. The phonon dispersions two nearest-neighbor F sites via the M site. Due to the of all these atomic structures including the pure 36-Ti symmetry of the diffusion curves, we plot here only the matrix are calculated and shown in Fig. 4. Obviously, 5 when a helium atom occupies the F site, the He-related A linear fit between 280 K and 850 K of the calculated phonon dispersions have the highest frequencies, which data gives Qxy=24.48 kJ/mol and Qz=23.16 kJ/mol, show a double degenerate dispersionless branch at 20.4 both somewhat higher than the molecular dynamic sim- THz. The F site locates at the geometric center of trian- ulation results of 17.34 kJ/mol [12] obtained using the gle formed by its three nearest Ti atoms, thus the inter- Einstein relation [29]. The diffusion coefficients of these stitial space is much more confined with respect to other two directions are very anisotropic, corresponding to the xy −6 2 z −7 states. At the T site, the highest He-related frequency pre-factors of D0 =5.06 ×10 m /s and D0=6.63×10 decreases to 19.1 THz and the degeneracy vanishes. In m2/s, respectively. This is consistent with the qualita- both cases, the frequencies are positive, indicating a true tive observation of helium migration enthalpy, where the minimum in the energy. By definition, the transition helium atom is easier to overcome the migration enthalpy state is characterized by the occurrence of one negative within the xy plane (cf. Fig. 3). Note that the theoretical eigenvalue in the dynamical matrix, showing the imag- diffusion coefficients obtained here should be the upper inary frequency in the phonon dispersion, which corre- bound of helium migration in α-Ti matrix. Due to the sponds to a motion of helium atom along the diffusion existence of vacancies or dislocations in a real crystal, as path. When helium locates at the transition state along well as the formation of helium bubble, the diffusion co- the F→T path, we note that a Ti-related phonon branch efficients will be reduced to some extent. To this end, we between 8.5 and 9.0 THz is obviously off from the bulk also made a relevant test to elaborate the effects of clus- of Ti-related phonon branches. This branch is related to tering behavior of two helium. We set two helium atoms motions of Ti atoms, which couple with motions of the He at the adjacent F sites along the z axis, and then per- impurity. At the octahedral site, the He-related branches form the relaxation of the system until the largest force are all negative, locating at the −6.8 THz nearby. The on any atom is smaller than 0.01 eV/A.˚ We find that two O-site is the point of maximum solution energy and the helium atoms form a dimer along the h001i direction cen- helium atom is unstable at this point, tending to dif- tered at the octahedral site. The formation energy is just fuse away from it. From the phonon dispersions, we can 2.05 eV per He atom, lower than that of one helium case evaluate the zero-point energy, the temperature depen- (cf. Table I). The formation of helium dimer will confine dent enthalpy, and free energy. The electronic energy and the diffusion of helium atom in Ti matrix, reducing the zero-point energy of these states are listed in Table III. diffusion coefficients. The zero-point energy at F site is 124.98 kJ/mol, which is the highest one among all the considered states due to the highest He-related branch. With the correction of zero-point energy, the most favorable interstitial site in energy for helium is still the F site at ground state. Finally, we can confirm that the helium atom diffuse in VI. CONCLUSION Ti matrix by Γ1 jump path along the z axis and Γ3 jump path within the xy plane, respectively. In summary, we have systematically studied the V. DIFFUSION COEFFICIENTS temperature-dependent diffusion coefficients of helium atom in α-Ti using transition state theory with accurate first-principles total energy and phonon calculations. It Under the assumption of the above jump paths and is found that the most stable solution site for He in Ti according to Eq. (2), the two diffusion coefficients in the matrix is the F site, with the lowest formation energy hcp-structured α-Ti can be expressed as of 2.385 eV. Two minimum energy pathways and associ- 2 ated saddle point structures are determined by using the Dxy = a Γ1 (11) CINEB method. Within the xy plane, the helium atom and diffuse between adjacent solution sites by crossing the T meta-stable site, while along the z axis it should pass the c 2 Dz = ( ) Γ3, (12) transition state of O site. The obtained diffusion coef- 2 ficients within the xy plane (Dxy) and along the z axis respectively, where a and c are the lattice parameters (Dz) show remarkable anisotropy, and the helium atom is of α-Ti. Diffusion coefficient presents more quantitative more easily to diffuse within the xy plane. The formation description on the diffusion features. In Fig. 5 we show of helium dimer centered at the octahedral site reduces the Arrhenius plot, i.e., ln(D) = ln(D0) − Q/kT , of the the total energy compared to the case of two isolate he- computed diffusion coefficients Dxy and Dz. As shown lium atoms, which confines the diffusion of helium and in Fig. 5, the diffusion of helium in α-Ti matrix shows further reduce the diffusion coefficients. Without consid- remarkably linear behavior and thus we can obtain the ering the vacancies or helium bubble effects, our theoret- activation energy and prefactor by linear fitting. In Table ical diffusion coefficients should be an upper bound on IV we listed our fitting results of pre-factors and activa- the helium diffusion in Ti matrix, which provides a good tion energies. reference for future experimental measurements on α-Ti. 6

Acknowledgments Science and Technology under Grant No. 2011A0301016.

This work was supported by NSFC under Grant No. 51071032 and by CAEP Foundations for Development of

[1] V. Sciani, P. Jung, Radiation Effects 78 87-99 (1983). 036105 (2011). [2] H. Schroeder and P. Batfalsky, J. Nucl. Mater. 103-104, [15] C. Wert and C. Zener, Phys. Rev. 76, 1169 (1949). 839 (1981). [16] H. Eyring, J. Chem. Phys. 3, 107 (1935). [3] H. Trinkaus, B.N. Singh, J. Nucl. Mater. 323 229 (2003). [17] G. H. Vineyard, J. Phys. Chem. Solids 3, 121 (1957). [4] R. Vassen, H. Trinkaus and P. Jung, J. Nucl. Mater. 183 [18] G. Kresse and D. Joubert, Phys. Rev. B 59, 1758 (1999). 1 (1991). [19] G. Kresse and J. Furthm¨uler, Phys. Rev. B 54, 11169 [5] G.J. Thomas, W.A. Swansiger and M.I. Baskes, J. Appl. (1996). Phys. 50 6942 (1979). [20] P. E. Bl¨ochl, Phys. Rev. B 50, 17953 (1994). [6] B.N. Singh, H. Trinkaus, J. Nucl. Mater. 186, 153-165 [21] H. J. Monkhorst and J. D. Pack, Phys. Rev. B 13, 5188 (1992). (1976). [7] R. Rajaraman, B. Viswanathan, M.C. Valsakumar, and [22] G. Henkelman, B. P. Uberuaga, and H. Jonsson, J. Chem. K.P. Gopinathan ,Phys. Rev. B 50, 597-600 (1994). Phys. 113, 9901 (2000). [8] B. Glam, S. Eliezer, D. Moreno, D. Eliezer, J. Nucl. [23] A. Togo, F. Oba, and I. Tanaka, Phys. Rev. B 78, 134106 Mater. 392 413-419 (2009). (2008). [9] W.D. Wilson, C.L. Bisson and M.I. Baskes, Phys. Rev. [24] W.D. Wilson, R.A. Johnson, in: P.C. Gehlen, J.R. B 24 5616 (1981). Beeler, R.I. Jaffee (Eds.), Interatomic Potentials und [10] E. Wimmer, W. Wolf, J. Sticht, and P. Saxe, Phys. Rev. Simulation of Lattice Defects, Plenum Press, New York, B 77 134305 (2008). p. 375, (1972). [11] M. Mantina, Y. Wang, R. Arroyave, L. Q. Chen, Z. K. [25] C.C. Fu, F. Willaime, Phys. Rev. B 72 064117 (2005). Liu, and C. Wolverton, Phys. Rev. Lett. 100 215901 [26] T. Seletskaia, Y. Osetsky, R. E. Stoller, and G. M. Stocks, (2008). Phys. Rev. Lett. 94, 046403 (2005). [12] M. Chen, Q. Hou, J. Wang, T.Y Sun, X.G. Long, S.Z. [27] Y. Lu, P. Zhang, J. Appl. Phys. 113, 193502 (2013). Luo, Solid State Commun. 148 178-181 (2008). [28] J.K. Nørskov, Phys. Rev. B 26 2875-2885 (1982). [13] M. Mantina, Y. Wang, L. Q. Chen, Z. K. Liu, and C. [29] G. Boisvert, L.J. Lewis, Phys. Rev. B 54 2880 (1996). Wolverton, Acta Mater. 57 4102-4108 (2009). [14] B.L Zhang, J. Wang, and Q. Hou, Chin. Phys. B 20