立方结构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. 关键词： 弹性系数 磁致伸缩材料 赝势平面波法 全电势线性缀加平面波法     

Al的微观组态与LaNi3.75Al1.25的结构和弹性第一性原理研究

在全电子水平上，采用广义梯度近似密度泛函理论和全势能线性缀加平面波方法并结合二维立方拟合方法，对LaNi_3.75Al_1.25合金的晶体结构与弹性性质进行了理论研究.计算结果给出合金的晶格常数a=b=0.5137　nm，c=0.4018　nm，Al原子在晶胞中的微观分布为同时占据部分3g和2c等价格位，弹性常数C₁₁+C₁₂=281.2，C₁₃=82.3，C₃₃=227.3，以及体弹性模量B＝124.5、切变模量G＝68.2　GPa.还对态密度、能带结构和电荷密度进行了计算分析，并给出材料LaNi_3.75Al_1.25的电子线性比热系数23.45　mJ/molK². 关键词： 3.75Al_1.25')" href="#">LaNi_3.75Al_1.25 储氢合金 全势能线性缀加平面波 弹性常数     

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

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

第一性原理研究硫化镉高压相变及其电子结构与弹性性质

采用第一性原理研究了CdS的六方纤锌矿(WZ), 立方闪锌矿(ZB) 和岩盐矿(RS)相在高压条件下的相稳定性、 相变点、电子结构以及弹性性能.WZ相与RS 相可以在相应的压强范围内稳定存在, 而ZB相不能稳定存在.压强大于2.18 GPa时, WZ相向RS相发生金属化相变.WZ相中S原子电负性大于Cd, 且电负性差值小于1.7, CdS的WZ相为共价晶体.高压作用下, S原子半径被强烈压缩, 有效核电荷增加, 对层外电子吸引能力提高, 电负性急剧增大, 导致S与Cd的电负性差值大于1.7, CdS的RS相以离子晶体存在. WZ相的C₄₄随压强增加呈下降趋势, 导致WZ相力学不稳定, 并向RS相转变.当压强大于2.18 GPa时, RS相C₁₁, C₁₂随压强增加而增大, 并且C₄₄保持稳定, 说明RS相具有良好的高压稳定性与力学性能. 关键词： 第一性原理 相变 电子结构 弹性性质     

Gruhn-Hess两体势修正模型向列相液晶微滴的Monte Carlo模拟

基于Gruhn-Hess两体势修正模型,用Monte Carlo方法模拟向列相液晶微滴.Gruhn-Hess两体势是空间各向异性的,修正模型的势参数依赖于液晶的弹性常数K₁₁、K₂₂、K₃₃及表面弹性常数K₁₃.假设向列相液晶微滴具有自由表面,在低温下计算分子二阶序参数和切向内禀强度在微滴内不同区域的变化,并与Gruhn-Hess两体势模型一作对比分析.结果表明:附加K₁₃项的修正模型在微滴表面产生内禀锚定,其大小决定了内禀易取向方向及内禀锚定强度的大小,且内禀强度的大小与K₃₃/K₁₁值有关;空间不完整的向列相液晶微滴由内层到外层有序度逐渐降低.     

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

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

L12型铝合金的结构、弹性和电子性质的第一性原理研究

运用第一性原理方法研究了L1₂型铝合金相Al₃Sc和Al₃Zr的晶体结构、电子结构和弹性.结合能和形成能的计算表明,两种合金具有较强的合金化能力,且Al₃Zr较Al₃Sc具有更强的结构稳定性.电子结构分析表明,费米能级以下较多的价电子数决定了Al₃Zr具有较强的结构稳定性.计算并分析比较了两种合金相的单晶弹性常数（C₁₁,C₁₂和C₄₄）以及多晶弹性模量（体弹性模量B、剪切模量G、杨氏模量Y、泊松比ν和各向异性因子A）.通过对比实验和其他理论计算结果,进一步分析和解释了两种合金相的力学性质. 关键词： 铝合金 第一性原理 结构和电子性质 弹性     

