钙钛矿结构ZrBeO3稳定性的第一性原理研究

温新竹 彭玉颜 刘明真

引用本文:
Citation:

钙钛矿结构ZrBeO3稳定性的第一性原理研究

    作者简介: 温新竹(1982-),男,副教授,主要从事材料学研究. E-mail: 64412900@qq.com;
  • 中图分类号: O521.2; TM911.4

First-Principles Study on Structural Stability of Perovskite ZrBeO3

  • CLC number: O521.2; TM911.4

  • 摘要: 基于密度泛函理论构建了钙钛矿结构ZrBeO3晶体模型,计算了该晶体模型结合能,表明了该构型热力学稳定性;计算出该结构在不同压力下的弹性常数,并据此计算了ZrBeO3的体积模量、剪切模量、杨氏模量、泊松比、BH/GH(体模量/剪切模量)等参数,结果表明该材料具有机械稳定性,随着等静压力增加,材料由脆性向韧性转变;计算了零压下ZrBeO3的硬度,为34.5 GPa,表明该结构晶体应为超硬材料;计算了ZrBeO3的声子能谱,结果表明ZrBeO3在低温零压下热动力学不稳定,为此分析比较了不同压力下的声子能谱、不同原子轨道及化学键布居值,研究表明随着压力增加,Be原子sp杂化后形成的Be-O共价键成分增强、Zr-O键离子键成分增强,晶格动力学趋于稳定。
  • 图 1  ZrBeO3晶格常数a随总能E变化曲线

    Figure 1.  Lattice constants (a) vs. total energy (E) of ZrBeO3

    图 2  ZrBeO3晶体模型

    Figure 2.  ZrBeO3 crystal model

    图 3  零压下的ZrBeO3声子能谱

    Figure 3.  ZrBeO3 phonon spectra at zero pressure

    图 4  不同压力下的ZrBeO3声子能谱

    Figure 4.  Phonon spectra of ZrBeO3 at different pressures

    图 5  不同压力下ZrBeO3电子总态密度及各原子分波态密度图

    Figure 5.  Total density of state of ZrBeO3 electrons and partial density of states of each atom at different pressures

    图 A1  ZrBeO3的电子能带结构

    Figure A1.  Electronic band structure of ZrBeO3

    图 A2  ZrBeO3的吸收谱

    Figure A2.  Optical absorption spectrum of ZrBeO3

    图 A3  ZrBeO3的反射谱

    Figure A3.  Optical reflectivity spectrum of ZrBeO3

    图 A4  ZrBeO3介电函数谱线

    Figure A4.  Dielectric function spectrum of ZrBeO3

    表 1  ZrBeO3模型晶格参数

    Table 1.  Lattice parameters of ZrBeO3 model

    CompoundLattice/nmWykoff coordinatesVolume/(10–3 nm3)Bond Length/nmBond population
    ZrBeO3 cubic0.346 1Zr(0.00, 0.00, 0.00)a
    Be(0.50, 0.50, 0.50)a
    O(0.50, 0.50, 0.00)a
    41.447 3Be-O 0.173 0
    Zr-O 0.244 7
    Be-O (3) 0.66
    Zr-O (3) 0.61
    下载: 导出CSV

    表 2  不同压力下的弹性常数、体积模量、剪切模量、杨氏模量、泊松比、BH/GH

    Table 2.  Elastic constants, bulk modulus, shear modulus, Young’s modulus, Poisson’s ratio, BH/GH under different pressures

    p/GPaC11C12C44BHGHE$v$BH/GH
    0 422.45149.33171.48 240.37157.55 385.860.2321.53
    30 582.04212.84222.98 335.91206.74 514.650.2451.62
    50 775.37315.82306.78 469.01273.25 686.430.2561.72
    1001 021.29401.25353.07 607.94335.17 850.860.2671.81
    1501 281.55513.85432.03 769.75412.761 049.020.2731.86
    2001 515.13614.39500.17 914.63479.621 224.670.2771.91
    2501 778.26731.79580.321 080.61556.761 425.470.2801.94
    3002 256.08940.80723.011 379.23696.111 787.600.2841.98
    下载: 导出CSV

    表 3  不同压力下原子轨道和化学键的布居值分布

    Table 3.  Atomic orbital and chemical bond population distribution at different pressures

    p/
    GPa
    AtomicspdBond Be-OBond Zr-O
    PopulationBond length/nmPopulationBond length/nm
    0Be/Zr/O(3)2.34/2.26/1.821.11/6.31/4.890/1.84/00.630.173 00.620.244 7
    50Be/Zr/O(3)2.31/2.18/1.791.24/6.24/4.920/1.91/00.690.162 50.530.229 8
    100Be/Zr/O(3)2.29/2.12/1.771.31/6.21/4.930/1.97/00.720.157 70.470.223 0
    150Be/Zr/O(3)2.27/2.08/1.761.37/6.18/4.940/2.01/00.750.154 20.400.218 0
    200Be/Zr/O(3)2.26/2.04/1.751.42/6.16/4.940/2.05/00.770.151 30.330.214 0
    下载: 导出CSV
  • [1] IKEGAMI K, LU M, OHON R, et al. Nonlinear electrical properties of thin films of a light-emitting perovskite type oxide Pr0.002(Ca0.6Sr0.4)0.997TiO3 [J]. Procedia Engineering, 2012, 36: 388–395. doi: 10.1016/j.proeng.2012.03.057
    [2] TAKASHIMA H, SHIMADA K, MIURA N, et al. Low-driving-voltage electroluminescence in perovskite films [J]. Advanced Materials, 2009, 21(36): 3699–3702. doi: 10.1002/adma.200900524
    [3] SAHA S, SINHA T P, MOOKERJEE A. Electric structure, chemical bonding, and optical properties of paraelectric BaTiO3 [J]. Physical Review B, 2000, 62(13): 699–702.
    [4] 赵国栋, 杨亚利, 任伟. 钙钛矿型氧化物非常规铁电研究进展 [J]. 物理学报, 2018, 67(15): 60–72.
    ZHAO G D, YANG Y L, REN W. Progress in unconventional ferroelectricity of perovskite-type oxides [J]. Acta Physica Sinica, 2018, 67(15): 60–72.
    [5] PARK S Y, KUMAR A, RABE K. Carbon fibers from polyacrylonitrile/cellulose [J]. 2016 APS Meeting, 2016, 6(17): 160–172.
    [6] FENNIE C J. Ferroelectrically induced weak ferromagnetism by design [J]. Physical Review Letters, 2008, 100(16): 167203. doi: 10.1103/PhysRevLett.100.167203
    [7] DOLAN D H, AO T. Cubic zirconia as a dynamic compression window [J]. Applied Physics Letters, 2008, 93(2): 021908. doi: 10.1063/1.2957996
    [8] RESTANI R, MARTIN M, KIVEL N, et al. Analytical investigations of irradiated inert matrix fuel [J]. Journal of Nuclear Materials, 2009, 385(2): 435–442. doi: 10.1016/j.jnucmat.2008.12.030
    [9] CONRADSON S D, DEGUELDRE C A, ESPINOSA-FALLER F J, et al. Complex behavior in quaternary zirconias for inert matrix fuel: what do these materials look like at the nanometer scale? [J]. Progress in Nuclear Energy, 2001, 38(3/4): 221–230.
    [10] WANG S J, ONG C K, XU S Y, et al. Crystalline zirconia oxide on silicon as alternative gate dielectrics [J]. Applied Physics Letters, 2001, 78(11): 1604–1606. doi: 10.1063/1.1354161
    [11] LIN Y S, PUTHENKOVILAKAM R, CHANG J P, et al. Interfacial properties of ZrO2 on silicon [J]. Journal of Applied Physics, 2003, 93(10): 5945–5952. doi: 10.1063/1.1563844
    [12] WILK G D, WALLACE R M, ANTHONY J M. High-κ gate dielectrics: current status and materials properties considerations [J]. Journal of Applied Physics, 2001, 89(10): 5243–5275. doi: 10.1063/1.1361065
    [13] MCEVOY A. Thin SOFC electrolytes and their interfacesâ: a near-term research strategy [J]. Solid State Ionics, 2000, 132(3/4): 159–165.
    [14] BILIĆ A T, GALE J D. Ground state structure of BaZrO3: a comparative first-principles study [J]. Physical Review B, 2009, 79(17): 174107. doi: 10.1103/PhysRevB.79.174107
    [15] ZHANG H W, FU X Y, NIU S Y, et al. Synthesis and photoluminescence properties of Eu3+-doped AZrO3 (A = Ca, Sr, Ba) perovskite [J]. Journal of Alloys and Compounds, 2008, 459(1/2): 103–106.
    [16] DUBEY V, TIWARI N. Structural and optical analysis on europium doped AZrO3 (A = Ba, Ca, Sr) phosphor for display devices application [C]//Bikaner, India. Author(s), 2016, 1728(1): 15-32.
    [17] RANDALL C A, BHALLA A S, SHROUT T R, et al. Classification and consequences of complex lead perovskite ferroelectrics with regard to B-site cation order [J]. Journal of Materials Research, 1990, 5(4): 829–834. doi: 10.1557/JMR.1990.0829
    [18] BRADHA M, HUSSAIN S, CHAKRAVARTY S, et al. Total conductivity in Sc-doped LaTiO3+δ perovskites [J]. Ionics, 2014, 20(9): 1343–1350. doi: 10.1007/s11581-014-1216-y
    [19] GEPPERT B, GROENEVELD D, LOBODA V, et al. Finite-element simulations of a thermoelectric generator and their experimental validation [J]. Energy Harvesting and Systems, 2015, 2(1/2): 97–103.
    [20] MUHAMMAD I D, AWANG M, MAMAT O, et al. First-principles calculations of the structural, mechanical and thermodynamics properties of cubic zirconia [J]. World Journal of Nano Science and Engineering, 2014, 4(2): 97–103. doi: 10.4236/wjnse.2014.42013
    [21] BAERENDS E J. Perspective on “Self-consistent equations including exchange and correlation effects” [J]. Theoretical Chemistry Accounts: Theory, Computation, and Modeling (Theoretica Chimica Acta), 2000, 103(3/4): 265–269.
    [22] SEGALL M D, LINDAN P J D, PROBERT M J, et al. First-principles simulation: ideas, illustrations and the CASTEP code [J]. Journal of Physics: Condensed Matter, 2002, 14(11): 2717–2744. doi: 10.1088/0953-8984/14/11/301
    [23] 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
    [24] JIA X F, HOU Q Y, XU Z C, et al. Effect of Ce doping on the magnetic and optical properties of ZnO by the first principle [J]. Journal of Magnetism and Magnetic Materials, 2018, 465: 128–135. doi: 10.1016/j.jmmm.2018.05.037
    [25] FISCHER T H, ALMLOF J. General methods for geometry and wave function optimization [J]. The Journal of Physical Chemistry, 1992, 96(24): 9768–9774. doi: 10.1021/j100203a036
    [26] DING Y C, CHEN M, WU W J. Phase stability, elasticity, hardness and the minimum thermal conductivity of Si2N2O polymorphs from first principles calculations [J]. Physica B: Condensed Matter, 2014, 449: 236–245. doi: 10.1016/j.physb.2014.05.042
    [27] GONZE X, BEUKEN J M, CARACAS R, et al. First-principles computation of material properties: the ABINIT software project [J]. Computational Materials Science, 2002, 25(3): 478–492. doi: 10.1016/S0927-0256(02)00325-7
    [28] DING J F, LI X M, CUI L L, et al. Electronic and optical properties of anion-doped c-ZrO2 from first-principles calculations [J]. Journal of Central South University, 2014, 21(7): 2584–2589. doi: 10.1007/s11771-014-2216-9
    [29] BORN M, HUANG K, LAX M. Dynamical theory of crystal lattices [J]. American Journal of Physics, 1955, 23(7): 474.
    [30] WU Z J, HAO X F, LIU X J, et al. Structures and elastic properties of OsN2 investigated via first-principles density functional calculations [J]. Physical Review B, 2007, 75(5): 054115. doi: 10.1103/PhysRevB.75.054115
    [31] ZHAO J J, WINEY J M, GUPTA Y M. First-principles calculations of second-and third-order elastic constants for single crystals of arbitrary symmetry [J]. Physical Review B, 2007, 75(9): 094105. doi: 10.1103/PhysRevB.75.094105
    [32] GAO F M. Theoretical model of intrinsic hardness [J]. Physical Review B, 2006, 73(13): 132104. doi: 10.1103/PhysRevB.73.132104
    [33] FAN C Z, ZENG S Y, LI L X, et al. Potential superhard osmium dinitride with fluorite and pyrite structure: first-principles calculations [J]. Physical Review B, 2006, 74(12): 125118. doi: 10.1103/PhysRevB.74.125118
    [34] 丁迎春, 肖冰. 一种超硬新材料BeP2N4的电子结构和力学性质及本征硬度 [J]. 物理化学学报, 2011, 27(7): 1621–1632. doi: 10.3866/PKU.WHXB20110730
    DING Y C, XIAO B. Electronic structure, mechanical properties and intrinsic hardness of a new superhard material BeP2N4 [J]. Acta Physico-Chimica Sinica, 2011, 27(7): 1621–1632. doi: 10.3866/PKU.WHXB20110730
    [35] SHARMA A D, SINHA M M. Lattice dynamics of protonic conductors AZrO3 (A = Ba, Sr & Pb): a comparative study [J]. Advanced Materials Research, 2013, 685: 191–194. doi: 10.4028/www.scientific.net/AMR.685.191
    [36] 刘哲, 李辉, 赵鹏. Ti5Al2C3与Ti2AlC、Ti3AlC2结构、弹性和电子性质的第一性原理对比研究 [J]. 人工晶体学报, 2019, 48(5): 834–839. doi: 10.3969/j.issn.1000-985X.2019.05.011
    LIU Z, LI H, ZHAO P. A first-principles comparative study of the structure, elasticity and electronic properties of Ti5Al2C3, Ti2AlC and Ti3AlC2 [J]. Journal of Artificial Lenses, 2019, 48(5): 834–839. doi: 10.3969/j.issn.1000-985X.2019.05.011
    [37] LAI J, JIA X, WANG D. Thermodynamically stable whilst kinetically labile coordination bonds lead to strong and tough self-healing polymers [J]. Nature Communications, 2019, 10(1): 155–167. doi: 10.1038/s41467-018-07819-1
    [38] GULL E, PARCOLLET O, MILLIS A J. Superconductivity and the pseudogap in the two-dimensional Hubbard model [J]. Physical Review Letters, 2013, 110(21): 256–298.
  • [1] 杨龙星刘雷刘红易丽顾小雨 . 高压下石榴子石结构和弹性的第一性原理研究. 高压物理学报, 2019, 33(6): 060104-1-060104-11. doi: 10.11858/gywlxb.20190785
    [2] 谭心贾亦超刘学杰 . 外压下ZrN结构相变的第一性原理研究. 高压物理学报, 2014, 28(2): 168-174. doi: 10.11858/gywlxb.2014.02.006
    [3] 韩林马麦宁徐志双周晓亚 . 辉石结构与相变的第一性原理研究. 高压物理学报, 2017, 31(2): 125-134. doi: 10.11858/gywlxb.2017.02.004
    [4] 丁迎春刘海军蒋孟衡陈敏陈勇明 . 高压下BeP2N4结构相变和电子结构的第一性原理计算. 高压物理学报, 2012, 26(6): 674-680. doi: 10.11858/gywlxb.2012.06.012
    [5] 宋功保梁敬魁饶光辉 . 固相线下Pr1-xBa2-yCax+yCu3O7体系的相关系和晶体结构. 高压物理学报, 2004, 18(4): 309-314 . doi: 10.11858/gywlxb.2004.04.004
    [6] 赵廷河阎学伟崔硕景陈久华刘立君牛伟赵伟 . 宝石翡翠的合成、热行为和热稳定性的研究. 高压物理学报, 1992, 6(4): 291-296 . doi: 10.11858/gywlxb.1992.04.008
    [7] 刘思远缪宇马雪姣李鑫高文泉程宇衡刘艳辉 . 基于第一性原理计算IrSb压力相变. 高压物理学报, 2019, 33(5): 052203-1-052203-7. doi: 10.11858/gywlxb.20190716
    [8] 李胜旨刘锦超杨向东蒋德琼 . Al、N共掺杂实现p-型ZnS的第一性原理研究. 高压物理学报, 2011, 25(6): 519-525. doi: 10.11858/gywlxb.2011.06.007
    [9] 李胜旨刘锦超杨向东郭艳锋许海全 . Mn、Fe掺杂ZnS的第一性原理计算. 高压物理学报, 2010, 24(6): 449-454 . doi: 10.11858/gywlxb.2010.06.008
    [10] 郝军华吴志强王铮金庆华李宝会丁大同 . 高压下SiO2的第一性原理计算. 高压物理学报, 2010, 24(4): 260-266 . doi: 10.11858/gywlxb.2010.04.004
    [11] 孙磊罗坤刘兵韩俏怡王小雨梁子太赵智胜 . 金属性硅同素异形体的第一性原理研究. 高压物理学报, 2019, 33(2): 020103-1-020103-8. doi: 10.11858/gywlxb.20190705
    [12] 宋功保韩苍穹梁敬魁饶光辉 . 基于组分Ca2Pr2Cu5O10三元固溶区的范围和晶体结构. 高压物理学报, 2004, 18(1): 27-35 . doi: 10.11858/gywlxb.2004.01.006
    [13] 孙莹钟鑫吕健马琰铭 . 氧族氢化物的压致金属化与奇异超导电性. 高压物理学报, 2018, 32(1): 010102-1-010102-10. doi: 10.11858/gywlxb.20170629
    [14] 李鑫马雪姣高文泉刘艳辉 . 高压下Ir2P晶体结构预测与物理性质. 高压物理学报, 2019, 33(1): 011103-1-011103-6. doi: 10.11858/gywlxb.20180645
    [15] 缪宇刘思远马雪姣金哲学 . 高压下MgN8晶体结构理论模拟与物性研究. 高压物理学报, 2020, 34(1): 011102-1-011102-8. doi: 10.11858/gywlxb.20190818
    [16] 张乐婷赵宇宏孙远洋邓世杰吉如意韩培德 . 压力作用下Mg2X(X=Si, Ge)相热力学性质的第一性原理研究. 高压物理学报, 2018, 32(3): 032201-1-032201-8. doi: 10.11858/gywlxb.20170630
    [17] 刘克伟于杰周晓龙胡明钰占建翔 . 高压下P42/mnmPnnm晶体结构SnO2密度泛函研究. 高压物理学报, 2017, 31(1): 81-88. doi: 10.11858/gywlxb.2017.01.012
    [18] 李全马琰铭 . 典型双原子分子晶体的高压解离和单原子相. 高压物理学报, 2013, 27(3): 313-324. doi: 10.11858/gywlxb.2013.03.001
    [19] 柳雷李晓东李延春唐玲云刘景毕延 . NiO的高压结构和等温状态方程研究. 高压物理学报, 2009, 23(3): 209-214 . doi: 10.11858/gywlxb.2009.03.008
    [20] 顾惠成陈良辰褚君浩鲍忠兴刘克岳王金义郑康立 . Hg1-xCdxTe在高压下的结构、状态方程与相变. 高压物理学报, 2000, 14(4): 298-301 . doi: 10.11858/gywlxb.2000.04.011
  • 2019-0802_A_origin.rar
  • 加载中
