BeO高压相变和声子谱的第一性原理计算

采用第一性原理方法计算了BeO在零温时的高压相变和三种结构在零温零压时的声子谱.相变的计算表明,在122GPa左右的压力下BeO会发生从纤锌矿(B4)结构到氯化钠(B1)结构的相变,而闪锌矿(B3)结构在零温零压下是一种可能的亚稳态结构.采用冷声子方法计算了这三种结构的BeO在零温零压下的声子谱.计算结果表明:B1结构在零温零压下是一种不稳定的结构;尽管B4结构和B3结构具有明显的相似性,仍然可以通过声子谱来很好的区分.最后根据准简谐近似理论计算得到了BeO的高温高压相图.     

高压下金属Ba的结构稳定性以及热动力学的第一原理研究

采用第一原理方法计算了高压下金属Ba的三个高压相 Ba-I, Ba-Ⅱ和Ba-V的稳定性及热动力学性质.结果表明, Ba的三个高压相在0 K时在其压力范围内都是动力学和力学稳定的;但随压力增加, Ba-I 和Ba-Ⅱ 的声子谱频率出现异常"软化",而Ba-V则出现"硬化".虽然 Ba-Ⅱ 和 Ba-V 同为六方密堆(hcp)结构,计算表明它们在高压下表现出了不同的弹性各向异性.计算同时发现 Ba-Ⅱ 在更高的压力下仍满足力学稳定条件,但声子谱有虚频存在, 表明动力学失稳是Ba-Ⅱ在压力下向Ba-I!V相转变的原因. 计算和比较了同为六方密堆(hcp)结构的Ba-Ⅱ和Ba-V在高压下的声速、 德拜温度、体模量、剪切模量等力学和热学性质, 展现了金属Ba在压力下的稳定机制和热动力学性质.     

高压下c-BN的力学、电子结构及光学性质研究

为了研究不同压力下c-BN的力学性质、电子结构以及光学性质的变化,基于密度泛函理论构建了不同压力下c-BN的晶体模型。发现c-BN在不同压力下均力学稳定,随着压力增加,c-BN的弹性常数逐渐增大,材料的可压缩性逐渐变差;c-BN在零压下的禁带宽度为4.367 eV,表明c-BN是间接宽带隙半导体,随着压力增加,c-BN的禁带宽度逐渐增大,态密度谱图变化不明显;对Milliken 布居分布在不同压力下进行分析,表明随着压力增加,B、N原子杂化后形成的B-N共价性增强。对c-BN的复介电函数、折射率、反射率等进行分析,发现随着压力增大,它们都产生一定蓝移,且在整个可见光谱范围以及红外与紫外（约从204nm开始）光谱的很大范围内都透明。研究结果对高压下c-BN的应用有一定参考价值。     

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

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

P掺杂YVO_4电子结构和弹性性能的第一性原理研究

基于密度泛函理论(DFT)的第一性原理平面波超软赝势平面波方法,在零压的条件下,建立了未掺杂以及P原子取代V原子的YVO_4与YV_(0.75)P_(0.25)O_4模型,然后对模型进行了电子结构和弹性性能的计算.计算结果表明:锆英石结构YVO_4与YV_(0.75)P_(0.25)O_4在零压下是稳定的;YV_(0.75)P_(0.25)O_4属于延展性材料,具有大的弹性各项异性特征,YV_(0.75)P_(0.25)O_4的体积模量、剪切模量、杨氏模量分别为136.61 GPa、62.14 GPa、161.87 GPa,泊松比为0.303,研究结果表明P掺杂能减少YVO_4的体积模量,并且增大其剪切模量与杨氏模量.     

过渡金属Tc及其氮化物TcN,TcN2,TcN3与TcN4低压缩性的第一性原理计算研究

采用基于密度泛函理论的第一性原理计算方法研究了Tc及其氮化物的弹性性质、电子结构、原子布局数等,并分析和计算了它的氮化物的理论硬度.结果表明随着Tc中氮的掺入量的渐增,出现层状结构,它们的弹性模量并非单调增加,其中TcN的体弹模量最大而剪切模量最小;层状结构的TcN3的剪切模量最大而体弹模量最小,TcN4的理论计算硬度最大.说明在Tc的氮化物中,其体弹模量与材料的晶体结构、平均每个原子上的电荷密度和材料的质量密度有关;化学键的共价性结构和氮元素的含量对理论计算硬度有正作用;而剪切模量的极大值则与其层状结构及体系中一定量的方向基本一致的N-N键相关.     