带有Dzyaloshinski-Mariya相互作用的两比特XXZ模型的纠缠量子热机

以带有Dzyaloshinski-Mariya 相互作用的两比特XXZ模型为工作物质构建纠缠量子热机. 在量子热力学平衡态下, 采用Kieu的形式描述了做功与热传递.对于不同的各向异性参数, 分析了热机循环中量子纠缠与热传递、做功以及机械效率等热力学量之间的关系. 结果表明: 在这个纠缠体系中, 热力学第二定律依然成立; 机械效率的等值线图是环状曲线; 当各向异性参数Δ较小时, 热机在C₁ > C₂ 和C₁ < C₂ 两区域运行, 当增大Δ值时, 热机只在C₁ > C₂ 区域运行. 关键词： 量子热机 并发度 机械效率     

水溶液中稀土离子与甘氨酸、丝氨酸、天冬氨酸作用的13C NMR研究

测定了酸性水溶液中甘氨酸、丝氨酸和天冬氨酸稀土络合物(Ln=La、Pr、Nd、Eu、Tb、Dy、Ho、Er、Tm和Yb)的¹³C诱导位移。对位移试剂的分析指出,三种氨基酸通过α-羧基以双齿形式配位于稀土,配位键长为0.23nm~0.25nm,天冬氨酸的γ-羧基也是配位基团。由本文与文献中已报道的各种氨基酸稀土络合物的¹³C诱导位移的系统分析表明,配体¹³C超精细偶合常数A值和结构因子G值有如下规律:(1)│A(C₀)│<│A(C_α)│;A(C₀)为正,A(C_α)为负;(2)│G(C₀)│>│G(C_α)│;配体碳核的G均为负值。     

高压下有序晶态合金Fe3Pt的低能声子不稳定性及磁性反常

有序晶态Fe₃Pt因瓦合金处于一种特殊的磁临界状态, 这种磁临界状态下体系的晶格动力学稳定性对压力极为敏感. 基于密度泛函理论的第一性原理的投影缀加平面波方法研究了不同晶态合金的Fe₃Pt的焓和磁性随压力的变化规律, 结果表明, 在压力小于18.54 GPa下, P4/mbm结构是热力学稳定的相. Pm3m结构、I4/mmm结构、DO₂₂结构的Fe₃Pt在铁磁性坍塌临界压力附近体系的总磁矩急剧下降并具有振荡现象, 且I4/mmm结构和DO₂₂结构的Fe₃Pt 在临界压力附近出现了Fe1原子磁矩反转现象. 在43 GPa下, DO₂₂结构的Fe₃Pt出现了亚铁磁微观磁特性突然增强且伴随着体积突然增大的现象. 在高压下, 对Pm3m结构Fe₃Pt的晶格动力学计算表明, 压力小于26.95 GPa的铁磁态下体系的自发磁化诱导了体系横向声学支声子软化, 表明体系中存在很强的自发体积磁致伸缩. 特别是在铁磁性坍塌临界压力41.9 GPa至磁性完全消失的57.25 GPa压力区间, 晶格动力学稳定性对压力更加敏感. 压力大于57.25 GPa时, 压力诱导了体系声子谱的稳定.     

First-principles study of structural and elastic properties of CrO2 at high pressure

The structural and elastic properties of CrO₂ in the rutile phase under high pressures have been investigated using pseudopotential plane-wave method based on density functional theory. The optimized lattice parameters and the bulk modulus at zero pressure agree well with available experimental and theoretical data. The elastic constants C ₁₁, C ₁₂, C ₄₄, C ₃₃, C ₁₃, and C ₆₆ at zero pressure are calculated to be 359.91, 264.69, 143.28, 309.45, 218.45, and 260.74 GPa, respectively. Elastic constants, bulk modulus, shear modulus, Young's modulus, and Poisson's ratio under pressures are obtained. Our results indicate that the rutile phase is mechanically stable below 11.99 GPa. The elastic anisotropy of rutile phase under pressures has also been predicted.     

