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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 所有人都存在会有的《收敛测试的困惑》
101/1返回列表
查看: 1649  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

ustbwgs

木虫 (正式写手)

[交流] 所有人都存在会有的《收敛测试的困惑》 已有5人参与

我们在算任何性质之前都要进行收敛测试，因为所有的数值计算，都是用的iteration integration, 所以只有完全收敛的那个值才具有意义。
比如K点，能量，energy-cut and so on. 但是在得到结果之后，经常会有人提出这样的问题，你看我的这图啊，我的K或energy cut点取这么大了，怎么还不收敛，很多时候问题很简单。画图的时候纵轴范围设得太小了，本来很早就收敛了，但你把纵轴这么一扩大，本来很平缓的让你放大了。可是这就提出一个问题，我们的纵轴究竟应该设多大的范围呢？因为这个没有个绝对得答案。。可是这个怎么去衡量呢？ Instead, 有的人把纵轴设得范围很大。比如0到10000，本来很本没有收敛，可是他这么一画图，就是一条很平缓的水平线。那么，究竟如何从一个图去衡量收敛性呢？
再者，比如你想算表面能，层错能，空位能，这些能量之类的计算，经常会遇到它根本就不收敛，也许把胞做得很大之后，这些能量可以收敛，但那样花费的计算时间会很大很大，你说它不收敛吧，它相差的还不是很大，它总是在线性很微弱减小。例如你做空位层测试，bulk层数测试，它就是收敛的很慢。所以平常找到的一些文献，大部分人都没有用很大的胞，但人家照样能发好文章。这些问题究竟如何处理？我们该如何面对这种问题？
回复此楼

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

模拟知识汇总 科研

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
doingwhatishoulddo!
1楼 2013-12-11 22:38:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

KalaShayminS

铁杆木虫 (著名写手)
★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
sunyang1988: 金币+1, 谢谢交流 2013-12-13 10:36:17
图不能表明收敛性。
收敛取决于你要算什么，例如如果要算表面能，需要精度是0.01eV/晶胞，那就得一直增加各种条件到满足这个为止，然后这个数大概就是error bar。如果只需要0.1eV/晶胞那条件就能取得低一些。如果无法在合理时间内收敛，那找个大家惯用的可接受的大小就行。
赞一下 回复此楼
2楼2013-12-11 22:53:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ustbwgs

木虫 (正式写手)
引用回帖:
2楼: Originally posted by KalaShayminS at 2013-12-11 22:53:40
图不能表明收敛性。
收敛取决于你要算什么，例如如果要算表面能，需要精度是0.01eV/晶胞，那就得一直增加各种条件到满足这个为止，然后这个数大概就是error bar。如果只需要0.1eV/晶胞那条件就能取得低一些。如果无 ...

Good, thank you for your consideration，根据你的建议：我们有表面能的收敛精度是0.01eV=10m eV。希望更多的朋友能补充出其它性质方面的收敛判据和精度。
赞一下 回复此楼
doingwhatishoulddo!
3楼2013-12-12 01:07:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangmomi

金虫 (小有名气)
★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
sunyang1988: 金币+1, 谢谢交流 2013-12-13 10:36:28
引用回帖:
3楼: Originally posted by ustbwgs at 2013-12-12 01:07:27
Good, thank you for your consideration，根据你的建议：我们有表面能的收敛精度是0.01eV=10m eV。希望更多的朋友能补充出其它性质方面的收敛判据和精度。...

楼主可以参考之前讨论收敛的标准。一般来讲，我们选择用力作为收敛标准而不是能量
赞一下 回复此楼
4楼2013-12-12 08:29:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

emilyoyang

木虫 (正式写手)
★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
sunyang1988: 金币+1, 谢谢交流 2013-12-13 10:36:36
收敛性测试时最最基本的问题，文章中不应该提及（至少大家默认你的收敛性性测试是对的，或者压根儿大家不会怀疑你连收敛性测试都做错了）
计算采用晶胞的大小跟能不能发好文章没有直接的逻辑吧
发好文章最重要的还是看idea以及还有老板的牛掰程度
赞一下 回复此楼
5楼2013-12-12 09:20:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ustbwgs

木虫 (正式写手)
引用回帖:
4楼: Originally posted by yangmomi at 2013-12-12 08:29:23
楼主可以参考之前讨论收敛的标准。一般来讲，我们选择用力作为收敛标准而不是能量...

首先感谢关注此话题。我这儿说的收敛标准不是指迭代计算时候的收敛标准，而是在测试各个参数时的收敛标准。
赞一下 回复此楼
doingwhatishoulddo!
6楼2013-12-12 10:30:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ustbwgs

木虫 (正式写手)
引用回帖:
5楼: Originally posted by emilyoyang at 2013-12-12 09:20:08
收敛性测试时最最基本的问题，文章中不应该提及（至少大家默认你的收敛性性测试是对的，或者压根儿大家不会怀疑你连收敛性测试都做错了）
计算采用晶胞的大小跟能不能发好文章没有直接的逻辑吧
发好文章最重要的还 ...

你说得很正确。但我们做计算得都知道，计算中用到得每个参数都必须是经过收敛测试得。只有这样，结果才可信。因为毕竟这是计算，和做实验不一样，必须有充分的理由说明你的结果是正确得才能让人相信。再者我认为不管发好文章也好，不好文章也好，所有的结果必须在自己心里有足够的可信度。
赞一下 回复此楼
doingwhatishoulddo!
7楼2013-12-12 10:34:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)
★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
sunyang1988: 金币+1, 专家考核, 谢谢交流 2013-12-13 10:36:48
其实取代简单的能量图，可以用能量对参数（比如cut-off energy）的导数图。也就是 dE / dEc，这样就不会出现上述的问题。
比如说有人的x轴(cut-off E）是一格5eV，另一个人是一格50eV，那么提出的0.1eV/atom的变化有什么意义呢？显然前面一个方法更快的到达了“收敛”。但是如果用dE / dEc，那么就会非常清晰了。
赞一下(2人) 回复此楼
思想重于技巧，内涵重于表象
8楼2013-12-12 11:49:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ustbwgs

木虫 (正式写手)
引用回帖:
8楼: Originally posted by souledge at 2013-12-12 11:49:49
其实取代简单的能量图，可以用能量对参数（比如cut-off energy）的导数图。也就是 dE / dEc，这样就不会出现上述的问题。
比如说有人的x轴(cut-off E）是一格5eV，另一个人是一格50eV，那么提出的0.1eV/atom的变化 ...

Perfect idea, thank you.
赞一下 回复此楼
doingwhatishoulddo!
9楼2013-12-12 19:07:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lq6865387

木虫 (著名写手)

有志成年
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
8楼: Originally posted by souledge at 2013-12-12 11:49:49
其实取代简单的能量图，可以用能量对参数（比如cut-off energy）的导数图。也就是 dE / dEc，这样就不会出现上述的问题。
比如说有人的x轴(cut-off E）是一格5eV，另一个人是一格50eV，那么提出的0.1eV/atom的变化 ...

在收集到几组阶段能和能量值之后，导数dE / dEc 的数值可以通过什么方法得出呢？？  用origin？？如果用的话如何操作呢？？谢谢赐教
赞一下 回复此楼
Better Late Than Never
10楼2013-12-16 19:39:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ustbwgs 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