CASTEP概述及一PPt详细教程.ppt

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,CASTEP,概述,CASTEP,是特别为固体材料学而设计的一个现代的量子力学基本程序，其使用了密度泛函,(DFT),平面波赝势方法，进行第一原理量子力学计算，以探索如半导体，陶瓷，金属，矿物和沸石等材料的晶体和表面性质典型的应用包括表面化学、键结构、态密度和光学性质等研究，,CASTEP,也可用于研究体系的电荷密度和波函数的,3D,形式此外，,CASTEP,可用于有效研究点缺陷,(,空位，间隙和置换杂质,),和扩展缺陷,(,如晶界和位错,),的性质Material Studio,使用组件对话框中的,CASTEP,选项允许准备、启动、分析和监测,CASTEP,服役工作计算：允许选择计算选项,(,如基集，交换关联势和收敛判据,),，作业控制和文档控制分析：允许处理和演示,CASTEP,计算结果这一工具提供加速整体直观化以及键结构图、态密度图形和光学性质图形关于,CASTEP,CASTEP,的任务,CASTEP,计算是要进行的三个任务中的一个，即单个点的能量计算，几何优化或分子动力学可提供这些计算中的每一个以便产生特定的物理性能。

性质为一种附加的任务，允许重新开始已完成的计算以便产生最初没有提出的额外性能在,CASTEP,计算中有很多运行步骤，可分为如下几组：,结构定义：必须规定包含所感兴趣结构的周期性的,3D,模型文件，有大量方法规定一种结构：可使用构建晶体,(Build Crystal),或构建真空板,(Build Vacuum Stab),来构建，也可从已经存在的的结构文档中引入，还可修正已存在的结构注意：,CASTEP,仅,能在,3D,周期模型文件基础上进行计算，必须构建超单胞，以便研究分子体系提示：,CASTEP,计算所需时间随原子数平方的增加而增加因此，建议是用最小的初晶胞来描述体系，可使用,BuildSymmetryPrimitive Cell,菜单选项来转换成初晶胞计算设置：合适的,3D,模型文件一旦确定，必须选择计算类型和相关参数，例如，对于动力学计算必须确定系综和参数，包括温度，时间步长和步数选择运行计算的磁盘并开始,CASTEP,作业结果分析：计算完成后，相关于,CASTEP,作业的文档返回用户，在项目面板适当位置显示这些文档的一些进一步处理要求获得可观察量如光学性质CASTEP,中选择一项任务,1,从模块面板,(Module Explorer),选择,CASTEPCalculation,。

2,选择设置表3,从任务列表中选择所要求的任务CASTEP,能量任务,CASTEP,能量任务允许计算特定体系的总能量以及物理性质除了总能量之外，在计算之后还可报告作用于原子上的力；也能创建电荷密度文件；利用材料观测仪,(Material,Visualizer,),允许目测电荷密度的立体分布；还能报告计算中使用的,Monkhorst,-Park,的,k,点的电子能量，因此在,CASTEP,分析中可生成态密度图对于能够得到可靠结构信息的体系的电子性质的研究，能量任务是有用的只要给定应力性质，也可用于计算没有内部自由度的高对称性体系的状态方程,(,即压力,-,体积，能量,-,体积关系,),注意：具有内部自由度的体系中，利用几何优化,(Geometry Optimization),任务可获得状态方程CASTEP,中能量的默认单位是电子伏特,(,eV,),，,各种能量单位的换算关系见,Mohr.P.J(2000).,1,eV,=0.036749308 Ha=23.0605 kcal/mole=96.4853 kJ/mole,CASTEP,几何优化任务,CASTEP,几何优化任务允许改善结构的几何，获得稳定结构或多晶型物。

通过一个迭代过程来完成这项任务，迭代过程中调整原子坐标和晶胞参数使结构的总能量最小化CASTEP,几何优化是基于减小计算力和应力的数量级，直到小于规定的收敛误差也可能给定外部应力张量来对拉应力，压应力和切应力等作用下的体系行为模型化在这些情况下反复迭代内部,几何优化处理产生的模型结构与真实结构紧密相似利用,CASTEP,计算的晶格参数精度列于右图应力张量直到与所施加的外部应力相等状态方程计算,在所,施加静压力下几何优化可用于确定材料的体模量,B,和对压力的导数,B,=,dB/,dP,过程包括计算理论状态方程,(EOS),，该方程描述单胞体积于外部静压力的关系工艺非常类似于真实实验：使用几何优化对话框中的应力列表将外部压力固定通过进行几何优化可以找到在此压力下的单胞体积随后的,P-V,数据分析与实验研究精确一致描述,EOS,选择分析表达式，其参数适于计算数据点最流行的,EOS,形式是三阶,Birch-,Murnaghan,方程：,式,中,V,0,为平衡体积Cohen,等,进行了,EOS,各种解析式的的详细比较研究注意：从相应实验中获得的,B,和,B,值,依赖于计算使用的压力值范围利用金刚石压砧获得的实验值通常在,0-30GPa,范围内，因此推荐理论研究也在这个范围内。

