24小时热门版块排行榜    

查看: 3971  |  回复: 13

奋力拼搏

至尊木虫 (文坛精英)

[求助] 能量的收敛精度是什么意思?如何得到?谢谢!!

用CASTEP计算单点能进行收敛性测试时,能量的收敛精度(单位是eV/atom)是什么意思?它是如何得到的?如果体系含有多个原子呢?请不吝赐教!!十分感谢!!!!
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

天道酬勤
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
奋力拼搏: 金币+5, 有帮助, 谢谢您的指点,如果能够再深入一点就更好了 2013-08-14 08:58:56
fegg7502: 金币+1, 3ks 2013-08-14 10:33:45
fegg7502: 应助指数+1, 专家考核, 鼓励交流 2013-08-14 10:33:57
两次迭代,系统(包含所有原子)总能量的差值,再除以系统的总原子数,该然后将数值的绝对值与设置的精度相比较,若小于预设精度,则迭代停止,计算完成。

» 本帖已获得的红花(最新10朵)

MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
2楼2013-08-13 18:05:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhenwuhuang

至尊木虫 (文学泰斗)

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +1
奋力拼搏: 金币+2, 有帮助, 谢谢 2013-08-14 08:59:45
fegg7502: 金币+1, 鼓励交流 2013-08-14 10:34:09
castep 结构优化不收敛
http://muchong.com/html/200905/1338065.html
3楼2013-08-13 18:46:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

飞行鸟

至尊木虫 (文坛精英)

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +1
奋力拼搏: 金币+2, 有帮助, 谢谢 2013-08-14 08:59:55
fegg7502: 金币+1 2013-08-14 10:34:18
一般是10的-5次方,但是如果要更好的话-6次方。 一般默认设置。
快乐吧
4楼2013-08-13 18:56:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

奋力拼搏

至尊木虫 (文坛精英)

送红花一朵
引用回帖:
2楼: Originally posted by 月只蓝 at 2013-08-13 18:05:50
两次迭代,系统(包含所有原子)总能量的差值,再除以系统的总原子数,该然后将数值的绝对值与设置的精度相比较,若小于预设精度,则迭代停止,计算完成。

首先感谢您的回答!确定截断能时,一般以相邻两个截断能计算出的能量之差小于0.001eV作为标准,请问“能量之差小于0.001eV中的能量差是体系的总能量差还是体系中每个原子的平均能量(体系的总能量除以体系中原子的个数之后的值)差”?渴望得到您的指点,谢谢!!
天道酬勤
5楼2013-08-14 08:53:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

奋力拼搏

至尊木虫 (文坛精英)

引用回帖:
4楼: Originally posted by 飞行鸟 at 2013-08-13 18:56:26
一般是10的-5次方,但是如果要更好的话-6次方。 一般默认设置。

用CASTEP进行收敛性测试时,为了确定截断能,一般以相邻两个截断能计算出的能量之差小于0.001eV作为标准,请问“能量之差小于0.001eV中的能量差是体系的总能量差?还是体系的总能量差除以体系的原子数目后得到的值”?渴望得到您的指点,谢谢!!
天道酬勤
6楼2013-08-14 08:57:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖


