天然橄榄石单晶的压缩性

秦霏; 王英; 巫翔; 秦善; 李晖; 李晓东; 杨科

doi:10.11858/gywlxb.2016.01.003

天然橄榄石单晶的压缩性

doi: 10.11858/gywlxb.2016.01.003

秦霏^1,,
王英¹,
巫翔^1,*, ,,
秦善¹,
李晖²,
李晓东³,
杨科⁴

1.
北京大学地球与空间科学学院，造山带与地壳演化教育部重点实验室，北京 100871
2.
北京工业大学固体微结构与性能研究所，北京 100124
3.
中国科学院高能物理研究所同步辐射实验室，北京 100049
4.
中国科学院上海应用物理研究所，上海 201024

基金项目:

国家自然科学基金 U1232204

国家自然科学基金 11075175

详细信息

作者简介:
秦霏(1991-), 女, 博士研究生, 主要从事地球高温高压晶体化学研究.E-mail:fei.qin@pku.edu.cn

通讯作者:
巫翔(1978-), 男, 博士, 研究员, 主要从事实验地球化学及矿物物理研究.E-mail:xiang.wu@pku.edu.cn

中图分类号: O521.2
计量
- 文章访问数: 7938
- HTML全文浏览量: 2539
- PDF下载量: 225
出版历程
- 收稿日期: 2015-02-01
- 修回日期: 2015-04-17

Compressibility of Natural Olivine Single-Crystals

QIN Fei^1
,,
WANG Ying¹,
WU Xiang^{1,*
, ,},
QIN Shan¹,
LI Hui²,
LI Xiao-Dong³,
YANG Ke⁴

1.
Key Laboratory of Orogenic Belts and Crustal Evolution, School of Earthand Space Sciences, Peking University, Beijing 100871, China
2.
Institute of Microstructure and Property of Advanced Materials, Beijing University of Technology, Beijing 100124, China
3.
Synchrotron Radiation Laboratory, Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 100049, China
4.
Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Shanghai 201024, China

摘要

摘要: 采用金刚石对顶砧高压装置和同步辐射X射线衍射技术，在常温下对两种天然橄榄石单晶(含铁的铁镁橄榄石(Mg_1.83Fe_0.17)SiO₄和不含铁的纯镁橄榄石Mg₂SiO₄)进行了高压结构行为研究。实验结果表明，两种天然橄榄石单晶在最高压力为9.9GPa的条件下没有发生相变。采用二阶Birch-Murnaghan状态方程描述压力-晶胞体积(p-V)关系，获得了铁镁橄榄石和纯镁橄榄石在零压下的晶胞体积(0.2929(3)、0.2899(1)nm³)和体积模量(140(3)、151(2)GPa)。两种镁橄榄石呈高度各向异性，轴压缩性表现为：b轴的压缩率最大，c轴次之，a轴的压缩率最小。随着铁组分的加入，橄榄石的体波速减小，两种橄榄石的体波速相差5.2%，该结果对于研究上地幔的物质构成和弹性波速具有一定的意义。
- 橄榄石 /
- 压缩性 /
- 各向异性 /
- 体波速
Abstract: In situ X-ray diffraction measurements for natural Fe-bearing and Fe-free Mg₂SiO₄ forsterite single-crystals ((Mg_1.83Fe_0.17)SiO₄ and Mg₂SiO₄) were performed with diamond anvil cell high-pressure device and synchrotron radiation X-ray diffraction at room temperature.The experimental results showed that no phase transition of two natural olivine single-crystals was observed within the upper highest pressure of 9.9GPa.The isothermal relationship of pressure and volume in the two kinds of crystals is well described by the second-order Birch-Murnaghan equation of state with the zero-pressure unit cell volume and zero-pressure bulk modulus of 0.2929(3)nm³ and 140(3)GPa for Fe-bearing forsterite and 0.2899(1)nm³ and 151(2)GPa for Fe-free forsterite, respectively.The linear compressibilities along a, b and c axes (k_a, k_b and k_c) are elastically anisotropic in the two kinds of forsterite and the anisotropy can be described as k_b>k_c>k_a.With the addition of iron into the crystals, the bulk velocity of olivine decreases and the difference between the two kinds of forsterite single-crystals is 5.2%.This result will be significant for the study of the composition of the earth's upper mantle and its elastic wave velocities.
- olivine /
- compressibility /
- anisotropy /
- bulk velocity

HTML全文

1. 引言

