合金元素对 α-Fe中H原子溶解及扩散影响的第一性原理研究
1.
2.
3.
First-Principles Study on the Influence of Alloying Elements on the Dissolution and Diffusion of Hydrogen Atoms in α-Fe
1.
2.
3.
通讯作者: 马 会,hma@imr.ac.cn,主要从事材料计算及相关研究
责任编辑: 李海兰
收稿日期: 2024-07-17 修回日期: 2024-10-23
| 基金资助: |
|
Corresponding authors: MA Hui, associate professor, Tel:
Received: 2024-07-17 Revised: 2024-10-23
| Fund supported: |
|
作者简介 About authors
宝尔根,男,1993年生,博士生
氢脆会导致暴露在含H环境中的金属变脆并开裂,对能源、制造、运输和航空航天等行业构成严重风险,理解H与钢的相互作用在控制氢脆问题中起着至关重要的作用。本工作利用基于密度泛函理论的第一性原理方法,计算了不同合金元素对H原子溶入完整的及有空位缺陷的α-Fe基固溶体的影响,并根据晶体轨道Hamilton布居分析了H原子与合金元素、Fe原子的相互作用,并用爬坡图像弹性带(CI-NEB)方法计算了合金元素对H原子在α-Fe基固溶体中扩散的影响。结果表明,H原子均倾向溶于合金原子的第二、第三近邻四面体间隙,且溶解焓取决于H与最近邻的合金原子或Fe原子之间的键强;有空位缺陷时,H原子在合金化后Sc、V、Cr、Mn、Co、Ni、Cu的第一近邻空位处的溶解焓比在完整固溶体的四面体间隙位的溶解焓低,而固溶Al、Si、Ti、Zr、Nb、Eu、W的效果则相反;Ti、Cu、Zr、Nb及稀土元素的固溶会增加H原子从合金原子的次近邻至第三近邻四面体间隙的扩散能垒,同时降低其反方向的能垒,而固溶Si、V、Cr、Mn、Co、Ni、Mo、W的效果则相反。
关键词:
Hydrogen embrittlement causes metals exposed to H-containing environments to become brittle and crack, posing serious risks to industries such as energy, manufacturing, transportation, and aerospace. Understanding the interaction between hydrogen and steel is crucial for addressing hydrogen embrittlement challenges. This study uses first-principles methods based on the density functional theory to examine how different alloying elements influence the dissolution of H atoms in α-Fe-based solid solutions, both with and without vacancy defects. The interactions between H atoms and alloying elements, as well as Fe atoms are analyzed using crystal orbital Hamilton populations. Additionally, the climbing image-nudged elastic band (CI-NEB) method was used to calculate the influence of the alloying elements on the diffusion of H atoms in α-Fe-based solid solutions. The results indicate that H atoms preferentially dissolve in the second and third nearest tetrahedral interstitial sites of alloy atoms. The solution enthalpy of H atoms is determined by the bond strength between the H atom and the nearest alloy atom as well as Fe atoms. For systems with vacancy defects after alloying, the solution enthalpies of H atoms at the first nearest vacancy for elements, such as Sc, V, Cr, Mn, Co, Ni, and Cu are lower than those at the tetrahedral interstitial sites of the complete α-Fe-based solid solutions. However, the opposite effect is observed for systems containing Al, Si, Ti, Zr, Nb, Eu, and W. Furthermore, alloying with Ti, Cu, Zr, Nb, and rare-earth elements increases the diffusion energy barriers for H atoms moving from the second nearest tetrahedral interstitial site to the third nearest site, while reducing diffusion energy barriers in the opposite direction. Conversely, the opposite effects are observed in α-Fe-based solid solutions containing Si, V, Cr, Mn, Co, Ni, Mo, and W.
Keywords:
本文引用格式
宝尔根, 王建韬, 徐文静, 马会, 陈星秋.
BAO Ergen, WANG Jiantao, XU Wenjing, MA Hui, CHEN Xing-qiu.
H2是现在以及将来许多工业所用的重要材料,是一种多用途、可持续的能源。这种可再生能源越来越多地用于各个行业,特别是在电力生产、运输和国内能源供应方面[1,2]。利用现有的天然气管道输送H2是一种节约能源和减少排放的方法。然而,由于管道钢的材料特性,即使是最低浓度的H也会导致金属和合金的脆化[3,4],这大大降低了管道的力学性能,严重影响了金属的断裂韧性[5,6],从而降低管道的耐久性和使用寿命,对管道的结构完整性和安全性构成威胁[7,8]。另一方面,在轴承服役过程中,有机润滑油的裂解产生许多H元素,这些H原子吸附在轴承表面,随后扩散到轴承钢中,导致氢脆,对轴承的使用性能产生负面影响[9,10]。输氢管道和轴承钢并不完全由单一金属Fe组成,在钢材料中掺杂Cr、Mn、Ni、Mo等合金元素可以改善钢的化学性能和力学性能[11~13],改性的合金成分可以提高合金钢的抗氢脆性能。因此,有必要进一步深入理解和揭示合金元素对H原子溶入钢内及其在钢中扩散的影响。
研究H原子的吸附和扩散过程是开展氢致开裂研究的基础。Yu等[14]利用密度泛函理论(density functional theory,DFT)和从头算原子热力学证明了在α-Fe表面H2分子在动力学和热力学上都倾向于解离,H在表面的稳定覆盖高度率依赖于温度和H2分压。Dadfarnia等[15]提出了一种修正的H输运模型,认为位错有助于提高裂纹起裂前的局部H浓度。Yanachkov等[16]通过对纯Fe进行应力松弛实验和自洽动力学Monte Carlo模拟,发现充氢影响位错的迁移,导致冷轧和退火样品中更容易形成位错胞结构。Zhang等[17]采用第一性原理计算方法,从吸附能和态密度两方面分析了合金元素对α-Fe (110)表面吸附H的影响。Omura等[18]通过氢渗透实验研究了充氢条件下合金元素对铁素体的微观结构、氢扩散率和拉伸性能的影响,发现氢扩散系数随合金浓度的升高而降低。然而,通过实验并不能详细地了解氢扩散的机制。Liu等[19]用理论和实验检测碳化物的H捕获行为,证明添加Mo元素增强了氢捕获能力,使H能够进入碳化物内的C空位陷阱,但这是从碳化物的角度所进行的研究而非钢中的基体α-Fe。
1 计算方法与模型
本工作使用基于第一性原理的计算机程序——Vienna从头算模拟包(Vienna ab initio simulation package,VASP)[22~24]对材料进行原子尺度模拟。使用广义梯度近似(GGA)-Perdew-Burke-Ernzerhof (PBE)泛函[25]及以平面波为基函数的投影缀加波(PAW)方法[26]。其中,截断能取为520 eV,能量的收敛标准为4 × 10-6 eV,力的收敛标准为0.15 eV/nm。使用一阶Methfessel-Paxton方法[27]确定展宽为0.1 eV的电子占据率。在结构优化的计算中,k空间由Monkhorst-Pack方法[28]分成3 × 3 × 3的网格。对于电子态密度(DOS)的计算,采用4 × 4 × 4的k空间网格,并用Blöchl修正的四面体法[29]来计算k空间中的积分。用VASP做完自洽计算后,用LOBSTER程序[30]读取电子波函数分析晶体轨道Hamilton布居(COHP)[31]。对于扩散路径中过渡态的搜索,采用VTST编译后的VASP中的爬坡图像弹性带(climbing image-nudged elastic band,CI-NEB)方法[32]。
合金原子M固溶于α-Fe中的固溶焓(ΔHM )为:
式中,
间隙H原子的溶解焓(ΔHH)为:
式中,
空位生成焓(ΔHvac)表达式为:
式中,
钢中常用的合金元素有B、C、N、Al、Si、Ti、V、Cr、Mn、Co、Ni、Cu、Zr、Nb、Mo、W、稀土元素等。当合金元素的含量较少时,合金元素会以溶质原子的形式固溶在α-Fe、γ-Fe和Fe3C中,形成以α-Fe、γ-Fe和Fe3C为基的固溶体[33]。其中,B、C、N原子(由于原子尺寸小)间隙固溶于α-Fe、γ-Fe和Fe3C中,而Al、Si、Ti、V、Cr、Mn、Co、Ni、Cu、Zr、Nb、Mo、W和稀土原子则是置换固溶于α-Fe、γ-Fe和Fe3C中。本工作考虑了常见合金元素Al、Si、Ti、V、Cr、Mn、Co、Ni、Cu、Zr、Nb、Mo、W及稀土元素置换固溶于α-Fe相中的情况。为兼顾计算的时间成本及合金元素含量较少这个要求,选取了对α-Fe单胞以3 × 3 × 3扩胞得到的超胞,共54个Fe原子,并置换其中一个Fe原子为合金原子(图1a)。
图1
图1
α-Fe基固溶体的晶体结构
Fig.1
Crystal structures of α-Fe-based solid solutions
(a) the structure with a supercell and a substituted alloy atom (Purple and green balls represent Fe atoms and an alloy atom, respectively)
(b) the H atom dissolved in the tetrahedral interstitial site (Yellow ball represents the H atom)
(c) the H atom dissolved in the first, second, and third nearest neighbor tetrahedral interstitial sites around the alloy atom (View parallel to c-axis; the distances from the alloy atom to the H atom in the first, second, and third nearest neighbor positions are 0.167, 0.259, and 0.328 nm, respectively)
(d) the H atom dissolved in the first, second, and third nearest neighbor tetrahedral interstitial sites around the alloy atom (View parallel to b-axis)
(e) vacancy defects near the alloy atom (Two yellow balls represent the first and second nearest neighbor vacancies around the alloy atom, respectively)
(f) diffusion paths of the H atom (From the second nearest tetrahedral interstitial site to the third nearest tetrahedral interstitial site, and from the third nearest tetrahedral interstitial site to the fourth nearest tetrahedral interstitial site)
2 计算结果与分析
2.1 合金元素置换固溶于 α-Fe中的难易程度
首先计算了合金元素Al、Si、Ti、V、Cr、Mn、Co、Ni、Cu、Zr、Nb、Mo、W和稀土元素固溶于α-Fe中的固溶焓,结果如图2所示。可见,Al、Si及过渡元素中的Ti、V、Cr、Co、Nb、Mo、W的固溶焓为负值,能自发地溶入α-Fe中;过渡元素中的Sc、Mn、Ni、Cu、Y、Zr及所有稀土元素的固溶焓为正值,不能自发溶入α-Fe中,其中Mn、Ni溶入α-Fe相对容易。整体而言,稀土元素固溶焓随着原子半径增大而降低,但其中Ce元素的固溶焓低于其他La系稀土元素。这与文献[33]中V、Cr能无限固溶,Al、Si、Mn、Co、Ni、Mo、W能以较宽的固溶度有限固溶,Ti、Cu以较窄的固溶度有限固溶,Zr、Nb以很小的固溶度固溶的结论大体一致。各组元的固溶度取决于点阵类型、原子尺寸和电子结构,与文献[33]有一定差别,其原因是第一性原理计算的是在0 K下的情况,而实际可能是在高温下固溶的合金元素。
图2
图2
α-Fe基固溶体中合金元素的固溶焓
Fig.2
Solid solution enthalpies of alloying elements in α-Fe-based solid solutions
2.2 合金元素对H原子溶入 α-Fe基固溶体的影响
文献[34,35]表明,相较于α-Fe中的八面体间隙位置,H原子倾向溶于四面体间隙位置,如图1b所示。计算的H原子在α-Fe四面体间隙位置的溶解焓为0.42 eV,与成应晋等[35]的计算值(0.39 eV)基本吻合。由图3a可以看出,溶入的H原子与周围的Fe原子成键。图3b1和b2分别为溶H后α-Fe的DOS图及投影晶体轨道Hamilton布居(pCOHP)图,-pCOHP图标示了键合(-pCOHP > 0)和反键合(-pCOHP < 0)对能带结构能量的贡献。可以看出,在(-8.8 eV + EF ~ -8.4 eV + EF,EF为Fermi能级)能量区间,Fe的4s、3d轨道电子与H的1s轨道电子成键。
图3
图3
溶H后α-Fe的电荷密度、电子态密度及投影晶体轨道Hamilton布居(pCOHP)图
Fig.3
Charge density, electronic density of states (DOS), and projected crystal orbital Hamilton populations (pCOHP) of α-Fe containing H
(a) charge density on the (100) plane (The iso-surface value is set to 0.04 electrons/a
(b1, b2) total electronic density of states (TDOS) and projected density of states of α-Fe-based solid solution containing H (b1) and negative pCOHP (-pCOHP) between H and the nearest, next-nearest Fe atoms (b2) (E—electron energy, EF—Fermi energy; 1NN and 2NN represent the nearest neighbors and next-nearest neighbors, respectively; “↑, ↓” represent spin-up and spin-down, respectively)
在固溶合金元素后,分别计算了H原子溶入合金原子的第一近邻、第二近邻及第三近邻四面体间隙位置(图1c和d)时的溶解焓。计算结果表明,H原子在大多数合金原子的第一近邻四面体间隙处不稳定,除了Cr、Mn、Co、Ni、Cu这几个原子序数与Fe接近的合金元素外,结构优化后H原子都会从合金原子的第一近邻四面体间隙移动到第二近邻四面体间隙位置。如图4所示,本工作所考虑的合金元素置换固溶在α-Fe基体时,H原子都更倾向溶于合金原子的第二、第三近邻四面体间隙,尤其是稀土元素固溶时,H原子倾向溶于合金原子的第二近邻四面体间隙位置。并且,这些合金元素的固溶都会降低H溶于第二、第三近邻四面体间隙的难度,尤其是Zr及稀土元素。
图4
图4
H原子在α-Fe基固溶体中不同四面体间隙位置的溶解焓
Fig.4
Solution enthalpies of H atoms in different tetrahedral interstitial sites within α-Fe-based solid solution (Gray dots, green triangles, and pink stars represent the cases where H atoms dissolve in the first, second, and third nearest tetrahedral interstitial sites of the alloy atoms, respectively; the dark green dashed line corresponds to the situation where H atoms dissolve in the interstitial sites of pure α-Fe, with an solid solution enthalpy of 0.42 eV)
这是由于,H原子溶于合金原子的第二近邻四面体间隙时,Zr及稀土元素与H原子之间形成了较强的化学键。使用-IpCOHP表示-pCOHP对能量在(-∞, EF)上的积分,它代表了特定键的强度。由补充材料(图S1)可以看出,合金原子的固溶后,H与最近邻Fe (键长约0.167 nm)的键强没有较大的变化,而H与合金原子(键长约0.259 nm)的键强变化较大,且IpCOHP的变化趋势与固溶不同合金后H的溶解焓的变化一致(固溶Yb、Eu除外),见图5中的绿色折线。这说明H溶于合金原子的第二近邻四面体间隙的难易取决于不同合金原子与H原子的成键强弱,合金原子与H原子成键越强(-IpCOHP越大),H原子就越容易溶入(ΔHH越低)。
图5
图5
固溶不同的合金元素对于H与合金原子、H与最近邻Fe之间的IpCOHP (pCOHP对能量在(-∞, EF)上的积分)的影响
Fig.5
Influence of alloying different elements on the integrated pCOHP (IpCOHP) between H and alloy atoms, as well as between H and the nearest Fe atoms (The green line represents the IpCOHP between H and the alloy atoms when H is dissolved in the second nearest tetrahedral interstitial site of the alloy atoms; the blue line represents the IpCOHP between H and the nearest Fe atoms when H is dissolved in the third nearest tetrahedral interstitial site of the alloy atoms; the purple dashed line corresponds to the IpCOHP between H and the next nearest Fe atoms when H is interstitially dissolved in pure α-Fe, which is -0.061 eV; the yellow dashed line corresponds to the IpCOHP between H and the nearest Fe atoms when H is interstitially dissolved in pure α-Fe, which is -0.478 eV)
2.3 合金元素对H原子溶于有空位缺陷的 α-Fe基固溶体的影响
图6
图6
在α-Fe基固溶体中合金原子的第一、第二近邻原子处产生空位的生成焓
Fig.6
Vacancy formation enthalpies at the first and second nearest sites of alloy atoms in α-Fe-based solid solutions (The black line and green line represent the first and second nearest vacancies of different alloy atoms, respectively; the dark green dashed line corresponds to the case of a single vacancy in pure Fe, with a formation enthalpy of 2.19 eV)
然后计算了H原子溶于合金原子的第一近邻位空位的情况,结果见图7所示。由于空位所具有的空间更大,所以H原子在α-Fe空位处的溶解焓(0.36 eV)比在完整α-Fe的四面体间隙位的溶解焓(0.42 eV)低。在Sc、V、Cr、Mn、Co、Ni、Cu合金化后,H原子在合金原子第一近邻空位处的溶解焓比在纯α-Fe空位处的溶解焓低,而且稀土原子的大尺寸会挤压第一近邻空位的空间,导致H原子进入稀土原子旁空位的溶解焓比其他合金元素略高;而H原子在Al、Si、Ti、Zr、Nb、Eu、W合金原子第一近邻空位处的溶解焓比在纯α-Fe空位处的溶解焓高。
图7
图7
H原子溶于合金原子的第一近邻空位处的溶解焓
Fig.7
Solid solution enthalpies of H atoms at the first nearest vacancy of alloy atoms (Light blue dashed line corresponds to the solid solution enthalpy of H in the vacancy of pure α-Fe, which is 0.36 eV; dark green dashed line corresponds to the solid solution enthalpy of H in the tetrahedral interstitial site of complete pure α-Fe, which is 0.42 eV)
2.4 合金元素对H原子在 α-Fe基固溶体中扩散的影响
图8
图8
H原子从各固溶合金原子的次近邻扩散至第三近邻四面体间隙位置的过渡态搜索
Fig.8
Transition state search for the diffusion of the H atom from the second nearest tetrahedral interstitial site to the third nearest tetrahedral interstitial site of different alloy atoms (Taking the energy of H atoms in the tetrahedral interstitial site of pure α-Fe as the reference)
(a) non-rare earth elements in solid solution
(b) rare earth elements in solid solution
由图8可以看出,Ti、Cu、Zr、Nb及所有稀土元素的固溶使得H原子从合金原子的次近邻四面体间隙至第三近邻四面体间隙的扩散能垒升高(0.01~0.13 eV),同时使反方向的扩散能垒降低(0.01~0.04 eV)。这意味着,Ti、Cu、Zr、Nb及稀土元素的次近邻四面体间隙位容易捕获H原子。与之相反,Si、V、Cr、Mn、Co、Ni、Mo、W的固溶使得H从次近邻至第三近邻四面体间隙的扩散能垒降低(0~0.03 eV),同时反方向的扩散能垒升高(0~0.02 eV)。这意味着,Si、V、Cr、Mn、Co、Ni、Mo、W的固溶会阻碍H被捕获。
综合考虑合金元素对于H原子在α-Fe基固溶体间隙处及空位处的溶解焓、空位的形成焓以及H原子从合金原子的次近邻至第三近邻四面体间隙位(及其反方向)扩散的影响,Zr和稀土元素合金化后的α-Fe基固溶体使得H原子更容易溶于合金原子附近,且不易扩散出去;Al、Si、Cr、Ni、Mo、W合金化后的α-Fe基固溶体可能能够抵抗氢脆,这与较高的Ni含量有利于抵抗氢脆[36]一致。
3 结论
(1) 在α-Fe基固溶体中,H原子都倾向溶于合金原子的第二、第三近邻四面体间隙。H在合金原子的次近邻四面体间隙的溶解焓取决于H与合金原子之间的成键强弱,而H在合金原子的第三近邻四面体间隙的溶解焓取决于H与最近邻的Fe的成键强弱。
(2) 在α-Fe基固溶体中,在合金元素的第一近邻位处更加容易产生空位,并且合金元素(除了V和Co)的固溶都会降低空位的生成焓。H原子在纯α-Fe空位处的溶解焓比在完整α-Fe的四面体间隙位的低,且H原子在合金化后Sc、V、Cr、Mn、Co、Ni、Cu的第一近邻空位处的溶解焓比在完整固溶体的四面体间隙位的低,而H原子在合金化后Al、Si、Ti、Zr、Nb、Eu、W的第一近邻空位处的溶解焓比在完整固溶体的四面体间隙位的高。
(3) 在α-Fe中固溶Ti、Cu、Zr、Nb或稀土元素,会使H原子从合金原子的次近邻至第三近邻四面体间隙的扩散能垒升高,而降低反方向的扩散能垒;而固溶Si、V、Cr、Mn、Co、Ni、Mo、W会使H从合金原子的次近邻至第三近邻四面体间隙的扩散能垒降低,同时增加反方向的扩散能垒。
(4) Zr和稀土元素合金化后的α-Fe基固溶体使得H原子更容易溶于合金原子附近,且不易扩散;Al、Si、Cr、Ni、Mo、W合金化后的α-Fe基固溶体可能能够抵抗氢脆。
参考文献
Estimation of the influence of compressed hydrogen on the mechanical properties of pipeline steels
[J].Consideration of the possibility of transporting compressed hydrogen through existing gas pipelines leads to the need to study the regularities of the effect of hydrogen on the mechanical properties of steels in relation to the conditions of their operation in pipelines (operating pressure range, stress state of the pipe metal, etc.). This article provides an overview of the types of influence of hydrogen on the mechanical properties of steels, including those used for the manufacture of pipelines. The effect of elastic and plastic deformations on the intensity of hydrogen saturation of steels and changes in their strength and plastic deformations is analyzed. An assessment of the potential losses of transported hydrogen through the pipeline wall as a result of diffusion has been made. The main issues that need to be solved for the development of a scientifically grounded conclusion on the possibility of using existing gas pipelines for the transportation of compressed hydrogen are outlined.
Ensuring natural gas infrastructure is suitable for hydrogen service
[J].
Hydrogen effects on X80 pipeline steel in high-pressure natural gas/hydrogen mixtures
[J].
Damage associated with interactions between microstructural characteristics and hydrogen/methane gas mixtures of pipeline steels
[J].
Hydrogen transport and metal embrittlement risk in storage and industrial applications
[J].
Effect of hydrogen gas impurities on the hydrogen dissociation on iron surface
[J].
Hydrogen embrittlement of medium Mn steels
[J].Recent research efforts to develop advanced–/ultrahigh–strength medium-Mn steels have led to the development of a variety of alloying concepts, thermo-mechanical processing routes, and microstructural variants for these steel grades. However, certain grades of advanced–/ultrahigh–strength steels (A/UHSS) are known to be highly susceptible to hydrogen embrittlement, due to their high strength levels. Hydrogen embrittlement characteristics of medium–Mn steels are less understood compared to other classes of A/UHSS, such as high Mn twinning–induced plasticity steel, because of the relatively short history of the development of this steel class and the complex nature of multiphase, fine-grained microstructures that are present in medium–Mn steels. The motivation of this paper is to review the current understanding of the hydrogen embrittlement characteristics of medium or intermediate Mn (4 to 15 wt pct) multiphase steels and to address various alloying and processing strategies that are available to enhance the hydrogen-resistance of these steel grades.
Hydrogen stress cracking behaviour in dissimilar welded joints of duplex stainless steel and carbon steel
[J].As the need for duplex stainless steel (DSS) increases, it is necessary to evaluate hydrogen stress cracking (HSC) in dissimilar welded joints (WJs) of DSS and carbon steel. This study aims to investigate the effect of the weld microstructure on the HSC behaviour of dissimilar gas-tungsten arc welds of DSS and carbon steel. In situ slow-strain rate testing (SSRT) with hydrogen charging was conducted for transverse WJs, which fractured in the softened heat-affected zone of the carbon steel under hydrogen-free conditions. However, HSC occurred at the martensite band and the interface of the austenite and martensite bands in the type-II boundary. The band acted as an HSC initiation site because of the presence of a large amount of trapped hydrogen and a high strain concentration during the SSRT with hydrogen charging. Even though some weld microstructures such as the austenite and martensite bands in type-II boundaries were harmless under normal hydrogen-free conditions, they had a negative effect in a hydrogen atmosphere, resulting in the premature rupture of the weld. Eventually, a premature fracture occurred during the in situ SSRT in the type-II boundary because of the hydrogen-enhanced strain-induced void (HESIV) and hydrogen-enhanced localised plasticity (HELP) mechanisms.
Study on mechanism of hydrogen generation from lubricants
[J].
Hydrogen permeation through steel during cathodic polarization of lubricating oils in a modified Devanathan-Stachurski cell
[J].In lubricated tribo-contacts, hydrogen ingress in steel is possible due to chemical reactions of lubricant components like base oils or additives, and/or contamination upon service particularly water, and/or corrosion processes, and/or electrostatic fields or current flow. Absorbed by the metal, atomic hydrogen may cause serious deleterious effects on the physical-chemical and mechanical properties, reducing the material's ability to withstand the design loads. The present research work is focused on analyzing the influence of electric field on lubricating oils in contact with steel surface. In order to evaluate the possibility of atomic hydrogen generation and permeation into the steel under cathodic polarization of lubricating oils the electrochemical permeation technique developed by Devanathan and Stachurski is used. The input cell of a Devanathan-Stachurski set up is appropriately modified by realizing a very close distance between the working electrode (steel membrane) and a Pt counter electrode with the oil between. This significantly increases the sensibility of the set up and allows the application of larger voltage and higher temperature to enable hydrogen generation from lubricating oils. The complex effects of cathodic polarization, temperature, additives and presence of water in model lubricating oils on atomic hydrogen permeation into steel is discussed.© 2022. The Author(s).
Alloying element additions to Ni3Al: Site preferences and effects on elastic properties from first-principles calculations
[J].
Characterization of hot deformation behavior of Ti-3Al-5Mo-4.5V alloy with a martensitic starting microstructure
[J].The hot deformation behavior of Ti–3Al–5Mo–4.5V alloy with an [Formula: see text] martensitic microstructure was studied in the temperature range of 700–1000[Formula: see text]C and strain rates range of 0.001–[Formula: see text] up to a height reduction of 50%. The results show that an ultrafine equiaxed microstructure with an average grain size of [Formula: see text] were successfully produced through thermomechanical processing of a martensitic starting microstructure. A processing map was successfully constructed and provides appropriate processing parameters for hot deformation which are located in the temperature range of 700–800[Formula: see text]C and the low strain rate range of 0.001–[Formula: see text]. The instability domain mainly occurs at the strain rate higher than [Formula: see text] for the whole deformation temperature, which should be avoided during practice. The flow softening mechanism of the alloy is determined to be continuous dynamic recrystallization and the unstable flow is caused by the macro-fracture and flow localization.
Influence of titanium on the microstructure and mechanical properties of foundry zinc alloy
[J].
High coverage H2 adsorption and dissociation on fcc Co surfaces from DFT and thermodynamics
[J].
Modeling hydrogen transport by dislocations
[J].
Effect of microstructure on the mechanical response of hydrogen-charged pure iron
[J].In this paper, we investigate how two different microstructures in pure iron affect the dislocation mobility in hydrogen-charged and non-charged samples by conducting stress-relaxation tests. The effective activation volume of the pure iron for both types of microstructures (cold-rolled and annealed samples) has been determined for both H-charged and uncharged material. Information about the dislocation structures formed during stress relaxation is provided by conducting TEM analysis. We employ a self-consistent kinetic Monte-Carlo (SCkMC) model of the ½ [111] screw dislocation in Fe to investigate how hydrogen affects the mobility and behavior of the dominant mobile dislocation in Fe at different stresses and H concentrations. The results from our simulations show the following: (i) at low stresses the deviation from the primary slip plane in the presence of H is lower than the deviation in the uncharged Fe. The deviation angle decreases with increasing H concentration; (ii) at higher shear stresses, the higher probability for kink-pair formation in the secondary (110) planes in the presence of H, leads to an enhanced deviation from the primary slip plane, which increases with increasing H concentration. We use the results of stress-relaxation tests and SCkMC simulations to propose an explanation for the formation of dislocation cell structures in pure and hydrogen charged Fe in the cold-rolled and annealed samples.
Effects of alloying element on hydrogen adsorption and diffusion on α-Fe (110) surfaces: First principles study
[J].Based on first principles density functional theory (DFT) methods, this study employed the Cambridge Serial Total Energy Package (CASTEP) module within Materials Studio (MS) software under the generalized gradient approximation to investigate the adsorption, diffusion behavior, and electronic properties of hydrogen atoms on α-Fe(110) and α-Fe(110)-Me (Mn, Cr, Ni, Mo) surfaces, including calculations of their adsorption energies and density of states (DOS). The results demonstrated that doping with alloy atoms Me increased the physical adsorption energy of H2 molecules on the surface. Specifically, Mo doping elevated the adsorption energy from −1.00825 eV to −0.70226 eV, with the largest relative change being 30.35%. After doping with Me, the chemical adsorption energy of two hydrogen atoms does not change significantly, among which doping with Cr results in a decrease in the chemical adsorption energy. Building on this, further analysis of the chemical adsorption of single atoms on the surface was conducted. By comparing the adsorption energy and the bond length between a hydrogen atom and iron/dopant metal atom, it was found that Mo doping has the greatest impact, increasing the bond length by 58.58%. Analysis of the DOS functions under different doping conditions validated the interaction between different alloy elements and H atoms. Simultaneously, simulations were carried out on the energy barrier crossed by H atoms diffusing into the metal interior. The results indicate that Ni doping facilitates the diffusion of H atoms, while Cr, Mn, and Mo hinder their diffusion, with Mo having the most significant effect, where its barrier is 21.88 times that of the undoped surface. This conclusion offers deep insights into the impact of different doping elements on hydrogen adsorption and diffusion, aiding in the design of materials resistant to hydrogen embrittlement.
Effects of alloying elements on hydrogen diffusion in iron
[J].
Engineering metal-carbide hydrogen traps in steels
[J].Hydrogen embrittlement reduces the durability of the structural steels required for the hydrogen economy. Understanding how hydrogen interacts with the materials plays a crucial role in managing the embrittlement problems. Theoretical models have indicated that carbon vacancies in metal carbide precipitates are effective hydrogen traps in steels. Increasing the number of carbon vacancies in individual metal carbides is important since the overall hydrogen trapping capacity can be leveraged by introducing abundant metal carbides in steels. To verify this concept, we compare a reference steel containing titanium carbides (TiCs), which lack carbon vacancies, with an experimental steel added with molybdenum (Mo), which form Ti-Mo carbides comprising more carbon vacancies than TiCs. We employ theoretical and experimental techniques to examine the hydrogen trapping behavior of the carbides, demonstrating adding Mo alters the hydrogen trapping mechanism, enabling hydrogen to access carbon vacancy traps within the carbides, leading to an increase in trapping capacity.
Effects of alloying elements on the solution and diffusion of oxygen at iron grain boundary investigated by first-principles study
[J].The effects of alloying elements (Si, Cr, Mo) on the solution and diffusion of oxygen (O) atoms at the grain boundary of iron (Fe) Σ5(310)/[001] are investigated by the simulations of ab initio density functional theory (DFT). It is found that Si, Mo and Cr prefer to segregate to the grain boundary, and further affect the solution and diffusion of O atoms at Fe grain boundaries. The segregated Cr promotes the solution of O, while Si and Mo inhibit the solution of O at the grain boundary. Meanwhile, Cr and Si accelerate the diffusion of O, and Mo retards the diffusion of O in the grain boundary. Further analysis indicates that the effects are closely related to the interactions between the alloying elements and O atoms, which are determined by the competition between the distortion of local structure and the charge transfer between local atoms. Finally, the effects of alloying elements on the O concentration distribution near the grain boundary are explored by employing the Langmuir–McLean models. This work not only provides insights into the effects of alloying elements on the solution and diffusion of O at grain boundaries, but also provides parameters of the atomic interactions for the initial oxidation simulation on a large scale, which relates to the growth of oxide in polycrystalline systems with various grain sizes at experimental temperatures.
Investigating the impact of alloying elements on hydrogen diffusion in Ti-based alloys via first-principles calculations
[J].
Ab initio molecular dynamics for open-shell transition metals
[J].
Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set
[J].
Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set
[J].
Generalized gradient approximation made simple
[J].
High-precision sampling for Brillouin-zone integration in metals
[J].
Special points for Brillouin-zone integrations
[J].
Improved tetrahedron method for Brillouin-zone integrations
[J].
LOBSTER: Local orbital projections, atomic charges, and chemical-bonding analysis from projector-augmented-wave-based density-functional theory
[J].We present an update on recently developed methodology and functionality in the computer program Local Orbital Basis Suite Toward Electronic-Structure Reconstruction (LOBSTER) for chemical-bonding analysis in periodic systems. LOBSTER is based on an analytic projection from projector-augmented wave (PAW) density-functional theory (DFT) computations (Maintz et al., J. Comput. Chem. 2013, 34, 2557), reconstructing chemical information in terms of local, auxiliary atomic orbitals and thereby opening the output of PAW-based DFT codes to chemical interpretation. We demonstrate how LOBSTER has been improved by taking into account time-reversal symmetry, thereby speeding up the DFT and LOBSTER calculations by a factor of 2. Over the recent years, the functionalities have also been continually expanded, including accurate projected densities of states (DOSs), crystal orbital Hamilton population (COHP) analysis, atomic and orbital charges, gross populations, and the recently introduced k-dependent COHP. The software is offered free-of-charge for non-commercial research.© 2020 Wiley Periodicals LLC.
Crystal orbital Hamilton populations (COHP): Energy-resolved visualization of chemical bonding in solids based on density-functional calculations
[J].
A climbing image nudged elastic band method for finding saddle points and minimum energy paths
[J].
Hydrogen in α-iron: Stress and diffusion
[J].
Investigation of interaction between α-Fe metal and H atom by ab-initio method
[J].
基于第一性原理的含空位α-Fe和H原子相互作用研究
[J].
/
| 〈 |
|
〉 |