AsI3电子结构与弹性性质的第一性原理研究

采用基于密度泛函理论的第一性原理赝势方法,对AsI3的平衡态晶格常数、弹性常数和电子结构进行了研究。研究结果表明,R-3结构的AsI3在零压下是稳定的,优化得到的平衡结构参数与实验值符合的很好。AsI3是脆性材料,具有大的弹性各向异性特征。R-3相AsI3的块体模量、剪切模量和杨氏模量分别为14.2 GPa, 9.8 GPa和23.9 GPa,泊松比为0.22,德拜温度是163 K。能带结构计算表明,AsI3是带隙为2.34 eV的间接带隙半导体。AsI3的化学键是弱共价键和离子键的混合。     

钍基碳化物ThC的高压结构稳定性研究

钍基碳化物ThC是重要的核材料,基于第一性原理计算,我们研究了ThC的结构稳定性及相变.声子色散曲线表明基态结构的晶格能够稳定到42 GPa附近,当压力继续升高时虚频的出现表明了动力学的不稳定.经比较基态Fm-3m及高压相Pnma结构的自由能,发现零温下,相变压力为22.4 GPa,符合声子色散曲线的结果,同时预测了两种结构的相变边界.     

利用第一性原理方法研究CuXSe_2(X=B,Al,Ga,In,Tl)力学性质

在广义梯度近似(GGA)和局域密度近似(LDA)下,采用第一性原理方法研究CuXSe2(X=B,Al,Ga,In,Tl)晶体结构的稳定性和力学性质.分析CuXSe2(X=B,Al,Ga,In,Tl)结构的晶格常数,弹性常数,体积模量,剪切模量,杨氏模量,泊松比,对比各材料的力学性质变换规律.计算结果表明,根据力学稳定判据,在零温零压下,CuXSe2(X=B,Al,Ga,In,Tl)的晶体结构是力学稳定的.经过同主族替换,发现晶格常数越大,体弹性模量就越小.这可以解释为X离子半径逐渐增大,晶格常数逐渐增加,晶体的可压缩性也增加.另外,CuTlSe2的剪切模量最小,不容易发生剪切形变.CuBSe2的杨氏模量最大,其刚度最高.由Pugh经验关系可知CuXSe2(X=B,Al,Ga,In,Tl)均属于韧性材料.     

利用第一性原理方法研究CuXSe (X=B，Al，Ga，In，Tl)力学性质

在广义梯度近似（GGA）和局域密度近似（LDA）下,采用第一性原理方法研究CuXSe2 (X=B,Al,Ga,In,Tl)晶体结构的稳定性和力学性质。分析CuXSe2 (X=B,Al,Ga,In,Tl)结构的晶格常数,弹性常数,体积模量,剪切模量,杨氏模量,泊松比,对比各材料的力学性质变换规律。计算结果表明,根据力学稳定判据,在零温零压下,CuXSe2 (X=B,Al,Ga,In,Tl)的晶体结构是力学稳定的。经过同主族替换,发现晶格常数越大, 体弹性模量就越小。 这可以解释为X 离子半径逐渐增大, 晶格常数逐渐增加, 晶体的可压缩性也增加。另外,CuTlSe2的剪切模量最小,不容易发生剪切形变。CuBSe2的杨氏模量最大,其刚度最高。由Pugh经验关系可知CuXSe2 ( X=B,Al,Ga,In,Tl)均属于韧性材料。     

First principle study of structural,electronic, mechanical,dynamic and optical properties of half-Heusler compound LiScSi under pressure

In this study, several physical properties of LiScSi compound with MgAgAs phase were investigated via the plane-wave pseudo-potential technique in density functional theory (DFT). The calculated total energy-atomic volume was fitted to the Murnaghan equation of state in order to obtain bulk modulus, their first derivatives and the lattice constant. These results were compared to findings of recent literature. Afterwards, the partial density of states (PDOS) and charge density differences were used to evaluate the electronic band structure of LiScSi under pressure. By analysing elastic properties (shear modulus, Poisson ratio, Young’s modulus, etc.) of the material, it has been shown that MgAgAs phase of the compound is mechanically stable under pressure. Moreover, the dynamical stability of this compound is calculated by means of the phonon dispersion curves and one-phonon DOS. Finally, the optical properties and related parameters (refractive index, dielectric function, and loss function) of LiScSi were examined with subject to different pressures.     

