CaF2熔化温度的分子动力学模拟

利用壳层模型分子动力学方法,研究了高温高压条件下CaF2的熔化温度,同时计算了温度为300K、压强上升到100GPa时CaF2 的状态方程.研究中考虑了分子动力学模拟的过热熔化,通过晶体的现代熔化理论,对CaF2 的分子动力学模拟熔化温度进行了修正, 获得了高温高压下CaF2的熔化温度.因此,常压下壳层模型分子动力学方法为研究物质熔化提供了一个很好的方法.     

壳层模型分子动力学在CaO研究中的应用

利用壳层分子动力学方法结合有效的对势,研究了高压条件下CaO的熔化曲线。研究表明,分子动力学模拟结果精确地再现了广泛压强范围内CaO的状态方程。研究中考虑了分子动力学模拟熔化存在的过热现象,通过晶体的现代熔化理论,对CaO的分子动力学模拟熔化温度进行了修正,获得了高温高压下CaO正确的熔化温度。因此,常压下引入壳层模型的分子动力学为研究物质熔化提供了一个很好的方法,这种方法可进一步推广到其它物质的高压熔化研究中。     

下地幔条件下钙钛矿结构MgSiO3熔化行为的研究

利用分子动力学方法结合有效的对势,模拟了下地幔条件下钙钛矿结构MgSiO3的熔化曲线．研究表明,分子动力学模拟结果精确地再现了广泛压强范围内钙钛矿结构MgSiO3的状态方程,并且熔化曲线与最新的实验结果也符合的很好．在压强上升到下地幔压强范围内,压强低于60 GPa时的钙钛矿结构MgSiO3熔化曲线比较陡,接着变得平缓．在核幔边界压强135 GPa时,钙钛矿结构MgSiO3的熔化温度是6500 K,明显低于Zerr和Boehler实验结果的外推结果．     

方镁石高压结构预测和高温结构稳定性研究

方镁石是镁方铁矿的终端组分,化学组成为氧化镁(MgO).理论预测的MgO熔化线和高压下实验测量结果之间存在巨大的分歧,为澄清歧见人们展开了对MgO高压结构的进一步研究,方镁石MgO高压结构预测及温度对结构稳定性的影响一直是高压凝聚态物理和地球物理领域的重要研究内容.本文利用基于密度泛函理论的第一性原理计算方法,对MgO实验结构、各种可能存在的结构及基于粒子群优化算法预测的晶体结构进行了系统深入的研究,发现MgO在0—580 GPa的压力范围内一直以稳定岩盐结构存在,580—800 GPa压力范围内的稳定结构为氯化铯结构.尽管NiAs六角密堆结构和纤锌矿结构能合理解释冲击压缩实验中沿MgO的P-V雨贡纽线在(170±10) GPa存在体积不连续的原因(Zhang L, Fei Y W 2008 Geophys.Res.Lett. 35 L13302)和高压下理论计算的熔化线与实验结果相差很大的原因(Aguado A, Madden P A 2005Phys.Rev.Lett.94 068501),但这两种结构连同闪锌矿结构及基于粒子群优化算法预测的晶体结构B8_2和P3m1仅为其亚稳结构.在MgO高压结构稳定性预测的基础上,本文利用经典分子动力学方法,分别引入能很好描述离子极化特性的壳层模型和离子压缩效应的呼吸壳层模型,对MgO岩盐结构的高温稳定性进行了模拟研究,给出了压力达150 GPa的高压熔化相图.     

Molecular Dynamics Simulations for Melting Temperatures of SrF2 and BaF2

利用壳层模型分子动力学方法,考虑萤石结构分子中的预熔化现象,对SrF₂和BaF₂的分子动力学模拟熔化温度进行修正,获得了高压下SrF₂和BaF₂的熔化温度. 同时给出了300 K、0.1 Mpa～7GPa和0.1 Mpa～3 GPa时SrF₂和BaF₂的状态方程,与已有研究结果的最大误差分别为0.3%和2.2%. 计算所得SrF₂和BaF₂常压下的熔点与已有的实验结果符合较好. 对于SrF₂和BaF₂分子体积变化和已有的熔化模拟的差别也做了比较和讨论.     

新的相互作用势对MD模拟钙钛矿结构MgSiO3热力学性质的影响

通过分析势能曲线解释了钙钛矿结构MgSiO3熔化模拟过程中模拟熔化温度存在较大差异的原因，并进一步研究了对势参数在分子动力学模拟中的影响. 通过调整已有的经验势得到了一组新的势参数，以此来进行分子动力学研究，得到的常温常压下摩尔体积与Belonoshko和Dubrovinsky的结果符合较好，并且其状态方程、常压下热容和常压下热膨胀系数与他人的实验值都较好地吻合. 另外，所得到的熔化温度也与以前的研究进行了比较.     

自由能方法与零压下Al的熔化温度

