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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » MS » CASTEP表面能收敛计算求助
171/2返回列表12下一页
查看: 3215  |  回复: 16
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

zhangxx93

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

» 收录本帖的淘帖专辑推荐

CSATEP-表面能

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2017-04-29 11:10:42
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

卡开发发

专家顾问 (著名写手)

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

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur
引用回帖:
3楼: Originally posted by zhangxx93 at 2017-04-30 17:53:51
work function逸出功，第三个回答我没有弄懂，能具体点么？或者有没有这方面的参考文献，谢谢
...

work function或者系统总能实际上都会受到真空大小的影响，而且没有表面偶极的情况（有表面偶极的情形需要修正，不修正的话真空无穷高也不会收敛），只要真空足够大这些物理量就能够收敛，选取一个指标作为参考就行，收敛的情况可能不见的是完全单调的趋势。你可以选择足够的原子层高（不必太大也不必太小），然后在足够真空层中优化完再把那些原子贴到一系列晶格常数ab相同c不同的格子以此作为初始结构，这样每次优化结构可以提升效率，或者这样的情况不优化结构也成。
赞一下(1人) 回复此楼

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

zhangxx93
不一定挂在论坛，计算问题问题欢迎留言。
5楼2017-05-02 05:26:05
已阅   关注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的回帖

zhangxx93

禁虫 (初入文坛)
送红花一朵
本帖内容被屏蔽
6楼2017-05-02 08:36:15
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

zhangxx93

禁虫 (初入文坛)
送红花一朵
本帖内容被屏蔽
7楼2017-05-02 09:29:51
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

M.E-小小子

铁虫 (小有名气)
引用回帖:
4楼: Originally posted by hyde1022 at 2017-05-02 04:47:47
1.优化前用束缚原子么？我分别进行了束缚和不束缚原子，得到的能量相差不大。
取决于你的计算方法和你怎么定义表面能。你给的信息不够，不足以判断。表面能的计算中为了抵消surface dipole，有时候会两个表面同时r ...

前辈能知道我一下有关ms的运算吗，最近用dmol3优化一类似双分子共聚物，始终没法收敛~新手

发自小木虫Android客户端
赞一下 回复此楼
进击的皮卡丘~~~
8楼2017-05-02 09:50:53
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

hyde1022

铜虫 (小有名气)

【答案】应助回帖

引用回帖:
6楼: Originally posted by zhangxx93 at 2017-05-02 08:36:15
我是用图片一中公式2-6进行计算，但实际计算是参考图片二和三进行计算，也就是用图一中的2-4，用tools计算得到的表面积，进行计算时没有除以2。
1.这种计算方法对么？
2.我进行4-10层表面能收敛计算，应该束缚原子 ...

那个计算方法没有问题。对于slab厚度没有一个绝对的值，任何计算都是以收敛为准。你先试试从13层开始算吧，然后再计算15层，不用束缚原子。由于对称性，中间层的原子怎么也是不动的。如果两次计算出来的表面能结果相差零点零几eV/A-2, 就差不多了。另外可以测量一下收敛后结构中，中间几层原子的间距，看是不是接近bulk数值。
你的slab只有七个原子，太少了。厚度不够，表面积也不够。表面积小的时候要取很多K点，不知道你设置够了没有。我建议你把表面积也设置大一点。计算的时候打开symmetry，用gamma-centered K mesh。多做几次测试，这种计算不太花时间
赞一下 回复此楼
9楼2017-05-02 13:06:30
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

zhangxx93

禁虫 (初入文坛)
本帖内容被屏蔽
10楼2017-05-02 15:32:39
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhangxx93 的主题更新
171/2返回列表12下一页
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求材料调剂 +10 隔壁陈先生 2026-03-12 10/500 2026-03-19 09:15 by 无际的草原
[考研] 能源材料化学课题组招收硕士研究生8-10名 +4 脱颖而出 2026-03-16 11/550 2026-03-19 08:11 by 脱颖而出
[考研] 0854可跨调剂，一作一项核心论文五项专利，省、国级证书40+数一英一287 +8 小李0854 2026-03-16 8/400 2026-03-18 14:35 by 搏击518
[考研] 298-一志愿中国农业大学-求调剂 +7 手机用户 2026-03-17 7/350 2026-03-18 14:34 by vgtyfty
[考研] 331求调剂（0703有机化学 +7 ZY-05 2026-03-13 8/400 2026-03-18 14:13 by 007_lilei
[考研] 材料，纺织，生物（0856、0710），化学招生啦 +3 Eember. 2026-03-17 9/450 2026-03-18 10:28 by Eember.
[考研] 301求调剂 +9 yy要上岸呀 2026-03-17 9/450 2026-03-18 08:58 by 无际的草原
[考研] 296求调剂 +5 大口吃饭 身体健 2026-03-13 5/250 2026-03-17 21:05 by 不惑可乐
[考研] 268求调剂 +8 一定有学上- 2026-03-14 9/450 2026-03-17 17:47 by laoshidan
[考研] 211本，11408一志愿中科院277分，曾在中科院自动化所实习 +6 Losir 2026-03-12 7/350 2026-03-17 12:09 by danranxie
[考研] 274求调剂 +5 时间点 2026-03-13 5/250 2026-03-17 07:34 by 热情沙漠
[考研] 机械专硕325，寻找调剂院校 +3 y9999 2026-03-15 5/250 2026-03-16 19:58 by y9999
[考研] 085600调剂 +5 漾漾123sun 2026-03-12 6/300 2026-03-16 15:58 by 漾漾123sun
[考研] 285求调剂 +6 ytter 2026-03-12 6/300 2026-03-16 15:05 by njzyff
[考博] 东华理工大学化材专业26届硕士博士申请 +6 zlingli 2026-03-13 6/300 2026-03-15 20:00 by ryzcf
[考研] 22408总分284求调剂 +3 InAspic 2026-03-13 3/150 2026-03-15 11:10 by zhq0425
[考研] 本科南京大学一志愿川大药学327 +3 麦田耕者 2026-03-14 3/150 2026-03-14 20:04 by 外星文明
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
[考研] 材料080500调剂求收留 +3 一颗meteor 2026-03-13 3/150 2026-03-14 10:54 by peike
[考研] 070303一志愿西北大学学硕310找调剂 +3 d如愿上岸 2026-03-13 3/150 2026-03-13 10:43 by houyaoxu
信息提示
关闭
请填处理意见
关闭