First principles study of structural,phonon, optical,elastic and electronic properties of Y3Al5O12

Structural, phonon, optical, elastic and electronic properties of Y₃Al₅O₁₂ have been investigated by means of the first principles method with the Cambridge Serial Total Energy Package (CASTEP) code based on the density functional theory. The calculated lattice parameters, valence charge density, bond length and single crystal elastic properties at zero pressure are in good agreement with the available experimental data. The close agreement with the experimental values provides a good confirmation of the reliability of the calculations. Optical, elastic and phonon properties of Y₃Al₅O₁₂ under pressures are performed. The results that are obtained show the changes of optical and elastic properties under the influence of applied pressure, and proving the dynamical stability of YAG are destructed when applied pressure up to 7 GPa. Moreover, polycrystalline elastic moduli are deduced according to the Reuss assumption. Those elastic constants provide important parameters that describe reliability of both physical model and engineering application at the atomistic level. The result of the density of states explains the nature of the electronic band structure.     

高温高压下立方氮化硼和六方氮化硼的结构、力学、热力学、电学以及光学性质的第一性原理研究

本文采用基于密度泛函理论的第一性原理平面波赝势和局域密度近似方法,优化了立方和六方氮化硼的几何结构,系统地研究了零温高压下立方和六方氮化硼的几何结构、力学、电学以及光学性质.结构与力学性质研究表明:立方氮化硼的结构更加稳定,两种结构的氮化硼均表现出一定的脆性,而六方氮化硼的热稳定性则相对较差;电学性质研究表明:立方氮化硼和六方氮化硼均为间接带隙半导体,且立方氮化硼比六方氮化硼局域性更强;光学性质结果显示:立方氮化硼和六方氮化硼对入射光的通过性都很好,在高能区立方氮化硼对入射光的表现更加敏感.此外,还研究了高温高压下立方氮化硼的热力学性质,并得到其热膨胀系数、热容、德拜温度和格林艾森系数随温度和压力的变化关系.本文的理论研究阐述了高压下立方氮化硼和六方氮化硼的相关性质,为今后的实验研究提供了比较可靠的理论依据.     

Prediction of Pressure-Induced Structural Transition and Mechanical Properties of Mg Y from First-Principles Calculations

Using the particle swarm optimization algorithm on crystal structure prediction,we first predict that Mg Y alloy undergoes a first-order phase transition from Cs Cl phase to P4/NMM phase at about 55 GPa with a small volume collapse of 2.63%.The dynamical stability of P4/NMM phase at 55 GPa is evaluated by the phonon spectrum calculation and the electronic structure is discussed.The elastic constants are calculated,after which the bulk moduli,shear moduli,Young's modui,and Debye temperature are derived.The brittleness/ductile behavior,and anisotropy of two phases under pressure are discussed in details.Our results show that external pressure can change the brittle behavior to ductile at10 GPa for Cs Cl phase and improve the ductility of Mg Y alloy.As pressure increases,the elastic anisotropy in shear of Cs Cl phase decreases,while that of P4/NMM phase remains nearly constant.The elastic anisotropic constructions of the directional dependences of reciprocals of bulk modulus and Young's modulus are also calculated and discussed.     

高温高压下碳化钨晶体的结构、力学、电子、光学以及热力学性能的第一性原理计算

本文利用密度泛函理论中的广义梯度近似对碳化钨晶体的三种结构(碳化钨相、闪锌矿相以及纤锌矿相)进行了优化,得到能量最低的稳定构型,并在此基础上计算了它的力学、电子、光学和高温高压下的热力学性质.研究表明:在0～300 GPa压力范围内,碳化钨相具有最高的稳定性.同时,高压下碳化钨相的弹性常数满足Born-Huang准则,且0 GPa和300 GPa下的声子色散没有虚频,证明了高压下碳化钨相的静力学稳定性和动力学稳定性.电子性质表明了碳化钨的金属性.光学性质表明碳化钨在高能区很难吸收光.热力学性质的研究表明:体积比V/V_0对压强的变化更敏感;高温时C_V曲线近似一条直线;给定压强下热膨胀系数α在600 K温度以上增长非常缓慢;压强对德拜温度Θ_D的影响较大;在低压下格林艾森系数γ的变化较大.     

First-principles investigations on structural,elastic, electronic properties and Debye temperature of orthorhombic Ni3Ta under pressure

