搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

三元Nb系和Ta系硼碳化物稳定性和物理性能的第一性原理研究

胡前库 秦双红 吴庆华 李丹丹 张斌 袁文凤 王李波 周爱国

三元Nb系和Ta系硼碳化物稳定性和物理性能的第一性原理研究

胡前库, 秦双红, 吴庆华, 李丹丹, 张斌, 袁文凤, 王李波, 周爱国
PDF
HTML
导出引用
  • 过渡金属轻元素化合物是高硬度材料的潜在候选. 以往研究多集中在二元过渡金属硼化物、碳化物和氮化物, 三元相的研究则相对较少. 本文通过提炼和堆垛已知相Nb3B3C (Ta3B3C)和Nb4B3C2 (Ta4B3C2)的结构基元, 构建不同组分的Nb-B-C和Ta-B-C三元相结构模型, 采用第一性原理计算方法, 计算所建结构的形成焓、声子谱和弹性常数, 通过判断其热力学、动力学和力学稳定性, 绘出了三元Nb-B-C和Ta-B-C相图, 成功预测了5种Nb-B-C和6种Ta-B-C三元稳定相. 力学和电学性能计算结果显示Nb-B-C和Ta-B-C三元稳定相均为高硬度导电材料, 硬度大约为25 GPa.
      通信作者: 周爱国, zhouag@hpu.edu.cn
    • 基金项目: 国家级-国家自然科学基金(51472075,51772077)
    [1]

    Tian Y J, Xu B, Zhao Z S 2012 Int. J. Refract. Met. Hard Mater. 33 93

    [2]

    包括, 马帅领, 徐春红, 崔田 2017 物理学报 66 036104

    Bao K, Ma S L, Xu C H, Cui T 2017 Acta Phys. Sin. 66 036104

    [3]

    Zhou X F, Sun J, Fan Y X, Chen J, Wang H T, Guo X J, He J L, Tian Y J 2007 Phys. Rev. B 76 100101

    [4]

    Wu Q H, Hu Q K, Hou Y M, Wang H Y, Zhou A G, Wang L B 2018 J. Phys. Condens. Matter 30 385402

    [5]

    Tian Y J, Xu B, Yu D L, Ma Y M, Wang Y B, Jiang Y B, Hu W T, Tang C C, Gao Y F, Luo K, Zhao Z S, Wang L M, Wen B, He J L, Liu Z Y 2013 Nature 493 385

    [6]

    Huang Q, Yu D L, Xu B, Hu W T, Ma Y M, Wang Y B, Zhao Z S, Wen B, He J L, Liu Z Y, Tian Y J 2014 Nature 510 250

    [7]

    徐波, 田永君 2017 物理学报 66 036201

    Xu B, Tian Y J 2017 Acta Phys. Sin. 66 036201

    [8]

    Wu Q H, Hu Q K, Hou Y M, Wang H Y, Zhou A G, Wang L B, Cao G H 2018 Mater. Des. 140 45

    [9]

    Cumberland R W, Weinberger M B, Gilman J J, Clark S M, Tolbert S H, Kaner R B 2005 J. Am. Chem. Soc. 127 7264

    [10]

    Chung H Y, Weinberger M B, Levine J B, Kavner A, Yang J M, Tolbert S H, Kaner R B 2007 Science 316 436

    [11]

    Gregoryanz E, Sanloup C, Somayazulu M, Badro J, Fiquet G, Mao H K, Hemley R J 2004 Nat. Mater. 3 294

    [12]

    Young A F, Sanloup C, Gregoryanz E, Scandolo S, Hemley R J, Mao H K 2006 Phys. Rev. Lett. 96 155501

    [13]

    Ivanovskii A L 2012 Prog. Mater. Sci. 57 184

    [14]

    陶强, 马帅领, 崔田, 朱品文 2017 物理学报 66 036103

    Tao Q, Ma S L, Cui T, Zhu P W 2017 Acta Phys. Sin. 66 036103

    [15]

    Hillebrecht H, Gebhardt K 2001 Angew. Chem. Int. Ed. 40 1445

    [16]

    胡前库, 侯一鸣, 吴庆华, 秦双红, 王李波, 周爱国 2019 物理学报 68 096201

    Hu Q K, Hou Y M, Wu Q H, Qin S H, Wang L B, Zhou A G 2019 Acta Phys. Sin. 68 096201

    [17]

    Wang P F, Weng M Y, Xiao Y, Hu Z X, Li Q H, Li M, Wang Y D, Chen X, Yang X N, Wen Y R, Yin Y X, Yu X Q, Xiao Y G, Zheng J X, Wan L J, Pan F, Guo Y G 2019 Adv. Mater. 31 1903483

    [18]

    Xiao W J, Xin C, Li S B, Jie J S, Gu Y, Zheng J X, Pan F 2018 J. Mater. Chem. A 6 9893

    [19]

    Kresse G, Furthmüller J 1996 Phys. Rev. B 54 11169

    [20]

    Perdew J P, Burke K, Ernzerhof M 1996 Phys. Rev. Lett. 77 3865

    [21]

    Le Page Y, Saxe P 2002 Phys. Rev. B 65 104104

    [22]

    Togo A, Tanaka I 2015 Scr. Mater. 108 1

    [23]

    Togo A, Chaput L, Tanaka I, Hug G 2010 Phys. Rev. B 81 174301

    [24]

    Mouhat F, Coudert F X 2014 Phys. Rev. B 90 224104

    [25]

    Wu Z J, Zhao E J, Xiang H P, Hao X F, Liu X J, Meng J 2007 Phys. Rev. B 76 054115

    [26]

    Pugh S F 1954 Philos. Mag. 45 823

    [27]

    Chen X Q, Niu H Y, Li D Z, Li Y Y 2011 Intermetallics 19 1275

  • 图 1  (a), (b) Ta3B3C; (c) Ta4B3C2; (d) Ta3BC2; (e) Ta6B4C3; (f) Ta7B4C4; (g) Ta7B6C3的晶体结构. 棕球: Ta原子; 蓝球: B原子; 粉球: C原子. Ta6B三棱柱和Ta6C八面体分别用绿色和褐色表示

    Fig. 1.  The crystal structures of (a), (b) Ta3B3C; (c) Ta4B3C2; (d) Ta3BC2; (e) Ta6B4C3; (f) Ta7B4C4; (g) Ta7B6C3. The light brown, blue and pink spheres represent Ta, B, and C atoms, respectively. The Ta6B triangular prisms and Ta6C octahedrons are painted green and dark brown.

    图 2  (a) Nb-B-C和(b) Ta-B-C三元相图. 红色, 稳定相; 蓝色, 亚稳相; 绿色, 不稳定相

    Fig. 2.  Ternary phase diagrams of (a) Nb-B-C and (b) Ta-B-C. Red, stable; blue, metastable; green, unstable.

    图 3  不同温度下 (a) Nb-B-C和(b) Ta-B-C三元相分别和其相应最稳定竞争组合相的自由能之差

    Fig. 3.  Energy differences of (a) Nb-B-C and (b) Ta-B-C ternary phases with respect to their most competing phases as a function of temperature.

    图 4  Nb-B-C和Ta-B-C三元相的声子色散曲线

    Fig. 4.  Phonon dispersion curves of Nb-B-C and Ta-B-C ternary phases.

    图 5  Nb-B-C和Ta-B-C三元相的态密度图

    Fig. 5.  Density of states of Nb-B-C and Ta-B-C ternary phases.

    表 1  不同成分Nb(m + n + 2)B(2m + 2)Cn和Ta(m + n + 2)B(2m + 2)Cn晶体的结构参数

    Table 1.  Structural parameters of Nb(m + n + 2)B(2m + 2)Cn and Ta(m + n + 2)B(2m + 2)Cn crystals.

    mn空间群模型晶格参数/Å模型晶格参数/Å
    abcabc
    01CmmmNb3B2C3.25413.8083.141Ta3B2C3.24013.6973.127
    02CmcmNb2BC3.23518.3303.153Ta2BC3.22018.1653.140
    03CmmmNb5B2C33.22522.9033.153Ta5B2C33.19922.6593.138
    04CmcmNb3BC23.21427.3763.156Ta3BC23.19827.1323.150
    11PmmmNb4B4C3.29018.9943.145Ta4B4C3.27718.8783.127
    12ImmmNb5B4C23.26723.6003.150Ta5B4C23.24823.3773.138
    13PmmmNb6B4C33.24328.0283.154Ta6B4C33.22527.8723.141
    14ImmmNb7B4C43.24232.5453.158Ta7B4C43.22432.3153.147
    21CmmmNb5B6C3.30224.4143.134Ta5B6C3.28924.2083.122
    22CmcmNb3B3C3.28428.8773.144Ta3B3C3.26728.6883.133
    23CmmmNb7B6C33.26433.3643.148Ta7B6C33.24633.1643.136
    24CmcmNb4B3C23.25737.8743.153Ta4B3C23.24337.6093.141
    31PmmmNb6B8C3.30914.8893.137Ta6B8C3.29814.7883.122
    32ImmmNb7B8C23.29034.2473.144Ta7B8C23.27634.0073.131
    33PmmmNb8B8C33.27619.3503.148Ta8B8C33.25819.2353.135
    34ImmmNb9B8C43.26843.2553.151Ta9B8C43.25242.9773.138
    41CmmmNb7B10C3.31235.1923.131Ta7B10C3.29934.9903.116
    42CmcmNb4B5C3.29639.6943.139Ta4B5C3.28039.4413.125
    43CmmmNb9B10C33.28144.2063.142Ta9B10C33.26343.9243.130
    44CmcmNb5B5C23.27348.7293.145Ta5B5C23.25748.4003.134
    下载: 导出CSV

    表 2  不同成分Nb-B-C相和Ta-B-C相的形成焓 (单位: eV/atom), $ \Delta{H}_{\rm{elements}} $表示单质为反应物, $ \Delta{H}_{\rm{comp}} $表示最稳定竞争组合为反应物

    Table 2.  Calculated formation enthalpies of different Nb-B-C and Ta-B-C phases (in eV/atom).$ \Delta{H}_{\rm{elements}} $ represents the elements as the reactants, and $\Delta{H}_{\rm{comp}}$ indicates the most stable composite as the reactants.

    Phases$ \Delta{H}_{\rm{elements}} $$ \Delta{H}_{\rm{comp}} $最稳定竞争组合Phases$ \Delta{H}_{\rm{elements}} $$ \Delta{H}_{\rm{comp}} $最稳定竞争组合
    Nb3B2C–0.6200.070Nb3B4 + 6NbB + Nb6C5 = 5Nb3B2CTa3B2C–0.6510.086Ta3BC2 + 3TaB = 2Ta3B2C
    Nb2BC–0.6190.029Nb3B4 + NbB + Nb6C5 = 5Nb2BCTa2BC–0.6640.035Ta3BC2 + TaB = 2Ta2BC
    Nb5B2C3–0.5860.0363Nb3B4 + Nb7B4C4 + 4Nb6C5 = 8Nb5B2C3Ta5B2C3–0.6550.0213Ta3BC2 + TaB = 2Ta5B2C3
    Nb3BC2–0.5860.019Nb3B4 + 3Nb7B4C4 + 4Nb6C5 = 16Nb3BC2Ta3BC2–0.660–0.002TaB + 2TaC = Ta3BC2
    Nb4B4C–0.6790.0303Nb3B4 + Nb7B4C4 = 4Nb4B4CTa4B4C–0.6910.044Ta7B4C4 + 3Ta3B4 = 4Ta4B4C
    Nb5B4C2–0.6680.006Nb3B4 + Nb7B4C4 = 2Nb5B4C2Ta5B4C2–0.6940.019Ta7B4C4 + Ta3B4 = 2Ta5B4C2
    Nb6B4C3–0.6450.005Nb3B4 + 3Nb7B4C4 = 4Nb6B4C3Ta6B4C3–0.6930.0043Ta7B4C4 + Ta3B4 = 4Ta6B4C3
    Nb7B4C4–0.632–0.0063Nb3B4 + 2C + 2Nb6C5 = 3Nb7B4C4Ta7B4C4–0.685–0.0173Ta3B4 + 4TaC = Ta7B4C4
    Nb5B6C–0.6970.0153Nb3B4 + C + 2Nb3B3C = 3Nb5B6CTa5B6C–0.6970.024C + Ta5B6 = Ta5B6C
    Nb3B3C–0.685–0.0013Nb3B4 + C + 3Nb4B3C2 = 7Nb3B3CTa3B3C–0.6990.0103Ta7B4C4 + 9Ta3B4 + 4C = 16Ta3B3C
    Nb7B6C3–0.6640.0005Nb3B3C + Nb4B3C2 = Nb7B6C3Ta7B6C3–0.6950.00085Ta7B4C4 + 7Ta3B4 + 4C = 8Ta7B6C3
    Nb4B3C2–0.648–0.0015Nb3B4 + 4C + 7Nb7B4C4 = 16Nb4B3C2Ta4B3C2–0.6840.0027Ta7B4C4 + 5Ta3B4 + 4C = 16Ta4B3C2
    Nb6B8C–0.6950.0192Nb3B4 + C = Nb6B8CTa6B8C–0.6850.0342Ta3B4 + C = Ta6B8C
    Nb7B8C2–0.6830.0083Nb3B4 + 2C + 4Nb3B3C = 3Nb7B8C2Ta7B8C2–0.6860.020Ta7B4C4 + 7Ta3B4 + 4C = 4Ta7B8C2
    Nb8B8C3–0.6650.008C + 8Nb3B3C = 3Nb8B8C3Ta8B8C3–0.6840.012Ta7B4C4 + 3Ta3B4 + 2C = 2Ta8B8C3
    Nb9B8C4–0.6510.008C + 5Nb3B3C + 3Nb4B3C2 = 3Nb9B8C4Ta9B8C4–0.6750.0133Ta7B4C4 + 5Ta3B4 + 4C = 4Ta9B8C4
    Nb7B10C–0.6930.021C + 2Nb2B3 + Nb3B4 = Nb7B10CTa7B10C–0.6770.030TaB2 + 2Ta3B4 + C = Ta7B10C
    Nb4B5C–0.6840.0112C + Nb3B3C + 3Nb3B4 = 3Nb4B5CTa4B5C–0.6790.026Ta7B4C4 + 19Ta3B4 + 12C = 16Ta4B5C
    Nb9B10C3–0.6680.012C + 2Nb3B3C + Nb3B4 = Nb9B10C3Ta9B10C3–0.6770.0193Ta7B4C4 + 17Ta3B4 + 12C = 8Ta9B10C3
    Nb5B5C2–0.6550.011C + 5Nb3B3C = 3Nb5B5C2Ta5B5C2–0.6700.0185Ta7B4C4 + 15Ta3B4 + 12C = 16Ta5B5C2
    下载: 导出CSV

    表 3  Nb-B-C和Ta-B-C三元相的弹性常数Cij、体模量B、剪切模量 G和维氏硬度Hv (单位: GPa)

    Table 3.  Elastic constants Cij, bulk modulus B, shear modulus G, Vickers hardness Hv of Nb-B-C and Ta-B-C ternary phases (in GPa).

    结构弹性常数力学性能a硬度
    C11C22C33C44C55C66C12C13C23BGB/GHChenHTian
    Nb3B3C544.3479.8522.8181.5171.9245.3170.9132.9162.2275.3189.71.4524.824.7
    Nb4B3C2551.5499.2548.5184.0175.1257.1183.2132.7157.8282.9195.81.4425.525.4
    Nb6B4C3533.3493.8548.1174.9161.3255.2175.4138.9151.7278.5189.51.4724.424.3
    Nb7B4C4535.9505.9526.4172.2161.3259.1184.0142.8152.6280.6188.31.4923.923.8
    Nb7B6C3553.1494.5563.2188.7179.6255.6176.4132.1157.7282.5198.91.4226.326.2
    Ta3B3C569.6514.4563.5194.1180.0261.8187.1147.3173.9295.9200.81.4725.325.3
    Ta4B3C2581.1535.3602.1197.3185.1275.8200.3146.0170.2305.7209.01.4626.226.2
    Ta3BC2550.0547.7550.0159.8159.5292.1216.7160.0149.2299.6191.81.5622.722.9
    Ta6B4C3584.7539.6614.2203.0189.9279.9195.5168.0144.1305.9213.91.4327.427.3
    Ta7B4C4563.1547.5571.5183.6170.4281.4200.2162.0164.3303.9200.81.5124.424.5
    Ta7B6C3584.7540.0614.2203.0190.0280.0195.5168.0144.1305.9213.91.4327.427.3
    TaB23022001.5124.424.5
    NbB22871951.4724.824.8
    TaC3242151.5125.625.9
    NbC2391611.4821.621.4
    SiC2131871.1433.632.2
    Al2O32321471.5818.718.7
    TiN2591801.4424.324.0
    注: a二元相力学性能数据来自Materials Project网站.
    下载: 导出CSV
  • [1]

    Tian Y J, Xu B, Zhao Z S 2012 Int. J. Refract. Met. Hard Mater. 33 93

    [2]

    包括, 马帅领, 徐春红, 崔田 2017 物理学报 66 036104

    Bao K, Ma S L, Xu C H, Cui T 2017 Acta Phys. Sin. 66 036104

    [3]

    Zhou X F, Sun J, Fan Y X, Chen J, Wang H T, Guo X J, He J L, Tian Y J 2007 Phys. Rev. B 76 100101

    [4]

    Wu Q H, Hu Q K, Hou Y M, Wang H Y, Zhou A G, Wang L B 2018 J. Phys. Condens. Matter 30 385402

    [5]

    Tian Y J, Xu B, Yu D L, Ma Y M, Wang Y B, Jiang Y B, Hu W T, Tang C C, Gao Y F, Luo K, Zhao Z S, Wang L M, Wen B, He J L, Liu Z Y 2013 Nature 493 385

    [6]

    Huang Q, Yu D L, Xu B, Hu W T, Ma Y M, Wang Y B, Zhao Z S, Wen B, He J L, Liu Z Y, Tian Y J 2014 Nature 510 250

    [7]

    徐波, 田永君 2017 物理学报 66 036201

    Xu B, Tian Y J 2017 Acta Phys. Sin. 66 036201

    [8]

    Wu Q H, Hu Q K, Hou Y M, Wang H Y, Zhou A G, Wang L B, Cao G H 2018 Mater. Des. 140 45

    [9]

    Cumberland R W, Weinberger M B, Gilman J J, Clark S M, Tolbert S H, Kaner R B 2005 J. Am. Chem. Soc. 127 7264

    [10]

    Chung H Y, Weinberger M B, Levine J B, Kavner A, Yang J M, Tolbert S H, Kaner R B 2007 Science 316 436

    [11]

    Gregoryanz E, Sanloup C, Somayazulu M, Badro J, Fiquet G, Mao H K, Hemley R J 2004 Nat. Mater. 3 294

    [12]

    Young A F, Sanloup C, Gregoryanz E, Scandolo S, Hemley R J, Mao H K 2006 Phys. Rev. Lett. 96 155501

    [13]

    Ivanovskii A L 2012 Prog. Mater. Sci. 57 184

    [14]

    陶强, 马帅领, 崔田, 朱品文 2017 物理学报 66 036103

    Tao Q, Ma S L, Cui T, Zhu P W 2017 Acta Phys. Sin. 66 036103

    [15]

    Hillebrecht H, Gebhardt K 2001 Angew. Chem. Int. Ed. 40 1445

    [16]

    胡前库, 侯一鸣, 吴庆华, 秦双红, 王李波, 周爱国 2019 物理学报 68 096201

    Hu Q K, Hou Y M, Wu Q H, Qin S H, Wang L B, Zhou A G 2019 Acta Phys. Sin. 68 096201

    [17]

    Wang P F, Weng M Y, Xiao Y, Hu Z X, Li Q H, Li M, Wang Y D, Chen X, Yang X N, Wen Y R, Yin Y X, Yu X Q, Xiao Y G, Zheng J X, Wan L J, Pan F, Guo Y G 2019 Adv. Mater. 31 1903483

    [18]

    Xiao W J, Xin C, Li S B, Jie J S, Gu Y, Zheng J X, Pan F 2018 J. Mater. Chem. A 6 9893

    [19]

    Kresse G, Furthmüller J 1996 Phys. Rev. B 54 11169

    [20]

    Perdew J P, Burke K, Ernzerhof M 1996 Phys. Rev. Lett. 77 3865

    [21]

    Le Page Y, Saxe P 2002 Phys. Rev. B 65 104104

    [22]

    Togo A, Tanaka I 2015 Scr. Mater. 108 1

    [23]

    Togo A, Chaput L, Tanaka I, Hug G 2010 Phys. Rev. B 81 174301

    [24]

    Mouhat F, Coudert F X 2014 Phys. Rev. B 90 224104

    [25]

    Wu Z J, Zhao E J, Xiang H P, Hao X F, Liu X J, Meng J 2007 Phys. Rev. B 76 054115

    [26]

    Pugh S F 1954 Philos. Mag. 45 823

    [27]

    Chen X Q, Niu H Y, Li D Z, Li Y Y 2011 Intermetallics 19 1275

  • [1] 胡前库, 侯一鸣, 吴庆华, 秦双红, 王李波, 周爱国. 过渡金属硼碳化物TM3B3C和TM4B3C2稳定性和性能的理论计算. 物理学报, 2019, 68(9): 096201. doi: 10.7498/aps.68.20190158
    [2] 赵红霞, 赵晖, 陈宇光, 鄢永红. 一维扩展离子Hubbard模型的相图研究. 物理学报, 2015, 64(10): 107101. doi: 10.7498/aps.64.107101
    [3] 汪志刚, 张杨, 文玉华, 朱梓忠. ZnO原子链的结构稳定性和电子性质的第一性原理研究. 物理学报, 2010, 59(3): 2051-2056. doi: 10.7498/aps.59.2051
    [4] 刘子媛(Ziyuan Liu), 潘金波, 张余洋, 杜世萱. 原子尺度构建二维材料的第一性原理计算研究. 物理学报, 2021, (): . doi: 10.7498/aps.70.20201636
    [5] 陈献, 程梅娟, 吴顺情, 朱梓忠. 石墨炔衍生物结构稳定性和电子结构的第一性原理研究. 物理学报, 2017, 66(10): 107102. doi: 10.7498/aps.66.107102
    [6] 黄炳铨, 周铁戈, 吴道雄, 张召富, 李百奎. 空位及氮掺杂二维ZnO单层材料性质:第一性原理计算与分子轨道分析. 物理学报, 2019, 68(24): 246301. doi: 10.7498/aps.68.20191258
    [7] 王 岩, 韩晓艳, 任慧志, 侯国付, 郭群超, 朱 锋, 张德坤, 孙 建, 薛俊明, 赵 颖, 耿新华. 相变域硅薄膜材料的光稳定性. 物理学报, 2006, 55(2): 947-951. doi: 10.7498/aps.55.947
    [8] 宋庆功, 姜恩永, 裴海林, 康建海, 郭 英. 插层化合物LixTiS2中Li离子-空位二维有序结构稳定性的第一性原理研究. 物理学报, 2007, 56(8): 4817-4822. doi: 10.7498/aps.56.4817
    [9] 白静, 王晓书, 俎启睿, 赵骧, 左良. Ni-X-In(X=Mn,Fe和Co)合金的缺陷稳定性和磁性能的第一性原理研究. 物理学报, 2016, 65(9): 096103. doi: 10.7498/aps.65.096103
    [10] 郑路敏, 钟淑英, 徐波, 欧阳楚英. 锂离子电池正极材料Li2MnO3稀土掺杂的第一性原理研究. 物理学报, 2019, 68(13): 138201. doi: 10.7498/aps.68.20190509
    [11] 王丹, 邹娟, 唐黎明. 氢化二维过渡金属硫化物的稳定性和电子特性: 第一性原理研究. 物理学报, 2019, 68(3): 037102. doi: 10.7498/aps.68.20181597
    [12] 张 凯, 冯 俊. 相对论Birkhoff系统的对称性与稳定性. 物理学报, 2005, 54(7): 2985-2989. doi: 10.7498/aps.54.2985
    [13] 王作雷. 一类简化Lang-Kobayashi方程的Hopf分岔及其稳定性. 物理学报, 2008, 57(8): 4771-4776. doi: 10.7498/aps.57.4771
    [14] 陈海军, 李高清, 薛具奎. 变分法研究一维Bose-Fermi系统的稳定性. 物理学报, 2011, 60(4): 040304. doi: 10.7498/aps.60.040304.1
    [15] 闫小童, 侯育花, 郑寿红, 黄有林, 陶小马. Ga, Ge, As掺杂对锂离子电池正极材料Li2CoSiO4的电化学特性和电子结构影响的第一性原理研究. 物理学报, 2019, 68(18): 187101. doi: 10.7498/aps.68.20190503
    [16] 谢 莉, 雷银照. 线性瞬态涡流电磁场定解问题解的唯一性和稳定性. 物理学报, 2006, 55(9): 4397-4406. doi: 10.7498/aps.55.4397
    [17] 时培明, 刘 彬, 刘 爽. 一类谐波激励相对转动非线性动力系统的稳定性与近似解. 物理学报, 2008, 57(8): 4675-4684. doi: 10.7498/aps.57.4675
    [18] 刘爽, 刘彬, 张业宽, 闻岩. 一类时滞非线性相对转动系统的Hopf分岔与周期解的稳定性. 物理学报, 2010, 59(1): 38-43. doi: 10.7498/aps.59.38
    [19] 刘浩然, 朱占龙, 时培明. 一类相对转动时滞非线性动力系统的稳定性分析. 物理学报, 2010, 59(10): 6770-6777. doi: 10.7498/aps.59.6770
    [20] 杨秀峰, 刘谋斌. 光滑粒子动力学SPH方法应力不稳定性的一种改进方案 . 物理学报, 2012, 61(22): 224701. doi: 10.7498/aps.61.224701
  • 引用本文:
    Citation:
