materialstudio一些基础设置问题

1、-. z.CASTEP的任务1. CASTEP能量任务CASTEP能量任务允许您计算指定系统的总能量，以及它的物理性质。除了总能量，原子上的力也会在计算完毕时报告。还创立了一个电荷密度文件，允许使用可视化工具直接观察电荷密度的空间分布。还报告了在计算中使用的monkhorst - pack k点的电子能量，以便在CASTEP分析过程中生成态密度图。能量任务对于研究可靠的构造信息体系的电子特性是非常有用的。只要指定了应力特性，它也可以用来计算没有部自由度的高对称系统的状态方程(即压力体积和/或能量-体积依赖)。注意：在具有部自由度的系统中，可以利用几何优化任务得到状态方程。CASTEP的能量的默

2、认单位是电子伏特(eV)。1 eV= 0.036749308 Ha =23.0605 kcal/mole =96.4853 kJ/mole2. CASTEP几何优化任务CASTEP几何优化任务允许优化几何构造，以获得一个稳定的构造或多态性。这是通过执行一个迭代的过程来完成的，在这个过程中，原子的坐标和可能的原胞参数被调整，从而使构造的总能量是最小的。CASTEP几何优化是基于减小计算力和应力的大小，直到它们变得小于定义的收敛误差。此外，还可以指定一个外部应力量，来模拟在力、压缩、剪切等情况下系统的行为。在这些情况下，部应力量是迭代的，直到它等于施加的外部应力。几何优化的过程一般会产生一个与实际

3、构造相似的模型构造。用CASTEP计算的晶格参数的准确性如图1所示(Milman等，2000)。Figure 1. E*perimental vs. CASTEP calculated lattice parameters状态方程的计算应用流体静压法的几何优化可用于确定材料的体积模量，B，压力导数、B = dB / dP。这个过程包括计算状态方程EOS，它描述了细胞体积对外部流体静压的依赖。该方法与实际实验非常相似:在几何优化对话框中使用最小化选项键确定外部压力，通过对CASTEP进展几何优化来确定压力的单元体积。对P- V数据集的后续分析与实验研究是完全一样的(参见Winkler,1999了

4、解更多细节)。选择一个解析表达式来描述EOS，并将其参数与计算的数据相适应。EOS最受欢送的形式是三阶Birch - Murnaghan方程:V0是平衡体积，Cohen et al.(2000). 等人对各种分析型EOS进展了详细的比拟研究。几何优化方法默认情况下，CASTEP使用BFGS几何优化方法。这种方法通常提供了找到最低能量构造的最快方法，这是在CASTEP支持细胞优化的唯一方案。阻尼分子动力学方法是一种可以和BFGS一样有效的方法，用于具有平面势能面的系统，例如在外表的分子晶体或分子。3CASTEP动力学任务CASTEP动力学任务允许你模拟一个构造中的原子在计算力的影响下如何移动。在

5、执行CASTEP动态计算之前，您可以选择一个热力学集合和相关的参数，定义模拟时间和模拟温度。选择热力学集合结合牛顿运动方程，你可以探索一个系统的恒能外表(NVE动力学)。然而，大多数自然现象发生在系统与环境交换热量的条件下。这些条件可以使用NVT集合(决定性的NOSE集合或随机的Langevin集合)进展模拟。外部压力可以通过使用NPH(常数焓)或NPT(恒定温度)集合来计算动力学计算。定义时间步长集成算法中的一个重要参数是时间步长。为了充分利用计算机时间，应该使用较大的时间步长。但是，如果时间步太大，可能导致积分过程中不稳定和不准确。通常情况下，这是一种系统的漂移运动。注意：量子力学分子动力

6、学的计算通常比基于力场的动力学应用需要更小的时间。然而，在CASTEP中，这不是一个问题，因为使用了波函数和密度外推法。约束在动力学CASTEP在分子动力学模拟中支持线性约束。但是，您只能通过MS界面应用两个最根本的约束:1质心可以是固定的2单个原子位置可以是固定的。更复杂的约束可以使用seedname.cell file.4. CASTEP弹性常量的任务CASTEP弹性常量任务允许你运行一组CASTEP计算，该计算提供所有必要的信息来获取具有任何对称的周期性构造的全部的66量弹性常量。CASTEP会自动生成一组缺陷的构造，它需要对称来产生完整的弹性常数集，然后为理想的构造和所有缺陷的构造提交

