两侧同时换到之前的修订记录前一修订版后一修订版 | 前一修订版后一修订版两侧同时换到之后的修订记录 |
atk:创建分子-表面体系_苯分子在au_111_表面 [2018/05/08 17:12] – [表] xie.congwei | atk:创建分子-表面体系_苯分子在au_111_表面 [2018/05/11 17:06] – [参考文献] xie.congwei |
---|
====== 创建分子-表面体系:苯分子在Au(111)表面 ====== | ====== 创建分子-表面体系:苯分子在 Au(111) 表面 ====== |
| |
版本:2017.0 | 版本:2017.0 |
| |
在本教程中,您将学习如何使用**Virtual NanoLab**(**VNL**)构建苯分子吸附在金(111)表面上两种不同的吸附构型,此处将其表示为Bz@Au(111)。特别地,您将: | 在本教程中,您将学习如何使用 **Virtual NanoLab**(**VNL**) 构建苯分子吸附在金 (111) 表面上两种不同的吸附构型,此处将其表示为 |
| Bz@Au(111)。特别地,您将: |
| |
* 使用Lui等人在参考文献中<color #ed1c24>[LRZ+13]</color>报道过的晶格常数,创建一个Au块体。 | * 使用 Lui 等人在参考文献中 <color #ed1c24>[LRZ+13]</color> 报道过的晶格常数,创建一个 Au 块体。 |
| |
* 参照<color #ed1c24>[LRZ+13]</color>中的参数构建Au(111)面的平板模型。 | * 参照 <color #ed1c24>[LRZ+13]</color> 中的参数构建 Au(111) 面的平板模型。 |
| |
* 在Au(111)面的平板模型上附着上苯分子。 | * 在 Au(111) 面的平板模型上附着上苯分子。 |
| |
* 使用VNL **Builder**中的高级构建工具修改苯分子的原子坐标以构建Bz@Au(111)构形。 | * 使用 VNL **Builder** 中的高级构建工具修改苯分子的原子坐标以构建 Bz@Au(111) 构形。 |
| |
<WRAP center round tip 100%> | <WRAP center tip 100%> |
=== 提示 === | === 提示 === |
构建完这些几何结构后,您可以在ATK中使用其中一种计算方法对它们进行优化,以执行进一步的分析。 | 构建完这些几何结构后,您可以在 ATK 中使用其中一种计算方法对它们进行优化,以执行进一步的分析。 |
</WRAP> | </WRAP> |
| |
<WRAP center round important 100%> | <WRAP center important 100%> |
=== 注意 === | === 注意 === |
本教程中描述的很多VNL功能都相当通用的,除了表面吸附分子的结构外,还可以构建很多复杂构形,在教程[[https://docs.quantumwise.com/tutorials/molecular_junction/molecular_junction.html#molecular-junction|Building a molecular junction]]中可以看到更多举例。 | 本教程中描述的很多VNL功能都相当通用的,除了表面吸附分子的结构外,还可以构建很多复杂构形,在教程 [[https://docs.quantumwise.com/tutorials/molecular_junction/molecular_junction.html#molecular-junction|Building a molecular junction]]中可以看到更多举例。 |
</WRAP> | </WRAP> |
| |
===== 工作流程概要 ===== | ===== 工作流程概要 ===== |
| |
以下是构建Bz@Au(111)构形所需主要步骤的概括。 | 以下是构建 Bz@Au(111) 构形所需主要步骤的概括。 |
| |
1.创建新的VNL项目,打开 {:atk:builder.png?direct&25|}} Builder。 | - 创建新的 VNL 项目,打开 {{:atk:builder.png?direct&25|}} **Builder**。 |
| - 从 **Builder Database** 中导入块体金的原始晶胞。 |
2.从**Builder Database**中导入块体金的原始晶胞。 | - 参照 <color #ed1c24>[LRZ+13]</color> 的数据,使用 **Lattice Parameters** 插件修改晶格常数 a 为 4.159 Å,采用 DFT 和 PBE 密度泛函。 |
| - 打开 **Surface(Cleave)** 插件,构建包含 6 个金属层的 3×3 Au(111) 平板,上方留出 20 Å 的空间。 |
3.参照<color #ed1c24>[LRZ+13]</color>的数据,使用**Lattice Parameters**插件修改金的晶格常数a为4.159 Å,采用DFT和PBE密度泛函。 | - 用 {{:atk:molecular_builder_2.png?direct&|30}} **Molecular Builder** 工具在平板构形上添加苯分子。 |
| - 选定苯原子,用 **Rotate** 插件绕y轴旋转分子 90° 使分子可以平行地吸附于平板表面。 |
4.打开**Surface(Cleave)**插件,构建包含6个金属层的3 × 3 Au(111)平板,上方留出20 Å的空间。 | - 点击 {{:atk:alignmentpoint02a.png?direct&25|}} 图标在分子的几何中心添加一个新原子。 |
| - 使用 {{:atk:move1.png?direct&25|}} **Move** 工具将分子移入 hcp-30° 吸附位。上一步中添加的额外原子用于将分子“吸入”既定位置,随后会将其从构形中删除。 |
5.用 {{:atk:molecular_builder_2.png?direct&|30}} **Molecular Builder**工具在平板构形上添加苯分子。 | - 给构形**重命名**,并在 Stash 中创建它的**副本**。 |
| - 给副本**重命名**,并使用 **Rotate** 插件绕z轴旋转苯 30°,创建 hcp-0° 吸附构形。 |
6.选定苯原子,用**Rotate**插件绕y轴旋转分子90°使分子可以平行地吸附于平板表面。 | |
| |
7.点击 {{:atk:alignmentpoint02a.png?direct&25|}} 图标在分子的几何中心添加一个新原子。 | |
| |
8.使用 {{:atk:move1.png?direct&25|}} **Move**工具将分子移入hcp-30°吸附位。上一步中添加的额外原子用于将分子“吸入”既定位置,随后会将其从构形中删除。 | |
| |
9.给构形**重命名**,并在Stash中创建它的**副本**。 | |
| |
10.给副本**重命名**,并使用**Rotate**插件绕z轴旋转苯30°,创建hcp-0°吸附构形。 | |
===== 详细说明 ===== | ===== 详细说明 ===== |
| |
==== 构建Au(111)表面 ==== | ==== 构建 Au(111) 表面 ==== |
| |
您将在此处由金的原始晶胞开始创建Au(111)的平板模型。 | 您将在此处由金的原始晶胞开始创建 Au(111) 的平板模型。 |
| |
=== 金块体 === | === 金块体 === |
| |
打开**VNL**,创建新项目。设置标题(此处为"//Bz_Au111//"),选择存放的文件夹路径(脚本、数据文件等),点击**OK**完成创建项目。 | 打开 **VNL**,创建新项目。设置标题(此处为 "//Bz_Au111//" ),选择存放的文件夹路径(脚本、数据文件等),点击 **OK** 完成创建项目。 |
| |
点击**OPEN**打开项目,开启您的VNL练习。 | 点击 **OPEN** 打开项目,开启您的 VNL 练习。 |
| |
{{ :atk:1.png?direct&550 |}} | {{ :atk:1.png?direct&550 |}} |
| |
| 在 VNL 的主窗口,打开 {{:atk:builder.png?direct&25|}} **Builder**,点击 **Stash** 旁边的 Add {{:atk:arrow.png?direct&5|}} From Database,打开常用材料的实验结构数据库。 |
| |
在VNL的主窗口,打开 {{:atk:builder.png?direct&25|}} **Builder**,点击**Stash**旁边的Add {{:atk:arrow.png?direct&5|}} From Database,打开常用材料的实验结构数据库。 | 在搜索栏搜索 “gold”,选中结果后点击 {{:atk:add_to.png?direct&30|}} 按钮,将结构导入 **Stash**。 |
| |
在搜索栏搜索“gold”,选中结果后点击 {{:atk:add_to.png?direct&30|}} 按钮,将结构导入**Stash**。 | |
| |
{{ :atk:2.png?direct&550 |}} | {{ :atk:2.png?direct&550 |}} |
| |
金的原始晶胞已经添加在**Stash**区域,现在可以用**Builder plugins**修改参数。 | 金的原始晶胞已经添加在 **Stash** 区域,现在可以用 **Builder plugins** 修改参数。 |
| |
本例中,依次选择Bulk Tools {{:atk:arrow.png?direct&5|}} Lattice Parameters插件,参照文献<color #ed1c24>[LRZ+13]</color>中采用DFT和PBE泛函情况下,将晶格常数a设为4.159 Å。 | 本例中,依次选择 Bulk Tools {{:atk:arrow.png?direct&5|}} Lattice Parameters 插件,参照文献 <color #ed1c24>[LRZ+13]</color> 中采用 DFT 和 PBE 泛函情况下,将晶格常数 a 设为 4.159 Å。 |
| |
{{ :atk:3.png?direct&900 |}} | {{ :atk:3.png?direct&900 |}} |
{{ :atk:4.png?direct&450 |}} | {{ :atk:4.png?direct&450 |}} |
| |
完成以上操作后,关闭插件小工具,返回到Builder。 | 完成以上操作后,关闭插件小工具,返回到 Builder。 |
| |
| |
| |
| |
=== Au(111)平板 === | === Au(111) 平板 === |
| |
下一步是金块体沿(111)方向切开,得到Au(111)表面的平板模型。 | 下一步是金块体沿 (111) 方向切开,得到 Au(111) 表面的平板模型。 |
| |
在 {{:atk:builder.png?direct&25|}} **Builder**里,Builders Surface(Cleave)插件切开金晶体,密勒指数$(h,k,l)=(1,1,1)$: | 在 {{:atk:builder.png?direct&25|}} **Builder** 里, Builders Surface(Cleave) 插件切开金晶体,密勒指数 $(h,k,l)=(1,1,1)$: |
| |
{{ :atk:5.png?direct&800 |}} | {{ :atk:5.png?direct&800 |}} |
{{ :atk:6.png?direct&650 |}} | {{ :atk:6.png?direct&650 |}} |
| |
默认的超胞横向大小为1×1。有些情况下,因为该默认值的存在会非常高效。但在当前的例子中,我们需要3×3的平板结构使相邻苯分子的平面间作用力最小化。 | 默认的超胞横向大小为 1×1。有些情况下,因为该默认值的存在会非常高效。但在当前的例子中,我们需要 3×3 的平板结构使相邻苯分子的平面间作用力最小化。 |
| |
因此,增加晶格矢量$\mathbf{v}_1$和$\mathbf{v}_2$的长度,如此 | 因此,增加晶格矢量 $\mathbf{v}_1$ 和 $\mathbf{v}_2$ 的长度,如此 |
| |
$$v_1 = 3u_1 + 0u_2$$ | $$v_1 = 3u_1 + 0u_2$$ |
{{ :atk:7.png?direct&700 |}} | {{ :atk:7.png?direct&700 |}} |
| |
<WRAP center round tip 100%> | <WRAP center tip 100%> |
=== 提示 === | === 提示 === |
您还可以通过使用鼠标手动地在右边窗口内显示的平面晶胞上移动红色和蓝色矢量实现**交互式地**修改晶格矢量。 | 您还可以通过使用鼠标手动地在右边窗口内显示的平面晶胞上移动红色和蓝色矢量实现**交互式地**修改晶格矢量。 |
</WRAP> | </WRAP> |
| |
最后,点击out-of-plane cell vector $\mathbf{v}_3$ is下方的下拉菜单选择**Non-periodic and slab-like**使其形成具有上下方真空的表面。增加**TOP vacuum**值到20 Å,平板的**Thickness**为6层,如下图所示。 | 最后,点击 out-of-plane cell vector $\mathbf{v}_3$ is下方的下拉菜单选择 **Non-periodic and slab-like** 使其形成具有上下方真空的表面。增加 **TOP vacuum** 值到 20 Å,平板的 **Thickness** 为 6 层,如下图所示。 |
| |
{{ :atk:8.png?direct&600 |}} | {{ :atk:8.png?direct&600 |}} |
| |
完成以上操作后,点击**Finish**。添加Au(111)平板构形到**Stash**区,关闭插件小工具。 | 完成以上操作后,点击 **Finish**。添加 Au(111) 平板构形到 **Stash** 区,关闭插件小工具。 |
| |
{{ :atk:9.png?direct&700 |}} | {{ :atk:9.png?direct&700 |}} |
| |
| |
==== 苯吸附在Au(111)面 ==== | ==== 苯吸附在 Au(111) 面 ==== |
| |
当然,最后一步就是将苯分子吸附在Au(111)平板上。我们的第一个目标就是文献<color #ed1c24>[LRZ+13]</color>中报道的hcp-30°构形。 | 当然,最后一步就是将苯分子吸附在 Au(111) 平板上。我们的第一个目标就是文献 <color #ed1c24>[LRZ+13]</color> 中报道的 hcp-30° 构形。 |
| |
选中**Stash**区的“//Gold(111)//”,点击鼠标右键或按下键盘上F2键将其重命名为“//Au(111)//”。复制该结构,再次用鼠标右键重命名为“//hcp30//”。 | 选中 **Stash** 区的 “//Gold(111)//”,点击鼠标右键或按下键盘上 F2 键将其重命名为 “//Au(111)//”。复制该结构,再次用鼠标右键重命名为 “//hcp30//”。 |
| |
在 {{:atk:builder.png?direct&25|}} **Builder**里,点击 {{:atk:molecular_builder_2.png?direct&|30}} 按钮打开**Molecular Builder**。 | 在 {{:atk:builder.png?direct&25|}} **Builder** 里,点击 {{:atk:molecular_builder_2.png?direct&|30}} 按钮打开 **Molecular Builder**。 |
| |
{{ :atk:12.png?direct&700 |}} | {{ :atk:12.png?direct&700 |}} |
| |
在**Molecular Builder**里,从下拉菜单中点击Fragments {{:atk:arrow.png?direct&5|}} Simple rings {{:atk:arrow.png?direct&5|}} Benzene。 | 在 **Molecular Builder** 里,从下拉菜单中点击 Fragments {{:atk:arrow.png?direct&5|}} Simple rings {{:atk:arrow.png?direct&5|}} Benzene。 |
| |
{{ :atk:13.png?direct&550 |}} | {{ :atk:13.png?direct&550 |}} |
| |
您现在可以将苯分子插入Au(111)构形中:点击Au(111)真空区域的某处将分子放置在那里。然后关闭Molecular Builder窗口。 | 您现在可以将苯分子插入 Au(111) 构形中:点击 Au(111) 真空区域的某处将分子放置在那里。然后关闭 Molecular Builder 窗口。 |
| |
<WRAP center round alert 100%> | <WRAP center alert 100%> |
=== 警告 === | === 警告 === |
需要注意的是,只能点击一次添加苯分子到Au(111)。因为每点击一次,就有一个分子被添加到同一个位置,新的原子会与那些已经存在的分子重叠! | 需要注意的是,只能点击一次添加苯分子到 Au(111)。因为每点击一次,就有一个分子被添加到同一个位置,新的原子会与那些已经存在的分子重叠! |
| |
谨记,您可以随时检查是否有重叠原子:点击Selection Tools {{:atk:arrow.png?direct&5|}} Close Neighbors插件,选择一个相对较小的阈值(例如,0.001 Å)。 | 谨记,您可以随时检查是否有重叠原子:点击 Selection Tools {{:atk:arrow.png?direct&5|}} Close Neighbors 插件,选择一个相对较小的阈值(例如,0.001 Å)。 |
</WRAP> | </WRAP> |
| |
然后,您需要将分子旋转90°。首先用鼠标在分子周围画一个矩形选中所有苯原子。打开Coordinate Tools {{:atk:arrow.png?direct&5|}} Rotate插件,将旋转轴设为“y”,旋转角度为90°,并确保勾选了Rotate around selection center的选项框。点击**Apply**执行旋转操作。 | 然后,您需要将分子旋转 90°。首先用鼠标在分子周围画一个矩形选中所有苯原子。打开 Coordinate Tools {{:atk:arrow.png?direct&5|}} Rotate 插件,将旋转轴设为 “y”,旋转角度为 90°,并确保勾选了 Rotate around selection center 的选项框。点击 **Apply** 执行旋转操作。 |
| |
{{ :atk:rotation.png?direct&700 |}} | {{ :atk:rotation.png?direct&700 |}} |
{{ :atk:rotated.png?direct&700 |}} | {{ :atk:rotated.png?direct&700 |}} |
| |
下一步,移动苯分子到吸附位置hcp-30°。首先,确保已选中所有的苯原子,然后点击 {{:atk:alignmentpoint02a.png?direct&25|}} 图标在分子的几何中心添加一个原子。您将利用这个原子为锚点将分子“吸入”表面的上方。 | 下一步,移动苯分子到吸附位置 hcp-30°。首先,确保已选中所有的苯原子,然后点击 {{:atk:alignmentpoint02a.png?direct&25|}} 图标在分子的几何中心添加一个原子。您将利用这个原子为锚点将分子“吸入”表面的上方。 |
| |
{{ :atk:snapatom.png?direct&450 |}} | {{ :atk:snapatom.png?direct&450 |}} |
再次全部选中包括额外原子在内的所有苯原子,旋转视图为构形的俯视图。 | 再次全部选中包括额外原子在内的所有苯原子,旋转视图为构形的俯视图。 |
| |
<WRAP center round tip 100%> | <WRAP center tip 100%> |
=== 提示 === | === 提示 === |
{{:atk:camera.png?direct&25|}} Camera工具会非常有用:选XY视图平面以获得沿z轴的完美视图。 | {{:atk:camera.png?direct&25|}} Camera 工具会非常有用:选 XY 视图平面以获得沿 z 轴的完美视图。 |
</WRAP> | </WRAP> |
| |
点击 {{:atk:move1.png?direct&25|}} 图标打开**Move**工具,然后点击苯分子中间的额外“吸入原子”,选中其为锚点原子,目标原子周围的阴影应该会变为红色。 | 点击 {{:atk:move1.png?direct&25|}} 图标打开 **Move** 工具,然后点击苯分子中间的额外“吸入原子”,选中其为锚点原子,目标原子周围的阴影应该会变为红色。 |
| |
{{ :atk:move2.png?direct&700 |}} | {{ :atk:move2.png?direct&700 |}} |
| |
在**Move**面板,确保您勾选了**Snap**的选项框。用鼠标移动锚点原子以拖拽分子将其放在第二个平板层的金原子上。 | 在 **Move** 面板,确保您勾选了 **Snap** 的选项框。用鼠标移动锚点原子以拖拽分子将其放在第二个平板层的金原子上。 |
| |
{{ :atk:snapped.png?direct&700 |}} | {{ :atk:snapped.png?direct&700 |}} |
| |
最后,参照文献<color #ed1c24>[LRZ+13]</color>调整分子-表面间的距离为3.62 Å:在**Move**面板,如下图所示在锚点原子位置$r_0$的z坐标处输入18.0272 Å。 | 最后,参照文献 <color #ed1c24>[LRZ+13]</color> 调整分子-表面间的距离为 3.62 Å:在 **Move** 面板,如下图所示在锚点原子位置 |
| $r_0$ 的 z 坐标处输入 18.0272 Å。 |
| |
{{ :atk:zcoordinate.png?direct&700 |}} | {{ :atk:zcoordinate.png?direct&700 |}} |
| |
<WRAP center round important 100%> | <WRAP center important 100%> |
=== 注意 === | === 注意 === |
一般情况下,为了设置所需的分子-表面距离,您一定要先知道金最顶层的笛卡尔坐标Z。这可以通过把指针悬停在Au(111)最顶层的一个原子上方而轻易获得。 | 一般情况下,为了设置所需的分子-表面距离,您一定要先知道金最顶层的笛卡尔坐标 Z。这可以通过把指针悬停在 Au(111) 最顶层的一个原子上方而轻易获得。 |
</WRAP> | </WRAP> |
| |
最后,选中锚点原子并删除。至此,您的“//hcp30//”构形应如下图: | 最后,选中锚点原子并删除。至此,您的 “//hcp30//” 构形应如下图: |
| |
{{ :atk:hcp30.png?direct&400 |}} | {{ :atk:hcp30.png?direct&400 |}} |
| |
<WRAP center round tip 100%> | <WRAP center tip 100%> |
=== 提示 === | === 提示 === |
您可以把构形保存为ATKPython脚本,或者将其**导出**为一些主流的文件格式。 | 您可以把构形保存为 ATKPython 脚本,或者将其**导出**为一些主流的文件格式。 |
</WRAP> | </WRAP> |
| |
| |
| |
=== 创建hcp-0°构形 === | === 创建 hcp-0° 构形 === |
| |
hcp-0°吸附构形与hcp-30°吸附构形相关,只要简单的将后者的吸附的苯分子旋转30°(假设Bz-Au(111)的距离恒定)就可以得到前者。 | hcp-0° 吸附构形与 hcp-30° 吸附构形相关,只要简单的将后者的吸附的苯分子旋转 30° (假设Bz-Au(111)的距离恒定)就可以得到前者。 |
| |
首先,在Stash区域**复制**“//hcp30//”,**重命名**为“//hcp0//”。然后选中新创建副本中的所有苯原子,再次打开Coordinate Tools {{:atk:arrow.png?direct&5|}} Rotate。 | 首先,在 Stash 区域**复制 **“//hcp30//”,**重命名**为 “//hcp0//”。然后选中新创建副本中的所有苯原子,再次打开 Coordinate Tools {{:atk:arrow.png?direct&5|}} Rotate。 |
| |
将构形沿表面所在平面的法线z轴旋转30°,“//hcp0//”构形应如下图: | 将构形沿表面所在平面的法线z轴旋转 30°, “//hcp0//” 构形应如下图: |
| |
{{ :atk:hcp00.png?direct&400 |}} | {{ :atk:hcp00.png?direct&400 |}} |
| |
| |
===== 参考文献 ===== | ===== 参考 ===== |
| |
| <color #ed1c24>[LRZ+13] </color> | (1, 2, 3, 4, 5, 6) Wei Liu, Victor G Ruiz, Guo-Xu Zhang, Biswajit Santra, Xinguo Ren, Matthias Scheffler, and Alexandre Tkatchenko. Structure and energetics of benzene adsorbed on transition-metal surfaces: density-functional theory with van der Waals interactions including collective substrate response. New Journal of Physics, 15(5):53046, 2013. [[http://stacks.iop.org/1367-2630/15/i=5/a=053046|URL: http://stacks.iop.org/1367-2630/15/i=5/a=053046]]. | | * <color #ed1c24>[LRZ+13] </color> | (1, 2, 3, 4, 5, 6) Wei Liu, Victor G Ruiz, Guo-Xu Zhang, Biswajit Santra, Xinguo Ren, Matthias Scheffler, and Alexandre Tkatchenko. Structure and energetics of benzene adsorbed on transition-metal surfaces: density-functional theory with van der Waals interactions including collective substrate response. New Journal of Physics, 15(5):53046, 2013. [[http://stacks.iop.org/1367-2630/15/i=5/a=053046|URL: http://stacks.iop.org/1367-2630/15/i=5/a=053046]]. |
| * 英文原文:[[https://docs.quantumwise.com/tutorials/benzene_au111/benzene_au111.html|https://docs.quantumwise.com/tutorials/benzene_au111/benzene_au111.html]] |