【答案】应助回帖

11楼2014-11-30 20:56:07

valenhou001

至尊木虫 (职业作家)

1ST强帖: 13
应助: 241 (大学生)
金币: 25701.7
散金: 602
红花: 166
帖子: 3782
在线: 873.8小时
虫号: 1007127
注册: 2010-04-27
专业: 凝聚态物性 II ：电子结构

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
whywhy01: 金币+10, ★★★很有帮助 2014-12-04 09:54:12

建议你做这样的测试：

１）增加电场强度的大小，看真空层区域的静电势的变化
２）将基矢的a23, a32改为0.0000后再计算测试一下。它们的值很小，接近０。可改为0。

12楼2014-11-30 20:59:24

whywhy01

金虫 (小有名气)

应助: 35 (小学生)
金币: 725.8
红花: 4
帖子: 113
在线: 46小时
虫号: 1515708
注册: 2011-11-29
性别: GG
专业: 凝聚态物性 II ：电子结构

引用回帖:

10楼: Originally posted by valenhou001 at 2014-11-30 20:53:02
原子层放在中间，对结果没有影响。

为啥你的基矢a23, a32不是为零？　你计算的什么类型的表面（表面方向和表面大小）？
难道你的POSCAR是通过ISIF=3优化出来的吗？...

初始的超胞是四方的，结构优化时ISIF=3产生了一点畸变，考虑到这个变化量不大就直接拿过来了。Ni和Co都是FCC结构，底面是(111)面，z轴沿<111>方向

13楼2014-11-30 21:00:45

magicmonk

至尊木虫 (著名写手)

应助: 172 (高中生)
金币: 11451.3
红花: 24
帖子: 1124
在线: 3108.8小时
虫号: 1191315
注册: 2011-01-17
专业: 理论和计算化学

引用回帖:

11楼: Originally posted by valenhou001 at 2014-11-30 20:56:07
如果你的POSCAR是通过ISIF=3优化出来的话，那你的这个POSCAR就有点问题了。
对于slab的超原胞，除非你修改过vasp源代码，用ISIF=3来优化这样体系的结构是有问题的，它会调整真空层的厚度。

请问如果要修改源代码的话，大致思路是如何的？

我们最深的恐惧不是我们能力不够，我们最深的恐惧是能力超越了极限。

14楼2014-12-01 09:25:12

valenhou001

至尊木虫 (职业作家)

1ST强帖: 13
应助: 241 (大学生)
金币: 25701.7
散金: 602
红花: 166
帖子: 3782
在线: 873.8小时
虫号: 1007127
注册: 2010-04-27
专业: 凝聚态物性 II ：电子结构

【答案】应助回帖

引用回帖:

13楼: Originally posted by whywhy01 at 2014-11-30 22:00:45
初始的超胞是四方的，结构优化时ISIF=3产生了一点畸变，考虑到这个变化量不大就直接拿过来了。Ni和Co都是FCC结构，底面是(111)面，z轴沿<111>方向...

如我前面说的，你这样得到的POSCAR有问题。重新计算。

对于表面体系的超原胞的计算，不要采用ISIF=3.
正确或合理的步骤：先优化和计算好块体材料的结构参数（晶格常数，原子坐标），再切表面构造出表面体系的超原胞，对这样得到的表面体系的超原胞进行只优化原子位置（也就是ISIF=2），平行于表面的晶格常数固定，表面方向上的晶格常数根据真空层厚度以及原子层厚度来取。

15楼2014-12-01 18:42:03

valenhou001

至尊木虫 (职业作家)

1ST强帖: 13
应助: 241 (大学生)
金币: 25701.7
散金: 602
红花: 166
帖子: 3782
在线: 873.8小时
虫号: 1007127
注册: 2010-04-27
专业: 凝聚态物性 II ：电子结构

更正："表面方向上的晶格常数根据真空层厚度以及原子层厚度来取"-----"表面法向方向上的晶格常数根据真空层厚度以及原子层厚度来取"

16楼2014-12-01 18:43:32

whywhy01

金虫 (小有名气)

应助: 35 (小学生)
金币: 725.8
红花: 4
帖子: 113
在线: 46小时
虫号: 1515708
注册: 2011-11-29
性别: GG
专业: 凝聚态物性 II ：电子结构

引用回帖:

16楼: Originally posted by valenhou001 at 2014-12-01 18:43:32
更正："表面方向上的晶格常数根据真空层厚度以及原子层厚度来取"-----"表面法向方向上的晶格常数根据真空层厚度以及原子层厚度来取"

我明白您的意思了，我的这个结构是放在多层薄膜中用ISIF=3优化过的，取出中间两层直接做静态计算，不知是否可以？
另外电场的作用应该如何调整参数使其能够体现出来呢？会不会是在纯金属薄膜里用DIPOL的方法加电场效果本来就不明显？

17楼2014-12-01 20:25:02

whywhy01

金虫 (小有名气)

应助: 35 (小学生)
金币: 725.8
红花: 4
帖子: 113
在线: 46小时
虫号: 1515708
注册: 2011-11-29
性别: GG
专业: 凝聚态物性 II ：电子结构

引用回帖:

15楼: Originally posted by valenhou001 at 2014-12-01 18:42:03
如我前面说的，你这样得到的POSCAR有问题。重新计算。

对于表面体系的超原胞的计算，不要采用ISIF=3.
正确或合理的步骤：先优化和计算好块体材料的结构参数（晶格常数，原子坐标），再切表面构造出表面体系的超 ...

我用您的方法试过了，晶体结构符合要求，但是电场还是没有明显作用。于是用加减电子的方法模拟电场了。还是很感谢您！

18楼2014-12-04 09:53:55