金属Fe与间隙H原子相互作用的密度泛函研究

采用基于密度泛函理论的第一性原理方法, 研究了不同摩尔比下H在α-Fe和γ-Fe晶格中的间隙占位情况, 计算了稳态晶体的总能量、结合能、溶解热、电子态密度、电荷差分密度和电荷布居, 分析了间隙H原子和Fe晶格之间的相互作用, 讨论了H溶解对α-Fe和γ -Fe晶体电子结构的影响. 结果表明: H溶解引起α-Fe和γ-Fe晶体点阵晶格畸变, 体积膨胀率随着溶氢量的增加而增加. 从能量角度分析发现, H优先占据α-Fe的四面体间隙位, 而在γ -Fe中优先 占据八面体间隙位. 态密度、电荷差分密度以及电荷布居分析发现, 间隙H原子与Fe晶格的相互作用仅由H的1s轨道电子和Fe的4s轨道电子所贡献, 二者作用力相对较弱, 这是造成H在Fe晶格中固溶度较低的主要原因之一. 关键词： 金属Fe 间隙H原子 第一性原理 溶解热     

α-Fe和γ-Fe中氢扩散行为的第一性原理计算

采用基于密度泛函理论的第一性原理方法,研究了α-Fe和γ-Fe中氢扩散行为,建立了两相的H原子扩散系数与温度的关系式,对比了H原子在两相中扩散的难易程度,并结合2205双相不锈钢的氢致开裂试验对计算结果进行了实验验证.研究结果表明:H原子在α-Fe中所需扩散激活能较γ-Fe低,其容易在α-Fe的缺陷处聚集,并最终造成氢致开裂等严重后果;而2205双相不锈钢的HIC实验表明,裂纹在铁素体相中表现为穿晶型开裂,并且在传播过程中由于受到奥氏体相的阻碍,裂纹沿两相界面发展,表现出混合型开裂方式;综上所述,铁素体相更容易产生氢致裂纹,这进一步证实了本文计算结果的可靠性.     

γ-Fe中氢扩散行为的第一性原理研究

不锈钢中氢脆问题一直是研究者们关注热点,然而至今为止,对于钢中合金元素及空位与H原子之间的作用机理依然不清晰.采用第一性原理方法研究了H原子在不同体系下γ-Fe(Fe₈H、Fe₇Cr/MoH、Fe₇H)的扩散激活能、扩散系数以及扩散过渡态,对比了H原子在不同体系γ-Fe中的扩散难易程度,并从电子结构层次解释其相互作用.研究结果表明,空位的存在会改变H原子在γ-Fe中的扩散路径.在含空位的晶胞中H原子稳定存在于空位附近的近四面体间隙中,并沿空位附近的四面体间隙或八面体到四面体间隙路径扩散. γ-Fe中空位对H原子的影响是局域性的,一方面是捕获H原子的陷阱;另一方面降低了H原子的扩散激活能,促进H原子的扩散.而Cr/Mo的掺杂会降低空位对H原子的束缚,也会抑制H原子的扩散.     

稀土La在α-Fe中占位倾向及对晶界影响的第一性原理研究

本文采用重合位置点阵理论构建了α-Fe的Σ3[110](112)对称倾转晶界模型,通过基于密度泛函理论的平面波超软赝势方法研究了稀土La元素在α-Fe中的占位倾向.结果表明,La在α-Fe晶界的杂质形成能最低,因而La原子倾向于占据晶界区;掺杂La前后的α-Fe晶界电子结构计算结果显示,La占位于α-Fe晶界会使体系中的电荷发生重新分配,将提供更多电子用于晶界区成键,使得Fe原子得到更多的电子,这将导致掺杂区原子间结合有离子化趋势,从而使La与晶界区相邻Fe原子之间的相互作用加强,也使晶界原子与晶界两侧Fe原子的键合加强,从能量角度解释了材料宏观力学性能变化的原因;计算同时发现,La加入后,也使晶界上的原子成键区态密度左移,降低了体系的总能量,使晶界结构更为稳定.     

