-
氢致裂纹是制约超高强度钢应用的关键问题, 掌握扩散氢的分布行为有助于弄清氢致裂纹的形成机理. 本文采用第一性原理方法计算了H原子占据 α-Fe晶格间隙和空位时的情况, 得到了晶体的稳定构型及能量, 并据此分析了H原子在晶格间隙和空位中的溶解倾向; 从Mulliken布居、电子密度分布、态密度分布等角度分析了H原子与 α-Fe晶体间隙和空位之间的相互作用. 结果表明: 间隙H原子倾向占据 α-Fe四面体间隙位, 其1s轨道电子与Fe的4s轨道电子呈微弱共轭杂化; 空位是强氢陷阱, H原子倾向占据空位内壁附近的等电荷面. 在真空0 K条件下单空位最多稳定溶解3个H原子, 且H原子之间未表现出自发形成H 2的倾向; 间隙和空位中的H原子溶入改变了Fe晶格内电子分布导致原子结合力弱化, 并在局部区域形成反键. 基于第一性原理能量计算结果开展热力学分析, 分析结果表明大多数情况下间隙H原子都是H主要的固溶形式, H平衡溶解度计算结果与实际符合良好.Hydrogen-induced cracking (HIC) is a key problem restricting the application of ultra-high strength steel. It is necessary to analyze the distribution of diffusible hydrogen to reveal the mechanism of HIC. The site occupation tendency of H in interstitial and vacancy positions are investigated by the ab-initio method, and the stable configuration and steady state energy are obtained. The solution tendency of H atom in interstitial and vacancy positions is analyzed based on the aforementioned results. Specifically, the Mulliken population, density of states, charge density difference are calculated and used to analyze the interaction between α-Fe metal and H atom. The results show that the dissolved H tends to occupy the interstitial sites of the body-centered cubic, the weak hybridization interaction between the interstitial hydrogen and its nearest neighbour Fe atom is contributed by the H 1s orbital and Fe 4s orbital. Vacancies can capture H atoms easily and H atoms tend to occupy the isoelectric surface near the inwall of the vacancies. A vacancy defect can hold up to three H atoms which are difficult to combine with each other to form H 2molecule by covalent bond. H atoms in vacancies and at interstitial positions change the charge distribution of the Fe lattice, which weakens the binding force of the atoms and forms anti-bonding orbital in local area. The proposed thermodynamical model allows the determining of the equilibrium vacancy and the dissolved H concentration for a given temperature and H chemical potential in the reservoir, and the calculated results are in good agreement with the actual results.
[1] [2] [3] [4] [5] [6] [7] [8] [9] [10] [11] [12] [13] [14] [15] [16] [17] [18] [19] [20] [21] [22] [23] [24] [25] [26] [27] [28] [29] [30] [31] -
晶体类型 a/Å b/Å c/Å α/(°) β/(°) γ/(°) V/Å 空间群 α-Fe+H(T-site) 5.6843 5.6612 5.6804 89.999 89.999 90.000 5.6843 115(P-4m2) α-Fe+H(O-site) 5.8019 5.6112 5.6112 90.000 90.000 90.000 5.8019 123(P4/mmm) 晶体类型 ${E_{{\rm{crystal}}}}$/eV ${E_{{\rm{ZP}}}}$/eV ${E_{{\rm{bind}}}}$/eV $E_{{\rm{form}}}$/eV $\Delta H_{{ {\rm{sol} } } }^{\rm{H} }$/eV α-Fe+H(T-site) –13861.050 0.246 5.182 0.390 0.390 α-Fe –13845.344 — 5.530 — 晶体类型 原子 轨道电荷占据数 总布居 净布居 s p d α-Fe-H(T-site) H 1.34 0 0 1.34 –0.34 Fe2, Fe4 0.62 0.66 6.65 7.93 0.07 Fe3, Fe11 0.62 0.67 6.65 7.94 0.06 Fe12, Fe10 0.65 0.72 6.62 7.99 0.01 Fe1, Fe9, Fe13, Fe14, Fe16 0.65 0.72 6.62 7.99 0.01 Fe8, Fe6 0.65 0.74 6.61 8.01 –0.01 α-Fe Fe 0.68 0.70 6.62 8.00 0 自由态 H 1.00 0 0 1.00 0 晶体类型 原子对 距离/Å 键布居 α-Fe-H (T-site) Fe2-H 1.6494 0.16 Fe3-H 1.6507 0.16 Fe2-Fe3 2.5558 –0.09 Fe2-Fe4 2.7286 –0.14 Fe8-Fe11 2.4783 0.17 Fe3-Fe12 2.4471 0.18 Fe8-Fe16 2.8401 0.05 Fe7-Fe8 2.4601 0.15 α-Fe-H (O-site) Fe2-Fe4 2.6287 –0.28 α-Fe Fe-Fe 2.4400 0.14 Fe-Fe 2.8174 0.06 晶体类型 a/Å b/Å c/Å α/(°) β/(°) γ/(°) V/Å 空间群 α-Fe+Vac 5.6033 5.6033 5.6033 90.000 90.000 90.000 175.923 221${\rm{(}}Pm\overline {{\rm{3}}m} )$ α-Fe+(Vac-1H) 5.6321 5.6103 5.6103 90.000 90.001 89.999 177.270 99${\rm{(}}P{\rm{4}}MM)$ α-Fe+(Vac-2H) 5.6285 5.6285 5.6484 90.000 90.000 90.000 178.940 123(P4/MMM) α-Fe+(Vac-3H) 5.6297 5.6598 5.6853 90.004 90.011 90.002 181.154 25(PMM2) α-Fe+(Vac-4H) 5.6727 5.6943 5.6723 89.966 90.540 89.973 183.221 38(AMM2) α-Fe+(Vac-5H) 5.6905 5.7086 5.7093 90.000 90.004 90.002 185.467 99${\rm{(}}P{\rm{4}}MM)$ α-Fe+(Vac-6H) 5.7407 5.7270 5.7208 89.433 89.691 89.692 188.064 5(C2) 晶体类型 ${E_{{\rm{crystal}}}}$/eV ${E_{{\rm{ZP}}}}$/eV ${E_{{\rm{bind}}}}$/eV $E_{{\rm{form}}}$/eV $\Delta H_{_{ {\rm{sol} } } }^{\rm{H} }$/eV α–Fe+Vac –12977.593 — 5.369 2.416 — α–Fe+(Vac-1H) –12993.933 0.141 5.055 1.928 –0.347 α–Fe+(Vac-2H) –13010.262 0.295 4.777 1.450 –0.324 α–Fe+(Vac-3H) –13026.330 0.478 4.513 1.234 –0.034 α–Fe+(Vac-4H) –13042.359 0.670 4.275 1.056 0.014 α–Fe+(Vac-5H) –13058.299 0.889 4.055 0.968 0.131 α–Fe+(Vac-6H) –13073.995 1.149 3.842 1.123 0.438 晶体类型 $E_{{\rm{trap}}}^{\rm{H}}$/eV 不考虑${E_{{\rm{ZP}}}}$ 考虑${E_{{\rm{ZP}}}}$ α-Fe+(Vac-1H) 0.633 0.778 α-Fe+(Vac-2H) 0.623 0.627 α-Fe+(Vac-3H) 0.361 0.211 α-Fe+(Vac-4H) 0.322 –0.011 α-Fe+(Vac-5H) 0.227 –0.297 α-Fe+(Vac-6H) –0.028 –0.772 晶体类型 原子 轨道电荷占据数 总布居 净布居 s p d α-Fe+Vac Fe7, Fe11, Fe13 0.74 0.72 6.68 8.13 –0.13 Fe3, Fe5, Fe9 0.68 0.73 6.63 8.04 –0.04 其余Fe原子 0.66 0.67 6.60 7.93 0.07 α-Fe+(Vac-2H) Fe13 0.71 0.71 6.73 8.15 –0.15 Fe9 0.66 0.70 6.62 7.98 0.02 Fe11 0.71 0.71 6.65 8.07 –0.07 α-Fe+Vac—α-Fe+(Vac-6H) H1 1.20—1.22 0 0 1.20—1.22 –0.20— –0.22 自由态 H 1.00 0 0 1.00 0 晶体类型 ${E_{{\rm{crystal}}}}$/eV ${E_{{\rm{ZP}}}}$/eV ${E_{{\rm{bind}}}}$/eV $E_{{\rm{form}}}$/eV $\Delta H_{{ {\rm{sol} } } }^{\rm{H} }$/eV α–Fe+Vac+H(T-site) –12993.354 — –12993.353 0.248—0.250 5.012 2.755—2.756 0.339—0.340 温度/K μH/eV cH/% 计算值 实验值 298.15 –0.239 2.08 × 10–2 4.41 × 10–2
2.88 × 10–2, 其中晶格溶H占总扩散H含量的43%[2] -
[1] [2] [3] [4] [5] [6] [7] [8] [9] [10] [11] [12] [13] [14] [15] [16] [17] [18] [19] [20] [21] [22] [23] [24] [25] [26] [27] [28] [29] [30] [31]
计量
- 文章访问数:8518
- PDF下载量:136
- 被引次数:0