外压下VN相变和力学性质的第一性原理研究

运用基于赝势平面波基组的密度泛函程序VASP并结合Quantum ESPRESSO,Phonopy软件包对压力下VN的结构、力学性质、声子色散关系进行了第一性原理的研究.分别对Na Cl型(B1),Cs Cl型(B2),WC型(Bh)三种构型的VN进行了计算,三种结构的体积能量曲线、焓压关系和声子谱表明在常压下六角WC结构与立方结构相比更稳定.随着压力增加VN由Bh结构到B1结构的相变发生在30 GPa左右,而B1结构到B2结构的相变可能发生在150 GPa左右.常压下三种结构的VN是力学稳定的,其弹性常数和弹性模量都有随压强的增大而增加的趋势,三者都是脆性材料.B1结构和B2结构坐标基矢方向上的杨氏模量数值与体对角线方向上的差距较大,体现出明显的各向异性.随压力的增加B1结构各向异性程度增大而B2结构各向异性程度减小.     

外压条件下TaN及Ta-Si-N结构相变的第一性原理研究

利用基于第一性原理密度泛函理论的赝势平面波方法,对TaN、Ta-Si-N的B1型(NaCl)和B2型(CsCl)复合结构在压力下的的力学常数、体积-能量关系、焓-压关系以及声子色散关系图进行了计算和研究,并分析了他们的相对稳定性.根据计算结果,推断在90—110 GPa左右压力时,TaN薄膜将会发生由B1型向B2型的相变,且材料脆性增加;当压力在190—200GPa的范围内,Ta-Si-N薄膜也将发生B1→B2型的结构转变,由延展性材料压变为脆性材料.     

钝感高能炸药三氨基三硝基苯高温高压下热力学性质的分子动力学模拟研究

热力学性质是钝感高能炸药1, 3, 5-三氨基-2, 4, 6-三硝基苯(TATB)爆轰性质和安全性评估分析的重要参数. 由于结构的复杂性, TATB炸药尚缺乏系统的实验和理论计算结果. 结合全原子力场和分子动力学的方法, 本文系统研究了不同温度和压力条件下TATB的力学性质和热力学参数, 得到了弹性模量、德拜温度等随温度、压力的变化情况, 并与实验进行了对比分析. 结果表明: 在 0-50 GPa外部压力下, TATB晶体保持力学稳定, 弹性常数和弹性模量随压力升高而增大, 各向异性程度随压力升高而减小, 泊松比和延展性则受压力的影响较小; 随温度的升高, TATB的力学稳定性逐渐下降, 有发生力学失稳的可能, 各弹性常数随温度升高而逐渐减小, 各向异性程度也随之减小; TATB 的声速和德拜温度同样随着压力升高而增大, 平均声速从0 GPa下的1833 m/s, 增加到10 GPa 下的3143 m/s, 德拜温度由0 GPa下的254 K增加到10 GPa的587 K. TATB 热膨胀系数的计算表明, 在200-500 K 温度常压情况下, 其体热膨胀系数为35.9×10^-5 K^-1, 与实验数据符合较好.     

氮化铂在高压下的相变及弹性性质的第一原理计算

用PBE形式下的广义梯度近似（GGA）赝势平面波方法研究了氮化铂的结构相变以及弹性性质,计算了氮化铂的氯化钠（B1）、氯化铯（B2）、闪锌矿（B3）、纤维矿（B4）等四种结构并应用高压下的焓与压强的关系,得出在常温常压下B4结构是最稳定的结构,这与Yu 等人得的结果一致,且 B4→B1及B1→B2的相变压强分别发生在36.7 GPa和 185.4 GPa,同时,研究了B4结构在高压的弹性性质,发现弹性常数、体模量、剪切模量、压缩波速、剪切波速以及德拜温度均随着压强的增大而单调增大     

SmN晶体的电子结构和高压相变的第一性原理

利用基于密度泛函理论的第一性原理,研究了SmN晶体的电子结构和高压相变. SmN晶体的电子结构具有半金属特征,多数自旋电子显示金属导电性,少数自旋电子显示半导体导电性. 高压相变的结果显示,SmN晶体经历从NaCl型(B1)到CsCl型(B2)结构转变的压致结构相变,相变压力117 GPa. 弹性系数的结果显示,在环境压力下SmN晶体的弹性系数满足玻恩稳定条件,标志着B1相是力学稳定结构. 声子谱结果显示,在环境压力下B1相是热力学稳定结构,与弹性系数的计算结果一致.     

第一性原理研究压力对Fe2B的结构、机械性质的影响