7、CASTEP任务。根据初始应变值和计算应力，对弹性常量进展实际计算是在分析阶段进展。弹性常量的精度，特别是剪切常量的精度，在很大程度上取决于SCF计算的质量，特别是关于布里渊区采样的质量和波函数的收敛程度。因此，您应该在SCF选项中选用Fine设置，以及k点取样和Fine 的FFT网格。在计算弹性常量之前没有必要进展几何优化，因此可以为实验观测的构造生成Cij数据。然而，如果你执行完整的几何优化，包括原胞优化，然后计算与理论基态相对应的构造的弹性常数，则会得到更一致的结果。CASTEP为每个有缺陷的构造的部坐标进展最小化处理，同时保持晶格参数不变。可以使用BFGS或damped MD几何优化方

8、法。要选择不同的方法，您必须选择几何优化任务，改变相关设置，然后返回到弹性常数任务。注意：作为弹性常数任务的一局部，CASTEP对没有缺陷的构造进展单点能量计算。这允许您分析电子和容量属性。这个计算的能量和应力并不被用来估计弹性常数。5. CASTEP过渡状态搜索任务CASTEP TS搜索任务允许您优化过渡态。当一个分子或晶体构造被建造时，通常需要把它优化成稳定的几何构造。改良过程是一个迭代的过程，在这个过程中，原子的坐标被调整，使构造的能量被带到一个定点，即该点的原子受的力是零。过渡态是一个定点，它是一个方向上的能量最大值(反响坐标的方向)和其他所有方向上能量是最小的。在化学反响过程中，总能

9、量自然变化。从反响物开场，能量增加到最大值然后降低到生成物的能量。反响路径上的最大能量称为活化能。与这种能量相对应的构造称为过渡态。TS Search任务对于预测化学反响的势垒和决定反响途径特别有用。它还可以用于寻找固体扩散或外表扩散的扩散势垒。同步传输方法CASTEP过渡态搜索使用同步传输方法执行。当反响物和生成物的合理构造存在时，它们的效果最好，但过渡状态的位置是未知的。从反响物和产物开场，同步传输方法插入反响路径以找到过渡状态。线性同步传输(LST)方法对最大的能量执行一次插值。二次同步传输(QST)方法通过约束最小化来对能量最大化进展搜索，以将过渡状态细化到高程度。CASTEP可用的选

10、项包括:1LST Ma*imum执行一次LST最大化，在反响物和生成物之间的最大限度。这是最快但最不准确的选项。这种方法决定的TS构造通常需要进一步的优化。2LST/Optimization执行一次LST最大化，接着是向反响路径的共轭的能量最小化。这就产生了一个能量更低的构造，更接近于真实的TS，而不是一个简单的LST。继续最小化步骤，直到到达能量的最小值或使共轭方向的数量耗尽。3Halgren - Lipsb是一种有限的LST/Optimization，它的设计是为了再现由Halgren和Lipsb推广的算法。在确定了LST最大值后，该方法执行共轭梯度最小化，但只在一个方向上进展。4plet

11、e LST/QST通过执行LST/Optimization计算开场。以这种方式获得的TS近似用于执行QST最大化。从这一点开场，另一个共轭梯度最小化开场执行。循环重复，直到一个静止点被定位，或者被允许的QST步骤数被耗尽。对所使用算法的完整描述在 HYPERLINK mk:MSITStore:C:Program%20Files%20(*86)AccelrysMaterials%20Studio%204.3sharedocSMCastep.chm:/Html/refsCastep.htm l govind_2003 Govind et al. (2003).中给出。6. CASTEP属性CAST

12、EP属性任务允许您在完成单点能量、几何优化或动力学运行之后计算电子、构造和振动特性。注意：为了利用CASTEP属性任务来计算属性，必须在当前project中给出一个适宜的仿真结果。您还可以在CASTEP计算对话框中使用Properties选项卡，要求将这些属性计算为CASTEP运行的一局部。您可以使用CASTEP分析对话框查看结果。CASTEP可以生成的属性如下:1带构造:在布里渊区，在模拟过程中使用电子电荷密度和电位，沿高对称方向的电子特征值在价层和导带是非自洽的，。2态密度: 使用电子电荷密度和在模拟过程中产生的电势，在一个fine HYPERLINK mk:MSITStore:C:Pro

13、gram%20Files%20(*86)AccelrysMaterials%20Studio%204.3sharedocSMCastep.chm:/Html/thBrillouin.htm l MPpoints Monkhorst-Pack grid的价带和导带中的电子特征值是非自洽的。3电子密度差:对原子密度的线性组合或构造中包含的原子集合密度的线性组合的电子密度差。4NMR:计算了化学屏蔽量和电场梯度。注意:NMR的CASTEP是另一个许可模块NMR CASTEP的一局部。只有购置了这个模块，才能执行NMR计算。5光学性质:计算电子带间跃迁的矩阵元素。CASTEP分析对话框可以用来生成网格

