材料第一性原理操作和试验方法课件

1、第一性原理计算的软件实现及其在材料设计中应用报告人：曹勇 博士生导 师 ：朱景川 教授2011.03.09第一性原理计算的软件实现及其在材料设计中应用报告人：曹勇 1 概述2 计算软件CASTEP的简介及工作环境3 模型搭建4 几何优化5 计算结果分析6 应用实例1 概述1 概述1 概述第一原理计算指的是只依赖于基本的物理常数和原子质量，通过量子力学理论求解运动方程来确定电子行为，进而预测物质性质。由于不借助任何经验、半经验理论或试验结果，因此也被称为从头算方法（ab initio method)。仅需5个基本物理常数，即电子的静止质量m、电子电量e、普朗克常数h、光速c和波尔玆曼常数kB 第

2、一原理计算方法是从最基本和最简单的物理事实出发，计算结果往往能够具有很高的精确度和可信度，成为独立与试验之外研究物质性质的又一利器。 1第一原理的可靠性1第一原理的可靠性Total Energy, Cohesive EnergyMetalDFT-LDAExperimentNb0.550.56Mo0.490.50Rh0.450.42Pd0.270.28Total energyof atoms (Ry.)Cohesive energy (Ry.)Total Energy, Cohesive EnergyMStructural Properties: MoleculesStructural prop

3、erties of diborane abd1 ()d2 ()DFT-LDA121o98o1.321.20DFT-GGA122o97o1.331.19EXP122o97o1.331.19Diborne (3c-2e bond)Structural Properties: MoleculLattice constant (a.u.) LDA predicts latticeconstants to within 2%.Bulk modulus (Mbar)LDA predicts bulk modulusto within a few percents.Structural Properties

4、LatticeconstantLattice constant (a.u.) LDA prMagnetism: BCC Fe(DFT-GGA)-2-1 0 1 2-8-6-4-2 0 2 4Energy (eV)DOSMajority spinMinority spinMeasured: mB=2.22/atomCalculated:mB=2.25/atomMagnetism: BCC Fe(DFT-GGA)-2-2-1 0 1 2-8-6-4-2 0 2 4Magnetism: HCP Co (DFT-GGA)Energy (eV)DOSMajority spinMinority spinM

5、easured: mB=1.72/atomCalculated:mB=1.63/atom-2-1 0 1 2-8-6-4-2 0 2 4MagnetMagnetism: CCP Ni(DFT-GGA)Energy (eV)DOSMajority spinMinority spin-1.5-1-0.5 0 0.5 1 1.5-8-6-4-2 0 2 4Measured: mB=0.62/atomCalculated:mB=0.65/atomMagnetism: CCP Ni(DFT-GGA)EneMagnetic PropertiesMetalCrystalstructureatomic con

6、figurationMeaured (mB/atom)Calculated (DFT) (mB)FeBCC3d64s22.222.25CoHCP3d74s21.721.63NiCCP3d94s20.620.65Magnetic PropertiesMetalCrysta由于在处理多电子体系相互作用时遇到的困难, 第一原理在实际计算中也需要做一定程度的近似，不过这样的近似也都是基于物理上的简单假设和数学形式上的合理简化而得到的。由于在处理多电子体系相互作用时遇到的困难, 第一原理在实际计Can the materials properties be understood knowing only

7、 the identities (atomic numbers) of their constituent atoms?CSiGeNBStructural properties?Electronic properties?Elastic properties?optical properties ?Phonon dispersion?Magnetic properties?Superconductivity?2第一原理所能解决的问题Can the materials properties bHexagoal Mo2CMo2C 是合金钢中的二次硬化相，它均匀、细小、弥散的分布于基体中，具有极高的

8、硬度和良好的高温稳定性。合金钢中合金元素种类繁多，是否存在其他的碳化物相，如Cr2C、V2C、Ni2C来取代其作用。 Mo2C 是否具有其他更加稳定的结构 空位或其他原子置换对其稳定性和性能有何影响2 第一原理所能解决的问题Hexagoal Mo2CMo2C 是合金钢中的二次硬化Orthorhombic and Hexagonal Fe3C晶体结构的稳定性判定虚拟相、亚稳相、稳定相。研究晶格常数和原子所处位置对晶体结构稳定性的影响。对于组成元素确定的体系，可能存在不同的晶体结构。分别对不同的晶体结构的能量极小值进行计算，比较不同晶体结构的能量极小值，便可确定稳定的晶体结构。2 第一原理所能解决

