-
过渡金属双氮化合物具有诸多优异的物理性质,如较高的热力学及化学稳定性、超高硬度及高电热导率等[1–3],因此激发了众多学者的研究兴趣。5d过渡金属(尤其是Re、Os、Ir等)具有相当大的体模量[4],但其中的金属键又使其剪切模量很低[5]。实验上,利用高压技术将N原子嵌入过渡金属中可以形成局域而定向的共价键,进而极大地增加了材料的硬度。近年来,实验上PtN2、OsN2、IrN2等[6–8]过渡金属双氮化合物相继被合成,但到目前为止,关于ReN2的实验报道还很少。为了给ReN2的实验合成提供技术支持,很多学者进行了相关理论研究。
2008年,Zhao等[9]基于ReN2相关化合物的研究成果,在密度泛函理论的基础上,采用第一性原理计算方法中的超软赝势平面波法计算了ReN2的Fm-3m、Pa-3、P42/mnm、Pmmn 4种晶体结构在零温零压下的平衡结构参数,并通过能量-体积曲线得出:P42/mnm-ReN2结构在零温零压下最稳定;2009年,Li等[10]提出ReN2可能存在一种新的正交晶系结构Pbcn-ReN2,并由焓差-压强图得出Pbcn-ReN2结构在常温常压下最稳定;2010年,Du等[11]调研发现P4/mmm-ReB2结构具有非常大的体模量与剪切模量,因此通过维也纳从头算模拟程序包VASP的全电子投影缀加平面波(PAW)方法对ReN2的类ReB2晶体结构P4/mmm进行了计算;2012年,Kawamura等[12]通过ReCl5、Li3N、NaCl在高压下的复分解反应,实验上得到了ReN2的MoS2结构;2013年,Wang等[13]从弹性、热力学、X射线衍射(XRD)谱及化学键等方面研究得出:MoS2-ReN2结构并不稳定,并利用CALYPSO(Crystal Structure Analysis by PSO Code)[14]程序对ReN2的最稳定结构进行了预测,结果发现ReN2的最稳定结构为C2/m结构,并预估当压强高于130 GPa时,C2/m结构可能会相变为P4/mbm结构。
另外,弹性常数是晶体非常重要的一个物理量。通过弹性常数可以得知晶体材料的诸多热力学参数,如热膨胀系数、比热容、德拜温度、热-弹性应力、声速及断裂韧性等,弹性常数还能够解释原子间成键、晶体声子谱等固态现象[15–16]。然而,关于ReN2在高压下的弹性常数至今还没有人研究。因此,本工作将着重研究C2/m-ReN2在高压下的弹性性质。
全文HTML
-
采用基于密度泛函理论的Materials Studio软件包中的CASTEP模块[17]计算电子结构的总能量。其中计算方法为Vanderbilt超软赝势平面波法[18],交换相关势为Perdew-Burke-Ernzerhof(PBE)形式的广义梯度近似(General Gradient Approximate,GGA)[19–20]。原子赝势的电子组态为:Re 5d56s2,N 2s22p3,几何优化为 Brodyden-Fletcher-Goldfarb-Shanno(BFGS)算法[21]。为增加计算精确度,平面波基函数的截断能设置为500 eV,根据Monkhorst-Pack方法[22]测试得出的布里渊空间的k点取样值设置为11×11×9,总能自洽收敛精度设置为1 μeV/cell,应变前后的总能差值设置为小于10 μeV,离子最大Hellmanne-Feynman力设置为小于0.3 eV/nm,离子位移最大值设置为小于0.1 pm,应力最大值设置为小于50 MPa。以上参数相对不同体积都能得到完全自洽的总能,并经过测试后给出了相对精确的结果。
-
弹性常数的计算采用了同时适用于声速计算的非体积守恒张力法。有限张力变量Cijkl,即弹性常数,定义为[23–24]
式中:
${\sigma _{ij}}$ 、ekl分别指施加的压力张量、欧拉应变张量。X、x则代表发生应变前后的坐标。各向同性压力相关公式为[24–25]式中:cijkl表示对无穷小欧拉应变张力的二阶导数。21个独立组元才能构成四阶张量,但由于晶体具有对称性,因此一般晶体的独立组元(即弹性常数的个数)常小于21个。单斜晶系有13个独立组元,分别为C11、C22、C33、C44、C55、C66、C12、C13、C15、C23、C25、C35和C46。
1.1. 电子结构总能计算
1.2. 弹性常数
-
寻找零温零压下晶体的平衡结构并了解其结构参数对于探索晶体原子间的相互作用具有重要意义。本研究先对ReN2的单斜晶系结构C2/m进行加压几何优化,从0~80 GPa每间隔5 GPa优化一次,优化过程中晶体结构的原子坐标、键长、键角等结构参数会发生弛豫,由此得到不同原胞体积下的总能量。然后将优化得到的各E-V数据输入Birch-Murnaghan状态方程[26],进而拟合得到C2/m-ReN2在平衡状态下的各个晶格参数(a、b、c及
$\beta$ )、平衡体积V0、体模量B0以及体模量对压强的一阶导数$B_0'$ 。表1列出了计算得到的结果及其他理论计算值[13, 27]。对比发现,计算结果与其他理论参考值符合得很好。 -
晶体结构的弹性性质与晶体的力学、动力学性质及热力学性质等息息相关,了解高压下晶体的弹性性质对于晶体基础性质的研究及探索其工业应用具有重要价值。将零温零压下计算得到的弹性常数代入单斜晶系结构的力学稳定性判定条件
通过验证发现C2/m-ReN2满足力学稳定性条件。另外,高压下的理论值Cij也符合上面所有条件,进而也保证了C2/m-ReN2的高压力学稳定性。图1为C2/m-ReN2的各个弹性常数随压强p的变化关系。由图可知:随着压强的增大,C11、C22、C33、C44、C55、C66、C12、C13、C23和C46单调递增,其中C11、C22和C33的增长最显著,C12、C13和C23增长相对缓慢,C44、C55、C66和C46增长不太明显,C15、C25和C35则随压强的增加而缓慢单调递减。
将各个Cij代入Watt[28]给出的计算BV、BR以及GV、GR的方程,再结合Voigt-Reuss-Hill均值方法[29],便可计算得到C2/m-ReN2的体模量B和剪切模量G,相关公式为
式中:
$ a \,=\, {C_{33}}{C_{55}} \,-\, C_{35}^2,\,\,\,\, b = {C_{23}}{C_{55}} \,-\, {C_{25}}{C_{35}},\,\,\,\, c = {C_{13}}{C_{35}} - {C_{15}}{C_{33}},\,\,\,\, d = {C_{13}}{C_{55}} - {C_{15}}{C_{35}},\,\,\,\, e = {C_{13}}{C_{25}} - {C_{15}}{C_{23}},\,\, f =$ ${C_{11}}\left( {{C_{22}}{C_{55}} - C_{25}^2} \right) - {C_{12}}\left( {{C_{12}}{C_{55}} - {C_{15}}{C_{25}}} \right)+ {C_{15}}\left( {{C_{12}}{C_{25}} - {C_{15}}{C_{22}}} \right) + {C_{25}}\left( {{C_{23}}{C_{35}} - {C_{25}}{C_{33}}} \right),\,g = {C_{11}}{C_{22}}{C_{33}} -{C_{11}}C_{23}^2 - {C_{22}}C_{13}^2 -$ $ {C_{33}}C_{12}^2 \,\,+\,\, 2{C_{12}}{C_{13}}{C_{23}},\,\, \varOmega \,\,=\,\, 2[{C_{15}}{C_{25}}\left( {{C_{33}}{C_{12}} \,\,-\,\,{C_{13}}{C_{23}}} \right) \,\,+\,\, {C_{15}}{C_{35}}\left( {{C_{22}}{C_{13}} \,\,-\,\, {C_{12}}{C_{23}}} \right) \,\,+\,\, {C_{25}}{C_{35}}\left( {{C_{11}}{C_{23}} \,\,-\, \,{C_{12}}{C_{13}}}\,\right)]\, - $ $\,[C_{15}^2\left( {{C_{22}}{C_{33}} - C_{23}^2} \right)+C_{25}^2\left( {{C_{11}}{C_{33}} - C_{13}^2} \right) + C_{35}^2\left( {{C_{11}}{C_{22}} - C_{12}^2} \right)] + g{C_{55}}$ 。将计算得到的体模量B和剪切模量G代入杨氏模量E及泊松比
$\sigma $ 的定义式计算得到零温零压下C2/m-ReN2结构的体模量B、剪切模量G、杨氏模量E及泊松比
$\sigma $ ,见表2,对比发现,与其他理论值[13, 27]相吻合。图2给出了高压下C2/m-ReN2的体模量B、剪切模量G、杨氏模量E随压强p的变化关系。由图2可知,这3个物理量在高压0~80 GPa范围内均呈线性增长趋势。
据Pugh[30]的理论研究,B/G可作为鉴别晶体韧脆性的标准:当B/G小于1.75时,晶体材料表现为脆性;当B/G超过1.75时,晶体材料表现为韧性。零温零压下计算得到C2/m-ReN2的B/G约为1.53,因此零温零压下C2/m-ReN2表现为脆性。另外,Frantsevich等[31]提出利用泊松比
$\sigma $ 亦可辨别材料的韧脆性:当$\sigma $ <1/3时,材料为脆性;当$\sigma $ >1/3时为韧性。由计算得到的泊松比为0.23可知,这两种方法对C2/m-ReN2在零温零压下的韧脆性预测一致,并且C2/m-ReN2在0~80 GPa范围内仍表现为脆性(见图3)。另外,还可以得到C2/m-ReN2的压缩波声速vp、剪切波声速vs及平均声速vm,所用公式为
图 4为0 K条件下C2/m-ReN2的各个声速(vp、vs、vm)与压强p的变化关系。由图4可知,vp、vs和vm都随压强的增加单调递增。
根据德拜温度
$\varTheta $ 可以得知固体的比热容、热膨胀系数、电导率、热导率等诸多物理性质。将计算得到的平均声速vm代入式中:h为普朗克常数,kB为玻尔兹曼常数,n为原子数,NA为阿伏加德罗常数,
$\rho $ 为密度,M为相对分子质量。由图2可知,德拜温度$\varTheta $ 随压强的增长缓慢递增。 -
硬度是判断晶体材料抵抗渗透、磨损、变形程度的物理量,与晶体的弹性模量密切相关[32]。由零温零压下C2/m-ReN2的剪切模量(242 GPa)可估测其具有较大的硬度值。2012年,Tian等[33]得到了修正后的维氏硬度公式
式中:k=G/B。根据(14)式,计算得到C2/m-ReN2的硬度值为27.66 GPa(见表2),由此可预见C2/m-ReN2的超硬不可压缩性。
2.1. 晶体结构参数
2.2. 高压弹性
2.3. 硬 度
-
应用第一性原理计算方法,理论上计算并分析了零温零压下晶体C2/m-ReN2的结构性质,并首次研究了高压下C2/m-ReN2的力学稳定性及弹性性质。其中,具体分析了C2/m-ReN2的弹性常数Cij、弹性模量(B、G、E)、德拜温度
$\varTheta $ 、声速(vp、vs、vm)随压强的变化关系。研究发现,除个别弹性常数外,这些物理量都随压强的增加而增加。此外,还根据比值B/G及泊松比$\sigma $ 估测了C2/m-ReN2的韧脆性,并理论计算了C2/m-ReN2的维氏硬度。