资源描述
,单击此处编辑母版标题样式,*,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,计算材料学主要内容：,密度泛函理论基础,计算模型的建立,晶体结构优化,表面吸附,电荷密度分布,化学反应中的过渡态搜索,气体扩散,X-ray,结构精修,教员：,郝绿原,hly 3600834,傅正平,fuzp,李震宇（理论）,zyli,张文华（理论）,whzhang,资料：, Studio,使用组件对话框中的,CASTAP,选项允许准备，启动，分析和监测,CASTAP,服役工作。,计算：允许选择计算选项（如基集，交换关联势和收敛判据），作业控制和文档控制。,分析：允许处理和演示,CASTAP,计算结果。这一工具提供加速整体直观化以及键结构图，态密度图形和光学性质图形。,CASTAP,的任务,CASTAP,计算是要进行的三个任务中的一个，即单个点的能量计算，几何优化或分子动力学。可提供这些计算中的每一个以便产生特定的物理性能。性质为一种附加的任务，允许重新开始已完成的计算以便产生最初没有提出的额外性能。,在,CASTAP,计算中有很多运行步骤，可分为如下几组：,* 结构定义：必须规定包含所感兴趣结构的周期性的,3D,模型文件，有大量方法规定一种结构：可使用构建晶体（,Build Crystal),或构建真空板,(Build Vacuum Stab),来构建，也可从已经存在的的结构文档中引入，还可修正已存在的结构。,注意：,CASTAP,仅能在,3D,周期模型文件基础上进行计算，必须构建超单胞，以便研究分子体系。,提示：,CASTAP,计算所需时间随原子数平方的增加而增加。因此，建议用最小的原胞来描述体系，可使用,BuildSymmetryPrimitive Cell,菜单选项来转换成原胞。,CASTAP,中选择一项任务,1,从模块面板（,Module Explorer),选择,CASTAPCalculation,。,2,选择设置表。,3,从任务列表中选择所要求的任务。,*,计算设置：合适的,3D,模型文件一旦确定，必须选择计算类型和相关参数，例如，对于动力学计算必须确定系综和参数，包括温度，时间步长和步数。选择运行计算的磁盘并开始,CASTAP,作业。,* 结果分析：计算完成后，相关于,CASTAP,作业的文档返回用户，在项目面板适当位置显示。这些文档的一些进一步处理要求获得可观察量如光学性质。,CASTAP,能量任务,CASTAP,能量任务允许计算特定体系的总能量以及物理性质。,除了总能量之外，在计算之后还可报告作用于原子上的力；也能创建电荷密度文件；利用材料观测仪（,Material Visualizer),允许目测电荷密度的立体分布；还能报告计算中使用的,Monkhorst-Park,的,k,点的电子能量，因此在,CASTAP,分析中可生成态密度图。,对于能够得到可靠结构信息的体系的电子性质的研究，能量任务是有用的。只要给定应力性质，也可用于计算没有内部自由度的高对称性体系的状态方程（即压力,-,体积，能量,-,体积关系）。,注意：具有内部自由度的体系中，利用几何优化（,Geometry Optimization),任务可获得状态方程。,CASTAP,中能量的默认单位是电子伏特,(eV),，各种能量单位的换算关系见,Mohr.P.J(2000).,1 eV=0.036749308 Ha=23.0605 kcal/mole=96.4853 kJ/mole,CASTAP,几何优化任务,CASTAP,几何优化任务允许改善结构的几何，获得稳定结构或多晶型物。通过一个迭代过程来完成这项任务，迭代过程中调整原子坐标和晶胞参数使结构的总能量最小化。,CASTAP,几何优化是基于减小计算力和应力的数量级，直到小于规定的收敛误差。也可能给定外部应力张量来对拉应力、压应力和切应力等作用下的体系行为模型化。在这些情况下,反,复迭,代内部应力张量直到与所施加的外部应力相等。,几何优化处理产生的模型结构与真实结构紧密相似。利用,CASTAP,计算的晶格参数精度列于右图。,状态方程计算,在所施加静压力下几何优化可用于确定材料的体模量,B,和对压力的导数,B,=,dB/dP,。,过程包括计算理论状态方程（,EOS,），该方程描述单胞体积与外部静压力的关系。工艺非常类似于真实实验：使用几何优化对话框中的应力列表将外部压力固定。通过进行几何优化可以找到在此压力下的单胞体积。随后的,P-V,数据分析与实验研究精确一致。描述,EOS,选择分析表达式，其参数适于计算数据点。最流行的,EOS,形式是三阶,Birch-Murnaghan,方程：,式中,V,0,为平衡体积。,Cohen,等进行了,EOS,各种解析式的的详细比较研究。,注意：从相应实验中获得的,B,和,B,值依赖于计算使用的压力值范围。利用金刚石压砧获得的实验值通常在,0-30GPa,范围内，因此推荐理论研究也在这个范围内。在研究中避免使用负压力值也很重要。此外，用于生成,P-V,数据序列的压力值可能是不均匀的，在低压力范围要求更精确采样以便获得体模量精确值。,P-V,几何优化方法,在默认条件下，,CASTAP,使用,BFGS,几何优化方法。该方法通常提供了寻找最低能量结构的最快途径，这是支持,CASTAP,单胞优化的唯一模式。,衰减分子动力学（,Damped molecular dynamics,）方法是另一种可以选择的方法，该方法对具有平滑势能表面的体系如,分子晶体或表面分子,与,BFGS,同样有效。,CASTAP,动力学任务,CASTAP,动力学任务允许模拟结构中原子在计算力的影响下将如何移动。,在进行,CASTAP,动力学计算以前，可以选择热力学系综和相应参数，定义模拟时间和模拟温度。,选择热力学系综,对牛顿运动定律积分允许探索体系恒值能量表面（,NVE,动力学）。然而，在体系与环境进行热交换条件下发生最本质的现象。使用,NVT,系综（或者是确定性的,Nos,系综或者是随机性的,Langevin,系综）可模拟该条件。,定义时间步长（,timestep,）,在积分算法中重要参数是时间步长。为更好利用计算时间，应使用大的时间步长。然而，如果时间步长过大，则可导致积分过程的不稳定和不精确。典型地，这表示为运动常数的系统偏差。,注意：量子力学分子动力学计算要求比力场动力学使用更小的时间步长。,动力学过程的约束,CASTAP,支持,Langevin NVT,或,NVE,动力学过程的线性约束。然而，借助,Material Studio,界面可以近似使用以下两种更基本的约束： 质心固定，单个原子固定。,使用,seedname,.cell,文档可以利用更复杂的约束。,CASTAP,性质任务,CASTAP,性质任务允许在完成能量，几何优化或动力学运行之后求出电子和结构性质。可以产生的性质如下：,* 态密度（,DOS,）：利用原始模拟中产生的电荷密度和势能，非自恰计算价带和导带的精细,Monkhorst-Pack,网格上的电子本征值。,* 带结构：利用原始模拟中产生的电荷密度和势能，非自恰计算价带和导带的布里渊区高对称性方向电子本征值。,* 光学性质：计算电子能带间转变的矩阵元素。,CASTAP,分析对话可用于生成包含可以测得的光学性质的网格和图形文件。,* 布局数分析：进行,Mulliken,分析。计算决定原子电荷的键总数和角动量（以及自旋极化计算所需的磁矩）。任旋地，可产生态密度微分计算所要求的分量。,* 应力：计算应力张量，并写入,seedname.castep,文档。,如果要进行单胞参数固定时进行几何优化运行和要检查点阵偏离平衡的程度，这些信息是有用的。例如，可进行符合于给定体系理论基态的固定单胞的点缺陷的超晶胞研究。几何优化后的应力值显示了与超单胞近似相关联的弹性效应。,注意：为计算某种性质，从适当模拟得到的结果文档必须以当前的文件夹形式出现。,目的,:,介绍,CASTEP,中的结构优化， 使用立体可视化工具显示等值面模块,: Materials Visualizer, CASTEP,前提,:,使用晶体建模工具,用第一原理预测,AlAs,的晶格参数,内容,1.,构建,AlAs,的晶体结构,2.,设置并进行,CASTEP,计算,3.,分析结果,4.,比较计算的结构参数和实验数据,(1),图示电荷密度,(2),图示态密度和带结构,引言,本指南介绍了,CASTEP,是如何使用量子力学方法来确定材料的晶体结构，使用者将学会如何构建晶体结构，设定一个,CASTEP,几何优化任务，然后分析计算结果。,背景,密度泛函理论,(DFT),在周期性大尺度材料上应用的进展，对材料设计和加工越来越重要。该理论使得研究者能对实验数据进行解释；并从未知晶体的结构性质、结合能和表面分子的活动性确定材料的本征性质。这些理论工具可用于指导新材料的设计，帮助研究者了解内在的化学和物理过程。,注意,:,如果你的服务器没有足够快的,CPU,，请慎用,CASTEP,进行几何优化计算，因为它会占用相当长的时间 。,Al,的分数坐标：,(0 0 0,）,(1/2 1/2 0) (1/2 0 1/2),（,0 1/2 1/2,）,As,的分数坐标：,(3/4 3/4 3/4,）,(1/4 1/4 3/4) (1/4 3/4 1/4),（,3/4 1/4 1/4,）,(0 0 0,）,(1/2 0 1/2),(1/2 1/2 0),（,0 1/2 1/2,）,(1/4 1/4 3/4),（,3/4 1/4 1/4,）,（,1/4 3/4 1/4,）,（,3/4 3/4 3/4,）,As: (3/4 3/4 3/4,）,= (1/4 1/4 1/4,）,1.,构建,AlAs,的晶体结构,空间群是,F-43m,构建一个晶体结构，需要知道该晶体的空间群、晶格参数和晶体的内坐标。对,AlAs,来说，空间群是,F-43m,，空间群代号为,216,。原胞有两个原子，,Al,和,As,的分数坐标分别为,(0, 0, 0),和,(0.25, 0.25, 0.25),，晶格参数为,5.6622,.,。,第一步是建立晶格。,在,D disk,上建立英文目录,D:class3,。按下面步骤,在,Project Explorer,内，建立,AlAs,根目录。,从“开始”或快捷图标 打开,MS,。,找到,class3,，按“打开”按钮,输入,AlAs,，这将是新的,Project,的名字。,在,Project Explorer,中，右击根目录,AlAs,，选择,New | 3D Atomistic Document,。,右击,3D Atomistic document,文件，将该文件重新命名为,AlAs.xsd,。,从菜单栏里选择,Build / Crystals / Build Crystal,。,Build Crystal,对话框显示出来。,点击,Enter group,输入,216,，按下,TAB,按钮,(,或,在,Enter group,中选择,F-43m),，,空间群信息更新为,F-43m,空间群。,空间群信息框中的信息也随着,F-43m,空间群的信息而发生变化 。,杨碚芳课,选择,Lattice Parameters,标签，把,a,值从,10.00,变为,5.662,。点击,Build,按钮。,单击此图标，然后可旋转晶格，显示其立体结构。,一个没有原子的,3D,格子显示在,3D Atomistic,文件里。,AsAl?,Atom # OX SITE x y z SOF H,Al 1 +3 4 a 0 0 0 1. 0,As 2 -3 4 c 0.25 0.25 0.25 1. 0,*end for ICSD #67784,Al,的分数坐标：,(0 0 0,）,(1/2 1/2 0) (1/2 0 1/2),（,0 1/2 1/2,）,As,的分数坐标：,(3/4 3/4 3/4,）,(1/4 1/4 3/4) (1/4 3/4 1/4),（,3/4 1/4 1/4,）,输入几个原子？,Ba,3,Si,6,O,9,N,4,作业,1,：解释符号和群,F-43m,的意思，给出参考书。,需输入,3,个,Ba,2,个,Si,2,个,N,3,个,O,*,从菜单栏中选择,Build / Add Atoms,。通过,Add Atoms,对话框，我们可以把原子添加到指定的位置，其对话框如下：,在,Add Atoms,对话框中选择,Options,标签，确定,Coordinate system,为,Fractional,。,*,如上所示，选择,Atoms,标签，通过周期表，在,Element,文本框中输入,Al,， 再输入,Al,的分数坐标,(0, 0, 0),，,然后按下,Add,按钮，铝原子就添加到结构中了。,*,从菜单栏中选择,View / Display Style,。在打开的对话框中，选择,Ball and stick,显示方式。,*,同前，在,Element,文本框中键入,As,。在,a, b, c,文本框中键入,0.25,。按,Add,按钮，这样,As,也加入到晶格中。关闭对话框。,说明：上面操作虽然只加入一个,Al,、一个,As,，但群的对称操作在晶体中补充了剩余的,Al,原子和,As,原子,(,等位原子,),。,如果新加原子以,line,方式显示，可重复前面步骤，改显示方式为,Ball and Stick,。,单击此图标，出现添加原子,Add Atoms,对话框。选择原子名称，输入分数坐标，按,Add,，则原子添加到晶体结构图中。重复操作，直到添加完晶胞中的所有原子。关闭,Add Atoms,框。,前面的添加原子操作也可用下面图标来实现。这里不再重复。,*,从上面的,AlAs,晶体结构看出，近邻晶胞中的原子也显示出来。这种显示表示了,AlAs,晶体中键的拓扑结构。当然，可以通过重新建造晶体结构来移去这些,近邻晶胞中的,原子。,从菜单栏中选择,Build / Crystals / Rebuild Crystal.,，打开对话框，按下,Rebuild,按钮。在显示出的晶体结构中那些原子就被移走了。已经把显示方式定为,Ball and Stick,。,按下面图示步骤，保存结果。,3D Viewer,内的晶体结构是传统的,晶胞,conventional (unit) cell,，显示了晶格的立方对称性。,CASTEP,利用了晶格的对称性，可以使用只包含两个原子的,原胞,primitive (unit) cell,来进行计算，而晶胞包含了,8,个原子。电荷密度、键长和每个原子的总能量将是一样的，而不管这个,unit cell,是如何被定义的。这样，使用原胞，原子数较少，计算量大大减小，计算时间将被缩短。,Note,：在计算磁性体系时，使用了自旋极化。这时要注意，电荷密度自旋波的周期是原胞的数倍。,*,选择菜单栏里的,Build / Symmetry / Primitive Cell,。,模型文件,(3D Viewer),显示为原胞,(primitive cell),。,AlAs,的原胞,不同角度,在晶体图上按右键，选,Label,，在出现的对话框中选,ElementSymbol,。,按,apply,，晶胞上显示元素符号。,2.,设置,CASTEP,计算任务,从工具栏中选择,CASTEP,工具，再选择,Calculation,或从菜单栏中选择,Modules | CASTEP |Calculation,。,CASTEP Calculation,对话框如下：,下面我们分两步，先优化,AlAs,的几何结构，再计算,Band Structure,和,Density of states,。,(1),优化,AlAs,晶体结构,* 把,Task,改为,Geometry Optimization,，计算精度,Quality,设置为,coarse,。,*,结构优化的默认设置是优化原子坐标。在本例中，我们不仅要优化原子坐标，同时也要优化晶格常数。,* 按下与,Task,相关的,More,按钮，勾选上,Optimize Cell,， 关闭此对话框。,当改变计算精度的时候，其它的参数也会自动作相应的变化。,*,选择,Electronic,标签栏，按下,More,按钮。在,SCF,对话框里作如下设置，将,Charge,由,0.5,改为,0.15,，,钩上,Fix occupancy,。,*,选择,Properties,标签栏，里面的计算任务都,不,要选。,*,选择,Job Control,标签栏，按下,More,按钮。在,CASTEP Job Control Options,对话框里，改变,Update interval,为,30.0 s,，关闭此对话框。,按下,Run,按钮，关闭对话框。,注意，此时,3D,结构为激活窗口。如果激活窗口是文本，,则,Run,为灰色。,几秒钟后，一个新文件夹出现在,Project Explorer,内，该文件夹包含了所有的计算结果。,如果使用客户端服务器模式，当工作结束时，文件会被传回到客户端。数据传输过程需要一定的时间，与文件的大小有关。,Job Explorer,显示了所有正在运行的工作的状态。它显示了很多有用的信息，包括服务器和工作代码。如果需要，也可以通过,Job Explorer,来中止运行工作。,在工作运行过程中，四个文件打开了。这些文件包含了晶体结构、结构优化过程中模型的更新、工作参数的设置和运行状态的信息，以及关于总能量、能量变化、应力、压力和位移随迭代次数变化的图表。,计算过程中出现的两个表示能量收敛的图框。,查看计算设置,几分钟后，计算结束，出现,Job Completed,提示，表示计算成功。,输出文本文档为,AlAs.castep,，包含优化信息，在,AlAs CASTEP GeomOpt,文件夹中。,按下面图示操作，关闭工作窗口中的文件。,(2),接下来，利用优化过的,AlAs,结构，计算,AlAs,的,Band Structure,和,Density of states,。,*,在结构优化的文件中，双击,AlAs.xsd,、,AlAs.castep,两个文件。这两个文件出现在工作窗口中。按下面图示设定计算任务，直至,Run,。,出现新的文件夹,*,打开新文件夹,AlAs CASTEP Properties,，双击,AlAs_BandStr.castep,，这此文件出现在工作窗口中。,*,按图示操作，显示,AlAs,的,Band Structure,。,(3),计算结束后，查看,AlAs,的,Band Structure,和,Density of states,的计算结果。,可用工具 放大、缩小能带图。,能带图下方的字母表示布里渊区的高对称点。单击,AlAs.xsd,文件，使其为当前活动窗口。再按图示点击计算设置 ，使当前状态为,Band Structure,。点击下方的,More,，则显示出对称点的坐标。,1,2,*,在新文件夹,AlAs CASTEP Properties,中双击,AlAs_DOS.castep,，这此文件出现在工作窗口中。,*,按图示操作，显示,AlAs,的,DOS,。,可用工具 放大、缩小能带图。,(4),前面的结构优化显示，,AlAs,的总能量随迭代次数振荡。下面改计算方法进行结构优化，避免振荡。计算步骤与前面类似，图示如下。,结构优化过程，可见随着迭代次数的增加，,AlAs,的总能量逐渐减小。,*,按下面图示步骤，保存计算结果，关闭工作窗口中的文件。,(5),同样，利用优化过的,AlAs,结构，计算,AlAs,的,Band Structure,和,Density of states,。,*,在新出现的文件夹中，打开,AlAs.castep,、,AlAs.xsd,两个文件。这两个文件都出现在工作窗口中。,作业,2,：在,.CASTEP,文件中找到并记下总电子能量，与前面优化后的总电子能量比较，哪个应该小？,*,设置计算任务，按,Run,进行计算，关闭对话框。,(6),计算结束后，查看,AlAs,的,Band Structure,和,Density of states,的计算结果。,*,计算结束后，出现提示，关闭。,*,保存计算结果，清理工作窗口。,*,在新出现的,AlAs CASTEP Properties,文件夹中，双击,AlAs_BandStr.castep,，该文件出现在工作窗口中。,*,按图示操作，显示带结构。,*,同样，双击,AlAs_DOS.castep,，显示态密度。,说明：,分析工具可以用来显示态密度,(DOS),和能带结构。,能带结构图显示了布里渊区内沿着高对称方向电子能量对,k,矢的依赖性。这些图提供了一个对材料的电子结构进行定性分析的非常有用的工具。譬如，与近自由的,s,、,p,电子构成的能带相比，很容易鉴别出,d,、,f,电子构成的窄带。,DOS,和,PDOS,图给出了材料的电子结构的一个快速定性图像，有时候它们可以直接和实验光谱结果相关联。,CASTEP,的主要输出结果文件,AlAs.castep,包含了有限的能带结构和,DOS,信息，更多的详细信息包含在,AlAs_BandStr.castep,文件内。,打开,Analysis,对话框，选上,Band structure,。从这个对话框可以看出，可以把能带结构和态密度信息显示在同一个图中。在,DOS,部分，选上,Show DOS,，单击,View,，出现的图包含了带结构和,DOS,两种信息。当然，可以分别显示能带结构和态密度。,可以按图片、数据格式输出图文件，数据可由,Excel,等软件读取。,还可以借助,CASTEP,来计算很多其他性质，比如反射率和介电函数等等。,DFT,计算带隙,Eg,，数值偏小。比较下图可知，计算精度高，,Eg,大。,(7),比较两次计算的结果,1,2,4,3,1.296,1.762,作业,2,：比较两次计算出的总电子能量,3,分析结果,当结果文件被传输回来，会得到下列数个文件：,AlAs.xsd,最后的优化结构,AlAs Trajectory.xtd -,一个轨迹文件，,包含了,每一个优化步骤后,的结构,AlAs.castep,包含了优化信息的输出,文本文件,AlAs.param,模拟所用的输入参数,计算任何一个性质，都会产生,.param,和,.castep,文件。,在,AlAs,结构中，由于对称性的存在，受力为,0,，但是应力的大小取决于晶格参数。这样，,CASTEP,就会努力去最小化系统的总能量和应力。因此，为保证计算能够合适地完成，检查压力收敛是非常重要的。,在,Project Explorer,内，双击,AlAs.castep,，将其激活为当前工作文件。选择菜单栏里的,Edit | Find.,，在文本框中输入“,completed successfully”,，按下,Find Next,按钮，,AlAs.castep,文件滚动。,看到一个含有两行的表格，最后一列的每一行都显示为,Yes,，这表明计算成功地结束。,开始创建晶胞时，就知道晶格长度为,5.6622 ,。因此，可以把能量最小化后的晶格长度与初始的实验数据相比较。实验的晶格长度基于晶胞结构，而不是原胞，因此需要将现在的原胞转化为晶胞，再与实验数据比较。,4.,比较,AlAs,晶体结构的计算数据与实验数据,*,双击,AlAs.xsd,使其为当前工作文件,*,从菜单栏里选择,Build / Symmetry / Conventional Cell,，晶胞显示出来。,*,有数种方法看到晶格长度，一种就是打开,Lattice Parameters,对话框。在模型文件,(3D Viewer),上右击，选择,Lattice Parameters,。格矢大约为,5.721128 ,。,*,另一种简单的方法是在左侧,Properties,中选择,Lattice 3D,，其中显示晶格常数为,5.72113 ,。,误差大约是,-0.5%,。这在,1-2%,典型误差范围内，这个误差值是赝势平面波方法与实验结果比较的期望误差。使用,GGA,函数，晶格常数的计算值偏大；使用,LDA,函数，晶格常数的计算值偏小；,继续之前，需要保存工作，并关闭所有窗口。,选择菜单栏上的,File | Save Project,，然后是,Window | Close All,。,*,从工具栏选择,CASTEP,，然后选择,Analysis(,或从菜单栏选择,Modules / CASTEP / Analysis),，再选中,Electron density,选项。,此时，会有一条信息,“,no results file is available,”,，说没有什么结果文件，所以需要指定结果文件。,在,Project Explorer,内，双击,AlAs.castep,。这将把结果文件和分析对话框关联起来。,5.,可视化电荷密度,可以用,CASTEP,分析工具得到电荷密度。,*,需要指定一个,3D,文件来显示等密度面。,在,Project Explorer,内，双击优化后的,AlAs.xsd,文件。,*,选择菜单栏里的,Build / Symmetry /,Primitive Cell,* CASTEP Analysis,对话框上的,Import,按钮现在是激活状态。,按下,Import,按钮。,AlAs,的电子密度等值面,作业,3,：电子密度有无负值？电子如何分布？,可以通过,Display Style,对话框来改变等密度面的设置。,*在等值面的,3D,文件上右击，选择,Display Style,，选择,Isosurface,标签栏。,Isosurface,标签栏显示出来：,*,这里可以改变不同的设置。,在,Isovalue,文本框里，输入,0.1,，然后按,TAB,键。,注意等密度面是如何改变的。,*,把,Transparency,滑条向右移动。,向右移动,Transparency,滑条的时候，等密度面变得越来越透明。,在,3D,文件上按住鼠标右键，滚动模型。,当模型滚动的时候，等密度面变成点状显示以提高滚动速度。,*,通过,Display Style,对话框可以移走等密度面。,取消选择,Visible,选项，关闭,Display Style,对话框。,可以在任何时候通过勾选上,Isosurface checkbox,来显示等密度面。,用,CPK,方式显示晶体结构，清楚地看到原子。,调整电子云的透明度，清楚地显示电子云与原子的相对位置。,破缺，积分小。,里面的球面,s,电子，弥散,.,As,As,Al,结构优化改变原子位置和晶格参数，得到一种电子密度分布，使得电子的总能量最小。,电子密度是电子电荷乘上波函数的平方。图中蓝色表示等值面的外面，灰色表示等值面的里面。由等值面的颜色、形状,(,包裹,As),、数值可以看出，电荷集中在,As,上离子晶体。,This is the end of the tutorial.,
展开阅读全文