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

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

Molecular dynamics of MgSiO3 perovskite melting

The melting curve of MgSiO分子动力学 MgSiO3钙钛矿 熔化温度 高压melting temperature, molecular dynamics, high pressureProject supported by the National Natural Science Foundation of China (Grant Nos 10274055 and 10376021),the Natural Science Foundation of Gansu Province, China (Grant No 3ZS051-A25-027) and the Scientific Research Foundation of Education Bureau of Gansu Province, China (Grant No 0410-01).2005-01-125/8/2005 12:00:00 AMThe melting curve of MgSiO3 perovskite is simulated using molecular dynamics simulations method at high pressure. It is shown that the simulated equation of state of MgSiO3 perovskite is very successful in reproducing accurately the experimental data. The pressure dependence of the simulated melting temperature of MgSiO3 perovskite reproduces the stability of the orthorhombic perovskite phase up to high pressure of 130GPa at ambient temperature, consistent with the theoretical data of the other calculations. It is shown that its transformation to the cubic phase and melting at high pressure and high temperature are in agreement with recent experiments.     

MgSiO3钙钛矿高温高压特性的第一性原理研究

本文采用第一性原理方法,计算了MgSiO3钙钛矿在零温和0 ~ 150 GPa静水压范围内的晶体结构和弹性模量,并利用准简谐近似Debye模型,拟合三阶Birch-Murnaghan物态方程得到了其高温高压下的热力学性质。通过与现有的理论和实验的结果数据比较,确认在0 ~ 2000 K的温度区间内,第一性原理计算结合Debye模型能够较可靠地模拟在下地幔压力范围内MgSiO3钙钛矿的热力学性质。     

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

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

高压下MgO熔化特性的分子动力学模拟

利用分子动力学方法,模拟研究了高压下MgO的熔化特性．通过晶体的现代熔化理论,对MgO的分子动力学模拟熔化温度进行了修正,得到了高温高压下MgO的熔化温度．计算得到的MgO熔化曲线和已有的实验及其它理论结果在0-135 GPa进行了比较,发现修正得到的MgO熔化温度和由Lindemann熔化方程及两相方法得到的结果在压力低于15 GPa时符合很好．同时,MgO熔化模拟有效解释了一阶相变分子动力学过程中出现的过热熔化现象．     

MgSiO3状态方程的分子动力学模拟

利用分子动力学方法,研究了高温高压下钙钛矿结构MgSiO_3的状态方程。研究表明,分子动力学模拟结果精确地再现了广泛温度和压强范围内MgSiO_3的摩尔体积。在300 K压强上升到140 GPa模拟的MgSiO_3状态方程和有效的实验值、他人的拟合值以及基于局域密度近似的第一原理计算结果基本一致。并且更高温度和更高压强下模拟的MgSiO_3状态方程和他人的计算值吻合的很好。另外,还分别计算了300、900、2000和3000 K压强上升到120 GPa时MgSiO_3的体积压缩率。     

The melting curve of CaSiO3 perovskite under lower mantle pressures

The high pressure melting curve of CaSiO₃ perovskite is simulated by using the constant temperature and pressure molecular dynamics method combined with effective pair potentials for the first time. The simulated results for the partial radial distribution function all compare well with experiment. The calculated equation of state is very successful in accurately reproducing the recent experimental data over a wide pressure range. The predicted high pressure melting curve is in good agreement with the experimental ones, and the melting curve up to the core–mantle boundary pressure, being very steep at lower pressures, rapidly flattens on increasing pressure. The present results also suggest the validity of the experimental data of Zerr and Boehler.     

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.     

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.     

Melting of the rare earth metals and f-electron delocalization

Melting curves for Pr, Nd, Sm, Gd, and Y were measured in a diamond-anvil-cell to nearly 100 GPa and 4000 K. f-electron volume collapses are observed as triple points for Pr (24 GPa and 1400 K) and Gd (65 GPa and 3100 K). These pressures coincide with the volume collapses observed at room temperature. For Nd and Sm, the f-electron volume collapse has not been observed at room temperature but appears at approximately 2000-2500 K as a broad minimum in the melting curve, similar to that of Ce, near 50 GPa (Nd) and 70 GPa (Sm). The melting curve of Y goes smoothly along the entire rare earth sequence.     

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

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

