第一性原理对金的高压相变和零温物态方程的计算

在密度泛函理论下，用缀加平面波加局域轨道方法，分别采用广义梯度近似(GGA)和局域密度近似(LDA)对金的面心立方晶格结构(fcc)、体心立方晶格结构(bcc)和六角密堆积结构(hcp)的结构能量进行了计算.在GGA下，计算得出fcc向hcp和hcp向bcc的相变分别发生在380GPa 和1250GPa；而LDA下相变分别发生在490GPa和790GPa.当计算压强达到2TPa时，bcc在这两种近似下仍然保持稳定的结构.根据不同体积下不同结构的电子态密度的特征，对发生相变的物理原因进行了定性的分析，在此基础上得到了金的零温状态方程. 关键词： 缀加平面波方法 固态相变 电子态密度 物态方程     

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

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

高压下TATB压缩性质的LDA和GGA比较研究

本文采用第一性原理密度泛函理论,结合平面波赝势方法,采用局域密度近似(LDA)和广义梯度近似(GGA)两种方法计算了TATB晶体在高压(0~7 GPa)压缩下的结构和物理性质,并与实验数据进行了比较.详细讨论了TATB晶体的晶体结构和分子构型随压力的变化.     

从头算方法研究面心立方铝在高温高压下的热力学状态方程

用从头算方法优化计算了面心立方铝的电子结构和总能,得到了它在零温下的状态方程和弹性性质.将得到的总能与晶格体积拟合到Debye模型,获得了非平衡态下的Gibbs自由能与温度、压力之间的关系,在此基础上计算了相应的热状态方程,利用Burakovsky-Preston-Silbar （BPS） 熔化模型计算了铝的熔化曲线.所有的电子结构和总能计算都是基于局域密度近似（LDA）和广义梯度近似（GGA）的平均得到的.计算得到的铝在高温、高压下的状态方程与一些热力学性质和熔化曲线同冲击波和静高压实验数据在225 G 关键词： 铝 热力学状态方程 从头算 熔化曲线     

金红石→萤石TiO2相变的第一性原理计算

利用第一性原理平面波密度泛函理论并结合准谐德拜模型,对TiO₂从金红石结构到萤石结构的相变进行理论研究.依据焓相等原理,得到TiO₂从金红石结构到萤石结构的相变压强为45.32 GPa,并通过吉布斯能的计算发现相变压强随温度的增加而增加,运用拟合的方法估计实验条件下的相变压强(温度为1900-2100K时,P_t=47.604~47.756 GPa).运用等吉布斯能原理,得到零压下从金红石到萤石结构TiO₂的相变温度为2029 K.计算结果与实验及其他理论值相符,并通过准谐德拜模型成功获得相变前后的结构和热力学性质,如相对体积、热膨胀系数、热容和德拜温度等.     

Elastic properties and phonon dispersions of rhenium in hexagonal-close-packed structure under pressure from first principles

The elastic properties of the hexagonal-close-packed (hcp) structure rhenium (Re) and their behavior under pressure are investigated using the local density approximation (LDA) and the generalized gradient approximation (GGA). The obtained high pressure elastic constants are well consistent with previous theoretical date, while large discrepancies are found between theory and the high pressure experiments. The calculated isothermal bulk modulus B₀ (376 GPa for GGA and 389 GPa for LDA) and its initial pressure derivative (4.52 for LDA and 4.58 for GGA) compare favorably with the experimental values. Moreover, it is found that the value of c/a, B/G, Poisson's ratio, and B_c/B_a are virtually independent of pressure. We also performed calculation for phonon dispersions at high pressure. GGA in our calculation exhibits a same trend as the high pressure experimental curve.     

LiH的热力学性质及分子内部运动对体系热力学性质的影响

利用基于密度泛函理论的第一性原理平面波超软赝势方法,计算了LiH在零温零压下的晶格常数、体弹模量,计算结果与实验值和其他理论计算值符合得较好.通过准谐德拜模型计算了LiH在压强为0-80GPa、温度为0-2000K范围内,体积膨胀率、热涨系数、德拜温度及定容热容随压强和温度的变化关系.最后,以代数方法（AM）和势能变分法(PVM)为基础,运用统计热力学理论计算了LiH分子内部运动对体系热力学性质的影响.     

Structural phase stability,elastic parameters and thermal properties of YN from first-principles calculation

The structural phase stability, elastic parameters and thermodynamic properties of YN at normal and under high pressure are reported. The calculations are mainly performed using the full-potential linearized augmented plane wave method within the density functional theory. Both local density approximation (LDA) and generalized gradient approximation (GGA) are used to model the correlation-exchange potential. The calculated equilibrium lattice parameter and the bulk modulus show good accordance with the experimental and previous theoretical reports. The phase transition from the NaCl (B1) structure to the CsCl (B2) structure is found to occur at 131?GPa within GGA and 115?GPa within LDA. The linear pressure coefficients of the different elastic moduli being addressed here are also determined along with the mechanical and dynamical stability criteria which are shown to be satisfied for YN with B1 phase under normal conditions. Besides, the heat capacity and other thermodynamic parameters are examined and discussed versus temperature.     

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.     

氧化镍（绿镍矿）等温压缩及高压相变研究

 本文采用DAC（金刚石压砧高压腔）装置，对氧化镍进行了静水压、非静水压、电导率测量等系统高压实验，获取了氧化镍等温压缩、高压相变及电导率压力效应的新结果，并在实验数据的基础上，对其高压相变与电性及磁性变化关系及体弹性模量作了分析讨论。     

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

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

