ZnS结构相变、电子结构和光学性质的研究

运用第一性原理平面波赝势和广义梯度近似方法, 对闪锌矿结构(ZB)和氯化钠结构(RS) ZnS的状态方程及其在高压下的相变进行计算研究, 分析相变点附近的电子态密度、能带结构和光学性质的变化机理. 结果表明: 通过状态方程得到ZB相到RS相的相变压强值为18.1 GPa, 而利用焓相等原理得到的相变压强值为18.0 GPa; 在结构相变过程中, sp³轨道杂化现象并未消失, RS相ZnS的金属性明显增强; 与ZB相ZnS相比, RS相ZnS的介电常数主峰明显增强, 并向低能方向出现了明显偏移, 使得介电峰向低能方向拓展, 在低能区电子跃迁大大增强. 关键词： 硫化锌 相变 电子结构 光学性质     

压力及杂质氢对金属锂弹性特性的影响

针对六角密堆金属锂16个原子超晶胞(supercell)、填隙一个氢原子的周期单元，采用基于密度泛函理论的平面波-赝势方法，研究了零温条件下压力及填隙氢掺杂对体系弹性性质的影响.结果表明氢掺杂导致体系的体模量增加.常压下掺杂体系的弹性常数C₁₁，C₃₃，C₆₆和C₁₂高于单质体系，剪切模量C₄₄有所下降，而C₁₃则与单质体系持平.压力作用下C₁₁，C₃₃和C₆₆一直大于单质体系，但C₁₂的值低于单质体系.在2GPa—4GPa压力区间内，弹性常数C₁₃呈反常变化，小于单质体系；在高压区掺杂体系的C₄₄和C₁₃则高于单质体系的相应量值，压力导致掺杂体系和单质体系之间剪切模的偏离加剧.掺杂体系在压力作用下依然保持压缩模的各向同性，具有和单质体系相似的特性. 关键词： 第一性原理 压力效应 弹性常数 金属锂     

ZnSe相变、电子结构的第一性原理计算

基于第一性原理平面波赝势(PWP)和广义梯度近似(GGA)方法,对闪锌矿结构(ZB)和岩盐结构(RS)的ZnSe在0—20GPa高压下的几何结构、态密度、能带结构进行了计算研究,分析了闪锌矿结构ZnSe和岩盐结构ZnSe的几何结构.在此基础上,研究了ZnSe的结构相变、弹性常数、成键情况以及相变压强下电子结构的变化机理.结果发现:通过焓相等原理得到的ZB相到RS相的相变压强为15.3GPa,而由弹性常数判据得到的相变压强为11.52GPa,但在9.5GPa左右并没有发现简单立方相的出现;在结构相变过程中,sp3轨道杂化现象并未消除,Zn原子的4s电子在RS相ZnSe的导电性中起主要贡献.     

立方结构Fe基磁性材料弹性系数第一性原理计算

通过赝势平面波法(CASTEP)及全电势线性缀加平面波法(FLAPW)，以bcc-Fe为对象，研究第一性原理计算立方结构Fe基磁性材料弹性系数的方法，分析影响计算立方结构Fe基磁性材料弹性系数准确度的各项因素. 结果表明，在第一性原理弹性系数计算中，晶格常数是决定弹性系数计算准确度的关键因素；势函数的选择也会影响计算准确度. 使用全电势基矢的FLAPW法可以得到更为精准的弹性系数计算结果. 计算得到bcc-Fe的弹性系数C₁₁，C₁₂，C₄₄分别为246 GPa，121 GPa，113 GPa，与实验值基本一致. 利用本方法，计算了新型Fe-Ga磁致伸缩材料的弹性系数C₁₁，C₁₂，C₄₄分别为207 GPa，166 GPa及108 GPa. 关键词： 弹性系数 磁致伸缩材料 赝势平面波法 全电势线性缀加平面波法     

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的关系. 关键词： 相变 热力学性质 第一性原理     

Investigations of phase transition, elastic and thermodynamic properties of GaP by using the density functional theory