Evidence of maximum in the melting curve of hydrogen at megabar pressures

Hydrogen at high pressures of ∼400 GPa might be in a zero-temperature liquid ground state (N. Ashcroft, J. Phys.: Condens. Matter A 12, 129 (2000), E. G. Brovrnan et al., Sov. Phys. JETP 35, 783 (1972)). If metallic hydrogen is liquid, the melting T _melt(P) line should possess a maximum. Here we report on the experimental evaluation of the melting curve of hydrogen in the megabar pressure range. The melting curve of hydrogen has been shown to reach a maximum with T _melt = 1050 ± 60 K at P = 106 GPa and the melting temperature of hydrogen decreases at higher pressures so that T _melt = 880 ± 50 K at P = 146 GPa. The data were acquired with the aid of a laser heating technique where diamond anvils were not deteriorated by the hot hydrogen. Our experimental observations are in agreement with the theoretical prediction of unusual behavior of the melted hydrogen [S. Bonev et al., Nature 481, 669 (2004)]. The article is published in the original.     

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.     

Melting of cubic boron nitride at extreme pressures

Due to its large pressure range of stability and inert nature, cubic boron nitride has been proposed as a potential pressure standard for high pressure experiments. It is extremely refractive upon compression, although its melting temperature is not known beyond 10 GPa. We apply first-principles molecular dynamics to evaluate the thermodynamics of zincblende structured (cubic) and liquid boron nitride at extreme temperatures and pressures, and compute the melting curve up to 1 TPa by integration of the Clapeyron equation. The resulting equations of state reveal that liquid boron nitride becomes denser than the solid phase at pressures of around 0.5 TPa. This is expressed as a turnover in the melting curve, which reaches a maximum at 510 GPa and 6550 ± 700 K. The origin of this density crossover is explained in terms of the underlying liquid structure, which diverges from that of the zincblende structured solid as the phases are compressed.     

Structure of liquid iron at pressures up to 58 GPa

We report structural data on liquid iron at pressures up to 58 GPa measured by x-ray scattering in a laser heated diamond anvil cell. The determined structure factor preserves essentially the same shape along the melting curve. Our data demonstrate that liquid iron at high pressures is a close-packed hard-sphere liquid. The results place important constraints on the thermodynamic and transport properties of liquid iron and the melting curve of iron.     

Equations of State and the Phase Diagram of Dielectric and Metallic Helium-4

A new form of the semiempirical equation of state for the liquid phase of helium-4 is proposed that is based on the assumption that the structure of this phase consists of a mixture of dielectric and metallic components. It is postulated that solid dielectric helium with density higher than 5.3 g/cm³ becomes a metal. The values of the parameters of the equations of state for both solid phases and the liquid phase of helium are calculated. The unknown values of the initial data for helium are taken by analogy with the parameters for deuterium. The phase diagram, shock adiabat, isentropes, isotherms, and electrical conductivity in these processes are calculated with the use of the equations of state of solid and liquid phases of helium-4. The results of calculation are compared with experimental data in the range of pressures of up to 35 GPa for an isotherm, up to 150 GPa for a shock adiabat, up to 42 GPa for the melting curve, and up to 2000 GPa for isentropes, and showed quite satisfactory agreement. Numerical extrapolation of the melting curve is performed to a range of ultrahigh pressures of up to 8000 GPa.     

Phase transitions in uranium dioxide at high pressures and temperatures

The melting curve is plotted for uranium dioxide with fluorite structure in a pressure range from ?2.5 to +100 GPa. This curve has a peak at the point 3348 K, 6 GPa, and has a negative derivative at high pressures. The pressure corresponding to a polymorphic transition of uranium dioxide (37 GPa) at a temperature of 1015 K is determined. The slope of the equilibrium curve of the polymorphic transition in UO₂ in the temperature range 300–1000 K is ? 56 K/GPa.     