稀土La对bcc-Fe中Cu扩散行为影响的第一性原理研究

采用基于密度泛函理论的第一性原理,研究了稀土La对bcc-Fe中Cu析出行为的影响.计算了α-Fe中La原子和Cu原子与空位之间,以及La原子和Cu原子之间的点缺陷结合能;在此基础上,讨论了α-Fe中La对Cu扩散激活能的关系.结果表明:La原子与空位之间有较强的相互吸引作用,且对近邻Cu原子也有一定的束缚.此外,La的加入使Cu原子近邻的空位形成能显著升高,这表明La,Cu偏聚区形成空位较为困难.与此同时,由于La原子对近邻空位和Cu原子的吸引作用,使Cu原子向近邻空位跳跃的迁移能有所升高.迁移能与空位形成能变化的计算结果显示,La原子的加入能够使α-Fe中Cu的扩散激活能显著升高,从而延缓了铁素体区富铜相的偏聚和析出.     

合金元素对α-Fe氢扩散行为影响的第一性原理研究

本文采用基于密度泛函理论的第一性原理方法,研究了H原子在α-Fe和α-Fe-Me(Me=Cr、Mo、Ni、Cu、Mn)晶胞中的扩散行为,建立了α-Fe和α-Fe-Me中氢扩散系数与温度间的关系表达式,探讨了合金元素Me对α-Fe中H扩散行为的影响作用.研究结果表明:合金原子Me与H原子间存在一定的键合作用,并且这种键合作用主要由Me价电子的d、s轨道与H的1s轨道所贡献;Cr、Ni、Cu、Mn元素将降低H原子扩散激活能,而Mo元素会提高氢扩散激活能,显著降低氢扩散系数,这说明添加适量Mo元素将有利于防止铁素体氢致开裂,这将对抗H₂S应力腐蚀新材料的设计提供一定理论依据.     

单轴应变对H在α-Fe中占位及扩散的影响

采用基于密度泛函理论的第一性原理方法,研究了H在不同单轴应变下α-Fe中的间隙占位,计算了H原子的溶解能、态密度、电荷差分密度和电荷布居.结果表明:不同单轴拉压应变作用下,H原子优先占据四面体间隙(Ts)位,且随着压应变减小、拉应变增加,H原子越易溶于α-Fe.压应变使得Ts位的H获得更多的电子,而拉应变减少了这种电荷转移.应用LST/QST过渡态搜索计算垂直应变方向的扩散.八面体间隙位是邻近Ts位H的扩散过渡态.扩散激活能与应变呈线性关系,且随着压应变的增加,扩散激活能降低,扩散更容易.     

NaZn13型Fe基化合物的结构和热力学性质研究

利用Chen-Mbius晶格反演获得的原子间相互作用势,对NaZn13型Fe基金属间化合物进行原子级模拟研究.计算结果表明,Si原子和Co原子均优先占据96i晶位,Si原子和Co原子替代Fe原子后晶体平均结合能降低.随着Co含量的增加,LaFe13-x-yCoySix和NdFe13-x-yCoySix的晶格参数逐渐降低.声子态密度中,稀土原子主要激发低频模,Si原子主要激发高频模.LaFe11.5-yCoySi1.5化合物的德拜温度随Co含量的增加而增高.     

H,O,N,He原子掺杂Al原子的第一性原理研究

采用密度泛函理论框架下的第一性原理平面波赝势方法,对Al中分别加入H,O,N和He原子后的晶体状态进行了研究.通过晶体结构和形成能分析比较了杂质原子占据不同位置的难易程度及对晶体稳定性的影响,并从态密度、电荷密度和电荷布居的角度,分析了其电子结构.结果表明:H、O和N原子占据金属Al的四面体间隙最稳定,而He原子主要占据金属Al的八面体间隙. O和N原子与Al原子具有强烈的共价作用,H原子与Al原子存在共价作用但相对较弱,而He原子与Al原子的相互作用以范德华力为主.进一步揭示了四种原子在金属Al中不同行为的电子机制.     

