GaP相变及热力学性质的理论研究

利用基于密度泛函理论的第一性原理方法研究了闪锌矿和氯化钠结构的GaP的相变及热力学性质.对两种结构的能量体积曲线做公切线,得到了从闪锌矿到氯化钠结构的相变压力约为26.2GPa,与实验结果一致.通过准谐德拜模型得到了不同温度下体积和热膨胀系数与压强的关系,以及不同压强下热容与温度的关系.     

高压下CaPo弹性性质和热力学性质的第一性原理研究

利用密度泛函理论的平面波赝势方法预测研究了CaPo从岩盐结构(B1结构)到氯化铯结构(B2结构)的相变以及B1结构CaPo高压下的弹性性质以及热力学性质等.通过等焓原理发现B1→B2的相变压力为22.8GPa. 同时计算了B1结构CaPo高压下的弹性常数以及剪切模量、杨氏模量等相关弹性参数,结果发现当压力超过20GPa时,B1结构CaPo开始不稳定了,这和等焓原理所得结果相符合. 最后通过Debye模型成功获取了B1结构C 关键词： 相变 弹性性质 热力学性质 CaPo     

AlAs相变及热动力学性质的第一性原理研究

利用基于密度泛函理论的全势能线性糕模轨函法研究了闪锌矿（B3）,NiAs（B8）和岩盐（B1）结构的AlAs的相变、结构性质以及热动力学性质.对B3-B8和B3-B1结构的能量体积曲线做公切线,得到了B3→B8相变压力为5.44 GPa,并预测到B3→B1相变压力为6.46 GPa.同时计算了高压下B8相的结构性质,结果显示当V/V₀≈0.7—1.05时,c/a基本保持恒定（仅有约 0.2%的波动）;当V/V₀≈0.4—0.7,c/a随着V/V₀的减小而近似线性地增大.通过状态方程拟合,得到了AlAs的相对体积V/V₀与压强P的函数关系,B8相的状态方程与实验结果符合很好.最后利用准谐德拜模型得到了AlAs的体弹模量B随压力P的变化关系以及不同压强下热容C_V与温度T的关系. 关键词： 相变 热力学性质 第一性原理     

硝酸丙酯结构、振动频率和热力学性质的密度泛函理论研究

用密度泛函理论(DFT)研究硝酸丙酯化合物的分子结构、振动光谱和热力学等基本性质.取BLYP、B3LYP方法和6-31G*、6-31G**、6-311G*、6-311G**基组,对硝酸丙酯分子的几何构型进行全优化计算并分析其电子结构性质.和考虑了二级相关能校正的MP2/6-311G*计算结果比较表明,B3LYP/6-31G*是研究许多较大体系化合物卓有成效和颇有前途的方法.在B3LYP/6-31G*的水平上对优化后的结构进行了正则振动频率分析,用因子0.95校正后的振动光谱和实验结果比较,符合较好.进一步     

高压下单晶LiF的光学及热力学性质的密度泛函理论研究

采用平面波赝势密度泛函方法,对单晶氟化锂(LiF)在0~500 GPa静水压下的光学性质进行了理论研究,并利用Vinet状态方程和准简谐Debye模型得到了其热力学性质.理论计算结果表明单晶氟化锂(LiF)在0~500 GPa静水压范围内具有良好的透明性,吸收波段随压强的增加而出现了蓝移.计算所得晶格常数、体积模量及其对压强的一阶导数与实验值相符合.     

NbSi2奇异高压相及其热力学性质的第一性原理研究

采用基于粒子群优化算法的结构预测程序CALYPSO, 并结合第一性原理的VASP程序, 在175 GPa发现NbSi₂的奇异立方高压相. 在此结构中, Nb原子形成金刚石结构, 而Si原子则形成正四面体镶嵌在金刚石结构中. 声子谱计算结果表明该结构是动力学稳定的. 电子结构分析表明, 六角相和立方相NbSi₂均为金属, 对金属性贡献较大的是Nb原子, 而且Nb和Si原子之间存在明显的p-d杂化现象, 电荷更多地聚集在Si四面体中. 利用“应力应变”方法, 计算了NbSi₂的弹性常数, 分析了其体积模量、剪切模量、杨氏模量和德拜温度等热动力学性质随压力的变化并进行了详细的讨论. 根据剪切模量和体积模量的比值分析了NbSi₂两种相结构的脆性和延展性, 发现压力会导致六角相NbSi₂的延展性增加, 但对立方相结构的延展性影响较小; 采用经验算法计算了NbSi₂两种相结构硬度变化情况, 结合这一比值进行了详细的分析. 弹性各向异性计算结果表明, 随着压力增加, 六角结构的各向异性增强, 而立方结构的各向异性减小.     