9、的问题Orthorhombic and Hexagonal Fe3Material Properties from First-PrinciplesFrom first principles!Predict new behaviors/properties of existing materialsDesign materials with desired propertiesUnderstand and explain materials propertiesBecoming realityMaterial Properties from First目前已有一些利用第一原理对具有周期性结构的

10、材料进行能带计算的商业软件，如VASP，ABINIT，Wien2K，CASTEP。目前已有一些利用第一原理对具有周期性结构的材料进行能带计算的2 计算软件CASTEP的简介及运行环境2 计算软件CASTEP的简介及运行环境计算软件CASTEP (Cambridge Sequential Total Energy Package) 是基于密度泛函理论的第一性原理计算量子力学程序。其使用的是总能量平面波赝势的方法。密度泛函理论 (Density Function Theory) 巧妙的将电子之间的交换势表示为密度泛函，从而使得薛定谔方程在考虑了电子之间的复杂相互作用之后，依然可以用自洽的方法求解。

11、计算软件CASTEP (Cambridge Sequenti关于CASTAP CASTAP是特别为固体材料学而设计的一个现代的量子力学基本程序，其使用了密度泛函(DFT）平面波赝势方法，进行第一原理量子力学计算，以探索如半导体，陶瓷，金属，矿物和沸石等材料的晶体和表面性质。 典型的应用包括表面化学，键结构，态密度和光学性质等研究， CASTAP也可用于研究体系的电荷密度和波函数的3D形式。此外， CASTAP可用于有效研究点缺陷（空位，间隙和置换杂质）和扩展缺陷（如晶界和位错）的性质。 Material Studio使用组件对话框中的CASTAP选项允许准备，启动，分析和监测CASTAP服役工

12、作。计算：允许选择计算选项（如基集，交换关联势和收敛判据），作业控制和文档控制。分析：允许处理和演示CASTAP计算结果。这一工具提供加速整体直观化以及键结构图，态密度图形和光学性质图形。关于CASTAP CASTAP计算 CASTAP计算是要进行的三个任务中的一个，即单个点的能量计算，几何优化或分子动力学。可提供这些计算中的每一个以便产生特定的物理性能。性质为一种附加的任务，允许重新开始已完成的计算以便产生最初没有提出的额外性能。在CASTAP计算中有很多运行步骤，可分为如下几组：* 结构定义：必须规定包含所感兴趣结构的周期性的3D模型文件，有大量方法规定一种结构：可使用构建晶体（Build

13、 Crystal)或构建真空板(Build Vacuum Stab)来构建，也可从已经存在的的结构文档中引入，还可修正已存在的结构。提示： CASTAP计算所需时间随原子数平方的增加而增加。因此，建议是用最小的初晶胞来描述体系，可使用BuildSymmetryPrimitive Cell菜单选项来转换成初晶胞。CASTAP计算 CASTAP中选择一项任务1 从模块面板（Module Explorer)选择CASTAPCalculation。2 选择设置表。3 从任务列表中选择所要求的任务。 计算设置：合适的3D模型文件一旦确定，必须选择计算类型和相关参数，例如，对于动力学计算必须确定系综和参数

14、，包括温度，时间步长和步数。选择运行计算的磁盘并开始CASTAP作业。* 结果分析：计算完成后，相关于CASTAP作业的文档返回用户，在项目面板适当位置显示。这些文档的一些进一步处理要求获得可观察量如电子结构图、光学性质。 CASTAP中选择一项任务 计算设置：合适的3CASTAP能量任务CASTAP能量任务允许计算特定体系的总能量以及物理性质。除了总能量之外，在计算之后还可报告作用于原子上的力；也能创建电荷密度文件；利用材料观测仪（Material Visualizer)允许目测电荷密度的立体分布；还能报告计算中使用的Monkhorst-Park的k点的电子能量，因此在CASTAP分析中可生

15、成态密度图。对于能够得到可靠结构信息的体系的电子性质的研究，能量任务是有用的。只要给定应力性质，也可用于计算没有内部自由度的高对称性体系的状态方程（即压力-体积，能量-体积关系）。CASTAP中能量的缺损单位是电子伏特(eV)，各种能量单位的换算关系见Mohr.P.J(2000).1 eV=0.036749308 Ha=23.0605 kcal/mole=96.4853 kJ/moleCASTAP能量任务CASTAP几何优化任务 CASTAP几何优化任务允许改善结构的几何，获得稳定结构或多晶型物。通过一个迭代过程来完成这项任务，迭代过程中调整原子坐标和晶胞参数使结构的总能量最小化。 CASTA

