LMTO方法在β－Si_3N_4能带结构和静态性质研究中的应用

用ＬＭＴＯ－ＡＳＡ方法和局域密度泛函理论，从第一原理计算了晶体结构复杂的β－Ｓｉ＿３Ｎ＿４的电子能带结构、晶体的总能－体积曲线，由总能－体积曲线求得晶体的平衡晶格常数ａ＿０、体模量Ｂ＿０和结合能Ｅ＿（ｃｏｈ），引用Ｌｏｗｄｉｎ微扰法减少计算工作量，计算结果与准确性较高的从头算赝势法的计算结果基本相符，本文计算量比从头算赝势法少很多，仅用７４阶久期方程，便获得合理的结果。     

MgO能带结构的LMTO—ASA计算

来用LMTO-ASA能带从头计算方法计算MgO晶体的能带结构,根据添加空球与否和对d态处理方法的不同,分成四种不同的计算方案,结果表明,其中添加空球并利用Lowdin微扰法计入空d态的计算方案所得的结果与从头计算的赝势法的计算结果相当符合,该方案仅用16×16阶久期方程,计算量小,准确性高,可以给出MgO的合理的能带结构。     

由Al12C原子集团构造晶体的理论探讨

使用基于混合基表示的第一原理赝势法,研究了由Ａｌ１２Ｃ原子集团立方密堆积构造而成的晶体的电子和几何结构．对所构成的晶体的晶格常数和集团中各原子在晶体场下的位置改变都进行了没有可调参数的从头计算．结果表明,这样构成的晶体将是金属而非半导体,Ａｌ１２Ｃ原子集团之间在构成晶体时有较强的相互作用,晶场中的集团不再具有类似惰性原子的性质．     

求解多电子体系问题的赝势方法

文章介绍了在从头计算方法中加入赝势来求解多电子体系问题的方法。实践证明,用这种方法可以简化计 算,其精度与精确算法相比,误差小于10-3。     

硼酸盐晶体Li4CaB2O6的从头计算研究

本文通过总能赝势方法理论研究了八面体晶体材料Li4CaB2O6的晶体和电子结构以及Mulliken布居数；本文使用的总能赝势方法包括局域密度近似和广义梯度近似两种密度泛涵理论.所得的结果与实验提供的数据吻合较好.通过分析电子能带结构，可以推断多晶材料Li4CaB2O6 适用于半导体光学材料.     

有机桂中间体（H2SiNaF）稳定性的理论分析

在３００－１５００Ｋ范围内，用从头算分子轨道法计算了Ｈ２ＳｉＮａＦ异构化的热力学函数（热焓Ｈ，自由能Ｇ，熵Ｓ），平衡常数Ｋ（Ｔ），转化率Ｃ（Ｔ），Ａ因子和速度常数ｋ（Ｔ）动力学的研究是基于含Ｗｉｇｎｅｒ校正的过渡态理论和高斯－８６程序计算结果之上，得到构型１和２是稳定的，３和４将快速、完全地转化为１和２，并指出，在作稳定性分析时，除作量化计算外，必须考虑温度效应。     

Ge,Si的形变势及其应变层异质结的带阶

采用从头赝势能带计算方法研究了（００１）界面应变对Ｇｅ、Ｓｉ能带结构、平均键能和带阶参数的影响．计算了相应的形变势，并利用形变势的研究结果，采用平均键能方法计算了Ｇｅ／Ｓｉ应变层异质结在不同生长厚度ｈ（或平行晶格常数ａ∥）情况下的价带带阶和导带带阶．     

用第一性原理研究金属铜的电子结构

基于局域密度泛函理论,采用第一性原理赝势法,首先计算了晶体铜晶格参数,在此基础上计算并分析了铜的电子结构。计算机所得结果与其它计算及实验一致。     

镁基碱金属二元合金热力学性质的理论计算

本文运用Ｍｉｅｄｅｍａ理论，系统地计算了镁与碱金属组成的二元合金系统的液态混和焓和固态合金形成焓。根据理想溶液理论，还计算了Ｌｉ－Ｍｇ系统在１０００Ｋ时的熵、以及自由能－成分曲线，所得的结果同已有的实验结果和第一原理赝势方法计算的结果符合得比较好。     

利用两种不同赝势对铝热力学特性的研究

利用在局域密度近似下的第一原理平面波赝势方法以及应用TM和HGH两种模守恒赝势,计算了在一定晶格常数范围内铝的面心立方结构的能量与体积关系．在准谐近似下引入德拜模型,用两种赝势计算了不同的热力学特性,如体积弹性模量、物态方程、热膨胀等,比较了两种赝势计算出的结果,在通常情况下,两种赝势计算的结果差别不大,在有些情况下两者的差别比较明显,但计算结果都与有效的实验值相符合。     

THEORETICALANALYSISONTHESTABILITYOFORGANOSILICONINTERMEDIATES：H2SiNaF