橄榄石(Mg, Fe)₂SiO₄作为上地幔最主要的代表性矿物之一，也是地幔橄榄岩和辉石岩的主要组成矿物，约占上地幔的60%~65%^[1]，其高压多形为β相(瓦兹利石)和γ相(林伍德石)。一般的橄榄石由镁橄榄石(Mg₂SiO₄)和铁橄榄石(Fe₂SiO₄)两种端元组分以不同的比例组成。橄榄石的高压相变和分解被认为是导致地幔转换带中地震波不连续面形成的直接原因。现有的研究普遍认为：橄榄石的α相→β相相变导致了410km地震波速间断面的形成，β相→γ相相变是520km波速间断面形成的原因^[1-2]，而660km处橄榄石的后尖晶石相变(即γ相橄榄石分解成钙钛矿和方镁铁矿)则标志着下地幔的开始^[3]。因此，研究橄榄石在高压条件下的结构性质，有助于了解地幔不连续面的成因，对于探索整个地幔的物质组成和演化、俯冲板片深源地震等地球深部问题具有重要的意义。

橄榄石属于岛状硅酸盐矿物，常压下为斜方晶系，空间群为Pbnm；其结构中有2种不同的八面体位置和3种不同的O原子位置，O做近似六方最紧密堆积，Mg和Fe填充于八面体空隙中，每个O原子与3个八面体阳离子和1个四面体阳离子相连接。国内外学者通过实验和理论计算等方法研究发现，镁铁比例、含水量等因素都会对橄榄石的各向异性产生显著的影响^[4-9]。Suzuki^[10]研究了具有不同端元组分的橄榄石的线膨胀系数，发现纯镁橄榄石的膨胀系数最大，并且体积热膨胀与晶胞体积呈反比关系。Hazen^[11-12]研究了镁、铁橄榄石的压缩特性，结果表明：铁橄榄石的体积模量(113GPa)小于镁橄榄石(132GPa)，更易于压缩。很多学者还对橄榄石在高温高压下的弹性性质进行了研究，讨论并分析了Fe和OH组分对橄榄石弹性波速的影响。Zha等人^[8]利用布里渊散射技术研究了镁橄榄石在16GPa压力范围内的弹性性质，得到了体积模量(K，GPa)和剪切模量(G，GPa)与压力(p，GPa)的关系：K=128.8(5)+4.2(2)p，G=81.6(2)+1.4(1)p。Speziale等人^[13]的研究表明，在12.1GPa的压力范围内铁橄榄石的零压体积模量(K₀)和零压剪切模量(G₀)分别为136.3(2)和51.2(2)GPa，∂K/∂p=4.9(1)。Liu等人^[14]通过第一性原理计算，研究了不同的镁铁比例对橄榄石弹性波速的影响，发现Fe含量的增加会使体积模量增大，使剪切模量减小。Ono等人^[15]给出了铁橄榄石到尖晶石相的相变边界，其压力-温度(p-T)关系为p(GPa)=0.5+0.0034T(K)。

尽管前人对橄榄石的各向异性和弹性性质进行了诸多探讨，但是对于具有不同含铁量的橄榄石的压缩性仍存在一些争议。为此，本工作拟采用两种天然的含铁和不含铁的镁橄榄石，结合金刚石对顶砧(Diamond Anvil Cell，DAC)装置和同步辐射X射线衍射(X-Ray Diffraction, XRD)技术，对两种镁橄榄石单晶在高压下的压缩性、各向异性进行对比研究，进一步探讨不同的橄榄石铁组分对弹性波速的影响。

2. 样品与实验

实验所用样品为天然的宝石级单晶橄榄石，其中黄绿色样品产自河北省张家口市大麻坪汉诺坝，无色样品的产地不详。利用JXA-8100型电子探针显微分析仪，测得两种样品的化学式分别为(Mg_1.83Fe_0.17)SiO₄和Mg₂SiO₄，即铁镁橄榄石和纯镁橄榄石。Mo靶的单晶X射线衍射数据表明：两种样品的空间群均为Pbnm；铁镁橄榄石样品的晶胞参数a=0.47956(2)nm，b=1.0178(1)nm，c=0.60369(5)nm，纯镁橄榄石样品的晶胞参数a=0.47382(6)nm，b=1.0218(2)nm，c=0.59817(6)nm。采用3M公司生产的金刚石薄膜砂纸，对两种单晶样品进行双面抛光，直至厚度约为15μm。选取尺寸约为40μm×50μm×15μm的样品作为压缩实验的测试对象。