16、P几何优化是基于减小计算力和应力的数量级，直到小于规定的收敛误差。也可能给定外部应力张量来对拉应力、压应力和切应力等作用下的体系行为模型化。在这些情况下反复迭代内部应力张量直到与所施加的外部应力相等。 几何优化处理产生的模型结构与真实结构紧密相似。利用CASTAP计算的晶格参数精度列于右图。CASTAP几何优化任务 CASTAP几何优CASTAP性质任务 CASTAP性质任务允许在完成能量，几何优化或动力学运行之后求出电子和结构性质。可以产生的性质如下：* 态密度（DOS）：利用原始模拟中产生的电荷密度和势能，非自恰计算价带和导带的精细Monkhorst-Pack 网格上的电子本征值。* 带结

17、构：利用原始模拟中产生的电荷密度和势能，非自恰计算价带和导带的布里渊区高对称性方向电子本征值。* 光学性质：计算电子能带间转变的矩阵元素。CASTAP分析对话可用于生成包含可以测得的光学性质的网格和图形文件。* 布局数分析：进行Mulliken 分析。计算决定原子电荷的键总数和角动量（以及自旋极化计算所需的磁矩）。任旋地，可产生态密度微分计算所要求的分量。* 应力：计算应力张量，并写入seedname.castep 文档。CASTAP性质任务启动程序启动程序CASTAP的运行环境CASTAP的运行环境材料第一性原理操作和试验方法课件材料第一性原理操作和试验方法课件材料第一性原理操作和试验方法课

18、件材料第一性原理操作和试验方法课件材料第一性原理操作和试验方法课件2.8.3 模型搭建2.8.3 模型搭建搭建晶胞模型能带计算过程输出结果分析包括电子结构和物理性质通过计算能量极小值进行几何优化确立稳定结构搭建晶胞模型能带计算过程输出结果分析包括电子结构和通过计算能MoSi2以MoSi2为例， 其空间群为14/mmm； 晶格常数为a=3.202，c=7.851； Mo的原子位置为(0,0,0)和(1/2,1/2,1/2)； Si的原子位置为(0,0,1/3),(0,0,2/3),(1/2,1/2,1/6)和(1/2,1/2,5/6)MoSi2以MoSi2为例，空间群是点阵、平移群(滑移面和螺旋

19、轴)和点群的组合。 230个空间群是由14个Bravais点阵与32个晶体点群系统组合而成。能使无限大晶体(从微观角度来看)自身重复的几何对称操作的集合。 原胞初级单（晶）胞原始单（晶）胞primitive (unit) cell惯用晶胞Conventional (unit) cell空间群是点阵、平移群(滑移面和螺旋轴)和点群的组合。 230由一组基矢所决定的平行六面体所围起来的最小重复单元就叫原胞（或初基单胞Primitive Cell。(固体物理常用)原胞代表的是几何形状、体积的最小重复单元，不具有物理内涵。原胞往往不能反映晶体的对称性，因而，习惯上常选择能反映晶体对称性的重复单元，这种

20、重复单元就叫惯用晶胞（conventional cell）或非初基单胞（nonprimitive cell)。晶胞一般不是最小的重复单元。其体积（面积）可以是原胞的数倍。 由一组基矢所决定的平行六面体所围起来的最小重复单元就叫原胞（寻找已经存在的常见晶体模型寻找已经存在的常见晶体模型寻找已经存在的常见晶体模型寻找已经存在的常见晶体模型寻找已经存在的常见晶体模型寻找已经存在的常见晶体模型寻找已经存在的常见晶体模型寻找已经存在的常见晶体模型建立所研究材料的结构模型 空间群 晶格参数 晶体的原子内坐标以VC为例， 其空间群为FM-3M； 晶格常数为a=4.160； V的原子位置为(0,0,0)； C