过渡金属氮化物中单原子固溶和力学性能的第一性原理研究

基于密度泛函理论(DFT)的第一性原理,计算Si原子在Ti族和V族氮化物中以及B、C和Ge原子在TiN晶体中固溶的稳定结构,讨论置换型和间隙型固溶的低能量稳定结构与晶体间距的关系,研究金属氮化物和固溶原子固溶结构的力学性能.结果表明:Si原子在TiN、ZrN、HfN和TaN晶体中固溶以及Ge原子在TiN晶体中固溶情况为,单原子不进入对应过渡金属氮化物晶体中形成间隙固溶或置换固溶,随着晶体间距离变化单原子可以在晶体之间形成间隙固溶或置换固溶;Si原子在NbN以及B原子在TiN晶体中可以实现间隙固溶,而不能形成置换固溶;Si原子在VN和C原子在TiN晶体中固溶结构形式均为置换固溶.单原子固溶形成低能量置换型固溶体和间隙型固溶体的弹性常数、体模量和剪切模量均低于原过渡金属氮化物的对应值.     

Vacancy trapping mechanism for multiple hydrogen and helium in beryllium: a first-principles study

The microscopic mechanism for H and He trapping by vacancy defects and bubble formation in a Be host lattice is investigated using first-principles calculations. A single He atom prefers to occupy a vacancy centre while H does not. He can segregate towards the vacancy from the interstitial site much more easily than H. Both H and He exhibit lower diffusion barriers from a remote interstitial to a vacancy with regard to their diffusion barriers inside a perfect Be solid. Up to five H or 12 He atoms can be accommodated into the monovacancy space, and the Be-He interaction is much weaker than Be-H. The physical origin for aggregation of multiple H or He atoms in a vacancy is further discussed. The strong tendency of H and He trapping at vacancies provides an explanation for why H and He bubbles were experimentally observed at vacancy defects in materials. We therefore argue that vacancies provide a primary nucleation site for bubbles of H and He gases inside Be materials.     

Effects of Cr on H and He trapping and vacancy complexes in V in a fusion environment: a first-principles study

First-principles density functional theory was used to investigate the interaction between hydrogen (H) and helium (He) in V–Cr alloy, which is a potential structural material for use in fusion reactors. When vacancies are present in the V–Cr alloy, a single He atom prefers to occupy the octahedral site near the Cr atom rather than vacancy centre, which differs from the cases of iron and tungsten. Because of the decrease of the electron density around the He atom, there was a strong interaction between He and H. In the vicinity of He-vac complexes, H atoms tend to stay in the tetrahedral site rather than occupy the octahedral-interstitial site. A single He-vac complex can trap as many as six H atoms, which is more than can be trapped by an empty vacancy in the V–Cr alloy because of the electronic density redistribution of vacancy vicinity. This strong attraction explains the enhanced retention of H and He observed near the surface of V and V-based alloys under both sequential and simultaneous bombardments. The results provide useful insight into the application of the V-based alloys as candidate structural materials in fusion Tokamaks.     

The multiple hydrogen location near a α-Fe vacancy

The interaction between 10 hydrogen atoms and a α-Fe structure having a vacancy (V) has been studied using a cluster model and a semi-empirical theoretical method. The energy of the system was calculated by the atom superposition and electron delocalization molecular orbital method. The electronic structure was studied using the concept of density of states and crystal orbital overlap population curves.For the study of a sequential absorption, the hydrogen atoms were positioned in their energy minima configurations, near to the tetrahedral sites neighboring the vacancy, except the last H atom that was located far from the vacancy. The energy difference for H agglomeration was also computed. The vacancy-H_n complexes become less stable than VH species for more than three hydrogen's.The changes in the electronic structure of Fe atoms near the vacancy were also analyzed. The interactions mainly involve Fe 3d and 4s atomic orbitals. The contribution of Fe p orbitals is much less important. The Fe-Fe bond weakened as new Fe-H and H-H pairs were formed. The effect of H atoms is limited to its first Fe neighbors. The detrimental effect of H atoms on the Fe-Fe bonds can be related to the mechanism for embrittlement in α-Fe.     

