Fe-3.24%Si的状态方程和声速的第一性原理计算:地球内核Si元素的约束

李佩芸 黄海军 李艳丽

引用本文:
Citation:

Fe-3.24%Si的状态方程和声速的第一性原理计算:地球内核Si元素的约束

    作者简介: 李佩芸(1994-),女,硕士研究生,主要从事第一性原理计算研究. E-mail:lipeiyun@whut.edu.cn;
    通讯作者: 李艳丽, liyanli128@whut.edu.cn
  • 中图分类号: O521.2

First-Principles Calculations of the Equation of State and Sound Velocity of Fe-3.24%Si: Implications for the Composition of Earth’s Inner Core

    Corresponding author: LI Yanli, liyanli128@whut.edu.cn
  • CLC number: O521.2

  • 摘要: 硅(Si)被认为是地球内核的主要轻元素,但其在内地核中的含量仍然存在争议。为了探索内地核中Si的含量,应用第一性原理方法对Fe-3.24%Si(Si的质量分数为3.24%)进行了研究。构造了4种Fe-3.24%Si的超晶胞,研究了不同的晶胞大小和自旋对优化结构的影响。结果表明:在100 GPa以上,自旋对Fe-3.24%Si的密度无影响;而在100 GPa以下,考虑自旋时的计算结果更接近实验值。基于0 K下的声速、状态方程和相关热力学参数,计算了Fe-3.24%Si在内地核条件下的密度和声速。研究发现:Fe-3.24%Si的密度低于纯铁的密度,略高于内地核的密度;纵波声速及剪切波声速与纯铁的声速很接近,但均明显高于内地核声速,因此排除了内地核含有大量Si元素的可能性。
  • 图 1  hcp-Fe-3.24%Si的结构

    Figure 1.  Structure of hcp-Fe-3.24%Si

    图 2  单位体积内的能量e随压力的变化

    Figure 2.  Variation of the energy per unit volume e versus pressure

    图 3  Fe-3.24%Si的压力随密度的变化

    Figure 3.  Calculated density vs. pressure for Fe-3.24%Si

    图 4  0 K状态下Fe-3.24%Si的声速与密度的关系

    Figure 4.  Sound velocity of Fe-3.24%Si versus density at 0 K

    图 5  内地核条件下Fe-3.24%Si的密度、纵波声速vp、体波声速vb和横波声速vs与铁和内地核的数据比较

    Figure 5.  Comparison of the density, longitudinal sound velocity vp, bulk sound velocity vb and shear velocity vs of Fe-3.24%Si with those of Fe and the data of the inner core

  • [1] RINGWOOD A E. On the chemical evolution and densities of the planets [J]. Geochimica et Cosmochimica Acta, 1959, 15(4): 257–283. doi: 10.1016/0016-7037(59)90062-6
    [2] BIRCH F. Density and composition of mantle and core [J]. Journal of Geophysical Research, 1964, 69(20): 4377–4388. doi: 10.1029/JZ069i020p04377
    [3] TAKAFUJI N, HIROSE K, MITOME M, et al. Solubilities of O and Si in liquid iron in equilibrium with (Mg, Fe)SiO3 perovskite and the light elements in the core [J]. Geophysical Research Letters, 2005, 32(6).
    [4] FISCHER R A, CAMPBELL A J, REAMAN D M, et al. Phase relations in the Fe-FeSi system at high pressures and temperatures [J]. Earth and Planetary Science Letters, 2013, 373: 54–64. doi: 10.1016/j.jpgl.2013.04.035
    [5] FISCHER R A, CAMPBELL A J, CARACAS R, et al. Equations of state in the Fe-FeSi system at high pressures and temperatures [J]. Journal of Geophysical Research: Solid Earth, 2014, 119(4): 2810–2827. doi: 10.1002/2013JB010898
    [6] TATENO S, KUWAYAMA Y, HIROSE K, et al. The structure of Fe-Si alloy in Earth’s inner core [J]. Earth and Planetary Science Letters, 2015, 418: 11–19. doi: 10.1016/j.jpgl.2015.02.008
    [7] OZAWA H, HIROSE K, YONEMITSU K, et al. High-pressure melting experiments on Fe-Si alloys and implications for silicon as a light element in the core [J]. Earth and Planetary Science Letters, 2016, 456: 47–54. doi: 10.1016/j.jpgl.2016.08.042
    [8] KNITTLE E, JEANLOZ R. Earth’s core-mantle boundary: results of experiments at high pressures and temperatures [J]. Science, 1991, 251(5000): 1438–1443. doi: 10.1126/science.251.5000.1438
    [9] DUBROVINSKY L, DUBROVINSKAIA N, LANGENHORST F, et al. Iron-silica interaction at extreme conditions and the electrically conducting layer at the base of Earth’s mantle [J]. Nature, 2003, 422(6927): 58. doi: 10.1038/nature01422
    [10] LIN J F, CAMPBELL A J, HEINZ D L, et al. Static compression of iron-silicon alloys: implications for silicon in the Earth’s core [J]. Journal of Geophysical Research: Solid Earth, 2003, 108(B1).
    [11] ASANUMA H, OHTANI E, SAKAI T, et al. Static compression of Fe0.83Ni0.09Si0.08 alloy to 374 GPa and Fe0.93Si0.07 alloy to 252 GPa: implications for the Earth’s inner core [J]. Earth and Planetary Science Letters, 2011, 310(1/2): 113–118. doi: 10.1016/j.jpgl.2011.06.034
    [12] BADRO J, FIQUET G, GUYOT F, et al. Effect of light elements on the sound velocities in solid iron: implications for the composition of Earth’s core [J]. Earth and Planetary Science Letters, 2007, 254(1/2): 233–238.
    [13] ANTONANGELI D, SIEBERT J, BADRO J, et al. Composition of the Earth’s inner core from high-pressure sound velocity measurements in Fe-Ni-Si alloys [J]. Earth and Planetary Science Letters, 2010, 295(1/2): 292–296.
    [14] MAO Z, LIN J F, LIU J, et al. Sound velocities of Fe and Fe-Si alloy in the Earth’s core [J]. Proceedings of the National Academy of Sciences, 2012, 109(26): 10239–10244. doi: 10.1073/pnas.1207086109
    [15] LIU J, LIN J F, ALATAS A, et al. Seismic parameters of hcp-Fe alloyed with Ni and Si in the Earth’s inner core [J]. Journal of Geophysical Research: Solid Earth, 2016, 121(2): 610–623. doi: 10.1002/2015JB012625
    [16] SAKAIRI T, SAKAMAKI T, OHTANI E, et al. Sound velocity measurements of hcp Fe-Si alloy at high pressure and high temperature by inelastic X-ray scattering [J]. American Mineralogist, 2018, 103(1): 85–90. doi: 10.2138/am-2018-6072
    [17] ANTONANGELI D, MORARD G, PAOLASINI L, et al. Sound velocities and density measurements of solid hcp-Fe and hcp-Fe-Si (9 wt.%) alloy at high pressure: constraints on the Si abundance in the Earth’s inner core [J]. Earth and Planetary Science Letters, 2018, 482: 446–453. doi: 10.1016/j.jpgl.2017.11.043
    [18] TSUCHIYA T, FUJIBUCHI M. Effects of Si on the elastic property of Fe at Earth’s inner core pressures: first principles study [J]. Physics of the Earth and Planetary Interiors, 2009, 174(1): 212–219.
    [19] CÔTÉ A S, VOČADLO L, DOBSON D P, et al. Ab initio lattice dynamics calculations on the combined effect of temperature and silicon on the stability of different iron phases in the Earth’s inner core [J]. Physics of the Earth and Planetary Interiors, 2010, 178(1/2): 2–7.
    [20] MARTORELL B, WOOD I G, BRODHOLT J, et al. The elastic properties of hcp-Fe1− xSi x at Earth’s inner-core conditions [J]. Earth and Planetary Science Letters, 2016, 451: 89–96. doi: 10.1016/j.jpgl.2016.07.018
    [21] HOHENBERG P, KOHN W. Inhomogeneous electron gas [J]. Physical Review, 1964, 136(3B): B864. doi: 10.1103/PhysRev.136.B864
    [22] PERDEW J P. Exchange and correlation in atoms, molecules, and solids: the density functional picture [M]//Electron Correlations and Materials Properties. Boston: Springer, 1999: 287–298.
    [23] GROSS E K U, DREIZLER R M. Density functional theory: an approach to the quantum many-body problem [M]. Berlin: Springer, 1990.
    [24] KOHN W, SHAM L J. Quantum density oscillations in an inhomogeneous electron gas [J]. Physical Review, 1965, 137(6A): A1697. doi: 10.1103/PhysRev.137.A1697
    [25] LANGRETH D C, PERDEW J P. Theory of nonuniform electronic systems. I. analysis of the gradient approximation and a generalization that works [J]. Physical Review B, 1980, 21(12): 5469. doi: 10.1103/PhysRevB.21.5469
    [26] PERDEW J P, CHEVARY J A, VOSKO S H, et al. Atoms, molecules, solids, and surfaces: applications of the generalized gradient approximation for exchange and correlation [J]. Physical Review B, 1992, 46(11): 6671. doi: 10.1103/PhysRevB.46.6671
    [27] 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. doi: 10.1088/0953-8984/14/11/301
    [28] VOIGT W. The relation between the two elastic moduli of isotropic materials [J]. Annals of Physics (Leipzig), 1889, 33: 573.
    [29] REUSS A. Calculation of the flow limits of mixed crystals on the basis of the plasticity of monocrystals [J]. Zeitschrift für Angewandte Mathematik und Mechanik, 1929, 9: 49–58. doi: 10.1002/zamm.19290090104
    [30] HILL R. The elastic behaviour of a crystalline aggregate [J]. Proceedings of the Physical Society Section A, 1952, 65(5): 349. doi: 10.1088/0370-1298/65/5/307
    [31] WANG C S, KLEIN B M, KRAKAUER H. Theory of magnetic and structural ordering in iron [J]. Physical Review Letters, 1985, 54(16): 1852. doi: 10.1103/PhysRevLett.54.1852
    [32] ASADA T, TERAKURA K. Cohesive properties of iron obtained by use of the generalized gradient approximation [J]. Physical Review B, 1992, 46(20): 13599. doi: 10.1103/PhysRevB.46.13599
    [33] COHEN R E, MUKHERJEE S. Non-collinear magnetism in iron at high pressures [J]. Physics of the Earth and Planetary Interiors, 2004, 143: 445–453.
    [34] BROWN J M, FRITZ J N, HIXSON R S. Hugoniot data for iron [J]. Journal of Applied Physics, 2000, 88(9): 5496–5498. doi: 10.1063/1.1319320
    [35] 冯磊. 高压下温度对Fe-8.6Si声速的影响 [D]. 武汉: 武汉理工大学, 2017: 72–83.
    FENG L. Effect of temperature on Fe-8.6Si sound velocity at high pressure [D]. Wuhan: Wuhan University of Technology, 2017: 72–83.
    [36] 经福谦. 实验物态方程导引 [M]. 2版. 北京: 科学出版社, 1999: 188–197.
    JING F Q. Introduction to experimental equation of state [M]. 2nd ed. Beijing: Science Press, 1999: 188–197.
    [37] BROWN J M, MCQUEEN R G. Phase transitions, Grüneisen parameter, and elasticity for shocked iron between 77 GPa and 400 GPa [J]. Journal of Geophysical Research: Solid Earth, 1986, 91(B7): 7485–7494. doi: 10.1029/JB091iB07p07485
    [38] BONESS D A, BROWN J M, MCMAHAN A K. The electronic thermodynamics of iron under Earth core conditions [J]. Physics of the Earth and Planetary Interiors, 1986, 42(4): 227–240. doi: 10.1016/0031-9201(86)90025-7
    [39] FEI Y, MURPHY C, SHIBAZAKI Y, et al. Thermal equation of state of hcp-iron: constraint on the density deficit of Earth’s solid inner core [J]. Geophysical Research Letters, 2016, 43(13): 6837–6843. doi: 10.1002/2016GL069456
    [40] ANDERSON O L. The power balance at the core-mantle boundary [J]. Physics of the Earth and Planetary Interiors, 2002, 131(1): 1–17. doi: 10.1016/S0031-9201(02)00009-2
    [41] BIRCH F. Elasticity and constitution of the Earth’s interior [J]. Journal of Geophysical Research, 1952, 57(2): 227–286. doi: 10.1029/JZ057i002p00227
    [42] HIROSE K, LABROSSE S, HERNLUND J. Composition and state of the core [J]. Annual Review of Earth and Planetary Sciences, 2013, 41: 657–691. doi: 10.1146/annurev-earth-050212-124007
    [43] ZHANG Y, SEKINE T, LIN J F, et al. Shock compression and melting of an Fe-Ni-Si alloy: implications for the temperature profile of the Earth’s core and the heat flux across the core-mantle boundary [J]. Journal of Geophysical Research: Solid Earth, 2018, 123(2): 1314–1327. doi: 10.1002/2017JB014723
    [44] ANTONANGELI D, KOMABAYASHI T, OCCELLI F, et al. Simultaneous sound velocity and density measurements of hcp iron up to 93 GPa and 1100 K: an experimental test of the Birch’s law at high temperature [J]. Earth and Planetary Science Letters, 2012, 331: 210–214.
    [45] ANTONANGELI D, OHTANI E. Sound velocity of hcp-Fe at high pressure: experimental constraints, extrapolations and comparison with seismic models [J]. Progress in Earth and Planetary Science, 2015, 2(1): 3. doi: 10.1186/s40645-015-0034-9
    [46] LIN J F, STURHAHN W, ZHAO J, et al. Sound velocities of hot dense iron: Birch’s law revisited [J]. Science, 2005, 308(5730): 1892–1894. doi: 10.1126/science.1111724
    [47] SAKAMAKI T, OHTANI E, FUKUI H, et al. Constraints on Earth’s inner core composition inferred from measurements of the sound velocity of hcp-iron in extreme conditions [J]. Science Advances, 2016, 2(2): e1500802. doi: 10.1126/sciadv.1500802
    [48] CHEN B, LAI X, LI J, et al. Experimental constraints on the sound velocities of cementite Fe3C to core pressures [J]. Earth and Planetary Science Letters, 2018, 494: 164–171. doi: 10.1016/j.jpgl.2018.05.002
    [49] GAO L, CHEN B, WANG J, et al. Pressure-induced magnetic transition and sound velocities of Fe3C: implications for carbon in the Earth’s inner core [J]. Geophysical Research Letters, 2008, 35(17).
  • [1] 李芳菲崔启良李敏周强邹广田 . 高温高压下液态水声速的研究不同状态方程的准确性验证. 高压物理学报, doi: 10.11858/gywlxb.2008.03.010
    [2] 龚自正华劲松经福谦谢鸿森郭捷 . 高压下声速温度系数的一种新算法. 高压物理学报, doi: 10.11858/gywlxb.2000.04.001
    [3] 顾惠成李凤英 . 有磁和无磁状态下YIG和YAG晶体高压下的声速和衰减. 高压物理学报, doi: 10.11858/gywlxb.1993.02.010
    [4] 王贵朝 . 高压下材料弹性声速的一种经验表达式. 高压物理学报, doi: 10.11858/gywlxb.1988.01.014
    [5] 张其黎赵艳红马桂存张弓木 . 铀高压状态方程的第一原理研究. 高压物理学报, doi: 10.11858/gywlxb.2016.01.005
    [6] 郝爱民周铁军朱岩刘鑫 . 高压下TiN的结构转变、弹性和热力学性质的第一性原理计算. 高压物理学报, doi: 10.11858/gywlxb.2012.04.006
    [7] 耿华运孙毅向士凯 . 复杂材料高压性质的计算与模拟:基于第一性原理方法的部分进展. 高压物理学报, doi: 10.11858/gywlxb.20190710
    [8] 杨晓翠赵玉伟高忠明刘鑫张立新王晓明郝爱民 . 高压下CaF2结构相变和光学性质的第一性原理计算. 高压物理学报, doi: 10.11858/gywlxb.2010.03.011
    [9] 赵锋 . 两种状态方程参数的计算. 高压物理学报, doi: 10.11858/gywlxb.2014.05.004
    [10] 陈俊祥耿华运 . 多孔材料温压状态方程计算简要评述. 高压物理学报, doi: 10.11858/gywlxb.20190767
    [11] 杨晓翠赵衍辉罗香怡肖俊平 . 镧系氮化物压致结构相变的第一性原理计算. 高压物理学报, doi: 10.11858/gywlxb.2014.05.006
    [12] 徐济安 . 一个等温状态方程(Ⅴ)不同温度的高压状态方程. 高压物理学报, doi: 10.11858/gywlxb.1988.03.003
    [13] 徐济安 . 一个等温状态方程(Ⅳ)理想固体和绝对压强坐标. 高压物理学报, doi: 10.11858/gywlxb.1988.01.004
    [14] 薛富增李招宁潘守甫 . 铁热稠密物质状态方程的相对论Hartree-Fock-Slater计算. 高压物理学报, doi: 10.11858/gywlxb.1989.01.008
    [15] 杨晶巫翔秦善 . (Fe0.03Ni0.97)8(Si0.79P0.21)3的等温状态方程研究. 高压物理学报, doi: 10.11858/gywlxb.2011.03.013
    [16] 陈湘涛李茂生陈栋泉 . 一种实用的弹塑性流体状态方程(Ⅰ). 高压物理学报, doi: 10.11858/gywlxb.1988.01.006
    [17] 李茂生陈栋泉 . 一种实用的弹塑性流体状态方程(Ⅱ). 高压物理学报, doi: 10.11858/gywlxb.1988.02.009
    [18] 杨金科龚自正邓力维张莉费英伟 . 顽火辉石(Mg0.92,Fe0.08)SiO3的冲击相变和高压状态方程及其地球物理意义. 高压物理学报, doi: 10.11858/gywlxb.2007.01.008
    [19] 张其黎赵艳红马桂存 . 金高压物态方程的第一原理研究. 高压物理学报, doi: 10.11858/gywlxb.2014.01.003
    [20] 冯磊黄海军冷春蔚杨刚 . 高压下Fe92.5O2.2S5.3的熔化温度. 高压物理学报, doi: 10.11858/gywlxb.2017.06.004
  • 加载中