21、的原子位置(1/2,1/2,1/2).建立所研究材料的结构模型搭建模型搭建模型选择空间群、设定晶格常数选择空间群、设定晶格常数选择空间群、设定晶格常数选择空间群、设定晶格常数材料第一性原理操作和试验方法课件选择空间群、设定晶格常数选择空间群、设定晶格常数添加原子，输入原子内坐标添加原子，输入原子内坐标材料第一性原理操作和试验方法课件材料第一性原理操作和试验方法课件材料第一性原理操作和试验方法课件材料第一性原理操作和试验方法课件材料第一性原理操作和试验方法课件材料第一性原理操作和试验方法课件建立超晶胞建立超晶胞转换原胞转换原胞设立P1空间群，降低对称性设立P1空间群，降低对称性问题1 说出上述晶

22、胞的晶胞类型，并指出后两个晶胞所含化学式VC的个数。问题2 P1空间群有何特点，什么情况下会将晶胞设置为P1空间群？问题1 说出上述晶胞的晶胞类型，并指出后两个晶胞所含化学式V2.8.4. 几何优化2.8.4. 几何优化能带计算的过程晶体是一个具有周期性结构的体系，输入时只能给出一个体积有限的晶体结构模型，需要利用周期性边界条件，才能得到整个晶体的能带结构。在选择了计算任务之后，需要对计算参数进行设置，如自洽场计算的精度、基组大小、k的取值等。能带计算的过程晶体是一个具有周期性结构的体系，输入时只能给出在能带计算中，基组是指用于线性组合成单电子波函数的基函数集合，在平面波赝势方法中，给出能量的

23、截断值Ecut。由于计算量的的原因计算时k值只取简约布里渊区中的有限值，以了l1 l2 l3的形式表示所取的k在b1、b2和b3方向上取值间隔分别为b1 / l1 , b2 / l2 , b3 / l3 , 。在其他条件不变的情况下，k的取值数目增加，得到的En(k) 函数的精确度增加，但计算量显著增加。自洽场计算的收敛用电子数密度或晶体总能量的收敛来标志，自洽场计算的精度是指收敛的标准。能带计算的过程在能带计算中，基组是指用于线性组合成单电子波函数的基函数集合 利用自洽场方法求解Kohn-Sham方程，得到所设结构的晶体总能量。在计算晶体的物理性质之前必须对所设的晶体结构模型进行几何优化，根

24、据关于能量、力、应力、位移的判据来判断晶体结构是否为稳定结构（总能量最小）。如果晶体结构不是稳定结构，重新设置晶格参数进行计算，直至得到稳定的晶体结构。然后对结构优化后的晶体进行物理性质计算，输出计算结果。能带计算的过程 利用自洽场方法求解Kohn-Sham方程，得到所设结构的晶几何优化几何优化是通过调节结构模型的几何参数来获得稳定结构的过程，其结果是使模型结构尽可能的接近真实结构。进行几何优化的判据可以根据研究的需要而定，一般是几个判据的组合使用。常用的判据有以下几个：自洽场收敛判据。对给定的结构模型进行自洽场计算时，相继两次自洽计算得到的晶体总能量之差足够小，即相继两次自洽计算的晶体总能量

25、之差小于设定值。力判据。每个原子所受的晶体内作用力足够小，即单个原子受力小于设定最大值。应力判据。每个结构模型单元中的应力足够小，即应力小于设定的最大值。位移判据。相继两次结构参数变化引起的原子位移的分量足够小，即原子位移分量小于设定的最大值。几何优化几何优化是通过调节结构模型的几何参数来获得稳定结构的材料第一性原理操作和试验方法课件设置参数设置参数设置参数设置参数设置参数设置参数设置参数设置参数设置电子结构参数设置电子结构参数几何优化几何优化几何优化(a)(b)MoMoCCFig. 1 The crystal structures of Mo2C with 12 atoms basis in

26、 (a) Ortho-Mo2C, (b) Hexa-Mo2C. 几何优化CFig. 1 The crystal struct材料第一性原理操作和试验方法课件5 计算结果分析5.1 能量分析5.2 电子结构分析5.3 弹性性质分析5.4 其他物理性质分析5 计算结果分析5.1 能量分析5.1 能量分析5.1 能量分析晶体总能量晶体总能量（不包括核的动能部分）可分为两部分：一部分是原子核与内层电子组成的离子实能量，这部分能量基本上与晶体结构无关，是一个常数，赝势方法中常把总能量中这部分不变的能量设为零；另一部分是总能量与离子实能量之差，包括离子实与价电子的相互作用，离子实之间的相互作用以及价电子之

