CASTEP及一PPt详细教程.pptVIP

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 2005年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对话框如下： 。 精选文档 在Setup标签中，把Task设置为Geometry Optimization，把Quality 设置为Fine，并且把Functional设置为GGA and PW91。 选择Electronic标签，按下More...按钮以得到CASTEP Electronic Options对话框。把Derived grid的设置从Standard改为Fine。关闭CASTEP Electronic Options对话框。 选择Job Control标签，选中你想要运行CASTEP工作的Gateway。 按下CASTEP Calculation对话框中的Run按钮。 优化之后，此构造的晶胞参数应为a=b=c=3.60556?。现在我们可以继续计算优化构造的弹性常数。 2 计算BN的弹性常数 选择CASTEP Calculation对话框中的Setup标签。从Task的下拉清单中选择Elastic Constants。按下More...按钮。 CASTEP Elastic Constants对话框如下： 。 精选文档 把Number of steps for each strain从4增加到6。关闭CASTEP Elastic Constants对话框。按下CASTEP Calculation对话框中的Run按钮。 CASTEP的弹性常数计算任务的结果以一批.castep输出文件的形式给出。这些文件中的每一个文件都代表确定的晶胞在假设的应变模式和应变振幅下的几何优化运行结果。对这些文件的习惯命名方式为：seedname_cij__m__n。对于给定的模式来说，m代表当前的应变模式，n代表当前的应变振幅。 CASTEP可以使用这些结果来分析 每一个运行计算出来的压力张量和产生一个有关弹性属性信息的文件。 从工具栏中选择CASTEP 工具，然后选择Analysis或者从菜单栏中选择Modules | CASTEP | Analysis。 附属性清单中选择Elastic constants，从BN的弹性常数计算工作 。 精选文档 中得到的结果文件应自动显示在Results file选框中。按下Calculate按钮。 在结果文件夹中创立了一个新的文档BN Elastic Constants.txt。 此文档中的信息包括输入的应变摘要和计算出的应力；每一种应变模式的线性适配结果包括适配的质量；在计算出的应力和给定对称性的弹性常数之间的对应;弹性常数Cij和弹性柔量Sij的表格以及最后推知的属性如体积模量和于之相反的可压缩性杨氏模量以及三唯方向上的Poisson比率和建造一个物质为各向同性媒介的模型所需要的Lame 常数。 3 弹性常数文件的描述 对于这种点阵类型需要两种应变模式。对于每一种应变模式，这里都有一个从各自的.castep文件中获取的计算出的压力的摘要。 。 精选文档 =============================================== Elastic constants from Materials Studio: CASTEP =============================================== Summary of the