24小时热门版块排行榜    

查看: 2311  |  回复: 11
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

Gina88

木虫 (正式写手)

[求助] ISMEAR=1的测试问题

一个体系有70多个原子,
我计算的时候ISMEAR=1,然后测试sigma,当sigma=0.2的时候,得到TS为-0.036318

但是几何优化的时候sigma=0.2时,得到TS=-0.089262,然后采用sigma=0.15,得到TS=-0.063869

那我sigma应该选0.2呢还是应该选0.15呢?
大家有遇到静态测试盒几何优化时,ISMEAR和sigma取相同的值,却得到不同数值TS的情况吗?

还有大家测试ENCUT和Kpoints的时候,最终精度是不是都达到了delta-total energy小于1meV呢?我觉得我自己测试的时候好像很难到达1meV的精度啊,即使到达1meV的精度,那ENCUT和KPOINTs取值也比较大,因而计算量很大。但是手册上及侯老师的确这么说的。
然后查看文献的时候,发现有的人用的delta-total energy的精度是10meV。
所以对着这个参数测试的精度犯糊涂了,
不知道诸位在具体计算的时候都是怎么测试的,还望不吝赐教!
先谢谢各位大侠了!
回复此楼

» 猜你喜欢

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

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

wzbhit

至尊木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
Gina88: 金币+7, ★★★★★最佳答案 2012-11-05 09:40:33
引用回帖:
5楼: Originally posted by Gina88 at 2012-11-04 23:43:04
您说的意思是不一定需要TS小于1meV/atom,只需要量级在1meV/atom?
还有,如果按照测试,sigma要比0.2大好多,如果后面的计算sigma也设置得比0.2大好多,得到的TS却比1meV/atom大好多啊,遇到这种情况该怎么办呢? ...

问题1. 是的,量级相同即可。
问题2。sigma值在计算中要一样,sigma一般经验取0.2 0.1 0.05等值,太多了,就会可能就会出现你说的后面大好多的情况。建议尽可能小,注意太小不好收敛。
问题3. 你的理解正确
从明天开始坚持吃早饭呀!
7楼2012-11-05 09:20:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答

colinchen2

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★
感谢参与,应助指数 +1
liliangfang: 金币+1, 谢谢指教 2012-11-04 19:47:06
Gina88: 金币+3, ★★★很有帮助 2012-11-04 23:33:59
是应该不同的吧.......

你可以看一下侯大师的教程啊,上面有提到
ISMEAR & SIGMA部分
选一个绝对值最小的......
2楼2012-11-04 10:21:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wzbhit

至尊木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
liliangfang: 金币+2, 谢谢提示 2012-11-04 19:47:13
Gina88: 金币+5, ★★★★★最佳答案 2012-11-04 23:34:25
收敛标准是1meV/atom。sigma=0.2时,得到TS=-0.089262 除以70=-0.001275已经在可接受的范围了。cutoff和kpoint测试一样的收敛收敛标准,没必要超胞能量变化在1meV,精度太高,因为你的体系70atom,属于比较大的体系。
从明天开始坚持吃早饭呀!
3楼2012-11-04 11:54:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Gina88

木虫 (正式写手)

引用回帖:
2楼: Originally posted by colinchen2 at 2012-11-04 10:21:38
是应该不同的吧.......

你可以看一下侯大师的教程啊,上面有提到
ISMEAR & SIGMA部分
选一个绝对值最小的......

非常感谢您的回答!
选一个绝对值最小的是什么意思呢,能够具体一点啊?
4楼2012-11-04 23:33:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 13:35 by C3vABopYbZJS
[教师之家] 某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意! +10 zju2000 2026-09-26 13/650 2026-09-30 10:26 by yexuqing
[有机交流] 同一个分子,一条来自文献,一条来自AI——不告诉你答案,你会选哪条? +4 tianxiaoxian 2026-09-28 7/350 2026-09-30 08:34 by 88817753
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:54 by yCO1Ll7aHtsw
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:47 by yCO1Ll7aHtsw
[博后之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:30 by etmYJ6d2rquH
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:03 by etmYJ6d2rquH
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:04 by etmYJ6d2rquH
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 14:58 by etmYJ6d2rquH
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 GDBe8tDZqE8z 2026-09-28 6/300 2026-09-29 14:41 by etmYJ6d2rquH
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-29 14:30 by etmYJ6d2rquH
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-29 14:24 by etmYJ6d2rquH
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 IDs3scOF0tjC 2026-09-28 7/350 2026-09-29 14:07 by etmYJ6d2rquH
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:42 by ZvyPK8n6Nfki
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-28 23:46 by b8fCuTHqckEt
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-28 23:12 by b8fCuTHqckEt
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[找工作] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 5/250 2026-09-28 22:32 by ez6fGg9abYaj
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 6/300 2026-09-28 22:15 by ez6fGg9abYaj
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-28 19:08 by 9lS3ad5oOymn
信息提示
请填处理意见