在３００～１５００Ｋ范围内，用从头算分子轨道法计算了Ｈ２ＳｉＮａＦ异构化的热力学函数（热焓Ｈ，自由能Ｇ，熵Ｓ），平衡常数Ｋ（Ｔ），转化率Ｃ（Ｔ），Ａ因子和速度常数ｋ（Ｔ），动力学的研究是基于含Ｗｉｇｎｅｒ校正的过渡态理论和高斯－８６程序（６－３１Ｇ＊基组）计算结果之上，得到构型１和２是稳定的，３和４将快速、完全地转化为１和２，并指出，在作稳定性分析时，除作量化计算外，必须考虑温度效应。     

高压下BeSe和BeTe相变和热力学性质的第一性原理计算

用从头算平面赝势密度泛函理论计算了BeSe 和BeTe从闪锌矿结构(B3)到NiAs(B8)结构的相变压强, 所得的结果与其他理论值和实验值相吻合; 通过准谐德拜模型研究了两种结构下的德拜温度和热容随压强以及温度的变化关系.     

船体的极限承载能力

用摄动－迦辽金法求得矩形板的单位轴向缩矩－压缩载荷之间的关系．按照Ｓｍｉｔｈ法得到各单元的平均应力－平均应变曲线和船体梁总纵弯矩与曲率的关系曲线，编制了ＦＵＳＲＡ程序．船舶模型的数值计算表明，结果与试验值相吻合，具有简便、省时等优点．     

磷化铟从闪锌矿结构到氯化钠结构的相变

用两种方法研究了InP从闪锌矿结构到氯化钠结构的相变，第一种是基于等焓的方法，另外一种是根据柯西平衡判据.利用从头算平面波赝势密度泛函理论研究了InP在一系列高压下的电子结构和弹性常数，同时根据等焓条件发现InP从闪锌矿结构到氯化钠结构的相变压强是12.1 GPa，与根据柯西平衡判据得到的结果一致，从而说明用这两种方法来研究晶体的结构相变是可行的.     

求解多电子体系问题的赝势方法

章介绍了在从头计算方法中加入赝势来求解多电子体系问题的方法。实践证明，用这种方法可以简化计算，其精度与精确算法相比，误差小于10^-3。     

ZnO晶体结构参数的第一性原理计算

采用ABINIT软件计算了在Troullier-Martins赝势条件下不同晶格常数的纤锌矿结构ZnO晶体总能,获得晶体总能与晶格常数关系的数据,根据能量最小原理,获得ZnO稳定结构的晶格常数为a0=3.201,c0=5.161,c/a=1.612;对数据通过Muranghan状态方程进行拟和得到弹性模量B0=168GPa,计算结果比较接近实验值.     

从第一原理反演计算同类原子间的对势

在晶向族表象中，以三晶向族的Ｍｂｉｕｓ反演公式为零级近似，其余晶向族的贡献为微扰给出了陈－任微扰公式。内聚能函数对晶格间距的关系由第一原理的ＬＭＴＯ－ＡＳＡ方法计算得到。原子间对势以改进的陈－任公式计算。文中在没有考虑体积效应的情况下给出金和铜的同类原子对势的计算结果。     

钛原子位置对立方CaTiO3电子结构的影响

基于广义剃度近似（GGA）的从头算平面波超软赝势方法,计算了CaTiO3晶体的晶格常数、总能量以及结合能,得到了与实验相一致的结果。分析了Ti原子Z坐标位置对CaTiO3能带结构、态密度（DOS）等的影响。Ti原子位置强烈地影响CaTiO3的电子结构。随着Z坐标的增大,Ti的部分3d电子能量降低,O的部分2p电子能量升高,使价带与导带重叠,禁带消失。     

从头算研究SinHn（n=4,6,8）簇合物的电子结构

用ＧＡＵＳＳＩＡＮ９２从头算程序采用ＬＡＮＬＩＤＺ基组研究了ＳｉｎＨｎ（ｎ＝４,６,８）簇合物的电子结构。不管它们的几何构型和对称性,都具有相对较高的Ｓｉ－Ｓｉσ占据轨道和相对较低的Ｓｉ－Ｓｉσ＾＊空轨道。它们的ＨＯＭＯ－ＬＵＭＯ能隙接近,同一簇合物的不同低位空轨道之间的能隙也很小。     

量子化学从头计算FORTRAN IV通用程序MQAB-80

本文介绍一个用ＦＯＲＴＲＡＮⅣ语言编写的ＳＣＦ－ＬＣＡＯ－ＭＯ法从头计算通用程序。它可用于计算闭壳和开壳的、包含ｓ，ｐ，ｄ轨道的各种中小分子、离子和自由基。程序配备了ＲＨＦ和ＵＨＦ两种方法。藉助于局部数学模型的简化和计算方案的优化提高了运算速度。在本程序上采用不同的基集计算了若干分子和离子，结果与国外著名文献报导相符。程序对ＵＨＦ法作了尝试性改进并用于最低激发态的计算，取得较好的结果。本程序使用灵活方便，易读易懂，适于普及。