摘要:Al-Mg-Si系铝合金在氢环境服役时, 因遭受氢损伤而导致力学性能退化. Mg₂Si是Al-Mg-Si系铝合金主要的热处理强化相, 其力学性能直接决定着Al-Mg-Si系铝合金的强度. 本文采用基于密度泛函数理论的第一性原理计算方法, 研究了间隙H原子对Mg₂Si力学性能的影响. 首先计算了Mg₂Si的单晶体弹性常数C₁₁, C₁₂和C₄₄, 在通过Hill模型计算了多晶体的弹性模量、泊松比及硬度, 并进一步计算了H对Mg₂Si晶粒拉伸性能的影响, 最后对H原子掺杂引起的Mg₂Si晶体态密度变化进行了分析. 结果表明: H原子的引入显著降低Mg₂Si晶体的剪切模量和弹性模量, 从而使得体系强度、硬度降低, 韧性提高. 拉伸性能计算也表明H原子引起Mg₂Si晶体断裂强度降低、断裂延伸率提高. 态密度分析表明, Mg₂Si晶体中引入氢原子, 将引起Mg₂Si晶体由半导体性质向金属性质的转变, 从而造成体系强度、硬度的降低, 韧性增加. 本文计算结果可以为揭示以Mg₂Si为增强相的材料在氢环境中强度降低机理提供参考依据.
关键词: 第一性原理/
Mg₂Si/
弹性性质/
态密度

English Abstract

全文HTML

--> --> -->

Mg₂Si金属间化合物是Al-Mg-Si系铝合金主要的热处理强化相, 高弹性模量、高硬度的Mg₂Si颗粒弥散分布于铝合金基体中, 通过阻碍位错运动在Al-Mg-Si系铝合金形成弥散强化作用, 使铝合金强度硬度大幅度提高^[1,2]. 由于Mg₂Si的高熔点(1085 ℃)^[3], 被用于替代SiC陶瓷制备铝基复合材料, 且具有良好的切削加工性和成形性. 因此, Mg₂Si金属间化合物的研究备受关注^[4].
对于Mg₂Si金属间化合物的理论研究, 主要采用基于密度泛函数理论的第一性原理计算方法. 任玉艳等^[5]对Mg₂Si金属间化合物的结构稳定性、热力学和力学性能进行了系统的第一性原理计算研究, 理论计算值与实验值十分接近. 余本海等^[6]采用第一性原理的方法对Mg₂Si的同质异相体的晶体结构、电子结构及弹性性质进行了全面研究.
氢由于其原子半径极小, 容易在材料中扩散, 并引起材料性能的变化. 氢原子对晶体结构和力学性质影响的第一性原理计算已经成为目前研究的热点. 韩秀丽等^[7]采用第一性原理的方法对Ti-H体系的晶体结构和能量进行了研究, 氢原子的加入引起Ti晶体的晶格畸变, 降低了晶体的结合能. 氢在α-Fe和γ-Fe中掺杂均使晶体结构产生不同程度的畸变, 并降低了材料的塑性^[8,9]. 姚宝殿等^[10]的研究显示氢对Ti₃SiC₂的力学性质影响不大. 对α-Ti—H键重叠集居数的研究发现^[11], 氢减弱了近邻原子间的成键强度, 存在弱键效应; 而在β-Ti—H中, 氢增加了近邻钛原子间的键重叠集居数, 存在强键效应. 氢的加入也会使Nb晶体发生类似的电子结构变化^[12]. α-Fe₂O₃—H中, H和O之间成键, 费米面附近出现了新的杂化能级, 价带顶变宽, 同时带隙由于导带底下移而减小, 表明氢可以用来调节能带结构^[13].
本文采用基于密度泛函数理论的第一性原理方法, 构建Mg₂Si-H体系, 研究了氢原子浓度变化对Mg₂Si金属间化合物弹性性质的影响, 并对其力学性能的变化进行了深入分析.

