MgB2状态方程的第一原理计算

利用基于局域密度近似框架下的第一原理平面波方法, 结合HGH型相对论分析赝势, 对MgB2的晶格参数和状态方程进行了计算,所得结果与实验值及其他理论得到的计算值相符合.     

SrS相变和弹性性质的第一原理计算

利用平面波赝势密度泛函方法和准谐德拜模型研究了SrS从NaCl结构到CsCl结构的相变以及弹性性质.在零温下,我们计算的相变压强为17.9 GPa,这与实验值和其他作者的计算值符合很好.研究还表明:相变压强随温度增加而非线性地增加,然而力学不稳性的压强随温度增加而线性地增加.     

超导体MgB2弹性常数的第一原理计算

利用全势线性muffin-tin轨道(FP-LMTO)方法, 结合在密度泛函理论框架下的广义梯度近似, 研究了六角密堆积结构超导体MgB2的晶格参数, 弹性常数, 以及体积模量及其对压强的微分. 计算结果显示当晶格参数c和a的比率c/a大约为1.138时, MgB2的结构最稳定.本文所得到的计算结果与实验值及其他作者利用不同方法得到的计算值相符合.     

纤维锌矿结构GaN热力学性质的第一性原理计算(英文)

利用平面波赝势密度泛函方法,结合广义梯度近似,对纤维锌矿GaN的结构和热力学性质进行了计算,发现纤维锌矿结构GaN的最稳定的结构对应的结构参数为:u=0.3769,c/a=1.628,a=3.204,andc=5.216,所得结果与实验和其他理论的结果相符.通过准谐德拜模型,我们还成功得到了w-GaN的相对体积、热膨胀系数、热容和德拜温度等热力学量.     

Y（Fe，Si）12化合物结构与磁性的第一原理计算

R（Fe，Si）12（R=Y，Nd）型稀土永磁材料具用重要的实际应用价值，然而对于V（Fe，Si）12化合物微观机理的理论计算工作至今尚未开展。采用新近发展的全势能线性缀加平面波（（L）APW）＋局域轨道（10）和广义梯度近似（GGA）密度泛函方法对其结构与磁性进行了研究。     

铀高压状态方程的第一原理研究北大核心CSCD

利用基于密度泛函理论的赝势平面波方法,研究了铀的状态方程。基态计算结果表明:在100GPa的压强范围内,α-U结构是最稳定的相。利用量子分子动力学方法(QMD),计算了有限温度下的状态方程,在此基础上计算了冲击Hugoniot点,并与实验数据及其他理论计算结果进行了比较,结果表明:在100GPa以下的压力区域内,QMD计算结果与实验结果符合得很好,而用快速计算状态方程(QEOS)方法计算的结果偏软。     

氘、氦及其混合物状态方程第一原理研究

利用量子分子动力学方法对氘、氦及其混合物的状态方程进行了研究. 计算了氦在密度0.32-5 g/cm³, 温度1000-50000 K范围内的状态方程, 并与化学模型的结果进行了比较; 同时计算了冲击Hugoniot曲线, 与气炮实验的结果符合得很好. 通过计算对分布函数及态密度, 探讨了氦在高温高压下发生金属-绝缘体转变的机理. 计算了氘在密度0.19-0.84 g/cm³, 温度20-50000 K范围内的状态方程; 并计算了理论Hugoniot状态, 由于没有考虑原子的零点运动, 在低温区, 理论结果比实验值小. 对氘氦混合物的状态方程进行了研究, 计算了温度和密度区间为100-50000 K, 0.19-0.84 g/cm³, 不同混合度下的293个状态点的状态方程. 对线性混合近似进行了考查, 结果表明线性混合近似是一个粗略的近似.     

第一原理研究界面弛豫对InAs/GaSb超晶格界面结构、能带结构和光学性质的影响

通过广义梯度近似的第一原理全电子相对论计算, 研究了不同界面类型InAs/GaSb超晶格的界面结构、电子和光吸收特性. 由于四原子界面的复杂性和低对称性, 通过对InAs/GaSb超晶格进行电子总能量和应力最小化来确定弛豫界面的结构参数. 计算了InSb, GaAs型界面和非特殊界面(二者交替)超晶格的能带结构和光吸收谱, 考察了超晶格界面层原子发生弛豫的影响.为了证实能带结构的计算结果, 用局域密度近似和Hartree-Fock泛函的平面波方法进行了计算. 对不同界面类型InAs/GaSb超晶格的能带结构计算结果进行了比较, 发现界面Sb原子的化学键和离子性对InAs/GaSb超晶格的界面结构、 能带结构和光学特性起着至关重要的作用.     

Ti,Nb, Al及其二元合金状态方程的第一原理计算及其应用

金属状态方程对于探究金属及合金原子的相互作用中起到了至关重要的作用.本文使用第一原理计算了Ti, Nb, Al及其二元合金能量与原子间距关系(E-r关系),并使用得到的E-r关系计算了金属及合金的体弹性模量,结果与实验值吻合的很好.计算结果表明,使用不同的赝势,计算不同金属表现出不同的适用性;二元合金的E-r曲线处于对应合金化元素曲线之间;并且合金的E-r曲线会更靠近合金内含量较高的元素的E-r曲线;并发现合金E-r关系可通过组成合金的纯金属的E-r关系计算获得,且用该方法计算材料的体弹性模量与实验值非常符合.     

混合物状态方程的计算

