高压下Fe92.5O2.2S5.3的熔化温度

冯磊 黄海军 冷春蔚 杨刚

引用本文:
Citation:

高压下Fe92.5O2.2S5.3的熔化温度

    作者简介: 冯磊(1987-), 男, 硕士, 主要从事固体地球物理研究.E-mail:feng_lei@whut.edu.cn;
    通讯作者: 黄海军, hjhuang@whut.edu.cn
  • 基金项目: 国家重点基础研究发展计划(973计划) 2014CB845904
    中央高校基本科研业务费专项资金 2015Ⅲ035
    中央高校基本科研业务费专项资金 2014-Ⅶ-006
    国家自然科学基金优秀青年科学基金 41322028

  • 中图分类号: O521.2;P3

Melting Temperatures of Fe92.5O2.2S5.3 under High Pressure

    Corresponding author: HUANG Hai-Jun, hjhuang@whut.edu.cn ;
  • CLC number: O521.2;P3

  • 摘要: 采用反向碰撞法与光分析技术,测量了Fe92.5O2.2S5.3(质量分数比)在208 GPa下的声速,发现固态Fe92.5O2.2S5.3的纵波声速在144 GPa下开始减小,直到165 GPa完全转变为液态体波声速,表明样品的完全熔化温度为(3 500±400)K。将该熔化温度作为参考点,应用Lindeman定律并外推至地球内外核边界可知,Fe92.5O2.2S5.3的熔化温度为(5 000±400)K。通过比较Fe、Fe-O、Fe-S以及Fe-O-S的熔化温度,发现O元素对Fe熔化温度的影响很小,S元素对Fe熔化温度的降低与其含量成正比。如果外地核中S的质量分数为2%~6%,则地球内外核界面温度为5 000~5 400 K。
  • 图 1  反向碰撞法实验装置图(a)和反向碰撞法波系作用图(b)

    Figure 1.  Sample configuration (a) and wave interaction (b) for sound velocity measurements in the reverse impact method

    图 2  实验No.081225和No.081226的界面粒子速度历史

    Figure 2.  Particle velocity histories at the impact interface for experiments No.081225 and No.081226

    图 3  拉格朗日声速与粒子速度的关系

    Figure 3.  Lagrangian sound velocity vs.particle velocity

    图 4  采用光分析技术测量声速的实验装置

    Figure 4.  Experimental configuration for sound velocity measurements using optical analyzer technique

    图 5  光分析技术测量声速实验中波的传播

    Figure 5.  Wave interaction for sound velocity measurements in optical analyzer technique

    图 6  实验No.110505在198 GPa下的粒子速度历史

    Figure 6.  Particle velocity history for experiment No.110505 at 198 GPa

    图 7  实验No.110505在198 GPa下的时间差Δt与样品厚度的关系

    Figure 7.  Time intervals Δt vs.sample thickness for experiment No.110505 at 198 GPa

    图 8  Fe92.5O2.2S5.3的声速与冲击压强的关系

    Figure 8.  Sound velocity of Fe92.5O2.2S5.3 vs.shock pressure

    图 9  Fe、Fe-O、Fe-S、Fe-O-S体系在高压下的熔化温度和液相线温度随压强的变化(红色实线为本研究得到的Fe92.5O2.2S5.3的熔化温度。Fe:Ma等[40],标记为M04;Anzellini等[39],标记为A13。Fe-O-S体系:Terasaki等[28],标记为T11;Huang等[14],标记为H10。Fe-S体系:Kamada等[26],标记为K12。)

    Figure 9.  Melting temperature of Fe-O-S system compared with those of Fe, Fe-O and Fe-S (The red solid line shows the melting temperature of Fe92.5O2.2S5.3 from this study.The lines labeled as M04 and A13 represent the melting temperatures of Fe from Ref.[40] and Ref.[39] respectively.The dashed lines labeled as H10 and T11 represent the melting temperature of Fe90O8S2[14] and the liquidus and solidus temperature of Fe-O-S[28].K12 represent the melting relationships in the Fe-Fe3S system up to the outer core conditions[26].)

    图 10  高压下Fe-S与Fe-O-S体系液相线温度的比较

    Figure 10.  Comparison of liquid-phase temperature between Fe-S and Fe-O-S system under high pressure

    表 1  Fe、FeO、FeS和Fe-O-S体系参数

    Table 1.  Parameters for Fe, FeO, FeS and Fe-O-S system

    Material ρ0/(g/cm3) C0/(km/s) λ γ0 q β0/(J·kg-1·K-2) κ E0/(kJ)
    α-Fe 7.85[43] 3.935[43] 1.578[43]
    ε-Fe 8.298[44] 4.720* 1.523* 1.76* 0.76* 0.091[45] 1.34[45] 76.268 4[44]
    FeO(B1) 5.71[46] 5.83* 0.99* 57.796 7[47]
    FeO(B8) 6.05[46, 48] 4.486* 1.699* 1.8[48] 1[48] 0[49] 0[49]
    FeS 4.602[50] 2.947[50] 1.578[50] 1.54[50] 1[50] 0.25[51] 1.34[51] 33.707 5[52]
    Fe90O8S2 6.69[14] 3.97[14] 1.58[14] 1.85[14] 0.87[14] 0.075[14] 1.39[14] 67.315
    Fe92.5O2.2S5.3 6.88 3.71 1.61 1.82 0.85 0.119 1.407 68.258
    Note:The asterisk * represents the fitted results from the Hugoniot data.
    下载: 导出CSV
  • [1] BIRCH F.Density and composition of mantle and core[J]. J Geophys Res, 1964, 69(20):4377-4388.
    [2] POIRIER J P.Light elements in the Earth's outer core:a critical review[J]. Phys Earth Planet Inter, 1994, 85(3/4):319-337.
    [3] HILLGREN V J, GESSMANN C K, LI J.Origin of the Earth and Moon[M]. New York:Springer, 2000:245-263.
    [4] HOLLAND H D, TUREKIAN K.Treatise on geochemistry[M]. Oxford:Elsevier, 2003:527-557.
    [5] MURTHY V R, HALL H T.The chemical composition of the Earth's core:possibility of sulphur in the core[J]. Phys Earth Planet Inter, 1970, 2(4):276-282.
    [6] USSELMAN T M.Experimental approach to the state of the core:Part Ⅰ, the liquidus relations of the Fe-rich portion of the Fe-Ni-S system from 30 to 100 kb[J]. Am J Sci, 1975, 275(3):278-290. doi: 10.2475/ajs.275.3.278
    [7] FEI Y W, BERTKA C M, FINGER L W.High-pressure iron-sulfur compound, Fe3S2, and melting relations in the Fe-FeS system[J]. Science, 1997, 275(5306):1621-1623. doi: 10.1126/science.275.5306.1621
    [8] FEI Y W, LI J, BERTKA C M, et al.Structure type and bulk modulus of Fe3S, a new iron-sulfur compound[J]. Am Mineral, 2000, 85(11/12):1830-1833.
    [9] LI J, FEI Y W, MAO H K, et al.Sulfur in the Earth's inner core[J]. Earth Planet Sci Lett, 2001, 193(3):509-514.
    [10] RINGWOOD A E.Composition of the core and implications for origin of the Earth[J]. Geochem J, 1977, 11(3):111-135.
    [11] URAKAWA S, KATO M, KUMAZAWA M.Experimental study on the phase relations in the system Fe-Ni-O-S up to 15 GPa[J]. High-Pressure Research in Mineral Physics:A Volume in Honor of Syun-iti Akimoto, 1987:95-111.
    [12] HELFFRICH G, KANESHIMA S.Seismological constraints on core composition from Fe-O-S liquid immiscibility[J]. Science, 2004, 306(5705):2239-2242.
    [13] HUANG H J, FEI Y W, CAI L C, et al.Evidence for an oxygen-depleted liquid outer core of the Earth[J]. Nature, 2011, 479(7374):513-516. doi: 10.1038/nature10621
    [14] HUANG H J, HU X J, JING F Q, et al.Melting behavior of Fe-O-S at high pressure:a discussion on the melting depression induced by O and S[J]. J Geophys Res, 2010, 115(B5).
    [15] OHTANI E, RINGWOOD A E, HIBBERSON W.Composition of the core, Ⅱ.effect of high pressure on solubility of FeO in molten iron[J]. Earth Planet Sci Lett, 1984, 71(1):94-103.
    [16] KATO T, RINGWOOD A E.Melting relationships in the system Fe-FeO at high pressures:implications for the composition and formation of the Earth's core[J]. Phys Chem Miner, 1989, 16(6):524-538.
    [17] RINGWOOD A E, HIBBERSON W.The system Fe-FeO revisited[J]. Phys Chem Miner, 1990, 17(4):313-319.
    [18] SEAGLE C T, HEINZ D L, CAMPBELL A J, et al.Melting and thermal expansion in the Fe-FeO system at high pressure[J]. Earth Planet Sci Lett, 2008, 265(3/4):655-665.
    [19] ANDRAULT D, BOLFAN-CASANOVA N, OHTAKA O, et al.Melting diagrams of Fe-rich alloys determined from synchrotron in situ measurements in the 15-23 GPa pressure range[J]. Phys Earth Planet Inter, 2009, 174(1):181-191.
    [20] WILLIAMS Q, JEANLOZ R.Melting relations in the iron-sulfur system at ultra-high pressures:implications for the thermal state of the Earth[J]. J Geophys Res B, 1990, 95(12):19299-19310.
    [21] BOEHLER R.Fe-FeS eutectic temperatures to 620 kbar[J]. Phys Earth Planet Inter, 1996, 96(2):181-186.
    [22] CHUDINOVSKIKH L, BOEHLER R.Eutectic melting in the system Fe-S to 44 GPa[J]. Phys Earth Planet Inter, 2007, 257(1/2):97-103.
    [23] CAMPBELL A J, SEAGLE C T, HEINZ D L, et al.Partial melting in the iron-sulfur system at high pressure:a synchrotron X-ray diffraction study[J]. Phys Earth Planet Inter, 2007, 162(1):119-128.
    [24] MORARD G, ANDRAULT D, GUIGNOT N, et al.In situ determination of Fe-Fe3S phase diagram and liquid structural properties up to 65 GPa[J]. Earth Planet Sci Lett, 2008, 272(3):620-626.
    [25] KAMADA S, TERASAKI H, OHTANI E, et al.Phase relationships of the Fe-FeS system in conditions up to the Earth's outer core[J]. Earth Planet Sci Lett, 2010, 294(1):94-100.
    [26] KAMADA S, OHTANI E, TERASAKI H, et al.Melting relationships in the Fe-Fe3S system up to the outer core conditions[J]. Earth Planet Sci Lett, 2012, 359:26-33.
    [27] TSUNO K, OHTANI E.Eutectic temperatures and melting relations in the Fe-O-S system at high pressures and temperatures[J]. Phys Chem Miner, 2008, 36(1):9-17.
    [28] TERASAKI H, KAMADA S, SAKAI T, et al.Liquidus and solidus temperatures of a Fe-O-S alloy up to the pressures of the outer core:implication for the thermal structure of the Earth's core[J]. Earth Planet Sci Lett, 2011, 304(3):559-564.
    [29] DUFFY T S, AHRENS T J.Compressional sound velocity, equation of state, and constitutive response of shock-compressed magnesium oxide[J]. J Geophys Res B, 1995, 100(1):529-542.
    [30] HU J B, ZHOU X M, TAN H, et al.Successive phase transitions of tin under shock compression[J]. Appl Phys Lett, 2008, 92(11):111905. doi: 10.1063/1.2898891
    [31] WENG J D, TAN H, WANG X, et al.Optical-fiber interferometer for velocity measurements with picosecond resolution[J]. Appl Phys Lett, 2006, 89(11):111101.
    [32] ASAY J R, CHHABILDAS L C.Determination of the shear strength of shock compressed 6061-T6 aluminum[M]//Shock Waves and High-Strain-Rate Phenomena in Metals.Springer, 1981: 417-431.
    [33] MITCHELL A C, NELLIS W J.Shock compression of aluminum, copper, and tantalum[J]. J Appl Phys, 1981, 52(5):3363-3374.
    [34] BROWN J M, SHANER J W.Rarefaction velocities in shocked tantalum and the high pressure melting point[M]//Shock Waves in Condensed Matter-1983.North Holland: Springer Science & Business Media, 1984: 91-94.
    [35] HUANG H J, JING F Q, CAI L C, et al.Grüneisen parameter along Hugoniot and melting temperature of ε-iron:a result from thermodynamic calculations[J]. Chin Phys Lett, 2005, 22(4):836-838.
    [36] WALLACE D C.Melting of elements[J]. Proc Roy Soc A, 1991, 433(1889):631-661. doi: 10.1098/rspa.1991.0068
    [37] BROWN J M, MCQUEEN R G.Phase transitions, Grüneisen parameter, and elasticity for shocked iron between 77 GPa and 400 GPa[J]. J Geophys Res B, 1986, 91(7):7485-7494.
    [38] ANDERSON O L.The power balance at the core-mantle boundary[J]. Phys Earth Planet Inter, 2002, 131(1):1-17.
    [39] ANZELLINI S, DEWAELE A, MEZOUAR M, et al.Melting of iron at Earth's inner core boundary based on fast X-ray diffraction[J]. Science, 2013, 340(6131):464-466.
    [40] MA Y, SOMAYAZULU M, SHEN G, et al.In situ X-ray diffraction studies of iron to Earth-core conditions[J]. Phys Earth Planet Inter, 2004, 143(1):455-467.
    [41] SHEN G, MAO H K, HEMLEY R J, et al.Melting and crystal structure of iron at high pressures and temperatures[J]. Geophys Res Lett, 1998, 25(3):373-376.
    [42] HUANG H J, JING F Q, CAI L C.Grüneisen parameter along Hugoniot and melting temperature of iron:a thermodynamic computational method[J]. Chin Phys Lett, 2005, 22(4):836-838.
    [43] BROWN J M, FRITZ J N, HIXSON R S.Hugoniot data for iron[J]. J Appl Phys, 2000, 88(9):5496-5498. doi: 10.1063/1.1319320
    [44] ANDERSON O L, DUBROVINSKY L, SAXENA S K, et al.Experimental vibrational Grüneisen ratio values for ε-iron up to 330 GPa at 300 K[J]. Geophys Res Lett, 2001, 28(2):399-402.
    [45] BONESS D A, BROWN J M, MCMAHAN A K.The electronic thermodynamics of iron under Earth core conditions[J]. Phys Earth Planet Inter, 1986, 42(4):227-240. doi: 10.1016/0031-9201(86)90025-7
    [46] JACKSON I, KHANNA S K, REVCOLEVSCHI A, et al.Elasticity, shear-mode softening and high-pressure polymorphism of wüstite (Fe1-xO)[J]. J Geophys Res B, 1990, 95(13):21671-21685.
    [47] ANDERSON O L.The Grüneisen parameter for iron at outer core conditions and the resulting conductive heat and power in the core[J]. Phys Earth Planet Inter, 1998, 109(3):179-197.
    [48] JEANLOZ R, AHRENS T J.Equations of state of FeO and CaO[J]. Geophys J Int, 1980, 62(3):505-528. doi: 10.1111/gji.1980.62.issue-3
    [49] STIXRUDE L, WASSERMAN E, COHEN R E.Composition and temperature of Earth's inner core[J]. J Geophys Res B, 1997, 102(11):24729-24739.
    [50] BROWN J M, AHRENS T J, SHAMPINE D L.Hugoniot data for pyrhotite and the Earth's core[J]. J Geophys Res B, 1984, 89(7):6041-6048.
    [51] ANDERSON W W, AHRENS T J.Shock temperature and melting in iron sulfides at core pressures[J]. J Geophys Res B, 1996, 101(3):5627-5642.
    [52] STEVENSON D J.Models of the Earth's core[J]. Science, 1981, 214(4521):611-619. doi: 10.1126/science.214.4521.611
    [53] DREIBUS G, PALME H.Cosmochemical constraints on the sulfur content in the Earth's core[J]. Geochimica et Cosmochimica Acta, 1996, 60(7):1125-1130.
    [54] MAHAN B, SIEBERT J, PRINGLE E A, et al.Elemental partitioning and isotopic fractionation of Zn between metal and silicate and geochemical estimation of the S content of the Earth's core[J]. Geochimica et Cosmochimica Acta, 2017, 196:252-270. doi: 10.1016/j.gca.2016.09.013
    [55] HIROSE K, LABROSSE S, HERNLUND J.Composition and state of the core[J]. Annu Rev Earth Planet Sci, 2013, 41(1):657-691.
  • [1] 龚自正华劲松经福谦谢鸿森郭捷 . 高压下声速温度系数的一种新算法. 高压物理学报, 2000, 14(4): 241-246 . doi: 10.11858/gywlxb.2000.04.001
    [2] 金柯吴强李加波周显明叶素华李俊 . 冲击加载下NaCl单晶高压声速与温度的同步测量. 高压物理学报, 2017, 31(6): 707-717. doi: 10.11858/gywlxb.2017.06.005
    [3] 张凌云戴诚达许灿华谭华 . 无氧铜冲击熔化温度的辐射法测量. 高压物理学报, 2005, 19(4): 365-370 . doi: 10.11858/gywlxb.2005.04.015
    [4] 戴诚达谭华 . 金属冲击温度的辐射法测量问题. 高压物理学报, 2006, 20(2): 113-121 . doi: 10.11858/gywlxb.2006.02.001
    [5] 谭叶肖元陆薛桃李俊金柯 . 镁铝合金的冲击熔化行为实验研究. 高压物理学报, 2019, 33(2): 020106-1-020106-6. doi: 10.11858/gywlxb.20190729
    [6] 王汝菊王积方 . 高压下Gd3Ga5O12晶体的弹性性能. 高压物理学报, 1987, 1(1): 32-38 . doi: 10.11858/gywlxb.1987.01.005
    [7] 马云张毅李加波蔡灵仓李剑峰柳雷敬秋民王志刚王翔翁继东 . DAC加载下材料高压声速的激光超声测量技术. 高压物理学报, 2011, 25(5): 416-420 . doi: 10.11858/gywlxb.2011.05.006
    [8] 李敏崔启良李芳菲周强武晓鑫邹广田 . 高温高压条件下液态氧的声速及弹性系数研究. 高压物理学报, 2008, 22(3): 286-290 . doi: 10.11858/gywlxb.2008.03.011
    [9] 马云张毅李加波蔡灵仓李剑峰王翔翁继东李军敬秋民李俊 . 利用激光超声技术测量介观尺度材料的声速. 高压物理学报, 2011, 25(2): 128-132 . doi: 10.11858/gywlxb.2011.02.006
    [10] 俞宇颖谭华胡建波戴诚达马云陈大年 . 利用VISAR测量LY12铝在冲击压缩下的声速. 高压物理学报, 2006, 20(2): 145-152 . doi: 10.11858/gywlxb.2006.02.006
    [11] 陈军曾代朋谭多望王荣波 . 超压爆轰产物声速以及热力学CJ点的实验研究. 高压物理学报, 2010, 24(3): 201-205 . doi: 10.11858/gywlxb.2010.03.007
    [12] 顾惠成李凤英 . 有磁和无磁状态下YIG和YAG晶体高压下的声速和衰减. 高压物理学报, 1993, 7(2): 138-142 . doi: 10.11858/gywlxb.1993.02.010
    [13] 李芳菲崔启良李敏周强邹广田 . 高温高压下液态水声速的研究不同状态方程的准确性验证. 高压物理学报, 2008, 22(3): 281-285 . doi: 10.11858/gywlxb.2008.03.010
    [14] 王贵朝 . 高压下材料弹性声速的一种经验表达式. 高压物理学报, 1988, 2(1): 92-95 . doi: 10.11858/gywlxb.1988.01.014
    [15] 王志刚徐亮李绪海王海阔贺端威孟川民 . 碳化硅-金刚石超硬复合材料的弹性性质. 高压物理学报, 2015, 29(4): 263-267. doi: 10.11858/gywlxb.2015.04.004
    [16] 马春丽武晓鑫黄凤仙李敏王晓霞周强李芳菲崔启良 . 小分子液体的高温布里渊散射研究. 高压物理学报, 2015, 29(1): 35-41. doi: 10.11858/gywlxb.2015.01.006
    [17] 胡金彪经福谦程菊鑫 . 铜的高压声速和冲击熔化. 高压物理学报, 1989, 3(3): 187-197 . doi: 10.11858/gywlxb.1989.03.003
    [18] 李西军周显明经福谦 . 一种多孔铁的高压声速和冲击熔化. 高压物理学报, 2001, 15(2): 117-121 . doi: 10.11858/gywlxb.2001.02.007
    [19] 池元斌王立中刘慎薪王秋平 . Gd2O2S:Eu荧光谱及能级的高压研究. 高压物理学报, 1994, 8(2): 87-98 . doi: 10.11858/gywlxb.1994.02.002
    [20] 杨金科龚自正邓力维张莉费英伟 . 顽火辉石(Mg0.92,Fe0.08)SiO3的冲击相变和高压状态方程及其地球物理意义. 高压物理学报, 2007, 21(1): 45-54 . doi: 10.11858/gywlxb.2007.01.008
  • 加载中
