ZnS电子结构的第一性原理.doc

______________________________________________________________________________________________________________精品资料纯ZnS材料的模型建立阿衣扎提1251624一、软件介绍用MaterialsStudio软件中的CASTEP计算程序下完成纯ZnS材料的模型。CASTEP是一种从头算量子力学程序，它基于密度泛函理论，利用总能量平面波赝势方法，用赝势代替离子势，通过平面波基组展开电子波函数，采用局域密度近似(LocalDensityApproximation，LDA)或广义梯度近似(GeneralizedGradientApproximation，GGA)对电子-电子相互作用的交换关联能进行校正。CASTEP计算程序是目前较为准确的电子结构计算方法。CASTEP适用于计算周期性的晶体结构，对于非周期性晶体结构一般采用特定的部分作为周期性晶体结构，建立单位晶胞后再进行计算。故CASTEP软件一般采用超晶胞模型，在周期系统中进行计算。CASTEP软件可以对晶体结构进行几何优化，得到晶胞稳定时的结构参数，通过计算超晶胞的总能和各孤立原子的总能可以得出晶体的结合能，通过计算电子分布密度和键布居来了解电荷转移情况、原子间的成键情况等。另外，还可以计算晶体及其原子的光学光谱、态密度和分态密度，对其相应的发光机制和电子机制进行分析研究。总的来说，它可以实现：计算体系的总能；进行结构优化；执行动力学任务；在设置的温度和关联参数下，研究体系中原子的运动行为；计算周期体系的弹性常数；化学反应的过度态搜索等。除此之外，计算一些晶体的性质，如能带结构、态密度、聚居数分析、声子色散关系、声子态密度、光学性质、应力等。量子力学计算精确度高但计算密集。直到最近，表征固体和表面所需的扩展体系的量子力学模拟对大多数研究者来说才切实可行。然而，不断发展的计算机功能和算法的进步使这种计算越来越容易实现。二、软件的使用1．计算任务的设置在CASTEP软件中行任务设置，主要是通过Visualizer应用窗口中的工具条之一“Calculation”来进行。我们可以更改工具框中的相应选项，来配置诸如：“电子选项”“结构优化选项”和“电子和结构性质选项”等。这几个选项是我们在运用CASTEP计算研究中非常重要的几个技术参数。其中，“电子选项”是很多其它计算任务也要涉及的。在CASTEP中还有如动力学、结构优化、弹性常数、过渡态等计算的设置。在程序运行之前，从研究的问题出发，要将软件中关键的一些任务参数设置成符合计算需要的值，我们才能得到所期望的运算结果。在利用CASTEP做有关能量、动力学、结构优化、弹性常数、过渡态等计算时，必须对电子选项进行设置。在电子选项中主要有精度设置、交换-关联函数的设置、赝势的设置、截断能的设置K点的设置。2．结构优化任务的设置结构优化是CASTEP计算中经常要进行的计算任务，特别是想要计算所关注体系的各种性质的时候，必须首先进行结构优化的计算，在得到结构优化结果文件以后，才能进行性质的计算。所以，正确的设置结构优化的参数是非常重要的。在CASTEP软件中，有四个参数来控制结构优化的收敛参数，第一个是能量的收敛精度，单位为eV/atom，是体系中每个原子的能量值；第二个是作用在每个原子上的最大力收敛精度；第三个是最大应变收敛精度，单位为GPa；第四个是最大位移收敛精度，。这些收敛精度指的是两次迭代求解之间的差，只有当某次计算的值与上一次计算的值相比小于设置的值时，计算才停止。3．计算体系性质的设置在CASTEP中可以计算体系的性质，如能带结构、态密度、聚居数分析、声子色散关系、声子态密度、光学性质、应力等。在计算能带和态密度这两项的计算设置之前，必需先进行自洽计算得到基态能量，而结构优化能够做到这一点，所以要在计算能带和态密度之前对体系进行结构优化。理论模型的构造闪锌矿ZnS结构的空间群为F-43M(216)，晶格常数a=b=c=0.54093nm，α=β=γ=90°，体对角线l／4处为S原子，八个顶角和六个面心为Zn原子。其2×1×1超晶胞结构如下图所示，共16个原子(考虑到计算效率，所有计算均采用1×1×1的超晶胞模型)。图1.1ZnS的2×1×1超晶胞结构图，(黄色球代表O原子，灰色球代表Zn原子)首先运行MSModeling软件，根据其默认提示创建一个新文档，得到图1.1所示的界面：图1.1MSModeling软件打开界面接着点击File中NewProject，将Untitled的文件名改为ZnS，选中Projectexplorer的目录上显示的改后的文件名，右键单击，选中New-3DAtomisticDocument，接着在3DA