博士论文答辩报告-演讲版.pptVIP

热烈欢迎各位老师、同学光临 下地幔主要矿物的熔化、弹性和热力学特性研究 　答辩人：　刘子江　　　　　　 　指导老师：杨向东　教授 报　告　内　容 ? 研究背景 ? 计算方法 ? 从头算方法 ? 分子动力学方法 ? 研究内容 ? MgSiO3钙钛矿和MgO的高压熔化曲线 ? 下地幔主要矿物的高压弹性行为 ? CaSiO3钙钛矿的热力学特性 ? 结论 研究背景之研究意义 研究背景：地球的演化及其内部结构 研究背景：地球的演化及其内部结构 研究背景：球形地震模型　　PREM (preliminary reference Earth model) 研究背景：地球下地幔化学元素的丰度 研究背景：地球下地幔的主要成分 研究背景：矿物学与其他学科之间的关系 研究背景：下地幔主要矿物的性质 矿物高压熔化曲线的研究意义 熔化研究进展 ? 实验研究的难度 ? 熔化曲线的理论研究 ? 熔化曲线的模拟 实验的熔化判据 ? 静高压（金刚石压砧技术－DAC）：表面织构判剧－用激光对样品进行加温加压，用显微镜观察晶体表面是否发生对流或变形来判断材料是否熔化。 ? 动高压（冲击波加载技术）：用冲击波加载的方法测量样品在不同冲击压力下的冲击温度，从冲击温度与压力曲线的变化趋势来确定样品的冲击熔化线。 两种实验方法之间的误差 　 。 静高压实验结果 晶体熔化的理论模型 ? 熔化的两相理论(J. R. Morris, Phys. Rev. B, V49 (1994) 3109) ? Lindemann熔化判据(Lindemann F A, Z. Physik. 11 (1910) 609) 1910年Lindemann指出，晶体中原子振动的均方根位移达到原子最近邻距离的某一临界比时熔化发生, 据此判剧和Einstein模型计算的单质的熔化温度场Debye温度有关: Tm= CmV1/2Θ2其中C为常数, m为质量, V为体积, Θ为晶体的Debye温度. ? 力学不稳定判据(Born　M, J. Chem. Phys. 7 (1939) 591)：晶体具有抵抗剪切应力的弹性阻力，而液体则没有，所以切变模量（或其中之一）的消失导致晶体的熔化。（Tallon and Born）固体在它的模量能够等温地转变为零的状态时熔化。 ? 热弹性失稳判据(Boyer L L, Phase Tran. 5 (1985) 1)：固体随温度的升高引起切变力的阻抗降低, 当固体熔化时, 热膨胀和等温压缩率趋于无穷大. ? 缺陷熔化判据(Gorecki T, Z. Metall, 65 (1974) 426)：缺陷的浓度随着温度升高而达到某一临界值时导致熔化. 理论模型、静高压和动高压熔化曲线 CaSiO3 钙钛矿的结构 MgO的相图和相变压 计算方法 第一性原理计算方法 从头算程序的结构 分子动力学程序的结构 研究内容：MgO的高压熔化曲线 研究内容：MgSiO3的高压熔化曲线 0压下计算的绝热弹性模量与实验的比较 MgO的相转变压强 MgO的弹性模量 MgO的弹性各向异性与压强的关系 CaSiO3钙钛矿的各向同性压缩波速和剪切波速 后钙钛矿相 Motohiko Murakami, Kei Hirose et al (Science,May 2004) announce new post-perovskite phase. Also found by Oganov Ono (Nature, July 2004) Based on “CaIrO3” or “UFeS3” structure. Stable at P120 GPa. D” will not be made of perovskite at all! MgSiO3钙钛矿的弹性常数与压强的变化关系 MgSiO3钙钛矿与后钙钛矿弹性模量的比较 MgSiO3钙钛矿和后钙钛矿相的各向同性压缩波速和剪切波速 ?利用经典分子动力学和Buckingham对势模型，确定了MgSiO3钙钛矿和MgO的高压熔化曲线。 ?利用第一性原理平面波赝势方法，研究了MgSiO3钙钛矿同质异构体、MgO和CaSiO3钙钛矿的高压弹性行为。 ?利用准谐近似Debye模型，系统地预测了CaSiO3钙钛矿高温高压下的热力学特性。 展望…… ?本文利用计算机模拟方法成功地研究了极端条件下地球下地幔三个主要矿物MgSiO3钙钛矿、MgO和CaSiO3钙钛矿的高压行为。而实际上富Mg硅酸盐还含有少量的铁，MgO也只是镁方铁矿（Mg,Fe）O的端元。因此，为了更好地理解地球下地幔主要成分的行为，就需要对（Mg,Fe）SiO3钙钛矿和（Mg,Fe）O的