高压装置采用Symmetry型DAC，砧面直径为300μm；封垫材料为T301不锈钢和Re，封垫厚度分别为200和250μm，预压后封垫厚度约为40μm；样品腔的直径约为150μm。将预处理好的单晶样品与压标物质Au箔片^[16]一起放入样品腔中，传压介质采用体积比为4:1的甲醇-乙醇混合溶液。铁镁橄榄石单晶的原位高压X射线衍射实验在北京同步辐射装置(BSRF)高压站(4W2束线)完成，X射线的波长为0.06199nm，聚焦光斑尺寸为30μm×40μm，每条谱的采谱时间约为300s。纯镁橄榄石单晶的原位高压X射线衍射实验在上海同步辐射装置(SSRF)高压站(15U束线)完成，X射线的波长为0.06199nm，聚焦光斑尺寸为4μm×3μm，每条谱的采谱时间为15s。两组实验的最高实验压力为9.9GPa。在测试过程中，为了得到更多的单晶衍射点，连续旋转DAC，使DAC与入射X射线的夹角在15°~20°区间内连续变化。采用GSE_ADA/RSV软件提取单晶衍射图谱中各衍射点对应的面网间距^[17]，进行晶面指标化，最后通过UnitCell软件拟合得到样品的晶胞参数^[18]。橄榄石样品的p-V关系通过二阶Birch-Murnaghan状态方程(Equation of State, EOS)进行描述^[19]

$p\left( V \right) = \frac{3}{2}{K_0}\left[ {{{\left( {\frac{{{V_0}}}{V}} \right)}^{7/3}} - {{\left( {\frac{{{V_0}}}{V}} \right)}^{5/3}}} \right]$

(1)

式中：p为压力，V为晶胞体积，V₀和K₀分别为零压下的晶胞体积和体积模量。

3. 实验结果分析和讨论

在本实验所测得的每张衍射图谱中至少可观测到25个单晶衍射点，除了来自金刚石压砧的衍射点外，其他所有的衍射点都可以用橄榄石晶体结构模型进行指标化。图 1为纯镁橄榄石和铁镁橄榄石的单晶衍射图谱。随着压力增加至9.9GPa，两种样品的单晶图谱中均没有出现新的衍射点，且衍射点未出现拉伸、形变、劈裂等现象，说明在此静水压力范围内两种镁橄榄石的结构保持稳定，没有发生相变。在数据处理过程中，最少使用了12个单晶衍射点进行晶胞参数拟合。例如：对于铁镁橄榄石，选取了(112)、(113)、(132)、(133)、(134)、(123)、(232)、(024)、(240)、(111)、(121)和(122)面网指数进行晶胞参数拟合(见表 1)。两种橄榄石单晶的晶胞体积随压力变化的关系如图 2所示，可见，随着压力的增加，晶胞体积逐渐被压缩。采用二阶Birch-Murnaghan状态方程对两种镁橄榄石单晶的实验结果进行p-V状态方程拟合，得到两种单晶样品的状态方程参数：对于铁镁橄榄石，V₀=0.2929(3)nm³，K₀=140(3)GPa；对于纯镁橄榄石，V₀=0.2899(1)nm³，K₀=151(2)GPa。

图 1 纯镁橄榄石和铁镁橄榄石单晶的XRD谱(方框所示为实验所得的单晶衍射点，数字表示对应的指标化结果；“D”表示金刚石的衍射点，衍射环为压标Au的衍射环)

Figure 1. Single-crystal XRD patterns of Fe-free forsterite and Fe-bearing forsterite (Squares show single-crystal diffraction peaksand figures indicate indexed results.'D' represents diamond reflections, while the diffraction rings in the patterns are from Au)

下载: 全尺寸图片幻灯片

表 1 不同压力下铁镁橄榄石和纯镁橄榄石的晶胞参数

Table 1. Unit-cell parameters of Fe-bearing and Fe-free forsterite at various pressures