27、间的相互作用。晶体总能量晶体总能量（不包括核的动能部分）可分为两部分：一部材料第一性原理操作和试验方法课件晶体总能量晶体总能量M. T. Yin和M. L. Cohen成功地利用第一性原理从头算法研究了Si和Ge的晶体结构，指出常压下金刚石结构在能量上最稳定，算出的晶格常数与实验值吻合的非常好。他们还从单点能的体积变化率预言了压力导致的相变，即当压力增加体积缩小时，其他相的单点能可能比常压下稳定相的单点能更低，从而发生相变。计算结果表明，虽然六角金刚石是个次稳定的结构，其单点能极小值仅次于金刚石结构，但在稳定结构附近的单点能体积变化率比金刚石结构大，而-锡结构的单点能体积变化与金刚石结构一致，

28、因此，他们预言，金刚石结构在静压力下将发生向-锡结构而不是六角金刚石结构相变，这一预言已被试验证实。M. T. Yin和M. L. Cohen成功地利用第一性原 晶体的稳定性与其结合能密切相关，结合能定义为 晶体分解为单个原子所需要的能量。因此结合能越大，晶体结构越稳定。其中，x，y分别为晶胞中A，B原子的个数；Etot为晶胞的总能量；EatomA,EatomB分别为自由状态下原子A，B的能量。结合能其中，x，y分别为晶胞中A，B原子的个数；Etot为晶胞的总FeCrMoVNbCEsolid （eV）865.732468.931936.861979.071551.63155.83EsolidA

29、,EsolidB分别为固态单质单个原子A，B的能量。形成能FeCrMoVNbCEsolid （eV）865.732Hexagoal and orthorhombic Mo2C形成能Hexagoal and orthorhombic Mo2C设定任务设定任务5.2 电子结构分析5.2 电子结构分析能带结构能态密度布局分析电荷密度差异能带结构性能参数性能参数能带结构单电子方程的本征值En(k)对每个n是一个对k准连续的、可区分（非简并情况）的函数，称为能带。所有的能带称为能带结构。人们常常关心由原子的价电子形成的能量相对较高的几个能带，其中在绝对零度(0K)时，被价电子填满的能带称为价带，而未被填

30、满或全空的能带称为导带。导带底与价带顶之间的能量区间称为禁带，导带底与价带顶的能量之差称为禁带宽度。能带结构单电子方程的本征值En(k)对每个n是一个对k准连续能带结构下图是CASTEP计算的闪锌矿结构的BN能带结构图，总能量的计算精度是1.010-5，k的取值为888，能量的截断值为250.0eV。图中虚线表示费米能级(EF=0),费米能级以下的能带是价带，费米能级以上的能带是导带。导带底与价带顶位于相同点的能带结构，称为直接跃迁型能带结构。导带底与价带顶位于不同点的能带结构，称为间接跃迁型能带结构。从图中可以看出，导带底与价带顶位于布里渊区的原点（点，图中以G点表示），这表明BN具有直接型

31、能带结构。能带结构下图是CASTEP计算的闪锌矿结构的BN能带结构图，能带结构由于三维晶体的波矢k也是三维的，En(k)需要四维空间，因此，一般使波矢k沿选定的方向取值，画出二维的En(k)图。所选定的直线方向一般是晶体倒易点阵的高对称方向，如立方晶体倒易点阵的轴（100方向）、轴（110方向）和轴（111方向）。全部En(k)函数给出能带结构，实际上人们最关心的是费米能级附近的一系列能带，因此，一般只计算有限个能带的En(k)函数。从能带结构图上可以直观地看到在指定方向上各能带En(k)函数随k的变化、导带底与价带顶的位置、禁带宽度及禁带能隙随k的变化。能带结构由于三维晶体的波矢k也是三维的

32、，En(k)需要四维空能态密度在原子中电子的本征态形成一系列分立的能级，可以具体标明各能级的能量，说明它们的分布情况。然而，在固体中，电子能级非常密集，形成准连续分布，标明其中每个能级的能量是没有意义的。为了描述固体能带中电子能级的分布，引入“能态密度”的概念。若能量在E(E+E)范围内的电子能态数目为Z，则N(E)称为能态密度（或态密度）。能态密度在原子中电子的本征态形成一系列分立的能级，可以具体标能态密度能态密度描述了电子态的能量分布，由态密度的定义式可推出，第n个能带的态密度Nn(E)为式中，为原胞体积，N为晶体中原胞数目，积分在布里渊区内进行。总态密度N(E)为所有能带的态密度之和，即