图(9)表(3)
计量
  • 文章访问数:  267
  • 阅读全文浏览量:  250
  • PDF下载量:  3
出版历程
  • 收稿日期:  2019-07-02
  • 录用日期:  2019-07-25
  • 网络出版日期:  2020-01-17
  • 刊出日期:  2020-02-01

钙钛矿结构ZrBeO3稳定性的第一性原理研究

    作者简介:温新竹(1982-),男,副教授,主要从事材料学研究. E-mail: 64412900@qq.com
  • 仰恩大学工程技术学院,福建 泉州 362014

摘要: 基于密度泛函理论构建了钙钛矿结构ZrBeO3晶体模型,计算了该晶体模型结合能,表明了该构型热力学稳定性;计算出该结构在不同压力下的弹性常数,并据此计算了ZrBeO3的体积模量、剪切模量、杨氏模量、泊松比、BH/GH(体模量/剪切模量)等参数,结果表明该材料具有机械稳定性,随着等静压力增加,材料由脆性向韧性转变;计算了零压下ZrBeO3的硬度,为34.5 GPa,表明该结构晶体应为超硬材料;计算了ZrBeO3的声子能谱,结果表明ZrBeO3在低温零压下热动力学不稳定,为此分析比较了不同压力下的声子能谱、不同原子轨道及化学键布居值,研究表明随着压力增加,Be原子sp杂化后形成的Be-O共价键成分增强、Zr-O键离子键成分增强,晶格动力学趋于稳定。