The phase transition of gallium phosphide (GaP) from zinc-blende (ZB) to a rocksalt (RS) structure is investigated by the plane-wave pseudopotential density functional theory (DFT). Lattice constant a₀, elastic constants c_ij, bulk modulus B₀ and the pressure derivative of bulk modulus B₀' are calculated. The results are in good agreement with numerous experimental and theoretical data. From the usual condition of equal enthalpies, the phase transition from the ZB to the RS structure occurs at 21.9 GPa, which is close to the experimental value of 22.0 GPa. The elastic properties of GaP with the ZB structure in a pressure range from 0 GPa to 21.9 GPa and those of the RS structure in a pressure range of pressures from 21.9 GPa to 40 GPa are obtained. According to the quasi-harmonic Debye model, in which the phononic effects are considered, the normalized volume V/V₀, the Debye temperature θ, the heat capacity C_v and the thermal expansion coefficient α are also discussed in a pressure range from 0 GPa to 40 GPa and a temperature range from 0 K to 1500 K.     

压力对有机半导体均苯四甲酸晶体结构转变和电子性质的影响

研究高压条件下均苯四甲酸(C₁₀H₆O₈)材料的结构和性质对探索有机半导体材料的应用有积极意义.基于密度泛函理论的第一性原理赝势平面波方法,开展了0-300 GPa压强下C₁₀H₆O₈晶体的结构、电子和光学性质的研究.晶格常数在压强20 GPa和150 GPa下出现了明显跳变,且原子之间随着压强变化反复地出现成键/断键现象,表明压强可诱导晶体结构变化.电子结构的性质表明,0 GPa的C₁₀H₆O₈晶体是带隙为3.1 eV的直接带隙半导体,而压强增加到150 GPa时,带隙突变为0 eV,表明了晶体由半导体转变为导体.当压强为160 GPa时,晶体又变成了能隙约为1eV的间接带隙半导体,这可能是费米能级附近仅受O-2p轨道电子影响所导致.通过对C₁₀H₆O₈晶体介电函数的分析,再次验证了晶体在150 GPa时发生了结构相变.同时...     

一维应变加载下单晶铁结构转变的微观研究

采用嵌入原子势和分子动力学方法,模拟了单晶铁在一维应变条件下由体心立方（bcc）转变为六角密排（hcp）结构的微观过程. 当应变加载至相变临界值时,hcp相开始均匀形核并沿{011}晶面长大为薄片状体系.弹性常数C₃₁和C₃₂在相变前被逐渐硬化,C₃₃则在相变前出现软化行为;当体系完全相变后,上述各弹性常数显示开始随体积压缩而迅速硬化,温度效应对晶格具有软化作用,可削弱C₃₃的硬化和软化过程;样品在压缩过程可出现孪晶结构,孪晶结构使晶格发生剪切变形.混合相中,hcp相势能比bcc相高,最大剪应力方向与bcc相反向;系统的偏应力与hcp相质量分数近似呈线性关系. 关键词： 结构转变 分子动力学 一维应变     

BaF2高压相变行为的第一性原理研究

通过理论计算研究了BaF₂在高压下的晶体结构及物理性质.结果表明,在3.5和18.3 GPa,BaF₂依次经历了Fm■m-Pnma-P6₃/mmc两次结构相变,相变过程伴随着体积的塌缩,均为一级相变.约15 GP时,Pnma相晶轴压缩性出现异常,表现为随压强增大,晶轴b_o轻微增加,a_o略微减小.对其电子态密度进行分析发现,在16 GPa以后,由于F1原子的p_y+p_z与p_x轨道电子离域,导致其带隙随压强增加而降低.在约20 GPa时,Pnma相完全转变为P6₃/mmc相,相变完成.对BaF₂的拉曼峰位随压强变化进行了计算,为其高压拉曼光谱行为提供了相应的理论依据.计算了P6₃/mmc相在不同压强下的声子色散曲线,揭示了其卸压过程中的滞后机制,计算结果还预测该物相至少可以稳定到80 GPa.     