利用自由能方法的分子动力学模拟,计算了零压下Al的熔化温度.在计算液相自由能的过程中,采用勒纳-琼斯(LJ)液体作为参考系统,同时将计算结果与Mei和Davenport等人的计算结果进行了比较,计算结果表明：1)选用LJ参考系统使液相自由能的计算时间节省一半,并且不影响熔化温度的计算结果;2)采用不同的埋入原子势(EAM)的分子动力学模拟计算得到的熔化温度与实验值都存在偏差,而就金属Al而言,采用Cai等人的EAM势的熔化温度的计算结果比Mei和Davenport及Morris等人采用的势模型的结果略有改 关键词： 熔化温度 自由能方法 分子动力学模拟     

氦热力学性质的分子动力学研究

 利用经典分子动力学和第一性原理分子动力学,研究了氦在高压下的熔化曲线、状态方程和非金属-金属转变。得到了氦在温度小于4.5 eV、 密度0.3~5.0 g/cm³范围内的状态方程,并把氦的熔化曲线的压强范围拓展到了50 GPa。氦的能隙宽度曲线表明,温度大大降低了氦的金属化密度。     

Mo高压熔化的分子动力学模拟

 采用第一原理方法计算了钼在零温下的结构,表明钼在500 GPa以下一直保持bcc结构（常温）,与实验一致。在零压附近计算了E-V关系,利用Murnaghan物态方程拟合得到了零压体积及其模量,与实验结果符合得很好。采用第一原理分子动力学模拟了钼的高压熔化性质。采用NVT系综计算了128个原子的系统,初始构形为bcc结构,体积分别为0.015 48、0.012 19、0.010 98、0.009 84、0.009 10 nm³/atom,计算了几个温度点,拟合得到了熔化曲线,熔化温度明显高于金刚石压砧（DAC）实验结果;将初始构形改变为fcc结构,模拟其熔化特性,得到的熔化温度明显下降,与激光加载DAC实验结果一致,认为可能的原因是钼熔化后形成的液体结构类似于fcc结构,而不是常态时的bcc结构。     

高压下钙钛矿结构MgSiO3的分子动力学研究

利用分子动力学方法,研究了高温高压下钙钛矿结构MgSiO3的状态方程.研究表明,分子动力学模拟结果很好地再现了广泛温度和压强范围内钙钛矿结构MgSiO3的摩尔体积.温度300 K压强上升到120 GPa模拟的钙钛矿结构MgSiO3状态方程和有效的实验结果基本一致.在更高温度和更高压强下模拟的钙钛矿结构MgSiO3状态方程和他人的计算值吻合的很好.另外,还分别计算了温度300 K,900 K,1500 K和2500 K压强上升到120 GPa时MgSiO3的体积压缩率.     

Molecular dynamics study for the melting curve of MgO at high pressure

Shell-model molecular dynamics method is used to study the melting temperatures of MgO at elevated temperatures and high pressures using interaction potentials. Equations of state for MgO simulated by molecular dynamics are in good agreement with available experimental data. The pressure dependence of the melting curve of MgO has been calculated. The surface melting and superheating are considered in the correction of experimental data and the calculated values, respectively. The results of corrections are compared with those of previous work. The corrected melting temperature of MgO is consistent with corrected experimental measurements. The melting temperature of MgO up to 140GPa is calculated.     

Properties of MgO at high pressures: Shell-model molecular dynamics simulation

Shell-model molecular dynamics (MD) simulation has been performed to investigate the melting of the major Earth-forming mineral: periclase (MgO), at elevated temperatures and high pressures, based on the thermal instability analysis. The interatomic potential is taken to be the sum of pair-wise additive Coulomb, van der Waals attraction, and repulsive interactions. The MD simulation with selected Lewis–Catlow (LC) potential parameters is found to be very successful in describing the melting behavior for MgO, by taking account of the overheating of a crystalline solid at ambient pressure. The thermodynamic melting curve is estimated on the basis of the thermal instability MD simulations and compared with the available experimental data and other theoretical results in the pressure ranges 0–150 GPa. Our simulated melting curve of MgO is consistent with results obtained from Lindemann melting equation and two-phase simulated data at constant pressure by Belonoshko and Dubrovinsky, in the pressure below 20 GPa. The extrapolated melting temperatures in the lower mantle are in good agreement with the results obtained from Wang's empirical model up to 100 GPa. Compared with experimental measurements, our results are substantially higher than that determined by Zerr and Boehler, and the discrepancy between DAC and MD melting temperatures may be well explained with different melting mechanisms. Meanwhile, the radial distribution functions (RDFs) of Mg–Mg, O–Mg, and O–O ion pairs near the melting temperature have been investigated.     

Melting curve of MgO from first-principles simulations

