24小时热门版块排行榜    

查看: 6422  |  回复: 30

valenhou001

至尊木虫 (职业作家)

如果你的POSCAR是通过ISIF=3优化出来的话,那你的这个POSCAR就有点问题了。
对于slab的超原胞,除非你修改过vasp源代码,用ISIF=3来优化这样体系的结构是有问题的,它会调整真空层的厚度。
11楼2014-11-30 20:56:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
whywhy01: 金币+10, ★★★很有帮助 2014-12-04 09:54:12
建议你做这样的测试:

1)增加电场强度的大小,看真空层区域的静电势的变化
2)将基矢的a23, a32改为0.0000后再计算测试一下。它们的值很小,接近0。可改为0。
12楼2014-11-30 20:59:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

whywhy01

金虫 (小有名气)

引用回帖:
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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

magicmonk

至尊木虫 (著名写手)

引用回帖:
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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

【答案】应助回帖

引用回帖:
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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

更正:"表面方向上的晶格常数根据真空层厚度以及原子层厚度来取"-----"表面法向方向上的晶格常数根据真空层厚度以及原子层厚度来取"
16楼2014-12-01 18:43:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

whywhy01

金虫 (小有名气)

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

我明白您的意思了,我的这个结构是放在多层薄膜中用ISIF=3优化过的,取出中间两层直接做静态计算,不知是否可以?
另外电场的作用应该如何调整参数使其能够体现出来呢?会不会是在纯金属薄膜里用DIPOL的方法加电场效果本来就不明显?
17楼2014-12-01 20:25:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

whywhy01

金虫 (小有名气)

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

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

我用您的方法试过了,晶体结构符合要求,但是电场还是没有明显作用。于是用加减电子的方法模拟电场了。还是很感谢您!
18楼2014-12-04 09:53:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

【答案】应助回帖

引用回帖:
18楼: Originally posted by whywhy01 at 2014-12-04 10:53:55
我用您的方法试过了,晶体结构符合要求,但是电场还是没有明显作用。于是用加减电子的方法模拟电场了。还是很感谢您!...

加减电子数,还不是模拟加静电场。 通过偶极矩、锯齿型势、平行板电容器加静电场会改变上下两个表面上的电荷分布,但是不改变体系的总电荷。

你试过把电场强度增大一些,比如0.2, 0.3或更大一点。
19楼2014-12-04 10:24:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

天空般明朗

金虫 (小有名气)

为什么设置IBRION = -1 呢?
人生就是不停的战斗!
20楼2014-12-04 11:27:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 whywhy01 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +20 木水思豆 2026-08-25 23/1150 2026-08-26 06:35 by echo8914667
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 有没有大神帮我看看基金代码 28+4 1234567wang 2026-08-24 10/500 2026-08-25 19:15 by lfy8008
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 19/950 2026-08-25 17:57 by chick875
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 2026国自然函评费到账 +20 羊腰板 2026-08-21 23/1150 2026-08-25 16:34 by zsna
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 22/1100 2026-08-25 15:36 by 医学老男孩
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +3 苏东坡二世 2026-08-23 3/150 2026-08-25 10:35 by shisan1313
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +5 process2012 2026-08-23 7/350 2026-08-25 09:42 by huixian257
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 只有每年这种时候来逛逛小木虫 +25 yaoyewhu2008 2026-08-20 27/1350 2026-08-25 08:01 by Equinoxhua
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
信息提示
请填处理意见