搜索

x

留言板

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

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

氟化镁高压萤石结构稳定性及热物性的数值模拟

孙小伟 宋婷 刘子江 万桂新 张磊 常文利

引用本文:
Citation:

氟化镁高压萤石结构稳定性及热物性的数值模拟

孙小伟, 宋婷, 刘子江, 万桂新, 张磊, 常文利

Numerical prediction of structural stability and thermodynamic properties for MgF2 with fluorite- type structure under high pressure

Sun Xiao-Wei, Song Ting, Liu Zi-Jiang, Wan Gui-Xin, Zhang Lei, Chang Wen-Li
PDF
HTML
导出引用
  • 氟化镁(MgF2)是工业用途广泛的重要碱土金属氟化物, 也是矿物质氟镁石的主要成分, 相比于电子结构和光学特性的研究, 人们从地球物理学的角度给予MgF2高压热物性的研究明显不够, 而组成地球下地幔矿物的高压熔化、热膨胀等热物性预测对理解地球的结构、动力学、演化及起源至关重要. 本文利用基于密度泛函理论的第一性原理方法, 通过热力学、动力学、力学稳定性计算表明萤石结构为MgF2高压结构, 根据等焓原理, 分别结合广义梯度近似和局域密度近似确定出了零温下MgF2晶体从稳定金红石结构到高压萤石结构的相转变压力为19.26 GPa和18.15 GPa, 且萤石结构至少稳定到135 GPa (相当于下地幔压力); 利用基于有效势参数模型的经典分子动力学方法, 通过模拟特定压力下MgF2体系的摩尔体积、总能随温度的变化确认了MgF2萤石结构在300—6000 K温度范围内的高温稳定性. 在此基础上, 考虑选用能够提高密堆固体平衡特性的交换相关泛函形式的广义梯度近似方法且结合准谐德拜模型, 以及利用根据从头算Hartree-Fock方法获得的数据拟合得到的可靠经验势参数结合经典分子动力学方法, 共同预测了萤石结构的MgF2在300—1500 K和0—135 GPa的温度和压力范围内的体积热膨胀系数、等温体模量、热弹性参数等重要热力学参量. 研究表明: MgF2萤石结构基于体积热膨胀系数和等温体模量得到的热弹性参数并非物态方程研究中通常假定的常数, 但在高温高压条件下, 其值接近于常数.
    MgF2 is an important member of alkaline-earth fluorides and has a wide range of applications in industry. Meanwhile, MgF2 occurs naturally as a mineral sellaite. Compared with the study of its electronic structure and optical properties, the researches of the behavior under high pressure of MgF2, especially the thermodynamic properties are still limited. The high-pressure melting, volume thermal expansion coefficient, and thermoelastic parameter of the Earth’s lower mantle mineral, like MgF2, are of interest and importance for understanding the physical nature of the functional material and for recognizing the structural compositions, dynamics, evolution and origin of the earth. Using the first-principles calculations based on density functional theory, the thermodynamic, mechanical, and dynamic stability of the fluorite-type structure for MgF2 are systematically studied. The calculations indicate that the fluorite-type structure is a high-pressure phase and it is stable at least up to 135 GPa. According to the principle of equal enthalpies, the phase transition pressure of MgF2 crystal from stable rutile structure to high pressure fluorite structure is determined to be 19.26 GPa and 18.15 GPa based on the the generalized gradient approximation and local density approximation calculations, respectively. The high-temperature structural stability of MgF2 with the fluorite-type structure is investigated and confirmed by using the classical molecular dynamics (MD) simulations by taking into account the molar volume and total energy change behavior in a temperature range from 300 to 6000 K. On the basis of previous research, the volume thermal expansion coefficient, isothermal bulk modulus, and thermoelastic parameter of MgF2 with the CaF2-type fluorite structure are predicted systematically in a temperature range from 300 to 1500 K and in a pressure range from 0 to 135 GPa with the help of the generalized gradient approximation of the revised Perdew-Burke-Ernzerhof form combined with quasiharmonic Debye model calculations and the molecular dynamics method combined with reliable interatomic potentials. An important discovery is that the thermoelastic parameter of this material under low temperature and low pressure is not a constant as assumed usually in previous studies of the equation of states, but it approaches to a constant under both high temperature and high pressure.
      通信作者: 孙小伟, sunxw_lzjtu@yeah.net
    • 基金项目: 国家级-国家自然科学基金(5156202)
      Corresponding author: Sun Xiao-Wei, sunxw_lzjtu@yeah.net
    [1]

    Appel R, Dyer C D, Lockwood J N 2002 Appl. Opt. 41 2470Google Scholar

    [2]

    Arroussi A, Ghezali M 2018 Optik 164 16

    [3]

    Wojciechowska M, Zieliński M, Pietrowski M 2003 J. Fluorine Chem. 120 1Google Scholar

    [4]

    Sun X W, Liu Z J, Song T, Quan W L, Chen Q F 2012 Phys. Scr. 85 065707Google Scholar

    [5]

    Haines J, Léger J M, Gorelli F, Klug D D, Tse J S, Li Z Q 2001 Phys. Rev. B 64 134110Google Scholar

    [6]

    Ming L C, Manghani M H 1979 Geophys. Res. Lett. 6 13Google Scholar

    [7]

    Öztürk H, Kürkçü C, Kürkçü C 2014 J. Alloys Compd. 609 185Google Scholar

    [8]

    Nelson J R, Needs R J, Pickard C J 2017 Phys. Rev. B 95 054118Google Scholar

    [9]

    Allan N L, Hines R I, Towler M D, Mackrodt W C 1994 J. Chem. Phys. 100 4710Google Scholar

    [10]

    Nishidate K, Baba M, Sato T, Nishikawa K 1995 Phys. Rev. B 52 3170

    [11]

    Catti M, Pavese A, Dovesi R, Roetti C, Causà M, 1991 Phys. Rev. B 44 3509Google Scholar

    [12]

    Nga Y A, Ong C K, 1993 J. Chem. Phys. 98 3240Google Scholar

    [13]

    Barrera G D, Taylor M B, Allan N L, Barron T H K, Kantorovich L N, Mackrodt W C 1997 J. Chem. Phys. 107 4337Google Scholar

    [14]

    Tian J H, Song T, Sun X W, Liu Z J, Quan W L, Guo P 2012 Physica B 407 551Google Scholar

    [15]

    Sun X W, Song T, Wei X P, Quan W L, Liu X B, Su W F 2014 Mater. Res. Bull. 52 151Google Scholar

    [16]

    Lin J F, Speziale S, Mao Z, Marquardt H 2013 Rev. Geophys. 51 244Google Scholar

    [17]

    Segall M D, Lindan P J, Probert M J, Pickard1C J, Hasnip P J, Clark S J, Payne M C 2002 J. Phys. Condens. Matter 14 2717Google Scholar

    [18]

    Ceperley D M, Alder B 1980 Phys. Rev. Lett. 45 566Google Scholar

    [19]

    Perdew J P, Zunger A 1981 Phys. Rev. B 23 5048Google Scholar

    [20]

    Perdew J P, Ruzsinszky A, Csonka G I, Vydrov O A, Scuseria G E, Constantin L A, Zhou X, Burke K 2008 Phys. Rev. Lett. 100 136406Google Scholar

    [21]

    Vanderbilt D 1990 Phys. Rev. B 41 7892Google Scholar

    [22]

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

    [23]

    Fischer T H, Almlof J 1992 J. Phys. Chem. 96 9768Google Scholar

    [24]

    Gonze X, Lee C 1997 Phys. Rev. B 55 10355Google Scholar

    [25]

    Karki B B, Ackland G J, Crain J 1997 J. Phys. Condens. Matter 9 8579Google Scholar

    [26]

    Fincham D 1992 Mol. Simul. 8 165Google Scholar

    [27]

    宋婷, 孙小伟, 魏小平, 欧阳玉花, 张春林, 郭鹏, 赵炜 2019 物理学报 68 126201Google Scholar

    Song T, Sun X W, Wei X P, Ouyang Y H, Zhang C L, Guo P, Zhao W 2019 Acta Phys. Sin. 68 126201Google Scholar

    [28]

    Cazorla C, Errandonea D 2013 J. Phys. Chem. C 117 11292

    [29]

    Song T, Sun X W, Liu Z J, Li J F, Tian J H 2012 Chin. Phys. B 21 037103Google Scholar

    [30]

    孙小伟, 褚衍东, 刘子江, 刘玉孝, 王成伟, 刘维民 2005 物理学报 54 5830Google Scholar

    Sun X W, Chu Y D, Liu Z J, Liu Y X, Wang C W, Liu W M 2005 Acta Phys. Sin. 54 5830Google Scholar

    [31]

    张计划, 丁建文, 卢章辉 2009 物理学报 58 1901Google Scholar

    Zhang J H, Ding J W, Lu Z H 2009 Acta Phys. Sin. 58 1901Google Scholar

    [32]

    Simanovskii D M, Schwettman H A 2003 Phys. Rev. Lett. 91 107601Google Scholar

    [33]

    Wang J, Yip S, Phillpot S R, Wolf D 1993 Phys. Rev. Lett. 71 4182Google Scholar

    [34]

    Blanco M, Francisco E, Luana V 2004 Comput. Phys. Commun. 158 57Google Scholar

    [35]

    Liu M, Lee C, Kaneko M, Nakahira K, Takano Y 2006 Appl. Opt. 45 1368Google Scholar

    [36]

    Sun X W, Liu Z J, Chen Q F, Quan W L, Chen Z G, Li Y H 2009 Mater. Res. Bull. 44 1729Google Scholar

  • 图 1  MgF2晶体(a) 金红石结构和(b) 立方萤石结构示意图, 其中大球代表Mg原子, 小球代表F原子

    Fig. 1.  Crystal structures of MgF2 with (a) the rutile-type phase and (b) the fluorite-type phase. The large and small spheres represent magnesium and fluorine atoms, respectively.

    图 2  利用(a) GGA和(b) LDA方法分别计算的MgF2晶体金红石结构和萤石结构零温下的焓随压力的关系, 内插图分别为两种结构的MgF2每个分子式的相对焓随压力的变化

    Fig. 2.  Calculated enthalpy as a function of pressure in the framework of (a) GGA and (b) LDA for MgF2 with the rutile-type and fluorite-type structures at zero temperature. In the inset, the relative enthalpy versus pressure is presented.

    图 3  利用GGA计算得到的MgF2晶体(a) 金红石结构在零压和相变压力为19.26 GPa下的能带结构以及(b) 萤石结构在相变压力为19.26和135 GPa下的能带结构

    Fig. 3.  Calculated band structures of MgF2 using GGA method: (a) The rutile-type phases at 0 and 19.26 GPa; (b) the fluorite-type phase at 19.26 and 135 GPa.

    图 4  利用GGA计算得到的MgF2晶体金红石结构在(a)零压、(b)相变压力为19.26 GPa下的声子谱和(c)萤石结构在135 GPa下的声子谱

    Fig. 4.  Calculated phonon spectra of MgF2 with the rutile-type phases at (a) 0 GPa and (b) 19.26 GPa and with (c) the fluorite-type phase at 135 GPa using GGA method.

    图 5  利用GGA近似计算的MgF2萤石结构的弹性常数Cij随外加压力的变化

    Fig. 5.  Pressure dependence of the elastic constants Cij for MgF2 with the fluorite-type structure.

    图 6  利用分子动力学模拟和第一性原理计算得到的MgF2萤石结构 (a) 在300 K下的体积比率随压力的变化和(b) 在50 GPa下的体积比率随温度的变化, 内插图为0.1 MPa下的模拟结果

    Fig. 6.  Volume ratios of MgF2 with the fluorite-type structure obtained from molecular dynamics simulations and first-principles calculations: (a) Volume ratios under different pressures at 300 K; (b) volume ratios under different temperatures at 50 GPa, where in the inset, the data at 0.1 MPa is presented.

    图 7  利用GGA和LDA计算得到的MgF2萤石结构的原胞总能随体积的变化

    Fig. 7.  Energy as a function of primitive cell volume for MgF2 with the fluorite-type structure using GGA and LDA calculations.

    图 8  利用分子动力学模拟得到不同压力下的MgF2萤石结构的(a) 摩尔体积随温度的变化和(b) 总能随温度的变化

    Fig. 8.  (a) Molar volume and (b) total energy of MgF2 with the fluorite-type structure as a function of temperature under different pressures calculated by molecular dynamics.

    图 9  分子动力学模拟得到MgF2晶体萤石结构的熔化相图及和MgO岩盐结构熔化相图的比较

    Fig. 9.  Melting phase diagram of MgF2 with the fluorite-type structure obtained from molecular dynamics (MD) simulations, in comparison with the calculated melting phase diagram of the rocksalt phase of MgO.

    图 10  模拟得到的300 K及其他不同高温(500, 1000和1500 K)下的MgF2萤石结构的体积热膨胀系数、等温体模量、热弹性参数随压力的变化

    Fig. 10.  Predicted volume thermal expansion coefficient α, isothermal bulk modulus KT, and thermoelastic parameter αKT of MgF2 with the fluorite-type structure as a function of pressure at 300 K and other different temperatures (500, 1000 and 1500 K).

    图 11  模拟得到的环境压力下及其他不同高压(50, 100和135 GPa)下的MgF2萤石结构的体积热膨胀系数、等温体模量、热弹性参数随温度的变化

    Fig. 11.  Predicted volume thermal expansion coefficient α, isothermal bulk modulus KT, and thermoelastic parameter αKT of MgF2 with the fluorite-type structure as a function of temperature at 0.1 MPa and other different pressures (50, 100 and 135 GPa).

  • [1]

    Appel R, Dyer C D, Lockwood J N 2002 Appl. Opt. 41 2470Google Scholar

    [2]

    Arroussi A, Ghezali M 2018 Optik 164 16

    [3]

    Wojciechowska M, Zieliński M, Pietrowski M 2003 J. Fluorine Chem. 120 1Google Scholar

    [4]

    Sun X W, Liu Z J, Song T, Quan W L, Chen Q F 2012 Phys. Scr. 85 065707Google Scholar

    [5]

    Haines J, Léger J M, Gorelli F, Klug D D, Tse J S, Li Z Q 2001 Phys. Rev. B 64 134110Google Scholar

    [6]

    Ming L C, Manghani M H 1979 Geophys. Res. Lett. 6 13Google Scholar

    [7]

    Öztürk H, Kürkçü C, Kürkçü C 2014 J. Alloys Compd. 609 185Google Scholar

    [8]

    Nelson J R, Needs R J, Pickard C J 2017 Phys. Rev. B 95 054118Google Scholar

    [9]

    Allan N L, Hines R I, Towler M D, Mackrodt W C 1994 J. Chem. Phys. 100 4710Google Scholar

    [10]

    Nishidate K, Baba M, Sato T, Nishikawa K 1995 Phys. Rev. B 52 3170

    [11]

    Catti M, Pavese A, Dovesi R, Roetti C, Causà M, 1991 Phys. Rev. B 44 3509Google Scholar

    [12]

    Nga Y A, Ong C K, 1993 J. Chem. Phys. 98 3240Google Scholar

    [13]

    Barrera G D, Taylor M B, Allan N L, Barron T H K, Kantorovich L N, Mackrodt W C 1997 J. Chem. Phys. 107 4337Google Scholar

    [14]

    Tian J H, Song T, Sun X W, Liu Z J, Quan W L, Guo P 2012 Physica B 407 551Google Scholar

    [15]

    Sun X W, Song T, Wei X P, Quan W L, Liu X B, Su W F 2014 Mater. Res. Bull. 52 151Google Scholar

    [16]

    Lin J F, Speziale S, Mao Z, Marquardt H 2013 Rev. Geophys. 51 244Google Scholar

    [17]

    Segall M D, Lindan P J, Probert M J, Pickard1C J, Hasnip P J, Clark S J, Payne M C 2002 J. Phys. Condens. Matter 14 2717Google Scholar

    [18]

    Ceperley D M, Alder B 1980 Phys. Rev. Lett. 45 566Google Scholar

    [19]

    Perdew J P, Zunger A 1981 Phys. Rev. B 23 5048Google Scholar

    [20]

    Perdew J P, Ruzsinszky A, Csonka G I, Vydrov O A, Scuseria G E, Constantin L A, Zhou X, Burke K 2008 Phys. Rev. Lett. 100 136406Google Scholar

    [21]

    Vanderbilt D 1990 Phys. Rev. B 41 7892Google Scholar

    [22]

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

    [23]

    Fischer T H, Almlof J 1992 J. Phys. Chem. 96 9768Google Scholar

    [24]

    Gonze X, Lee C 1997 Phys. Rev. B 55 10355Google Scholar

    [25]

    Karki B B, Ackland G J, Crain J 1997 J. Phys. Condens. Matter 9 8579Google Scholar

    [26]

    Fincham D 1992 Mol. Simul. 8 165Google Scholar

    [27]

    宋婷, 孙小伟, 魏小平, 欧阳玉花, 张春林, 郭鹏, 赵炜 2019 物理学报 68 126201Google Scholar

    Song T, Sun X W, Wei X P, Ouyang Y H, Zhang C L, Guo P, Zhao W 2019 Acta Phys. Sin. 68 126201Google Scholar

    [28]

    Cazorla C, Errandonea D 2013 J. Phys. Chem. C 117 11292

    [29]

    Song T, Sun X W, Liu Z J, Li J F, Tian J H 2012 Chin. Phys. B 21 037103Google Scholar

    [30]

    孙小伟, 褚衍东, 刘子江, 刘玉孝, 王成伟, 刘维民 2005 物理学报 54 5830Google Scholar

    Sun X W, Chu Y D, Liu Z J, Liu Y X, Wang C W, Liu W M 2005 Acta Phys. Sin. 54 5830Google Scholar

    [31]

    张计划, 丁建文, 卢章辉 2009 物理学报 58 1901Google Scholar

    Zhang J H, Ding J W, Lu Z H 2009 Acta Phys. Sin. 58 1901Google Scholar

    [32]

    Simanovskii D M, Schwettman H A 2003 Phys. Rev. Lett. 91 107601Google Scholar

    [33]

    Wang J, Yip S, Phillpot S R, Wolf D 1993 Phys. Rev. Lett. 71 4182Google Scholar

    [34]

    Blanco M, Francisco E, Luana V 2004 Comput. Phys. Commun. 158 57Google Scholar

    [35]

    Liu M, Lee C, Kaneko M, Nakahira K, Takano Y 2006 Appl. Opt. 45 1368Google Scholar

    [36]

    Sun X W, Liu Z J, Chen Q F, Quan W L, Chen Z G, Li Y H 2009 Mater. Res. Bull. 44 1729Google Scholar

  • [1] 刘泽涛, 陈博, 令伟栋, 包南云, 康冬冬, 戴佳钰. 冲击压缩下金属钯的结构相变. 物理学报, 2022, 71(3): 037102. doi: 10.7498/aps.71.20211511
    [2] 王碧涵, 李冰, 刘旭强, 王毫, 蒋升, 林传龙, 杨文革. 毫秒时间分辨同步辐射X射线衍射和高压快速加载装置及应用. 物理学报, 2022, 71(10): 100702. doi: 10.7498/aps.71.20212360
    [3] 田春玲, 刘海燕, 王彪, 刘福生, 甘云丹. 稠密流体氮高温高压相变及物态方程. 物理学报, 2022, 71(15): 158701. doi: 10.7498/aps.71.20220124
    [4] 刘泽涛, 陈博, 令伟栋, 包南云, 康冬冬, 戴佳钰. 冲击压缩下金属钯的结构相变研究. 物理学报, 2021, (): . doi: 10.7498/aps.70.20211511
    [5] 宋婷, 孙小伟, 魏小平, 欧阳玉花, 张春林, 郭鹏, 赵炜. 方镁石高压结构预测和高温结构稳定性研究. 物理学报, 2019, 68(12): 126201. doi: 10.7498/aps.68.20190204
    [6] 陶强, 马帅领, 崔田, 朱品文. 过渡金属硼化物的结构与性质. 物理学报, 2017, 66(3): 036103. doi: 10.7498/aps.66.036103
    [7] 肖宏宇, 刘利娜, 秦玉琨, 张东梅, 张永胜, 隋永明, 梁中翥. B2O3添加宝石级金刚石单晶的生长特性. 物理学报, 2016, 65(5): 050701. doi: 10.7498/aps.65.050701
    [8] 胡翠娥, 曾召益, 蔡灵仓. 极端条件下锆的动力学稳定性研究. 物理学报, 2015, 64(4): 046401. doi: 10.7498/aps.64.046401
    [9] 胡永金, 吴云沛, 刘国营, 罗时军, 何开华. ZnTe结构相变、电子结构和光学性质的研究. 物理学报, 2015, 64(22): 227802. doi: 10.7498/aps.64.227802
    [10] 濮春英, 王丽, 吕林霞, 于荣梅, 何朝政, 卢志文, 周大伟. NbSi2奇异高压相及其热力学性质的第一性原理研究. 物理学报, 2015, 64(8): 087103. doi: 10.7498/aps.64.087103
    [11] 肖宏宇, 李尚升, 秦玉琨, 梁中翥, 张永胜, 张东梅, 张义顺. 高温高压下掺硼宝石级金刚石单晶生长特性的研究. 物理学报, 2014, 63(19): 198101. doi: 10.7498/aps.63.198101
    [12] 王金荣, 朱俊, 郝彦军, 姬广富, 向钢, 邹洋春. 高压下RhB的相变、弹性性质、电子结构及硬度的第一性原理计算. 物理学报, 2014, 63(18): 186401. doi: 10.7498/aps.63.186401
    [13] 卢志文, 仲志国, 刘克涛, 宋海珍, 李根全. 高温高压下Ag-Mg-Zn合金中金属间化合物的微观结构与热动力学性质的第一性原理计算. 物理学报, 2013, 62(1): 016106. doi: 10.7498/aps.62.016106
    [14] 秦杰明, 王皓, 曾繁明, 李建利, 万玉春, 刘景和. 高温高压下MgxZn1-xO固溶体的制备. 物理学报, 2010, 59(12): 8910-8914. doi: 10.7498/aps.59.8910
    [15] 李晓兵, 赵祥永, 汪尧进, 王飞飞, 陈超, 罗豪甦. 由BaTiO3晶体结构相变时的介电特性研究其电场作用下的偶极子偏转路径. 物理学报, 2009, 58(6): 4225-4229. doi: 10.7498/aps.58.4225
    [16] 孔令刚, 康晋锋, 王 漪, 刘力锋, 刘晓彦, 张 兴, 韩汝琦. CoxTi1-xO2-δ体材中氢退火引起的铁磁性及结构相变. 物理学报, 2006, 55(3): 1453-1457. doi: 10.7498/aps.55.1453
    [17] 白锁柱, 姚 斌, 郑大方, 邢国忠, 苏文辉. 新型BCN化合物的结构表征和相转变. 物理学报, 2006, 55(11): 5740-5744. doi: 10.7498/aps.55.5740
    [18] 李成斌, 贾天卿, 孙海轶, 李晓溪, 徐世珍, 冯东海, 王晓峰, 葛晓春, 徐至展. 飞秒激光对氟化镁烧蚀机理研究. 物理学报, 2006, 55(1): 217-220. doi: 10.7498/aps.55.217
    [19] 孙小伟, 褚衍东, 刘子江, 刘玉孝, 王成伟, 刘维民. 高温高压下闪锌矿相GaN结构和热力学特性的分子动力学研究. 物理学报, 2005, 54(12): 5830-5836. doi: 10.7498/aps.54.5830
    [20] 胡林华, 戴松元, 王孔嘉. 溶胶-凝胶法制备的纳米TiO2结构相变及晶体生长动力学. 物理学报, 2003, 52(9): 2135-2139. doi: 10.7498/aps.52.2135
计量
  • 文章访问数:  9366
  • PDF下载量:  143
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-02-26
  • 修回日期:  2020-04-29
  • 上网日期:  2020-05-13
  • 刊出日期:  2020-08-05

/

返回文章
返回