论文部分内容阅读
（Mg，Fe）SiO3钙钛矿作为下地幔最主要的组成矿物，其高压物性（如状态方程，热力学稳定性等）的研究对于准确限定地幔矿物组成，了解地幔动力学过程（如地幔对流模式，地震火山形成机理等）具有十分重要的意义。本学位论文采用静高压与冲击压缩技术相结合的手段，并辅以计算机模拟技术，对（Mg，Fe）SiO3钙钛矿开展了高温高压状态方程及相稳定性的理论、实验研究。这些研究为建立准确、合理的下地幔矿物学组成模型提供了确凿的数据。提升了对于地球深部物理的了解与认识。
本学位论文还利用静高压大腔体技术对地核的可能组成元素体系Fe—S—C开展了一系列高温高压熔化实验研究。基于对淬火样品的定量定性分析，初步掌握该三元合金体系熔化关系随压力、温度、初始组份的变化规律。这对于推测星体演化过程中核元素的分布及建立行星核矿物组成模型有极为重要的意义。
本文的主要研究内容和认识如下：
（1）通过对大腔体压机样品装配腔的设计与改进，利用活塞圆筒压机和大压机成功合成了适合冲击波实验尺寸要求的大块钙钛矿MgSiO3样品。通过微区电子探针、拉曼光谱和X射线衍射分析证实合成材料确为钙钛矿MgSiO3，为下一步的冲击波高压实验研究提供了所必须的初始样品，使得在高压下对钙钛矿MgSiO3开展直接冲击压缩研究成为可能。
（2）在47～107 GPa冲击压力（估算温度为600～3300K）范围内进行了初始样品为钙钛矿MgSiO3的冲击压缩状态方程测量实验。这是国际上首次开展的以钙钛矿MgSiO3为冲击样品的动高压实验，是本文的创新点之一。得到冲击波速度（Us）与粒子速度（up）关系： Us=6．42+1．48up。通过Rankin—Hugoniot方程拟合得到了Gruneisen参数γ0=1．33（q=1）。利用三阶欧拉有限应变方程拟合得到的绝热体积模量Kos及其对压力的一阶偏导数Kos’分别是254（±10） GPa和3．9（±0．17）。通过与以顽火辉石为初始样品的冲击数据的对比，发现：MgSiO3钙钛矿在下地幔温压条件下无化学分解反应发生，并且90GPa冲击压力以上时MgSiO3顽火辉石才彻底转变为钙钛矿相。本文的数据完全排除了低压相顽火辉石的影响，首次提供了通过直接冲击压缩钙钛矿所得到的状态方程参数，为今后构建地幔矿物学组成模型提供了更准确的热力学约束。也为进一步利用冲击手段研究地幔高压物性做了有益的探索。
（3）在3．5—20GPa，1100—1700℃范围内开展了Fe—C—S三元体系的高温高压熔化实验研究，通过对淬火样品纹理的仔细观察，并结合扫描电镜、电子探针等测试分析手段，考察了压力、温度、组份对该体系高压相关系的影响。低压下，该体系将出现富S熔体与富C熔体的不互溶现象。随着压力的升高，不互溶液体渐渐融合，并最终在4．9GPa成为一个化学均匀的熔体。压力的改变还促使了一系列复杂中间相的生成，共熔温度随之变化。在等压线上随着温度的降低，Fe—C合金首先结晶，继而是熔点较低的Fe—S合金结晶出溶。使用两种不同的初始样品导致固相线下形成的矿物相迥异。实验结果对于了解在行星演化过程中，（富含C，S）行星核的分化历史及对核元素组成的限定具有非常重要的意义：当行星核分化过程中核压小于SGPa时，由于熔体的不互融，将出现分层的核，即富S的液体外核及富C的液体内核。随着核的冷却，Fe—C合金将慢慢析出形成固体内核。而对于核压大于5GPa的行星体，在其核内将无分层现象。特别地，基于本文的研究推论，本文提出了含C的固体内地核及富含S（贫C）的液态外地核模型。
（4）采用密度泛函理论和平面波赝势方法针对MgSiO3钙钛矿开展了第一性原理计算研究。模拟得到了0—120GPa范围内MgSiO3钙钛矿零温昌格常数、密度、声速等弹性数据，模拟得到的状态方程参数（K0=264GPa，K0’=3．83）与实验数据较好符合，证实了本文所采用的计算方法的可靠性。考虑到真实下地幔的高温环境，利用热力学方法对模拟得到的声速、密度、体模量进行了温度修正，并将修正后的数据与地球初步参考模型的相应剖面比较，对比结果很好的支持了以MgSiO3钙钛矿为主要组成矿物的下地幔矿物组成模型。
（5）针对MgSiO3和（Mgx，Fe1—x）SiO3钙钛矿在高压下可能存在的微结构相变，采用密度泛函理论和平面波赝势方法分别对MgSiO3（Pbnm—P4mbm—Pm3m）及（Mgx，Fe1—x）SiO3（Pm—Pmmm—P4mmm）的可能相变序列在0—150GPa内开展了第一性原理计算模拟研究。基于在同一压力下不同相的焓差的比较，钙钛矿MgSiO3在高压下的相稳定序列（由稳定到次稳定）是：Pbnm（orthorhombic）—P4mbm（tetragonal）—Pm3m（cubic）．而对于（Mg0．75， Fe0．25）SiO3来说，相应的相稳定序列是pm（orthorhombic）—pmmm（tetragonal）—p4mmm（cubic）。计算结果支持具有低对称性的正交MgSiO3较其他相结构在下地幔压力条件下更为稳定，此外，在MgSiO3加入少量Fe将使得其稳定性降低。