14、和包含可测量的光学属性的图表文档。6轨道:提供关于电子波函数的信息。这允许您可视化各种电子状态(轨道)的3D分布。此信息还需要对STM配置文件进展可视化。7声子:对于声子色散的运行，声子频率和沿着布里渊区的高对称方向的特征向量进展计算。在声子密度的计算中，声子频率和特征向量是在Monkhorst-Pack grid计算的。在分析过程中需要这些信息，以显示所有状态的总声子密度。它也被用来计算热力学性质。注意：Phonon计算考虑到现有的固定原子约束，而不考虑整个任务设置。这些固定原子被排除在振动性质的计算之外，这符合partial Hessian方法。小窍门：您必须检查Edit Constrai

15、nts对话框的Atom选项上的Fi*分数位置复选框，可以从Modify菜单中访问。8极化率和红外光谱: 随着红外强度(反响电场在红外围)变化的光学(=)和直流(= 0)介电常数或光学(=)和静态(= 0)分子极化率被计算。介电常数与固体材料有关，而极化率和红外强度与用超胞方法构建的分子有关。9Population analysis:Mulliken analysis被执行。计算了键数和角动量解析的原子电荷(以及自旋极化计算的磁矩)。另外，还生成了状态局部密度(PDOS)计算所需的权重。10应力:应力量的计算和写入到seedname.castep文件中。这个信息是有用的，例如，你执行一个几何优化

16、运行，其中的单元参数是固定的，可以检查晶格到平衡位置的距离是多少。例如，一个点缺陷的超细胞研究应该与给定系统的理论基态相对应的固定原胞一块进展。在几何优化后的应力值给出了与超级细胞近似有关的弹性效应的大小的提示。7.在固体中模拟无序许多晶体构造具有静态的位置无序。MS允许你通过定义混合原子来模拟无序。1混合原子无序晶体中的原子位置可以用一个混合原子来描述，它由两个或多个元素类型组成。相对浓度可以设置为任意数量的原子，其中总浓度必须是100%。混合原子描述是固体溶液、金属合金、无序矿物等的最常用表现形式。2技术上的限制目前，CASTEP的分子动力学(MD)算法的实现与虚拟晶体近似(VCA)还不相

17、容，VCA是用来在固体中建立无序模型的。在目前的版本中，没有任何一个MD(NVE,NPT，等等)为无序晶体工作。以下属性对于无序系统是不可用的:数量分析振动性质(声子色散，声子态密度)光学性质此外，基于dampedMD的几何优化是不可用的(因为MD本身是不允许的)。Atom占据态不能用于基于虚拟晶体近似(VCA)的当前CASTEP实现。局部原子占用可以被看作是混合型原子的一个特例，其中一局部的混合物是一个虚无的原子。这种情况不受一些技术和根本原因的支持。上面列出的限制分为三类:对混合原子(数量分析、振动特性)的性质的定义不明确或根本没有意义。有可能进展计算的特性，但这需要更多的理论工作(光学，

18、局部原子占据)。由于纯粹的技术原因而被限制的特性(分子动力学)。这些限制是由用户界面施加的;如果您尝试使用不允许的设置组合运行，则当您选择运行或保存文件时，会显示一个消息框，解释为什么无法执行运行。CASTEP计算的设置1电子选项设置电子选项必须设置为能量、几何优化、动力学、弹性常数和TS搜索计算。用于属性计算的设置从将要重新启动的计算的输入文件中复制，并且不受任何通过界面进展的更改的影响。用于定义计算的根本设置可以在CASTEP计算对话框的设置和电子设置选项中得到。您可能需要修改的重要设置如下:Quality计算质量：决定了计算的准确度和本钱System parameters系统参数：描述系

19、统的电荷和自旋状态E*change-correlation functional交换关联泛函：确定DFT描述的水平Pseudopotentials赝势：影响计算的准确性和本钱Basis set基组设置指定平面波基组的质量k-points ：确定布里渊区抽样的准确性SCF：指定用于电子最小化的方案Quality为了帮助您选择CASTEP参数的适当值，可以使用四个预设的参数。这些预设的参数在质量上有Coarse粗糙，适合快速近似计算；到Ultra-fine超精细，用于计算需要高精度的计算。当您选择质量设置时，设置以下参数与Coarse, Medium, 和 Fine设置相对应的截断能的值是对应于特

