高压下MgN8晶体结构理论模拟与物性研究

缪宇 刘思远 马雪姣 金哲学

引用本文:
Citation:

高压下MgN8晶体结构理论模拟与物性研究

    作者简介: 缪 宇(1999-),男,本科,主要从事高压结构计算研究. E-mail: 15053219411@163.com;
    通讯作者: 金哲学, jinzhexue@ybu.edu.cn
  • 中图分类号: O521.2

Theoretical Simulation and Physical Properties of MgN8 Crystal Structure under High Pressure

    Corresponding author: JIN Zhexue, jinzhexue@ybu.edu.cn
  • CLC number: O521.2

  • 摘要: 基于密度泛函理论第一性原理的方法,使用CALYPSO结构搜索技术结合VASP软件,在0~100 GPa压强范围内对MgN8的晶体结构进行预测,并对预测的结构进行系统研究。结果表明:在常压下,空间群为P4/mbmα-MgN8晶体结构的焓值最低;当压强达到24.3 GPa和68.3 GPa时发生相变,分别相变成空间群为P4/mncβ-MgN8相和空间群为Cmcmγ-MgN8相,两次相变均为对应体积坍塌的一级相变。电子性质计算结果表明,α-MgN8相的导带与价带之间具有3.09 eV的带隙,表明该结构具有非金属性;β相和γ相具有明显的金属特征。Bader电荷转移计算表明,随着压力的增加,Mg原子向N原子转移的电荷逐渐增多。
  • 图 1  MgN8的焓差曲线以及α-MgN8β-MgN8γ-MgN8相体积随压强的变化关系

    Figure 1.  Calculated enthalpies per formula unit of pressure of MgN8 and the calculated volume versus pressure of α-MgN8, β-MgN8 and γ-MgN8

    图 2  α-MgN8β-MgN8γ-MgN8相的晶体结构

    Figure 2.  Crystal structures of α-MgN8, β-MgN8 and γ-MgN8

    图 3  α-MgN8β-MgN8γ-MgN8相声子色散关系和声子态密度

    Figure 3.  Phonon-dispersion curves and the phonon density of states of α-MgN8, β-MgN8 and γ-MgN8

    图 4  α-MgN8β-MgN8γ-MgN8相的能带结构和电子态密度

    Figure 4.  Band structure and partial density of states of α-MgN8, β-MgN8 and γ-MgN8

    图 5  α-MgN8β-MgN8γ-MgN8相的电子局域函数

    Figure 5.  Electron localization function of α-MgN8, β-MgN8 and γ-MgN8

    表 1  α-MgN8β-MgN8γ-MgN8相的平衡态晶格常数

    Table 1.  Lattice parameters of α-MgN8, β-MgN8 and γ-MgN8 in the equilibrium state

    PhasePressure/GPaabcα/(°)β/(°)γ/(°)
    α-MgN8 (P4/mbm)0 5.9135.9136.57290.090.090.0
    β-MgN8 (P4/mnc)24.36.2196.2193.72490.090.090.0
    γ-MgN8 (Cmcm)68.34.1674.1678.68090.090.051.1
    下载: 导出CSV

    表 2  α-MgN8β-MgN8γ-MgN8相的平衡态晶格常数和原子位置

    Table 2.  Lattice parameters and atomic coordinate of α-MgN8, β-MgN8 and γ-MgN8 in the equilibrium state

    PhasePressure/GPaAtomsWyckoff positionxyz
    α-MgN8 (P4/mbm)0 Mg 2b0.0000.0000.500
    N1 2a0.0000.0000.000
    N2 2c0.0000.5000.500
    N3 4e0.0000.0000.821
    N4 4f0.0000.5000.915
    N5 4h0.8590.3590.500
    β-MgN8 (P4/mnc)24.3Mg 2a0.0000.0001.000
    N16i0.4160.2190.823
    γ-MgN8 (Cmcm)68.3Mg 8g0.624–0.376 1.250
    N116h0.226–0.135 0.617
    N216h1.1150.4711.386
    下载: 导出CSV

    表 3  α-MgN8β-MgN8γ-MgN8相电荷转移

    Table 3.  Calculated Bader charges of α-MgN8, β-MgN8 and γ-MgN8

    PhasePressure/GPaAtomsNumberCharge value/eCharge transfer/e
    α-MgN8 (P4/mbm)0Mg16.26 1.74
    N114.84 0.16
    N214.86 0.13
    N325.04–0.04
    N425.46–0.46
    N525.52–0.52
    β-MgN8 (P4/mnc)24.3Mg16.06 1.94
    N85.24–0.24
    γ-MgN8 (Cmcm)68.3Mg15.80 2.20
    N145.28–0.28
    N245.27–0.27
    下载: 导出CSV
  • [1] GROCHALA W, HOFFMANN R, FENG J, et al. The chemical imagination at work in very tight places [J]. Angewandte Chemie International Edition, 2007, 46(20): 3620–3642. doi: 10.1002/anie.200602485
    [2] HEMLEY R J. Effects of high pressure on molecules [J]. Annual Review of Physical Chemistry, 2000, 51(1): 763–800. doi: 10.1146/annurev.physchem.51.1.763
    [3] SCHETTINO V, BINI R. Molecules under extreme conditions: chemical reactions at high pressure [J]. Physical Chemistry Chemical Physics, 2003, 5(10): 1951–1965. doi: 10.1039/b301381b
    [4] 李鑫, 马雪姣, 高文泉, 等. 高压下Ir2P晶体结构预测与物理性质 [J]. 高压物理学报, 2019, 33(1): 011103. doi: 10.11858/gywlxb.20180645
    LI X, MA X J, GAO W Q, et al. Evolution of crystal structures and electronic properties for Ir2P under high pressure [J]. Chinese Journal of High Pressure Physics, 2019, 33(1): 011103. doi: 10.11858/gywlxb.20180645
    [5] MA X J, LI X, ZHOU D, et al. Phase diagram and bonding states of Ir-P binary compounds at high pressures [J]. Journal of Alloys and Compounds, 2019, 791: 1257–1262. doi: 10.1016/j.jallcom.2019.03.051
    [6] BRULS R J, HINTZEN H T, METSELAAR R. Preparation and characterisation of MgSiN2 powders [J]. Journal of Materials Science, 1999, 34(18): 4519–4531. doi: 10.1023/A:1004645407523
    [7] PARKIN I P, NARTOWSKI A M. Solid state metathesis routes to Group IIIa nitrides: comparison of Li3N, NaN3, Ca3N2 and Mg3N2 as nitriding agents [J]. Polyhedron, 1998, 17(16): 2617–2622. doi: 10.1016/S0277-5387(97)00454-3
    [8] KOBASHI M, OKAYAMA N, CHOH T. Synthesis of AlN/Al alloy composites by in situ reaction between Mg3N2 and aluminum [J]. Materials Transactions, JIM, 1997, 38(3): 260–265. doi: 10.2320/matertrans1989.38.260
    [9] EREMETS M I, GAVRILIUK A G, TROJAN I A, et al. Single-bonded cubic form of nitrogen [J]. Nature Materials, 2004, 3(8): 558–563. doi: 10.1038/nmat1146
    [10] LORENZ H, KÜHNE U, HOHLFELD C, et al. Influence of MgO on the growth of cubic boron nitride using the catalyst Mg3N2 [J]. Journal of Materials Science Letters, 1988, 7(1): 23–24. doi: 10.1007/BF01729904
    [11] PAPACONSTANTOPOULOS D A, PICKETT W E, KLEIN B M, et al. Electronic properties of transition-metal nitrides: the group-V and group-VI nitrides VN, NbN, TaN, CrN, MoN, and WN [J]. Physical Review B, 1985, 31(2): 752–761. doi: 10.1103/PhysRevB.31.752
    [12] MENG Y, MAO H K, ENG P J, et al. The formation of sp3 bonding in compressed BN [J]. Nature Materials, 2004, 3(2): 111–114. doi: 10.1038/nmat1060
    [13] RAZA Z, PICKARD C J, PINILLA C, et al. High energy density mixed polymeric phase from carbon monoxide and nitrogen [J]. Physical Review Letters, 2013, 111(23): 235501. doi: 10.1103/PhysRevLett.111.235501
    [14] KNITTLE E, WENTZCOVITCH R M, JEANLOZ R, et al. Experimental and theoretical equation of state of cubic boron nitride [J]. Nature, 1989, 337(6205): 349–352. doi: 10.1038/337349a0
    [15] XIA Y, LI Q, MA Y M. Novel superhard polymorphs of Be3N2 predicted by first-principles [J]. Computational Materials Science, 2010, 49(1): S76–S79. doi: 10.1016/j.commatsci.2010.01.045
    [16] ZHU S S, PENG F, LIU H Y, et al. Stable calcium nitrides at ambient and high pressures [J]. Inorganic Chemistry, 2016, 55(15): 7550–7555. doi: 10.1021/acs.inorgchem.6b00948
    [17] JI D P, CHONG X Y, FENG J. Electronic, mechanical and hydrogen storage properties of novel Mg3N2 [J]. Journal of Alloys and Compounds, 2019, 800: 8–15. doi: 10.1016/j.jallcom.2019.06.021
    [18] WANG Y C, LV J, ZHU L, et al. Crystal structure prediction via particle-swarm optimization [J]. Physical Review B, 2010, 82(9): 094116. doi: 10.1103/PhysRevB.82.094116
    [19] MONKHORST H J, PACK J D. Special points for Brillouin-zone integrations [J]. Physical Review B, 1976, 13(12): 5188. doi: 10.1103/PhysRevB.13.5188
    [20] PERDEW J P, BURKE K, ERNZERHOF M. Generalized gradient approximation made simple [J]. Physical Review Letters, 1996, 77(18): 3865–3868. doi: 10.1103/PhysRevLett.77.3865
    [21] BORN M, HUANG K. Dynamical theory of crystal lattices [J]. American Journal of Physics, 1954, 39(2): 113–127.
    [22] BECKE A D, EDGECOMBE K E. A simple measure of electron localization in atomic and molecular systems [J]. The Journal of Chemical Physics, 1990, 92(9): 5397–5403. doi: 10.1063/1.458517
    [23] BADER R F W. Atoms in molecules [J]. Accounts of Chemical Research, 1985, 18(1): 9–15. doi: 10.1021/ar00109a003
    [24] OGANOV A R, CHEN J H, GATTI C, et al. Ionic high-pressure form of elemental boron [J]. Nature, 2009, 457(7231): 863–867. doi: 10.1038/nature07736
  • [1] 刘思远缪宇马雪姣李鑫高文泉程宇衡刘艳辉 . 基于第一性原理计算IrSb压力相变. 高压物理学报, 2019, 33(5): 052203-1-052203-7. doi: 10.11858/gywlxb.20190716
    [2] 李鑫马雪姣高文泉刘艳辉 . 高压下Ir2P晶体结构预测与物理性质. 高压物理学报, 2019, 33(1): 011103-1-011103-6. doi: 10.11858/gywlxb.20180645
    [3] 丁迎春刘海军蒋孟衡陈敏陈勇明 . 高压下BeP2N4结构相变和电子结构的第一性原理计算. 高压物理学报, 2012, 26(6): 674-680. doi: 10.11858/gywlxb.2012.06.012
    [4] 韩林马麦宁徐志双周晓亚 . 辉石结构与相变的第一性原理研究. 高压物理学报, 2017, 31(2): 125-134. doi: 10.11858/gywlxb.2017.02.004
    [5] 杨龙星刘雷刘红易丽顾小雨 . 高压下石榴子石结构和弹性的第一性原理研究. 高压物理学报, 2019, 33(6): 060104-1-060104-11. doi: 10.11858/gywlxb.20190785
    [6] 郝军华吴志强王铮金庆华李宝会丁大同 . 高压下SiO2的第一性原理计算. 高压物理学报, 2010, 24(4): 260-266 . doi: 10.11858/gywlxb.2010.04.004
    [7] 谭心贾亦超刘学杰 . 外压下ZrN结构相变的第一性原理研究. 高压物理学报, 2014, 28(2): 168-174. doi: 10.11858/gywlxb.2014.02.006
    [8] 杨晓翠赵玉伟高忠明刘鑫张立新王晓明郝爱民 . 高压下CaF2结构相变和光学性质的第一性原理计算. 高压物理学报, 2010, 24(3): 225-230 . doi: 10.11858/gywlxb.2010.03.011
    [9] 刘克伟于杰周晓龙胡明钰占建翔 . 高压下P42/mnmPnnm晶体结构SnO2密度泛函研究. 高压物理学报, 2017, 31(1): 81-88. doi: 10.11858/gywlxb.2017.01.012
    [10] 郝爱民杨晓翠赵玉伟刘鑫宋爱君张卫国辛伟 . 高压下钡的硫化物的结构相变和光学性质研究. 高压物理学报, 2010, 24(2): 113-119 . doi: 10.11858/gywlxb.2010.02.006
    [11] 唐旭东张增明赵智王中平丁泽军 . 高压诱致的反式联苯乙烯酮结构相变与化学反应. 高压物理学报, 2009, 23(6): 401-406 . doi: 10.11858/gywlxb.2009.06.001
    [12] 温新竹彭玉颜刘明真 . 钙钛矿结构ZrBeO3稳定性的第一性原理研究. 高压物理学报, 2020, 34(1): 011202-1-011202-9. doi: 10.11858/gywlxb.20190802
    [13] 武琪彭放李庆华雷力李荣祺 . LiAl5O8高压相变研究. 高压物理学报, 2012, 26(3): 338-342. doi: 10.11858/gywlxb.2012.03.015
    [14] 张程祥吴绍曾吴姝妍吴雪原 . NaCl、KCl和NaF晶体高压结构相变点的计算. 高压物理学报, 1993, 7(3): 232-237 . doi: 10.11858/gywlxb.1993.03.011
    [15] 何林尹君 . 空位点缺陷对金刚石在高压下的电子结构和光学性质的影响. 高压物理学报, 2013, 27(6): 802-806. doi: 10.11858/gywlxb.2013.06.002
    [16] 任晓光崔雪菡吴宝嘉顾广瑞 . 高压下金属间化合物CaAlSi的电子结构和晶格动力学性质. 高压物理学报, 2014, 28(2): 161-167. doi: 10.11858/gywlxb.2014.02.005
    [17] 刘玉孝伍丽莎齐云姚勇 . ZnSe高压相变和声子谱的第一性原理研究. 高压物理学报, 2017, 31(2): 97-102. doi: 10.11858/gywlxb.2017.02.001
    [18] 王汝菊李凤英姚玉书王积方林学仕 . 用超声方法研究掺丙氨酸氘化TGS晶体的相变. 高压物理学报, 1989, 3(2): 156-158 . doi: 10.11858/gywlxb.1989.02.009
    [19] 王霖刘冰冰王卉侯元元艾晓雷潘跃武崔启良邹广田刘洪江倪永红 . 纳米硫化锌球壳的高压相变研究. 高压物理学报, 2005, 19(4): 357-360 . doi: 10.11858/gywlxb.2005.04.013
    [20] 靳锡联崔田 . 高压下典型小分子晶体的结构和行为. 高压物理学报, 2013, 27(2): 188-198. doi: 10.11858/gywlxb.2013.02.004
  • 加载中