Ab initio thermodynamic characteristics of the formation of oxygen vacancies,and boron,carbon, and nitrogen impurity centers in anatase

Using the ab initio projector augmented wave (PAW) method, calculations are performed for the electronic energy-band structure of titanium dioxide having the structure of anatase doped with boron, nitrogen, and carbon. Thermodynamic characteristics are determined for the formation of impurity centers, such as the preference energy for the interstitial position, the energy of impurity oxidation, and the energy of oxygen vacancy formation. It is shown that under the conditions of thermodynamic equilibrium the interstitial position of boron atoms is stable, whereas carbon atoms, depending on the oxygen pressure, can occupy both interstitial positions and substitutional positions of oxygen atoms, and nitrogen atoms replace oxygen atoms. It is shown that the presence of oxygen vacancies promotes the thermodynamic stability of carbon and nitrogen atoms. The obtained densities of electronic states correspond to ESR spectroscopy data, which indicates the presence of spin-polarized electrons in the states of the oxygen vacancy.     

Energetics,electronic structure,and positron annihilation studies of carbon-vacancy complexes in iron

Using the density functional theory and the embedded cluster model, we have calculated the energetics of carbon-vacancy complexes in Fe as a function of distance separating the carbon atom(s) and the vacancy along different crystallographic directions. Carbon is found to prefer off-center sites from the vacancy and maintains a constant distance of approximately 3.35a ₀ from the nearest Fe atom(s) independent of its direction from the vacancy center. This distance agrees closely with that in stoichiometric Fe₃C. Our results suggest that the equilibrium site of carbon is one where it is coordinated to three Fe atoms in its nearest neighbor shell. Further decoration of vacancies by more than one carbon atom has been found to be energetically favorable. The binding of carbon to iron atoms is not caused by charge transfer between the atoms, but rather has a magnetic origin. Positron lifetimes at vacancies and vacancy-carbon complexes have also been calculated for various configurations. Our results are in general agreement with experiment.     

The unambiguous identification of vacancy migration stage in resistivity recovery of irradiated metals based on trapping of vacancies at probe impurity atoms

The method of resistivity recovery is a powerful tool for studying point defects in irradiated metals. However, the method is nonspecific with respect to a type (vacancy or interstitial) of investigated defects. To overcome this shortcoming, we made use of opposite signs of excess electric charges of the vacancies and self-interstitial atoms in a lattice. Resistivity loss takes place on trapping of vacancies at impurity atoms if excess charges of the impurity atom and defect are opposite in sign. A way of selecting the impurity atoms with the excess charge opposite in sign to that of vacancies is proposed. The specific evolution of resistivity recovery spectrum induced by vacancy trapping at the selected impurity atoms (probe traps) allows one to unambiguously identify the stage of free long-range vacancy migration.     

Na在B掺杂的空位石墨烯上吸附性能的第一性原理研究

采用基于密度泛函理论的第一性原理方法,研究了本征石墨烯和B掺杂的空位石墨烯吸附Na原子的电荷密度、吸附能、态密度、储存量以及电极电压.结果表明,两种石墨烯中,Na原子的最佳吸附位置都是H位.B掺杂的空位石墨烯对Na原子的吸附能是-2.08 eV,比本征石墨烯对Na原子的吸附能(-0.71eV)低很多.B掺杂的空位石墨烯中Na原子与B原子发生轨道杂化,本征石墨烯中没有杂化现象.B掺杂的空位石墨烯能够吸附12个Na原子,较本征石墨烯多.因此,B掺杂的空位石墨烯更适合储钠.