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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » MS » CASTEP表面能收敛计算求助
171/2返回列表12下一页
查看: 3463  |  回复: 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下一页
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的？ +16 爱看书的可乐 2026-08-26 18/900 2026-08-28 11:12 by winsaint
[基金申请] 为什么资助数各大高校都创新高，自己申请怎么就这么难 +10 Kittylucky 2026-08-27 11/550 2026-08-28 10:45 by Kittylucky
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜？ 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解！ +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 系统查不到 +10 董八千 2026-08-26 10/500 2026-08-26 16:30 by Equinoxhua
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 2026国自然函评费到账 +22 羊腰板 2026-08-21 25/1250 2026-08-26 15:00 by zuocuiping
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 今天务委会开完了，明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 国合现在查不到了吗？ +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
信息提示
关闭
请填处理意见
关闭