24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1958  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 290调剂生物0860 +27 哇哈哈,。 2026-04-11 30/1500 2026-04-12 21:40 by 学zh
[考研] 一志愿085802 323分求调剂 +10 drizzle_9 2026-04-12 11/550 2026-04-12 20:24 by bljnqdcc
[考研] 086003调剂求助 +21 苏弋万 2026-04-09 22/1100 2026-04-11 20:25 by dongdian1
[考研] 调剂求助 +6 果然有我 2026-04-11 7/350 2026-04-11 16:22 by 明月此时有
[考研] 298求调剂 +9 钉叮咚冬瓜 2026-04-07 11/550 2026-04-11 09:35 by zhq0425
[考研] 材料与化工调剂 +12 否极泰来2026 2026-04-10 13/650 2026-04-11 00:28 by wangjihu
[考研] 一志愿华南理工大学331分材料求调剂 +9 天下ww 2026-04-09 9/450 2026-04-10 22:58 by Ftglcn90
[考研] 265求调剂 +12 风说她早忘了 2026-04-10 13/650 2026-04-10 18:56 by chemisry
[考研] 材料专业344求调剂 +16 hualkop 2026-04-10 21/1050 2026-04-10 17:28 by laoshidan
[考研] 301求调剂 +5 149. 2026-04-10 5/250 2026-04-10 15:45 by 柴小白
[考研] 293调剂 +25 yj1221 2026-04-08 26/1300 2026-04-10 15:02 by 柴小白
[考研] 材料复试求调剂 +20 xhhdjdjsjks 2026-04-09 20/1000 2026-04-10 10:25 by 孙小小12457
[考研] 一志愿鲁东大学071000生物学学硕初试分数276求调剂 +20 慕绝cc 2026-04-09 24/1200 2026-04-10 10:22 by xujun0624
[考研] 284求调剂 +7 让我上岸吧阿西 2026-04-09 7/350 2026-04-09 18:59 by haironglove
[考研] 生物学调剂,一志愿西南大学348,Top期刊一区二作、二区三作,三等奖学金三次 +4 candyyyi 2026-04-09 4/200 2026-04-09 18:39 by l_paradox
[考研] 材料专硕初试分332一志愿西北工业大学, +12 故人?? 2026-04-09 12/600 2026-04-09 18:34 by Ccclqqq
[考研] 328求调剂 +17 lftmya 2026-04-07 18/900 2026-04-09 08:05 by 5268321
[考研] 307求调剂 +14 超级伊昂大王 2026-04-06 14/700 2026-04-08 07:03 by 无际的草原
[考研] 求考研材料调剂 +3 材化李可 2026-04-07 3/150 2026-04-08 00:21 by JourneyLucky
[考研] 071000生物学,一志愿深圳大学296分,求调剂 +12 TIckLw 2026-04-06 13/650 2026-04-07 20:34 by lijunpoly
信息提示
请填处理意见