Material	p/(GPa)	a/(nm)	b/(nm)	c/(nm)	V/(nm³)
Fe-free forsterite	0.0001	0.47382(6)	1.0218(2)	0.59817(5)	0.28961(6)
	0.5	0.47283(4)	1.0236(1)	0.59743(4)	0.28917(3)
	1.2	0.47279(4)	1.0206(2)	0.59637(4)	0.28776(6)
	2.9	0.47124(4)	1.0168(2)	0.59438(4)	0.28481(6)
	4.0	0.47015(4)	1.0149(2)	0.59275(4)	0.28284(6)
	5.3	0.46926(4)	1.0126(2)	0.59106(4)	0.28060(6)
	6.2	0.46872(4)	1.0089(2)	0.58968(4)	0.27885(6)
	7.0	0.46810(4)	1.0085(2)	0.58839(4)	0.27777(6)
	7.7	0.46761(4)	1.0077(2)	0.58746(4)	0.27667(6)
	8.7	0.46677(4)	1.0043(2)	0.58692(4)	0.27514(6)
	9.9	0.46609(4)	1.0013(2)	0.58564(4)	0.27240(5)
Fe-bearing forsterite	0.0001	0.47832(2)	1.0150(1)	0.60201(5)	0.29230(4)
	2.6	0.47847(7)	1.0114(1)	0.59514(3)	0.28801(3)
	4.1	0.47840(7)	1.0070(2)	0.59284(2)	0.28561(3)
	4.8	0.47682(7)	1.0058(2)	0.59165(3)	0.28376(3)
	5.8	0.47538(8)	1.0047(2)	0.59022(2)	0.28200(3)
	7.0	0.47413(8)	1.0023(2)	0.58876(3)	0.27979(4)
	7.6	0.47407(8)	1.0000(2)	0.58805(3)	0.27880(3)
	8.4	0.47306(8)	0.9996(2)	0.58687(2)	0.27753(3)
	9.1	0.47199(9)	0.9984(3)	0.58630(2)	0.27628(4)
	9.7	0.47122(8)	0.9982(3)	0.58553(4)	0.27542(3)

下载: 导出CSV

| 显示表格

图 2 铁镁橄榄石和纯镁橄榄石在常温下的p-V关系

Figure 2. p-V relationship for Fe-bearing forsteriteand Fe-free forsterite at 300K

下载: 全尺寸图片幻灯片

图 3 晶胞参数(a、b和c)随压力的变化

Figure 3. Lattice parameters a, b and c as a function of pressure

下载: 全尺寸图片幻灯片

图 3显示了两种橄榄石单晶的晶格常数(a、b、c)随压力(p)变化的关系。通过对实测晶轴大小(d)进行线性拟合，得到轴压缩率K_d，其拟合公式为

$d = {d_0} + {K_{\rm{d}}}\mathit{p}$

(2)

式中：d₀为零压下的晶轴大小。从图 3可以清楚地看到，随着压力的增加，两种镁橄榄石样品的各个轴长均基本呈线性减小。由此可以得到：对于铁镁橄榄石，a、b、c轴的轴压缩率分别为(-1.077、-1.848、-1.338)×10^-21m/Pa；对于纯镁橄榄石，a、b、c轴的轴压缩率分别为(-0.749、-2.248、-1.296)×10^-21m/Pa。两种镁橄榄石存在相同的压缩规律：b轴的压缩率最大，c轴次之，a轴的压缩率最小，即a轴最难压缩。

本研究获得的两种天然镁橄榄石的零压体积模量K₀为140GPa(铁镁橄榄石)和151GPa(纯镁橄榄石)。零压体积模量K₀的实测值不仅与实验细节和方案(样品成分、单晶/粉晶、压力标定、传压介质、表征方法等)相关，还与状态方程的拟合过程有关，即K₀与K′₀(K′₀=-(∂K/∂p)|_p=0)呈负相关性。为了更合理地比较K₀，采用统一的处理模式，即在固定K′₀=4的条件下，对前人报道的p-V数据进行重新拟合，获得K₀值，结果列于表 2，其中p_max为最高压力。从表 2可以看出，本研究得到的零压体积模量略高于前人报道的研究结果。其主要原因可能是实验所采用的压标物质不同：本实验采用内标法，通过Au的状态方程确定压力；Downs^[20]和Andrault^[21]等人通过测定红宝石的荧光峰来定压；而Nestola等人^{[4, 22]}则采用石英来定压。此外，在传压介质的使用上也不尽相同：本实验采用体积比为4:1的甲醇-乙醇混合溶液，Nestola等人^{[4, 22]}采用体积比为16:3:1的甲醇-乙醇-水混合溶液，而Andrault等人^[21]则未使用传压介质。在理论计算方面，Liu等人^[14]通过第一性原理计算得到Fe组分的增加会使镁橄榄石体积模量增大的结论，但是他们没有考虑Fe的3d电子的强关联作用，理论计算时未加入On-Site Coulomb作用能(U)，计算结构模型中也没有考虑Fe/Mg有序占位等问题，这些因素都会对最终结果产生一定的影响。橄榄石的晶体结构表现为：孤立的[SiO₄]四面体由金属阳离子Mg²⁺和Fe²⁺连接起来。结构中，O原子做近似六方最紧密堆积，Mg和Fe填充1/2的八面体空隙，Si填充1/8的四面体空隙。对比两种橄榄石的压缩性不难发现，随着含Fe量的增加，镁橄榄石的体积模量出现减小趋势。由于Fe原子的半径大于Mg原子，并且Fe—O键长大于Mg—O键长，在镁橄榄石八面体位置上Fe对Mg的置换，使得铁镁橄榄石更容易压缩，因此其弹性模量要小于纯镁橄榄石。在镁橄榄石的结构中，[(Mg, Fe)O₆]八面体沿c轴共棱连接呈锯齿状链，沿b轴方向的孔隙度更大些，所以b轴的轴压缩率最大。在铁镁橄榄石中，Fe的存在也使b轴的轴压缩率增加，因此b轴更容易压缩。