Mg₂Si金属间化合物具有反萤石结构, 空间群号为Fm3m (No. 255), 为面心立方 (FCC)结构, 晶格常数a = b = c = 0.676 nm, α = β = γ = 90° ^[14]. Si原子占据顶角和面心位置, Mg原子位于Si原子组成的四面体间隙中形成边长为0.5a的简单立方结构, 余留体心位置的八面体间隙, 成为H等小半径原子优先固溶位置.
本文采用基于密度泛函数理论(DFT)和平面波赝势的CASTEP(Cambridge sequential total energy package)模块, 对不同H原子浓度的Mg₂Si-H体系进行模拟计算. 首先构建Mg₂Si晶胞, 几何优化后构建1 × 1 × 1, 1 × 1 × 2, 1 × 2 × 2的Mg₂Si超晶胞, 在体心八面体间隙中掺杂H原子, 形成H原子掺杂浓度分别为7.69 at.%, 4 at.%和2.05 at.%的Mg₂Si-H体系, 如图1所示. 对掺杂H原子后的超晶胞再进行几何优化, 得到稳定的Mg₂Si-H体系, 然后对Mg₂Si-H体系进行弹性性质计算. 为更直观地反映H原子掺杂对Mg₂Si力学性的影响, 分别计算了1 × 1 × 2的Mg₂Si超晶胞和1 × 1 × 2的Mg₂Si-H (4 at.%)超晶胞沿[001]方向进行拉伸的应力-应变曲线. 拉伸应变通过改变晶胞[001]方向的尺寸进行施加, 几何优化后获得对应的应力值.
图 1 H原子位于八面体间隙的Mg₂Si-H体系超晶胞结构模型
Figure1. Super-cell structural model of Mg₂Si-H system.

计算过程采用周期性边界条件, 结构优化和弹性性质计算过程中, 采用广义梯度近似(genaral gradient approximation, GGA)中的质子平衡方程(proton balance equation, PEB)来处理电子间的交互关联能, 势函数选用倒异空间中表述的超软赝势(ultrasoft pseudopotential). 计算中平面波截断能量(cut off energy)取400 eV, 最大内应力为1 × 10^–3 GPa, 原子最大位移为2 × 10^–4 nm, 能量收敛标准为1 × 10^–6 eV/atom. 第一布里渊区k-point矢量分别取9 × 9 × 9, 9 × 9 × 4和9 × 4 × 4.
在三向应力下, 材料共有36个弹性常数C_ij (i = 1—6, j = 1—6). Mg₂Si晶体为立方结构, 因此存在三个独立的弹性常数, 即C₁₁, C₁₂和C₄₄, 分别体现体系在[100], [111]和[011]三个方向的抗变形能力. 当C₁₁, C₁₂和C₄₄同时满足C₁₁–C₁₂ > 0, C₁₁ > 0, C₄₄ > 0和C₁₁+2C₁₂ > 0四个条件时, 合金体系才可以稳定存在, 满足力学性能稳定性条件^[15]. 立方晶系剪切模量C'可由C₁₁和C₁₂计算:

$C' = \frac{{{C_{11}} - {C_{12}}}}{2}.$

材料的力学性质在微观上可以通过弹性常数、体模量、杨氏模量、剪切模量、泊松比和体模量/剪切模量之比等体现. 多晶弹性模量通常可由单晶弹性常数计算获得, 计算方法可用Voigt模型和Reuss模型描述. 经验证, Voigt模型和Reuss模型所得模量通常为多晶模量的上限和下限, 因此Hill^[16,17]对两者进行数学平均, 发现所得值与实测值一致性更好. 因此, 实际应用中Hill模型更为常用. 因此本文中Mg₂Si体模量B_H, 剪切模量G_H和杨氏模量E的计算如下:

${B_{\rm{V}}} = {B_{\rm{R}}} = \frac{1}{3}\left( {{C_{11}} + 2{C_{12}}} \right),$

${G_{\rm{V}}} = \frac{1}{5}\left( {{C_{11}} - {C_{12}} + {C_{44}}} \right),$

${G_{\rm{R}}} = \frac{{5{C_{44}}\left( {{C_{11}} - {C_{12}}} \right)}}{{4{C_{44}} + 3\left( {{C_{11}} - {C_{12}}} \right)}},$