一种高密度玻璃的多形性高压相变和物态方程研究

在二级轻气炮上用无氧铜飞片直接撞击重玻璃平板样品(密度为4.817 g/cm³,材料牌号:ZF6)开展了冲击压缩实验研究,压力范围为52.1—167.8 GPa,并采用多通道瞬态辐射高温计和光分析技术测量了其雨贡纽线、高压声速和冲击波温度等动态特性.实验结果显示,上述性质在三个不同压力区间出现不连续性变化,表明冲击压缩下该样品材料存在多形性高压相变,相变起始压力分别为23,78和120 GPa.实测声速先是随冲击压力的增高而增加,并在78 GPa附近出现急剧下降,之后又随压力增长,并在120 GPa之后下降到体波声速,表明材料进入高压熔化相.温度数据同样在78和120 GPa处出现明显的不连续变化,并在120 GPa之后变化趋于平缓与计算的Lindeman熔化线相符,进一步印证了上述相变行为.实测雨贡纽数据与LASL数据库中的重玻璃数据相符,结果显示除23 GPa附近有一明显的突变外,高压区数据几近线性变化,表明重玻璃的两个高压相变均为二级相变.本文报道的重玻璃材料高压物性数据和序列相变认识对于发展反向加载技术、提高材料声速测量精度和适用压力范围具有实用价值. 关键词： 重玻璃 冲击温度 卸载声速 冲击相变     

顽火辉石(Mg0.92,Fe0.08)SiO3的冲击相变和高压状态方程及其地球物理意义

 用阻抗匹配法和电探针技术在48~140 GPa冲击压力范围内对化学组分为(Mg_0.92, Fe_0.08)SiO₃、初始密度为3.06 g/cm³的天然顽火辉石进行了冲击压缩实验。根据本工作13发实验数据，结合McQueen等人的数据可以看出，(Mg_0.92, Fe_0.08)SiO₃顽火辉石在冲击压缩过程中，大约经历三个明显区域：低压相区，压力范围为0~40 GPa；混合相区，压力范围为40～67 GPa；高压相区，压力范围为68～140 GPa。在低压相区，D-u关系已由McQueen给出；而在高压相区（68～140 GPa），可由本实验数据得到。由叠加原理计算得到的混合物(Mg_0.92, Fe_0.08)O(Mw)+SiO₂(St)的D-u关系及p-ρ关系曲线明显偏离了实验数据的拟合曲线，从而排除了在高达140 GPa冲击压力下，钙钛矿结构的(Mg_0.92, Fe_0.08)SiO₃发生向氧化物化学分解相变的可能性。对高压相区的实验数据进行拟合，可以得到(Mg_0.92, Fe_0.08)SiO₃钙钛矿的Grüneisen参数γ。通过三阶Birch-Murnaghan有限应变状态方程，由冲击波实验数据得到了零压等熵体积模量K_0S=259.6(9) GPa及其对压力的一阶偏导数K′_0S=4.20(5)，其ρ₀=4.19 g/cm³。(Mg_0.92, Fe_0.08)SiO₃钙钛矿冲击压缩下的密度数据与PREM密度剖面吻合很好，支持钙钛矿为主要成分的下地幔模型。     

High-pressure melting curve of nitrogen and the liquid-liquid phase transition

The melting curve of nitrogen was measured up to 71 GPa, a fourfold increase in pressure over previous measurements. The measurements were made using the laser-heated diamond anvil cell and melting was detected in situ by the laser speckle method. The melting temperature rises linearly up to a maximum at 50 GPa and 1920 K, and with increasing pressure suddenly decreases linearly to 1400 K at 71 GPa. This sharp drop in the melting slope (dT/dP) above 50 GPa indicates the appearance of a liquid denser than the solid and of a liquid-liquid phase transition. The sharpness of the changes suggests that the transition is first order and is a liquid-liquid polymer transition. This conclusion is consistent with earlier theoretical studies and experimental evidence that pressure transforms molecular nitrogen into a chainlike polymeric form.