表 2 橄榄石的状态方程参数

Table 2. Equation of state parameters for olivines

Composition	V₀/(nm³)	K₀/(GPa)	p_max/(GPa)	Pressurecalibration	Pressuremedium	Method	Reference
Mg₂SiO₄	289.9(1)	151(2)^a	9.9	Au	ME	XRD(SC)	This study
Mg₂SiO₄	289.7(3)	132(6)^a	5.0	Olivine	─	XRD(SC)	Ref.[11]
Mg₂SiO₄	255	128(8)^a	41.8	Ruby	─	XRD(PD)	Ref.[21]
Mg₂SiO₄	290.14(9)	125(2)^a	17.2	Ruby	ME	XRD(SC)	Ref.[20]
Mg_1.84Fe_0.16SiO₄	292.5(1)	128.8(5)^a	8.0	Quartz	MEW	XRD(SC)	Ref.[22]
Mg_1.83Fe_0.17SiO₄	292.9(3)	140(3)^a	9.7	Au	ME	XRD(SC)	This study
Mg_1.8Fe_0.2SiO₄	290.4(5)	131(1)^a	9.7	Ruby	ME	XRD(SC)	Ref.[23]
Mg_1.66Fe_0.34SiO₄	269	134(10)^a	34	Ruby	─	XRD(PD)	Ref.[21]
Mg_1.6Fe_0.4SiO₄	294.6(1)	129.4(5)^a	7.4	Quartz	MEW	XRD(SC)	Ref.[4]
Mg_1.42Fe_0.58SiO₄	296.6(2)	129.4(5)^a	8.2	Quartz	MEW	XRD(SC)	Ref.[4]
Mg₂SiO₄		128.8(5)^b	16.2	Ruby	Ar/He/ME	BS	Ref.[8]
Mg_1.8Fe_0.2SiO₄	292.0(1)	131(2)^c	32.4	Quartz	ME	BS	Ref.[24]
Note:(1) Superscript letter a represents refitting K₀ with fixed K′₀=4; (2) Superscript letter b represents that K₀ was based on Brillouin scattering, and K′₀ was 4.2(2); (3) Superscript letter c represents that K₀ was based on Brillouin scattering, and K′₀ was 3.8(2); (4) ME:4:1 methanol-ethanol mixture; MEW:16:3:1 methanol-ethanol-water mixture; SC:single crystal; PD:powder; BS:Brillouin scattering.

下载: 导出CSV

| 显示表格

根据体积模量K和对应的密度ρ，通过v_Φ=(K/ρ)^1/2，可以得到样品的体波速v_Φ，如图 4所示。从图 4可以看出：镁橄榄石在零压下的体波速随Fe含量的增加呈现减小趋势；纯镁橄榄石样品在零压下的体波速(6.84km/s)大于铁镁橄榄石的体波速(6.50km/s)，两者相差约5.2%。Zha等人^[24]的研究结果表明，San Carlos橄榄石(Mg_1.8Fe_0.2SiO₄)在2.5GPa压力下的体波速为6.42km/s，在8.1GPa压力下为6.87km/s，与本研究得到的体波速6.68km/s(2.6GPa)和7.02km/s(8.4GPa)相差不大。由此可以得到：在镁橄榄石中，随着Fe组分的加入，体积模量和体波速均减小。另外，Andrault等人^[21]指出，橄榄石的弹性模量和弹性波速随着含水量的增加而逐渐减小。由于高压下橄榄石保持着高度的各向异性(见图 3)，因此在地球深部动力学(如地幔对流、洋壳深俯冲、地幔柱等)的作用下，上地幔的主要矿物橄榄石在特定区域内可能会沿着一定的方向择优排布，引起地震波的异常行为。本研究的结果显示，铁镁橄榄石比纯镁橄榄石更容易压缩，并且伴随着铁组分的加入，橄榄石的体波速减小。该结论对进一步了解上地幔的物质组成、蠕变机制及其动力学过程有一定的帮助。

