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

Application of shell model in molecular dynamics simulation to MgO

The P－V－T equation of state of MgO has been simulated under high pressure and elevated temperature using the molecular dynamics (MD) method with the breathing shell model (BSM). It is found that the MD simulation with BSM is very successful in reproducing accurately the measured molar volumes of MgO over a wide range of temperature and pressure. In addition, the MD simulation reproduces accurately the measured volume compression data of MgO up to 100GPa at 300K. It is demonstrated that the MD simulated P－V－T equation of state of MgO could be applied as a useful internal pressure calibration standard at elevated temperatures and high pressures.     

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 and Grüneisen parameters of NaCl at high pressure

The Buckingham potential has been employed to simulate the melting and thermodynamic parameters of sodium chloride (NaCl) using the molecular dynamics (MD) method. The constant-volume heat capacity and Grüneisen parameters have been obtained in a wide range of temperatures. The calculated thermodynamic parameters are found to be in good agreement with the available experimental data. The NaCl melting simulations appear to validate the interpretation of superheating of the solid in the one-phase MD simulations. The melting curve of NaCl is compared with the experiments and other calculations at pressure 0－30GPa range.     

Effect of pressure on the thermal expansion of MgO up to 200~GPa

Constant temperature and pressure molecular dynamics (MD) simulations are performed to investigate the thermal expansivity of MgO at high pressure, by using effective pair-wise potentials which consist of Coulomb, dispersion, and repulsion interactions that include polarization effects through the shell model (SM). In order to take into account non-central forces in crystals, the breathing shell model (BSM) is also introduced into the MD simulation. We present a comparison between the volume thermal expansion coefficient α dependences of pressure P at 300 and 2000~K that are obtained from the SM and BSM potentials and those derived from other experimental and theoretical methods in the case of MgO. Compared with the results obtained by using the SM potentials, the MD results obtained by using BSM potentials are more compressible. In an extended pressure and temperature range, the α value is also predicted. The properties of MgO in a pressure range of 0--200~GPa at temperatures up to 3500~K are summarized.     

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

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

Melting of dense sodium

High-pressure high-temperature synchrotron diffraction measurements reveal a maximum on the melting curve of Na in the bcc phase at approximately 31 GPa and 1000 K and a steep decrease in melting temperature in its fcc phase. The results extend the melting curve by an order of magnitude up to 130 GPa. Above 103 GPa, Na crystallizes in a sequence of phases with complex structures with unusually low melting temperatures, reaching 300 K at 118 GPa, and an increased melting temperature is observed with further increases in pressure.     

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

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

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₂分子体积变化和已有的熔化模拟的差别也做了比较和讨论.     

相变材料热物理性质的分子动力学模拟

为从微观尺度探寻相变材料的热物性变化机理, 本文采用分子动力学的方法, 构建了由正二十二烷组成的无定形结构的相变材料体系, 采用周期性边界条件以及COMPASS力场对相变材料的比热以及导热系数进行了模拟, 并对纯正二十二烷进行了DSC测试. 结果表明, 模拟所得的相变材料热容与文献实验值的偏差是6.5%, 熔点与DSC实验值的偏差是0.98%. 当温度为288–318 K时, 相变材料的导热系数在0.1–0.4 W·m^-1·K^-1 范围内波动, 且随着压力增大略呈下降趋势. 关键词： 扩散系数 比热 导热系数 分子动力学     

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

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

Simulated equation of state of CaF2 with fluorite-type structure at high temperature and high pressure

The pressure-volume-temperature (P-V-T) equation of state (EOS), isothermal bulk modulus, and thermal expansivity of CaF₂ with cubic fluorite-type structure are investigated using the constant temperature and pressure shell model molecular dynamics (MD) method with effective pair potentials which consist of the Coulomb, dispersion, and repulsion interaction. It was shown that MD simulation is very successful in accurately reproducing the measured volumes of the CaF₂ over a wide range of pressures. The simulated P-V data matched X-ray diffraction experimental results up to 9.5 GPa at 300 K. In addition, volume thermal-expansion coefficient and isothermal bulk modulus were also calculated and compared with available experimental data and the latest theoretical results at ambient condition. At extended temperature and pressure ranges, The P-V EOS under different isotherms at selected temperatures, T-V EOS under different isobars at selected pressures, thermal expansivity, and isothermal bulk modulus were predicted up to 1500 K and 10 GPa. The detailed knowledge of thermodynamic behavior and EOS at extreme conditions are of fundamental importance to the understanding of the physical properties of CaF₂.     

