CASTEP模块计算表面上及吸附能.ppt

显示出CO分子PDOS的图表文档，表示不同原子壳层对态密度的贡献。 选中某一原子，重复上面步骤，可显示该原子对DOS的贡献。见下图。 PDOS of CO molecule CO PDOS O PDOS C PDOS 对(2x1) CO on Pd(110).xsd重复上面的操作。 DOS PDOS PDOS of (2x1) CO on Pd(110) 按Shift键，选中C和O。从工具栏中选择CASTEP工具，然后选择Analysis，选中Density of states，选上Partial，只选sum，按下View按钮 ，结果见后面标有CO的图。此图表示CO对DOS的贡献。同样可选中表面的Pd，或基体的Pd，计算其对DOS的贡献。结果分别见后面标有相应字母的图。 Total DOS Pd (surface) CO Pd (matrix) 体Pd的DOS表现出金属特性。 PDOS of (2x1) CO on Pd(110) 孤立CO的电子态集中在－20、－5、－2.5eV处。与孤立的CO分子的电子态相比，CO约束在表面时，能带有展宽，并向低能方向移动。 PDOS of CO molecule 本次课程结束 作业1：将自己计算的能量发给老师。 作业2：结合总DOS，分析表面Pd、体Pd对材料性质的贡献。 a0 d0 1 1 2 3 4 1 2 3 4 5 6 7 8 a0 5 6 1 2 4 3 zPd-Pd可以由晶格参数a0除以√2得到。 现在我们把距离改为分数长度，可以通过晶格参数（Lattice parameters）工具得到。 在3D模型文档中单击右键，选择Lattice parameters。注意c的值。 (000) 为了计算z的分数坐标，我们仅需要用晶格参数c除以zC（结果为0.382）。 从菜单栏中选择Build | Add Atoms，然后选中Options标签。确保Coordinate system为Fractional。选择Atoms标签，选中C原子，把a改为0，b为0.5，c为0.382, 按下Add按钮。 在3D Viewer上按右键，选Label。 如果我们想确认我们已经正确的设置了模型，可以使用Measure/Change工具。 单击工具栏中Measure/Change工具 的选项箭头，然后选择Distance。在C原子上单击左键，松开后C原子上出现十字叉丝。将鼠标移至Pd原子，再单击左键，则如图所示，显示Pd-C的键长为1.929 ?，数字呈红色 。 将鼠标恢复为 ， 键长数字变为绿色。 鼠标单击键长数字，其颜色变为黄色。这时按Delete，可删除该数字。 下一步是添加氧原子。 *在Add Atoms对话框中，把 Element改为O。 在实验中，C-O键的长度为1.15 ?。 ZO=ZC+ZC-O=4.12+1.15=5.27 O的分数坐标为 ZO/C=5.27/10.78=0.489， * 把这个值写到氧的z分数坐标 上，氧的z坐标值为0.489。按下Add按钮。关闭对话框。 同样用 查看C－O的键长，正确。 计算Pd表面结构时，使用的对称性是P1。但是即使添加了CO，体系仍有更高的对称性。我们可以通过Find Symmetry工具来找到其对称性，并强加对称性（Impose Symmetry），加快计算。 在工具栏中选择Find Symmetry工具 ，按下Find Symmetry按钮，随后按下Impose Symmetry按钮。对称性为PMM2。 PMM2 在3D模型文档中单击右键，选择Display Style。选中Lattice标签，把Style改为Default。 结构如下所示： 在优化几何结构之前，我们先用save project保存，再把它保存到(2x1) CO on Pd(110)文件夹中。 从菜单栏中选择File | Save As...，引导到(2x1) CO on Pd(110)文件。把文档保存为(2x1) CO on Pd(110).xsd。 现在可以优化结构 。 从菜单栏中选择File | Save Project ，然后选择Window | Close All。在Project Explorer中，打开(1x1)CO on Pd(110)文件夹中的(1x1)CO on Pd(110) .xsd。从工具栏中选