这里会显示出您选择的修订版和当前版本之间的差别。
两侧同时换到之前的修订记录前一修订版后一修订版 | 前一修订版 | ||
atk:分子器件模拟 [2017/03/05 18:51] – [分子器件模拟] dong.dong | atk:分子器件模拟 [2021/12/16 11:23] (当前版本) – ↷ 链接因页面移动而自动修正 59.51.24.40 | ||
---|---|---|---|
行 3: | 行 3: | ||
本教程专门介绍分子器件的计算和分析,研究的经典器件由一个对苯二硫酚(DTB)与两个金电极表面构成,这个结构在很多文献里都有报道。 | 本教程专门介绍分子器件的计算和分析,研究的经典器件由一个对苯二硫酚(DTB)与两个金电极表面构成,这个结构在很多文献里都有报道。 | ||
- | 本教程中的计算时间很长,作为入门的ATK量子输运计算教程,请参考: | + | 本教程中的计算时间很长,作为入门的QuantumATK量子输运计算教程,请参考: |
* [[atk: | * [[atk: | ||
行 9: | 行 9: | ||
* [[atk: | * [[atk: | ||
- | 本教程中使用的ATK-DFT方法参见ATK的使用手册。 | + | 本教程中使用的ATK-DFT方法参见QuantumATK的使用手册。 |
行 19: | 行 19: | ||
<WRAP center important 100%> | <WRAP center important 100%> | ||
=== 注意 === | === 注意 === | ||
- | 这个几何结构未经优化。这对本教程中要讨论的概念没有影响,但对于实际研究中的计算,结构优化至关重要。要了解如何弛豫器件几何结构,请参见: | + | 这个几何结构未经优化。这对本教程中要讨论的概念没有影响,但对于实际研究中的计算,结构优化至关重要。要了解如何弛豫器件几何结构,请参见:[[atk: |
</ | </ | ||
行 51: | 行 51: | ||
<WRAP center info 100%> | <WRAP center info 100%> | ||
- | 并行计算可以大幅度提升计算速度,要了解如何并行计算,请参考:[[atk: | + | 并行计算可以大幅度提升计算速度,要了解如何并行计算,请参考:[[atk: |
</ | </ | ||
行 136: | 行 136: | ||
要同时将结构显示在图上,只需要将 DeviceConfiguration 数据对象拖动到上面已经显示了本征态的Viewer窗口里。 | 要同时将结构显示在图上,只需要将 DeviceConfiguration 数据对象拖动到上面已经显示了本征态的Viewer窗口里。 | ||
- | {{: | + | <WRAP center box 100%> |
- | {{: | + | {{: |
+ | {{: | ||
+ | </ | ||
注意到,所有的态都是关于苯分子平面反对称的,这说明这几个态都是苯分子的 $\pi$-轨道,由与平面垂直的碳原子 6 个 pz 轨道线性组合而成的,因此这样的 $\pi$-轨道共有 6 个。 | 注意到,所有的态都是关于苯分子平面反对称的,这说明这几个态都是苯分子的 $\pi$-轨道,由与平面垂直的碳原子 6 个 pz 轨道线性组合而成的,因此这样的 $\pi$-轨道共有 6 个。 | ||
行 166: | 行 169: | ||
<WRAP center info 100%> | <WRAP center info 100%> | ||
- | 前面两个能量处只有一个透射本征态主导,最后一个能量处则有两个主导的透射本征值。后面会讨论为什么他们没有加在一起给出总的透射谱 T(E)。 | + | 前面两个能量处只有一个透射本征态主导,最后一个能量处则有两个主导的透射本征值。后面会讨论为什么他们加在一起没有得到总的透射谱 T(E)。 |
</ | </ | ||
行 214: | 行 217: | ||
{{ : | {{ : | ||
- | * 选中第一个本征值,点击 Eigenstates 按钮。稍等片刻,第一个本征值对应的本征态就计算好了,Viewer 窗口自动弹出。适当调整 isovalue 数值,得到漂亮的图。 | + | * 选中第一个本征值,点击 Eigenstates 按钮。稍等片刻,对应的本征态就计算好了,Viewer 窗口自动弹出。适当调整 isovalue 数值,得到漂亮的图。 |
{{ : | {{ : | ||
行 331: | 行 334: | ||
为使用 1.6V 和 0V 的 ElectrostaticDifferencePotential 和 ElectronDifferenceDensity 结果计算电压降,我们使用右侧的GridOperations工具。 | 为使用 1.6V 和 0V 的 ElectrostaticDifferencePotential 和 ElectronDifferenceDensity 结果计算电压降,我们使用右侧的GridOperations工具。 | ||
- | * 按住 '' | + | * 按住 '' |
* 点击打开 **GridOperations** 工具; | * 点击打开 **GridOperations** 工具; | ||
* 设置函数为 A1 - B1; | * 设置函数为 A1 - B1; | ||
行 340: | 行 343: | ||
重复上述步骤,计算电压诱导的电荷密度。 | 重复上述步骤,计算电压诱导的电荷密度。 | ||
- | {{: | + | {{ : |
现在计算结果分别保存于 '' | 现在计算结果分别保存于 '' | ||
行 351: | 行 354: | ||
在 Property 里详细调节图示的细节。 | 在 Property 里详细调节图示的细节。 | ||
- | {{: | + | {{ : |
- | {{: | + | {{ : |
如果设置如上图,则得到下面的图: | 如果设置如上图,则得到下面的图: | ||
- | {{: | + | {{ : |
将电压降进行一维作图通常比较有用:选中 EDP 中的 $f(x)$ 图标,点击右侧的 **1D Projector** 工具。投影方法改为 // | 将电压降进行一维作图通常比较有用:选中 EDP 中的 $f(x)$ 图标,点击右侧的 **1D Projector** 工具。投影方法改为 // | ||
- | {{: | + | {{ : |
此处的一维投影作图,C轴取值范围为 0 ~ 1,对应了真实的 z 方向长度。 | 此处的一维投影作图,C轴取值范围为 0 ~ 1,对应了真实的 z 方向长度。 |