搜索

x
中国物理学会期刊

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

CSTR: 32037.14.aps.69.20200289

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

CSTR: 32037.14.aps.69.20200289
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.

     

    目录

    /

    返回文章
    返回