24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 1974  |  回复: 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的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

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的回帖
普通回帖

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的回帖

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的回帖

Gina88

木虫 (正式写手)

引用回帖:
3楼: Originally posted by wzbhit at 2012-11-04 11:54:45
收敛标准是1meV/atom。sigma=0.2时,得到TS=-0.089262 除以70=-0.001275已经在可接受的范围了。cutoff和kpoint测试一样的收敛收敛标准,没必要超胞能量变化在1meV,精度太高,因为你的体系70atom,属于比较大的体系 ...

您说的意思是不一定需要TS小于1meV/atom,只需要量级在1meV/atom?
还有,如果按照测试,sigma要比0.2大好多,如果后面的计算sigma也设置得比0.2大好多,得到的TS却比1meV/atom大好多啊,遇到这种情况该怎么办呢?

还有,您说cutoff和kpoint测试一样的收敛标准,指的是delta-total energy达到1meV/atom即可认为收敛,不必一定达到手册上说的delta-total energy达到1meV?比如:对于我的体系,delta-total energy只要小于70meV即可认为收敛,而不必小于1meV?
5楼2012-11-04 23:43:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

colinchen2

金虫 (小有名气)

引用回帖:
4楼: Originally posted by Gina88 at 2012-11-04 23:33:41
非常感谢您的回答!
选一个绝对值最小的是什么意思呢,能够具体一点啊?...

...我楼下那位大侠讲的很详细了嗯~
6楼2012-11-05 07:26:43
已阅   回复此楼   关注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的回帖

Gina88

木虫 (正式写手)

引用回帖:
7楼: Originally posted by wzbhit at 2012-11-05 09:20:31
问题1. 是的,量级相同即可。
问题2。sigma值在计算中要一样,sigma一般经验取0.2 0.1 0.05等值,太多了,就会可能就会出现你说的后面大好多的情况。建议尽可能小,注意太小不好收敛。
问题3. 你的理解正确...

非常感谢!

那这么说来,前面测试就没什么用啦?
还是测试标准更高,以保证后面的计算满足TS在1meV/atom量级?
大侠帮我看看我的sigma测试结果,怎么不是sigma越大,TS(EENTRO)越大啊?
0.10    -0.021321
0.12    -0.036190
0.14    -0.048559
0.16    -0.052802
0.18    -0.047653
0.20    -0.036058
0.22    -0.022676
0.24    -0.011552
0.26    -0.005187
8楼2012-11-05 09:48:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wzbhit

至尊木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★
franch: 金币+2, 谢谢回帖交流, 2012-11-05 14:55:49
Gina88: 金币+3, ★★★★★最佳答案 2012-11-05 21:35:05
引用回帖:
8楼: Originally posted by Gina88 at 2012-11-05 09:48:31
非常感谢!

那这么说来,前面测试就没什么用啦?
还是测试标准更高,以保证后面的计算满足TS在1meV/atom量级?
大侠帮我看看我的sigma测试结果,怎么不是sigma越大,TS(EENTRO)越大啊?
0.10    -0.021321 ...

变化趋势比较诡异,看不出来原因,把INCAR贴出来看看
从明天开始坚持吃早饭呀!
9楼2012-11-05 11:44:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Gina88

木虫 (正式写手)

引用回帖:
9楼: Originally posted by wzbhit at 2012-11-05 11:44:51
变化趋势比较诡异,看不出来原因,把INCAR贴出来看看...

for i in 0.10 0.12 0.14 0.16 0.18 0.20 0.22 0.24 0.26 0.28 0.30
do
cat >INCAR< SYSTEM=Si111_relax
ISTART=0;ICHARG=2

PREC=Accurate
ENCUT=300
GGA=91;VOSKOWN=1
ISMEAR=1;SIGMA=$i
EDIFF=1E-6

LREAL=A
ALGO=Veryfast;IALGO=48

NGX=170
NGY=170
NGZ=80

NBANDS=200

!ISIF=4
!ISYM=2

!NSW=200
!IBRION=2
!POTIM=0.5

LWAVE=.FALSE.
LCHARG=.FALSE.
!
10楼2012-11-05 21:34:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Gina88 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 复试调剂 +20 积极向上; 2026-04-10 22/1100 2026-04-15 09:07 by lanmenghu
[考研] 085404 22408 309分求调剂 +7 lzmk 2026-04-14 7/350 2026-04-14 17:02 by lihaoda1994
[考研] 药学305求调剂 +10 玛卡巴卡boom 2026-04-10 10/500 2026-04-14 15:55 by zs92450
[考研] 人工智能320调剂08工类还有机会吗 +18 振—TZ 2026-04-10 19/950 2026-04-14 10:34 by screening
[考研] 2026硕士调剂_能动_河南农业大学 +4 河南农业大学-能 2026-04-12 4/200 2026-04-13 22:01 by bljnqdcc
[基金申请] 有爆料,一个青年教师卖房得400万,然后换了一个四青帽子 +11 babu2015 2026-04-08 11/550 2026-04-13 16:33 by probebill
[考研] 一志愿西交机械专硕求调剂 +9 求上岸的小王 2026-04-10 9/450 2026-04-13 16:08 by jiangguiquan11
[考研] +10 李多米lee. 2026-04-12 11/550 2026-04-12 22:58 by yuyin1233
[考研] 一志愿郑州大学 22408 305分求调剂 +5 安小满zzz 2026-04-08 5/250 2026-04-12 00:41 by 蓝云思雨
[考研] 0860004 求调剂 309分 +9 Yin DY 2026-04-08 9/450 2026-04-11 22:55 by dongdian1
[考研] 求调剂,一志愿大连理工大学354分 +5 雨声余生 2026-04-11 6/300 2026-04-11 16:12 by 雨声余生
[考研] 085600材料与化工329分求调剂 +16 叶zilin 2026-04-10 16/800 2026-04-11 11:04 by may_新宇
[考研] 085410-273求调剂 +6 X1999 2026-04-10 6/300 2026-04-11 10:32 by Delta2012
[考研] 调剂 +12 卷卷卷心菜_ 2026-04-09 13/650 2026-04-10 22:36 by Ftglcn90
[考研] 中科院总分315求调剂 +8 lallalh 2026-04-09 8/400 2026-04-10 19:30 by dick_runner
[论文投稿] mdpi小修rvr时间四五天了 20+3 哈哈high 2026-04-08 5/250 2026-04-10 16:02 by 北京莱茵润色
[考研] 江苏大学 工科调剂 捡漏 +3 Evan_Liu 2026-04-09 5/250 2026-04-10 10:22 by Evan_Liu
[考研] 0703化学求调剂 +21 不知名的小卅 2026-04-08 21/1050 2026-04-09 18:55 by l_paradox
[考研] 二次调剂求老师收留 +3 笑笑袁 2026-04-08 3/150 2026-04-08 23:50 by 醉在风里
[考研] 296求调剂 +3 汪!?! 2026-04-08 3/150 2026-04-08 22:00 by zhouyuwinner
信息提示
请填处理意见