铁的自旋转变对菱镁矿热力学性质的影响

马超杰 吴潇 马阳阳 何开华 姬广富

引用本文:
Citation:

铁的自旋转变对菱镁矿热力学性质的影响

    作者简介: 马超杰(1994-),男,硕士研究生,主要从事高温高压材料的模拟计算研究.E-mail: machaojie1994@126.com;
    通讯作者: 何开华, khhe@cug.edu.cn
  • 中图分类号: O521.2

Effect of Spin Transition of Iron on Thermodynamic Properties of Magnesiosiderite

    Corresponding author: HE Kaihua, khhe@cug.edu.cn ;
  • CLC number: O521.2

  • 摘要: 含铁菱镁矿(Mg,Fe)CO3是碳进入地球深部的主要载体之一,铁的进入会引起矿物物理性质的变化。采用第一性原理计算方法,研究了菱镁矿含铁及铁的自旋转变对菱镁矿热力学性质的影响。含铁菱镁矿的低自旋态体积比不含铁菱镁矿小;高自旋态在低温端的体积比不含铁菱镁矿略微增大,在高温端却减小;在所研究的温压范围内,低自旋态的体积始终比高自旋态的体积小。含铁菱镁矿高自旋态的热膨胀系数减小,而自旋转变会导致热膨胀系数增加。考虑高低两种自旋态共存时的热力学性质时,计算结果表明:自旋态共存时的热膨胀系数、速度在自旋共存区间内分别呈现异常增大峰和异常减小峰,并且这些异常变化峰随着温度的升高向高压方向移动。
  • 图 1  (Mg0.5Fe0.5)CO3的LS态与HS态的焓差ΔHHHSHLS)随压力的变化

    Figure 1.  Pressure dependence of the enthalpy difference (ΔH) between LS state and HS state (HHSHLS) for (Mg0.5Fe0.5)CO3

    图 2  (Mg0.5Fe0.5)CO3的HS态与LS态的体积V(a)和热膨胀系数α(b)随温度的变化关系

    Figure 2.  Temperature dependence of the (a) volume V, (b) thermal expansion coefficient α for both the HS and LS state of (Mg0.5Fe0.5)CO3

    图 3  MgCO3与(Mg0.5Fe0.5)CO3的定容比热容${C_{V,\rm{m}}}$(a)、格临爱森常数γ(b)、体变模量KT(c)与温度的关系

    Figure 3.  Temperature dependence of the (a) specific heat capacity of constant volume ${C_{V,\rm{m}}}$, (b) Grüneisen parameter γ, (c) bulk modulus KT of MgCO3 and (Mg0.5Fe0.5)CO3

    图 4  n随温度和压力的变化关系

    Figure 4.  Ratio n as the function of pressure at different temperature

    图 5  ${\left. {\dfrac{{\partial n}}{{\partial T}}} \right|_p}$随温度的变化关系(a) 及${\left. {\dfrac{{\partial n}}{{\partial p}}} \right|_T}$随压力的变化关系(b)

    Figure 5.  ${\left. {\dfrac{{\partial n}}{{\partial T}}} \right|_p}$as the function of temperature (a) and ${\left. {\dfrac{{\partial n}}{{\partial p}}} \right|_T}$as the function of pressure (b)

    图 6  MS态的体积V(a)、热膨胀系数α(b)、速度v(c)随压力的变化关系

    Figure 6.  Volume V (a), thermal expansion coefficient α (b), velocity v (c) of MS state as the function of pressure

  • [1] MANNING C E, SHOCK E L, SVERJENSKY D A. The chemistry of carbon in aqueous fluids at crustal and upper-mantle conditions: experimental and theoretical constraints [J]. Reviews in Mineralogy and Geochemistry, 2013, 75(1): 109–148. doi: 10.2138/rmg.2013.75.5
    [2] DASGUPTA R, HIRSCHMANN M M. Melting in the Earth’s deep upper mantle caused by carbon dioxide [J]. Nature, 2006, 440(7084): 659–662. doi: 10.1038/nature04612
    [3] ROHRBACH A, SCHMIDT M W. Redox freezing and melting in the Earth’s deep mantle resulting from carbon-iron redox coupling [J]. Nature, 2011, 472(7342): 209–212. doi: 10.1038/nature09899
    [4] WALTER M J, KOHN S C, ARAUJO D, et al. Deep mantle cycling of oceanic crust: evidence from diamonds and their mineral inclusions [J]. Science, 2011, 334(6052): 54–57. doi: 10.1126/science.1209300
    [5] LAVINA B, DERA P, DOWNS R T, et al. Siderite at lower mantle conditions and the effects of the pressure-induced spin-pairing transition [J]. Geophysical Research Letters, 2009, 36(23): 23306. doi: 10.1029/2009GL039652
    [6] LAVINA B, DERA P, DOWNS R T, et al. Effect of dilution on the spin pairing transition in rhombohedral carbonates [J]. High Pressure Research, 2010, 30(2): 224–229. doi: 10.1080/08957959.2010.485391
    [7] LAVINA B, DERA P, DOWNS R T, et al. Structure of siderite FeCO3 to 56 GPa and hysteresis of its spin-pairing transition [J]. Physical Review B, 2010, 82(6): 064110. doi: 10.1103/PhysRevB.82.064110
    [8] NAGAI T, ISHIDO T, SETO Y, et al. Pressure-induced spin transition in FeCO3-siderite studied by X-ray diffraction measurements [J]. Journal of Physics: Conference Series, 2010, 215(1): 012002.
    [9] FARFAN G, WANG S B, MA H W, et al. Bonding and structural changes in siderite at high pressure [J]. American Mineralogist, 2012, 97(8/9): 1421–1426.
    [10] LIN J F, LIU J, JACOBS C, et al. Vibrational and elastic properties of ferromagnesite across the electronic spin-pairing transition of iron [J]. American Mineralogist, 2012, 97(4): 583–591. doi: 10.2138/am.2012.3961
    [11] LIU J, LIN J F, MAO Z, et al. Thermal equation of state and spin transition of magnesiosiderite at high pressure and temperature [J]. American Mineralogist, 2014, 99(1): 84–93. doi: 10.2138/am.2014.4553
    [12] LIU J, LIN J F, PRAKAPENKA V B. High-pressure orthorhombic ferromagnesite as a potential deep-mantle carbon carrier [J]. Scientific Reports, 2015, 5(1): 07640. doi: 10.1038/srep07640
    [13] HAZEN R M, JONES A P, BAROSS J A. Carbon in Earth [J]. Reviews in Mineralogy and Geochemistry, 2013, 75(1).
    [14] ISSHIKI M, IRIFUNE T, HIROSE K, et al. Stability of magnesite and its high-pressure form in the lowermost mantle [J]. Nature, 2004, 427(6969): 60–63. doi: 10.1038/nature02181
    [15] OGANOV A R, ONO S, MA Y M, et al. Novel high-pressure structures of MgCO3, CaCO3 and CO2 and their role in Earth’s lower mantle [J]. Earth and Planetary Science Letters, 2008, 273(1/2): 38–47.
    [16] LIN J F, SPEZIALE S, MAO Z, et al. Effects of the electronic spin transitions of iron in lower mantle minerals: implications for deep mantle geophysics and geochemistry [J]. Reviews of Geophysics, 2013, 51(2): 244–275. doi: 10.1002/rog.20010
    [17] SPEZIALE S, MILNER A, LEE V E, et al. Iron spin transition in Earth’s mantle [J]. Proceedings of the National Academy of Sciences of the United States of America, 2005, 102(50): 17918–17922. doi: 10.1073/pnas.0508919102
    [18] TSUCHIYA T, WENTZCOVITCH R M, DA SILVAC R S, et al. Spin transition in magnesiowüstite in Earth’s lower mantle [J]. Physical Review Letters, 2006, 96(19): 198501. doi: 10.1103/PhysRevLett.96.198501
    [19] GONCHAROV A F, STRUZHKIN V V, JACOBSEN S D. Reduced radiative conductivity of low-spin (Mg, Fe)O in the lower mantle [J]. Science, 2006, 312(5777): 1205–1208. doi: 10.1126/science.1125622
    [20] LIN J F, VANKÓ G, JACOBSEN S D, et al. Spin transition zone in Earth’s lower mantle [J]. Science, 2007, 317(5845): 1740–1743. doi: 10.1126/science.1144997
    [21] CROWHURST J C, BROWN J M, GONCHAROV A F, et al. Elasticity of (Mg, Fe)O through the spin transition of iron in the lower mantle [J]. Science, 2008, 319(5862): 451–453. doi: 10.1126/science.1149606
    [22] WENTZCOVITCH R M, JUSSTO J F, WU Z, et al. Anomalous compressibility of ferropericlase throughout the iron spin cross-over [J]. Proceedings of the National Academy of Sciences of the United States of America, 2009, 106(21): 8447–8452. doi: 10.1073/pnas.0812150106
    [23] MARQUARDT H, SPEZIALE S, REICHMANN H J, et al. Elastic shear anisotropy of ferropericlase in Earth’s lower mantle [J]. Science, 2009, 324(5924): 224–226. doi: 10.1126/science.1169365
    [24] HSU H, UMEMOTO K, WU Z Q, et al. Spin-state crossover of iron in lower-mantle minerals: results of DFT+U investigations [J]. Reviews in Mineralogy and Geochemistry, 2010, 71(1): 169–199. doi: 10.2138/rmg.2010.71.09
    [25] ANTONANGELI D, SIEBERT J, ARACNE C M, et al. Spin crossover in ferropericlase at high pressure: aseismologically transparent transition? [J]. Science, 2011, 331(6013): 64–67. doi: 10.1126/science.1198429
    [26] WU Z Q, JUSTO J F, WENTZCOVITCH R M. Elastic anomalies in a spin-crossover system: ferropericlase at lower mantle conditions [J]. Physical Review Letters, 2013, 110(22): 228501. doi: 10.1103/PhysRevLett.110.228501
    [27] HOLMSTRÖM E, STIXRUDE L. Spin crossover in ferropericlase from first-principles molecular dynamics [J]. Physical Review Letters, 2015, 114(11): 117202. doi: 10.1103/PhysRevLett.114.117202
    [28] SONG Y L, HE K H, SUN J, et al. Effects of iron spin transition on the electronic structure, thermal expansivity and lattice thermal conductivity of ferropericlase: a first principles study [J]. Scientific Reports, 2019, 9(1): 4172. doi: 10.1038/s41598-019-40454-4
    [29] CERANTOLA V, MCCAMMON C, KUPENKO I, et al. High-pressure spectroscopic study of siderite (FeCO3) with a focus on spin crossover [J]. American Mineralogist, 2015, 100(11/12): 2670–2681.
    [30] MATTILA A, PYLKKÄNEN T, RUEFF J-P, et al. Pressure induced magnetic transition in siderite FeCO3 studied by X-ray emission spectroscopy [J]. Journal of Physics: Condensed Matter, 2007, 19(38): 386206. doi: 10.1088/0953-8984/19/38/386206
    [31] SHI H, LUO W L, JOHANSSON B, et al. First-principles calculations of the electronic structure and pressure-induced magnetic transition in siderite FeCO3 [J]. Physical Review B, 2008, 78(15): 155119. doi: 10.1103/PhysRevB.78.155119
    [32] HSU H, HUANG S C. Spin crossover and hyperfine interactions of iron in (Mg, Fe)CO3 ferromagnesite [J]. Physical Review B, 2016, 94(6): 060404. doi: 10.1103/PhysRevB.94.060404
    [33] FU S, YANG J, LIN J F. Abnormal elasticity of single-crystal magnesiosiderite across the spin transition in Earth’s lower mantle [J]. Physical Review Letters, 2017, 118(3): 036402. doi: 10.1103/PhysRevLett.118.036402
    [34] KRESSE G, FUTHMÜLLER J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set [J]. Physical Review B, 1996, 54(16): 11169–11186. doi: 10.1103/PhysRevB.54.11169
    [35] KRUKAU A V, VYDROV O A, IZMAYLOV A F, et al. Influence of the exchange screening parameter on the performance of screened hybrid functional [J]. The Journal of Chemical Physics, 2006, 125(22): 224106. doi: 10.1063/1.2404663
    [36] MONKHORST H J, PACK J D. Special points for Brillouin-zone integrations [J]. Physical Review B, 1976, 13(12): 5188–5192. doi: 10.1103/PhysRevB.13.5188
    [37] TOGO A, OBA F, TANAKA I. First-principles calculations of the ferroelastic transition between rutile-type and CaCl2-type SiO2 at high pressures [J]. Physical Review B, 2008, 78(13): 4106.
    [38] WU Z, JUSTO J F, DA SILVAC R S, et al. Anomalous thermodynamic properties in ferropericlase throughout its spin crossover [J]. Physical Review B, 2009, 80(1): 014409. doi: 10.1103/PhysRevB.80.014409
    [39] KOMABAYASHI T, HIROSE K, NAGAYA Y, et al. High-temperature compression of ferropericlase and the effect of temperature on iron spin transition [J]. Earth and Planetary Science Letters, 2010, 297(3/4): 691–699.
    [40] MAO Z, LIN J F, LIU J, et al. Thermal equation of state of lower-mantle ferropericlase across the spin crossover [J]. Geophysical Research Letters, 2011, 38(23): 23308.
  • [1] 张乐婷赵宇宏孙远洋邓世杰吉如意韩培德 . 压力作用下Mg2X(X=Si, Ge)相热力学性质的第一性原理研究. 高压物理学报, 2018, 32(3): 032201-1-032201-8. doi: 10.11858/gywlxb.20170630
    [2] 常少梅 . 高压下立方BC3的力学和热力学性质. 高压物理学报, 2018, 32(2): 021101-1-021101-8. doi: 10.11858/gywlxb.20170640
    [3] 郝爱民周铁军朱岩刘鑫 . 高压下TiN的结构转变、弹性和热力学性质的第一性原理计算. 高压物理学报, 2012, 26(4): 395-401. doi: 10.11858/gywlxb.2012.04.006
    [4] 李胜旨刘锦超杨向东郭艳锋许海全 . Mn、Fe掺杂ZnS的第一性原理计算. 高压物理学报, 2010, 24(6): 449-454 . doi: 10.11858/gywlxb.2010.06.008
    [5] 温新竹彭玉颜刘明真 . 钙钛矿结构ZrBeO3稳定性的第一性原理研究. 高压物理学报, 2020, 34(1): 011202-1-011202-9. doi: 10.11858/gywlxb.20190802
    [6] 赵凯蔡灵仓张修路罗雰 . 过渡金属钼高压物性的第一原理研究. 高压物理学报, 2016, 30(2): 94-100. doi: 10.11858/gywlxb.2016.02.002
    [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] 谭心贾亦超刘学杰 . 外压下ZrN结构相变的第一性原理研究. 高压物理学报, 2014, 28(2): 168-174. doi: 10.11858/gywlxb.2014.02.006
    [10] 韩林马麦宁徐志双周晓亚 . 辉石结构与相变的第一性原理研究. 高压物理学报, 2017, 31(2): 125-134. doi: 10.11858/gywlxb.2017.02.004
    [11] 郝军华吴志强王铮金庆华李宝会丁大同 . 高压下SiO2的第一性原理计算. 高压物理学报, 2010, 24(4): 260-266 . doi: 10.11858/gywlxb.2010.04.004
    [12] 孙磊罗坤刘兵韩俏怡王小雨梁子太赵智胜 . 金属性硅同素异形体的第一性原理研究. 高压物理学报, 2019, 33(2): 020103-1-020103-8. doi: 10.11858/gywlxb.20190705
    [13] 杨龙星刘雷刘红易丽顾小雨 . 高压下石榴子石结构和弹性的第一性原理研究. 高压物理学报, 2019, 33(6): 060104-1-060104-11. doi: 10.11858/gywlxb.20190785
    [14] 丁迎春刘海军蒋孟衡陈敏陈勇明 . 高压下BeP2N4结构相变和电子结构的第一性原理计算. 高压物理学报, 2012, 26(6): 674-680. doi: 10.11858/gywlxb.2012.06.012
    [15] 刘红魏冬青赵纪军郭永新龚自正 . 高压下液态硝基甲烷的分子动力学模拟. 高压物理学报, 2004, 18(4): 319-327 . doi: 10.11858/gywlxb.2004.04.006
    [16] 任晓光崔雪菡吴宝嘉顾广瑞 . 高压下金属间化合物CaAlSi的电子结构和晶格动力学性质. 高压物理学报, 2014, 28(2): 161-167. doi: 10.11858/gywlxb.2014.02.005
    [17] 郝爱民杨晓翠赵玉伟刘鑫宋爱君张卫国辛伟 . 高压下钡的硫化物的结构相变和光学性质研究. 高压物理学报, 2010, 24(2): 113-119 . doi: 10.11858/gywlxb.2010.02.006
    [18] 李鑫马雪姣高文泉刘艳辉 . 高压下Ir2P晶体结构预测与物理性质. 高压物理学报, 2019, 33(1): 011103-1-011103-6. doi: 10.11858/gywlxb.20180645
    [19] 何林尹君 . 空位点缺陷对金刚石在高压下的电子结构和光学性质的影响. 高压物理学报, 2013, 27(6): 802-806. doi: 10.11858/gywlxb.2013.06.002
    [20] 朱姗姗殷春浩徐振坤吴彩萍侯磊田 . 表面吸附氧原子的GaAs电子结构及光学性质的研究. 高压物理学报, 2013, 27(3): 468-472. doi: 10.11858/gywlxb.2013.03.023
  • 加载中