First-principles calculations based on density functional theory, both with the local density approximation (LDA) and with generalized gradient corrections (GGA), have been used to simulate solid and liquid MgO in direct coexistence in the range of pressure 0 < or = p < or = 135 GPa. The calculated LDA zero pressure melting temperature is T(LDA)m = 3110 +/- 50 K, in good agreement with the experimental data. The GGA zero pressure melting temperature T(GGA)m = 2575 +/- 100 K is significantly lower than the LDA one, but the difference between the GGA and the LDA is greatly reduced at high pressure. The LDA zero pressure melting slope is dT/dp approximately 100 K/GPa, which is more than 3 times higher than the currently available experimental one from Zerr and Boehler [Nature (London) 371, 506 (1994)]. At the core mantle boundary pressure of 135 GPa MgO melts at Tm = 8140 +/- 150 K.     

The high-pressure melting curve of CaO

Shell-model molecular dynamics simulation has been performed to investigate the melting of the major Earth-forming mineral CaO at elevated temperatures and high pressures, based on thermal instability analysis. The interatomic potential is taken to be the sum of effective pair-wise additive Coulomb, van der Waals attraction, and repulsive interactions. It is shown that the simulated molar volume of CaO is successful in reproducing recent experimental data and our DFT-GGA calculations up to the core–mantle boundary pressure of 135 GPa. The pressure dependence of the simulated high pressure melting temperature of CaO is in good agreement with the results obtained from the Lindemann melting equation at a pressure of below 7 GPa. The extrapolated melting temperatures are in good agreement with the results obtained from Wang’s empirical model up to 60 GPa. The predicted high pressure melting curve, being very steep at lower pressures, rapidly flattens on increasing pressure. The thermodynamic properties of the rocksalt phase of CaO are summarized in the 0–135 GPa pressure range and for temperatures up to 9300 K.     

高压下MgO的热弹性研究

利用第一原理平面波赝势方法和广义梯度近似研究了广泛温度和压强范围内MgO的热弹性特性．MgO从低压NaCl到高压CsCl结构的相变压强为397 GPa，表明它在地球内部不会发生相变．在压强上升到150 GPa时，MgO的绝热弹性模量跟0 GPa时的实验值和其他赝势在高压下的计算结果基本一致．从0 GPa上升到20 GPa时，MgO的各向异性逐渐减小；在20?150 GPa时绝对值逐渐增大．MgO明显地违背了Cauchy条件，反应出非中心多体力是重要的．另外，MgO的热力学参量与实验也符合的很好．     

Equations of state of MgO at high pressure and temperature

Abstract

A synchrotron X-ray diffraction study on MgO has been done at simultaneous high pressure and temperature. The lattice parameter of MgO has been measured up to a static pressure of 6 GPa and a temperature of 1273 K, using a large volume pressure cell and energy-dispersive synchrotron X-ray powder diffraction. The compression was made following six high-temperature isotherms. A Vinet equation of state was used to fit the experimental P-V-T data. The Vinet's model compares very well with the experimental data above the Debye temperature (760 K) and allows the use of MgO as an alternative internal pressure calibrant for experiments at high temperature.     

高压下块状Pd40Ni40P20非晶合金的晶化

 本文研究了Pd₄₀Ni₄₀P₂₀块状非晶在4 GPa及常压下的晶化过程。得到了时间-温度转变图。结果表明：高压下样品的晶化温度明显升高，压力对原子的长程扩散及相分离熔体的粘性流动均有抑制作用。在接近熔点进行高压退火时，获得了单相过饱和固溶体。其晶体结构为面心立方。     

无氧铜冲击熔化温度的辐射法测量

 在发生冲击熔化的情况下，金属样品/窗口界面压力下的熔化温度与卸载温度数值相等，且十分接近于界面温度值。根据这一结论，利用二级轻气炮加载手段和光辐射法测温技术，用氟化锂（LiF）单晶作透明窗口，获得了110~140 GPa压力范围内无氧铜的熔化温度。实验表明，无氧铜的高压熔化温度数据与文献发表的无氧铜高压声速实验结果是一致的，铜的高压熔化规律可用Lindemann熔化定律近似描述。采用的熔化温度测量方法不必反演出冲击温度，简化了冲击熔化温度的数据处理方法，为金属冲击熔化温度测量提供了一种潜在的技术途径。     

High P–T Brillouin scattering study of H2O melting to 26 GPa

The properties of solid and liquid phases of H₂O at high pressure and temperature remain an active area of research. In this study, Brillouin spectroscopy has been used to determine the temperature dependence of sound velocities in H₂O as a function of pressure up to 26 GPa through the phase field of ice VII and into the liquid to a maximum temperature of 1200 K. The Brillouin shift of the quasi-longitudinal acoustic mode moves to lower frequencies upon melting at each pressure. As a test of the method, measurements of the melting of Ar by Brillouin scattering at several pressures show a similar behavior for the acoustic mode, and measured melting points are consistent with previous results. The results of H₂O melting are consistent with previously reported melting curves below 20 GPa. The data at higher pressure indicate that ice melts at a higher temperature than a number of previous studies have indicated.