【www.arisingsemi.com--软件制图】

电子构型
纯ZnS材料的模型建立
阿衣扎提                      1251624
一、软件介绍
用Materials Studio软件中的CASTEP计算程序下完成纯ZnS材料的模型。CASTEP是一种从头算量子力学程序,它基于密度泛函理论,利用总能量平面波赝势方法,用赝势代替离子势,通过平面波基组展开电子波函数,采用局域密度近似(Local Density Approximation,LDA)或广义梯度近似(Generalized Gradient Approximation,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.1 ZnS的2×1×1超晶胞结构图,(黄色球代表O原子,灰色球代表Zn原子)
首先运行MS Modeling软件,根据其默认提示创建一个新文档,得到图1.1所示的界面:

接着点击File中New Project ,将Untitled的文件名改为ZnS,选中Project explorer的目录上显示的改后的文件名,右键单击,选中New-3D Atomistic Document,接着在3D Atomistic Document右键单击,将其更名(rename)为ZnS。(注意此次改的文件名必须与上一步骤中的文件名保持一致,只要保持一致即可,文件名可随意取),从菜单中选择Build-Crystals-Build Crystal,屏幕显示出Build Crystal的对话框,在Enter group中选择216F43-M,如下图1.2所示;接着选择Lattice Parameters标签,把a的数值从10.000改为5.405,如图1.3所示:


得到ZnS晶胞后就可以添加原子了。
从菜单栏中选择Build | Add Atoms。通过它我们可以把原子添加到指定的位置,在Add Atoms对话框中选择Options标签,确定Coordinate system为Fractional。如上所示。选择Atoms标签,在Element文本框中键入S在a,b,c文本框中分别键入0.250,0.250,0.250,如图1.4所示。然后按下Add按钮,硫原子就添加到结构中了。同理在Element文本框中键入Zn,(也可以调出元素周期表选中Zn元素)在a,b,c文本框中分别键入0.000,0.000,0.000,如图1.5,按Add按钮,锌原子就添加到结构中了,关闭对话框,便得到了图1.6所示的ZnS电子结构图。
当然,也可以通过重新建造晶体结构来移去这些原子。从菜单栏中选择Build | Crystals | Rebuild Crystal...,按下Rebuild按钮即可以重建ZnS的电子结构。


把显示方式变为Ball and Stick。
在模型文档中右键单击,选择Display Styles,按下Ball and stick按钮。
关闭对话框。建立了如图1.7的结构。在计算的过程中始终坚持实际计算结合理论分析,使理论和分析不脱节。

2.计算方法
我们的计算采用基于密度泛函理论(density functional theory,DFT)的广义梯度近似(generalized gradient approximation,GGA)修正方法。计算中不考虑自旋影响;采用超软赝势(ultra-soft pseudo-potential,USSP)近似描述价电子与离子实之间的相互作用;交换相干函数用GGA/PBE进行处理,在倒易的K空间中,通过平面波截止能量的选择可以提高计算精度,同时运用快速傅立叶变化技术,使能量、相互作用势等计算在实空间和倒易空间之间快速转换。
能量计算应用Pulay密度混合法,设置自治精度为5.0×10-6eV/atom,平面波截断能Ecut取为350eV,作用在每个原子上的力不大于0.01eV/ Å,内应力不大于0.02GPa,布里渊区的积分计算采用5×5×5的Monkhorst-pack[20]特殊k点对全布里渊区求和,能量计算都在倒易空间进行。同时为了得到稳定精确的计算结果,我们先对超晶胞进行了几何结构优化,得到晶胞参数后,再优化其内坐标,得到稳定的结构之后,再进行单点能的计算。
采用第一性原理的方法研究了本征ZnS电子结构的一些性质。计算了ZnS系统的能带结构、几何参数和电子态密度。结果表明,ZnS为直接禁带半导体材料,其带隙为3.68eV。导带部分主要分布在4.0~8.7eV,主要来源于锌4s态的贡献;价带可以分为处于-6.8~-4.9eV能量范围的下价带区和处于-4.9~0.0eV能量范围的上价带区,下价带主要来源于高度局域的锌3d态的贡献,上价带主要来源于硫3p态的贡献,在-4.0eV附近,锌4s态与硫3p态高度杂化。通过Mulliken布局分析得到ZnS晶体中Zn原子失去电子,为电子的给与体,S原子得到电子,是电子受主,且Zn原子与S原子形成的是共价键。
通过这些计算结果可以对实验中遇到的一些问题和现象从微观结构方面进行解释,并通过理论计算结果为实验提供理论依据,同时对制备高质量的ZnS材料也具有重要的指导意义。

本文来源:http://www.arisingsemi.com/it/60912/