图(10)表(1)
计量
  • 文章访问数:  343
  • 阅读全文浏览量:  61
  • PDF下载量:  86
出版历程
  • 收稿日期:  2017-05-08
  • 录用日期:  2017-05-17
  • 刊出日期:  2017-12-25

高压下Fe92.5O2.2S5.3的熔化温度

    作者简介:冯磊(1987-), 男, 硕士, 主要从事固体地球物理研究.E-mail:feng_lei@whut.edu.cn
    通讯作者: 黄海军, hjhuang@whut.edu.cn
  • 武汉理工大学理学院,湖北武汉 430070
基金项目:  国家重点基础研究发展计划(973计划) 2014CB845904中央高校基本科研业务费专项资金 2015Ⅲ035中央高校基本科研业务费专项资金 2014-Ⅶ-006国家自然科学基金优秀青年科学基金 41322028

摘要: 采用反向碰撞法与光分析技术,测量了Fe92.5O2.2S5.3(质量分数比)在208 GPa下的声速,发现固态Fe92.5O2.2S5.3的纵波声速在144 GPa下开始减小,直到165 GPa完全转变为液态体波声速,表明样品的完全熔化温度为(3 500±400)K。将该熔化温度作为参考点,应用Lindeman定律并外推至地球内外核边界可知,Fe92.5O2.2S5.3的熔化温度为(5 000±400)K。通过比较Fe、Fe-O、Fe-S以及Fe-O-S的熔化温度,发现O元素对Fe熔化温度的影响很小,S元素对Fe熔化温度的降低与其含量成正比。如果外地核中S的质量分数为2%~6%,则地球内外核界面温度为5 000~5 400 K。

