应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » MS » CASTEP表面能收敛计算求助
52/1返回列表
查看: 3052  |  回复: 16
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

zhangxx93

禁虫 (初入文坛)
本帖内容被屏蔽

» 猜你喜欢
1楼 2017-04-29 11:10:42
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur
引用回帖:
13楼: Originally posted by zhangxx93 at 2017-05-02 16:02:27
我进行了Al（111）4层，真空层分别为10，17，20，30A，得到的能量基本没有变化Eslab=225.112～225.113eV
...

这样一般足够了。
回复此楼
不一定挂在论坛，计算问题问题欢迎留言。
14楼2017-05-02 16:59:28
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 17 个回答

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur
1、要计算表面能只能固定中间层，否则下表面不弛豫是没什么意义的，如果层高足够的情况，中心的性质和bulk的差不多，此时固定与否差别不大。
2、经验上看 ，层高越薄，越不稳定，表面能越大，你的结果合理。
3、10A不一定够，你实际参考work function的收敛或者能量随着真空的变化。
剩下两问我没啥经验。
赞一下(1人) 回复此楼
不一定挂在论坛，计算问题问题欢迎留言。
2楼2017-04-30 10:49:16
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

zhangxx93

禁虫 (初入文坛)
本帖内容被屏蔽
3楼2017-04-30 17:53:51
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

hyde1022

铜虫 (小有名气)

【答案】应助回帖

感谢参与，应助指数 +1
1.优化前用束缚原子么？我分别进行了束缚和不束缚原子，得到的能量相差不大。
取决于你的计算方法和你怎么定义表面能。你给的信息不够，不足以判断。表面能的计算中为了抵消surface dipole，有时候会两个表面同时relax。这就需要很多原子层，你可以选择束缚中间的几层；如果层数足够多的话不束缚也可以。而且由于对称性，如果你有奇数层原子，中间那层怎么都是不动的。第二种方法，你也可以选择束缚底层原子，只relax一个表面，但这个时候要注意你的能量的reference，因为如果你用Al金属的能量做reference的话，这个方法就不准了。

2.随着层数N的增大，表面能是越来越大还是越来越小？我得到的结果是越来越小
同上，我假设你是用Al 的bulk结构做reference来计算表面能，那么我建议你用一个较厚的slab，然后relax两个表面。若果层数足够的话，得到的表面能是颖收敛的，再加更多层也不会有变化。

3有进行过Al（111）表面能收敛计算的么，我设置的真空层为10A，合理么？
用2中的方法，由于对称性，两个面的surface dipole相互抵消，所以不用太高的真空层。10A也许够，但是多加点更保险。

4.我的参考文献收敛值为1.05，我得到的值为0.62，差距大么？
不知道你们的计算参数和模型有什么不同，很难判断。注意如果relax两个表面的话，最后的结果要除以2。如果你们的模型一样，计算方法一样，用的reference state一样，那么这两个结果差距相当大

5.表面能前后数据相差百分之多少算收敛？
不知道什么意思，是两步SCF计算的能量差吗？
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

zhangxx93
4楼2017-05-02 04:47:47
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 17 个回答
信息提示
关闭
请填处理意见
关闭