图(5)
计量
  • 文章访问数:  186
  • 阅读全文浏览量:  207
  • PDF下载量:  3
出版历程
  • 收稿日期:  2019-05-22
  • 录用日期:  2019-06-26
  • 网络出版日期:  2019-09-25

Fe-3.24%Si的状态方程和声速的第一性原理计算:地球内核Si元素的约束

    作者简介:李佩芸(1994-),女,硕士研究生,主要从事第一性原理计算研究. E-mail:lipeiyun@whut.edu.cn
    通讯作者: 李艳丽, liyanli128@whut.edu.cn
  • 武汉理工大学理学院,湖北 武汉 430070

摘要: 硅(Si)被认为是地球内核的主要轻元素,但其在内地核中的含量仍然存在争议。为了探索内地核中Si的含量,应用第一性原理方法对Fe-3.24%Si(Si的质量分数为3.24%)进行了研究。构造了4种Fe-3.24%Si的超晶胞,研究了不同的晶胞大小和自旋对优化结构的影响。结果表明:在100 GPa以上,自旋对Fe-3.24%Si的密度无影响;而在100 GPa以下,考虑自旋时的计算结果更接近实验值。基于0 K下的声速、状态方程和相关热力学参数,计算了Fe-3.24%Si在内地核条件下的密度和声速。研究发现:Fe-3.24%Si的密度低于纯铁的密度,略高于内地核的密度;纵波声速及剪切波声速与纯铁的声速很接近,但均明显高于内地核声速,因此排除了内地核含有大量Si元素的可能性。