English Abstract

    • 众所周知,地核是由固态内核和液态外核组成,因此地球内外核界面温度等于地核物质在此压强(330 GPa)下的熔化温度,而该温度也是构建地球温度剖面的重要参考温度。根据地震波数据和高压实验结果可知,外核物质由铁和轻元素(质量分数约10%)组成[1-4],其中可能的轻元素包括氧(O)、硫(S)、硅(Si)、碳(C)、氢(H),而O和S被认为是主要的候选轻元素[5-14]。因此,为了构建地核的温度剖面,十分有必要研究O和S元素在高压下对Fe熔化温度的影响。

      对于Fe-O体系,实验研究发现:在低于50 GPa的压强范围内,共晶点温度比Fe的熔化温度低50~350 K,并且当O的含量小于共晶含量时,Fe-O体系的液相温度与Fe的熔化温度非常接近[15-19],即O对Fe熔化温度的影响很小。此外,通过实验研究还发现:Fe-S体系的共晶温度比纯Fe的熔化温度低1 000~1 500 K,并且随着压强的增大,共晶体成分向S减少的方向偏移[6-9, 20-26];随着S含量的增大,Fe-S体系的液相温度急剧降低。对于Fe-O-S三元体系,Urakawa等[11]和Tsuno等[27]使用大腔体压机分别在15 GPa和15~27 GPa压强下实验研究了Fe-O-S体系的熔化行为。Urakawa等[11]发现Fe-O-S共晶体系中存在液体不互溶区,并且随着压强的增大,这种不互溶区逐渐消失;且在15 GPa时,Fe-O-S共析线(共晶点的连线)在Fe-S轴上的投影与Fe-S在18.5 GPa下的液相线非常一致[19],说明O元素对Fe-S体系液相线温度的影响可以忽略。在更高的压强下,Huang等[14]利用冲击波压缩技术在低于230 GPa的压强范围内研究了富氧体系的Fe-O-S(Fe、O、S的质量分数比w(Fe):w(O):w(S)=90:8:2)的熔化行为,其液相线外推到内外核边界的温度为5 400 K。利用金刚石对顶砧与X射线同步辐射技术,Terasaki等[28]在0~157 GPa压强范围内测量了富硫体系的Fe-O-S(w(Fe):w(O):w(S)=85:2:13)的固相线和液相线温度,在误差范围内其液相线温度与Huang等[14]的测量结果非常接近。如果在高压下O对Fe-S体系液相线的影响非常小,那么由Huang等[14]和Terasaki等[28]的研究结果可知高压下S的含量对Fe-S体系液相线温度的影响将变得很小,与低压下的实验结论[6-9, 20-26]不同。而此结论是否正确,还需要进行进一步的高压实验。为此,本研究选择富硫体系的Fe-O-S样品作为研究对象,利用冲击波压缩技术研究其在高压下的熔化行为。

    • 初始原料为高纯Fe、FeS和Fe2O3粉末,按质量比78.15:14.53:7.32在玛瑙研钵中研磨1 h。将干燥后的混合物压制成片,并密封于钼“杯”中,以杜绝Fe在高温下被氧化的可能性。将钼“杯”放置在大腔体压机的样品腔中进行烧结,压强约为2 GPa,温度为850 ℃。X射线衍射(XRD)分析发现,样品中不存在Fe2O3相,表明Fe2O3已经完全还原成FeO。根据混合物中Fe、FeS和Fe2O3的质量比,计算得到样品中Fe、O、S的质量分数比w(Fe):w(O):w(S)=92.5:2.2:5.3,以下将样品表示为Fe92.5O2.2S5.3。样品的平均密度为6.88(0.02)g/cm3,与理论密度6.89 g/cm3非常接近。

      为了判断样品在冲击压缩过程中是否完全熔融,即在冲击加载下纵波声速是否变为体波声速,必须测量样品在冲击加载过程中的声速。尽管Huang等[13]已经报道了Fe92.5O2.2S5.3在高压下的声速数据,但是受篇幅限制,对测量方法和数据处理过程只进行了简单介绍,并未详细说明,为此本研究将对此进行详细论述。在测量Fe92.5O2.2S5.3的声速时,本研究采用了反向碰撞技术[29-30]和光分析技术[14]。反向碰撞技术可以在较低的压强下同时测出纵波声速和体波声速,而光分析技术可在较高的压强下测量样品熔融前的纵波声速和熔融后的体波速度。图 1(a)为反向碰撞技术的原理图。Fe92.5O2.2S5.3样品作为飞片,直接高速碰撞镀有铝膜的LiF单晶窗口。为了避免前冲气体对测量信号产生影响,在铝膜表面用环氧树脂粘贴一层8 μm厚的铝箔。碰撞后,在样品和铝箔中分别产生冲击波DsL和DAlL(波速分别为DsLDAlL),如图 1(b)所示。当铝箔中的冲击波DAlLt0时刻到达铝膜/LiF单晶界面时,全光纤位移干涉仪(Displacement Interferometer System for Any Reflector,DISAR)[31]会观测到其粒子速度突然增加,并保持不变。样品中的冲击波DsL传播到样品后表面时,会反射稀疏波CsL。当铝箔中的稀疏波CAlLt1时刻到达铝膜/LiF窗口界面时,粒子速度会缓慢减小。上述粒子速度的变化过程如图 2所示。

      图  1  反向碰撞法实验装置图(a)和反向碰撞法波系作用图(b)

      Figure 1.  Sample configuration (a) and wave interaction (b) for sound velocity measurements in the reverse impact method

      图  2  实验No.081225和No.081226的界面粒子速度历史

      Figure 2.  Particle velocity histories at the impact interface for experiments No.081225 and No.081226

      在拉格朗日坐标系中,样品中的纵波声速CL采用如下公式计算

      $ {\mathit{C}_{\rm{L}}}{\rm{ = }}\frac{{{\mathit{d}_{\rm{s}}}}}{{{\mathit{t}_{\rm{1}}}{\rm{ - }}{\mathit{t}_{\rm{0}}}{\rm{ - }}{\mathit{d}_{\rm{s}}}{\rm{/}}\mathit{D}_{\rm{s}}^{\rm{L}}}}\frac{{{\mathit{\rho }_{\rm{0}}}}}{\mathit{\rho }} $

      式中:ds为样品厚度,ρ0ρ分别为样品的初始密度和压缩后密度。在卸载过程中,如果可以确定弹性-塑性转变点在t2时刻,则体波声速CB可由下式确定

      $ {\mathit{C}_{\rm{B}}}{\rm{ = }}\frac{{{\mathit{d}_{\rm{s}}}}}{{{\mathit{t}_{\rm{2}}}{\rm{ - }}{\mathit{t}_{\rm{0}}}{\rm{ - }}{\mathit{d}_{\rm{s}}}{\rm{/}}\mathit{D}_{\rm{s}}^{\rm{L}}}}\frac{{{\mathit{\rho }_{\rm{0}}}}}{\mathit{\rho }} $

      然而,根据粒子速度历史(见图 2),很难准确判定t2的值。Asay等[32]发现样品在弹塑性卸载过程中,拉格朗日声速与粒子速度呈线性变化,他们将这种线性关系外推到卸载前的粒子速度,由此得到样品在卸载前的体波声速(见图 3)。

      图  3  拉格朗日声速与粒子速度的关系

      Figure 3.  Lagrangian sound velocity vs.particle velocity

      光分析技术实验装置如图 4所示。样品靶由3个不同厚度的样品构成,样品后表面为LiF单晶窗口,为了避免样品与窗口间的空隙在冲击过程中发出强光对测量信号产生干扰,在LiF表面粘贴一层金属钽箔。光分析技术中波系的相互作用见图 5。当飞片与样品高速碰撞时,在飞片与样品中分别产生两个传播方向相反的冲击波(DfL和DsL,波速分别为DfLDsL),当样品中的冲击波DsLB1时刻到达样品/LiF窗口界面时,DISAR会监测到粒子速度突然增加,同时在界面向样品内部反射一个稀疏波Cs′L。飞片中的冲击波DfL传播到飞片后表面时,会反射一个向前传播的稀疏波,它们经过飞片和样品时的速度分别为CfLCsL。当稀疏波CsL与来自样品后表面的稀疏波Cs′L相遇时,其稀疏波速度变为Cs1L,并在E1时刻到达样品/LiF窗口界面,此时会观测到界面粒子速度缓慢减小。在其他样品/窗口界面也会观测到类似的现象。由于稀疏波在样品中传播的速度CsL大于冲击波速度DsL,如果样品足够厚,稀疏波总会在样品中追赶上冲击波,如图 5A点所示(所对应的厚度d为追赶厚度)。由于△A1B1E1∽△A2B2E2,因此Δt1t2=(d-d1)/(d-d2),其中d1d2分别为样品1和样品2的厚度。如果在实验中可以布置更多厚度不同的样品,可以发现Δti与厚度di呈线性变化,当Δt=0即可得到追赶厚度。图 6显示了实验No.110505的粒子速度历史,样品厚度分别为1.535、2.336和3.065 mm,测量的时间间隔分别为579、536和500 ns,通过线性拟合得出样品的追赶厚度d=(12.74±0.06)mm,如图 7所示。

      图  4  采用光分析技术测量声速的实验装置

      Figure 4.  Experimental configuration for sound velocity measurements using optical analyzer technique

      图  5  光分析技术测量声速实验中波的传播

      Figure 5.  Wave interaction for sound velocity measurements in optical analyzer technique

      图  6  实验No.110505在198 GPa下的粒子速度历史

      Figure 6.  Particle velocity history for experiment No.110505 at 198 GPa

      图  7  实验No.110505在198 GPa下的时间差Δt与样品厚度的关系

      Figure 7.  Time intervals Δt vs.sample thickness for experiment No.110505 at 198 GPa

      在拉格朗日坐标下声速为

      $ \frac{\mathit{d}}{{\mathit{D}_{\rm{s}}^{\rm{L}}}}{\rm{ = }}\frac{{{\mathit{d}_{\rm{f}}}}}{{\mathit{D}_{\rm{f}}^{\rm{L}}}}{\rm{ + }}\frac{{{\mathit{d}_{\rm{f}}}}}{{\mathit{C}_{\rm{f}}^{\rm{L}}}}{\rm{ + }}\frac{\mathit{d}}{{\mathit{C}_{\rm{s}}^{\rm{L}}}} $

      式中:df为飞片厚度。在本实验过程中,为了在更高压力下获得声速,飞片材料采用Ta。在拉格朗日坐标系中,飞片和样品中的冲击波速度可采用以下公式计算

      $ \mathit{D}_{\rm{f}}^{\rm{L}}{\rm{ = }}{\mathit{C}_{{\rm{0, Ta}}}}{\rm{ + }}{\mathit{\lambda }_{{\rm{Ta}}}}\left( {\mathit{W}{\rm{ - }}\mathit{u}} \right) $

      $ \mathit{D}_{\rm{s}}^{\rm{L}}{\rm{ = }}{\mathit{C}_{{\rm{0, s}}}}{\rm{ + }}{\mathit{\lambda }_{\rm{s}}}\mathit{u} $

      式中:C0λ为Hugoniot参数,W为飞片速度,u为粒子速度。对于Fe92.5O2.2S5.3样品,C0,s=(3.71±0.12)km/s,λs=1.61±0.04[13];对于Ta飞片,C0,Ta=3.293 km/s,λTa=1.307[33]。粒子速度u由冲击阻抗匹配法得到,Ta飞片的CfL取Brown等[34]给出的欧拉声速Cf,样品的欧拉声速为CsLρ0/ρ的乘积。

    • Fe92.5O2.2S5.3的纵波声速和体波声速的测量结果如图 8所示。当冲击压强为142 GPa时,反向碰撞法测得的样品卸载信号不理想,根据实验数据计算得到的体波声速有很大误差,因此图 8中没有给出142 GPa下的体波声速。在144~165 GPa压强范围内,声速测量值沿Hugoniot线出现不连续性;而当压强高于165 GPa时,体波声速测量值与理论计算值[13]一致,表明样品在165 GPa下已经完全熔化,转变为液体。根据能量守恒关系,Fe92.5O2.2S5.3的液相线温度由(6)式和(7)式得到[35]

      $ \begin{array}{l} {\mathit{E}_{\rm{H}}}{\rm{ = }}\frac{{\rm{1}}}{{\rm{2}}}\left( {{\mathit{V}_{\rm{0}}}{\rm{ - }}\mathit{V}} \right){\mathit{p}_{\rm{H}}}{\rm{ + }}{\mathit{E}_{\rm{0}}}\\ \;\;\;\;\;\;{\rm{ = }}{\mathit{E}_\mathit{x}}{\rm{ + }}\int_0^{{\mathit{T}_{\rm{m}}}} {{\mathit{c}_\mathit{V}}{\rm{d}}\mathit{T}{\rm{ + }}{\mathit{T}_{\rm{m}}}{\rm{\Delta }}\mathit{S}} \end{array} $

      $ {\mathit{p}_{\rm{H}}}{\rm{ = }}\frac{{{\mathit{\rho }_{\rm{0}}}\mathit{C}_0^2{\rm{(1 - }}{\mathit{\rho }_{\rm{0}}}{\rm{/}}\mathit{\rho }{\rm{)}}}}{{{{{\rm{[1 - }}\mathit{\lambda }{\rm{(1 - }}{\mathit{\rho }_{\rm{0}}}{\rm{/}}\mathit{\rho }{\rm{)]}}}^{\rm{2}}}}} $

      图  8  Fe92.5O2.2S5.3的声速与冲击压强的关系

      Figure 8.  Sound velocity of Fe92.5O2.2S5.3 vs.shock pressure

      式中:E为比内能,V为比容,p为压强,T为温度,Tm为熔化温度;下标“H”代表沿Hugoniot线,“x”代表沿0 K等温线,“0”代表室温;ΔS为熔融熵,ΔS=0.19R[36],其中R为气体常数;cV为定容比热容,包括晶格比热容(cVl)和电子比热容(cVe)。Fe92.5O2.2S5.3cV可根据Fe、FeO、FeS的定容比热容,通过可加性原理得到,即

      $ {\mathit{c}_\mathit{V}}{\rm{ = }}\sum\limits_\mathit{i} {{\mathit{w}_\mathit{i}}{\mathit{c}_{\mathit{V}{\rm{l}}\mathit{i}}}} {\rm{ + }}\sum\limits_\mathit{i} {{\mathit{w}_\mathit{i}}{\mathit{\beta }_{{\rm{0}}\mathit{i}}}{{{\rm{(}}{\mathit{\rho }_{{\rm{0}}\mathit{i}}}{\rm{/}}{\mathit{\rho }_\mathit{i}}{\rm{)}}}^{{\mathit{\kappa }_\mathit{i}}}}\mathit{T}} $

      式中:wi为组分i的质量分数;β0iκi为组分i的电子热参数;ρ0iρi分别为组分i的初始密度和密度;cVli为组分i的晶格比热容贡献,在高温下满足Dulong-Petit定律,近似等于3R/μ(μ为摩尔质量)。组分Fe、FeO、FeS的详细参数如表 1所示。在(6)式中,固态下冷能Ex的计算方程为

      $ {\mathit{E}_\mathit{x}}{\rm{ = }}{\mathit{E}_{\rm{H}}}{\rm{ - }}\int_0^\mathit{T} {{\mathit{c}_\mathit{V}}{\rm{d}}\mathit{T}} $

      Material ρ0/(g/cm3) C0/(km/s) λ γ0 q β0/(J·kg-1·K-2) κ E0/(kJ)
      α-Fe 7.85[43] 3.935[43] 1.578[43]
      ε-Fe 8.298[44] 4.720* 1.523* 1.76* 0.76* 0.091[45] 1.34[45] 76.268 4[44]
      FeO(B1) 5.71[46] 5.83* 0.99* 57.796 7[47]
      FeO(B8) 6.05[46, 48] 4.486* 1.699* 1.8[48] 1[48] 0[49] 0[49]
      FeS 4.602[50] 2.947[50] 1.578[50] 1.54[50] 1[50] 0.25[51] 1.34[51] 33.707 5[52]
      Fe90O8S2 6.69[14] 3.97[14] 1.58[14] 1.85[14] 0.87[14] 0.075[14] 1.39[14] 67.315
      Fe92.5O2.2S5.3 6.88 3.71 1.61 1.82 0.85 0.119 1.407 68.258
      Note:The asterisk * represents the fitted results from the Hugoniot data.

      表 1  Fe、FeO、FeS和Fe-O-S体系参数

      Table 1.  Parameters for Fe, FeO, FeS and Fe-O-S system

      固态样品的冲击温度可由下式给出[37]

      $ {\rm{d}}\mathit{T}{\rm{ = - }}\mathit{T}\left( {\frac{{{\mathit{\gamma }_{{\rm{eff}}}}}}{\mathit{V}}} \right){\rm{d}}\mathit{V}{\rm{ + }}\frac{{\rm{1}}}{{{\rm{2}}{\mathit{c}_\mathit{V}}}}{\rm{[(}}{\mathit{V}_{\rm{0}}}{\rm{ - }}\mathit{V}{\rm{)d}}\mathit{p}{\rm{ + (}}\mathit{p}{\rm{ - }}{\mathit{p}_{\rm{0}}}{\rm{)d}}\mathit{V}{\rm{]}} $

      式中:γeff为Grüneisen参数,包括晶格的贡献γl和电子的贡献γeγeff的计算公式为

      $ {\mathit{\gamma }_{{\rm{eff}}}}{\rm{ = (}}{\mathit{c}_{\mathit{V}{\rm{l}}}}{\mathit{\gamma }_{\rm{l}}}{\rm{ + }}{\mathit{c}_{\mathit{V}{\rm{e}}}}{\mathit{\gamma }_{\rm{e}}}{\rm{)/}}{\mathit{c}_\mathit{V}} $

      Fe92.5O2.2S5.3的晶格参数γl=γ0(ρ0/ρ)q,其中γ0为常温下的Grüneisen参数,q为参数;由可加性原理及表 1中Fe、FeO、FeS的晶格参数可计算得到γl。Fe92.5O2.2S5.3的电子参数γe与纯Fe的γe较为接近,γe=2[38]。将以上参数代入(6)式中,计算得出Fe92.5O2.2S5.3样品在165 GPa压强下的熔化温度为(3 500±400)K。基于熔化温度的计算值,应用Lindeman原理,得出Fe92.5O2.2S5.3样品的熔化曲线,如图 9所示。

      $ {\rm{\Delta }}{\mathit{T}_{\rm{m}}}{\rm{ = 2}}{\mathit{T}_{\rm{m}}}\left( {{\mathit{\gamma }_{{\rm{eff}}}}{\rm{ - }}\frac{{\rm{1}}}{{\rm{3}}}} \right)\frac{{{\rm{\Delta }}\mathit{\rho }}}{\mathit{\rho }} $

      图  9  Fe、Fe-O、Fe-S、Fe-O-S体系在高压下的熔化温度和液相线温度随压强的变化(红色实线为本研究得到的Fe92.5O2.2S5.3的熔化温度。Fe:Ma等[40],标记为M04;Anzellini等[39],标记为A13。Fe-O-S体系:Terasaki等[28],标记为T11;Huang等[14],标记为H10。Fe-S体系:Kamada等[26],标记为K12。)

      Figure 9.  Melting temperature of Fe-O-S system compared with those of Fe, Fe-O and Fe-S (The red solid line shows the melting temperature of Fe92.5O2.2S5.3 from this study.The lines labeled as M04 and A13 represent the melting temperatures of Fe from Ref.[40] and Ref.[39] respectively.The dashed lines labeled as H10 and T11 represent the melting temperature of Fe90O8S2[14] and the liquidus and solidus temperature of Fe-O-S[28].K12 represent the melting relationships in the Fe-Fe3S system up to the outer core conditions[26].)

      为了研究O元素与S元素对纯Fe熔化温度的影响,将Fe92.5O2.2S5.3的熔化温度与已发表的Fe、Fe-O-S、Fe-O、Fe-S体系的熔化温度和液相线温度进行比较,如图 9所示。对于Fe,最近Anzellini等[39]通过快速XRD研究了Fe在高压下的熔化温度,在实验误差范围内比Ma等[40]和Shen等[41]利用XRD得到的测量结果以及Huang等[42]利用动高压实验得到的计算结果高300 K。对于Fe-O体系,Seagle等[18]测量的Fe-FeO液相线温度非常接近,说明高压下O对Fe熔化温度的影响很小。对于Fe-O-S体系,由冲击波实验测得的Fe90O8S2[14]和Fe92.5O2.2S5.3的液相线温度均低于Fe的熔化温度,并且随着S含量的增大,Fe-O-S体系的液相线温度越来越低,但高于Fe-FeS的熔化温度[21-26]。在图 9中,对于Fe85O2S13,Terasaki等[28]用激光加热金刚石压砧和原位XRD技术测得的液相线温度比较高,由于其S含量超过共晶点S的含量,因此所测结果有可能是Fe3S在高压下的液相线温度。以上仅定性分析了O和S在高压下对Fe熔化温度的影响。为了进一步说明O和S对Fe熔化温度的影响,将Fe-O-S的液相线温度投影到Fe-S相图上,如图 10所示。Fe-O-S在15 GPa下的液相线[11]与Fe-S在18.5 GPa下的液相线[19]非常一致,说明O元素对Fe-S液相线的影响很小。Kamada等[26]测量了123 GPa下Fe-Fe3S的共晶点温度和液相温度,发现Fe90O8S2、Fe92.5O2.2S5.3的液相线温度恰好在Fe的熔化温度与Fe-Fe3S的共晶点温度连线之间,并低于Fe-Fe3S在液相区的温度。在330 GPa下,Fe90O8S2和Fe92.5O2.2S5.3的液相温度低于Fe的熔化温度,并高于Kamada等[26]依据实验数据估算的Fe92S8(质量分数比)的熔化温度。综合图 10所示相图可知,Fe-O-S与Fe-S的液相线在高压下基本一致;O对Fe熔化温度的影响很小,而S对Fe熔化温度的降低与其含量成正比。如果外地核是Fe-O-S体系,根据地球化学估算,地核中S的质量分数在2%~6%之间[53-54],则内外核边界温度为5 000~5 400 K,与根据核幔边界所推算的温度5 200~5 700 K[55]在误差范围内非常接近。

      图  10  高压下Fe-S与Fe-O-S体系液相线温度的比较

      Figure 10.  Comparison of liquid-phase temperature between Fe-S and Fe-O-S system under high pressure

    • 测量了Fe92.5O2.2S5.3在208 GPa压强下的声速,发现其纵波声速在144 GPa下开始逐渐减小,到165 GPa时完全转变为液态体波声速。通过计算得到样品完全熔化温度为(3 500±400)K。将该熔化温度作为参考点,应用Lindeman定律并外推至地球内外核边界可知,Fe92.5O2.2S5.3的熔化温度为(5 000±400)K。通过比较Fe、Fe-O、Fe-S以及Fe-O-S的熔化温度发现,Fe-O-S体系与Fe-S体系的液相线温度非常接近,说明O元素对Fe熔化温度的影响很小,而S元素对Fe熔化温度的降低与其含量成正比。如果外地核中S元素的质量分数在2%~6%之间,则内外核界面温度为5 000~5 400 K。

参考文献 (55)

目录

    /

    返回文章
    返回