第六讲第一原理计算方法简介及Materials Studio中Castep使用教学文稿.ppt

现在它的对称性是PMA2。下面我们来优化它的几何结构。 从工具栏中选择CASTEP工具，然后选择Calculation。 对于本次计算，我们需要改变k点的格子参数，这样我们可以比较本次计算和上次计算的能量值。 选中CASTEP Calculation对话框中的Electronic标签，按下More...按钮。选择k-points标签，把Custom grid parameters改为：a = 2, b = 3, c = 1。关闭对话框，按下Run按钮。 计算开始。计算结束后，在下面的内容中我们需要详细的摘录整个体系的能量。我们可以进行下一步，摘录先前计算的能量。 8．分析能量 在这一部分，我们将要计算化学吸收能DEchem，定义如下： 允许CO分子依着彼此倾斜，然后减低分子的自我排斥力，会导致能量的增加。排斥能可从下面的公式得到： 为计算这些属性，我们需要从CASTEP的文本输出文档中摘录每一次模拟的整个能量。 在Project Explorer中，打开CO molecule/CO CASTEP GeomOpt文件夹中的CO.castep。按下CTRL+F键，搜索Final Enthalpy。向下滑动数行，在下面的表格中记录下出现在“TOTAL ENERGY IS”此行之后的数值。重复此操作，找到其它体系的整个能量，完成下面的表格。 如果获取了所有的能量值，使用上面的等式很简单便可计算出DEchem and DErep。他们的数值大约分别为-1.688eV和0.037eV。 9．分析态密度 下面我们要检查态密度（DOS）的改变。这会使我们对CO在Pd(110).上的连接机制有更深入的了解。为了做到这一点，我们需要显示孤立的CO分子和(2x1) CO 在Pd(110)上的态密度。 在Project Explorer中，打开CO molecule/CO CASTEP GeomOpt文件夹中的CO.xsd。 从工具栏中选择CASTEP工具，然后选择Analysis，选中Density of states，选上Partial，不选f 和sum，但其他的选项都保持原先的状态。按下View按钮 。 Simulation Total Energy(ev) CO molecule -5102 Pd(110) -2103 11 CO on Pd(110) -2103 21 CO on Pd(110) -5.970036697*103 显示出CO分子PDOS的图表文档。如右： 对(2x1) CO on Pd(110).xsd重复上面的操作。 PDOS of (2x1) CO on Pd(110) PDOS of CO molecule 很明显，孤立的CO分子的电子态大约在-20eV, -5eV 和 -2.5eV，比CO约束在表面时的能量要降低一些。 计算BN的弹性常数 背景 当前可应用于大周期性体系的密度泛函理论(DFT)方法的发展已经成为解决材料设计和处理难题的重点。这个理论可以使我们阐释实验数据和预测材料潜在的属性，从一种未知晶体的结构属性到结合能和表面分子的反应都可以预测。这些工具可以用来指导和引导我们设计新材料，允许研究人员推定潜在的化学和物理过程。 绪论 在本指南中，我们将学习如何使用CASTEP来计算弹性常数和其他的力学性能。首先我们要优化BN立方晶体的结构，然后计算它的弹性常数。 本指南主要包括以下内容： 1 优化BN立方晶体的结构 2 计算BN的弹性常数 3 弹性常数文件的描述 1 优化BN立方晶体的结构 在计算弹性常数之前并不是必须要进行几何优化，所以对于实验观测到的结构我们可以产生Cij数据。尽管如此，如果我们完成几何优化包括晶胞优化，可以获得更多相容的结果，然后再计算与理论基态对应的结构的弹性常数。 弹性常数的精确度，尤其是切变常数的精确度，主要取决于SCF计算的品质，特别是布里渊区取样和波函数收敛程度的品质。所以我们应该使用SCF公差为Fine的设置和 k点取样以及来自FFT格子的Fine设置。 首先应导入BN结构。 在菜单栏中选择File | Import，从structures/semiconductors中选中BN.xsd。BN的晶体结构如右： 现在设置几何优化。 从工具栏中选择CASTEP 工具，然后选择Calculation或从菜单栏中选择Modules | CASTEP | Calculation。 CASTEP Calculation对话框如下：