紧接上文：

Materials Studio官方教程：CASTEP——模拟Pd(110)表面上CO分子体系的STM图【1】

3、创建STM等值面

当计算结束的时候，可以显示电荷差分密度。在此之前关闭所有窗口。

从菜单栏中选择Window | Close All。

现在并排打开 (1x1)和(2x1) 计算的输出结构。

打开(1x1) CO on Pd (1 1 0) CASTEP Energy文件夹中的(1x1) CO on Pd (1 1 0).xsd文件。并打开(2x1) CO on Pd (1 1 0) CASTEP Energy文件夹中的(2x1) CO on Pd (1 1 0).xsd文件。从菜单栏中选择Window | Tile Vertically。

使得(1x1) CO on Pd (110).xsd为当前文档，打开Display Style对话框，在Lattice 选项卡中，设置A方向范围的Max.值为2.00，B方向为2.00，C方向为0.90。选择(2x1) CO on Pd(110).xsd文档，更改B方向范围Max.值为2.00，C方向为0.90。

为了更好地区分表面和吸附质，可能需要更改所有Pd原子的原子显示样式。

对于每个文档，选择每个Pd层的一个原子。在Display Style对话框的Atom选项卡上，选择CPK单选按钮，并将CPK半径更改为0.9。关闭Display Style对话框。

将两个结构定向，使(110)表面排列的原子朝向相同的方向。

模型应与下图类似：

Pd(110)表面上吸附CO分子的1 × 1和1 × 2超晶胞

记下氧原子的z坐标，无论是哪种情况，它们都应约为5.4和5.3 Å。

4、向STM等值面添加额外信息

现在将体积信息导入到两个结构文档中。

单击工具条上的CASTEP按钮，从下拉列表中选择Analysis，打开CASTEP Analysis对话框。从列表中选择STM profile。

打开(1x1) CO on Pd (1 1 0).castep文件，确保其列于输入结果文件Results file中。设置STM偏压STM bias的值为1.0，取消勾选View isosurface on import复选框。使得(1x1) CO on Pd (110).xsd为当前文档，单击Import按钮。

对(2x1) CO on Pd (110).castep输出文件和(2x1) CO on Pd (110).xsd结构文件重复此操作。关闭CASTEP Analysis对话框。

下一步将创建模拟图像，以模拟两个结构在恒定高度模式下的STM实验。这可以通过模拟STM体积信息显示切片来实现。

在体积可视化Volume Visualization工具栏上，单击切片Slice按钮的下拉箭头，然后选择平行于A&B轴Parallel to A & B axis。对另一个文档重复此操作。

此过程将创建一个将成为STM图像的体积切片。在实验中，扫描隧道显微镜通常在实际原子上方进行，并且依赖于隧道。为了在模拟中模拟这种设置，必须将切片设置在氧原子上方相对较高的位置，但也要远离Pd原子底部的周期性镜像。O原子上方1.0 Å的值是一个很好的折衷方案，只要对同类原子进行比较，并为两个模拟设置相同的参数值。

选择切片。在Properties Explorer中，双击切片位置Slice Position属性并将Z值更改为6.4。对另一个文档重复此操作。

注意：对于研究问题中的高质量STM模拟，必须比最高原子高出至少3Å，并使用相应高度较大的单位晶胞。要模拟合理的分辨率，可能必须显著增加平面波截断能。

现在调整颜色尺度以显示相关信息。

选择切片并单击色条Color Maps按钮。使用右箭头按钮将From和To值分别设置为Properties Explorer中显示的SliceMappedMin和SliceMappedMax值。将Spectrum更改为Black-White，并将Bands设置为128。

对另一个文档重复此操作，然后关闭Color Maps对话框。

注：在复杂模拟中，如果色标相同，则可以获得不同吸附质的相对亮度估计值。但是，仅当所有计算中的所有设置完全相同时，此选项才有效。

可以更改表现形式，以使周期性表面的一部分不被切片覆盖。最终结果应与下图类似：

对Pd(110)表面上吸附CO分子的1 × 1和1 × 2超晶胞体系模拟的STM图像。每个图像左侧计算的体积切片对应于预测的STM对比度，而扩展的超晶胞显示STM图像下方原子的位置。

根据该模拟，在STM模拟中，应该很容易区分CO吸附在Pd(110)上的1×1和2×1结构。STM图像的形状表明，对比度主要是由于CO分子的π轨道导致的。在1×1超晶胞中整齐排列，但在2×1结构中偏移。体积切片是仅由与偏压相对应的轨道的密度构成的，即费米能级以上约1 eV的轨道。颜色对应于态的总密度，白色区域对应于最高DOS。

本教程到此结束。

【系列教程】

Materials Studio官方教程：CASTEP——Pd(110)表面上的CO分子的差分电荷密度【1】

Materials Studio官方教程：CASTEP——计算铁磁性铁的声子谱【1】

Materials Studio官方教程：CASTEP——预测锗的热力学性质【1】

Materials Studio是久负盛名计算模拟软件，问世20余年来，经过不断地迭代优化，使其功能异常强大，极易上手，初学者只需通过简单的参数设置和点击鼠标就能完成DFT计算。其计算可靠性久经考验，备受Nature、Science等顶级期刊认可。

华算科技和Materials Studio官方代理深圳浦华系统联合推出Materials Studio建模、计算、分析课程。课程专为零基础学员设计，沿着理论讲解、模型搭建、性质计算、结果分析层层递进讲解，带你快速入门DFT计算。课程极度注重实操，全程线上直播，提供无限回放，课程群在线答疑。（详情点击下方图片跳转）

识别下方二维码报名，或者联系手机13005427160。