基于密度泛函理论的第一性原理方法，系统研究了0～30 GPa压力下四方相Fe₂B的结构，稳定性和力学性质等.随着压力的增加，计算得到的晶格参数逐渐减小，所有结构均满足热力学稳定性，体积模量和剪切模量逐渐增大，韧性得到有效改善.硬度呈现先增大后减小的趋势，并在18 GPa压力时达到最大值，各向异性先增大后减小.德拜温度变化趋势与杨氏模量的变化趋势相同.同时还研究了相变结构——正交相的相关性质，结果表明其结构满足热力学稳定性但不满足力学稳定性，是否稳定存在还需要进行后续实验研究.     

氮化铂在高压下的相变及弹性性质的第一原理计算

用PBE形式下的广义梯度近似（GGA）赝势平面波方法研究了氮化铂的结构相变以及弹性性质，计算了氮化铂的氯化钠（B1）、氯化铯（B2）、闪锌矿（B3）、纤维矿（B4）等四种结构并应用高压下的焓与压强的关系，得出在常温常压下B4结构是最稳定的结构，这与Yu 等人得的结果一致，且 B4→B1及B1→B2的相变压强分别发生在36.7 GPa和 185.4 GPa，同时，研究了B4结构在高压的弹性性质，发现弹性常数、体模量、剪切模量、压缩波速、剪切波速以及德拜温度均随着压强的增大而单调增大     

压力下MoS_2的结构相变以及热动力学性质的第一原理研究

采用基于密度泛函理论的第一原理方法研究了层状MoS2在压力下的热动力学性质和相变机制.计算表明MoS2的2Hc结构在17.5GPa会相变到2 Ha结构,与此前理论结果20GPa基本一致.对比分析了两个结构在压力下的弹性常数、体模量、波速、德拜温度、线性体模量、热膨胀系数和定容热容等热动力学性质.研究表明MoS2的2 Hc和2 Ha结构在0~60GPa都满足力学稳定性条件,说明相变不是由于力学稳定性丧失导致,并且两个高压相在压力下呈现出较强的弹性各向异性,在0~50GPa内其a轴抗压缩均能力强于c轴.在相变机制上,Mulliken布居分析表明,随着压力增加,S原子向Mo原子转移电子以及Mo原子内s电子向d电子转移对MoS2从2 Hc结构相变到2 Ha结构起到重要作用.     

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

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

高压下两种8—羟基喹啉络合物的发光行为和结构变化