高压拉曼光谱方法研究联苯分子费米共振

测量了0—15 GPa压强下联苯分子的拉曼光谱. 结果表明,随压强增加,分子内和分子间π-π共轭和离域效应增强,谱线的绝对强度变大、蓝移. 联苯分子的两费米共振谱线强度比R_f/a减少,频率差Δ增加,当压强为8 GPa时,费米共振现象消失,利用Betran理论得出了固有频率差Δ₀和耦合系数ω随压强的变化关系,通过高压下相变进行了解释,并探讨了高压下费米共振耦合变弱的机理. 关键词： 联苯 费米共振 高压 拉曼光谱     

Structural transition behavior of ZnS nanotetrapods under high pressure

ZnS nanotetrapods synthesized via a solvothermal route have a octahedral core with a zincblende (ZB) structure and four hexprism-shaped arms consisting of alternately stacking ZB and wurtzite (WZ) phases, where the WZ phase has a higher volume percentage. In situ angular-dispersive X-ray diffraction (ADXRD) measurements were carried out to study the structural behavior of ZnS nanotetrapods under high pressure up to 41.3?GPa. The initial WZ structure exhibits a very high mechanical stability to ～11.3?GPa. Both the WZ and ZB structures transform to the rocksalt (RS) structure at ～15.4?GPa. The bulk moduli of the WZ (148.2?±?8.9?GPa) and RS (165.6?±?9.9?GPa) phases are both larger than the previously reported values. These phenomena are discussed based on the alternating epitaxial growth of the WZ and ZB phases in the arms of nanotetrapods. Our study suggests that the internal structure of nanomaterials could also greatly affect their stability and transition behavior.     

Elastic and Thermodynamic Properties of CdSe from First-Principles Calculations

The lattice parameters, bulk modulus, phase transition pressure, and temperature dependencies of the elastic constants c_ij of CdSe are investigated by using the Cambridge Serial Total Energy Package (CASTEP) program in the frame of Density Functional Theory (DFT). It is found that the phase transitions from the ZB structure to the RS structure and from WZ structure to RS structure are 2.2 GPa and 2.8 GPa, respectively. Our results agree well with the available experimental data and other theoretical results. The aggregate elastic modulus (B, G, E, A), the Poisson's ratio (υ), the Grüneisen parameter (γ), the Debye temperature Θ_D on pressure and temperature are also successfully obtained.     

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.     

Pressure-induced phase transition of B-type Y_2O_3

The synthesized monoclinic(B-type) phase of Y_2O_3 has been investigated by in situ angle-dispersive x-ray diffraction in a diamond anvil cell up to 44 GPa at room temperature. A phase transition occurs from monoclinic(B-type) to hexagonal(A-type) phase at 23.5 GPa and these two phases coexist even at the highest pressure. Parameters of isothermal equation of state are V_0= 69.0(1) ~3, K_0= 159(3) GPa, K_0= 4(fixed) for the B-type phase and V_0= 67.8(2) ~3, K_0= 156(3) GPa,K'_0= 4(fixed) for the A-type phase. The structural anisotropy increases with increasing pressure for both phases.     

Phase transitions in the superionic protonic conductor CsHSeO4 investigated by Brillouin spectroscopy

A Brillouin investigation in CsHSeO₄ has been performed over the temperature range 20–165 °C which includes two phase transitions, in particular the transition to the superionic phase near T_s = 129 °C. We observed strong discontinuities for elastic constants C₁₁, C₂₂ and C₃₃ at T_s and a broadening of the Brillouin lines above T_s. The results are discussed on the basis of a linear coupling between strains and mobile protons.     

High pressure polymorph of CdS predicted by first principles