图(6)
计量
  • 文章访问数:  2039
  • 阅读全文浏览量:  1089
  • PDF下载量:  21
出版历程
  • 收稿日期:  2019-12-06
  • 录用日期:  2019-12-26
  • 网络出版日期:  2020-04-03
  • 刊出日期:  2020-04-01

铁的自旋转变对菱镁矿热力学性质的影响

    作者简介:马超杰(1994-),男,硕士研究生,主要从事高温高压材料的模拟计算研究.E-mail: machaojie1994@126.com
    通讯作者: 何开华, khhe@cug.edu.cn
  • 1. 中国地质大学(武汉)数学与物理学院,湖北 武汉 430074
  • 2. 中国工程物理研究院流体物理研究所冲击波物理与爆轰物理重点实验室,四川 绵阳 621999

摘要: 含铁菱镁矿(Mg,Fe)CO3是碳进入地球深部的主要载体之一,铁的进入会引起矿物物理性质的变化。采用第一性原理计算方法,研究了菱镁矿含铁及铁的自旋转变对菱镁矿热力学性质的影响。含铁菱镁矿的低自旋态体积比不含铁菱镁矿小;高自旋态在低温端的体积比不含铁菱镁矿略微增大,在高温端却减小;在所研究的温压范围内,低自旋态的体积始终比高自旋态的体积小。含铁菱镁矿高自旋态的热膨胀系数减小,而自旋转变会导致热膨胀系数增加。考虑高低两种自旋态共存时的热力学性质时,计算结果表明:自旋态共存时的热膨胀系数、速度在自旋共存区间内分别呈现异常增大峰和异常减小峰,并且这些异常变化峰随着温度的升高向高压方向移动。