The structural, elastic, electronic properties and Debye temperature of Ni₃Ta under different pressures are investigated using the first-principles method based on density functional theory. Our calculated equilibrium lattice parameters at 0 GPa well agree with the experimental and previous theoretical results. The calculated negative formation enthalpies and elastic constants both indicate that Ni₃Ta is stable under different pressures. The bulk modulus B, shear modulus G, Young’s modulus E and Poisson’s ratio ν are calculated by the Voigt–Reuss–Hill method. The bigger ratio of B/G indicates Ni₃Ta is ductile and the pressure can improve the ductility of Ni₃Ta. In addition, the results of density of states and the charge density difference show that the stability of Ni₃Ta is improved by the increasing pressure. The Debye temperature Θ _D calculated from elastic modulus increases along with the pressure.     

新型石墨插层化合物HfC2结构与性质的第一性原理研究

在高压下,预测一种新型石墨插层化合物HfC₂.采用第一性原理方法对其在0 GPa下的结构和性质进行研究,分别采用GGA-PBESOL、GGA-PW91和LDA方法进行结构优化,得到的晶体学数据基本相同.弹性常数和声子谱计算证实其力学和晶格动力学稳定性,表明HfC₂在0 GPa下能够稳定存在.采用GGA-PBESOL方法计算得到HfC₂的体模量和剪切模量达到265 GPa和118 GPa,Pugh比k<0.57,是一种具有高体模量的韧性材料.HfC₂存在C-C、Hf-C共价作用,且具有金属特性和特殊层状结构,是其具有高体模量和韧性的原因.最后,对HfC₂在0~500 GPa内的键长、体模量、剪切模量、k值等进行研究,探索其力学性质随压力变化的规律.     

Electronic,elastic and dynamical properties of MgSe under pressure: rocksalt and iron silicide phase

The electronic, elastic and dynamical properties of MgSe in the rocksalt (B1) and iron silicide (B28) phase and the effects of pressure on these properties are investigated using first-principles method. The calculated electronic band structure indicates that the B1 phase of MgSe presents an indirect band-gap feature and the band gaps initially increase with pressure and subsequently decrease upon compression. Remarkably, an indirect-to-direct band-gap transition has been observed at the phase transition pressure. The elastic constants, bulk modulus, shear modulus, Young’s modulus, elastic anisotropy and B/G ratio of MgSe in the B1 and B28 phase at high pressure have also been investigated. The bulk modulus, shear modulus and Young’s modulus all increase monotonously with the increasing of pressure for the B1 and B28 phase of MgSe. The calculated phonon frequencies of the B1 phase at zero pressure agree well with available theoretical results. And the transverse acoustic phonon TA(X) mode of this phase completely softening to zero at 82 GPa. The phonon curves of the B28 phase under pressure have also been successfully investigated.     

Lattice stabilities,mechanical and thermodynamic properties of Al_3Tm and Al_3Lu intermetallics under high pressure from first-principles calculations

The effects of high pressure on lattice stability, mechanical and thermodynamic properties of L1_2 structure Al_3Tm and Al_3Lu are studied by first-principles calculations within the VASP code. The phonon dispersion curves and density of phonon states are calculated by using the PHONONPY code. Our results agree well with the available experimental and theoretical values. The vibrational properties indicate that Al_3Tm and A_3Lu keep their dynamical stabilities in L1_2 structure up to 100 GPa. The elastic properties and Debye temperatures for Al_3Tm and Al_3 Lu increase with the increase of pressure. The mechanical anisotropic properties are discussed by using anisotropic indices AG, AU, AZ, and the threedimensional(3D) curved surface of Young's modulus. The calculated results show that Al_3Tm and Al_3Lu are both isotropic at 0 GPa and anisotropic under high pressure. In the present work, the sound velocities in different directions for Al_3Tm and Al_3Lu are also predicted under high pressure. We also calculate the thermodynamic properties and provide the relationships between thermal parameters and temperature/pressure. These results can provide theoretical support for further experimental work and industrial applications.     

CaS电子结构和热力学性质的第一性原理计算

根据密度泛函理论,采用平面波赝势和广义梯度方法,计算了Ca S的晶体结构和电子结构.通过准谐徳拜模型预测了硫化钙的体积变化率、体弹模量、热膨胀系数分别与温度和压强的变化关系,以及热容和温度的变化关系.