Melting line of hydrogen at high pressures

The insulator to metal transition in solid hydrogen was predicted over 70 years ago but the demonstration of this transition remains a scientific challenge. In this regard, a peak in the temperature versus pressure melting line of hydrogen may be a possible precursor for metallization. However, previous measurements of the fusion curve of hydrogen have been limited in pressure and temperature by diffusion of hydrogen into the gasket or diamonds. To overcome this limitation we have used an innovative technique of pulsed laser heating of the sample and find a peak in the melting line at P=64.7+/-4 GPa and T=1055+/-20 K.     

Diamond under pressure: Ab-initio calculations of the equation of state and optical phonon frequency revisited

The pressure–volume relationship and the zone-center optical phonon frequency of cubic diamond at pressures up to 600?GPa have been calculated based on density functional theory (DFT). Different methods (ranging from pseudopotential to all-electron approaches) and different approximations [local density approximation (LDA), generalized gradient approximation (GGA)] for the exchange-correlation (XC) energy have been applied. The results are compared to recent precise measurements up to 140?GPa. Possible implications for the experimental pressure determination based on the 1986 calibration of the ruby luminescence method are discussed.     

First-principles investigation of the structural,optoelectronic and thermodynamic properties of InAsxP1-x ternary alloys under hydrostatic pressure effect

The main objective of our work is the study of structural, optoelectronic and thermodynamic properties of InAs_xP_1-x alloys in the zinc-blende structure using the full potential linearized augmented plane wave method (FP-LAPW) based on density functional theory (DFT). Different exchange correlation potentials were used, as well as the local density approximation (LDA) and the generalized gradient approximation (GGA) parameterized by Perdew–Burke–Ernzerhof (PBE-GGA) and PBE sol-GGA of Perdew, to estimate structural properties such as lattice parameters, the bulk modulus and its first pressure derivative. For electronic properties, the Tran-Blaha modified Becke–Johnson potential (TB-mBJ) was used for density of states (DOS) and band structure calculations. The results show that the compounds of interest are semiconductors with direct band gaps for the full range of x compositions and that the optical band gap decreases from 1.58 to 0.41 eV with increasing As concentrations. The obtained results show a good agreement with experimental and theoretical data found in the literature. In addition, we have investigated the dielectric function as well as the refractive index and the reflectivity. The electronic and optical properties were studied under hydrostatic pressure (P = 0, 5, 10, 15, 20, and 25 GPa), and it was found that the band gaps of the binary compounds change from a direct to an indirect harmonic Debye model was used, which takes into account the effect of pressure P and temperature T on the lattice parameter, to explore the heat capacity, the Debye temperature and the entropy under pressures ranging from 0 to 20 GPa and temperatures ranging from 0 to 1200 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.     

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.     

In situ X-Ray study of thermal expansion and phase transition of iron at multimegabar pressure

The density of varepsilon-iron has been calculated at pressures and temperatures up to 300 GPa and 1300 K, respectively. We observe varepsilon to beta phase transition at pressures between 135 and 300 GPa and temperature above 1350 K; the pattern can be interpreted in terms of double hexagonal close-packed structure. The density calculated at high pressure and temperature (330-360 GPa and 5000-7000 K) closely matches with preliminary reference Earth model density, thereby imposing constraint on the composition of the Earth's inner core.     

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.     

常温及常压至1.2GPa条件下异辛烷的拉曼光谱研究

利用碳化硅压腔在室温(25℃)下,研究了异辛烷(2,2,4-三甲基戊烷)在常压至1.2GPa条件下的拉曼光谱特征。研究结果表明,异辛烷CH2和CH3的碳氢伸缩振动的拉曼位移随着压力的增大均呈线性向高频方向移动,其拉曼位移与压力的函数关系为:ν2 873=0.002 8P+2 873.3;ν2 905=0.004 8P+2 905.4;ν2 935=0.002 7P+2 935.0;ν2 960=0.012P+2 960.9。在1.0GPa附近,异辛烷的拉曼位移出现突变,与显微镜下观察发生的异辛烷液-固相变一致。结合异辛烷在常压下的熔点数据,获得了异辛烷的液-固两相相图,并根据克拉贝龙方程获得了液-固相转变过程中的摩尔体积变化量ΔVm=4.46×10-6 m3.mol-1和熵变ΔS=-30.32J.K-1.mol-1。     

Thermodynamic properties of OsB under high temperature and high pressure

The energy-volume curves of OsB have been obtained using the first-principles plane-wave ultrasoft-pseudopotential density functional theory (DFT) within the generalized gradient approximation (GGA) and local density approximation (LDA). Using the quasi-harmonic Debye model we first analyze the specific heat, the coefficients of thermal expansion as well as the thermodynamic Grüneisen parameter of OsB in a wide temperature range at high pressure. At temperature 300 K, the coefficients of thermal expansion α_V by LDA and GGA calculations are 1.67×10⁻⁵ 1/K and 2.01×10⁻⁵ 1/K, respectively. The specific heat of OsB at constant pressure (volume) is also calculated. Meanwhile, we find that the Debye temperature of OsB increases monotonically with increasing pressure. The present study leads to a better understanding of how the OsB materials respond to pressure and temperature.