English Abstract

  • 含碳固溶体的存在会影响地球内部的物理和化学性质[1-4]。近年来,菱镁矿(MgCO3)被认为是碳进入地球深部的主要载体之一,因其在地球深部碳循环中的关键作用而引起广泛关注[5-13]。Isshiki等[14]、Oganov等[15]的研究表明,菱镁矿在下地幔的温压条件下能够稳定存在;Hazen等[13]、Oganov等[15]开展的高温高压实验研究揭示,菱镁矿和菱铁矿的固溶体[(Mg,Fe)CO3]可以在低于100 GPa的压力条件下稳定存在。

    铁是多价态的过渡金属,会对菱镁矿和菱铁矿固溶体的性质产生非常重要的影响。此外,铁的自旋在一定的压力和温度条件下可以发生转变,引起物理性质变化[16]。国内外学者对铁方镁石(Mg,Fe)O的自旋转变开展了大量研究,结果表明,铁方镁石中的铁(Fe2+)在40~50 GPa范围内会从高自旋(High spin, HS)态向低自旋(Low spin, LS)态转变,并伴随着结构、电子、光学、弹性和热力学性质的异常变化[16-28]。对于含铁的菱镁矿,已有的高压穆斯堡尔光谱[29]、X射线发射光谱[30]、激光拉曼光谱[10]、X射线衍射[5, 8, 11]等实验研究和第一性原理计算研究[31-32]都证实,(Mg,Fe)CO3中的铁在约45 GPa时从HS态向LS态转变,导致(Mg,Fe)CO3的体积比MgCO3减小6%~10%[11]。Liu等[11]、Hsu等[32]、Fu等[33]利用理论计算和实验相结合的方法,详细讨论了弹性和地震波速在自旋转变时的变化,得到了非常有意义的结果。其他热力学参数如热膨胀系数、格临爱森常数和比热容等在自旋转变条件下的性质尚缺少报道。

    本研究利用第一性原理计算方法,开展高温高压下(Mg,Fe)CO3在HS、LS及HS和LS态共存时的混合自旋(Mixed spin, MS)态的热力学性质(热膨胀系数、体变模量、体积、速度等热力学参数)研究,并与不含铁的MgCO3的相关性质进行对比,分析引起(Mg,Fe)CO3热力学性质变化的机制。研究结果可为研究地幔深部碳的行为以及地幔在全球碳循环中的作用提供制约因素。

    • 菱镁矿MgCO3属于三角晶系,其空间群为$R\bar 3c$,原胞中包含10个原子(2个Mg原子、2个C原子和6个O原子)。为了讨论铁及其在高温高压下的自旋转变对菱镁矿物理性质的影响,本研究选择用1个Fe原子替换1个Mg原子,得到Fe和Mg的摩尔比为1∶1,分子式为(Mg0.5Fe0.5)CO3

      几何结构优化和相关能量的计算[34]采用基于密度泛函理论(Density of functional theory, DFT)的第一性原理分子动力学计算软件VASP完成。由于结构中有铁存在,需要考虑强关联作用,因此计算中考虑了Hubbard参数U的影响,即LDA+U。已有的研究表明,U值会随自旋转变而改变,本研究选取Tsuchiya等[18]、Krukau等[35]通过线性响应理论计算得到的U值,分别为ULS = 5.3 eV,UHS = 4.0 eV。K点采用Monkhorst-Pack方法生成以Γ点为中心的15 × 15 × 15网格,截断能设置为1 000 eV[36]。总能量的收敛阈值设置为1 × 10−6 eV/cell,原子力的收敛阈值为1 × 10−3 eV/cell。采用基于准谐近似的PHONOPY软件计算得到二阶原子间力常数(IFCs)、声子谱和其他热力学参数[37]。构造求解高对称q点声子频率本征值的动力学矩阵,然后利用倒空间中动力学矩阵的傅里叶变换,计算其他一般q点的声子频率。所有计算中,以Γ为中心的q点网格选取为30 × 30 × 30。计算热力学参数之前,用VASP软件完成不同位移或不同体积下结构的能量计算,计算过程中的参数设置与上述几何优化参数设置相同。

    • 根据已有的研究结果,MS态时的吉布斯自由能可表示为

      $ G{\rm{(}}n{\rm{,}}p{\rm{,}}T{\rm{)}} = n{G_{{\rm{LS}}}}{\rm{(}}p{\rm{,}}T{\rm{)}} + {\rm{(}}1 - n{\rm{)}}{G_{{\rm{HS}}}}{\rm{(}}p{\rm{,}}T{\rm{)}} + {G_{\rm{mix}}} $

      式中:n为LS态在MS态中所占的百分比,p为压力,T为温度,GHSGLS分别为HS态和LS态的吉布斯自由能,Gmix为HS-LS混合态的吉布斯自由能。Gmix可以表示为

      $ {G_{{\rm{mix}}}}{\rm{(}}n{\rm{)}} = {k_{\rm{B}}}T{X_{{\rm{Fe}}}}{\rm{[}}n{\rm{ln}}\;n + {\rm{(1 - }}n{\rm{)ln(1 - }}n{\rm{)]}} $

      式中:XFe为铁在含铁菱镁矿中的摩尔分数,kB为玻尔兹曼常数。在给定的压力和温度条件下,通过求最小化MS态吉布斯自由能得到LS态的占比n

      $n = {[1 + m(2s + 1)\exp (\Delta {G_{{\rm{LS}}-{\rm{HS}}}}/{X_{{\rm{Fe}}}}{k_{\rm{B}}}T)]^{ - 1}}$

      式中:ΔGLS−HS为(Mg0.5Fe0.5)CO3在LS态和HS态的吉布斯自由能之差。对于Fe2+来说,自旋和轨道简并量子数分别为s = 2,m = 3。得到n后可以推导出MS态的热膨胀系数α(n)和等温体积模量KT(n)分别为

      $\alpha (n) = \frac{{nV_{\rm{LS}}\alpha _{\rm{LS}} + (1 - n)V_{\rm{HS}}\alpha _{\rm{HS}} + (V_{\rm{LS}} - V_{\rm{HS}}){{\left. {\dfrac{{\partial n}}{{\partial T}}} \right|}_p}}}{{V(n)}}$

      ${K_T}(n) = \frac{{V(n)}}{{\dfrac{{nV_{\rm{LS}}}}{{K_{T,{\rm LS}}}} + (1 - n)\dfrac{{V_{\rm{HS}}}}{{K_{T,{\rm HS}}}} - (V_{\rm{LS}} - V_{\rm{HS}}){{\left. {\dfrac{{\partial n}}{{\partial p}}} \right|}_T}}}$

      式中:${\left. {\dfrac{{\partial n}}{{\partial p}}} \right|_T}$${\left. {\dfrac{{\partial n}}{{\partial T}}} \right|_p}$为LS态占比的导数。根据式(4)和式(5)导出其他参数的表达式

      ${K_S} = {K_T}(1{\rm{ + }}{\alpha ^2}{K_T}V/{C_{V,\rm{m}}})$

      $\gamma = {{\alpha {K_T}V} / {{C_{V,\rm{m}}}}}$

      $v = \sqrt {{{{K_S}} / \rho }} $

      式中:KS为绝热体积模量,${C_{V,\rm{m}}}$为定容比热容,ρ为密度。根据上述公式可以计算出MS态下的体积V、速度v和格临爱森常数γ

    • 含铁菱镁矿(Mg0.5Fe0.5)CO3的LS态与HS态的焓差ΔHHHSHLS)随压力的变化趋势如图1所示。由图1可知:在低压时,ΔH < 0,HS态的焓值小,因此HS态更稳定;在高压时,ΔH > 0,LS态的焓值变小,LS态变得更稳定。本研究计算得到的自旋转变压力为44.5 GPa,与实验观察得到的结果(40~52 GPa)很好地符合[5, 10-11, 29-30],与Hsu等[32]的计算结果(48 GPa)也符合较好,较小的差异源于铁在含铁菱镁矿中的摩尔分数不同。

      图  1  (Mg0.5Fe0.5)CO3的LS态与HS态的焓差ΔHHHSHLS)随压力的变化

      Figure 1.  Pressure dependence of the enthalpy difference (ΔH) between LS state and HS state (HHSHLS) for (Mg0.5Fe0.5)CO3

    • 图2(a)给出了菱镁矿含铁前后的体积对比。由图2(a)可以看出,菱镁矿含铁后HS态对应的体积比LS态对应的体积大,在30 GPa、300 K的温压条件下,(Mg0.5Fe0.5)CO3在HS和LS态下的体积分别为77.138 10 Å3和73.366 07 Å3,自旋转变后体积减小5.0%左右。对比含铁菱镁矿与不含铁菱镁矿的体积可知,含铁菱镁矿LS态的体积减小,HS态的体积则在低温端增大,在高温端减小,说明含铁菱镁矿的HS态对应的热膨胀系数小于不含铁菱镁矿。本研究中自旋转变导致的体积变化幅度比Liu等[11]实验测量得到的变化幅度稍小,这是由于实验所用样品为(Mg0.35Fe0.65)CO3,体积变化幅度的差异源于测量样品中铁的摩尔分数不同。

      图  2  (Mg0.5Fe0.5)CO3的HS态与LS态的体积V(a)和热膨胀系数α(b)随温度的变化关系

      Figure 2.  Temperature dependence of the (a) volume V, (b) thermal expansion coefficient α for both the HS and LS state of (Mg0.5Fe0.5)CO3

      热膨胀系数是研究矿物热力学性质的重要参数,本研究计算了含铁对菱镁矿热膨胀系数的影响。图2(b)分别给出了含铁菱镁矿的HS、LS态在30 GPa和60 GPa时的热膨胀系数,为了便于比较,图2(b)中给出了对应压力下MgCO3的热膨胀系数以及部分前人的结果[11, 32]。由图2(b)可知,在相同的温压条件下,(Mg0.5Fe0.5)CO3的HS与LS两种自旋态的热膨胀系数均小于MgCO3,即含铁会导致菱镁矿的热膨胀系数急剧减小。在30 GPa、300 K的温压条件下,计算得到MgCO3与(Mg0.5Fe0.5)CO3的热膨胀系数分别为8.05 × 10−6 K−1和3.13 × 10−6 K−1,减小幅度为61.0%;(Mg0.5Fe0.5)CO3的HS和LS态的热膨胀系数分别为3.13 × 10−6 K−1和4.22 × 10−6 K−1,铁的自旋转变导致的热膨胀系数增加幅度为34.8%。此外,热膨胀系数会随着压力的增加而减小,并且对温度的依赖性减弱。这与含铁方镁石在低压下HS态的热膨胀系数比LS态大,而在高压下HS态的热膨胀系数反而比LS态小的变化趋势有所差别[28]。根据格临爱森定律

      $\alpha = {{\gamma {C_{V,\rm{m}}}} / {{K_T}V}}$

      由于影响热膨胀系数的主要物理量为格临爱森常数γ、定容比热容CV,m、体积模量ΚT及体积V,因此可以从上述热力学参数的变化来解释含铁菱镁矿热膨胀系数的变化机制。

      图3给出了(Mg0.5Fe0.5)CO3的热力学参数随温度和压力的变化关系。由图3可知,定容比热容在低温段先随温度升高而急剧增大,随后趋于饱和,MgCO3的定容比热容为231.75 J/(mol·K),(Mg0.5Fe0.5)CO3在HS、LS态时的定容比热容分别为232.99 J/(mol·K)和231.58 J/(mol·K),彼此之间的差异很小,因此定容比热容对热膨胀系数的影响可以忽略不计。由1.2节的研究结果可知,菱镁矿含铁及铁的自旋转变引起的体积变化在5.0%左右,无法成为热膨胀系数剧烈变化的决定性因素。由图3(b)图3(c)可知,菱镁矿含铁后格临爱森常数明显减小,由HS态转变为LS态后明显增大,与热膨胀系数的变化趋势很好地符合。需要注意的是,含铁菱镁矿LS态的格临爱森常数γ比不含铁时大,而热膨胀系数却比不含铁时小,由此可以推断,体变模量KT对热膨胀系数的变化也起到至关重要的作用。因此,菱镁矿含铁及铁的自旋转变导致的热膨胀系数变化是由格临爱森常数γ与体变模量KT共同决定的。

      图  3  MgCO3与(Mg0.5Fe0.5)CO3的定容比热容${C_{V,\rm{m}}}$(a)、格临爱森常数γ(b)、体变模量KT(c)与温度的关系

      Figure 3.  Temperature dependence of the (a) specific heat capacity of constant volume ${C_{V,\rm{m}}}$, (b) Grüneisen parameter γ, (c) bulk modulus KT of MgCO3 and (Mg0.5Fe0.5)CO3

    • Liu等[11]、Hsu等[32]、Fu等[33] 的研究表明,含铁菱镁矿的MS态在高温高压下会引起弹性及声速的异常变化。研究MS态热力学参数的前提是计算出自旋共存时LS态的占比n。根据式(1)~式(8)计算出MS态下LS的占比nnpT的一阶导数${\left. {\dfrac{{\partial n}}{{\partial T}}} \right|_p}$${\left. {\dfrac{{\partial n}}{{\partial p}}} \right|_T}$,以及αVv等热力学参数。

      图4给出了不同温度下n随压力的变化关系。由图4可知:在低温时,n急剧增大到1,即HS态与LS态共存的压力区间很窄;随着温度升高,HS态与LS态共存的压力区间明显增大;当温度为300 K时,HS态与LS态共存的压力区间约为5 GPa;温度为1 500 K时,两者共存的压力区间增大到约15 GPa,与Liu等[11]、Hsu等[32]通过实验和理论计算得到的趋势一致。同时,通过计算发现:自旋转变的压力随温度的升高而增加;温度为300 K时,自旋转变发生在约40 GPa;而温度为3 000 K时,自旋转变发生在约60 GPa。该自旋转变行为在铁方镁石的理论计算和实验研究中也有报道[18, 20, 22, 38-40]。HS态与LS态共存区间和自旋转变随压力的变化关系决定了n对温度和压力的一阶导数,如图5所示。计算表明:当压力为50 GPa时,${\left. {\dfrac{{\partial n}}{{\partial T}}} \right|_p}$异常变化的起始温度为500 K左右;70 GPa时,${\left. {\dfrac{{\partial n}}{{\partial T}}} \right|_p}$异常变化的起始温度为1 500 K左右。同时${\left. {\dfrac{{\partial n}}{{\partial T}}} \right|_p}$峰的宽度也由50 GPa时的1 000 K增大到70 GPa时的1 700 K,由此可知,自旋转变引起的${\left. {\dfrac{{\partial n}}{{\partial T}}} \right|_p}$峰随着压强的增加向高温方向移动,同时峰的宽度也随之增加。图5(b)给出了${\left. {\dfrac{{\partial n}}{{\partial p}}} \right|_T}$随压强的变化关系,其峰高和峰宽随压力的变化趋势与${\left. {\dfrac{{\partial n}}{{\partial T}}} \right|_p}$类似。

      图  4  n随温度和压力的变化关系

      Figure 4.  Ratio n as the function of pressure at different temperature

      图  5  ${\left. {\dfrac{{\partial n}}{{\partial T}}} \right|_p}$随温度的变化关系(a) 及${\left. {\dfrac{{\partial n}}{{\partial p}}} \right|_T}$随压力的变化关系(b)

      Figure 5.  ${\left. {\dfrac{{\partial n}}{{\partial T}}} \right|_p}$as the function of temperature (a) and ${\left. {\dfrac{{\partial n}}{{\partial p}}} \right|_T}$as the function of pressure (b)

      因为${\left. {\dfrac{{\partial n}}{{\partial T}}} \right|_p}$${\left. {\dfrac{{\partial n}}{{\partial p}}} \right|_T}$可以不等于零,由式(4)和式(5)可知,等温体弹模量KT(n)和热膨胀系数α(n)受${\left. {\dfrac{{\partial n}}{{\partial T}}} \right|_p}$${\left. {\dfrac{{\partial n}}{{\partial p}}} \right|_T}$的影响。本研究计算了(Mg0.5Fe0.5)CO3在MS态下的体积、热膨胀系数和声速等热力学参数,并与MgCO3、(Mg0.5Fe0.5)CO3的HS态的参数及前人的研究结果进行了对比。首先是体积的变化,由于$V(n) = {\left. {\dfrac{{\partial G(n)}}{{\partial p}}} \right|_T} = {\left. {\dfrac{{\partial G(n)}}{{\partial p}}} \right|_{T,n}} + {\left. {\dfrac{{\partial G(n)}}{{\partial n}}} \right|_{T,p}}{\left. {\dfrac{{\partial n}}{{\partial p}}} \right|_T}$,在平衡时,${\left. {\dfrac{{\partial G(n)}}{{\partial n}}} \right|_{T,p}} = 0$,因此,(Mg0.5Fe0.5)CO3在MS态的体积为V(n) = nVLS(p,T) + (1-n)VHS(p,T),其中VLSVHS分别为(Mg0.5Fe0.5)CO3在LS态与HS态时的体积。由于在相同压力下,(Mg0.5Fe0.5)CO3的HS态体积比LS态体积大,因此,由图6(a)可以看出,在压力转变区间内,MS态的体积比HS态体积小,与前人的理论和实验结果是相符的。

      图  6  MS态的体积V(a)、热膨胀系数α(b)、速度v(c)随压力的变化关系

      Figure 6.  Volume V (a), thermal expansion coefficient α (b), velocity v (c) of MS state as the function of pressure

      图6(b)图6(c)分别给出了热膨胀系数α和速度v随压力的变化关系。由MS态时热膨胀系数的计算公式(4)可知,α主要受到LS态的占比n对温度的一阶导数${\left. {\dfrac{{\partial n}}{{\partial T}}} \right|_p}$的影响,其中α${\left. {\dfrac{{\partial n}}{{\partial T}}} \right|_p}$成反比,由图5(a)可知,随着HS态向LS态发生自旋转变,在自旋共存区间内,${\left. {\dfrac{{\partial n}}{{\partial T}}} \right|_p}$的变化趋势为先突然减小然后增大,因此热膨胀系数α的变化趋势为先突然增大然后减小,这与体积的变化趋势不同;由式(5)、式(6)、式(8)可知,MS态的速度主要由LS态的占比n对压强的一阶导数${\left. {\dfrac{{\partial n}}{{\partial p}}} \right|_T}$决定,且v${\left. {\dfrac{{\partial n}}{{\partial p}}} \right|_T}$成反比,由图5(b)可知,在自旋共存区间内,${\left. {\dfrac{{\partial n}}{{\partial p}}} \right|_T}$的变化趋势为先突然增大后减小,因此速度v在自旋共存区间的变化趋势为先突然减小然后增大。在300 K时,MS态(Mg0.5Fe0.5)CO3在约40 GPa时产生异常变化,而在1 500 K时,这类异常变化出现在约50 GPa,该变化趋势与${\left. {\dfrac{{\partial n}}{{\partial T}}} \right|_p}$${\left. {\dfrac{{\partial n}}{{\partial p}}} \right|_T}$的变化一致,表明热力学参数的变化受n的影响。αVv由于自旋转变引起的异常随温度的增大向高压方向移动,并且变得更加平滑,但在2 500 K左右时,αVv的异常仍然比较明显。本研究计算结果与前人的实验研究均符合较好,在相同温度下,仅峰的位置存在较小的偏差,这是由于实验中Fe2+的浓度较高,其自旋转变的压力也随之相应地提高。(Mg0.5Fe0.5)CO3的MS态体积和速度的异常变化也会影响其在地球深部条件下的热传导特征。

    • 采用第一性原理计算方法研究了含铁菱镁矿(Mg0.5Fe0.5)CO3中铁的掺入以及自旋转变对其热力学性质的影响,得出以下主要结论。

      (1)含铁后菱镁矿晶体的体积发生明显变化,HS态在低温端体积增大,而在高温端体积比不含铁时减小,LS态的体积明显减小,表明铁的自旋转变会导致含铁菱镁矿晶体体积减小。

      (2)含铁后菱镁矿的热膨胀系数显著减小,但是铁的自旋转变导致热膨胀系数有所增大,但仍然比不含铁时要小;格临爱森常数与体变模量的变化是热膨胀系数变化的主要决定因素。

      (3)含铁菱镁矿的HS态与LS态两种自旋共存时,${\left. {\dfrac{{\partial n}}{{\partial T}}} \right|_p}$${\left. {\dfrac{{\partial n}}{{\partial p}}} \right|_T}$在自旋共存区间内不等于零,并在自旋过渡区域出现异常变化,此特征会导致含铁菱镁矿的体积在自旋共存压力区间突然减小至接近LS态时的体积;而热膨胀系数、速度也因受${\left. {\dfrac{{\partial n}}{{\partial T}}} \right|_p}$${\left. {\dfrac{{\partial n}}{{\partial p}}} \right|_T}$的影响,分别在自旋转变过程中出现异常增大峰和异常减小峰,并且这些异常变化峰随着温度的升高向高压方向移动。

      若要进一步明确含铁菱镁矿的热力学性质对地幔在碳循环中的作用,还需全面考虑地幔中其他矿物对菱镁矿的影响,如铁方镁石和钙钛矿等在地球深部同时存在时所表现出的物理化学性质。

参考文献 (40)

目录

    /

    返回文章
    返回