Study of first principles on anisotropy and elastic constants of Y3Al2 compound

In this work primarily structural parameters of the Y₃Al₂ was optimized. In the calculations, PBE type Exchange-Correlation function was selected in the GGA approximation and ultra-soft pseudopotential was used. Crystalline lattice parameters a = 15.581 Bohr, c / a = 0.923. Total energy was calculated by small deformations in volume. These values were fitted to the equation of state and the bulk modulus was calculated as 58.4 GPa. The working structure has elastic constants as C₁₁, C₁₂, C₁₃, C₃₃, C₄₄ and C₆₆. These constants have been calculated under ambient pressure. Using the elastic constants, bulk, shear, Young modulus and Poisson ratio (in Voigt approach) were obtained as 58.12, 34.31, 86.01 GPa and 0.25, respectively. Melting temperature of the material was guessed from the elastic constant. The results were compared with the current experimental and theoretical data.     

First-principles study of elastic and electronic properties of layered ternary nitride SrZrN2 under pressure

The structural, elastic, and electronic properties of SrZrN₂ under pressure up to 100?GPa have been carried out with first-principles calculations based on density functional theory. The calculated lattice parameters at 0?GPa and 0?K by using the GGA-PW91-ultrasoft method are in good agreement with the available experimental data and other previous theoretical calculations. The pressure dependence of the elastic constants and the elastic-dependent properties of SrZrN₂, such as bulk modulus B, shear modulus G, Young's modulus E, Debye temperature Θ, shear and longitudinal wave velocity V_S and V_L, are also successfully obtained. It is found that all elastic constants increase monotonically with pressure. When the pressure increases up to 140?GPa, the obtained elastic constants do not satisfy the mechanical stability criteria and a phase transition might has occurred. Moreover, the anisotropy of the directional-dependent Young's modulus and the linear compressibility under different pressures are analysed for the first time. Finally, the pressure dependence of the total and partial densities of states and the bonding property of SrZrN₂ are also investigated.     

First-principles study of structural,elastic, electronic and thermodynamic properties of topological insulator Bi2Se3 under pressure

The structural, elastic, electronic and thermodynamic properties of the rhombohedral topological insulator Bi₂Se₃ are investigated by the generalized gradient approximation (GGA) with the Wu–Cohen (WC) exchange-correlation functional. The calculated lattice constants agree well with the available experimental and other theoretical data. Our GGA calculations indicate that Bi₂Se₃ is a 3D topological insulator with a band gap of 0.287 eV, which are well consistent with the experimental value of 0.3 eV. The pressure dependence of the elastic constants C_ij, bulk modulus B, shear modulus G, Young’s modulus E, and Poisson’s ratio σ of Bi₂Se₃ are also obtained successfully. The bulk modulus obtained from elastic constants is 53.5 GPa, which agrees well with the experimental value of 53 GPa. We also investigate the shear sound velocity V_S, longitudinal sound velocity V_L, and Debye temperature Θ_E from our elastic constants, as well as the thermodynamic properties from quasi-harmonic Debye model. We obtain that the heat capacity C_v and the thermal expansion coefficient α at 0 GPa and 300 K are 120.78 J mol^?1 K^?1 and 4.70 × 10^?5 K^?1, respectively.     

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.     

First principle calculation of elastic and thermodynamic properties of stishovite