图 4 常温下橄榄石的体波速随压力的变化

Figure 4. Bulk velocity as a function of pressure

下载: 全尺寸图片幻灯片

4. 结论

对两种天然镁橄榄石单晶进行了原位高压同步辐射X射线衍射研究，结果表明：在实验压力范围内镁橄榄石稳定存在，未发生相变。采用二阶Birch-Murnaghan状态方程对实验数据进行拟合，得到了铁镁橄榄石和纯镁橄榄石在零压下的晶胞体积(0.2929(3)、0.2899(1)nm³)和体积模量(140(3)、151(2)GPa)。两种镁橄榄石呈现高度各向异性，并且具有相同的轴压缩性规律，即b轴的压缩率最大，c轴次之，a轴的压缩率最小，在上地幔可能会沿着一定的方向择优排布。Fe对Mg的替换会减小镁橄榄石的体积模量，对比两种橄榄石的弹性性质发现，Fe组分的加入使橄榄石的体波速减小，该结果对于研究上地幔的组成及动力学过程具有一定的意义。

图 1 纯镁橄榄石和铁镁橄榄石单晶的XRD谱(方框所示为实验所得的单晶衍射点，数字表示对应的指标化结果；“D”表示金刚石的衍射点，衍射环为压标Au的衍射环)

Figure 1. Single-crystal XRD patterns of Fe-free forsterite and Fe-bearing forsterite (Squares show single-crystal diffraction peaksand figures indicate indexed results.'D' represents diamond reflections, while the diffraction rings in the patterns are from Au)

下载: 全尺寸图片幻灯片

图 2 铁镁橄榄石和纯镁橄榄石在常温下的p-V关系

Figure 2. p-V relationship for Fe-bearing forsteriteand Fe-free forsterite at 300K

下载: 全尺寸图片幻灯片

图 3 晶胞参数(a、b和c)随压力的变化

Figure 3. Lattice parameters a, b and c as a function of pressure

下载: 全尺寸图片幻灯片

图 4 常温下橄榄石的体波速随压力的变化

Figure 4. Bulk velocity as a function of pressure

下载: 全尺寸图片幻灯片

表 1 不同压力下铁镁橄榄石和纯镁橄榄石的晶胞参数

Table 1. Unit-cell parameters of Fe-bearing and Fe-free forsterite at various pressures

Material	p/(GPa)	a/(nm)	b/(nm)	c/(nm)	V/(nm³)
Fe-free forsterite	0.0001	0.47382(6)	1.0218(2)	0.59817(5)	0.28961(6)
	0.5	0.47283(4)	1.0236(1)	0.59743(4)	0.28917(3)
	1.2	0.47279(4)	1.0206(2)	0.59637(4)	0.28776(6)
	2.9	0.47124(4)	1.0168(2)	0.59438(4)	0.28481(6)
	4.0	0.47015(4)	1.0149(2)	0.59275(4)	0.28284(6)
	5.3	0.46926(4)	1.0126(2)	0.59106(4)	0.28060(6)
	6.2	0.46872(4)	1.0089(2)	0.58968(4)	0.27885(6)
	7.0	0.46810(4)	1.0085(2)	0.58839(4)	0.27777(6)
	7.7	0.46761(4)	1.0077(2)	0.58746(4)	0.27667(6)
	8.7	0.46677(4)	1.0043(2)	0.58692(4)	0.27514(6)
	9.9	0.46609(4)	1.0013(2)	0.58564(4)	0.27240(5)