20、定的元素，它取自存储在数据库中赝势文件。这些值是由单个原子和双原子分子的收敛试验确定的，它们分别对应于2eV/atom、0.3V/atom和0.1V/atom的总能量收敛值。在3D构造文档中所选的元素中所选择的质量的最大值被承受。Quality设置也影响了诸如几何优化、弹性常数计算和过渡态搜索等任务的收敛误差值。System parameters系统电荷CASTEP允许在具有电荷的系统上执行计算。这对于研究半导体的缺陷(Milman et al .，1993和Hakala等，2000)是很有用的。注意：在这个版本的CASTEP中没有实现Makovl Payne关联，它改良了总能量对带电周期系统

21、的超胞大小的收敛性。这意味着具有原胞单元带电系统获得的结果不则可靠。设置系统电荷从菜单栏中选择模块| CASTEP |计算，以显示CASTEP计算对话框。选择Setup键。在电荷场中指定系统电荷。自旋极化通过CASTEP进展自旋极化的DFT计算，可以研究磁系统。该设置用于过渡金属氧化物，*些无机外表研究，以及含有磁性元素的金属系统(Fe,Co,Mn,Ni)。自旋极化CASTEP计算通常由一个可变磁矩进展。然而，这一过程并不能保证找到最低的能态。DFT的解可以收敛到各种局部极小值，从而找到亚稳态。最常见的解对应于高自旋和低自旋态。CASTEP生成的状态依赖于初始的磁场配置，因为解可能会收敛到最近

22、的局部最小值，而不是全局最小值。有两种方法来定义初始的磁配置:要么指定均匀地分布在每单元的总磁矩，或者提供关于单元中每个原子的自旋的绝对值和方向(向上或向下)的详细信息。前一种方法可以用于相对简单的系统，只有两种解决方案(磁性和非磁性)。后一种方法，它指定了系统中原子的自旋状态，更一般，并且提供了更多的灵活性。有可能建立铁磁或反铁磁计算，以获得不同的起始自旋排列。用均匀分布的初始磁矩建立自旋极化计算从菜单栏中选择模块| CASTEP |计算，以显示CASTEP计算对话框。选择Setup选项。检查自旋极化复选框，并将使用正式自旋作为初始复选框。指定初始的自旋值。为了优化旋转值，选择Electro

23、nic，点击More按钮|CASTEP电子选项对话框。选择SCF(自洽场)选项卡。取消修复占用复选框并检查优化的总自旋复选框。用原子矩的指定顺序建立自旋极化的计算注意：只有当电子构造的密度混合方案被用于电子构造的最小化时，CASTEP才正确处理原子的初始磁矩。CASTEP将忽略原子信息，并将对所有能带或EDFT最小化使用均匀分布的磁矩。总力矩的值是由系统中正式自旋的总和决定的。为原子设置Set formal spins正式的旋转。在设置了正式的自旋后，操作find and impose the symmetry 。从菜单栏中选择模块| CASTEP |计算，以显示CASTEP计算对话框。选择设

24、置选项卡。检查自旋极化并将正式自旋作为初始的复选框。选择电子标签，点击More|CASTEP电子选项对话框。选择SCF(自洽场)选项卡。选择密度混合作为电子最小化装置。取消修复占用复选框并检查优化的总自旋复选框。在优化自旋状态时包含足够数量的空带是很重要的。空带的默认数量是4，这意味着自旋值最大可增加到8。如果CASTEP发现的实际极小值对应于高自旋状态(尽管这不是一个计算从高自旋和收敛到低旋转时的情况)，则从低自旋状态开场的计算将会非常慢。注意:如果空带的数量缺乏以容纳高自旋状态，在SCF中CASTEP可能无法收敛，特别是在一个具有可变占用数(金属系统)的计算中。2.设置e*change-c