Ab initio phonon calculations on CdS are performed to probe the high pressure structural behaviors. We predicted an unstable transverse acoustic (TA) mode for NaCl-CdS (B1) and a phase transition of B1→Pmmn driven by this soft mode is thus identified, excluding probable high pressure Cmcm phase. Furthermore, a softening TA phonon mode at the zone boundary M point of CsCl-CdS (B2) is predicted, which results in the phase transition from Pmmn to tetrahedral P4/nmm (B10). Enthalpy calculation reveals that Pmmn phase becomes energetically more favorable than the B1 phase over 51.2 GPa, and B10 phase is stable in a pressure range of 80.3-85.5 GPa, above which B10 phase will decompose into Cd and S.     

β-PtO2: Phononic,thermodynamic, and elastic properties derived from first-principles calculations

β-PtO₂ is a useful transition metal dioxide, but its fundamental thermodynamic and elastic properties remain unexplored. Using first-principles calculations, we systematically studied the structure, phonon, thermodynamic and elastic properties of β-PtO₂. The lattice dynamics and structural stability of β-PtO₂ under pressure were studied using the phonon spectra and vibrational density of states. The vibrational frequencies of the optical modes of β-PtO₂ increase with elevating pressure; this result is comparable with the available experimental data. Then, the heat capacities and their pressure responses were determined based on the phonon calculations. The pressure dependence of the Debye temperature was studied, and the results were compared in two distinct aspects. The elastic moduli of β-PtO₂ were estimated through the Voigt–Reuss–Hill approximation. The bulk modulus of β-PtO₂ increases linearly with pressure, but the shear modulus is nearly independent of pressure. Our study revealed that the elastic stiffness coefficients C₄₄, C₅₅ and C₆₆ play a primary role in the slow variation of the shear modulus.     

{Cu3}单分子磁体在磁场中的热纠缠

本文研究单分子磁体Na₉[Cu₃Na₃(H₂O)₉ (α-AsW₉O₃₃)₂]·26H₂O中三角自旋 环在磁场作用下的热纠缠性质, 利用数值计算求出任意两个Cu²⁺离子量子比特之间的配对纠缠度, 分别记为C₁₂, C₂₃和C₁₃. 研究结果表明, 磁场的方向和大小以及温度对配对纠缠度具有重要影响, 而且参数的变化对C₁₂, C₂₃和C₁₃的影响也是各不相同. 给出外加三个不同方向的磁场时, 配对纠缠度C₁₂, C₂₃和C₁₃各自对应的临界温度T_c随磁场强度的变化图, 由此可以得到单分子磁体三角自旋环中存在纠缠态的参数范围. 通过选择适当的磁场方向和大小以及温度等实验参数, 可以有效地调节和提高单分子磁体中的配对纠缠度. 关键词： 配对纠缠 单分子磁体 三角自旋环     

First principles study of structural stability,electronic structure and mechanical properties of ReN and TcN

The crystal structure, structural stability, electronic and mechanical properties of ReN and TcN are investigated using first principles calculations. We have considered five different crystal structures: NaCl, zinc blende (ZB), NiAs, tungsten carbide (WC) and wurtzite (WZ). Among these ZB phase is found to be the lowest energy phase for ReN and TcN at normal pressure. Pressure induced structural phase transitions from ZB to WZ phase at 214 GPa in ReN and ZB to NiAs phase at 171 GPa in TcN are predicted. The electronic structure reveals that both ReN and TcN are metallic in nature. The computed elastic constants indicate that both the nitrides are mechanically stable. As ReN in NiAs phase has high bulk and shear moduli and low Poisson's ratio, it is found to be a potential ultra incompressible super hard material.     

Pressure induced phase transition in ZnS

The phase transition of ZnS from the zincblende (ZB) structure to the rocksalt (RS) structure is investigated by the ab initio plane-wave pseudopotential density functional theory method. It is found that the pressures for transition from the ZB structure to the RS structure are 17.5 GPa from total energy-volume data and 15.4 GPa from equal enthalpies, consistent with the experimental data. From the high pressure elastic constants obtained, we find that the ZB structure ZnS is unstable when the applied pressure is larger than 17 GPa. Moreover, the dependence of the normalized primitive cell volume V/V₀ on pressure P can also be successfully obtained.