Fe-bearing forsterite	0.0001	0.47832(2)	1.0150(1)	0.60201(5)	0.29230(4)
	2.6	0.47847(7)	1.0114(1)	0.59514(3)	0.28801(3)
	4.1	0.47840(7)	1.0070(2)	0.59284(2)	0.28561(3)
	4.8	0.47682(7)	1.0058(2)	0.59165(3)	0.28376(3)
	5.8	0.47538(8)	1.0047(2)	0.59022(2)	0.28200(3)
	7.0	0.47413(8)	1.0023(2)	0.58876(3)	0.27979(4)
	7.6	0.47407(8)	1.0000(2)	0.58805(3)	0.27880(3)
	8.4	0.47306(8)	0.9996(2)	0.58687(2)	0.27753(3)
	9.1	0.47199(9)	0.9984(3)	0.58630(2)	0.27628(4)
	9.7	0.47122(8)	0.9982(3)	0.58553(4)	0.27542(3)

下载: 导出CSV

表 2 橄榄石的状态方程参数

Table 2. Equation of state parameters for olivines

Composition	V₀/(nm³)	K₀/(GPa)	p_max/(GPa)	Pressurecalibration	Pressuremedium	Method	Reference
Mg₂SiO₄	289.9(1)	151(2)^a	9.9	Au	ME	XRD(SC)	This study
Mg₂SiO₄	289.7(3)	132(6)^a	5.0	Olivine	─	XRD(SC)	Ref.[11]
Mg₂SiO₄	255	128(8)^a	41.8	Ruby	─	XRD(PD)	Ref.[21]
Mg₂SiO₄	290.14(9)	125(2)^a	17.2	Ruby	ME	XRD(SC)	Ref.[20]
Mg_1.84Fe_0.16SiO₄	292.5(1)	128.8(5)^a	8.0	Quartz	MEW	XRD(SC)	Ref.[22]
Mg_1.83Fe_0.17SiO₄	292.9(3)	140(3)^a	9.7	Au	ME	XRD(SC)	This study
Mg_1.8Fe_0.2SiO₄	290.4(5)	131(1)^a	9.7	Ruby	ME	XRD(SC)	Ref.[23]
Mg_1.66Fe_0.34SiO₄	269	134(10)^a	34	Ruby	─	XRD(PD)	Ref.[21]
Mg_1.6Fe_0.4SiO₄	294.6(1)	129.4(5)^a	7.4	Quartz	MEW	XRD(SC)	Ref.[4]
Mg_1.42Fe_0.58SiO₄	296.6(2)	129.4(5)^a	8.2	Quartz	MEW	XRD(SC)	Ref.[4]
Mg₂SiO₄		128.8(5)^b	16.2	Ruby	Ar/He/ME	BS	Ref.[8]
Mg_1.8Fe_0.2SiO₄	292.0(1)	131(2)^c	32.4	Quartz	ME	BS	Ref.[24]
Note:(1) Superscript letter a represents refitting K₀ with fixed K′₀=4; (2) Superscript letter b represents that K₀ was based on Brillouin scattering, and K′₀ was 4.2(2); (3) Superscript letter c represents that K₀ was based on Brillouin scattering, and K′₀ was 3.8(2); (4) ME:4:1 methanol-ethanol mixture; MEW:16:3:1 methanol-ethanol-water mixture; SC:single crystal; PD:powder; BS:Brillouin scattering.

下载: 导出CSV

参考文献(24)