English Abstract

  • 基于地球物理、地球化学和冶金学的研究,现已确定地核由铁和少量轻元素组成,如氧(O)、硫(S)、硅(Si)、碳(C)和氢(H)[1]。由于Si在地球中的含量非常高[1-2],高温高压下在铁中的溶解度也很高[3-7],并且铁与硅酸盐反应会生成硅铁合金[8-9],因此Si被认为是地核中的主要轻元素。

    目前,有关Si元素在内地核中的含量(以下如无特别说明,均为质量分数)仍然存在争议。例如:Lin等[10]测量了300 K下Fe-Si体系的状态方程,认为内核中存在约4%的Si;Asanuma等[11]分别测量了hcp-Fe0.93Si0.07在374 GPa和hcp-Fe0.83Ni0.09Si0.08在252 GPa压力下的状态方程,并由此预测内地核中含有3.9%~5.8%的Si;基于实验数据和第一性原理计算结果,Fischer等[5]认为内地核中Si的含量为6%~8%;Tateno等[6]根据实验测量结果推测内核中含有约7%的Si。

    地核中的候选轻元素除了满足密度约束条件外,还应满足声速的约束条件。2007年,Badro等[12]采用非弹性X射线散射技术(Inelastic X-ray Scattering,IXS)对FeSi在78 GPa下的声速进行研究,预测内地核中Si的含量为2.8%。Antonangeli等[13]使用相同的技术对Fe0.89Ni0.04Si0.07在108 GPa下的声速进行研究,认为内地核中Si的含量为1%~2%。2012年,Mao等[14]采用IXS和X射线衍射(X-Ray Diffraction,XRD)技术对hcp-Fe和hcp-Fe0.85Si0.15的声速进行研究,估测内地核中Si的含量为8%。Liu等[15]分别测量了206 GPa以内hcp-Fe和136 GPa以内hcp-Fe86.8Ni8.6Si4.6的声速,认为内地核中Si的含量为2.3%。最近,Sakairi等[16]通过高压实验研究了hcp-Fe0.89Si0.11(Fe-6%Si)的声速,指出内核中存在3%~6%的Si。Antonangeli等[17]测量了300 K下hcp-Fe和hcp-Fe-9%Si在高压下的声速,认为内核的Si不超过(3±2)%。受实验技术的限制,上述声速测量所能达到的温度和压强条件远小于内地核的温压环境,因此在外推过程中不可避免地引入了误差。

    基于密度泛函的第一性原理计算方法将不受上述条件的约束。Tsuchiya等[18]应用第一性原理方法计算了Fe-3.24%Si在0 K时的弹性常数,发现fcc (面心立方)结构和hcp(六方密堆积)结构的密度-速度关系与内地核数据一致。同样,运用第一性原理计算方法,Côté等[19]发现Fe-Si合金的结构与Si的含量有关:当Si的含量小于7%时,Fe-Si合金的结构为hcp结构;当Si的含量超过7%且温度高于4 000 K时,Fe-Si合金为fcc结构。Martorell等[20]计算了360 GPa时不同温度下的声速,认为Si不可能是内地核中的主要轻元素。综上所述,内地核中Si元素的含量还存在较大的争议,需要进一步开展理论和实验研究。

    • 本研究采用第一性原理计算方法对Fe-3.24%Si(Si的质量分数为3.24%)的状态方程和声速进行研究。计算中,将密度泛函理论(DFT)[21]的第一性原理计算方法与赝势平面波法相结合,交换关联能(XC)[22-23]采用Perdew-Burke-Ernzerh(PBE)的广义梯度近似(GGA)[24-25]。采用超软赝势(USPP)[26],其中截断能为380 eV,k点的设定为7×7×4。能量收敛精度为1.0×10–5 eV/atom,原子间作用力的收敛精度为0.03 eV/Å,原子间应力收敛精度为0.05 GPa,原子的最大允许位移为0.001 Å。基于上述参数,计算了0 K时不同压力下hcp-Fe-3.24%Si的密度和弹性常数。

      本研究基于应力和应变的关系获得了弹性常数[27]。根据晶体弹性力学相关理论,应力和应变之间的关系为

      $ \left[ {\begin{array}{*{20}{c}} {{\sigma _1}} \\ {{\sigma _2}} \\ {{\sigma _3}} \\ {{\sigma _4}} \\ {{\sigma _5}} \\ {{\sigma _6}} \end{array}} \right] = \left[\begin{array}{ccccccccc} {{c_{11}}}&{{c_{12}}}&{{c_{13}}}&{{c_{14}}}&{{c_{15}}}&{{c_{16}}}\\ {{c_{21}}}&{{c_{22}}}&{{c_{23}}}&{{c_{24}}}&{{c_{25}}}&{{c_{26}}}\\ {{c_{31}}}&{{c_{32}}}&{{c_{33}}}&{{c_{34}}}&{{c_{35}}}&{{c_{36}}}\\ {{c_{41}}}&{{c_{42}}}&{{c_{43}}}&{{c_{44}}}&{{c_{45}}}&{{c_{46}}}\\ {{c_{51}}}&{{c_{52}}}&{{c_{53}}}&{{c_{54}}}&{{c_{55}}}&{{c_{56}}}\\ {{c_{61}}}&{{c_{62}}}&{{c_{63}}}&{{c_{64}}}&{{c_{65}}}&{{c_{66}}} \end{array}\right]\left[ {\begin{array}{*{20}{c}} {{\varepsilon _1}} \\ {{\varepsilon _2}} \\ {{\varepsilon _3}} \\ {{\varepsilon _4}} \\ {{\varepsilon _5}} \\ {{\varepsilon _6}} \end{array}} \right] $

      弹性常数矩阵有36个弹性常数。根据晶体的对称性,弹性常数矩阵可以简化为

      ${{E}} = \left[\begin{array}{ccccccccc} {{c_{11}}}&{{c_{12}}}&{{c_{13}}}&0&0&0\\ {{c_{21}}}&{{c_{22}}}&{{c_{23}}}&0&0&0\\ {{c_{31}}}&{{c_{32}}}&{{c_{33}}}&0&0&0\\ 0&0&0&{{c_{44}}}&0&0\\ 0&0&0&0&{{c_{55}}}&0\\ 0&0&0&0&0&{{c_{66}}} \end{array}\right] $

      设零压p0下的晶体体积为V0。在本工作中,铁硅合金为六角立方结构。Voigt[28]认为体积模量KV和剪切模量GV可表示为

      ${K_{\rm V}} = \frac{{2\left( {{c_{11}} + {c_{12}}} \right) + 4{c_{13}} + {c_{33}}}}{9}$

      ${G_{\rm V}}{\rm{ = }}\frac{{ {M + 3{c_{11}} - 3{c_{12}} + 12{c_{44}} + {c_{66}}} }}{{30}}$

      Reuss[29]提出的体积模量KR和剪切模量GR分别为

      ${K_{\rm R}} = {C^2}/M$

      ${G_{\rm R}} = \dfrac{{15}}{{\dfrac{{18{K_{\rm V}}}}{{{C^2}}} + \dfrac{6}{{{c_{11}} - {c_{12}}}} + \dfrac{6}{{{c_{44}}}} + \dfrac{3}{{{c_{66}}}}}}$

      式中:MC分别表示为

      $M = {c_{11}} + {c_{12}} - 4{c_{13}} + 2{c_{33}}$

      ${C^2} = \left( {{c_{11}} + {c_{12}}} \right){c_{33}} - 2c_{13}^2$

      而Hill[30]认为晶体中应力和应变分布不均匀,其体积模量KH和剪切模量GH应为Voigt和Reuss建议的平均值,即

      ${K_{\rm{H}}} = \frac{1}{2}\left( {{K_{\rm{V}}} + {K_{\rm{R}}}} \right)$

      ${G_{\rm{H}}} = \frac{1}{2}\left( {{G_{\rm{V}}} + {G_{\rm{R}}}} \right)$

      在固体中,横波声速vs、纵波声速vp及体波声速vb与体积模量和剪切模量的关系为

      $ {v_{\rm{p}}} = \sqrt {\frac{{K + \dfrac{4}{3}G}}{\rho }},\;\;\;\;{v_{\rm{s}}} = \sqrt {\frac{G}{\rho }},\;\;\;\;{v_{\rm{b}}} = \sqrt {v_{\rm{p}}^2 - \frac{4}{3}v_{\rm{s}}^2} $

      式中:$ \rho$为密度。

    • 优化了4种不同结构的hcp-Fe-3.24%Si,通过比较能量,得到最接近实际值的结构(II-64-atoms)。基于第一性原理计算的状态方程和声速,计算了内核温压条件下的密度和声速,以此限定内核中Si的含量。

    • 考虑4种不同结构的hcp-Fe-3.24%Si(Fe0.937 5Si0.062 5),即2×2×2(16-atoms)、2×2×4(32-atoms)、4×2×4(I-64-atoms)和4×2×4(II-64-atoms),如图1所示,其中:4×2×4(I-64-atoms)的超晶胞是4个Si原子均匀分散在64个Fe原子内,4×2×4(II-64-atoms)结构是4个Si原子处于64个Fe原子中最分散的位置。这4种Fe-3.24%Si结构均在基态下进行了优化,并考虑了自旋的影响。根据每一单位体积的能量,可以选出最稳定的结构,如图2所示。当压强低于100 GPa时,II-64-atoms具有最低能量;当压强高于100 GPa时,自旋对II-64-atoms结构的能量无影响,因此单位体积内II-64-atoms具有最低的能量,即此结构最稳定。

      图  1  hcp-Fe-3.24%Si的结构

      Figure 1.  Structure of hcp-Fe-3.24%Si

      图  2  单位体积内的能量e随压力的变化

      Figure 2.  Variation of the energy per unit volume e versus pressure

    • 沿0 K等温线,用第一性原理方法计算得到状态方程,如图3所示,其中红色实心圆和黑色实心方块分别代表Fe-3.24%Si有自旋和无自旋计算结果。当压强低于100 GPa时,Fe-3.24%Si有自旋时的密度高于无自旋时的密度,并且更接近实验数据,表明低压时必须考虑电子自旋的影响。其他研究工作[31-33]也表明,当压强低于约60 GPa时,有自旋状态比无自旋状态更稳定。当压强高于100 GPa时,与Asanuma等[11]通过实验测量的Fe-3.4%Si的实验数据一致,并且有无自旋对计算结果无影响。将不考虑自旋的Fe-3.24%Si的计算结果用三阶Birch-Murnaghan状态方程(其中初始密度$ {\rho _0}$=(8.87±0.09) g/cm3)拟合

      图  3  Fe-3.24%Si的压力随密度的变化

      Figure 3.  Calculated density vs. pressure for Fe-3.24%Si

      $ p = \frac{3}{2}{K_0}\left[ {{{\left( {\frac{\rho }{{{\rho _0}}}} \right)}^{\textstyle\frac{7}{3}}} - {{\left( {\frac{\rho }{{{\rho _0}}}} \right)}^{\textstyle\frac{5}{3}}}} \right] \times \left\{ {1 + \frac{3}{4}\left( {{K'_0} - 4} \right)\left[ {{{\left( {\frac{\rho }{{{\rho _0}}}} \right)}^{\textstyle\frac{2}{3}}} - 1} \right]} \right\} $

      得到零压下的体积模量K0=(368±38) GPa,其压力导数为${K'_0}$=3.72±0.31。将计算结果外推至330 GPa,与Martorell等[20]采用第一性原理分子动力学模拟的Fe-3.24%Si结果一致,但却略高于Tsuchiya等[18]用第一性原理计算方法得到的结果。

      图3也给出了根据Hugoniot数据计算的Fe-3.24%Si的300 K状态方程,计算方法如下。Mie-Grüneisen方程表示为

      ${p_{300}}(V) = {p_{\rm{H}}}(V) - \frac{1}{V}\int_{300}^{{T_{\rm{H}}}} {{C_V}{\gamma _{{\rm{eff}}}}} {\rm{d}}T$

      式中:${p_{\rm{H}}}(V)$为Hugoniot状态方程。

      ${p_{\rm{H}}} = {\rho _0}C_0^2{{\left( {1 - {{{\rho _0}} / \rho }} \right)} / {{{\left[ {1 - \lambda \left( {1 - {{{\rho _0}} / \rho }} \right)} \right]}^2}}}$

      式中:C0$ \lambda $为Hugoniot参数。到目前为止,还未见发表Fe-3.24%Si的Hugoniot数据。为此,由纯铁的Hugoniot参数C0=3.935 km/s和$ \lambda $=1.578,纯铁的初始密度$ {\rho _0}$=7.85 g/cm3[34],以及Fe-9%Si的Hugoniot参数C0=(4.603±0.101) km/s,$ \lambda $=1.505±0.037,Fe-9%Si的初始密度$ {\rho _0}$=7.386 g/cm3[35],根据可加性原理[36]得到Fe-3.24%Si的Hugoniot数据

      $V(p) = \sum\limits_i^{} {{w_i}V_{0i}^{}\left\{ {1 - \frac{1}{{2\lambda _i^2}}\left[ {\left( {2{\lambda _i} + \frac{{C_{0i}^2}}{{V_{0i}^{}p}}} \right) - \sqrt {{{\left( {2{\lambda _i} + \frac{{C_{0i}^2}}{{V_{0i}^{}p}}} \right)}^2} - 4\lambda _i^2} } \right]} \right\}} $

      式中:下标i代表i组分,wi是组分i的质量分数,V0为基态体积。根据(15)式,计算得到Fe-3.24%Si的Hugoniot参数为C0=4.182 km/s,$ \lambda $=1.551,以及0 K密度$ {\rho _0}$=7.668 g/cm3。(13)式中的冲击温度可以通过下式估算

      ${\rm{d}}T = - T\frac{{{\gamma _{{\rm{eff}}}}}}{V}{\rm{d}}V + \frac{1}{{2{C_V}}}\left[ {{p_{\rm{H}}}\left( V \right) + \left( {{V_0} - V} \right)\frac{{{\rm{d}}p_{\rm{H}}^{}\left( V \right)}}{{{\rm{d}}V}}} \right]{\rm{d}}V$

      ${C_V} = {C_{V{\rm{l}}}} + {C_{V{\rm{e}}}}$

      ${\gamma _{{\rm{eff}}}} = \left( {{\gamma _{\rm{l}}}{C_{V{\rm{l}}}} + {\gamma _{\rm{e}}}{C_{V{\rm{e}}}}} \right)/{C_V}$

      式中:CVl代表晶格对定容比热容的贡献,简称晶格比热,在高温下CVl=3R/$ \mu $,其中R为理想气体常数,摩尔质量$ \mu $=54.11 g/mol;$ {C_{V{\rm{e}}}} = {\beta _0}{\left( {{\rho _{0\varepsilon }}/\rho } \right)^\kappa }T$[37]为电子对定容比热容的贡献,简称电子比热。由于Fe-3.24%Si中Si的含量很小,因此假设其电子比热的参数和铁一致,即$ {\beta _0}$=0.091 J·kg–1·K–2$\kappa = 1.34$[38]。另外,在本研究中还假设Fe-3.24%Si晶格和电子的Grüneisen参数都近似等于铁的相关参数$ {\gamma _{\rm l}} = 1.74{({\rho _{0\varepsilon }}/\rho )^{0.78}}$[39]$ {\gamma _{\rm e}}=2$[40]。基于hcp-Fe-9%Si[6]在室温下的密度$ {{\rho _{0\varepsilon }}}$=7.578 g/cm3和hcp-Fe[39]的密度$ {{\rho _{0\varepsilon }}}$=8.269 g/cm3,通过内插法可估计出hcp-Fe-3.24%Si在室温下的密度$ {{\rho _{0\varepsilon }}}$=7.995 g/cm3。将这些参数代入(13)式,计算得到Fe-3.24%Si在室温下的状态方程,如图3中品红色虚线所示,与Fe-3.4%Si的静高压数据一致[11],也与100 GPa以上的第一性原理计算结果一致。

    • 由前面对比分析可知,当压强在100 GPa以上时,有无自旋对Fe-3.24%Si的状态方程无影响。考虑计算资源的问题,本研究假设自旋对声速亦无影响,在此仅计算II-64-atoms结构在高压零温下的声速。Fe-3.24%Si[18, 20]的纵波声速vp与密度呈线性关系(见图4

      图  4  0 K状态下Fe-3.24%Si的声速与密度的关系

      Figure 4.  Sound velocity of Fe-3.24%Si versus density at 0 K

      ${v_{\rm{p}}} = - \left( {2.08 \pm 0.43} \right) + \left( {1.10 \pm 0.04} \right)\rho $

      也就是说,Fe-3.24%Si的纵波声速vp在室温下满足Birch定律[41]。在低压下,所有Fe-3.24%Si的计算结果与实验测量的Fe-4.5%Ni-3.7%Si的纵波声速vp一致[13]。对于体波声vb,本研究通过第一性原理得到的计算结果与Tsuchiya等[18]和Martorell等[20]计算的体波速度vb比较吻合。基于之前计算Fe-3.24%Si状态方程的方法,还可以根据热力学方法估算体波速度

      $v_{\rm{b}}^2 = - {v^2}{\left( {\frac{{\partial p}}{{\partial V}}} \right)_T} + \gamma _{_{{\rm{eff}}}}^2{C_V}T$

      0 K下,根据(20)式计算的体波速度vb在低压下略高于根据(11)式得到的vb,并且两者的差异在误差范围内,高压下计算得到的vb一致。因此,可以用(20)式估算内核高温高压条件下的体波速度vb。由(19)式和(20)式,计算得到的横波声速vs如下

      $\frac{4}{3}v_{\rm{s}}^2 = v_{\rm{p}}^2 - v_{\rm{b}}^2$

      在误差范围内,由(21)式计算的横波声速与由(10)式得到的横波声速相近。根据状态方程和纵波声速的拟合结果,vbvs的相对误差分别约为11%和24%。

    • 根据0 K的状态方程和相应的热力学参数,可以通过下式计算内核条件下Fe-3.24%Si的密度

      ${p_T}(V) = {p_0}(V) + \frac{1}{V}\int_0^T {{C_V}{\gamma _{{\rm{eff}}}}} {\rm{d}}T$

      内地核温度剖面$ {T_{\rm{E}}} = {T_{{\rm{ICB}}}}{\left( {\rho /{\rho _{{\rm{ICB}}}}} \right)^{1.5}}$可根据内外核边界处的温度TICB=(5 400±300) K[42-43]计算。对于纵波声速,前面的计算结果表明Fe-Si体系的纵波声速在室温下满足Birch定律,在高温高压下是否依然呈线性关系还存在争论[44]。Antonangeli等[45]在低于93 GPa、不高于1 100 K的温压范围内测量了Fe的声速,结果表明Fe的纵波声速在1100 K下仍然满足Birch定律。然而,Lin等[46]、Mao等[14]和Sakamaki等[47]的研究结果表明Fe的纵波声速在高温下与密度不再呈线性关系,意味着Birch定律在高温条件下并不成立。2017年,冯磊[35]应用动高压实验技术研究了Fe-9%Si在216 GPa、4 420 K时的纵波声速,测量数据与300 K下Fe-9%Si[17]的纵波声速非常一致,表明Fe-9%Si的纵波声速在高温高压下仍然满足Birch定律。在内地核高温高压条件下,假设Fe-3.24%Si的纵波声速仍然满足Birch定律,则其体波声速和横波声速可由(20)式和(21)式计算,结果如图5所示。为了研究Si对Fe的密度和声速的影响,还计算了内核条件下铁的密度和声速。铁的密度可由(13)式计算,其体积模量为K0=(172.7±1.4) GPa,${K'_0}$=4.79±0.05,初始密度$ {\rho _0}$=8.269 g/cm3[39];铁的纵波声速$ {v_{\rm{p}}} = (1.206\rho - 4.000)$ km/s[17]。在内核条件下,Fe-3.24%Si的密度低于Fe的密度,略高于内核的密度,表明Si降低了铁的密度。Fe-3.24%Si的纵波声速与铁的纵波声速接近,表明Si对铁的纵波声速的影响很小。Martorell等[20]的计算结果也支持这一结论,尽管他们认为温度会降低Fe和Fe-Si合金的纵波声速。Fe-3.24%Si的体波速度与PREM数据一致,并且高于铁;而Fe-3.24%Si的横波声速略低于Fe,但是均明显高于内地核的横波声速。由此可见,Si对铁的纵波声速和横波声速的影响较小。由于纯铁的声速明显高于内地核的声速,因此内地核中不可能含有大量的Si,而应含有能明显降低Fe的密度和声速的轻元素,如C[48-49]

      图  5  内地核条件下Fe-3.24%Si的密度、纵波声速vp、体波声速vb和横波声速vs与铁和内地核的数据比较

      Figure 5.  Comparison of the density, longitudinal sound velocity vp, bulk sound velocity vb and shear velocity vs of Fe-3.24%Si with those of Fe and the data of the inner core

    • 研究了0 K下Fe-3.24%Si的状态方程,并将其与静态高压实验数据、理论计算和动态高压实验数据进行了比较,研究发现:当压强高于100 GPa时,自旋对计算结果没有影响;在低压条件下考虑自旋时,其计算结果更接近实验值。通过计算弹性常数,得到了Fe-3.24%Si的纵波声速、横波声速和体波声速,发现根据第一性原理计算的纵波声速与Fe-4.5%Ni-3.7%Si的实验测量结果一致,并且满足Birch定律。

      基于由第一性原理方法计算的状态方程和声速,利用热力学方法,计算了内地核温压环境下Fe-3.24%Si的密度和声速,并与内地核的密度和声速进行对比,发现:Fe-3.24%Si的密度高于PREM模型给出的内地核密度,但明显低于纯铁的密度;Fe-3.24%Si的纵波和横波声速接近Fe的声速,但是明显高于内地核的声速。由此排除了内地核含有大量Si元素的可能性。

参考文献 (49)

目录

    /

    返回文章
    返回