搜索

x

留言板

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

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

第一性原理研究Li2NH的晶格动力学和热力学性质

赵玉娜 丛红璐 成爽 于娜 高涛 马俊刚

引用本文:
Citation:

第一性原理研究Li2NH的晶格动力学和热力学性质

赵玉娜, 丛红璐, 成爽, 于娜, 高涛, 马俊刚

First-principles study of lattice dynamical and thermodynamic properties of Li2NH

Zhao Yu-Na, Cong Hong-Lu, Cheng Shuang, Yu Na, Gao Tao, Ma Jun-Gang
PDF
HTML
导出引用
  • 采用第一性原理的赝势平面波方法系统地研究了Li2NH的电子结构、晶格动力学和热力学性质. 计算得到的晶格常数与先前的理论和实验结果符合得很好. 运用线性响应理论计算了整个布里渊区高对称方向上的声子色散曲线和相应的声子态密度, 发现Li2NH(Pnma)声子色散曲线没有虚频, 动力学性能相对最稳定, 计算结果和先前实验及理论数据符合得很好. 最后, 利用得到的声子态密度进一步预测了Li2NH的热力学性质, 包括晶格振动对Helmholtz自由能、内能、熵和热容的贡献, 计算结果在一定程度上可为Li-N-H储氢体系的应用提供理论指导.
    One of the key issues for scale applications of hydrogen energy is the availability of safe, efficient and ecnomicical hydrogen storage technologies. In the past few years, light metal hydrides have attracted considerable attention due to their high hydrogen capacity. With a hydrogen capacity up to ~6.5 wt%, Li2NH is regarded as one of the most promising hydrogen storage materials. Although the hydrogen physical and thermodynamic properties of Li2NH have been studied, the electronic structure, phonon vibration mode and thermodynamic properties of Li2NH have not yet been resolved. In this paper, by using the first principles based on the density functional theory (DFT), we investigate the electronic structure, lattice dynamical and thermodynamic properties of Li2NH in detail.Firstly, the structure of Li2NH is optimized and the lattice parameters and total energy of the crystals are calculated. As shown by the calculation results, the lattice parameters are in good agreement with previous theoretical and experimental results. Our lowest-energy structure of Li2NH has orthorhombic Pnma symmetry at T=0 K for all of the proposed structures. Secondly, the electronic band-structure studies reveal that Li2NH has a small band gap of about 2.0 eV. The analysis of total and partial density of states of Li2NH show that the bonding between the N and H has a covalent character. Thirdly, the lattice dynamical properties of Li2NH are investgated at the corresponding equilibrium states. These results show that only the phonon dispersion curves of Li2NH (Pnma) without negative frequencies are calculated along the high-symmetry points. The optical modes of phonon frequencies at Γ point are assigned as Raman and Infrared-active modes. Based on the calculated phonon density of states, the thermodynamic properties are computed, such as the Helmholtz free energy, internal energy, entropy and the constant-volume specific heat versus temperature. The calculation results may explore the applications in areas of hydrogen storage for Li-N-H, which is of great importance forusing hydrogen in the future.
      通信作者: 马俊刚, jgma@bjtuhbxy.cn
    • 基金项目: 国家自然科学基金应急管理项目(理论物理专项)(批准号:11547224)和河北省高等教育教学改革研究与实践项目(批准号:2018GJJG649)资助的课题.
      Corresponding author: Ma Jun-Gang, jgma@bjtuhbxy.cn
    • Funds: Project supported by the Emergency Management Program of the National Natural Science Foundation of China(Special Funds for Theoretical Physics)(Grant No. 11547224) and the Research and Practice of the Higher Education reform in Hebei Province, China (Grant No. 2018GJJG649).
    [1]

    Chen P, Xiong Z, Luo J, Lin J, Tan K L 2002 Nature 420 302Google Scholar

    [2]

    Ohoyama K, Nakamori Y, Orimo S, Yamada K 2005 J. Phys. Soc. Jpn. 74 483Google Scholar

    [3]

    Noritake T, Nozaki H, Aoki M, Towata S, Kitahara G, Nakamori Y, Orimo S 2005 J. Alloys Compd. 393 264Google Scholar

    [4]

    Herbst J F, Hector Jr L G 2005 Phys. Rev. B 72 125120Google Scholar

    [5]

    Balogh M P, Jones C Y, Herbst J F, Hector Jr. L G, Kundrat M 2006 J. Alloys Compd. 420 326Google Scholar

    [6]

    Mueller T, Ceder G 2006 Phys. Rev. B 74 134104Google Scholar

    [7]

    Kojima Y, Kawai Y 2005 J. Alloys Compd. 395 236Google Scholar

    [8]

    Magyari-Köpe B, Ozoliņš V, Wolverton C 2006 Phys. Rev. B 73 220101(R)Google Scholar

    [9]

    Song Y, Guo Z X 2006 Phys. Rev. B 74 195120Google Scholar

    [10]

    Hector Jr L G, Herbst J F 2008 J. Phys. Condens. Matter 20 064229Google Scholar

    [11]

    Yang J, Lamsal J, Cai Q, Yelon W B, James W J 2008 MRS Proceedings 1098 1098 1098−HH03-06

    [12]

    Miceli G, Cucinotta C, Bernasconi M, Parrinello M 2010 J. Phys. Chem. C 114 15174Google Scholar

    [13]

    Miceli G, Ceriotti M, Angioletti-Uberti S, Bernasconi M, Parrinello M 2011 J. Phys. Chem. C 115 7076Google Scholar

    [14]

    Wolverton C, Siegel J D, Akbarazadeh R A, Ozolis V 2008 J. Phys. Condens. Matter 20 064228Google Scholar

    [15]

    陈玉红, 吕晓霞, 杜瑞, 董肖, 张材荣, 康龙, 罗永春 2013 稀有金属材料与工程 4 2

    Chen Y H, Lv X X, Du R, Dong X, Zhang C R, Kang L, Luo Y C 2013 Rare Metal Materials and Engineering 4 2

    [16]

    Rajeswarapalanichamy R, Santhosh M, Sudhapriyanga G, Kanagaprabha S, Iyakutti K 2015 Acta Metall. Sinica 28 550Google Scholar

    [17]

    Crivello J C, Gupta M, Černý R, Latroche M, Chandra D 2010 Phys. Rev. B 81 104113Google Scholar

    [18]

    Wang Q, ChenY G, Zheng X, Niu G, Wu C L, Tao M D T 2009 Physica B 404 3431Google Scholar

    [19]

    The ABINIT code is a common project of the Université Catholique de Louvain, and other contributors (URL http://www.abinit.org)

    [20]

    Troullier N, Martins J L 1991 Phys. Rev. B 43 1993Google Scholar

    [21]

    Monkhorst H J, Pack J D 1976 Phys. Rev. B 13 5188Google Scholar

    [22]

    Baroni S, Giannozzi P, Testa A 1987 Phys. Rev. Lett. 58 1861Google Scholar

    [23]

    Baroni S, de Gironcoli S, Dal Corso A, Giannozzi P 2001 Rev. Mod. Phys. 73 515Google Scholar

    [24]

    Lee C, Gonze X 1995 Phys. Rev. B 51 8610Google Scholar

    [25]

    Li W 2011 Ph. D. Dissertation (Zhejiang: Zhejiang University)

    [26]

    Aulbur W G, Jonsson L, Wilkins J W 2000 Solid States Phys. 54 1Google Scholar

    [27]

    Stampfl C, Van de Walle C G 1999 Phys. Rev. B 5 9

    [28]

    Gupta M, Gupta R P 2007 J. Alloys Compd. 319 446

    [29]

    Yao J H 2007 Ph. D. Dissertation (London: University of London) 0-239

    [30]

    陈玉红2008 博士学位论文(兰州: 兰州理工大学)

    Chen Y H 2008 Ph. D. Dissertation(Lanzhou: Lanzhou University of Technology)(in Chinese)

    [31]

    Born M, Huang K 1954 Dynamical Theory of Crystal Lattices (Oxford:Oxford University Press) p121

    [32]

    Maradudin A A, Montroll E W, Weiss C H, Ipatova I P 1971 Theory of Lattice Dynamics in the Harmonic Approximation (2nd Ed.) (New York: Academic Press)

  • 图 1  Li2NH的三种晶体结构图 (a) Li2NH ($F\bar 43m$); (b) Li2NH ($Fm\bar 3m$); (c) Li2NH ($Pnma$)

    Fig. 1.  Crystal structures of Li2NH: (a) Li2NH ($F\bar 43m$); (b) Li2NH ($Fm\bar 3m$); (c) Li2NH ($Pnma$).

    图 2  Li2NH的能带结构

    Fig. 2.  Band structure of Li2NH.

    图 3  Li2NH的总态密度和分波态密度

    Fig. 3.  Total and partial density of states of Li2NH(states/eV)

    图 4  Li2NH($Pnma$)的声子色散曲线

    Fig. 4.  The phonon dispersion curves of Li2NH ($Pnma$).

    图 5  Li2NH($F\bar 43m$)的声子色散曲线

    Fig. 5.  The phonon dispersion curves of Li2NH ($F\bar 43m$).

    图 6  Li2NH($Fm\bar 3m$)的声子色散曲线

    Fig. 6.  The phonon dispersion curves of Li2NH ($Fm\bar 3m$)

    图 7  Li2NH($Pnma$)的声子态密度

    Fig. 7.  The phonon density of states of Li2NH ($Pnma$).

    图 8  Li2NH的热力学性质 (a)晶格振动对自由能的贡献ΔF; (b)内能ΔE随温度的变化; (c)熵S随温度的变化; (d)热容Cv随温度的变化 (▼表示Debye的T3定律)

    Fig. 8.  Thermodynamic properties for Li2NH: (a) The phonon contribution to the free energies ΔF of Li2NH; (b)the internal energies ΔE verus temperature; (c)the entropy S verus temperature; (d) the constant-volume specific heats Cv verus temperature(▼ meansT3 law of Debye)

    表 1  Li2NH晶体结构参数

    Table 1.  Lattice parameters of Li2NH.

    Space groupLattice parameter Etot/eV
    PresentExperimentalOther CalculatedPresent
    $F\bar 43m$a = 5.6249a = 5.0769[2]a = 5.159[17],5.223[17],5.074[17],5.649[16],5.1076[18],5.2968[9]–1920.4448
    $Pnma$a = 7.781
    b = 3.623
    c = 4.902
    a = 7.733[8],7.704[16],7.742[17],7.753[17],7.775[17]
    b = 3.60[8],3.75[16],3.604[17],3.609[17],3.618[17]
    c = 4.872[8],5.074[16],4.883[17],4.890[17],4.881[17]
    –2702.3588
    $Fm\bar 3m$a = 5.0696a = 5.0742[3]a = 5.047[3],5.045[16]–1913.6655
    下载: 导出CSV

    表 2  Li2NH($Pnma$)在布里渊区中心Γ点的光学模声子频率(单位:cm–1)

    Table 2.  Phonon frequencies (unit:cm–1) at the Γ point of Li2NH.

    B2uB1uB3gB2gB3uAuB1gAg
    123.6137.9160161.5168198.5222223.5
    310195.8380232.0289243.0359287
    408318432300350.5446435335
    644450660350.5385625654380
    550412620542
    630580668598
    312264031233132
    3134
    下载: 导出CSV

    表 3  Li2NH的热力学函数

    Table 3.  Thermodynamic functions for Li2NH.

    Temprature/KS/J·(mol-c·K)–1CV/J·(mol-c·K)–1ΔH(HH298)/KJ·mol–1
    298248.9273.80
    300250.2274.92.9
    500405.4327.155.1
    700519.0347.075.0
    900607.7358.586.5
    1100680.5366.694.2
    1300742.3372.8100.8
    1500796.0377.5105.5
    下载: 导出CSV
  • [1]

    Chen P, Xiong Z, Luo J, Lin J, Tan K L 2002 Nature 420 302Google Scholar

    [2]

    Ohoyama K, Nakamori Y, Orimo S, Yamada K 2005 J. Phys. Soc. Jpn. 74 483Google Scholar

    [3]

    Noritake T, Nozaki H, Aoki M, Towata S, Kitahara G, Nakamori Y, Orimo S 2005 J. Alloys Compd. 393 264Google Scholar

    [4]

    Herbst J F, Hector Jr L G 2005 Phys. Rev. B 72 125120Google Scholar

    [5]

    Balogh M P, Jones C Y, Herbst J F, Hector Jr. L G, Kundrat M 2006 J. Alloys Compd. 420 326Google Scholar

    [6]

    Mueller T, Ceder G 2006 Phys. Rev. B 74 134104Google Scholar

    [7]

    Kojima Y, Kawai Y 2005 J. Alloys Compd. 395 236Google Scholar

    [8]

    Magyari-Köpe B, Ozoliņš V, Wolverton C 2006 Phys. Rev. B 73 220101(R)Google Scholar

    [9]

    Song Y, Guo Z X 2006 Phys. Rev. B 74 195120Google Scholar

    [10]

    Hector Jr L G, Herbst J F 2008 J. Phys. Condens. Matter 20 064229Google Scholar

    [11]

    Yang J, Lamsal J, Cai Q, Yelon W B, James W J 2008 MRS Proceedings 1098 1098 1098−HH03-06

    [12]

    Miceli G, Cucinotta C, Bernasconi M, Parrinello M 2010 J. Phys. Chem. C 114 15174Google Scholar

    [13]

    Miceli G, Ceriotti M, Angioletti-Uberti S, Bernasconi M, Parrinello M 2011 J. Phys. Chem. C 115 7076Google Scholar

    [14]

    Wolverton C, Siegel J D, Akbarazadeh R A, Ozolis V 2008 J. Phys. Condens. Matter 20 064228Google Scholar

    [15]

    陈玉红, 吕晓霞, 杜瑞, 董肖, 张材荣, 康龙, 罗永春 2013 稀有金属材料与工程 4 2

    Chen Y H, Lv X X, Du R, Dong X, Zhang C R, Kang L, Luo Y C 2013 Rare Metal Materials and Engineering 4 2

    [16]

    Rajeswarapalanichamy R, Santhosh M, Sudhapriyanga G, Kanagaprabha S, Iyakutti K 2015 Acta Metall. Sinica 28 550Google Scholar

    [17]

    Crivello J C, Gupta M, Černý R, Latroche M, Chandra D 2010 Phys. Rev. B 81 104113Google Scholar

    [18]

    Wang Q, ChenY G, Zheng X, Niu G, Wu C L, Tao M D T 2009 Physica B 404 3431Google Scholar

    [19]

    The ABINIT code is a common project of the Université Catholique de Louvain, and other contributors (URL http://www.abinit.org)

    [20]

    Troullier N, Martins J L 1991 Phys. Rev. B 43 1993Google Scholar

    [21]

    Monkhorst H J, Pack J D 1976 Phys. Rev. B 13 5188Google Scholar

    [22]

    Baroni S, Giannozzi P, Testa A 1987 Phys. Rev. Lett. 58 1861Google Scholar

    [23]

    Baroni S, de Gironcoli S, Dal Corso A, Giannozzi P 2001 Rev. Mod. Phys. 73 515Google Scholar

    [24]

    Lee C, Gonze X 1995 Phys. Rev. B 51 8610Google Scholar

    [25]

    Li W 2011 Ph. D. Dissertation (Zhejiang: Zhejiang University)

    [26]

    Aulbur W G, Jonsson L, Wilkins J W 2000 Solid States Phys. 54 1Google Scholar

    [27]

    Stampfl C, Van de Walle C G 1999 Phys. Rev. B 5 9

    [28]

    Gupta M, Gupta R P 2007 J. Alloys Compd. 319 446

    [29]

    Yao J H 2007 Ph. D. Dissertation (London: University of London) 0-239

    [30]

    陈玉红2008 博士学位论文(兰州: 兰州理工大学)

    Chen Y H 2008 Ph. D. Dissertation(Lanzhou: Lanzhou University of Technology)(in Chinese)

    [31]

    Born M, Huang K 1954 Dynamical Theory of Crystal Lattices (Oxford:Oxford University Press) p121

    [32]

    Maradudin A A, Montroll E W, Weiss C H, Ipatova I P 1971 Theory of Lattice Dynamics in the Harmonic Approximation (2nd Ed.) (New York: Academic Press)

  • [1] 黄鳌, 卢志鹏, 周梦, 周晓云, 陶应奇, 孙鹏, 张俊涛, 张廷波. Al和O间隙原子对-Al2O3热力学性质影响的第一性原理计算. 物理学报, 2017, 66(1): 016103. doi: 10.7498/aps.66.016103
    [2] 邓世杰, 赵宇宏, 侯华, 文志勤, 韩培德. 高压下Ti2AlX(X=C,N)的结构、力学性能及热力学性质. 物理学报, 2017, 66(14): 146101. doi: 10.7498/aps.66.146101
    [3] 吴若熙, 刘代俊, 于洋, 杨涛. CaS电子结构和热力学性质的第一性原理计算. 物理学报, 2016, 65(2): 027101. doi: 10.7498/aps.65.027101
    [4] 米尔阿里木江, 艾力, 买买提热夏提, 买买提, 亚森江, 吾甫尔. 非对易相空间中谐振子体系热力学性质的探讨. 物理学报, 2015, 64(14): 140201. doi: 10.7498/aps.64.140201
    [5] 胡翠娥, 曾召益, 蔡灵仓. 极端条件下锆的动力学稳定性研究. 物理学报, 2015, 64(4): 046401. doi: 10.7498/aps.64.046401
    [6] 濮春英, 王丽, 吕林霞, 于荣梅, 何朝政, 卢志文, 周大伟. NbSi2奇异高压相及其热力学性质的第一性原理研究. 物理学报, 2015, 64(8): 087103. doi: 10.7498/aps.64.087103
    [7] 颜小珍, 邝小渝, 毛爱杰, 匡芳光, 王振华, 盛晓伟. 高压下ErNi2B2C弹性性质、电子结构和热力学性质的第一性原理研究. 物理学报, 2013, 62(10): 107402. doi: 10.7498/aps.62.107402
    [8] 赵立凯, 赵二俊, 武志坚. 5d过渡金属二硼化物的结构和热、力学性质的第一性原理计算. 物理学报, 2013, 62(4): 046201. doi: 10.7498/aps.62.046201
    [9] 赵玉娜, 高涛, 吕金钟, 马俊刚. Li-N-H储氢体系热力学性质的第一性原理研究. 物理学报, 2013, 62(14): 143101. doi: 10.7498/aps.62.143101
    [10] 张炜, 陈文周, 王俊斐, 张小东, 姜振益. MnPd合金相变, 弹性和热力学性质的第一性原理研究. 物理学报, 2012, 61(24): 246201. doi: 10.7498/aps.61.246201
    [11] 王斌, 刘颖, 叶金文. 高压下TiC的弹性、电子结构及热力学性质的第一性原理计算. 物理学报, 2012, 61(18): 186501. doi: 10.7498/aps.61.186501
    [12] 杨则金, 令狐荣锋, 程新路, 杨向东. Cr2MC(M=Al, Ga)的电子结构、弹性和热力学性质的第一性原理研究. 物理学报, 2012, 61(4): 046301. doi: 10.7498/aps.61.046301
    [13] 李晓凤, 刘中利, 彭卫民, 赵阿可. 高压下CaPo弹性性质和热力学性质的第一性原理研究. 物理学报, 2011, 60(7): 076501. doi: 10.7498/aps.60.076501
    [14] 忻晓桂, 陈香, 周晶晶, 施思齐. LiFePO4 晶格动力学性质的第一性原理研究. 物理学报, 2011, 60(2): 028201. doi: 10.7498/aps.60.028201
    [15] 李雪梅, 韩会磊, 何光普. LiNH2 的晶格动力学、介电性质和热力学性质第一性原理研究. 物理学报, 2011, 60(8): 087104. doi: 10.7498/aps.60.087104
    [16] 李沛娟, 周薇薇, 唐元昊, 张华, 施思齐. CeO2的电子结构,光学和晶格动力学性质:第一性原理研究. 物理学报, 2010, 59(5): 3426-3431. doi: 10.7498/aps.59.3426
    [17] 李世娜, 刘永. Cu3N弹性和热力学性质的第一性原理研究. 物理学报, 2010, 59(10): 6882-6888. doi: 10.7498/aps.59.6882
    [18] 黄良锋, 李延龄, 倪美燕, 王贤龙, 张国仁, 曾雉. 氢掺杂单层石墨体系的晶格动力学研究. 物理学报, 2009, 58(13): 306-S312. doi: 10.7498/aps.58.306
    [19] 王海燕, 崔红保, 历长云, 李旭升, 王狂飞. AlAs相变及热动力学性质的第一性原理研究. 物理学报, 2009, 58(8): 5598-5603. doi: 10.7498/aps.58.5598
    [20] 刘娜娜, 宋仁伯, 孙翰英, 杜大伟. Mg2Sn电子结构及热力学性质的第一性原理计算. 物理学报, 2008, 57(11): 7145-7150. doi: 10.7498/aps.57.7145
计量
  • 文章访问数:  5339
  • PDF下载量:  71
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-01-24
  • 修回日期:  2019-04-26
  • 上网日期:  2019-07-01
  • 刊出日期:  2019-07-05

第一性原理研究Li2NH的晶格动力学和热力学性质

  • 1. 北京交通大学海滨学院, 黄骅 061199
  • 2. 四川大学原子与分子物理研究所, 成都 610065
  • 通信作者: 马俊刚, jgma@bjtuhbxy.cn
    基金项目: 国家自然科学基金应急管理项目(理论物理专项)(批准号:11547224)和河北省高等教育教学改革研究与实践项目(批准号:2018GJJG649)资助的课题.

摘要: 采用第一性原理的赝势平面波方法系统地研究了Li2NH的电子结构、晶格动力学和热力学性质. 计算得到的晶格常数与先前的理论和实验结果符合得很好. 运用线性响应理论计算了整个布里渊区高对称方向上的声子色散曲线和相应的声子态密度, 发现Li2NH(Pnma)声子色散曲线没有虚频, 动力学性能相对最稳定, 计算结果和先前实验及理论数据符合得很好. 最后, 利用得到的声子态密度进一步预测了Li2NH的热力学性质, 包括晶格振动对Helmholtz自由能、内能、熵和热容的贡献, 计算结果在一定程度上可为Li-N-H储氢体系的应用提供理论指导.

English Abstract

    • 2002年, Chen等[1]提出Li3N理论可逆储氢量10.4 wt%, 储氢反应分两步完成, 即

      $ {\rm{L}}{{\rm{i}}_{\rm{3}}}{\rm{N + }}{{\rm{H}}_2} \Leftrightarrow {\rm{L}}{{\rm{i}}_{\rm{2}}}{\rm{NH + LiH}}, $

      $ {\rm{L}}{{\rm{i}}_{\rm{2}}}{\rm{NH + }}{{\rm{H}}_2} \Leftrightarrow {\rm{LiN}}{{\rm{H}}_{\rm{2}}}{\rm{ + LiH}}, $

      其中第2步反应的可逆储氢高达6.5 wt%, 而且放氢温度低(约150 ℃), 因此在实际应用中主要利用第2步反应. Li2NH作为Li-N-H体系储氢材料的中间生成物, 被认为是非常有前景的储氢材料, 引起各国科研人员的广泛关注.

      据报道, 先前关于Li2NH的合成、晶体结构的确定及储氢反应能等方面的研究很多[2-10], 尤其是对Li2NH的晶体结构的研究不论实验[2,3,5]还是理论[4,6,8-10]仍存在争议. Ohoyama等[2]通过中子粉末衍射实验得出Li2NH是面心立方结构, 空间群为$F\bar 43m$, 4个H原子随机占据16e位置; Noritake等[3]利用同步加速器–射线衍射方法测到Li2NH晶体是反萤石型结构, 空间群为$Fm\bar 3m$, 晶格常数a = 5.0742Å, H原子随机占据N原子周围48 h位置. Kojima和Kawai[7]从实验上测定了Li2NH晶体的光谱,—NH基的特征谱位于3180和3250 cm–1. Mueller和Ceder[6]以及Magyari-Köpe等[8]分别基于密度泛函理论报道了Li2NH在温度T = 0K下具有最低的能量结构, 呈六角$Pnma$结构, H原子占位4c位置. Herbst 和Hector[4]使用基于密度泛函理论的VASP软件包, 分别应用LDA和GGA研究了LiNH2和Li2NH的晶体结构和电子结构, 并计算了反应(1)式的反应焓. Song和Guo[9]对Li-N-H体系相关的一系列化合物的生成焓和键能等进行了计算, 从理论上证明氨中间体机理的正确性, 并且指出LiNH2中两个氢的脱附能量不同. 最近利用第一性原理对Li2NH体系的研究有很多[11-18], 如Rajeswarapalanichamy等[16]基于密度泛函理论的第一性原理研究了Li-N-H体系的基态性质, 如晶格常数、体变模量和价电子密度等, 同时报道了高压下的结构相变; Crivello等[17]通过VASP软件包, 广义梯度近似(general gradient approximation, GGA) 泛函计算了Li-N-H相关化合物的晶体结构参数, 分析了Li2NH的电子结构、电荷密度和生成焓; Wang等[18]利用第一性原理报道了Li-N-H体系的电子结构和化学键, 并给出了化合物的生成焓和反应(2)的反应焓. 但是迄今为止, 关于Li2NH晶体的晶格动力学和热力学性质的研究还没有相关报道. 因此, 本文利用密度泛函理论(DFT)的第一性原理方法对Li2NH晶体的电子结构、晶格动力学和热力学性质进行报道, 得到布里渊区中心Γ的声子频率、声子色散曲线和相应的态密度, 最后在简谐近似内得到了Li2NH的热力学性质, 其中包括晶格振动对Helmholtz自由能、内能、熵和热容的贡献, 并和已有的实验及理论数据进行比较分析.

    • Li2NH公认的晶体结构模型有三种, 第一种是文献[3]中给出的Li2NH常温常压下为立方反萤石结构(空间群为$Fm\bar 3m$, NO.225), 晶格常数a = 5.0742Å, Li原子占据8c(0.25,0.25,0.25)Wyckoff 格位, N原子占据4a(0,0,0)Wyckoff 格位, H原子占据48h(0.11,0.11,0.0)Wyckoff 格位, 而且H原子随机占据N原子周围48h位置. 第二种为Ohoyama等[2]利用中子粉末衍射实验确定的Li2NH晶体室温下为立方反萤石结构(空间群为$F\bar 43m,$ NO.216), 晶格常数a = 5.0769Å, Li原子占据4c(0.25,0.25,0.25)和4d(0.75,0.75,0.75和)Wyckoff 格位, N原子占据4a(0,0,0)Wyckoff 格位, 而H原子倾向占据16e(0.093,0.093,0.093)Wyckoff 格位, H原子随机占据16e位置. 第三种是文献[8]基于密度泛函理论报道的Li2NH在0K下呈正交晶系(空间群为$Pnma$, NO.62), 为最低的能量结构, 晶格常数a = 7.733Å, b = 3.6Å, c = 4.872Å, 两个Li原子分别占据4c(0.102,0.25,0.464)和(0.201,0.25,0.914)Wyckoff 格位, N原子占据4c(0.354,0.25,0.218)Wyckoff 格位, 而H原子占位4c(0.424,0.25,0.399)Wyckoff 格位. 三种Li2NH晶体的晶胞结构如图1所示.

      图  1  Li2NH的三种晶体结构图 (a) Li2NH ($F\bar 43m$); (b) Li2NH ($Fm\bar 3m$); (c) Li2NH ($Pnma$)

      Figure 1.  Crystal structures of Li2NH: (a) Li2NH ($F\bar 43m$); (b) Li2NH ($Fm\bar 3m$); (c) Li2NH ($Pnma$).

    • 文中所有计算采用基于密度泛函理论平面波赝势方法的 ABINIT 程序包[19], 赝势选用Troullier-Martins模守恒赝势[20], 在局域密度近似(LDA)下, 采用Teter参数化的Ceperley-Alder型交换关联势, 电子波函数由平面波基组展开. 经收敛测试, 平面波截断动能选取40 Hartree, 布里渊区积分选择4 × 4 × 4的Monkhorst-Pack[21]网格, 计算结果具有很好的收敛性. 通过BFGS方法来优化原子位置和晶格常数, 原子受力收敛标准设置5.0d-4Ha/a.u.. 声子频率计算时, $q$点选择在倒易空间6 × 6 × 6网格上, 在线性响应的密度泛函微扰理论[22,23]框架下计算这些$q$点的动力学矩阵, 经傅里叶变换得到实空间的力常数矩阵, 再利用快速傅里叶变换方法就可求出倒易空间中任意$q$点的动力学矩阵, 对动力学矩阵进行对角化即可得到声子频率[24]. 根据声子频率的计算可得到声子态密度, 利用声子态密度和简谐近似计算热力学性质.

    • 首先对图1中的三种结构进行几何优化, 优化后的晶格常数见表1. 为了验证计算结果的可靠性, 把相应的实验结果和理论计算数据也列在表1中. 从表中可以看出, 三种结构的晶格参数和实验及理论结果符合得非常好, 尤其是空间群为$Pnma$的Li2NH的优化后晶格参数a = 7.781Å, b = 3.623Å, c = 4.902Å,和Magyari-Köpe等[8]理论预测的结果符合得最好, 误差仅在1%以内.利用优化好的晶格参数计算了三种结构的总能量, 结果见表1, 结果发现空间群为$Pnma$的Li2NH总能最低, 相对其他结构来说较稳定, 这与文献[8]预测给出的最低能结构是T=0K下H原子占据4c格位的正交晶系(空间群为$Pnma$)一致, 并且Wolverton等[14]也指出T = 0K时具有正交晶系$Pnma$结构的Li2NH总能较低, 同时Crivello等[17]报道了0K下Li2NH的正交$ Pnma $结构比立方$ F\bar 43m$结构稳定, 主要原因是$Pnma$结构中N-H层的排列分布导致了其能量最低. 上述这些结果充分说明本文计算方法和参数的选择对于体系来说是合理的. 下面有关电子结构的计算是在这个优化结构基础上进行的.

      Space groupLattice parameter Etot/eV
      PresentExperimentalOther CalculatedPresent
      $F\bar 43m$a = 5.6249a = 5.0769[2]a = 5.159[17],5.223[17],5.074[17],5.649[16],5.1076[18],5.2968[9]–1920.4448
      $Pnma$a = 7.781
      b = 3.623
      c = 4.902
      a = 7.733[8],7.704[16],7.742[17],7.753[17],7.775[17]
      b = 3.60[8],3.75[16],3.604[17],3.609[17],3.618[17]
      c = 4.872[8],5.074[16],4.883[17],4.890[17],4.881[17]
      –2702.3588
      $Fm\bar 3m$a = 5.0696a = 5.0742[3]a = 5.047[3],5.045[16]–1913.6655

      表 1  Li2NH晶体结构参数

      Table 1.  Lattice parameters of Li2NH.

    • 电子结构一般包括能带结构和态密度, 通过计算得到Li2NH($Pnma$)沿布里渊区高对称方向的能带结构如图2所示, 主要由低价带、价带与导带三部分组成, 最低能带–42 eV附近的能带在态密度图中未画出, 其中虚线代表费米能级. 导带最低点和价带最高点都位于Γ点, 表明Li2NH为直接带隙材料, 导带和价带间的禁带宽度约为2.0 eV, 这和理论得到的带隙2.65 eV[9], 2.8 eV[25]和约3 eV[4]基本符合一致. 总而言之, 本文理论计算的带隙值相对偏小, 一般认为是LDA忽略了交换关联作用的非局域性导致的. 为了克服LDA所带来的带隙估计过小问题, 很多科学家进行了多方面的努力. 最有代表性的是GW近似方法对半导体带隙进行修正[26], 可以得到大致与实验相符合的结果, 这也是课题组接下来要进行的研究工作. 不过文中LDA的计算结果并不影响对能带结构的分析[27].

      图  2  Li2NH的能带结构

      Figure 2.  Band structure of Li2NH.

      态密度对于分析材料中的原子成键和材料特性有重要的意义. 图3给出了Li2NH的总态密度和分态密度图. 从图中可以看出, Li2NH的价带主要分为两个区, –14 eV附近的低能段和–5.5—0 eV高能段, –5.5 eV以上的价带很显然被分成两部分, 其能隙约为1.5 eV, 与文献[28]给出的态密度结果一致; Li2NH的价带在低能段存在一个很高的成键峰,主要是由N s和H s态电子贡献, 并有少量的Li p态电子的贡献. 然而, 在–5.5—–2.5 eV的高能段主要是H s和N p杂化贡献的, 说明N原子和H原子有明显的共振, 从–2.5 eV到费米能级0 eV范围内形成的成键峰, 主要是Li原子和N原子的相互作用贡献的, 比H s和N p杂化形成的弱峰还要弱一些[17]. 进一步分析知, 价带中N和H原子的重叠比Li和N原子的重叠相对多一些, 表明Li2NH内N—H键有强烈的共价键, Li—N的共价键较弱[18].

      图  3  Li2NH的总态密度和分波态密度

      Figure 3.  Total and partial density of states of Li2NH(states/eV)

    • 晶格振动的理论已成为研究材料的宏观性质和微观过程的重要理论基础, 研究晶格振动性质不仅可以了解材料的热力学性质, 还可以深入理解电子结构、光学性质等. Li2NH三种结构的声子色散曲线分别如图4图6所示. 由图可见, Li2NH($F\bar 43m$)和Li2NH($Fm\bar 3m$)的声子色散曲线中都存在虚频, 说明这两种结构在动力学上是不稳定的, 而Li2NH($Pnma$)声子色散曲线没有虚频, 说明该结构的相对稳定性, 这与我们关于基态能量计算的结果相一致, 再一次证明了Li2NH($Pnma$)为较稳定的结构. 从图4中可以看到, Li2NH的频率主要集中在0—700 cm–1和3100—3200 cm–1两段, 振动分析其结构, 低频对应于N—H键的弯曲振动, 而高频对应于N—H键的伸缩振动, 该段为—NH基的特征谱.计算结果和文献[29]中实验上用原位红外光谱仪测量Li2NH在高达380 ℃以上的红外吸收峰3150 cm–1一致. 并且和文献[30]中利用B3LYP/6-31G*方法的频率修正因子0.9613修正后得到的N—H伸缩振动位于3198.18—3374.77 cm–1符合得很好. 同时发现本文理论计算和文献[7]从实验上测定了Li2NH晶体的光谱, —NH基的特征谱位于3180和3250 cm–1, 几乎完全一致. 进一步说明了—NH基在晶体中保持其完整性[30]. 可见, 我们的计算结果和实验及理论数据符合很好.

      图  4  Li2NH($Pnma$)的声子色散曲线

      Figure 4.  The phonon dispersion curves of Li2NH ($Pnma$).

      图  5  Li2NH($F\bar 43m$)的声子色散曲线

      Figure 5.  The phonon dispersion curves of Li2NH ($F\bar 43m$).

      图  6  Li2NH($Fm\bar 3m$)的声子色散曲线

      Figure 6.  The phonon dispersion curves of Li2NH ($Fm\bar 3m$)

      下面详细分析Li2NH($Pnma$)布里渊区中心Γ的光学模频率. 根据晶格振动理论, 由于其每个原胞中包含16个原子, 故在每个q点有3个声学振动模(Γaco)和45个光学振动模(Γopt). 通过施加声学模求和规则[31]可以使三个声学模的频率严格等于零而满足平移对称性要求. 由群理论分析, Γ点的光学模(Γopt)不可约表示为 Γopt=4 B2u+7 B1u+4 B3g+8 B2g+7 B3u+4 Au+4 B1g+7 Ag. 其中B2u, B1uB3u模是红外活性模, B3g, B2g, B1gAg为拉曼活性模, Au模既不属于拉曼活性模也不属于红外活性模, 称之为Silent 模. 计算得到的Li2NH($Pnma$)布里渊区中心Γ点的光学模频率见表2, 低频范围在100—700 cm–1, 而高频分布在3100—3200 cm–1, 中间存在一个很宽的声子禁带. 经过大量调研发现, 这是第一次对Li2NH ($Pnma$)布里渊区中心Γ点的光学模频率进行的研究, 故实验数据几乎没有. 希望本文工作对实验结果起到一定的预测作用.

      B2uB1uB3gB2gB3uAuB1gAg
      123.6137.9160161.5168198.5222223.5
      310195.8380232.0289243.0359287
      408318432300350.5446435335
      644450660350.5385625654380
      550412620542
      630580668598
      312264031233132
      3134

      表 2  Li2NH($Pnma$)在布里渊区中心Γ点的光学模声子频率(单位:cm–1)

      Table 2.  Phonon frequencies (unit:cm–1) at the Γ point of Li2NH.

    • 了解材料的热力学性质对于研究其稳定性和化学反应是至关重要的. 基于图7得到的声子态密度, 按照晶格振动理论很容易得到Li2NH的热力学性质. 在简谐近似内, 晶格振动对Helmholtz自由能ΔF的贡献为[32]

      图  7  Li2NH($Pnma$)的声子态密度

      Figure 7.  The phonon density of states of Li2NH ($Pnma$).

      $ \Delta F=3 n N k_{B} T \int_{0}^{\omega_{\max }} \ln \left\{2 \sin \frac{\hbar \omega}{2 k_{B} T}\right\} g(\omega) d \omega, $

      晶格振动对内能的贡献ΔE[32]

      $ \Delta E = 3nN\frac{\hbar }{2}\int_0^{{\omega _{\max }}} {\omega \coth \left( {\frac{{\hbar \omega }}{{2{k_{\rm{B}}}T}}} \right)} g\left( \omega \right){\rm{d}}\omega, $

      晶格振动对熵S的贡献为[32]

      $\begin{split} S =\; & 3nN{k_{\rm{B}}}\int_0^{{\omega _{\max }}}\left[ \frac{{\hbar \omega }}{{2{k_{\rm{B}}}T}}\coth \left( {\frac{{\hbar \omega }}{{2{k_{\rm{B}}}T}}} \right) \right.\\ &\left.- {\rm{ln}}\{ 2\sin \frac{{\hbar \omega }}{{2{k_{\rm{B}}}T}}\} \right]g\left( \omega \right){\rm{d}}\omega,\end{split}$

      晶格振动对热容的贡献为[32]

      ${C_V} = 3nN{k_{\rm{B}}}{\int_0^{{\omega _{\max }}} {\left( {\frac{{\hbar \omega }}{{2{k_{\rm{B}}}T}}} \right)} ^2}\csc {h^2}\left( {\frac{{\hbar \omega }}{{2{k_{\rm{B}}}T}}} \right)g(\omega ){\rm{d}}\omega,$

      其中kB是玻尔兹曼常数, n是原胞中原子的数目, N是晶胞的数目, $\omega $为声子频率, ${\omega _{\max }}$是声子频率最大值, $g\left( \omega \right)$是声子态密度, 满足归一化条件$\displaystyle\int_0^{{\omega _{\max }}} {g\left( \omega \right)} {\rm{d}}\omega = 1$. 图8给出了Li2NH晶体晶格振动对自由能的贡献、内能、熵及比热随温度的变化, 可以看到体系熵S、内能ΔE随温度增加呈非线性增加, 自由能ΔF下降, 比热容CV高温时达到稳定. 从图8(b)中可以看到, 在温度300 K以上, 内能ΔE随温度的升高几乎呈现线性关系,趋于kBT行为. 从图8(d)中可以看到, 低温时比热容CV几乎正比于T的三次方(见红色▼标注点), 即Debye的T3定律, 高温时接近于一常数$3R \times $16 = 398.88 J/mol-c.K(见虚线), 即经典的Dulong-Petit定律, 说明本文的计算是合理可靠的.

      图  8  Li2NH的热力学性质 (a)晶格振动对自由能的贡献ΔF; (b)内能ΔE随温度的变化; (c)熵S随温度的变化; (d)热容Cv随温度的变化 (▼表示Debye的T3定律)

      Figure 8.  Thermodynamic properties for Li2NH: (a) The phonon contribution to the free energies ΔF of Li2NH; (b)the internal energies ΔE verus temperature; (c)the entropy S verus temperature; (d) the constant-volume specific heats Cv verus temperature(▼ meansT3 law of Debye)

      表3给出了计算的S和热容CV数据, 同时以298.15 K时的焓H298为标准, 计算了Li2NH在不同温度下ΔH(H-H298)的值, 结果也列在表3中. 我们知道晶格振动并不是严格简谐的,存在热膨胀,实际的热膨胀是原子之间的非谐作用所引起的, 导致简谐近似内计算得到的S和ΔH数据在高温下或多或少会存在一点偏差. 虽然我们没有找到实验的热容CV值, 但是基于密度泛函微扰理论计算得到热力学函数是合理的, 理论模拟能够得到材料的热力学性质, 希望为实验上提供参考.

      Temprature/KS/J·(mol-c·K)–1CV/J·(mol-c·K)–1ΔH(HH298)/KJ·mol–1
      298248.9273.80
      300250.2274.92.9
      500405.4327.155.1
      700519.0347.075.0
      900607.7358.586.5
      1100680.5366.694.2
      1300742.3372.8100.8
      1500796.0377.5105.5

      表 3  Li2NH的热力学函数

      Table 3.  Thermodynamic functions for Li2NH.

    • 本文基于密度泛函理论的第一性原理方法计算了Li2NH的电子结构、晶格动力学和热力学性质. 优化后的晶格常数和实验及理论数据符合很好, 其中Li2NH($Pnma$)的计算结果与理论值相差最小, 该体系相对最稳定.能带结构和态密度显示Li2NH属于直接带隙材料, 带隙约2.0 eV, N—H键有强烈的共价性. 基于线性响应的密度泛函理论计算了Li2NH的声子色散曲线和声子态密度, 结果发现只有Li2NH($Pnma$)的声子色散曲线没有虚频, 动力学性能相对最稳定, —NH基的特征谱结果和先前实验及理论数据符合得很好. 同时对布里渊区中心Γ点的光学模频率进行了详细分析. 根据声子态密度计算了Li2NH的热力学性质, 结果发现熵、内能随温度增加呈非线性增加, Helmholtz自由能下降, 热容高温时趋于稳定. 遗憾的是Li2NH的声子频率数据未发现, 希望本文的计算结果能为Li-N-H体系的实验研究提供进一步的指导.

参考文献 (32)

目录

    /

    返回文章
    返回