fegg7502: 金币+1, 应助指数+1, 鼓励交流 2013-08-14 10:34:30
引用回帖:
5楼: Originally posted by 奋力拼搏 at 2013-08-14 08:53:57
首先感谢您的回答!确定截断能时,一般以相邻两个截断能计算出的能量之差小于0.001eV作为标准,请问“能量之差小于0.001eV中的能量差是体系的总能量差还是体系中每个原子的平均能量(体系的总能量除以体系中原子的 ...

一般认为是体系的总能量之差。
但实际上,对于原子数很多的体系,这一标准很难达到;对于大体系,采用0.01eV的总能量标准,或者平均原子能量之差。
收敛性测试的目的是证明,即便采用更高的精度设置,也不会改变我们要下的科学判断或结论。
选取不同的Ecut、k-points为x轴,以体系总能量为y轴,随着Ecut、k-points的增大,总能量减小,随后趋于平直。
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
7楼2013-08-14 09:59:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

奋力拼搏

至尊木虫 (文坛精英)

引用回帖:
7楼: Originally posted by 月只蓝 at 2013-08-14 09:59:36
一般认为是体系的总能量之差。
但实际上,对于原子数很多的体系,这一标准很难达到;对于大体系,采用0.01eV的总能量标准,或者平均原子能量之差。
收敛性测试的目的是证明,即便采用更高的精度设置,也不会改变 ...

什么是大体系?标准是什么?它一般含有多少个原子?谢谢
天道酬勤
8楼2013-08-14 10:02:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
奋力拼搏: 金币+10, ★★★很有帮助, 非常感谢 2013-08-14 10:19:21
fegg7502: 金币+2, 应助指数+1, 3ks 2013-08-14 10:34:51
引用回帖:
8楼: Originally posted by 奋力拼搏 at 2013-08-14 10:02:24
什么是大体系?标准是什么?它一般含有多少个原子?谢谢...

大体系是个定性的概念,也认为30个原子以上,也有50个原子以上,计算的精度要求越高,被定义的原子数往往越少。
做收敛性测试的时候,可以先按照体系总能量之差的标准来考察,对于原子数不多的体系,这一标准容易达到,如果难以达到,退而求其次。
其实,总能量之差0.001eV的标准,0.01eV的标准以及平均原子能量之差都是人为定性,如7楼说的,选取不同的Ecut、k-points为x轴,以体系总能量为y轴,随着Ecut、k-points的增大,总能量减小,随后趋于平直,选取趋于平直的第一个点,或者更安全的第二个点,就是计算精度与耗时的良好平衡。
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
9楼2013-08-14 10:11:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

奋力拼搏

至尊木虫 (文坛精英)

引用回帖:
9楼: Originally posted by 月只蓝 at 2013-08-14 10:11:54
大体系是个定性的概念,也认为30个原子以上,也有50个原子以上,计算的精度要求越高,被定义的原子数往往越少。
做收敛性测试的时候,可以先按照体系总能量之差的标准来考察,对于原子数不多的体系,这一标准容易 ...

十分感谢!!
天道酬勤
10楼2013-08-14 10:18:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 奋力拼搏 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 不要再数国自然申请书的 filecode 的分隔符个数了 +11 sadpencil 2026-07-12 16/800 2026-07-13 11:25 by Vivilian
[公派出国] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 obzynn3tuv 2026-07-11 4/200 2026-07-13 11:25 by emlgkvoijw
[有机交流] 轴手性化合物的确定 10+3 trimethoxy 2026-07-07 4/200 2026-07-13 11:18 by celan
[考研] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +4 obzynn3tuv 2026-07-11 5/250 2026-07-13 11:05 by emlgkvoijw
[基金申请] 明天E口面上会评 +4 布布和一二 2026-07-12 5/250 2026-07-13 10:44 by zzybjr
[公派出国] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 obzynn3tuv 2026-07-11 4/200 2026-07-13 10:35 by emlgkvoijw
[考博] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 obzynn3tuv 2026-07-11 4/200 2026-07-13 10:25 by emlgkvoijw
[论文投稿] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 obzynn3tuv 2026-07-11 4/200 2026-07-13 10:15 by emlgkvoijw
[硕博家园] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 obzynn3tuv 2026-07-11 3/150 2026-07-13 09:56 by emlgkvoijw
[考博] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 obzynn3tuv 2026-07-11 4/200 2026-07-13 00:24 by v2i90kmrcu
[基金申请] 祈祷青基必中 50+6 hadengry 2026-07-09 16/800 2026-07-12 20:31 by 登山虫虫
[基金申请] 生命口会评 +11 wuke100666 2026-07-07 12/600 2026-07-11 19:30 by hhs666
[基金申请] 关于如何从代码看上不上会 +17 且听虎啸 2026-07-09 23/1150 2026-07-10 19:51 by sdfapple719
[有机交流] chemdraw 20+5 ?慕鱼? 2026-07-08 6/300 2026-07-10 15:37 by 紫枫星域
[基金申请] b口会评 +6 当空照 2026-07-08 6/300 2026-07-09 20:49 by tonyliu0401
[考研] 在职考研 +6 wu5d 2026-07-08 7/350 2026-07-09 16:40 by 化工小霸王呀
[论文投稿] MSER送审了还被拒稿 +3 S778950906 2026-07-08 5/250 2026-07-09 15:46 by S778950906
[有机交流] 硝基还原 20+3 暮年飞雪 2026-07-07 5/250 2026-07-09 13:30 by tfang
[基金申请] 化学口会评?应用基础研究和基础研究,希望能分开 +3 Tide man 2026-07-08 3/150 2026-07-09 11:36 by dede5556
[有机交流] Suzuki偶联中硼酸酯原料掉下的硼酸酯是否可以自身偶联,生成B2Pin2杂质 500+3 天下大同 2026-07-06 3/150 2026-07-07 10:31 by martyr_oyx
信息提示
请填处理意见