class_3_2009_用第一原理预测AIAs晶格参数.ppt

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 结构优化－改变原子位置和晶格参数，得到