Theoretical study of the interaction between hydrogen and 4d alloying atom in nickel

Theoretical study of the interaction between hydrogen and 4d alloying atom in nickel 1 corresp first-author Corresponding authors, Tel.: +86 (021) 39194735, E-mail address: liuwei@sinap.ac.cn Corresponding authors, Tel.: +86 (021) 39194735, E-mail address: liuwei@sinap.ac.cn liuwei@sinap.ac.cn 1 1 1 1 1 1 1 1 1 corresp Corresponding authors, Tel.: +86 (021) 39194115, E-mail address: liuwenguan@sinap.ac.cn Corresponding authors, Tel.: +86 (021) 39194115, E-mail address: liuwenguan@sinap.ac.cn liuwenguan@sinap.ac.cn Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Shanghai 201800, China Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Shanghai 201800, China In order to investigate the tritium behaviors in Hastelloy N alloy in molten salt reactor, first-principles calculations are used to study the interaction between hydrogen and 4d transition metal alloying atom in nickel-based alloy. The interaction energies between 4d elements and H are calculated. Atomic size effects and electron distribution effects are analyzed. The hydrogen-4d interactions are compared with the hydrogen-3d interactions calculated in our previous work. hydrogen tritium transition metal nickel first-principles calculation 1. Introduction 1 Hastelloy N, a Ni-based alloy, is the special structural material for molten salt reactor (MSR) [ 1 1 , 2 2 ] with excellent corrosion resistance against molten salts. Tritium (T), produced by neutron reactions with FLiBe molten salt, can flow readily through the Hastelloy N walls of various systems in MSR. With radioactivity, T release to the environment must be controlled. Knowledge of H behavior in Ni-based alloy is essential for predicting T behavior in MSR and guiding T control processes [ 3 3 ]. Also, H behavior in alloys is a fundamental subject for understanding hydrogen embrittlement [ 4 4 , 5 5 ] suffered by alloys. While it is difficult to experimentally explore H behaviors in alloys, first-principles calculations can be used to study H behaviors in alloys and their underlying physical mechanisms at the atomic and electronic levels [ 6 6 , 7 7 ]. Since the alloying additives can directly modify H behavior in alloys, the interactions between them are very important. In our previous first-principles calculations [ 8 8 ], the interactions between H and 3d transition metal (TM) alloying atom in Ni were found to be influenced greatly by the electron distribution, which do not obey the elastic mechanism [ 9 9 ] or the “reverse stability” rule [ 10 10 ]. In this work, we further investigated the interaction between 4d TM atoms and H in Ni with first-principles calculations, and the differences between 4d and 3d elements in their interaction with H. 2. Computational details 1 The computational method is stated in detail in our previous work [ 8 8 ], and spin-polarized electronic state calculations are performed within the density functional theory [ 11 11 , 12 12 ] using Vienna ab-initio simulation package (VASP) [ 13 13 ]. A 2×2×2 Ni supercell containing 32 Ni atoms is used with a Monkhorst-Pack 7×7×7 k-point mesh. In the Ni-4d supercell, a Ni atom is replaced by a 4d atom. In the Ni-4d-H supercell, one H atom is placed in the octahedral or tetrahedral interstitial site (O-site or T-site) nearest to the substitutional 4d atom. Energies of H vibrations are taken into account by calculating the zero point energy (ZPE). To a first-order approximation, only H vibrations are calculated and network vibrations are neglected [ 14 14 ]. The computational method of the interaction energy between 4d element and H is the same as our previous work [ 8 8 ] according to Refs.[ 9 9 ] and [ 10 10 ]. A positive value of the interaction energy means the rejection to H in the Ni-TM-H system, and a negative value refers to the attraction to H. 3. Results and discussion 1 The calculated total energy of the Ni-4d-H supercell with H in O-site or T-site is shown in Fig. 1. Fig. 1. It can be seen that H in O-site is more energetically preferable than in T-site, since O-site is larger than T-site. This is the same as the case of 3d elements [ 8 8 ]. Furthermore, the T-site in the Ni-4d-H supercell for the cases of Zr, Nb, Mo, Ag and Cd is unstable for H atom. Thus, only the cases of H in O-site are taken into account in the following sections. The calculated 4d-H interaction energies are shown in Fig.2(a) Fig.2(a) . In order to validate the 32-atoms supercell, the interaction energies calculated in 108-atoms supercell (denoted by ‘108’ in parentheses) are also shown in Fig. 2(a) Fig. 2(a) . The change trends of interaction energies in 32-atoms and 108-atoms supercells are almost the same. However, the interaction energies calculated in 32-atoms are all positive, and those calculated in 108-atoms are negative (expect Cd). This difference can be attributed to the influence of H concentration (i.e. the supercell size), which is shown in the results of 3d-H interaction [ 8 8 ]. In later text, only the results calculated in 32-atoms supercell are analyzed, which represents the case of a relatively high H and 4d TM atoms concentrations and can be used to compare with our results of the 3d-H interaction conveniently [ 8 8 ]. The TM-H interaction can be decomposed into the mechanical and chemical effects. The mechanical effect is caused by the atomic size effect mechanism, and the chemical effect is induced by the electron distribution. Since the change trends of the bond length and binding energy in Fig. 3 Fig. 3 do not have evident correlation with that of the interaction energy shown in Fig. 2(a) Fig. 2(a) , neither the bond length nor the binding energy can explain the 4d-H interaction energy in Fig. 2(a) Fig. 2(a) alone, as displayed in Fig. 3 Fig. 3 , which is consistent with the results of the 3d-H interaction [ 8 8 ]. For the 3d-H interaction, the interaction energies are all (except Zn) negative ( Fig.3 Fig.3 in Ref. [ 8 8 ]). However, for the 4d-H interaction, the interaction energies are all positive ( Fig. 2a Fig. 2a ). Comparing Fig. 3 Fig. 3 in this work with Fig. 2 Fig. 2 in Ref. [ 8 8 ], the Ni-4d bonds in Ni-TM supercell are longer than the Ni-3d bonds due to the bigger atomic sizes of 4d elements. The bigger 4d atomic size induces a stronger extrusion in the supercell resulting in a stronger repulsive force to H atom. So, the effect of atomic size plays a dominant role in the 4d-H interaction, which is mainly different from the 3d-H interaction. To explore the chemical components, of the electro-negativities of the 4d elements [ 15 15 ] shown in Fig. 2(b) Fig. 2(b) , the calculated charge density and the charge density difference [ 16 16 ] are demonstrated in Figs. 4 Figs. 4 and 5 5 . For Zr, the atomic size is relatively big and close to those of Ag and Cd in Fig. 3 Fig. 3 , so we can expect that Zr should exert a strong repulsion to H. However, the repulsion of Zr is obviously weaker than those of Ag and Cd as displayed in Fig. 2(a) Fig. 2(a) , which can be attributed to the chemical component. Due to the lowest electronegativity of Zr in Fig. 2(b) Fig. 2(b) , Zr loses electrons easily, and the electrons are gained by the Ni-H bonds (e.g. Ni1-H bond in Fig. 4b Fig. 4b ) which become stronger than those in the Ni-H supercell (e.g. Ni1-H bond in Fig. 4a Fig. 4a ). The stronger Ni-H bonds reduce the repulsive interaction in the Ni-Zr-H supercell, which can explain the difference between Zr and Ag (or Cd). However, with the increase of electronegativity from Zr to Mo, the Ni-H bonds become weak as shown in Fig. 4(b–d) Fig. 4(b–d) , hence, the repulsive interaction to H in the Ni-4d-H supercell become strong. This agrees with the interaction energy in Fig. 2(a) Fig. 2(a) . For Tc, the reduced electronegativity ( Fig. 2b Fig. 2b ) and the smallest atomic size ( Fig. 3 Fig. 3 ) are both helpful to reduce the repulsive interaction to H ( Fig. 2a Fig. 2a ). For Ru, Rh and Pd, there are no obvious differences in the calculated charge density and the charge density difference ( Figs. 5a-5c, 5f-5h Figs. 5a-5c, 5f-5h ), thus the corresponding interaction energies change little Fig. 2a Fig. 2a ). However, for Ag and Cd, the repulsive interactions to H become strong as judged by the interaction energies in Fig. 2(a) Fig. 2(a) , which has a close association with the apparent charge depletion between the Ag (or Cd) and H atoms ( Figs. 5i and 5j Figs. 5i and 5j ) and the relatively big atomic size ( Fig. 3 Fig. 3 ). Especially for Cd, the biggest atomic size and the stable electronic configuration with completely filled 4d and 5s orbitals result in the strongest repulsive interaction to H. The structure of the alloys (e.g. Hastelloy N alloy) is complex and also involves many kinds of defects, which will be addressed in our future works. So far there has been little research on the direct 4d-H interaction in Ni-based alloys by the experimental studies. However, the study of the 4d-H interaction in this paper is helpful for further experimental studies. 4. Conclusion 1 In this work, we use first-principles calculations to investigate the interaction between H and 4d TM elements in Ni-based alloy. The interaction energies between 4d elements and H are calculated. The TM-H interaction has been decomposed into the atomic size effect and electron effects, and both of these effects are analyzed. The effect of atomic size plays a dominant role in the 4d-H interaction, which is mainly different from the 3d-H interaction. We believe that this work is very helpful for the future experimental studies. References LeBlanc D . Molten salt reactors: A new beginning for an old idea . Nucl. Eng. Des . 2010 , 240 : 1644 - 1656 . DOI: 10.1016/j.nucengdes.2009.12.033 http://doi.org/10.1016/j.nucengdes.2009.12.033 Molten salt reactors: A new beginning for an old idea A Technology Roadmap for Generation IV Nuclear Energy Systems, US DOE Nuclear Energy Research Advisory Committee and the Generation IV International Forum, December 2002,GIF-002-00 Available from: http://gifinelgov http://gifinelgov A Technology Roadmap for Generation IV Nuclear Energy Systems, US DOE Nuclear Energy Research Advisory Committee and the Generation IV International Forum, December 2002,GIF-002-00 II Noecker FF, and DuPont JN . Metallurgical investigation into ductility dip cracking in Ni-based alloys: Part II . Weld. J . 2009 , 88 : 62S - 77S . Metallurgical investigation into ductility dip cracking in Ni-based alloys: Part II Louthan MR , Rawl DE , Caskey GR , et al . Hydrogen embrittlement of metals . Mat. Sci. Eng . 1972 , 10 : 357 - 368 . DOI: 10.1016/0025-5416(72)90109-7 http://doi.org/10.1016/0025-5416(72)90109-7 Hydrogen embrittlement of metals Barnoush A and Vehoff H . Recent developments in the study of hydrogen embrittlement: Hydrogen effect on dislocation nucleation . Acta Mater . 2010 , 58 : 5274 - 5285 . DOI: 10.1016/j.actamat.2010.05.057 http://doi.org/10.1016/j.actamat.2010.05.057 Recent developments in the study of hydrogen embrittlement: Hydrogen effect on dislocation nucleation Novakovic N , Novakovic JG , Matovic L , et al . Ab initio calculations of MgH2, MgH2:Ti and MgH2:Co compounds . Int. J. Hydrogen Energ . 2010 , 35 : 598 - 608 . DOI: 10.1016/j.ijhydene.2009.11.003 http://doi.org/10.1016/j.ijhydene.2009.11.003 Ab initio calculations of MgH2, MgH2:Ti and MgH2:Co compounds Qin L and Jiang C . First-principles based modeling of hydrogen permeation through Pd-Cu alloys . Int. J. Hydrogen Energ . 2012 , 37 : 12760 - 12764 . DOI: 10.1016/j.ijhydene.2012.06.029 http://doi.org/10.1016/j.ijhydene.2012.06.029 First-principles based modeling of hydrogen permeation through Pd-Cu alloys Liu W , Qian Y , Zhang D , et al . First-principles calculations of the interaction between hydrogen and 3d alloying atom in nickel . J. Nucl. Mater . 2015 , 465 : 254 - 259 . DOI: 10.1016/j.jnucmat.2015.06.001 http://doi.org/10.1016/j.jnucmat.2015.06.001 First-principles calculations of the interaction between hydrogen and 3d alloying atom in nickel Hu QM , Xu DS , Yang R , et al . First-principles investigation of solute-hydrogen interaction in a alpha-Ti solid solution . Phys. Rev. B 2002 , 66 : 064201 . DOI: 10.1103/PhysRevB.66.064201 http://doi.org/10.1103/PhysRevB.66.064201 First-principles investigation of solute-hydrogen interaction in a alpha-Ti solid solution Li YJ , Kulkova SE , Hu QM , et al . Interaction between hydrogen and the alloying atom in palladium . Phys. Rev. B 2007 , 76 : 064110 . DOI: 10.1103/PhysRevB.76.064110 http://doi.org/10.1103/PhysRevB.76.064110 Interaction between hydrogen and the alloying atom in palladium Hohenberg P and Kohn W . Inhomogeneous Electron Gas . Phys. Rev. B 1964 , 136 : B864 . DOI: 10.1103/PhysRev.136.B864 http://doi.org/10.1103/PhysRev.136.B864 Inhomogeneous Electron Gas Kohn W and Sham LJ . Self-Consistent Equations Including Exchange and Correlation Effects . Phys. Rev . 1965 , 140 : 1133 . DOI: 10.1103/PhysRev.140.A1133 http://doi.org/10.1103/PhysRev.140.A1133 Self-Consistent Equations Including Exchange and Correlation Effects Kresse G and Furthmuller J . Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set . Phys. Rev. B 1996 , 54 : 11169 - 11186 . DOI: 10.1103/PhysRevB.54.11169 http://doi.org/10.1103/PhysRevB.54.11169 Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set Connetable D , Wang Y and Tanguy D . Segregation of hydrogen to defects in nickel using first-principles calculations: The case of self-interstitials and cavities . J. Alloy Compd . 2014 , 614 : 211 - 220 . DOI: 10.1016/j.jallcom.2014.05.094 http://doi.org/10.1016/j.jallcom.2014.05.094 Segregation of hydrogen to defects in nickel using first-principles calculations: The case of self-interstitials and cavities http://en.wikipedia.org/wiki/Electronegativity http://en.wikipedia.org/wiki/Electronegativity . Sun SN , Kioussis N , Lim SP , et al . Impurity effects on atomic bonding in NI3AL . Phys. Rev. B 1995 , 52 : 14421 - 14430 . DOI: 10.1103/PhysRevB.52.14421 http://doi.org/10.1103/PhysRevB.52.14421 Impurity effects on atomic bonding in NI3AL 10.1007/s41365-017-0235-6 This work was supported by the "Strategic Priority Research Program" of the Chinese Academy of Sciences (No. XDA0202020), and the National Natural Science Foundation of China (No. 51601212). 2016-10-21 2016-12-23 2017-01-19 2017-05-06 2017-06-01 2021-04-11 © China Science Publishing & Media Ltd. (Science Press), Shanghai Institute of Applied Physics, the Chinese Academy of Sciences, ChineseNuclear Society and Springer Nature Singapore Pte Ltd. 2017 This is a post-peer-review, pre-copyedit version of an article published in Nuclear Science and Techniques. The final authenticated version is available online at: http://dx.doi.org/10.1007/s41365-017-0235-6 http://dx.doi.org/10.1007/s41365-017-0235-6 . 2017 Nuclear Science and Techniques 06 28