Using ab initio plane-wave pseudo-potential density functional theory method,the elastic constants and band structures of stishovite were calculated.The calculated elastic constants under ambient conditions agree well with previous experimental and theoretical data.C13,C33,C44,and C66 increase nearly linearly with pressure while C 11 and C 12 show irregularly changes with pressure over 20 GPa.The shear modulus(C11-C12)/2 was observed to decrease drastically between 40 GPa and 50 GPa,indicating acoustic mode softening in consistency with the phase transition to CaCl 2-type structure around 50 GPa.The calculated band structures show no obvious difference at 0 and 80 GPa,being consistent with the high incompressibility of stishovite.With a quasi-harmonic Debye model,thermodynamic properties of stishovite were also calculated and the results are in good agreement with available experimental data.     

钒基固溶体储氢材料弹性性质第一性原理研究

采用基于密度泛函理论的第一性原理赝势平面波方法,计算了（59Cr-41Ti）_100－xV_x（x=5,15,30,60,80,100）六种钒基储氢合金的晶格常数、弹性性质和电子态密度,计算结果与实验值符合较好.发现当x=60时的钒基合金具有较好弹性性质,杨氏模量为14930 GPa,切变模量为5442 GPa及体弹模量为19396 GPa.结合实验循环性能分析认为在吸放氢过程中合金已经发生塑性变形,弹性 关键词： 钒基固溶体 储氢合金材料 密度泛函理论 弹性性质     

Characterisation of the high-pressure structural transition and elastic properties in boron arsenic

This paper carries out the First principles calculation of the crystal structures (zinc blende (B3) and rocksalt (B1)) and phase transition of boron arsenic (BAs) based on the density-functional theory. Using the relation between enthalpy and pressure, it finds that the transition phase from the B3 structural to the B1 structural occurs at the pressure of 113.42GPa. Then the elastic constants C₁₁, C₁₂, C₄₄, bulk modulus, shear modulus, Young modulus, anisotropy factor, Kleinman parameter and Poisson ratio are discussed in detail for two polymorphs of BAs. The results of the structural parameters and elastic properties in B3 structure are in good agreement with the available theoretical and experimental values.     

高压下TiC的弹性、电子结构及热力学性质的第一性原理计算

利用基于密度泛函理论的第一性原理平面波赝势方法 并结合准谐徳拜模型研究了NaCl结构的TiC在高压下的弹性性质、电子结构和热力学性质. 计算所得零温零压下的晶格常数、体弹模量及弹性常数与实验值符合得很好. 零温下弹性常数和弹性模量随压强增大而增大. 通过态密度和电荷密度的分析, Ti-C键随压强增大而增强. 运用准谐德拜模型, 成功计算了TiC在高温高压下的体弹模量、熵、热膨胀系数、徳拜温度、 Grüneisen参数和比热容. 结果表明压强对体弹模量、热膨胀系数和徳拜温度的影响大于温度对其的影响. 热容随着压强升高而减小, 在高温高压下, 热容接近Dulong-Petit极限.     

Structural,elastic and mechanical properties of orthorhombic SrHfO3 under pressure from first-principles calculations

Structural, elastic and mechanical properties of orthorhombic SrHfO₃ under pressure have been investigated using the plane-wave ultrasoft pseudopotential technique based on the first-principles density functional theory. The calculated equilibrium lattice parameters and elastic constants of orthorhombic SrHfO₃ at zero pressure are in good agreement with the available experimental and calculational values. The lattice parameters, total enthalpy, elastic constants and mechanical stability of orthorhombic SrHfO₃ as a function of pressure were studied. With the increasing pressure, the lattice parameters and volume of orthorhombic SrHfO₃ decrease whereas the total enthalpy increases. Orthorhombic SrHfO₃ is mechanically stable with low pressure (<52.9 GPa) whereas that is mechanically instable with high pressure (>52.9 GPa). The bulk modulus, shear modulus, Young's modulus and mechanical anisotropy of orthorhombic SrHfO₃ as a function of pressure were analyzed. It is found that orthorhombic SrHfO₃ under pressure has larger bulk modulus, better ductility and less mechanical anisotropy than orthorhombic SrHfO₃ at 0 GPa.