在研究中避免使用负压力值也很重要此外，用于生成,P-V,数据序列的压力值可能是不均匀的，在低压力范围要求更精确采样以便获得体模量精确值P-V,几何优化方法,在缺损条件下，,CASTEP,使用,BFGS,几何优化方法该方法通常提供了寻找最低能量结构的最快途径，这是支持,CASTEP,单胞优化的唯一模式衰减分子动力学,(,Damped molecular dynamics,),方法是另一种可以选择的方法，该方法对具有平滑势能表面的体系如分子晶体或表面分子与,BFGS,同样有效CASTEP,动力学任务,CASTEP,动力学任务允许模拟结构中原子在计算力的影响下将如何移动在进行,CASTEP,动力学计算以前，可以选择热力学系综和相应参数，定义模拟时间和模拟温度选择热力学系综,对,牛顿运动定律积分允许探索体系恒值能量表面,(NVE,动力学,),然而，在体系与环境进行热交换条件下发生最本质的现象使用,NVT,系综,(,或者是确定性的,Nos,系综或者是随机性的,Langevin,系综,),可模拟该条件定义时间步长,(,timestep,),在积分算法中重要参数是时间步长为更好利用计算时间，应使用大的时间步长。

然而，如果时间步长过大，则可导致积分过程的不稳定和不精确典型地，这表示为运动常数的系统偏差注意：量子力学分子动力学计算要求比力场动力学使用更小的时间步长动力学过程的约束,CASTEP,支持,Langevin,NVT,或,NVE,动力学过程的线性约束然而，借助,Material Studio,界面可以近似使用以下两种更基本的约束：质心固定，单个原子固定,使用,seedname,.cell,文档可以利用更复杂的约束CASTEP,性质任务,CASTEP,性质任务允许在完成能量，几何优化或动力学运行之后求出电子和结构性质可以产生的性质如下：,态密度,(DOS),：利用原始模拟中产生的电荷密度和势能，非自恰计算价带和导带的精细,Monkhorst,-Pack,网格上的电子本征值带结构：利用原始模拟中产生的电荷密度和势能，非自恰计算价带和导带的布里渊区高对称性方向电子本征值光学性质：计算电子能带间转变的矩阵元素CASTEP,分析对话可用于生成包含可以测得的光学性质的网格和图形文件布局数分析：进行,Mulliken,分析计算决定原子电荷的键总数和角动量,(,以及自旋极化计算所需的磁矩,),任旋地，可产生态密度微分计算所要求的分量。

应力：计算应力张量，并写入,seedname.castep,文档如果要进行单胞参数固定时进行几何优化运行和要检查点阵偏离平衡的程度，这些信息是有用的例如，可进行符合于给定体系理论基态的固定单胞的点缺陷的超晶胞研究几何优化后的应力值显示了与超单胞近似相关联的弹性效应注意：为计算某种性质，从适当模拟得到的结果文档必须以当前的文件夹形式出现用第一原理预测,AlAs,的晶格参数,本指南主要是阐明在,Materials Studio,当中如何运用量子力学来测定物质的晶体结构你将从中学到如何构建晶体结构以及如何设置,CASTEP,几何优化运行和分析结果本指南的内容如下：,1,构建,AlAs,的晶体结构,2,设置和运行,CASTEP,中的计算,3,分析结果,4,比较实验数据和结构,注意,如果你的服务器没有足够快的,CPU,，,本指南限制使用,CASTEP,进行,几何优化计算，因为它会占用相当长的时间1,构建,AlAs,的晶体结构,为了构建晶体结构，我们需要知道你想要构建的晶体的空间群信息，晶格参数以及它的内部坐标以,AlAs,为例，它的空间群是,F-43m,或空间群数字是,216,它有两种基本元素,Al,和,As,，,其分数坐标分别为,(0 0 0),和,(0.25 0.25,0.25,),。

它的晶格参数为,5.6622,埃第一步是构建晶格在,Project explorer,的跟目录上右键单击，选中,New|3D Atomistic Document,接着在,3D Atomistic Document,右键单击，把它更名为,AlAs,从菜单中选择,Build|Crystals|Build Crystal,，,然后显示出,Build Crystal,对话框，如下：,在,Enter group,中选择,F-43m,或在,Enter group,中单击，然后键入,216,，再按下,TAB,键,.,(,空间群信息框中的信息也随着,F-43m,空间群的信息而发生变化,),选择,Lattice Parameters,标签，把,a,的数值从,10.00,改为,5.662,单击,Build,按钮一个没有原子的晶格就在,3D model document,中显示出来现在我们就可以添加原子了从菜单栏中选择,Build|Add Atoms,通过它，我们可以把原子添加到指定的位置，其对话框如下：,在,Add Atoms,对话框中选择,Options,标签，确定,Coordinate system,为,Fractional,。

如上所示选择,Atoms,标签，在,Element,文,本框中键入,Al,，然后按下,Add,按钮铝原子就添加到结构中了在,Element,文本框中键入,As,在,a,b,c,文本框中键入,0.25,按,Add,按钮关闭对话框原子添加完毕，我们再使用对称操作工具来构建晶体结构当中剩余的原子这些原子也显示在邻近的单胞中当然，我们也可以通过重新建造晶体结构来移去这些原子从菜单栏中选择,Build|Crystals|Rebuild Crystal.,，,按下,Rebuild,按钮在显示出的晶体结构中那些原子就被移走了我们可以把显示方式变为,Ball and Stick,在模型文档中右键单击，选择,Display Styles,，,按下,Ball and stick,按钮关闭对话框在,3D,视窗中的晶体结构是传统的单胞，它显示的是格子的立方对称如果存在的话，,CASTEP,使用的则是格子的全部对称,.,既包含有两个原子的原胞和包含有,8,个原子的单胞是相对应的,.,不论单胞如何定义，电荷密度，键长，每一类原子的总体能量都是一样的，并且由于使用了较少的原子,使计算时间得以减少从菜单栏中选择,Build|Symmetry|Primitive Cell,。

在模型文档中显示如下：,AlAs,的原胞,2,设置和运行,CASTEP,中的计算,从工具栏中选择,CASTEP,工具，再选择,Calculation,或从菜单栏中选择,Modules|CASTEP|Calculat。