摘要:
采用基于密度泛函理论的第一性原理方法,研究了H在不同单轴应变下α-Fe中的间隙占位,计算了H原子的溶解能、态密度、电荷差分密度和电荷布居.结果表明:不同单轴拉压应变作用下,H原子优先占据四面体间隙(Ts)位,且随着压应变减小、拉应变增加,H原子越易溶于α-Fe.压应变使得Ts位的H获得更多的电子,而拉应变减少了这种电荷转移.应用LST/QST过渡态搜索计算垂直应变方向的扩散.八面体间隙位是邻近Ts位H的扩散过渡态.扩散激活能与应变呈线性关系,且随着压应变的增加,扩散激活能降低,扩散更容易.
Abstract:
As is well known,hydrogen plays an important role in degrading mechanical properties of steel due to hydrogen embrittlement behavior.Thus,much attention should be paid to the interaction between hydrogen atom and Fe matrix especially in theoretical calculation and mechanism study.In this paper,the site occupations of hydrogen atom under different single axis strains in interstitial of α-Fe atoms are studied by the first principles calculation based on the spin-polarized density functional theory.The Kohn-Sham equations are solved under periodic boundary conditions,by using revised Perdew-Burke-Ernzerhof version of the generalized gradient approximation to account for the electron exchange and correlation.The total energy of the steady state crystal,binding energy,solution energy,density of states,charge density difference and charge population are calculated.The analyses of solution energy and density of states indicate that the hydrogen atom preferentially occupies the tetrahedral interstitial of α-Fe atoms under different single axis strains.With increasing tensile strain or reducing compressive strain,hydrogen atom prefers to occupy the site of tetrahedral interstitial.The analyses of charge population and charge density difference reveal that the hydrogen atom collects charges from Fe atoms,leading to electron density redistribution.The tensile strain reduces the charge transfer slightly while the compressive stress promotes the transfer process.The LST/QST (linear synchronous transit/quadratic synchronous transit) transition state search method is used to investigate the diffusion of hydrogen atom between two tetrahedral interstitials along the direction perpendicular to strain.Diffusion of hydrogen atom goes through transition state where the hydrogen atom is coordinated at octahedral interstitial.The minimum energy pathway for hydrogen diffusion under strainless state indicates the diffusion activation energy with a value of 0.58 eV.It is noticeable that the diffusion activation energy and the strain conforms to linear relation and are consistent with the fitting formula,Q =0.508 + 2.6ε.The diffusion activation energy increases with reducing compressive strain or increasing tensile strain.According to the calculation process and analysis results,we introduce the interaction between hydrogen atom and α-Fe at a level of electronic structure systematically and figure out the diffusion of hydrogen influenced by different states of stress.