其中B_V和G_V分别为Voigt模型的体积模量和剪切模量, B_R和G_R分别为 Reuss模型的体积模量和剪切模量. Hill模型中体积模量B和剪切模量G分别为:

${B_{\rm{H}}} = \frac{{{B_{\rm{V}}} + {B_{\rm{R}}}}}{2},$

${G_{\rm{H}}} = \frac{{{G_{\rm{V}}} + {G_{\rm{R}}}}}{2}.$

杨氏模量E和泊松比ν可由B和G计算得到^[18,19]:

$E = \frac{{9{B_{\rm{H}}}{G_{\rm{H}}}}}{{3{B_{\rm{H}}} + {G_{\rm{H}}}}},$

$\nu = \frac{{3{B_{\rm{H}}} - 2{G_{\rm{H}}}}}{{2\left( {3{B_{\rm{H}}} + {G_{\rm{H}}}} \right)}}.$

通常情况下, B反映材料抗体积压缩能力; G反映材料抗剪切应变的能力, 且与硬度相关^[19]; E反映材料抗正应变的能力, 也反映材料的硬度, E越大, 材料的硬度越高^[20]. 剪切应力下, ν越小则稳定性越好. B与G的比值是材料工程应用中必须考虑的重要参数, 能反映材料的脆韧性, 以1.75作为临界值. 当 B/G > 1.75时材料呈韧性, 否则呈脆性^[20,21].
材料的硬度与剪切模量、体模量之间存在直接的对应关系^[22,23]. Miao等 ^[24]在通过第一性原理研究材料的力学性质中, 提出了硬度H与杨氏模量E和泊松比ν之间计算的经验公式:

$H = \frac{{\left( {1 - 2\nu } \right)E}}{{6\left( {1 + \nu } \right)}}.$

2

-->

3.1.晶体弹性常数

Mg₂Si-H体系晶体弹性常数C₁₁, C₁₂, C₄₄和C'计算结果如表1所列. 从表1中可以看出, 所有体系中C₁₁, C₁₂和C₄₄均同时满足C₁₁–C₁₂ > 0, C₁₁ > 0, C₄₄ > 0和C₁₁+2C₁₂ > 0四个条件, 说明本文所建立的Mg₂Si-H体系可以稳定存在, 满足力学性能稳定性条件^[15]. 本工作中的Mg₂Si的弹性常数C₁₁, C₁₂和C₄₄计算值分别为113.2, 22.2和45.2 GPa, 与其他****的理论计算结果^[6,20]相比, 最大计算误差分别为1.15%, 2.70%和0.88%, 计算精度较高; 与其他****的试验结果^[21]相比, 本工作中的, C₁₂和C₄₄的最大计算误差分别为10.16%, 14.61和6.80%, 说明本文计算方法具有较高的可靠性.

Phase	Source	Elastic constants/GPa
		C11	C12	C44	C'	C12–C44
Mg2Si	This work	113.2	22.2	45.2	45.5	–23.0
	Calculated [20]	113.7	22.8	43.5
	Calculated [6]	114.5	21.5	45.6
	Experimental [21]	126.0	26.0	48.5
Mg2Si-2.05 at.%H	This work	110.3	25.1	38.8	42.6	–13.7
Mg2Si-4 at.%H	This work	109.0	26.1	37.8	41.4	–11.7
Mg2Si-7.69 at.%H	This work	90.6	35.5	34.9	27.5	0.6

表1Mg₂Si-H体系晶体弹常数计算结果
Table1.Elastic constant of Mg₂Si-H system.