金属铍热力学性质的理论研究

使用第一性原理方法结合平均场模型研究了压力从0到150GPa、温度从0到1500K，金属铍六角密排结构(hcp)的热力学性质，包括铍的常态性质，等温高压物态方程，以及常压下平衡体积、体弹模量随温度的变化，Hugoniot曲线等.0K物态方程由广义梯度近似下的密度泛函理论计算，粒子热运动的贡献由平均场模型计算.由于铍的Debye温度比较高，计算自由能时考虑了零点振动能修正.计算结果与已有的静力学和冲击波实验数据符合得非常好. 关键词： 热力学性质 物态方程 第一原理计算     

二氧化碳热力学性质的理论计算

用HF、MP2、QCISD和B3LYP四种方法,在6-311G* *基集合水平上对气态CO2(X1Σg)分子全优化,可直接得到热力学性质熵S和定压热容CP,在这四种方法中,B3LYP方法算得的结果与实验值最吻合,其相对误差分别只有0.00%和-0.15%;CO2分子的生成焓包含电子焓(为-393.31 kJ.mol-1)和核部分焓(为1.653 kJ.mol-1)两部分,其生成焓相对误差为0.37%.结果表明,所用理论方法计算气态CO2分子在1atm下的热力学性质是可行的.     

常压下固态KNO2热力学性质的密度泛函理论研究

采用密度泛函理论结合准谐德拜模型研究常压下300~725 K间KNO₂立方结构的热力学性质,重点分析常压下定压热容、定容热容、熵、德拜温度、体膨胀系数、平衡体积和体弹模量随温度的变化.结果显示,常压下计算的定压热容随温度的变化与实验数据符合较好,而计算的熵与实验数据相差较大.计算得到KNO₂的平均体膨胀系数约为1.837 8×10^-5K^-1,常温下(300 K)KNO₂的德拜温度约为667.13K.     

First-principles calculations for transition phase and thermodynamic properties of GaAs

The transition phase of GaAs from the zincblende (ZB) structure to the rocksalt (RS) structure is investigated by ab initio plane-wave pseudopotential density functional theory method, and the thermodynamic properties of the ZB and RS structures are obtained through the quasi-harmonic Debye model. It is found that the transition from the ZB structure to the RS structure occurs at the pressure of about 16.3\,GPa, this fact is well consistent with the experimental data and other theoretical results. The dependences of the relative volume V/V₀ on the pressure P, the Debye temperature \Th and specific heat C_V on the pressure P, as well as the specific heat C_V on the temperature T are also obtained successfully.     

First-principles study on the phase transitions,crystal stabilities and thermodynamic properties of TiN under high pressure

The structural and thermodynamic properties of titanium nitride (TiN) have been investigated by merging first-principles calculations and particle-swarm algorithm. The three phases are identified for TiN, including the B1, the $P{6}_3/mmc$, and the B2 phases. A new phase of anti-TiP structure with the space group $P{6}_3/mmc$ has been predicted. The calculated phase transition from the B1 to the $P{6}_3/mmc$ occurs at 270 GPa. The vibrational, elastic, and thermodynamic properties for the three phases have been calculated and discussed.     

First-principle study of phase stability,electronic structure and thermodynamic properties of cadmium sulfide under high pressure

By employing first principle and a quasi-harmonic Debye model, we study the phase stability, phase transition, electronic structure and thermodynamic properties of cadmium sulfide (CdS). The results indicate that CdS is a typical ionic crystal and that the zinc-blende phase in CdS is thermodynamically unstable. Moreover, the heat capacity of the wurtzite and rocksalt phases of CdS decreases with pressure and increases with temperature, obeying the rule of the Debye T³ law at low temperature and the Dulong–Petit limit at high temperature.     