English Abstract

  • 钙钛矿氧化物(ABO3)由于具有铁电、铁磁、超导、压电等特性,在光波导、激光倍频、压电传感器、可调谐电容器、高容量存储单元等方面具有广泛的应用[1-4]。Park等[5]发现在电荷有序的短周期LaVO3/SrVO3钙钛矿超晶格中,通过调控钒离子价态可以形成拥有很大极化强度的亚稳态。Fennie[6]在寻求强磁电耦合材料的过程中发现PbTiO3/BiFeO3电磁耦合机制不同。Sonali等[3]计算了顺电BaTiO3的电子结构、化学键、光学性质,揭示出其应用的巨大潜力。常见锆酸盐类物质(AZrO3,A = Ca, Sr, Ba等)一般具有高熔点、高热膨胀系数、低热导率、高化学稳定性、高辐射稳定性,以及优良的离子导电性能,在热障涂层材料[7]、核材料[8]、耐火材料[9]、发光材料[10-12]、高温离子传感器件[13-14]等诸多领域得到广泛应用。目前关于ABO3结构与性能方面的计算预测研究众多,相对于锆基功能材料的广泛应用,ABO3型结构的锆酸盐晶体材料的理论研究较少[15-16],而对于Zr元素与B、O可能构成的两种ZrBeO3、BeZrO3结构几乎没有任何相关报道。本研究基于密度泛函理论(Density functional theory, DFT)的第一性原理,采用模守恒赝势平面波方法,通过对ZrBeO3的化学键Mulliken布居值以及声子谱、声子态密度的计算分析,从电子结构角度出发探讨ZrBeO3稳定性差的本质原因,研究结果对于进一步优化ZrBeO3的稳定性,缩短其实验研究及应用进程,以及开发Zr基多功能材料提供一定的理论参考。

    • ABO3型钙钛矿的结构为简立方晶格,空间群Pm$\overline 3$m(221),晶胞中有5个原子,B占据立方晶格体心位置,A占据8个顶点位置,晶胞6个面心位置由O离子占据,O离子与中心B位离子形成BO6八面体。整个晶体也可看成由BO6八面体共顶点连接而成,各八面体之间的空隙位置由A位离子占据,A、B离子与O离子的配位数分别为12和6,一般而言,若要组成ABO3型钙钛矿结构,A、B、O的离子半径应满足紧密堆积原则

      ${R_{\rm A}} + {R_{\rm O}} = 1.414t({R_{\rm B}} + {R_{\rm O}})$

      式中:RARBRO分别为A、B和O离子的半径;t为容忍因子,一般情况下,t介于0.78~1.05[17]之间时,可组成钙钛矿结构,当t = 1时,体系为理想钙钛矿结构。

      由于Be2+、Zr4+、O2–离子半径分别为0.17、0.72、1.21Å[18-19],若要满足空间几何结构稳定,根据式(1),几种元素只能组合成Zr占据晶胞顶点位置的ZrBeO3钙钛矿结构,其容忍因子为0.989,非常接近理想钙钛矿结构。经过多次模拟和优化,最终建立的ZrBeO3晶胞的晶格常数为0.346 1 nm。

    • 计算在Materials studio 2017软件的CASTEP[20]模块中完成。在密度泛函理论框架下[21],采用广义梯度近似(General gradient approximate, GGA)方法[22-23],其中交换关联函数选择Perdew-Burke-Ernzerhof(PBE)函数[24],几何结构优化采用Broyden-Fletcher-Goldfarb-Shanno (BFGS)算法[25-26],声子能谱计算采用线性响应方法,布里渊区的积分采用以G点为中心6 × 6 × 6的K点设置。考虑过渡金属元素Zr的计算精确性,计算中离子实与价电子之间的相互作用选用模守恒赝势[27]描述,具体的电子组态分别为Be(1s22s2)、Zr(4s24p64d25s2)、O(2s22p4),平面波截止能Ecut为898.0 eV。结构的自洽优化收敛标准设置为:最大位移5.0×10−4 nm,最大内应力收敛标准0.02 GPa,原子间的最大相互作用力0.01 eV/nm,结构的总体能量收敛小于5×10−6 eV/atom。

    • 至目前为止还没有形成稳定的ZrBeO3晶体化合物,因此未检索到相关实验数据,在理论计算方面Materials project平台数据库中有所有元素经过高通量计算的初步结构筛选数据,本研究在了解该平台晶体结构信息的基础上建立晶体模型,明确Be、Zr、O元素在晶胞中各自占位后,去除晶胞的所有对称性,对晶胞常数a、b、c在0.138~0.432 nm(晶胞中原子接近重叠至晶胞中各原子不能成键,取值间隔0.01 nm)大范围内分别取任意值进行结构充分弛豫,发现弛豫后结果均为晶格常数在0.34 nm附近的立方钙钛矿结构,说明钙钛矿构型是该三元化合物中最稳定相。进一步以a = 0.34 nm为中心作钙钛矿构型ZrBeO3的晶格常数a与总能的变化曲线,如图1所示,可得出最低能量点–2 627.83 eV处于晶格常数a = 0.347 nm附近,对此晶格常数建立模型,再次优化后得到最终模型晶格常数为0.346 1 nm的钙钛矿结构,如图2所示,晶胞相关参数列于表1

      图  1  ZrBeO3晶格常数a随总能E变化曲线

      Figure 1.  Lattice constants (a) vs. total energy (E) of ZrBeO3

      CompoundLattice/nmWykoff coordinatesVolume/(10–3 nm3)Bond Length/nmBond population
      ZrBeO3 cubic0.346 1Zr(0.00, 0.00, 0.00)a
      Be(0.50, 0.50, 0.50)a
      O(0.50, 0.50, 0.00)a
      41.447 3Be-O 0.173 0
      Zr-O 0.244 7
      Be-O (3) 0.66
      Zr-O (3) 0.61

      表 1  ZrBeO3模型晶格参数

      Table 1.  Lattice parameters of ZrBeO3 model

      图  2  ZrBeO3晶体模型

      Figure 2.  ZrBeO3 crystal model

      定义ZrBeO3结合能[28]Ef = EtEZrEBe – 3EO,其中Et为结构充分优化后的钙钛矿型ZrBeO3的总能量(–2 627.83 eV),EZrEBeEO分别为Zr、Be、O的单自由原子能量,Ef为–38.27eV(–7.65 eV/atom),说明ZrBeO3钙钛矿构型具有较强的热力学稳定性。

    • 立方晶系独立弹性常数C11C12C44描述了晶体对外加应变的响应刚度,计算了钙钛矿结构ZrBeO3在0~300 GPa不同压力下的弹性常数,如表2所示。结合玻恩-黄稳定性判据[29-30]

      p/GPaC11C12C44BHGHE$v$BH/GH
      0 422.45149.33171.48 240.37157.55 385.860.2321.53
      30 582.04212.84222.98 335.91206.74 514.650.2451.62
      50 775.37315.82306.78 469.01273.25 686.430.2561.72
      1001 021.29401.25353.07 607.94335.17 850.860.2671.81
      1501 281.55513.85432.03 769.75412.761 049.020.2731.86
      2001 515.13614.39500.17 914.63479.621 224.670.2771.91
      2501 778.26731.79580.321 080.61556.761 425.470.2801.94
      3002 256.08940.80723.011 379.23696.111 787.600.2841.98

      表 2  不同压力下的弹性常数、体积模量、剪切模量、杨氏模量、泊松比、BH/GH

      Table 2.  Elastic constants, bulk modulus, shear modulus, Young’s modulus, Poisson’s ratio, BH/GH under different pressures

      ${C_{11}} > 0,\quad {C_{44}} > 0,\quad{C_{11}} > \left| {{C_{12}}} \right|,\quad{C_{11}} + 2{C_{12}} > 0$

      进行判断,所计算的ZrBeO3在不同压力下的钙钛矿结构在力学方面均满足稳定性要求。

      根据Voigt-Reuss-Hill近似,立方晶系体模量BH、剪切模量GH、杨氏模量E、泊松比$v$可以根据式(3)~式(6)计算[31],结果如表2所示。

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

      $G_{\rm H} = \frac{1}{{10}}\left( {{C_{11}} - {C_{12}} + 3{C_{44}}} \right) + \frac{{5\left( {{C_{11}} - {C_{12}}} \right){C_{44}}}}{{8{C_{44}} + 6\left( {{C_{11}} - {C_{12}}} \right)}}$

      $E = 9{B_{\rm{H}}}{G_{\rm{H}}}/\left( {3{B_{\rm{H}}} + {G_{\rm{H}}}} \right)$

      $v = \left( {3{B_{\rm{H}}} - 2{G_{\rm{H}}}} \right)/\left[ {2\left( {3{B_{\rm{H}}} + {G_{\rm{H}}}} \right)} \right]$

      表2中数据可以看出,ZrBeO3体弹性模量、剪切模量都比较大,理论上物质的硬度与体弹性模型密切相关,按一般经验理论初步判断它可能具有超硬性。Gao[32]基于第一性原理提出一种基于共价轨道重叠布居数的固体本征硬度值计算方法

      ${H_{\rm v}} = {\left[ {{{({H^{{\rm{Be - O}}}})}^{n_1}}{{({H^{{\rm{Zr - O}}}})}^{n_2}}} \right]^{1/(n_1 + n_2)}}$

      $H_{\rm v}^\text {μ} = A{P^\text{μ}}{(v_{\rm{b}}^\text{μ} )^{{\rm{ - }}\tfrac{5}{3}}}$

      $\upsilon _{\rm{b}}^\text{μ} = \frac{{{{\left( {{d^\text{μ} }} \right)}^3}}}{{\sum\nolimits_{\rm v} {[{{({d'})}^3}N_{\rm b}^\text{μ} ]} }} = \frac{{{{\left( {{d^\text{μ}}} \right)}^3}}}{{\sum\nolimits_{\rm v} {[{{({d^{\rm v}})}^3}{{N^\text{μ}}}/\varOmega ]} }} = \frac{{{{\left( {{d^\text{μ}}} \right)}^3}}}{{\sum\nolimits_{\rm v} {[{{({d^{\rm v}})}^3}{N^\text{μ}}]} }}$

      式中:A为常量,取740,上标Be-O/Zr-O代表其化学键;P为布居值,Pμ为μ类型化学键的Mulliken重叠布居数;${\upsilon _{\rm{b}}}$为键体积(${\nu _{\rm{b}}}$ = V/nV为原胞体积,n为原胞中所含该类键的数量);$\nu _{\rm{b}}^{\text{μ}}$为μ类型化学键的体积;上标μ、v表示原胞中不同价键成分;dμ为键长;$N_{\rm{b}}^{\text{μ}}$为单位体积内μ键的个数;Ω为晶胞体积;$N_{}^{\text{μ}}$为晶胞内μ键的总数。代入表2中数据,计算得到钙钛矿结构的ZrBeO3常压下硬度约为34.5 GPa,超过SiC的硬度[32],在超硬应用材料方面应该有较大开发潜力。材料泊松比$v$为固体材料最大拉伸强度与最大剪切强度之比,表征材料受挤压或拉伸时的收缩率或膨胀率[33],按照断裂行为的判据,低泊松比材料属于脆性材料,这一点在BH/GH得到验证:通常认为BH/GH > 1.75代表材料呈韧性[34]BH/GH < 1.75代表材料呈脆性(金刚石B/G = 0.8)。从表2中可以明显看出,若常压下ZrBeO3能稳定存在,其应为脆性材料,随着等静压力增大,材料逐渐向韧性转变。

    • 图3给出了晶胞结构优化后的ZrBeO3声子能谱,每单位晶胞有5个原子,因此有15支晶格振动格波模式。由于钙钛矿构型的体心高对称性,在G点(体心Be原子位置)15支格波呈共5组简并态,而在Q点(两Zr原子连线中点)格波谱线开裂间距最大[35],由Be和Zr原子质量差别较大所致,很明显声学支的晶格振动格波存在较多虚频(软化模式[36]),说明晶格整体在低温零压下处于热动力学不稳定状态,这可能是到目前为止实验未合成出该化合物的原因。晶胞在不同压力下的声子能谱如图4所示。可以看出,随着压力增大,声子虚频逐渐减小,等静压加到200 GPa时声子虚频消失。在273 K的温度下,对晶胞做了常压和200 GPa的动力学模拟(模拟动画见期刊官网的资源附件),结果显示200 GPa压力下晶格振动稳定,所以可认为ZrBeO3为高压生成稳定相。由于生成压力条件极高,因此实验室条件无法人工合成,但深层地壳或地核中可能存在该化合物。根据热力学稳定但动力学不活泼特性物质一般能形成力学强度高的材料[37],可以预测该化合物的力学强度很高。

      图  3  零压下的ZrBeO3声子能谱

      Figure 3.  ZrBeO3 phonon spectra at zero pressure

      图  4  不同压力下的ZrBeO3声子能谱

      Figure 4.  Phonon spectra of ZrBeO3 at different pressures

      表3对比了不同压力下ZrBeO3原子轨道及化学键布居值,可以从化学成键上说明在加压情况下ZrBeO3如何趋向于稳定。Be原子的s轨道布居值随等静压力增加而逐渐减小,p轨道布居值从1.11开始随压力增大而增加,说明钙钛矿结构ZrBeO3若要形成,Be原子不能处于基态,需要有s轨道的电子跃迁到p轨道,形成s-p杂化后与O原子成键,这可能是低压下ZrBeO3不能稳定的根本原因;同时随着压力增大,Be-O键布居值逐渐增加,键长逐渐减小,说明共价键属性增强,Be-O之间结合更紧密,从而使得八面体结构更稳定。而对于Zr原子,其sp轨道布居值随压力增大而减小,d轨道布居值增加,Zr-O键布居值和键长均减小,说明压力增加过程中,Zr-O逐渐向离子键趋势转变,电荷更多集中到O原子附近以Be-O共价键形式存在,加强稳固Be、O形成的八面体结构,Zr原子更多地起到“填隙”作用,使得钙钛矿结构稳定。

      p/
      GPa
      AtomicspdBond Be-OBond Zr-O
      PopulationBond length/nmPopulationBond length/nm
      0Be/Zr/O(3)2.34/2.26/1.821.11/6.31/4.890/1.84/00.630.173 00.620.244 7
      50Be/Zr/O(3)2.31/2.18/1.791.24/6.24/4.920/1.91/00.690.162 50.530.229 8
      100Be/Zr/O(3)2.29/2.12/1.771.31/6.21/4.930/1.97/00.720.157 70.470.223 0
      150Be/Zr/O(3)2.27/2.08/1.761.37/6.18/4.940/2.01/00.750.154 20.400.218 0
      200Be/Zr/O(3)2.26/2.04/1.751.42/6.16/4.940/2.05/00.770.151 30.330.214 0

      表 3  不同压力下原子轨道和化学键的布居值分布

      Table 3.  Atomic orbital and chemical bond population distribution at different pressures

      为更清楚地显示压力对原子间相互作用的影响,给出不同压力下的电子总态密度及Be、Zr、O分波态密度图,如图5所示。在电子总态密度图中可以很容易看出,随着压力增大,电子总态密度向高能态方向增加,且每个能带宽度均逐渐增大,尤其在高能量区域(8~28 eV)的sp电子最为明显,说明压力增大过程中,电子态密度有整体向高能级轨道跃迁的趋势,电子之间的相互作用逐渐增强,尤其是s-p电子之间。

      图  5  不同压力下ZrBeO3电子总态密度及各原子分波态密度图

      Figure 5.  Total density of state of ZrBeO3 electrons and partial density of states of each atom at different pressures

      从Be原子的电子分波态密度图中,能观察出其高能量区域中电子s-p杂化随着压力增大越来越明显。O原子分波态密度总体随压力变化不特别明显,高能区域变化趋势与总态密度相同,而在低能区域有个小波峰逐渐变大,表明该能量区域的电子密度有所增加(总体来说成分较小),可认为是Zr-O逐渐向离子键趋势转变促成。从过渡金属Zr原子分波态密度图中发现,随着压力增大,s轨道和d轨道出现分裂,这与价壳层中存在未配对电子有关。Be-O、Zr-O化学键的状态可以用“赝能隙”[38](Pseudogap,即在费米能级两侧分别有两个尖峰,而两个尖峰之间的态密度并不为零)判断,一般认为赝能隙宽度直接反映共价键强弱,越宽说明共价性越强。Be-O键的赝能隙宽度取决于O原子价带顶的p电子到Be原子导带底的p电子分布的距离,Zr-O键宽度取决于O原子价带顶的p电子到Zr原子导带底的d电子分布的距离。从图5(b)图5(d)中数值标注可读出,Be-O键的共价性随压力增大而明显增强;Zr-O键共价性质却由于Zr的d轨道随着压力增大逐渐出现未配对电子情况而变得相对复杂,其赝能隙先增大后因轨道分裂出现子峰而减小,而轨道分裂意味着出现未配对电子,未配对电子的出现应是O原子吸收了Zr原子外围的价电子所致(即Zr-O键的离子性质随压力增大而逐渐增强)。

    • 计算了优化后晶胞的电子能带结构,能带结构图说明该化合物是具有间接能隙1.918 eV的半导体材料;优化后的晶胞的光学特性计算(主要关心可见光附近波段)结果表明,该材料在紫外段300 nm、可见光段500 nm、红外段700~900 nm有较强的吸收峰,可见光段600 nm附近处有最低的反射率。计算了该材料的介电常数随入射光能量的变化频率谱线,结果表明若能合成稳定化合物,其在半导体器件、光传感、光催化领域均具有较大潜力。考虑到材料的极不稳定性,相关计算结果的数据意义不大,暂作为支撑材料附上,见附录A。

    • 在满足钙钛矿几何结构稳定性、热力学机械稳定性的基础上,基于密度泛函理论(DFT)框架,使用广义梯度近似平面波方法(GGA)构建了钙钛矿构型ZrBeO3具体的晶体模型。结合能计算表明,该构型具有热力学稳定性;弹性常数、体积模量、剪切模量、杨氏模量的计算表明,该晶体具备较高的力学强度;硬度计算表明,该结构晶体硬度可与SiC相比拟;不同压力下的泊松比、BH/GH计算表明,材料在压力增加的情况下会由脆性向韧性转变。计算了ZrBeO3的声子能谱,结果表明ZrBeO3在低温零压下热动力学不稳定。为此对比分析了不同压力下的声子能谱、不同原子轨道和化学键布居值,晶胞的电子态密度及各个原子的电子分波态密度研究表明随着压力增加,Be原子sp杂化后形成的Be-O共价键成分增强、Zr-O键离子键成分增强,在200 GPa附近时晶格动力学趋于稳定。

      图  A1  ZrBeO3的电子能带结构

      Figure A1.  Electronic band structure of ZrBeO3

      图  A2  ZrBeO3的吸收谱

      Figure A2.  Optical absorption spectrum of ZrBeO3

      图  A3  ZrBeO3的反射谱

      Figure A3.  Optical reflectivity spectrum of ZrBeO3

      图  A4  ZrBeO3介电函数谱线

      Figure A4.  Dielectric function spectrum of ZrBeO3

      附录A

参考文献 (38)

目录

    /

    返回文章
    返回