第一原理对铝的静态结构和相变的计算

 在密度泛函理论（DFT）和广义梯度近似（GGA）下,用缀加平面波加局域轨道（APW+lo）方法对铝的晶格常数、体弹模量以及在静态高压下的固态相变进行了计算。计算得出面心立方晶格结构（fcc）向六角密堆积结构（hcp）和体心立方结构（bcc）的相变分别发生在220 GPa和330 GPa,hcp向体心结构bcc的相变发生在380 GPa。计算结果和实验数据以及其它理论计算符合较好。     

固氩高压物态方程和弹性性质的密度泛函理论计算

 用平面波赝势方法结合局域密度近似密度泛函理论（DFT-LDA）计算了零温下固态氩晶体在压力0~82 GPa的p-V关系和弹性性质,计算结果与静高压实验数据符合较好,计算结果表明局域密度近似方法能较好地描述固氩晶体高压下的性质,采取合理的方法和计算参数,惰性气体固态晶体高压下的力学性质可以比较准确地计算出来,这可为一些还不能通过实验进行研究的物态分析提供借鉴。     

高压下Al的零温电子结构和物态方程

 用LMTO法（线性Muffin-Tin轨道法），计算了金属铝的超高压电子结构及零温物态方程。Al的压缩比到10，压力达10 TPa。根据第一原理计算结果，在带结构方面，s带总是处于Fermi能以下，随着压力的增加，s、p和d的带宽增加，随后则杂化程度增加，所以s、d轨道的电子占据数连续变化。上述变化对压力的影响也是连续的，换言之，从s→d转变没有理由说明Al在0.5 TPa附近冲击Hugoniot的斜率的拐弯现象。而这一点则是与Altshuler的观点不同。物态方程的第一原理计算结果表明，bcc结构比fcc结构要软，但因差别不很大，即使发生fcc→bcc的转变，也不会引起Hugoniot的质的变化（拐弯）。Al的LMTO物态方程还表明，我们以前所采用的半经验的冷压误差可达30%。为此，根据LMTO结果，给出新的冷压表达式p=∑_i=0⁵a_iδ^i/3，δ=Ω₀/Ω，其中a₀=-7.327 79，a₁=18.754 3，a₂=10.209 7，a₃=-2.523 53，a₄=-4.787 2，a₅=6.065 94，拟合误差小于3%。     

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

高压下SiO_2的第一性原理计算

基于密度泛函理论(DFT)的第一性原理,采用Hartree-Fork(HF)方法,分别计算了Si O2的α-石英结构、金红石结构以及氯化钙结构的总能量随体积的变化关系。利用Murnaghan状态方程,通过能量和体积拟合,得到了3种结构的体变模量及其对压强的一阶导数。计算结果表明,随着压强的增加,Si O2会从α-石英结构转变为金红石结构,与实验结果和其它理论结果一致;金红石结构与氯化钙结构之间不存在相变,可以共存。此外,对具有α-石英结构的Si O2的晶格常数、电子态密度和带隙随压强的变化关系进行了计算和分析,结果表明:加压作用下,能带向高能方向移动,Si─O键缩短,电子数转移增加,带隙展宽,电荷发生重新分布。     

稠密流体氮高温高压相变及物态方程

氮的高温高压物态方程以及相图对于研究和制备高能量密度含能材料至关重要.本文采用基于密度泛函理论的分子动力学模拟方法,研究了液氮的高温高压行为,给出900—25000 K, 2—200 GPa区间流体氮的物态方程以及组分、相态变化.在上述相空间,观察到流体氮分子相-聚合物相以及聚合物-原子相的相变发生.获得的液氮Hugoniot理论曲线与实验结果吻合较好,发现30—60 GPa区间Hugoniot曲线的软化与分子-聚合物流体相的相变有关;在60 GPa后Hugoniot曲线变陡峭与流体氮进入聚合物相区有关.     

加压下ZnO结构相变和电子结构的第一性原理计算

基于密度泛函(DFT)理论的第一性原理,计算半导体ZnO纤锌矿结构和岩盐矿结构状态方程及其在高压下的相变,分析加压下体相ZnO的晶格常数、电子态密度和带隙随压力的变化关系,并将计算结果与文献中的理论和实验数据进行比较.验证在计算金属氧化物时,应用局域密度(LDA)近似计算出的相变压力普遍偏高,采用广义梯度(GGA)近似得到的结果与实验符合较好.     

混合物物态方程的计算

在采用体积相加原理计算混合物物态方程的基础上,建立了一种物理模型确定混合物温度。根据混合物中各组分温度和压强平衡条件,采用压强-密度迭代方法计算给出混合物物态方程,编制了两种组分的混合物物态方程计算程序。为检验建立的温度模型的合理性及程序的有效性,分析了不同密度、温度状态的氢（H2）和钨（W）组成的混合物状态参量,计算了以下情形及其组合情形的混合物物态方程：H2和W以不同质量比混合;质量比固定,单组分状态不同;温度区间和密度区间不同。研究表明:实际应用中在建立的混合物温度模型基础上确定的混合物物态方程是合理的。     

混合物物态方程的计算

 在采用体积相加原理计算混合物物态方程的基础上,建立了一种物理模型确定混合物温度。根据混合物中各组分温度和压强平衡条件,采用压强-密度迭代方法计算给出混合物物态方程,编制了两种组分的混合物物态方程计算程序。为检验建立的温度模型的合理性及程序的有效性,分析了不同密度、温度状态的氢（H2）和钨（W）组成的混合物状态参量,计算了以下情形及其组合情形的混合物物态方程：H₂和W以不同质量比混合;质量比固定,单组分状态不同;温度区间和密度区间不同。研究表明:实际应用中在建立的混合物温度模型基础上确定的混合物物态方程是合理的。     

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

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