测定了在高压条件下两种金属(钙和锌）的8-羟基喹啉络合物的晶体粉末样品的发光行为和原位X光衍射光谱．结果表明,压力对其发光性质产生极大的影响．随着压力的增加,8—羟基喹啉钙的发光强度在3GPa以内时大大增加,随后发光强度快速下降．到7GPa左右时几乎为零,而8-羟基喹啉锌的发光强度随压力的增加而逐渐降低,到7GPa左右时约为常压的10％。高压下的原位X光衍射结果表明,8—羟基喹啉锌的晶体在3—4GPa开始发生非品化相变,在7GPa时该非晶化相变完成,样品的x光衍射完全消失．而8—羟基喹啉锌在压力的作用下(至16GPa)没有发生明显的相变。     

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.     

Elastic and thermodynamic properties of vanadium nitride under pressure and the effect of metallic bonding on its hardness

By the particle-swarm optimization method, it is predicted that tetragonal P42mc, 141md, and orthorhombic Amm2 phases of vanadium nitride （VN） are energetically more stable than NaCl-type structure at 0 K. The enthalpies of the predicted three new VN phases, along with WC, NaC1, AsNi, CsCl type structures, are calculated each as a function of pressure. It is found that VN exhibits the WC-to-CsCl type phase transition at 256 GPa. For the considered seven crystal- lographic VN phases, the structures, elastic constants, bulk moduli, shear moduli, and Debye temperatures are investigated. Our calculated equilibrium structural parameters are in very good agreement with the available experimental results and the previous theoretical results for the NaC1 phase. The Debye temperatures of VN predicted three novel phases, which are all higher than those of the remaining structures. The elastic constants, thermodynamic properties, and elastic anisotropies of VN under pressure are obtained and the mechanical stabilities are analyzed in detail based on the mechanical stability criteria. Moreover, the effect of metallic bonding on the hardness of VN is also investigated, which shows that VNs in P42mc, 141md, and Amm2 phases are potential superhard phases. Further investigation on the experimental level is highly recommended to confirm our calculations presented in this paper.     

Full potential calculation of structural, elastic properties and high-pressure phase of binary noble metal carbide: ruthenium carbide

We present in this paper the results of an ab initio theoretical study within the local density approximation (LDA) to determine in rock-salt (B1), cesium chloride (B2), zinc-blende (B3), and tungsten carbide (WC) type structures, the structural, elastic constants, hardness properties and high-pressure phase of the noble metal carbide of ruthenium carbide (RuC).The ground state properties such as the equilibrium lattice constant, elastic constant, the bulk modulus, its pressure derivative, and the hardness in the four phases are determined and compared with available theoretical data. Only for the three phases B1, B3, and WC, is the RuC mechanically stable, while in the B2 phase it is unstable, but in B3 RuC is the most energetically favourable phase with the bulk modulus 263 GPa, and at sufficiently high pressure (P_t=19.2 GPa) the tungsten carbide (WC) structure would be favoured, where ReC-WC is meta-stable.The highest bulk modulus values in the B3, B2, and WC structures and the hardnesses of H(B3)=36.94 GPa, H(B1)=25.21 GPa, and H(WC)=25.30 GPa indicate that the RuC compound is a superhard material in B3, and is not superhard in B1 and WC structures compared with the H(diamond)=96 GPa.     

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

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

Mechanical,electronic and thermodynamic properties of TE-C36 under high pressure

The structural parameters, mechanical, electronic and thermodynamic properties of TE-C36 under high pressure were calculated via the density functional theory in combination with the quasi-harmonic Debye model. The results show that the pressure has significant effects on the equilibrium structure parameters, mechanical, electronic and thermodynamic properties of TE-C36. The obtained ground state structural parameters are in good agreement with previous theoretical results. The mechanically and dynamically stable under pressure were confirmed by the calculated elastic constants and phonon dispersion spectra. The elastic constants, elastic modulus, B/G ratio, Poisson’s ratio and Vicker’s hardness were determined in the pressure range of 0–100?GPa. The elastic anisotropy of TE-C36 under pressure are also determined in detail. The electronic structure calculations reveal that TE-C36 remains a direct band gap semiconductor when the pressure changes from 0 to 100?GPa, and the band gap decreases with increasing pressure. Furthermore, the pressure and temperature dependence of thermal expansion coefficient, heat capacity and Debye temperature are predicted in a wide pressure (0–90?GPa) and temperature (0–2500?K) ranges. The obtained results are expected to provide helpful guidance for the future synthesis and application of TE-C36.     

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.     

Lattice instability and martensitic transformation in LaAg predicted from first-principles theory

The electronic structure, elastic constants and lattice dynamics of the B(2) type intermetallic compound LaAg are studied by means of density functional theory calculations with the generalized gradient approximation for exchange and correlation. The calculated equilibrium properties and elastic constants agree well with available experimental data. From the ratio between the bulk and shear moduli, LaAg is found to be ductile, which is unusual for B(2) type intermetallics. The computed band structure shows a dominant contribution from La 5d states near the Fermi level. The phonon dispersion relations, calculated using density functional perturbation theory, are in good agreement with available inelastic neutron scattering data. Under pressure, the phonon dispersions develop imaginary frequencies, starting at around 2.3 GPa, in good accordance with the martensitic instability observed above 3.4 GPa. By structural optimization the high pressure phase is identified as orthorhombic B(19).     

高压下γ'-Fe4N晶态合金的声子稳定性与磁性

利用基于密度泛函理论的第一性原理研究了高压下有序晶态γ’-Fe₄N合金的晶格动力学稳定性与磁性. 对比没有考虑磁性的γ’-Fe₄N的声子谱, 得出压力小于1 GPa时, 自发磁化诱导了铁磁相γ’-Fe₄N基态晶格动力学稳定. 压力在1.03-31.5 GPa时, Σ线上的点(0.37, 0.37, 0)、对称点X和M 上相继出现了声子谱软化现象. 压力在31.5-60.8 GPa时, 压致效应与自发磁化对诸原子的作用达到了稳定平衡, 表现出了声子谱稳定. 压力大于61.3 GPa时, 随着压力的增大压力诱导体系动力学不稳定性越强. 通过软模相变理论对于γ’-Fe₄N, 在10 GPa下的声学支声子的M点处软化现象的处理, 发现了动力学稳定的高压新相P2/m-Fe₄N. 压力小于1 GPa时高压新相P2/m-Fe₄N 是热力学稳定的相, 且磁矩与γ’-Fe₄N的磁矩几乎相同. 2.9-19 GPa时, P2/m相的焓比γ’相的焓小, 基态结构更稳定. 大于20 GPa时, 两相磁矩几乎相同.     

First-principles calculations of structural stability and mechanical properties of tungsten carbide under high pressure

The structural stability and mechanical properties of WC in WC-, MoC- and NaCl-type structures under high pressure are investigated systematically by first-principles calculations. The calculated equilibrium lattice constants at zero pressure agree well with available experimental and theoretical results. The formation enthalpy indicates that the most stable WC is in WC-type, then MoC-type finally NaCl-type. By the elastic stability criteria, it is predicted that the three structures are all mechanically stable. The elastic constants C_ij, bulk modulus B, shear modulus G, Young?s modulus E and Poisson?s ratio ν of the three structures are studied in the pressure range from 0 to 100 GPa. Furthermore, by analyzing the B/G ratio, the brittle/ductile behavior under high pressure is assessed. Moreover, the elastic anisotropy of the three structures up to 100 GPa is also discussed in detail.