33、积分在等能面SE上进行。能态密度能态密度描述了电子态的能量分布，由态密度的定义式可推能态密度总态密度N(E)是各能带的态密度之和，总电子数N等于N(E)从负无穷到费米能级EF的积分，即态密度是一个十分有用的概念，使用态密度可以用对电子能量E的积分代替在布里渊区内对k的积分。态密度经常用于电子结构的快速可视分析，价带宽度、能隙及电子态密度N(E)的主要特征处强度和数目等特性有助于定性解释实验得到的光谱数据。态密度分析还有助于理解电子结构的变化，如外压引起的电子结构的变化。能态密度总态密度N(E)是各能带的态密度之和，总电子数N等于能态密度能态密度（DOS）可分为总态密度、分波态密度(PDOS)和

34、局域态密度(LDOS)。局域态密度和分波态密度是电子结构分析十分有用的半定量工具。LDOS显示系统中原子的电子态对能态密度谱的每个部分的贡献。PDOS根据电子态密度角动量来进一步分辨这些贡献，确定DOS的主要峰是否具有s, p或d电子的特征。LDOS和PDOS分析可对体系中电子杂化的本质和体系的XPS谱、光谱中的主要特征的来源提供定性解释。PDOS计算基于Mulliken布局分析，这种布局分析可将原子对每个能带的贡献归属于指定的原子轨道，对所有能带中某些原子轨道的贡献求和，可得到加权的态密度。可选择不同的加权方式，例如，将指定原子的所有原子轨道的对各能带的贡献加起来便得到LDOS。能态密度能态