Nd(Fe,Si)11Cx化合物的全电子计算(x=0,2)

采用基于第一原理的全势能线性缀加平面波加局域轨道((L)APW lo)方法对Nd(Fe,Si)11Cx化合物(x=0,2)的电子结构进行了计算,得到了化合物态密度和磁矩等信息.计算结果表明NdFe9Si2化合物中Si原子主要与4b和32i位Fe原子产生杂化,导致Fe原子磁矩减小.NdFe9Si2C2化合物C原子使32i位Fe原子磁矩进一步降低,同时减弱了Si原子的影响,使得4b位Fe原子磁矩增大.     

Theoretical calculation for the equations of state and phase transitions of sodium and potassium hydrides

Abstract

The ionic overlap-compression model is used to calculate the equations of state, as well as the equilibrium properties, of sodium and potassium hydrides (NaH and KH). The present results agree with the experimental ones well. The NaC1-to-CsC1 phase transition pressures for both crystals are also determined. The agreement of the theoretical pressures (23.0 GPa for NaH and 4.9 GPa for KH) with the experimental measurements (29.3 GPa and 4 GPa) is rather good. The calculation shows that the effect of the zero-point vibration to the equilibrium properties and the transition pressures should not be ignored.     

Theoretical calculation for the equation of state and phase transition of lithium hydride crystal

Abstract

From the point of view of overlapping interactions between the nearest neighbours, while considering the compression effect of each ion, an ionic overlap-compression model is founded and applied to lithium hydride. The repulsive potential and cohesive energy curves of the crystal are calculated by a one-parameter variational method. The obtained equilibrium lattice constant (3.865 a₀), cohesive energy (? 218.82 kcal/mol), and bulk modulus (353 kbar) agree with experimental values surprisingly well. The calculated values of the equation of state also reach a good agreement with the experimental ones available below 40 kbar. A phase transition from NaCl to CsCl structure is predicted to occur around 0.85 Mbar, with a volume jump of about 6%.     

Second-order phase transitions and Ehrenfest equation for strained solids

The equations describing the second-order phase transitions at the hydrostatic and non-hydrostatic pressures are considered. It is shown that the proportionality coefficient between an “effective” volume and the true one V^??=?AV is inversely proportional to the compressibility of the solid at a uniaxial pressure and has a jump at the second-order phase transition. In the case of the non-hydrostatic pressure the “effective” volume of the solid is not a continuous function of temperature and has a jump at the phase transition as well. The Ehrenfest equation is generalized to the solids with an arbitrary homogeneous elastic deformation accompanied by change of the solid volume, in particular, to the solid strained by the uniaxial, biaxial or triaxial pressure. It is shown that the sum of the derivatives of the phase transition temperature with respect to uniaxial pressures applied along axes a, b, c does not coincide with the derivative of the phase transition temperature with respect to the hydrostatic pressure.     

The role of elastic strains at non-ferroelastic displacive phase transitions

The role of elastic strains at structural phase transitions is illustrated within the Landau theory and its first-order corrections due to critical fluctuations and defects are described. The Landau theory is sufficient to demonstrate the impossibility of bulk nucleation in a supercooled symmetrical phase and the absence of heterophase fluctuations in solids. The critical fluctuations are known to convert a second-order transition in an Ising-like system in a solid to a first-order one. Close to the mean-field tricritical point the effect can be described, for displacive systems, within a first-order perturbation theory and takes place for the Heisenberg systems as well. The influence of defects on these transitions is mediated essentially by the elastic strains. Defects smear the transition. For the “random local field” defects and an incommensurate (Heisenberg-like) transition this effect is so strong that first-order perturbation theory leads to a divergence.     

Influence of the nuclear equation of state on the hadron-quark phase transition in neutron stars

We study the hadron-quark phase transition in the interior of neutron stars, and examine the influence of the nuclear equation of state on the phase transition and neutron star properties. The relativistic mean field theory with several parameter sets is used to construct the nuclear equation of state, while the     

Structure phase transformation and equation of state of cerium metal under pressures up to 51 GPa

This study presents high pressure phase transitions and equation of states of cerium under pressures up to 51 GPa at room temperature. The angle-dispersive x-ray diffraction experiments are carried out using a high energy synchrotron x-ray source. The bulk moduli of high pressure phases of cerium are calculated using the Birch–Murnaghan equation. We discuss and correct several previous controversial conclusions, which are caused by the measurement accuracy or personal explanation. The c/a axial ratio of ε-Ce has a maximum value at about 29 GPa, i.e., c/a ≈ 1.690.     

Effect of high pressure on polymorphic phase transition and electronic structure of XAs (X=Al,Ga, In)

The structural and electronic properties of XAs (X = Al, Ga, In) under pressure have been investigated using ab-initio pseudo-potential approach within local density approximation in B3→B1→B2 phases. The values of phase transition pressures show reasonably good agreement with the experimental data and better than others. The B1→B2 phase transition in InAs is not seen. The volume collapse computed from equation of state (EOS) is found to be in good agreement with the experimental values. Under ambient conditions, the energy of B3 phase is lowest as compared to other phases, while at high pressures beyond B1→B2 phase transition, the energy of B2 phase is found to be lower than that of B1 phase showing correct stability of the phases. There is relatively smaller enthalpy associated with B3→B1 transition as compared to B3→B2 transition. The electronic structures have also been computed at different pressures. We have also reported the effect of pressure on energy gap and valence band width.