Crystal-liquid phase relations in silicon at negative pressure

The stable and metastable melting relations for silicon in the diamond and Si136 clathrate-II structures at positive and negative pressures are calculated by molecular dynamics computer simulation. The simulated liquid and crystalline clathrates undergo cavitation at approximately -3 and -12 GPa. Between these limits a stretched crystal would transform directly to gas in response to a mechanical instability. Most importantly, the clathrate-II crystal becomes thermodynamically stable over the diamond at negative pressure below -1 GPa at the melting point. Si136 should then crystallize from a slightly stretched liquid, which would have the same volume as a diamond-structure crystal.     

A new differential thermal analysis setup for measuring high pressure phase transitions

A new setup has been recently developed in the toroidal opposed-anvil device ‘Conac 40’ to perform differential thermal analysis in a high pressure range (0–6 GPa). To evaluate the precision and the reliability of the setup, the high pressure melting curve of germanium and the transition points of α-iron have been investigated up to 5.3 GPa and compared with previous results.     

Investigation of the effect of pressure on thermodynamic properties and thermoelastic phase transformation of CuAlNi alloys: A molecular dynamics study

Thermoelastic phase transformations and thermodynamic properties of CuAlNi alloys at 0, 1, 2 and 3 GPa pressures were investigated by using MD simulation in this study. The interactions between atoms were modelled by Sutton-Chen type of embedded atom method (SCEAM) that is based on many-body interaction. It was observed that thermoelastic phase transformation in the ternary alloy system occurred at the end of thermal process. Radial distribution function (RDF) was used in order to analysis the structures obtained from MD simulation using the simulation techniques’ thermodynamic parameters. The transformation temperatures, enthalpy and entropy of the ternary alloy system have been observed to be changing with the applied pressure. In addition, it was found that the elastic energy has been decreased about 22% by applied pressure whereas Gibbs free energy has been increased about 60% by applied pressure. The values of the thermodynamical parameters obtained in this study were observed to be in close agreement with the experimental study.     

The p–T phase diagram of KNbO3 by a dielectric constant measurement

A dielectric constant measurement was carried out on perovskite-type ferroelectrics KNbO₃ over a wide range of temperature under high pressure. The temperature- and pressure-dependence of the dielectric constant clarified that all temperatures of the transitions from the ferroelectric rhombohedral to orthorhombic, to tetragonal and then to the paraelectric cubic phase, decrease with increasing pressure. These results indicate that the orthorhombic–tetragonal transition takes place at 8.5 GPa and the tetragonal–cubic transition at 11 GPa, at room temperature.     

High-pressure differential thermal analysis measurements of the melting curve of lithium

The pressure effect on the melting behavior of lithium has been measured by observing the latent heat signatures of melting and freezing using differential thermal analysis (DTA). Samples were repeatedly melted and recrystallized at selected pressures up to 15 GPa in a multianvil press. Despite the weak DTA signals due to small sample sizes at high pressure, the melting and freezing temperatures were clearly determined from the derivatives of the DTA traces. We measured a drop in the melting temperature between 9 and 12 GPa, yielding a maximum at 10(2) GPa and 245(2) ^°C. This work highlights both the successes and failures of recent theoretical models for the melting behavior of this and related systems.     

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.     

Melting curve of NaCl0.5Br0.5 up to 4.5 GPa

Abstract

The melting curve of NaCl_0.5Br_0.5 has been measured under pressure up to 4.5 GPa. The melting temperatures of Ag and NaCl have been used to determine the pressure in the sample at its melting temperature.