计量
  • 文章访问数:  944
  • PDF下载量:  78
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-02-18
  • 修回日期:  2020-04-01
  • 刊出日期:  2020-06-05

三元Nb系和Ta系硼碳化物稳定性和物理性能的第一性原理研究

  • 河南理工大学材料科学与工程学院, 河南省深地材料科学与技术重点实验室, 焦作 454000
  • 通信作者: 周爱国, zhouag@hpu.edu.cn
    基金项目: 国家级-国家自然科学基金(51472075,51772077)

摘要: 过渡金属轻元素化合物是高硬度材料的潜在候选. 以往研究多集中在二元过渡金属硼化物、碳化物和氮化物, 三元相的研究则相对较少. 本文通过提炼和堆垛已知相Nb3B3C (Ta3B3C)和Nb4B3C2 (Ta4B3C2)的结构基元, 构建不同组分的Nb-B-C和Ta-B-C三元相结构模型, 采用第一性原理计算方法, 计算所建结构的形成焓、声子谱和弹性常数, 通过判断其热力学、动力学和力学稳定性, 绘出了三元Nb-B-C和Ta-B-C相图, 成功预测了5种Nb-B-C和6种Ta-B-C三元稳定相. 力学和电学性能计算结果显示Nb-B-C和Ta-B-C三元稳定相均为高硬度导电材料, 硬度大约为25 GPa.