35、密度（DOS）可分为总态密度、分波态密度(PDO能态密度能态密度结果输出结果输出材料第一性原理操作和试验方法课件能态密度在自旋极化体系中，电子和电子具有不同的空间波函数，即占据不同的能态。可分别计算电子的态密度N(E)和电子的态密度N(E) ，它们的和给出总态密度，它们的差N(E)- N(E) 被称为自旋态密度（SDOS）。材料的磁性质与自旋态密度有关。Fe是典型的自旋极化体系，用CASTEP计算-Fe的能态密度，总能量的计算精度是1.010-5，k的取值为101010，能量的截断值为300.0eV，分别得到电子和电子的能态密度，如图所示。能态密度在自旋极化体系中，电子和电子具有不同的空间波函

36、数布局分析对电子电荷在各组分原子之间的分布情况进行计算，称为布局分析。有多种布局分析方法，其中被广泛采用的布局分析方法是Mulliken布局分析。布局分析可以给出原子上、原子轨道上、两原子间的电子电荷分布，依次称为原子布局、轨道布局和键布局。布局分析对电子电荷在各组分原子之间的分布情况进行计算，称为布布局分析布局分析为原子间成键提供了一个客观判据，并且两原子间的重叠布局还可用于评价一个键的共价性或离子性。键布局的值高表明键是共价的，键布局的值低表明键是一种离子相互作用。还可以用有效离子价来进一步评价键的离子性，有效离子价定义为阴离子物种上原来的离子电荷与Mulliken电荷之差，若这个值为零，

37、则表明该键是完全离子键，若这个值大于零，则表明该键的共价成分增加。布局分析布局分析为原子间成键提供了一个客观判据，并且两原子间电荷密度差异电荷密度差异利用能带理论除了计算En(k)函数、能态密度等描述能带结构的物理量，还可以对固体的一些物理性质，如弹性常数、光学性质、热力学函数等，进行计算。计算固体的物理性质时，应先进行几何优化。利用能带理论除了计算En(k)函数、能态密度等描述能带结构的5.3 弹性性质输出与分析5.3 弹性性质输出与分析弹性性质材料的弹性常数描述了它对所加应力的相应。或者反过来说，弹性常数描述了为维持一个给定形变所需的应力。应力 和应变 均为二阶对称张量，可分别用 (i=1

38、, 2,6)和 (j=1, 2,6)来表示，则线弹性常数 可以表示为一个66的对称矩阵对于小的应力应变，有知道晶体总能量ET可算出弹性常数 。利用计算得到的弹性常数 ，还可计算体弹模量、泊松比的性质。弹性性质材料的弹性常数描述了它对所加应力的相应。或者反过来说材料第一性原理操作和试验方法课件待优化结构完成之后待优化结构完成之后VC弹性常数VC弹性常数材料第一性原理操作和试验方法课件材料第一性原理操作和试验方法课件弹性常数弹性常数5.4 其他物理性质分析5.4 其他物理性质分析热力学性质对体系热力学性质的描述基于声子，声子是晶格震动的能量子。声子的角频率与波矢q的函数关系(q)称为声子谱或色散关

39、系。利用第一原理计算声子谱(q)的方法有两种：超胞法和线性响应法。由声子谱(q)可计算体系的焓H、熵S自由能F和晶格热容CV，它们都是温度的函数。下图是用CASTEP计算的立方ZrO2的焓(H)、功函数(F)、温度与熵的乘积(TS)和等容热容(CV)随温度的变化。热力学性质对体系热力学性质的描述基于声子，声子是晶格震动的能光学性质光学性质光学性质光学性质6 应用实例6 应用实例晶体结构稳定性对于组成元素确定的体系，可能存在不同的晶体结构。分别对不同的晶体结构的能量极小值进行计算，比较不同晶体结构的能量极小值，便可确定稳定的晶体结构。Orthorhombic and Hexagonal Fe3C

40、(a)(b)MoMoCCFig. 1 The crystal structures of Mo2C with 12 atoms basis in (a) Ortho-Mo2C, (b) Hexa-Mo2C. 晶体结构稳定性对于组成元素确定的体系，可能存在不同的晶体结晶体结构稳定性对于相同结构的晶体结构，研究空位或原子置换对晶体结构稳定性的影响。通过对能量极小值进行计算，比较不同处理后的能量极小值，便可确定稳定的晶体结构。B. Xiao等b还利用第一原理和试验的方法研究了W、Mo、B在Cr7C3中多元复合添加的影响，研究结果表明复合添加B和W或Mo有利于Fe16Cr12C12结构的化学稳定性。

41、(Chemical Physics Letters, 2008)晶体结构稳定性对于相同结构的晶体结构，研究空位或原子置换对晶体导电性及结合键的强弱直接跃迁型能带结构间接跃迁型能带结构能隙的宽度结合键的共价性和离子性结合键的强度材料的磁性质与自旋态密度有关。晶体导电性及结合键的强弱直接跃迁型能带结构间接跃迁型能带结计算晶体结构的物理性质，计算力学稳定性及材料的韧脆性C11 (GPa)C12 (GPa)C44 (GPa)G (GPa)B (GPa)V4C4578.22147.17176.27518.50290.85V3FeC4638.79135.60130.95592.23301.60V3CrC4

42、677.78128.92191.09633.59310.62V3MoC4689.40115.43203.90653.73305.57V3NbC4643.42129.27150.84599.92301.13Nb4C4721.0294.34188.75699.39303.28Nb3FeC4635.02127.90136.46592.88297.41Nb3CrC4653.19123.42171.72614.22299.61Nb3MoC4724.36111.81183.33696.80318.69Nb3VC4699.14102.10180.08672.65300.53表6 金属原子M (M=Fe, C

43、r, Mo)置换后VC与NbC碳化物的三个独立的弹性常数及体弹模量根据判据，以上碳化物机械稳定。发生置换后，VC的体弹模量升高。Nb的碳化物中，Nb3MoC的体弹模量最高。固体的弹性常数与材料的力学性能密切相关，可以用来表征晶体对外界应力的响应，可作为原子间相互作用力大小的判据，并且弹性常数在材料的稳定性、硬度和强度上扮演着重要的角色。计算晶体结构的物理性质，计算力学稳定性及材料的韧脆性C11G (GPa)E (GPa)vAB/GV4C4191.04470.180.2310.8181.523V3FeC4170.10429.550.2630.5211.773V3CrC4220.84535.590

44、.2130.6961.407V3MoC4233.73558.730.1950.7111.307V3NbC4186.70464.170.2430.5871.613Nb4C4231.04552.760.1960.6021.313Nb3FeC4174.87438.640.2540.5381.701Nb3CrC4204.21499.210.2220.6481.467Nb3MoC4225.03546.470.2140.5991.416Nb3VC4220.33531.180.2050.6031.364表7 金属原子M (M=Fe, Cr, Mo)置换后VC与NbC碳化物的剪切模量、杨氏模量、各向异性因子及B/G比根据参考文献，B/G大于1.75，材料表现为韧性；反之则表现为脆性。由