计算材料学实验

《计算材料学实验》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、实验四利用MaterialStudio分析研究晶体材料性能姓名：高峰学号：1309050307班级：材料物理1303日期：2016.3.29一、实验目的1、了解MaterialStudio(MS)软件中有关固体材料科学设计各个模块功能；2、掌握在MS软件MaterialsVisualizer子模块中创建晶体结构模型；3、掌握在MS材料计算软件中研究晶体材料性质的方法；4、掌握查看和分析晶体材料属性的方法。5、分析AlAs晶体的晶格常数、态密度、能带图等性质。二、实验原理及方法MaterialsStudio是世界上著名的计算材料软件，其拥有十分强大的材料计算和设计功能。该软件中的

2、MaterialsVisualizer子模块，可用于构建晶体模型。根据晶体学与本课程所学知识，每种晶体都有其特有的晶体结构排布，每个原子占据一定的原子位置，每种品体都具有相应的品格常数、对称性等品体学数据。MaterialsVisualizer子模块可给出晶体结构的直观模型，并可以分析出晶体的晶格参数、空间群、原子坐标的数据。CASTAP是特别为固体材料学而设计的一个现代的量子力学基本程序，其使用了密度泛函(DFT)平面波澳势方法，进行第一原理量子力学计算，以探索如半导体，陶瓷，金属，矿物和沸石等材料的晶体和表面性质。MaterialStudio使用组件对话框中的CASTAP选

3、项允许准备，启动，分析和监测CASTAP服役工作。CASTAP的任务：能量任务、几何优化任务、动力学任务、性质任务。三、实验仪器MaterialsStudio软件、计算机。四、实验步骤1•构建晶体模型(1)创建空的单胞：build--buildcrystal,选择合适的空间群和晶格常数。(2)添加原子：build—addatoms输入原子种类及原子坐标。若有原子在晶胞外，选择rebu订d。(3)构建原胞：build--symmetry-primitivecello2.晶体模型计算(1)结构优化：modules--CASTEP-calculation--setup--task--

4、geometryoptimization，计算精度根据计算的时间耍求及单胞中的原子个数，一般选为fine；函数选默认的GGA-PBE方法。(2)性质计算：在setup的同级对话框中选择properties,CASTEP可完成的性质计算及计算过程中的计算精度控制等都可以在本对话框中进行选择。然后点击run,开始进行计算。3.晶体模型分析根据Modules--CASTEP--analysis工具，可以分析岀已计算的晶体性质。五、实验结果分析l.AIAs晶体结构AlAs空间群是F-43m,它的晶格参数为5.6622埃，依照步骤对AlAs晶体计算。计算后优化的3D视图的晶胞以及其电子密

5、度等能面如下图所示：半导体（AlAs晶胞）的能态密度和能带结构图如下:CAST"Strurtur*CA5WOntrty"Energy(«Y)meGASTEPDwrnatyofSUtes从以上的图中可以看岀，上方区域的线条表示的是导带，下方区域的线条表示价带，中间空白区域是禁带。禁带宽度在半导体范围内，这表明AlAs晶体具有半导体的性质。在态密度曲线图中，虚线表示的是费米能级。能带图显示的是在布里渊区屮K矢量沿着高对称性方向上的电子能量依赖度。一个原子的一个原子轨道在一维周期条件下将产生一条能带，能带的带宽取决于这些原子轨道的在周期方向上的成键强度，强度越大，带宽越大，成键越弱,

6、带宽越小。介电函数与频率的关系曲线如下图所示：CASTEPOp^calPropertiesScissorsoperator=0eV,Instrumentalsmearing:0.5eVCalculationgeometry:Polarized,polanzatwndrectx>n:(1.0000.0000.000)LossfunctonCASTB>OpbcalPropertoesSdssorsoperator-0eV,Instrumentalsmeanng:O.SeYCalculabongeometry:Polanzed,polanzatondrecoon:(1.0000.00

7、00.000)能量变化、最大位移量、最大内力、最大应力的收敛变化规律如下图所示:CASTEPOptmizat»nConvergenceConvergeoze(log10).EnergyChange(eV/atom)Max.Displacement⑺Max.Force(eV/?)Max.Stress(GPa)CASTEPOptealPropertiesSossorsoperator=0eV,InstrumentalsmeaririQ:0.5eVCalaiabongeometry:Polariz