25、orrelation 泛函CASTEP支持一个局域交换关联泛函，LDA，三个梯度修正的交换关联泛函，GGA和一组非局域泛函，用于总能量自洽计算。1LDA泛函CA-PZ是CASTEP中唯一可用的局域(LDA)泛函。它基于Ceperley and Alder(1980)由Perdew and Zunger(1981)参数化的数据。2GGA泛函GGA泛函提供了比LDA泛函更好的电子子系统的总体描述。LDA描述倾向于低估原子，从而使键长和细胞体积通常被低估了几个百分点，而体积的模量则相应地被高估了。GGA纠正了这个错误，但可能会被高估，从而高估键长。PBE ( HYPERLINK mk:MSITSto

26、re:C:Program%20Files%20(*86)AccelrysMaterials%20Studio%204.3sharedocSMCastep.chm:/Html/refsCastep.htm l pbe_1996 Perdew et al., 1996)是默认的交换关联泛函。特别推荐用于研究与金属外表相互作用的分子，尽管对大块系统也相当可靠。PW91(Perdew et al .，1992)泛函应该用于与文献数据进展比拟，因为它是使用最广泛的GGA泛函。PW91泛函所获得的结果应该与PBE所获得的结果非常相似。另一个非局域泛函，RPBE ( HYPERLINK mk:MSITSto

27、re:C:Program%20Files%20(*86)AccelrysMaterials%20Studio%204.3sharedocSMCastep.chm:/Html/refsCastep.htm l hammer_1999 Hammer et al., 1999)，是PBE泛函的改良形式，旨在改良金属外表的描述。最新的GGA势是WC, HYPERLINK mk:MSITStore:C:Program%20Files%20(*86)AccelrysMaterials%20Studio%204.3sharedocSMCastep.chm:/Html/refsCastep.htm l wu_

28、cohen_2006 Wu and Cohen (2006).。作者声称对晶格常数、晶体构造和金属外表的能量有显著的改良，超过了最受欢送的Perdew-Burke-Ernzerhof (PBE) GGA.3非局域交换关联泛函非局部交换关联泛函是由广义的kohn sham方程(Seidl et al,1996)所产生的，与LDA或GGA的计算相比，它旨在改良绝缘体和半导体的带隙的描述。这种额外的准确性来自于更多的时间消耗计算。一些可用非局部交换关联泛函近似:HFHF-LDAs*s*-LDAPBE0推荐的泛函是s*-LDA 和PBE0;其余局部主要用于测试其他由Hartree-Fock方案实现的程

29、序包。注意：使用这种方法有一些限制: 非局部交换只能用于具有能量和几何优化任务(没有原胞优化)的绝缘体。它与应力、NMR、声子或极化计算不相容。这种公式不能与原子或超软势混用。选择交换关联泛函从菜单栏中选择模块| CASTEP |。选择Setup选项卡。选择一个局域(LDA)，一个梯度校正(GGA)或非局域泛函从与功能选项相关联的第一个下拉列表。从第二个下拉列表中选择具体的函数。3.设置赝势1Ultrasoftpseudopotentials超软赝势和norm-conserving pseudopotentials模守恒赝势CASTEP在默认情况下使用的是超软赝势(USP)。它通常比模守恒赝势

30、(NCP)更准确、更有效。数据库中的USP文件有一个.usp或.uspcc扩展，NCP文件有一个.recpot扩展。注意：对于LDA和GGA交换关联泛函，应使用不同的超软赝势。数据库中的命名约定是_.uspcc用于LDA泛函和_PBE.uspcc用于GGA泛函。在默认情况下，MS试图找到最适合于交换关联势的赝势，尽管这种行为可以被拒绝。注意：与NCP计算相比，使用USPs进展光学特性的计算不则准确(详见光学特性)。然而，由于误差预计很小，无法量化。因此，如果需要的话，您应该使用USPs来进展几何优化，然后使用NCPs进展额外的能量和光学特性的计算。与NCP相比照，USP使得CASTEP计算，可

31、以使用较小的截断能，所以在计算时间上产生明显的优势。然而，USP的形式更加复杂，对于诸如声子、NMR特性或局域交换关联泛函的线性响应的实现来说，这是非常棘手的问题。因此，CASTEP在处理一些任务和性质的计算上，只能使用norm-conserving potentials。MS可用的NCPs数据库相当旧，在*些情况下可能包含尚未经过充分测试的赝势。因此，可以使用开源软件OpiumThe Optimized Pseudopotential Interface / Unification Module产生NCPs。这个包允许您使用最新的赝势生成技术产生.recpot格式的文件。除了改变全局选择的赝势，还可以根据交换关联泛函的选择值来更改每个元素的选择。提示：个人选择使用其他的soft versions of the USPs版本，对诸如氧和硅等元素也很有用，在这里可以这需要更小的截断能，但可能比标准势的可转移性更小。2Pseudopotentials generated on the flyCASTEP代码允许您generate pseudopotentials on the fly，例如，您可以提供参数，而不是来自数据库的文件来govern the generation。这种方法有很多优点，例如，在原子和固态计算中使用了一样的交换关联泛函;通过改变核心的半径，可以产生较软