图(5)表(3)
计量
  • 文章访问数:  102
  • 阅读全文浏览量:  97
  • PDF下载量:  2
出版历程
  • 收稿日期:  2019-08-06
  • 录用日期:  2019-09-25
  • 网络出版日期:  2020-01-14
  • 刊出日期:  2020-02-01

高压下MgN8晶体结构理论模拟与物性研究

    作者简介:缪 宇(1999-),男,本科,主要从事高压结构计算研究. E-mail: 15053219411@163.com
    通讯作者: 金哲学, jinzhexue@ybu.edu.cn
  • 1. 延边大学理学院,吉林 延吉 133000
  • 2. 延边大学工学院,吉林 延吉 133000

摘要: 基于密度泛函理论第一性原理的方法,使用CALYPSO结构搜索技术结合VASP软件,在0~100 GPa压强范围内对MgN8的晶体结构进行预测,并对预测的结构进行系统研究。结果表明:在常压下,空间群为P4/mbmα-MgN8晶体结构的焓值最低;当压强达到24.3 GPa和68.3 GPa时发生相变,分别相变成空间群为P4/mncβ-MgN8相和空间群为Cmcmγ-MgN8相,两次相变均为对应体积坍塌的一级相变。电子性质计算结果表明,α-MgN8相的导带与价带之间具有3.09 eV的带隙,表明该结构具有非金属性;β相和γ相具有明显的金属特征。Bader电荷转移计算表明,随着压力的增加,Mg原子向N原子转移的电荷逐渐增多。