[1]	RINGWOOD A E.Composition and petrology of the Earth's mantle[M]. New York:McGraw-Hill Inc, 1975:138.
[2]	RINGWOOD A E, MAJOR A.Synthesis of Mg₂SiO₄-Fe₂SiO₄ spinel solid solutions[J]. Earth Planet Sci Lett, 1966, 1(4):241-245. doi: 10.1016/0012-821X(66)90077-X
[3]	AGEE C B.Phase transformations and seismic structure in the upper mantle and transition zone[J]. Rev Miner Geochem, 1998, 37(1):165-203.
[4]	NESTOLA F, PASQUAL D, SMYTH J R, et al.New accurate elastic parameters for the forsterite-fayalite solid solution[J]. Am Mineral, 2011, 96(11/12):1742-1747. http://adsabs.harvard.edu/abs/2011AmMin..96.1742N
[5]	ROUQUETTE J, KANTOR I, MCCAMMON C A, et al.High-pressure studies of (Mg_0.9Fe_0.1)₂SiO₄ olivine using Raman spectroscopy, X-ray diffraction, and Mössbauer spectroscopy[J]. Inorg Chem, 2008, 47(7):2668-2673. doi: 10.1021/ic701983w
[6]	SMYTH J R, FROST D J, NESTOLA F, et al.Olivine hydration in the deep upper mantle:effects of temperature and silica activity[J]. Geophys Res Lett, 2006, 33:L15301. doi: 10.1029/2006GL026194
[7]	IRIFUNE T, ISSHIKI M.Iron partitioning in a pyrolite mantle and the nature of the 410km seismic discontinuity[J]. Nature, 1998, 392:702-705. doi: 10.1038/33663
[8]	ZHA C S, DUFFY T S, DOWNS R T, et al.Sound velocity and elasticity of single-crystal forsterite to 16GPa[J]. J Geophys Res, 1996, 101(B8):17535-17545. doi: 10.1029/96JB01266
[9]	DUFFY T S, ZHA C S, DOWNS R T, et al.Elasticity of forsterite to 16GPa and the composition of the upper mantle[J]. Nature, 1995, 378:170-173. doi: 10.1038/378170a0
[10]	SUZUKI I.Thermal expansion of periclase and olivine, and their anharmonic properties[J]. J Phys Earth, 1975, 23(2):145-159. doi: 10.4294/jpe1952.23.145
[11]	HAZEN R M.Effects of temperature and pressure on the crystal structure of forsterite[J]. Am Mineral, 1976, 61:1280-1293. http://ndt.oxfordjournals.org/external-ref?access_num=10.1081/JDI-120005366&link_type=DOI
[12]	HAZEN R M.Effects of temperature and pressure on the crystal structure of ferromagnesian olivine[J]. Am Mineral, 1977, 62:286-295.
[13]	SPEZIALE S, DUFFY T S.Single-crystal elasticity of fayalite to 12GPa[J]. J Geophys Res, 2004, 109:B12202. doi: 10.1029/2004JB003162
[14]	LIU L, DU J G, LIU W, et al.Elastic behavior of (Mg_xFe_1-x)₂SiO₄ olivine at high pressure from first-principles simulations[J]. J Phys Chem Solids, 2010, 71(8):1094-1097. doi: 10.1016/j.jpcs.2010.03.013
[15]	ONO S, KIKEGAWA T, HIGO Y.In situ observation of a phase transition in Fe₂SiO₄ at high pressure and high temperature[J]. Phys Chem Miner, 2013, 40(10):811-816. doi: 10.1007/s00269-013-0615-3
[16]	FEI Y W, RICOLLEAU A, FRANK M, et al.Toward an internally consistent pressure scale[J]. Proc Natl Acad Sci, 2007, 104(22):9182-9186. doi: 10.1073/pnas.0609013104
[17]	DERA P.GSE-ADA data analysis program for monochromatic single crystal diffraction with area detector[M]. Argonne, IL:GeoSoilEnviro CARS, 2007.
[18]	HOLLAND T J B, REDFERN S A T.Unit cell refinement from powder diffraction data:the use of regression diagnostics[J]. Mineral Mag, 1997, 61(1):65-77. doi: 10.1180-minmag.1997.061.404.07/
[19]	BIRCH F.Finite elastic strain of cubic crystals[J]. Phys Rev, 1947, 71(11):809-824. doi: 10.1103/PhysRev.71.809
[20]	DOWNS R T, ZHA C S, DUFFY T S, et al.The equation of state of forsterite to 17.2GPa and effects of pressure media[J]. Am Mineral, 1996, 81:51-55.
[21]	ANDRAULT D, BOUHIFD M A, ITI J P, et al.Compression and amorphization of (Mg, Fe)₂SiO₄ olivines:an X-ray diffraction study up to 70GPa[J]. Phys Chem Miner, 1995, 22(2):99-107. doi: 10.1007/BF00202469
[22]	NESTOLA F, NIMIS P, ZIBERNA L, et al.First crystal-structure determination of olivine in diamond:composition and implications for provenance in the Earth's mantle[J]. Earth Planet Sci Lett, 2011, 305(1/2):249-255. http://www.sciencedirect.com/science/article/pii/S0012821X11001518
[23]	ZHANG L.Single crystal hydrostatic compression of (Mg, Mn, Fe, Co)₂SiO₄ olivines[J]. Phys Chem Miner, 1998, 25(4):308-312. doi: 10.1007/s002690050119
[24]	ZHA C S, DUFFY T S, DOWNS R T, et al.Brillouin scattering and X-ray diffraction of San Carlos olivine:direct pressure determination to 32GPa[J]. Earth Planet Sci Lett, 1998, 159(1/2):25-33. http://www.sciencedirect.com/science/article/pii/S0012821X98000636