多介质流体动力学过程的数值模拟往往涉及混合物状态方程的计算. 做图法和Newton 法是混合物状态方程计算常采用的方法, 前者虽直观精度却差, 后者计算效率高却只具有局部收敛性, 当解与其初始猜测值相差较远时Newton法不一定能够获得收敛解. 为此, 本文给出一种具有大范围收敛性的嵌入算法(imbedding method)求解混合物状态方程, 其基本思想是通过引入嵌入参数, 将待解的混合物状态方程和易解的混合物状态方程线性组合, 构成嵌入方程组, 当嵌入参数从0连续地变化到1 时, 嵌入方程组的解由易解的混合物状态方程的解连续地变化为待解的混合物状态方程的解. 嵌入方程组可由Newton法迭代求解, 也可转化为以嵌入参数为自变量的常微分方程组, 从而易于由成熟的计算方法如梯形法等进行求解. 进一步利用热力学基本关系, Maxwell形式的微分方程描述了压力和温度随嵌入参数的演化速率与应变速率和组分质量分数演化速率的关系. 对铅锡混合物热力学量的计算表明了本文算法的有效性.     

Thermodynamic properties of MgO under high pressure from first-principles calculations

The equations of state (EOS) and other thermodynamic properties of the rocksalt (RS) structure MgO are investigated by ab initio plane-wave pseudopotential density functional theory method. The obtained results are consistent with the experimental data and those calculated by others. Through the quasi-harmonic Debye model, in which the phononic effects are considered, the dependences of relative volume V/V₀ on pressure P, cell volume V on temperature T, and Debye temperature Θ and specific heat C_V on pressure P are successfully obtained. The variation of the thermal expansion with temperature and pressure is investigated, which shows the temperature has hardly any effect on the thermal expansion at higher pressure.     

First-principles calculations of thermodynamic properties of TiB2 at high pressure

The equations of state (EOS) and other thermodynamic properties of TiB₂ are investigated using ab initio plane-wave pseudopotential density functional theory method. The obtained results are in good agreement with the experimental measured data and other theoretical calculated ones. Through the quasi-harmonic Debye model, in which the phononic effects are considered, the dependences of relative volume V/V₀ on pressure P, cell volume V on temperature T, and Debye temperature Θ and specific heat C_V on pressure P are successfully obtained.     

First-principles calculations for structure and equation of state of MgB2 at high pressure

We investigate the structure and the equation of state of compound MgB₂ at high pressure using the full-potential linearized muffin-tin orbital scheme within the generalized gradient approximation correction in the frame of density functional theory. Through the quasi-harmonic Debye model, in which the phononic effects are considered, we have obtained successfully the bulk modulus and the thermal expansion of MgB₂.     

First-principles calculations of structure and high pressure phase transition in gallium nitride

The phase transitions of semiconductor GaN from the Wurtzite (WZ) structure and the zinc-blende (ZB) structure to the rocksalt (RS) structure are investigated by using the first-principles plane-wave pseudopotential density functional method combined with the ultrasoft pseudopotential scheme in the generalized gradient approximation (GGA) correction. It is found that the phase transitions from the WZ structure and the ZB structure to the RS structure occur at pressures of 46.1 GPa and 45.2 GPa, respectively. The lattice parameters, bulk moduli and their pressure derivatives of these structures of GaN are also calculated. Our results are consistent with available experimental and other theoretical results. The dependence of the normalized formula-unit volume $V/V_{0 }$ on pressure $P$ is also successfully obtained.     

Analysis of K-prime equations of state

In the present communication, we introduce a new K-prime equation of state (EoS), which can be used to understand the interior of the Earth. The newly developed EoS is found to yield similar results as given by Stacey K-prime EoS and Keane EoS, in reference to the seismological data. However, the zero pressure and infinite pressure extrapolation of higher pressure derivatives of bulk modulus are found different for different K-prime EoS.     

Three equations of state for solids considering thermal effect

 The Einstein model to consider thermal effect in universal equations of state (UEOS) is modified. It is proposed that the zero-point vibration term should be deleted in a thermal UEOS, and the parameters cannot be directly taken as experimental data at a reference temperature, V_R, B_R, and , but their values at absolute zero temperature, V₀, B₀, and . An approach is proposed to solve V₀, B₀, and from V_R, B_R, and . The approaches are applied to three typical universal EOSs, including the Baonza, mGLJ and Morse EOSs. The numerical results show that the solved values of parameters are almost identical for different EOSs. And the thermo-physical properties predicted through different EOSs are almost identical at zero- and low-pressure conditions, once the same approach and input experimental data are used to solve the parameters. It is concluded that the prediction of thermo-physical properties at zero- and low-pressure conditions cannot be taken as the criteria to judge the applicability of a universal EOS.     

A detonation stability formulation for arbitrary equations of state and multi-step reaction mechanisms

A general normal-mode linear stability formulation of steady planar detonation waves is presented that is valid both for an arbitrary equation of state and for multi-step, multi-species chemical kinetics. The general formulation can be used for many purposes, including an examination of gaseous detonation stability with complex reaction kinetics in which the individual reacting species have variable thermochemical properties. In the present paper, we consider two cases that could not be obtained by previous one-step chemistry, polytropic gas formulations: the first concerns the effect of a difference in heat capacities between product and fuel species, as well as a possible mole change, in a single-step irreversible reaction. The second examines the effects of exothermic or endothermic heat release/absorption in the chain-initiation stage of a model three-step reaction.