图2给出了弹性常数C₁₁, C₁₂, C₄₄和C'随氢原子浓度增加的变化趋势. C₁₁值反映了体系在[100]方向的抗变形能力. 从图2中可以看出, 随氢原子浓度从0增加到7.69 at.%过程中, C₁₁从113.2 GPa逐步减小到90.6 GPa, 降低了20%, 说明H原子的引入显著降低了Mg₂Si在[100]方向抵抗变形的能力. C₁₂反映了体系[111]方向的抗变形能力. 从图2中可看出, H原子浓度从0增加到7.69 at.%过程中, Mg₂Si的C₁₂值从22.2 GPa增加到35.5 GPa, 增加了59.9%, 说明氢原子的引入使得Mg₂Si抵抗[111]方向的抗变形能力明显提高. C₄₄反映了体系[011]方向的抗变形能力, 随H原子浓度从0增加到7.69 at.%, Mg₂Si的C₄₄从45.2 MPa降低到34.9 GPa, 降低了22.8%, 说明Mg₂Si在[011]方向抗变形能力随氢原子浓度的增加而有所降低. 经(1)式计算得到的单晶体剪切模量C'则随着氢浓度从0增加到7.69 at.%, 从45.5 GPa降低到27.5 GPa, 降低了39.5%, 说明H原子具有显著降低Mg₂Si单晶体的抗剪切变形能力. 此外, C₁₂–C₄₄值还可表征晶体的脆性, 当C₁₂–C₄₄ > 0时晶体为塑性, 当C₁₂–C₄₄ < 0晶体为脆性^[17]. 从表1可以看出, 当H原子浓度达到7.69 at.%时, C₁₂–C₄₄值从由负值增加到正值, Mg₂Si晶体发生了由脆性向塑性的转变.
图 2 弹性常数C₁₁, C₁₂, C₄₄和C'随氢原子浓度增加的变化趋势
Figure2. Relationships between C₁₁, C₁₂, C₄₄ and C' as a function of H atoms content.

2 -->

3.2.多晶体弹性模量

经(5)—(9)式计算的剪切模量(G_H)、体模量(B_H)、弹性模量(E)、泊松比(ν)、脆性(B_H/G_H)及硬度(H)如表2所列. 本工作中的G_H, B_H, E和ν值分别为45.3, 52.5, 105.6 GPa和0.16, 与其他****的理论计算结果^[6,20]相比, 本工作的最大计算误差分别为2.21%, 1.73%, 1.42%和4.29%, 计算精度较高; 与其他****试验结果^[21]相比, 本工作中的B_H和E的最大计算误差分别为12.38%和13.63%, 与试验结果同样十分接近, 说明本文计算方法具有较高的可靠性.

Phase	Source	Modulus
		GH/GPa	BH/GPa	E/GPa	ν	BH/GH	H/GPa
Mg2Si	This work	45.3	52.5	105.6	0.16	1.16	10.31
	Calculated [20]	46.2	52.5	107.1	0.16	—	—
	Calculated [6]	44.3	53.4	104.1	0.17	—	—
	Experimental [21]	—	59.0	120.0	—	—	—
Mg2Si-2.05 at.%H	This work	40.2	53.4	96.5	0.17	1.32	9.07
Mg2Si-4 at.%H	This work	39.2	53.7	94.7	0.18	1.37	8.56
Mg2Si-7.69 at.%H	This work	31.7	53.8	79.5	0.28	1.69	4.55

表2Mg₂Si-H体系模量
Table2.Modulus of Mg₂Si-H system.

图3给出了剪切模量(G_H)、体模量(B_H)、弹性模量(E)、泊松比(ν)、脆性(B_H/G_H)及硬度(H)随氢原子含量增加的变化趋势. 弹性模量E是表征材料刚度的重要参数, E值越高, 材料刚度越大. 本工作计算结果表明, 未引入氢原子时, Mg₂Si的弹性模量E 为105.6 GPa, 随着体系内H原子浓度从0增加到7.69 at.%, 弹性模量E逐步降低到79.5 GPa, 降低了24.7%. 切变模量G与弹性模量E类似, 随着体系氢原子浓度从0增加到7.69 at.%, 切变模量G从45.3 GPa降低到31.7 GPa, 降低了30%. 体模量B则随体系氢浓度的改变未表现出明显的变化. 弹性模量E和切变模量G的大幅度降低, 表明Mg₂Si-H体系内氢原子浓度的增加, 显著降低了材料的强度、硬度. 与α-Ti相似^[25], Mg₂Si也表现出H致软化现象. 硬度计算结果显示, 体系的硬度从不含H原子时的10.31 GPa逐步降低到7.69 at.% H浓度时的4.55 GPa, 降低了55.8%, 印证了弹性模量E和切变模量G的计算结果. 随着H原子浓度增加, 体系硬度降低, B/G值逐步增加, 泊松比ν也逐步增大. 当H原子浓度达到7.69 at.%时, 体系的泊松ν为0.28, 与普通低合金钢基本相当^[26]. 说明当大量H原子溶于Mg₂Si时, Mg₂Si将逐步失去其本身具有的高强度、高硬度和高脆性, 可能严重降低以Mg₂Si为增强相材料的强度.
图 3 剪切模量(G_H)、体模量(B_H)、弹性模量(E)、泊松比(ν)、脆性(B_H/G_H)及硬度(H)与体系H原子浓度的变化关系
Figure3. Shear moduli G, bulk moduli B, Young's moduli E, Poisson ratios ν, brittleness B_H/G_H and Hardness H as a function of H atoms content.