English Abstract

  • 随着材料科学研究的深入,研究者们对材料的研究不再局限于常压范围,开始更多地关注高压下材料的结构和物理性质。在高压作用下,晶体的体积被压缩,原子间的间距变小,电子轨道发生变化,甚至原子的周期性排列也会改变,使晶体结构发生相变[1-3]。晶体结构在高压环境下具有较为丰富的相变序列,许多常压下无法合成的材料可在高压环境下合成,因此利用高压技术,可以丰富人们对物质世界的认识[4-5]

    高压下聚合氮以及氮化物的物理、化学性质受到较多关注。将氮作为媒介来合成不同结构的含氮材料是研究聚合氮及氮化物的基础[6-8],如利用氮单质在高温高压环境下合成立方氮[9],利用氮元素和其他元素合成具有超硬性质的化合物等[10]。在氮的化合物中,金属氮化物的潜在应用价值巨大,如超导材料[11]、特殊化学制品[12]、热稳材料、高能量密度材料和超硬材料等[13-14]。针对金属氮化物的研究也愈发广泛。

    对于碱土金属氮化物的理论研究较多。Be3N2的晶体结构研究结果表明,R3m$P\bar 3m1$两个晶体结构的维氏硬度分别为51 GPa和54 GPa,说明两种结构的Be3N2是潜在的超硬材料[15]。在高压下发现Ca与N的多个配比结构,得出CaN3C2/c结构中N原子形成孤立的六环;进一步对键长和键级进行计算,发现该结构的CaN3是潜在的高能量密度材料[16]。Mg3N2在高压下具有丰富的物理性质。研究发现,H2分子很容易在τ-Mg3N2α-Mg3N2α′-Mg3N2相的表面分解为H原子并进入材料内部,说明Mg3N2是潜在的储氢材料[17]

    由上述研究可知,碱土金属氮化物的物理性质丰富,潜在应用价值巨大。为此,对高压下Mg的多氮化合物的晶体结构和物理性质进行研究具有十分重要的意义。本研究结合晶体结构预测技术和基于密度泛函理论的第一性原理方法,在0~100 GPa的压强范围内,对MgN8的晶体结构进行预测,并系统分析预测得到的α-MgN8β-MgN8γ-MgN8相,所得结论对进一步研究碱土金属氮化物具有一定的参考价值。

    • 采用基于粒子群优化算法的CALYPSO软件,对MgN8进行结构搜索[18],压强范围为0~100 GPa,模拟晶胞采用2倍胞和4倍胞。对于搜索出的结构,使用VASP软件包进行结构优化[19],其中描述电子间的交换关联势能采用广义梯度近似(GGA)下的Perdew-Burke-Ernzerh(PBE)交换关联泛函[19],赝势采用全电子投影缀加平面波方法[20],Mg原子的价电子为2p63s2,N原子的价电子为2s22p3。为了保证计算精确,使能量收敛精度小于1 meV/atom,经过能量收敛测试,得到平面波的截断能为600 eV;第一布里渊区积分采用Monkhorst-Pack网格方法,网格的间距为0.2 nm–1;优化过程中,在确保空间群不改变的情况下,选用10–5 eV为自洽能量收敛的最小值,优化应力收敛设置为0.001 eV/Å。对于优化后的结构,使用超晶胞的方法,利用PHONOPY软件计算声子色散关系和声子态密度,选取不同的扩胞数和k点网格:α-MgN8相为1 × 1 × 1和7 × 7 × 6;β-MgN8相为2 × 2 × 4和6 × 6 × 9;γ-MgN8相为3 × 1 × 1和9 × 5 × 4。使用VASP软件包计算电子局域函数。

    • 运用CALYPSO软件搜索,在0~100 GPa压强范围内,对预测出的MgN8晶体结构进行晶格常数和原子位置的优化,计算时设定的温度为0 K,根据G = HTSG为吉布斯自由能,H为焓,T为温度,S为熵),可以用焓代替系统的自由能。通过计算优化后结构的焓值随压强的变化关系,可以得到晶体结构的热力学稳定区间。绘制出的焓差随压强变化的曲线如图1所示。由图1(a)可知:在常压下,空间群为P4/mbm的MgN8晶体结构的焓值最低,记为α-MgN8相;当压强达到24.3 GPa时,空间群为P4/mnc的晶体结构具有更低的焓值,记为β-MgN8相;当压强达到68.3 GPa时,空间群为Cmcm的晶体结构能量更低,记为γ-MgN8相。图1(b)给出了α-MgN8β-MgN8γ-MgN8相的晶胞体积随压强的变化。分析可知:α-MgN8相的晶胞体积随着压强的增大而减小,晶胞体积随压强连续变化;当压强达到24.3 GPa时,MgN8发生体积坍缩,坍缩率为18.87%;在68.3 GPa时,发生结构相变,体积塌缩率为2.62%;两次相变均属于一级相变。

      图  1  MgN8的焓差曲线以及α-MgN8β-MgN8γ-MgN8相体积随压强的变化关系

      Figure 1.  Calculated enthalpies per formula unit of pressure of MgN8 and the calculated volume versus pressure of α-MgN8, β-MgN8 and γ-MgN8

      α-MgN8β-MgN8γ-MgN8相的晶体结构见图2,其平衡态晶格常数以及原子位置信息列于表1表2

      PhasePressure/GPaabcα/(°)β/(°)γ/(°)
      α-MgN8 (P4/mbm)0 5.9135.9136.57290.090.090.0
      β-MgN8 (P4/mnc)24.36.2196.2193.72490.090.090.0
      γ-MgN8 (Cmcm)68.34.1674.1678.68090.090.051.1

      表 1  α-MgN8β-MgN8γ-MgN8相的平衡态晶格常数

      Table 1.  Lattice parameters of α-MgN8, β-MgN8 and γ-MgN8 in the equilibrium state

      图  2  α-MgN8β-MgN8γ-MgN8相的晶体结构

      Figure 2.  Crystal structures of α-MgN8, β-MgN8 and γ-MgN8

      PhasePressure/GPaAtomsWyckoff positionxyz
      α-MgN8 (P4/mbm)0 Mg 2b0.0000.0000.500
      N1 2a0.0000.0000.000
      N2 2c0.0000.5000.500
      N3 4e0.0000.0000.821
      N4 4f0.0000.5000.915
      N5 4h0.8590.3590.500
      β-MgN8 (P4/mnc)24.3Mg 2a0.0000.0001.000
      N16i0.4160.2190.823
      γ-MgN8 (Cmcm)68.3Mg 8g0.624–0.376 1.250
      N116h0.226–0.135 0.617
      N216h1.1150.4711.386

      表 2  α-MgN8β-MgN8γ-MgN8相的平衡态晶格常数和原子位置

      Table 2.  Lattice parameters and atomic coordinate of α-MgN8, β-MgN8 and γ-MgN8 in the equilibrium state

      图2(a)所示的结构为0 GPa下优化后的α-MgN8相结构。在α-MgN8相结构中,每个Mg原子被6个N原子包围,构成八面体结构。其中,Mg-N的键长有2类:一类是中心Mg原子与八面体上下顶点的两个N相连,键长为2.280 Å;另一类为中心的Mg原子与八面体中其余4个N相连,键长为2.109 Å。α-MgN8相的晶格结构常数:a = b = 5.913 Å,c = 6.572 Å,α = β = γ = 90.0°。其中,Mg原子的Wyckoff占位为2b (0.000, 0.000, 0.500),N原子的Wyckoff占位分别为2a (0.000, 0.000, 0.000)、2c (0.000, 0.500, 0.500)、4e (0.000, 0.000, 0.821)、4f (0.000, 0.500, 0.915)和4h (0.859, 0.359, 0.500)。图2(b)所示的结构为24.3 GPa下优化后的β-MgN8相结构。在该结构中,每个Mg原子被8个N原子包围,8支Mg—N键的键长相同,为2.184 Å,构成正六面体。该结构的晶体结构常数:a = b = 6.219 Å,c = 3.724 Å,α = β = γ = 90.0°。其中,Mg原子的Wyckoff占位为2a (0.000, 0.000, 1.000),N原子Wyckoff占位为16i (0.416, 0.219, 0.823)。图2(c)所示的结构为68.3 GPa下优化后的γ-MgN8相结构。在该结构中,每个Mg原子与8个N原子相连,构成六面体结构。六面体中的8支Mg-N键有2种键长:4支2.081 Å和4支2.069 Å。该结构的晶体结构常数:a = 4.167 Å,b = 4.167 Å,c = 8.680 Å,α = β = 90.0°,γ = 51.1°。其中,Mg原子的Wyckoff占位为8g (0.624, –0.376, 1.250),N原子Wyckoff占位分别为16h (0.226, –0.135, 0.617)和16h (1.115, 0.471, 1.386)。

    • 通过计算声子色散关系,可以判断晶格结构是否具有动力学稳定性。晶格结构具有动力学稳定性的判断条件是,所有简正声子频率都为有限的实值[21],如果是虚值,就可以判断此结构出现了声子软化现象,晶格结构不稳定。计算了预测结构的声子谱和声子态密度,如图3所示。

      图  3  α-MgN8β-MgN8γ-MgN8相声子色散关系和声子态密度

      Figure 3.  Phonon-dispersion curves and the phonon density of states of α-MgN8, β-MgN8 and γ-MgN8

      图3(a)可知,0 GPa时α-MgN8相在整个布里渊区未出现声子软化现象,说明该相具有动力学稳定性,最大光学支的频率为71.35 THz。由图3(b)图3(c)可知,24.3 GPa和68.3 GPa时,β-MgN8γ-MgN8相具有动力学稳定性,最大光学支频率分别为44.64 THz和45.63 THz。

    • 为了探究预测得到的MgN8晶体结构的电子性质,分别研究了α-MgN8相在0 GPa、β-MgN8相在24.3 GPa和γ-MgN8相在68.3 GPa时的能带结构和电子态密度,如图4所示。

      图  4  α-MgN8β-MgN8γ-MgN8相的能带结构和电子态密度

      Figure 4.  Band structure and partial density of states of α-MgN8, β-MgN8 and γ-MgN8

      由能带图可知,压强为0 GPa时,α-MgN8相的导带与价带之间有3.093 eV的带隙,表明该相结构具有非金属性;当压强为24.3 GPa和68.3 GPa时,β-MgN8γ-MgN8相的导带与价带跨越费米面发生交叠,表明这两相结构具有金属性。

      由电子态密度的计算结果可知,α-MgN8β-MgN8γ-MgN8相费米能级处的电子态密度均主要由N原子的p轨道贡献,而N原子的s轨道和Mg原子的sp轨道对电子态密度的贡献相对较少。N原子的p轨道与s轨道之间存在杂化,形成了N—N共价键。

      为了研究α-MgN8β-MgN8γ-MgN8相晶体结构的化学键,计算了这3个相的电子局域函数(Electron localization function, ELF)[22],如图5所示。

      图  5  α-MgN8β-MgN8γ-MgN8相的电子局域函数

      Figure 5.  Electron localization function of α-MgN8, β-MgN8 and γ-MgN8

      图5(a)显示了α-MgN8结构的三维电子局域函数。α-MgN8相的N原子之间形成N2和N3,其中:N2的两个N原子之间存在电子局域,形成了共价键,键长为1.116 Å,两个N原子外围存在孤对电子;N3中间的N原子与相邻的两个N原子之间存在电子局域,形成两种N3共价键,其键长分别为1.178 Å和1.182 Å,N3两侧的N原子周围也存在孤对电子。这3种N-N键的键长介于N≡N三键(1.10 Å)和N=N双键(1.25 Å)之间。图5(b)显示了压强为24.3 GPa时β-MgN8结构的三维电子局域函数。该相N原子之间形成链状结构,每个N原子与相邻的两个N原子之间均有电子局域,形成共价键,键长为1.308 Å和1.316 Å,这两种N-N键的键长介于N=N双键(1.25 Å)和N—N单键(1.45 Å)之间。图5(c)显示68.3 GPa下γ-MgN8结构的三维电子局域函数。由图5(c)可知,N原子之间形成链状结构,每个N原子与相邻的两个N原子之间均有电子局域,形成共价键,键长为1.282、1.297和1.319 Å,3种N-N键的键长介于N=N双键(1.25 Å)和N—N单键(1.45 Å)之间。比较上述3相的N原子键长数据可以发现,β相和γ相的N-N键更容易断裂,释放出大量的能量,从而成为潜在的高能量密度材料。

      图5(d)图5(e)图5(f)为MgN8晶体3种相ELF的二维切面图,晶格切面分别选取(010)、(110)和(010),ELF的数值范围为0~1。ELF = 1表示电子完全局域化,ELF = 0表示电子完全离域化或该处没有电子。从图5(d)图5(e)图5(f)中可以看出:在Mg原子与N原子之间存在电子局域,并且N原子附近对应的ELF值接近1,为高度局域的电子分布,说明Mg原子与N原子之间存在极性共价键;N原子与N原子之间对应的ELF值接近1,表明N原子与N原子之间存在共价键。

      为了清晰地描述N原子与Mg原子之间的电子转移情况,计算了Bader电荷转移[23-24],计算结果见表3。N原子的电负性强于Mg原子,故N原子吸引电子的能力更强,Mg原子与N原子之间电荷转移的施主为Mg原子,受主为N原子。当压强为0 GPa时,α-MgN8相每个Mg原子向周围N原子的电荷转移为1.74e。由于该结构中的N原子不仅与Mg原子相连,还形成了N2和N3结构,N2与N3结构的存在导致N原子既可以作为受主,接收来自Mg原子和其他N原子的电荷,又可以作为施主,将电荷转移给其他N原子;当压强为24.3 GPa时,β-MgN8相中每个Mg原子向周围N原子的电荷转移为1.94e;当压强为68.3 GPa时,γ-MgN8相中每个Mg原子向周围N原子的电荷转移为2.20e。由此可见,随着压强的增加,电荷转移逐渐增多。

      PhasePressure/GPaAtomsNumberCharge value/eCharge transfer/e
      α-MgN8 (P4/mbm)0Mg16.26 1.74
      N114.84 0.16
      N214.86 0.13
      N325.04–0.04
      N425.46–0.46
      N525.52–0.52
      β-MgN8 (P4/mnc)24.3Mg16.06 1.94
      N85.24–0.24
      γ-MgN8 (Cmcm)68.3Mg15.80 2.20
      N145.28–0.28
      N245.27–0.27

      表 3  α-MgN8β-MgN8γ-MgN8相电荷转移

      Table 3.  Calculated Bader charges of α-MgN8, β-MgN8 and γ-MgN8

    • 运用CALYPASO软件并结合密度泛函理论第一性原理计算方法,对MgN8的晶体结构进行预测,同时对预测得出的晶体结构进行了结构优化。研究表明:在0~100 GPa范围内,3个相的相变序列为P4/mbmP4/mncCmcm;两次相变过程中,体积均发生坍塌,坍塌率分别为18.87%和2.62%。电子性质研究表明:α-MgN8相结构具有非金属性,带隙为3.093 eV;β-MgN8γ-MgN8相结构具有明显的金属性。电子局域函数显示,Mg原子与N原子之间形成极性共价键,N原子与N原子之间形成共价键。高压β-MgN8γ-MgN8相更适合用作储能材料,有望成为高能量密度材料的候选。Bader电荷转移计算结果表明,N原子具有更强的电负性,Mg原子与N原子之间电荷转移的施主为Mg原子,受主为N原子。在α-MgN8β-MgN8γ-MgN8相结构中,每个Mg原子向周围N原子的电荷转移分别为1.74e、1.94e和2.20e,即随着压强的增加,电荷转移增多。研究结果可以为碱土金属氮化物的进一步研究提供参考。

参考文献 (24)

目录

    /

    返回文章
    返回