Electronic and thermodynamic properties of B2-FeSi from first principles

The electronic and thermodynamic properties of B2-FeSi have been investigated using the first-principles method based on the plane-wave basis set. The calculated equilibrium lattice constant is in good agreement with available experimental and theoretical data. Our results have shown that B2-FeSi was a narrow gap semiconductor of above 0.055 eV and exhibited metallic characteristics. The density of states (DOS) can also describe orbital mixing. Using the quasi-harmonic Debye model, the thermodynamic properties of B2-FeSi have been analyzed. Variations of the Debye temperature Θ_D, thermal expansion α, heat capacity C_v, entropy S and the Grüneisen parameter γ on temperature T and pressure P were obtained successfully in the ranges of 0-2400 K and 0-140 GPa.     

Pressure-dependent mechanical and thermodynamic properties of newly discovered cubic Na2He

Recently for the first time, a stable compound of He and Na (Na₂He) is predicted at high pressure. We explore the pressure-dependent elastic, mechanical and thermodynamic properties of this newly discovered Na₂He by using ab initio technique. The calculation presents good accordance between the theoretical and experimental lattice parameters. Though the most stable structure of Na₂He is found at 300?GPa, present study ensures the mechanical stability of this compound up to 500?GPa. The study of Cauchy pressure, Pugh's ratio, and Poisson's ratio implies the ductile manner of Na₂He up to 500?GPa. According to the value of Poisson's ratio the bonding force exists in Na₂He is central. The study of Zener anisotropy factor indicates that Na₂He is an anisotropic material but near at 300?GPa approximately isotropic nature of Na₂He is revealed. The study of the bulk modulus, shear modulus, Young's modulus and Vickers hardness implies that the hardness of Na₂He can be improved by applying external pressure. However, the Debye temperature, melting temperature and minimum thermal conductivity of Na₂He are also calculated and discussed at different pressures.     

Structural and thermodynamic properties of WB at high pressure and high temperature

The structure parameters and electronic structures of tungsten boride (WB) have been investigated by using the density functional theory (DFT). Our calculating results display the bulk modulus of WB are 352±2 GPa (K′₀=4.29) and 322±3 GPa (K′₀=4.21) by LDA and GGA methods, respectively. We have analyzed the probable reason of the discrepancy from the bulk modulus between theoretical and experimental results. The compression behavior of the unit cell axes is anisotropic, with the c-axis being more compressible than the a-axis. By analyzing the bond lengths information, it also demonstrated that WB has a lower compressibility at high pressure. From the partial densities of states (PDOS) of WB, we found that the Fermi lever is mostly contributed by the d states of W atom and p states of B atom and that the contributions from the s, p states of W atom and s states of B atom are small. Moreover, using the Gibbs 2 program, the thermodynamic properties of WB are obtained in a wide temperature range at high pressure for the first time in this work.     

Structural and thermodynamic properties of Os from first-principles calculations

We investigate the structural, thermodynamic and electronic properties of Os by plane-wave pseudopotential density functional theory method. The obtained lattice constants, bulk modulus and cell volumes per formula unit are well consistent with the available experimental data. Especially, from our calculated bulk modulus, we conclude that Os is more compressible than diamond. Moreover, the temperature induced phase transition of Os from HCP structure to FCC structure has been obtained. It is found that the transition temperature of Os at zero pressure is 2702 K. However no transition pressure is found in our calculations. The effect of bulk modulus B as well as other thermodynamic properties of Os (including the thermal expansion α and the Grüneisen constant γ) on temperatures have also been studied. Our calculated thermal expansion α=1.510×10⁻⁵ K⁻¹ and the Grüneisen constant γ=2.227 for HCP structure at room temperature agree very well with the experimental data. The density of states for HCP structure at 0 K and FCC structure at transition temperature 2702 K are also investigated in our work.     

基于FORTRAN和VB的常用制冷剂热物性计算软件的开发

采用状态方程和多项式函数拟合两种方法,利用现有的物性计算,提出一些常用制冷剂的热力性质计算公式。利用FORTRAN和VB6.0建立了包含R12、R22、R134a的热物性计算软件。程序开发中采用层次模块的思想和开放式结构,便于增加新的制冷剂,具有良好的可扩充性。