2 -->

3.3.拉伸性能

图4给出了1 × 1 × 2的Mg₂Si超晶胞和1 × 1 × 2的Mg₂Si-H(4 at.%)超晶胞沿[001]方向进行拉伸的应力-应变曲线. 从图中可以看出, 当Mg₂Si中掺杂4 at.%的H原子时, 其抗拉强度从15.21 GPa降低到14.12 GPa, 延伸率则从0.3增加到0.36, 说明H原子的掺杂使Mg₂Si的强度降低, 塑性增大. 此外, 掺杂4 at.% H原子的Mg₂Si在弹性阶段的曲线斜率明显比不掺杂时要小, 说明H原子的掺杂使Mg₂Si的弹性模量降低. 拉伸过程的应力-应变曲线结果与弹性模量计算结果得到很好地相互应证.
图 4 Mg₂Si 1 × 1 × 2超晶胞和4 at.% Mg₂Si-H超晶胞沿[001]方向进行拉伸的应力-应变曲线
Figure4. Stress vs. strain curves of 1 × 1 × 2 super cell with 4 at.% and without H during stretching along [001] direction

2 -->

3.4.态密度分析

为探索H原子对Mg₂Si力学性能影响微观机制, 对体系的态密度和分波态密度进行了分析, 如图5所示. 图5(a)为Mg₂Si晶体的电子态密度和分波态密度. 由图可见, Mg₂Si晶体的能带主要密集分布于费米能级两侧, 具有典型的半导体特性. 分波态密度显示, 在–9 — –6 eV能量范围内, Mg₂Si的态密度主要由Mg的3s和Si的3s态电子贡献, Si的3p态电子有少量贡献. 在–5—0 eV能量范围内, Mg₂Si的态密度主要由Mg的3p和Si的3s态电子贡献, Si的3p态电子有少量贡献. 在0—20 eV能量范围内, Mg₂Si的态密度主要由Si的3p和3s态电子贡献, Mg的3p态电子有少量贡献. 当体系中引入7.79 at.%的H原子时, H原子显著增加了Mg和Si原子S和P态电子的贡献率, 总态密度明显提高; 同时, 体系在0—20 eV能量范围内的态密度峰缩减到0—6 eV范围内, 变现为体系的禁带变窄, 如图5(b)所示. 于此同时, H原子的引入, 还引起了费米能级处态密度的增加, 态密度跨过费米能级区域. 因此, H原子的引入使体系具有了一定的金属性质, 从而引起材料强度、硬度的降低, 韧性增加.
图 5 (a) Mg₂Si 和 (b) Mg₂Si-7.69 at.%H 总态密度和分波态密度
Figure5. DOS and PDOS of (a) Mg2Si and (b) Mg₂Si-7.69 at.% H.

1) Mg₂Si-H体系内H原子浓度从0增加到7.69 at.%, 剪切模量和弹性模量分别下降24.7%和31.7%, 从而造成体系强度硬度降低, 韧性提高.
2)态密度分析表明, Mg₂Si晶体中引入氢原子, 将引起Mg₂Si晶体由半导体性质向金属性质的转变, 从而造成体系强度、硬度的降低, 韧性增加.