English Abstract

    • 硬质材料拥有高硬度、抗高压、耐摩擦磨损等优异的力学性能, 同时又具有良好的电学、光学和热学性能, 在众多工业领域都发挥着重要的作用[1]. 高硬度材料通常可以分为两类: 第一类材料是轻元素B, C, N, O结合所形成的强共价键化合物, 比如金刚石、立方BN、立方BC2N、立方BC5、B6O、不同成分C-N材料等[2-8]. 这一类材料往往是超硬材料(维氏硬度Hv > 40 GPa), 硬度比较高, 但是合成条件比较苛刻, 一般需要高温高压才能合成. 另一类高硬度材料存在于过渡金属(transition metals, TM)的轻元素化合物, 比如硼化物、碳化物和氮化物. 过渡金属硼化物体系的研究主要集中在TMB2组分, 比如2005年 Hv超过20 GPa的OsB2被成功制备[9], 2007年Chung等[10]采用电弧方法成功合成出了RuB2和ReB2. 过渡金属碳化物体系的研究较多集中在TiC, ZrC, HfC, VC, NbC, TaC等TMC组分, 这些化合物往往通过常压合成方法均可制备得到. 在过渡金属氮化物体系中, OsN2, IrN2和PtN2化合物采用对顶砧高压技术被成功合成[11,12], 这些材料均为极难压缩材料, 但合成条件比较苛刻, 合成所需压力很高. 关于过渡金属硼化物、碳化物和氮化物的研究现状, 可参阅综述文献[2,13,14].

      过渡金属轻元素化合物的研究主要集中于二元过渡金属硼化物、碳化物和氮化物, 三元相的研究则相对较少. 由于硼、碳、氮三种元素在元素周期表中毗邻, 三种元素原子大小差不多, 其核外电子排布也相近, 因此可推理, 三元过渡金属轻元素化合物的存在是完全有可能的. 2001年德国科学家Hillebrecht和Gebhardt[15]以Nb, B和C三种单质作为反应原料, 以Al-Cu金属合金作为助熔剂, 成功制备了Nb3B3C和Nb4B3C2晶体, 并给出了其晶体结构. 考虑到类质同晶效应, 我们以Nb3B3C和Nb4B3C2晶体结构为模板, 用28种其他过渡金属元素分别替代两种结构中Nb原子位置, 理论计算了这些新构建化合物的结构稳定性, 找到了Ta3B3C和Ta4B3C2两种稳定新相[16]. 由于Hillebrecht和Gebhardt所合成的晶体尺寸较小, 所以实验未能测得其性能[15]. 采用第一性原理计算了实验合成的Nb3B3C和Nb4B3C2相、以及理论预测稳定的Ta3B3C和Ta4B3C2相的力学和电学性能, 显示这4种化合物都是导电的高硬度材料[16].

      根据材料基因组学, 晶体材料最基本的结构基元是中心原子与周围原子所组成的配位多面体. 结构基元按照一定的空间排列方式进行周期排列, 从而构成晶体. 结构基元内的成键方式决定了晶体材料的本征物理化学性质, 因此从这个角度说, 结构基元就是材料的基因. 例如, Wang等[17]和Xiao等[18]以“无机类苯环”或“无机类苯环芳香性”作为功能结构基元, 设计了一系列锂电池正极材料.

      本文通过分析Nb3B3C和Nb4B3C2 (Ta3B3C和Ta4B3C2)结构, 得出其结构基元, 通过裁剪和堆垛结构基元, 构建不同组分的Nb-B-C和Ta-B-C结构模型, 采用第一性原理计算方法, 计算所建结构的形成焓、声子谱和弹性常数, 判断其热力学、动力学和力学稳定性, 基于所得Nb-B-C和Ta-B-C三元相图, 得出稳定的Nb-B-C和Ta-B-C化合物, 并研究稳定相的力学和电学性能.

    • 计算采用基于密度泛函理论的第一性原理平面波赝势方法, 由VASP软件包[19]完成. 结构优化过程中, 采用全优化方式, 即原子位置和晶格常数未做任何限制. 电子间交换关联能采用广义梯度近似下的PBE泛函. 平面波截断能Ecut设置为600 eV. 赝势采用全电子投影缀加平面波方法[20]. Nb原子的价电子为4p64d45s1, Ta原子的价电子为5p65d36s2, B原子的价电子为2s22p1, C原子的价电子为2s22p2. 第一布里渊区积分采用Monkhorst-Pack形式的特殊K点方法, K点精度设置为2π × 0.03 Å–1. 优化过程中能量迭代收敛标准为1 × 10–5 eV/atom. 弹性常数的计算采用应力应变方法[21]. 声子谱的计算采用有限位移法, 由PHONOPY软件[22]计算完成.

    • 由于Nb3B3C和Ta3B3C, Nb4B3C2和Ta4B3C2结构相同, 所以在以下讨论中以Ta3B3C和Ta4B3C2为例进行说明. 图1给出了Ta3B3C和Ta4B3C2的晶体结构. 可以看出Ta3B3C和Ta4B3C2都具有层状的堆垛方式, 堆垛的单元为Ta-C区和Ta-B区. 在Ta-C区中, C原子(Ta原子)处于周围6个Ta原子(C原子)所形成的八面体的中心, 组成岩盐矿型晶体结构, 所以Ta原子和C原子的配位数均为6. 在Ta-B区中, Ta原子和B原子形成AlB2型结构, B原子位于6个Ta原子所形成的三棱柱中心. 每个B原子除了与这6个Ta原子成键以外, 还与相邻的3个B原子成键, 因此其配位数为9. 通过分析结构, 可以得知Ta6C八面体和Ta6B三棱柱就是这两个结构的结构基元. Ta3B3C和Ta4B3C2两相的区别就在于具有不同的结构基元层数.

      图  1  (a), (b) Ta3B3C; (c) Ta4B3C2; (d) Ta3BC2; (e) Ta6B4C3; (f) Ta7B4C4; (g) Ta7B6C3的晶体结构. 棕球: Ta原子; 蓝球: B原子; 粉球: C原子. Ta6B三棱柱和Ta6C八面体分别用绿色和褐色表示

      Figure 1.  The crystal structures of (a), (b) Ta3B3C; (c) Ta4B3C2; (d) Ta3BC2; (e) Ta6B4C3; (f) Ta7B4C4; (g) Ta7B6C3. The light brown, blue and pink spheres represent Ta, B, and C atoms, respectively. The Ta6B triangular prisms and Ta6C octahedrons are painted green and dark brown.

      通过调整Ta6C八面体和Ta6B三棱柱的数量, 可以构造一系列以Ta6C八面体和Ta6B三棱柱为结构基元的三元Ta-B-C化合物. 通过分析图1中所示Ta3B3C和Ta4B3C2结构, 可以得出该系列Ta-B-C化合物结构式为TaB(TaB2)mTaB(TaC)n, 化学式为Ta(m + n + 2)B(2m + 2)Cn. 对于Ta3B3C相来说, m = 2, n = 2; 对于Ta4B3C2相来说, m = 2, n = 4. 在本文中限定m取值范围为0 ≤ m ≤ 4, n取值范围为1 ≤ n ≤ 4. 因此通过mn数值的排列组合, 可以分别构建20种三元Nb-B-C化合物和20种三元Ta-B-C化合物, 如表1中所列. 部分Ta-B-C化合物晶体结构见图1中所示.

      mn空间群模型晶格参数/Å模型晶格参数/Å
      abcabc
      01CmmmNb3B2C3.25413.8083.141Ta3B2C3.24013.6973.127
      02CmcmNb2BC3.23518.3303.153Ta2BC3.22018.1653.140
      03CmmmNb5B2C33.22522.9033.153Ta5B2C33.19922.6593.138
      04CmcmNb3BC23.21427.3763.156Ta3BC23.19827.1323.150
      11PmmmNb4B4C3.29018.9943.145Ta4B4C3.27718.8783.127
      12ImmmNb5B4C23.26723.6003.150Ta5B4C23.24823.3773.138
      13PmmmNb6B4C33.24328.0283.154Ta6B4C33.22527.8723.141
      14ImmmNb7B4C43.24232.5453.158Ta7B4C43.22432.3153.147
      21CmmmNb5B6C3.30224.4143.134Ta5B6C3.28924.2083.122
      22CmcmNb3B3C3.28428.8773.144Ta3B3C3.26728.6883.133
      23CmmmNb7B6C33.26433.3643.148Ta7B6C33.24633.1643.136
      24CmcmNb4B3C23.25737.8743.153Ta4B3C23.24337.6093.141
      31PmmmNb6B8C3.30914.8893.137Ta6B8C3.29814.7883.122
      32ImmmNb7B8C23.29034.2473.144Ta7B8C23.27634.0073.131
      33PmmmNb8B8C33.27619.3503.148Ta8B8C33.25819.2353.135
      34ImmmNb9B8C43.26843.2553.151Ta9B8C43.25242.9773.138
      41CmmmNb7B10C3.31235.1923.131Ta7B10C3.29934.9903.116
      42CmcmNb4B5C3.29639.6943.139Ta4B5C3.28039.4413.125
      43CmmmNb9B10C33.28144.2063.142Ta9B10C33.26343.9243.130
      44CmcmNb5B5C23.27348.7293.145Ta5B5C23.25748.4003.134

      表 1  不同成分Nb(m + n + 2)B(2m + 2)Cn和Ta(m + n + 2)B(2m + 2)Cn晶体的结构参数

      Table 1.  Structural parameters of Nb(m + n + 2)B(2m + 2)Cn and Ta(m + n + 2)B(2m + 2)Cn crystals.

      表1给出了构建的20种三元Nb-B-C化合物和20种三元Ta-B-C化合物的空间群及其结构优化后的晶格参数. 全部结构均为正交晶系. 当m为偶数、n为奇数时, 对称性为Cmmm; m为偶数、n为偶数时, 对称性为Cmcm; m为奇数、n为奇数时, 对称性为Pmmm; m为奇数、n为偶数时, 对称性为Immm. 各个结构在ac轴方向排布相同, 只是沿b轴层数不同, 所以各个结构的晶格参数ac大小很接近.

    • 鉴于Nb3B3C, Nb4B3C2, Ta3B3C和Ta4B3C2皆为稳定相[16], 因此构造的其他组分三元Nb-B-C和Ta-B-C体系中也可能存在稳定结构. 一个热力学稳定的三元化合物, 意味着其不能分解成单质、二元相以及其他三元相的任意组合, 也就是其能量要小于相应单质、二元相以及其他三元相的任意组合. 首先判断所构建的三元Nb-B-C化合物和三元Ta-B-C化合物相对于单质的稳定性, 计算公式如下:

      $\begin{split} \Delta {H_{{\rm{elements}}}} =\; & [(x + y + z) {H_{{\rm{TM}} _x{\rm{B}}_y{\rm{C}}_z}} - x {H_{{\rm{TM}}}} \\ & - y {H_{\rm{B}}} - z {H_{\rm{C}}}]/(x + y + z).\end{split}$

      式中$ \Delta {H_{{\rm{elements}}}} $表示TMxByCz (TM = Ta或Nb)相对其三种单质的形成焓(eV/atom), ${H}_{{\rm{T}}{\rm{M}}_x{\rm{B}}_y{\rm{C}}_z}$, $ {H}_{{\rm{T}}{\rm{M}}} $, $ {H}_{{\rm{B}}} $$ {H}_{{\rm{C}}} $分别是三元TMxByCz相、单质TM、单质硼以及石墨的焓值(eV/atom). 如果$ \Delta{H}_{\rm{elements}} > 0$, 则意味着三元相是不稳定的; 否则就说明三元相相对其三种单质材料来说是稳定的. 三元Nb-B-C化合物和三元Ta-B-C化合物相对单质的形成焓列于表2中. 可见所有化合物相对单质的形成焓均为负值, 这说明相对单质, 所构建的Nb-B-C化合物和Ta-B-C化合物都是稳定的.

      Phases$ \Delta{H}_{\rm{elements}} $$ \Delta{H}_{\rm{comp}} $最稳定竞争组合Phases$ \Delta{H}_{\rm{elements}} $$ \Delta{H}_{\rm{comp}} $最稳定竞争组合
      Nb3B2C–0.6200.070Nb3B4 + 6NbB + Nb6C5 = 5Nb3B2CTa3B2C–0.6510.086Ta3BC2 + 3TaB = 2Ta3B2C
      Nb2BC–0.6190.029Nb3B4 + NbB + Nb6C5 = 5Nb2BCTa2BC–0.6640.035Ta3BC2 + TaB = 2Ta2BC
      Nb5B2C3–0.5860.0363Nb3B4 + Nb7B4C4 + 4Nb6C5 = 8Nb5B2C3Ta5B2C3–0.6550.0213Ta3BC2 + TaB = 2Ta5B2C3
      Nb3BC2–0.5860.019Nb3B4 + 3Nb7B4C4 + 4Nb6C5 = 16Nb3BC2Ta3BC2–0.660–0.002TaB + 2TaC = Ta3BC2
      Nb4B4C–0.6790.0303Nb3B4 + Nb7B4C4 = 4Nb4B4CTa4B4C–0.6910.044Ta7B4C4 + 3Ta3B4 = 4Ta4B4C
      Nb5B4C2–0.6680.006Nb3B4 + Nb7B4C4 = 2Nb5B4C2Ta5B4C2–0.6940.019Ta7B4C4 + Ta3B4 = 2Ta5B4C2
      Nb6B4C3–0.6450.005Nb3B4 + 3Nb7B4C4 = 4Nb6B4C3Ta6B4C3–0.6930.0043Ta7B4C4 + Ta3B4 = 4Ta6B4C3
      Nb7B4C4–0.632–0.0063Nb3B4 + 2C + 2Nb6C5 = 3Nb7B4C4Ta7B4C4–0.685–0.0173Ta3B4 + 4TaC = Ta7B4C4
      Nb5B6C–0.6970.0153Nb3B4 + C + 2Nb3B3C = 3Nb5B6CTa5B6C–0.6970.024C + Ta5B6 = Ta5B6C
      Nb3B3C–0.685–0.0013Nb3B4 + C + 3Nb4B3C2 = 7Nb3B3CTa3B3C–0.6990.0103Ta7B4C4 + 9Ta3B4 + 4C = 16Ta3B3C
      Nb7B6C3–0.6640.0005Nb3B3C + Nb4B3C2 = Nb7B6C3Ta7B6C3–0.6950.00085Ta7B4C4 + 7Ta3B4 + 4C = 8Ta7B6C3
      Nb4B3C2–0.648–0.0015Nb3B4 + 4C + 7Nb7B4C4 = 16Nb4B3C2Ta4B3C2–0.6840.0027Ta7B4C4 + 5Ta3B4 + 4C = 16Ta4B3C2
      Nb6B8C–0.6950.0192Nb3B4 + C = Nb6B8CTa6B8C–0.6850.0342Ta3B4 + C = Ta6B8C
      Nb7B8C2–0.6830.0083Nb3B4 + 2C + 4Nb3B3C = 3Nb7B8C2Ta7B8C2–0.6860.020Ta7B4C4 + 7Ta3B4 + 4C = 4Ta7B8C2
      Nb8B8C3–0.6650.008C + 8Nb3B3C = 3Nb8B8C3Ta8B8C3–0.6840.012Ta7B4C4 + 3Ta3B4 + 2C = 2Ta8B8C3
      Nb9B8C4–0.6510.008C + 5Nb3B3C + 3Nb4B3C2 = 3Nb9B8C4Ta9B8C4–0.6750.0133Ta7B4C4 + 5Ta3B4 + 4C = 4Ta9B8C4
      Nb7B10C–0.6930.021C + 2Nb2B3 + Nb3B4 = Nb7B10CTa7B10C–0.6770.030TaB2 + 2Ta3B4 + C = Ta7B10C
      Nb4B5C–0.6840.0112C + Nb3B3C + 3Nb3B4 = 3Nb4B5CTa4B5C–0.6790.026Ta7B4C4 + 19Ta3B4 + 12C = 16Ta4B5C
      Nb9B10C3–0.6680.012C + 2Nb3B3C + Nb3B4 = Nb9B10C3Ta9B10C3–0.6770.0193Ta7B4C4 + 17Ta3B4 + 12C = 8Ta9B10C3
      Nb5B5C2–0.6550.011C + 5Nb3B3C = 3Nb5B5C2Ta5B5C2–0.6700.0185Ta7B4C4 + 15Ta3B4 + 12C = 16Ta5B5C2

      表 2  不同成分Nb-B-C相和Ta-B-C相的形成焓 (单位: eV/atom), $ \Delta{H}_{\rm{elements}} $表示单质为反应物, $ \Delta{H}_{\rm{comp}} $表示最稳定竞争组合为反应物

      Table 2.  Calculated formation enthalpies of different Nb-B-C and Ta-B-C phases (in eV/atom).$ \Delta{H}_{\rm{elements}} $ represents the elements as the reactants, and $\Delta{H}_{\rm{comp}}$ indicates the most stable composite as the reactants.

      但是除了单质, 三元相还有可能分解成单质、二元相和其他三元相的任意组合. 在所有各种可能任意组合中, 焓值总和最低的那一组合定义为该三元相的最稳定竞争组合. 对Nb-B-C和Ta-B-C来说, 通过在ICSD数据库搜索其各种单质、二元相和三元相结构, 计算其焓值. 三元相相对最稳定竞争组合的热力学稳定性可用下列公式进行计算:

      $ \Delta{H}_{\rm{comp}}={H}_{{\rm{T}}{\rm{M}}x{\rm{B}}y{\rm{C}}z}{-H}_{\rm{comp}}, $

      式中$ \Delta{H}_{\rm{comp}} $是TMxByCz (TM = Ta或Nb)相对最稳定竞争组合的形成焓(eV/atom), $ {H}_{{\rm{T}}{\rm{M}}x{\rm{B}}y{\rm{C}}z} $$ {H}_{\rm{comp}} $分别是三元相TMxByCz和最稳定竞争组合的焓值(eV/atom). 如果$ \Delta{H}_{\rm{comp}} > 0$, 则意味着三元相TMxByCz会分解成最稳定竞争组合; 否则就说明相对所有可能分解产物, 三元TMxByCz化合物都能稳定存在.

      基于Nb-B-C和Ta-B-C体系中二元相和三元相相对单质的形成焓, 分别绘出了Nb-B-C和Ta-B-C的三元成分相图, 如图2所示. 表2列出了三元Nb-B-C相和三元Ta-B-C相的形成焓及其最稳定竞争组合. 从图2表2中, 可以清楚地看出各相的稳定性. 相对最稳定竞争组合来说, Nb系三元相中有Nb3B3C, Nb4B3C2和Nb7B4C4是稳定的, 这与这三相已被成功合成的实验结果[15]是一致的, 这也证明了计算结果的可靠性. Ta系中有Ta3BC2和Ta7B4C4两相是稳定的. 之前的计算文献[16]报道Ta3B3C和Ta4B3C2两相是稳定的, 是因为没有考虑Ta3BC2和Ta7B4C4的原因, 这一点从Ta3B3C和Ta4B3C2最稳定竞争组合中均包括Ta7B4C4相可看出.

      图  2  (a) Nb-B-C和(b) Ta-B-C三元相图. 红色, 稳定相; 蓝色, 亚稳相; 绿色, 不稳定相

      Figure 2.  Ternary phase diagrams of (a) Nb-B-C and (b) Ta-B-C. Red, stable; blue, metastable; green, unstable.

      除了上述4种稳定相以外, 表2中还有一些化合物的形成焓虽为正值, 但却很接近零值. 把$0 < \Delta {H}_{\rm{comp}} \leqslant 0.005\;{\rm{eV}}/{\rm{atom}}$的化合物定义为亚稳相. 基于此标准, Nb-B-C体系中的Nb6B4C3, Nb7B6C3和Ta-B-C体系中的Ta6B4C3, Ta7B6C3, Ta3B3C和Ta4B3C2皆为亚稳相. 第一性原理计算所得能量为绝对零度下物质的基态能量, 因此如果加以考虑温度作用, 亚稳相的形成能是有可能变成负值的. 采用准谐近似方法[22,23]计算了Nb-B-C和Ta-B-C体系中三元稳定相和三元亚稳相及其最稳定竞争组合相在0—2000 K下的自由能, 计算结果如图3所示. 在Nb-B-C体系中, 考虑零点振动能之后, 亚稳相Nb7B6C3变成了稳定相, 这与之前实验合成得到该相的结果[15]相一致. 在1730 K, 亚稳相Nb6B4C3形成能由正转负, 意味着该相变成了稳定相. 在Ta-B-C体系中, 考虑零点振动能之后, 亚稳相Ta7B6C3变成了稳定相. 亚稳相Ta6B4C3, Ta4B3C2和Ta3B3C分别在210, 360和1100 K转变成了稳定相. 但同时注意到, 基态稳定相Ta3BC2在130 K转变成了不稳定相. 综上所述, 在温度高于室温条件下, 在Nb-B-C和Ta-B-C体系中存在Nb3B3C, Nb4B3C2, Nb6B4C3, Nb7B4C4, Nb7B6C3和Ta3B3C, Ta4B3C2, Ta6B4C3, Ta7B4C4, Ta7B6C3 10种稳定相. 除此之外, 还预测了一种基态稳定相Ta3BC2. 目前除了Nb3B3C, Nb4B3C2, Nb7B4C4, Nb7B6C3已经在实验上合成以外, 我们期待其他新相在未来可以实验合成.

      图  3  不同温度下 (a) Nb-B-C和(b) Ta-B-C三元相分别和其相应最稳定竞争组合相的自由能之差

      Figure 3.  Energy differences of (a) Nb-B-C and (b) Ta-B-C ternary phases with respect to their most competing phases as a function of temperature.

    • 通过计算声子谱可以判断一个结构的动力学稳定性. 图4给出了本文预测的Nb-B-C和Ta-B-C体系中三元稳定相的声子谱. 可以看出所有结构都没有声子虚频存在, 从而表明这些结构在动力学上是稳定的.

      图  4  Nb-B-C和Ta-B-C三元相的声子色散曲线

      Figure 4.  Phonon dispersion curves of Nb-B-C and Ta-B-C ternary phases.

      材料结构的力学稳定性可以通过检验其弹性常数是否满足“波恩-黄”稳定判据[24]来判断. 所预测的热力学稳定相都属于正交晶系. 对于正交晶系晶体来说, 其稳定结构的9个独立弹性常数需满足如下条件[24]: $ {C}_{11}>0, {C}_{44}>0, {C}_{55}>0, {C}_{66}>0, $ ${C}_{11}{C}_{22}>{C}_{12}^{2}, {C}_{11}{C}_{22}{C}_{33}+2{C}_{12}{C}_{13}{C}_{23}-{C}_{11}{C}_{23}^{2}- $ ${C}_{22}{C}_{13}^{2}-{C}_{33}{C}_{12}^{2}>0 $. 预测的Nb-B-C和Ta-B-C体系中三元稳定相的弹性常数列于表3中. 经过计算, 表3中所列结构的弹性常数都满足上述判据, 说明这些相都是力学上的稳定相.

      结构弹性常数力学性能a硬度
      C11C22C33C44C55C66C12C13C23BGB/GHChenHTian
      Nb3B3C544.3479.8522.8181.5171.9245.3170.9132.9162.2275.3189.71.4524.824.7
      Nb4B3C2551.5499.2548.5184.0175.1257.1183.2132.7157.8282.9195.81.4425.525.4
      Nb6B4C3533.3493.8548.1174.9161.3255.2175.4138.9151.7278.5189.51.4724.424.3
      Nb7B4C4535.9505.9526.4172.2161.3259.1184.0142.8152.6280.6188.31.4923.923.8
      Nb7B6C3553.1494.5563.2188.7179.6255.6176.4132.1157.7282.5198.91.4226.326.2
      Ta3B3C569.6514.4563.5194.1180.0261.8187.1147.3173.9295.9200.81.4725.325.3
      Ta4B3C2581.1535.3602.1197.3185.1275.8200.3146.0170.2305.7209.01.4626.226.2
      Ta3BC2550.0547.7550.0159.8159.5292.1216.7160.0149.2299.6191.81.5622.722.9
      Ta6B4C3584.7539.6614.2203.0189.9279.9195.5168.0144.1305.9213.91.4327.427.3
      Ta7B4C4563.1547.5571.5183.6170.4281.4200.2162.0164.3303.9200.81.5124.424.5
      Ta7B6C3584.7540.0614.2203.0190.0280.0195.5168.0144.1305.9213.91.4327.427.3
      TaB23022001.5124.424.5
      NbB22871951.4724.824.8
      TaC3242151.5125.625.9
      NbC2391611.4821.621.4
      SiC2131871.1433.632.2
      Al2O32321471.5818.718.7
      TiN2591801.4424.324.0
      注: a二元相力学性能数据来自Materials Project网站.

      表 3  Nb-B-C和Ta-B-C三元相的弹性常数Cij、体模量B、剪切模量 G和维氏硬度Hv (单位: GPa)

      Table 3.  Elastic constants Cij, bulk modulus B, shear modulus G, Vickers hardness Hv of Nb-B-C and Ta-B-C ternary phases (in GPa).

    • 多晶材料的体模量B和剪切模量G可以通过以下公式[25]计算:

      $ B = 1/2\left( {{B_{\rm{V}}} + {B_{\rm{R}}}} \right),G = 1/2\left( {{G_{\rm{V}}} + {G_{\rm{R}}}} \right), $

      其中下标V和R分别代表Voigt和Reuss方法.

      正交晶系晶体的BV, BR, GVGR可由弹性常数通过以下公式计算得出:

      $\begin{split} {B_{\rm{V}}} =\;& 1/9\left[ {{C_{11}} + {C_{22}} + {C_{33}} + 2\left( {{C_{12}} + {C_{13}} + {C_{23}}} \right)} \right],\\ {G_{\rm{V}}} =\;& 1/15[{C_{11}} + {C_{22}} + {C_{33}} + 3({C_{44}} + {C_{55}} + {C_{66}}) - ({C_{12}} + {C_{13}} + {C_{23}})],\\ {B_{\rm{R}}} =\;& [{C_{11}}\left( {{C_{22}} + {C_{33}} - 2{C_{23}}} \right) + {C_{22}}\left( {{C_{33}} - 2{C_{13}}} \right) - 2{C_{33}}{C_{12}} + {C_{12}}\left( {2{C_{23}} - {C_{12}}} \right) + {C_{13}}\left( {2{C_{12}} - {C_{13}}} \right) \\ & +{C_{23}}\left( {2{C_{13}} - {C_{23}}} \right){]^{ - 1}} \times [{C_{13}}\left( {{C_{12}}{C_{23}} - {C_{13}}{C_{22}}} \right) + {C_{23}}\left( {{C_{12}}{C_{13}} - {C_{23}}{C_{11}}} \right) + {C_{33}}\left( {{C_{11}}{C_{22}} - {C^2}_{12}} \right)],\\ {B_{\rm{R}}} =\;& [{C_{11}}\left( {{C_{22}} + {C_{33}} - 2{C_{23}}} \right) + {C_{22}}\left( {{C_{33}} - 2{C_{13}}} \right) - 2{C_{33}}{C_{12}} + {C_{12}}\left( {2{C_{23}} - {C_{12}}} \right) + {C_{13}}\left( {2{C_{12}} - {C_{13}}} \right) \\ & +{C_{23}}\left( {2{C_{13}} - {C_{23}}} \right){]^{ - 1}} \times \varDelta ,\\ {G_{\rm{R}}} =\;& 15\{ 4[{C_{11}}\left( {{C_{22}} + {C_{33}} + {C_{23}}} \right) + {C_{22}}\left( {{C_{33}} + {C_{13}}} \right) + {C_{33}}{C_{12}} - {C_{12}}\left( {{C_{23}} + {C_{12}}} \right) - {C_{13}}\left( {{C_{12}} + {C_{13}}} \right) \\ & -{C_{23}}\left( {{C_{13}} + {C_{23}}} \right)]/\varDelta + 3[\left( {1/{C_{44}}} \right) + \left( {1/{C_{55}}} \right) + \left( {1/{C_{66}}} \right)]{\} ^{ - 1}},\\[-10pt] \end{split} $

      其中Δ = C13(C12C23C13C22) + C23(C12C13C23C11) + C33(C11C22C122).

      计算的Nb-B-C和Ta-B-C体系中三元稳定相的体模量B和剪切模量G列于表3中. 所得三元相体模量B位于275—306 GPa的狭小数值范围内, 剪切模量G位于190—214 GPa的狭小数值范围内. 这说明无论是Nb-B-C三元相, 还是Ta-B-C三元相, 它们的力学性能很相似, 这与它们具有相同的结构基元有关. Ta-B-C化合物的BG值均稍大于同一成分的Nb-B-C相, 这是因为Ta原子半径稍小于Nb原子, Ta—B和Ta—C键长要比Nb—B和Nb—C键长稍短. 这一点也可以从表1中所列两体系的晶格常数看出来.

      体模量和剪切模量的比值可以用来判断一个材料是脆性还是韧性. 根据Pugh经验判据[26], 如果B/G > 1.75, 则该材料属于韧性; 反之则为脆性. 根据表3B/G比值, 可知Nb-B-C和Ta-B-C体系中所有三元相均为脆性材料.

      硬度是本文所计算材料的重要性能. 硬度计算采用陈星秋公式[27]和田永君公式[1]:

      $ {H}_{\rm{Chen}}=2{\left({k}^{2}G\right)}^{0.585}-3, $

      $ {H}_{\rm{Tian}}=0.92{k}^{1.137}{G}^{0.708}. $

      (5)式和(6)式中, HChenHTian是维氏硬度, kG/B的比值.

      Nb-B-C和Ta-B-C体系中三元稳定相的硬度数值列于表3中. 可看出, 两种计算公式所得硬度数值很接近, 误差小于0.2 GPa. Ta-B-C三元相的硬度稍高于Nb-B-C三元相. 所有三元相的硬度值很接近, 位于23.8—27.4 GPa范围内, 均小于40 GPa, 这说明均不是超硬材料, 但也属于高硬度材料. 这些三元相和已知常见的二元相硬质材料硬度相当, 因此未来可被作为硬质材料来使用.

      为了探索Nb-B-C和Ta-B-C体系中三元稳定相的电学性能, 计算了其态密度, 结果如图5所示. 由于结构的相似性, 各个结构的态密度也很相似. 为了表示更加清楚, 把Ta3BC2相的态密度图放大加以显示. 所示态密度图中, 虚线所示费米面处没有带隙存在, 电子态不为零, 这说明这些三元相均为导体. 通过图中所示不同电子的部分态密度可知, 费米面处的电子态主要来源于d电子, 这说明这些三元相的导电性主要来源于结构中过渡金属原子的d电子, 这和之前文献[16]中对Ta3B3C和Ta4B3C2态密度的分析结果是一致的.

      图  5  Nb-B-C和Ta-B-C三元相的态密度图

      Figure 5.  Density of states of Nb-B-C and Ta-B-C ternary phases.

    • 本文通过分析实验合成的Nb3B3C和Nb4B3C2结构, 得出其结构基元为Nb6C八面体和Nb6B三棱柱. 通过调整Nb6C (Ta6C)八面体和Nb6B (Ta6B)三棱柱的数量, 构建了20种三元Nb(m + n + 2)B(2m + 2)Cn化合物和20种三元Ta(m + n + 2)B(2m + 2)Cn化合物, 采用基于密度泛函理论第一性原理计算方法, 通过分析这些化合物的热力学、动力学和力学稳定性, 成功预测了Nb3B3C, Nb4B3C2, Nb6B4C3, Nb7B4C4, Nb7B6C3, Ta3B3C, Ta4B3C2, Ta6B4C3, Ta7B4C4, Ta7B6C3, Ta3BC2等11种三元稳定相. 其中Nb6B4C3, Ta6B4C3, Ta4B3C2和Ta3B3C分别在1730, 210, 360和1100 K温度以上才能稳定存在. 而Ta3BC2相只能在130 K温度以下才能稳定存在. 这11种稳定相皆为导电的脆性硬质材料, 其硬度值位于23.8—27.4 GPa范围内. 目前除了Nb3B3C, Nb4B3C2, Nb7B4C4, Nb7B6C3已经在实验上合成以外, 我们期待其他7种新相在未来的实验合成.

参考文献 (27)

目录

    /

    返回文章
    返回