版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請盡快對帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

Gina88

木蟲 (正式寫手)

1ST強(qiáng)帖: 1
應(yīng)助: 9 (幼兒園)
金幣: 3328.8
散金: 124
紅花: 8
帖子: 395
在線: 113小時(shí)
蟲號(hào): 1047684
注冊: 2010-06-26
專業(yè): 理論和計(jì)算化學(xué)

[求助] ISMEAR=1的測試問題

一個(gè)體系有70多個(gè)原子，
我計(jì)算的時(shí)候ISMEAR=1，然后測試sigma，當(dāng)sigma=0.2的時(shí)候，得到TS為-0.036318

但是幾何優(yōu)化的時(shí)候sigma=0.2時(shí)，得到TS=-0.089262，然后采用sigma=0.15，得到TS=-0.063869

那我sigma應(yīng)該選0.2呢還是應(yīng)該選0.15呢？
大家有遇到靜態(tài)測試盒幾何優(yōu)化時(shí)，ISMEAR和sigma取相同的值，卻得到不同數(shù)值TS的情況嗎？

還有大家測試ENCUT和Kpoints的時(shí)候，最終精度是不是都達(dá)到了delta-total energy小于1meV呢？我覺得我自己測試的時(shí)候好像很難到達(dá)1meV的精度啊，即使到達(dá)1meV的精度，那ENCUT和KPOINTs取值也比較大，因而計(jì)算量很大。但是手冊上及侯老師的確這么說的。
然后查看文獻(xiàn)的時(shí)候，發(fā)現(xiàn)有的人用的delta-total energy的精度是10meV。
所以對著這個(gè)參數(shù)測試的精度犯糊涂了，
不知道諸位在具體計(jì)算的時(shí)候都是怎么測試的，還望不吝賜教！
先謝謝各位大俠了！

回復(fù)此樓

» 猜你喜歡

化學(xué)工程321分求調(diào)劑已經(jīng)有22人回復(fù)
322求調(diào)劑已經(jīng)有4人回復(fù)
考研調(diào)劑已經(jīng)有3人回復(fù)
326求調(diào)劑已經(jīng)有8人回復(fù)
333求調(diào)劑已經(jīng)有7人回復(fù)
一志愿東華大學(xué)控制學(xué)碩320求調(diào)劑已經(jīng)有3人回復(fù)
考研化學(xué)學(xué)碩調(diào)劑，一志愿985 已經(jīng)有7人回復(fù)
0703化學(xué)調(diào)劑已經(jīng)有15人回復(fù)
0703化學(xué)調(diào)劑，六級(jí)已過，有科研經(jīng)歷已經(jīng)有14人回復(fù)
326求調(diào)劑已經(jīng)有4人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:

用VASP-DFPT+phonopy計(jì)算石墨烯聲子譜求助已經(jīng)有37人回復(fù)
HSE計(jì)算的一些問題，希望大家?guī)臀医鉀Q謝謝已經(jīng)有13人回復(fù)
關(guān)于Si算例結(jié)果分析已經(jīng)有5人回復(fù)
超晶格結(jié)構(gòu)原子弛豫優(yōu)化出問題已經(jīng)有5人回復(fù)
求助一下如何超算中心提交vasp任務(wù) 已經(jīng)有8人回復(fù)
求大神幫忙解釋截?cái)嗄軠y試中的問題！已經(jīng)有7人回復(fù)
sigma測試問題已經(jīng)有5人回復(fù)
NEB過渡態(tài)問題，搜索的結(jié)果為什么和插入結(jié)果差距這么大？已經(jīng)有6人回復(fù)
關(guān)于VASP中對ENCUT收斂性測試的疑問已經(jīng)有5人回復(fù)
vasp階段能測試求助，謝謝��！急求。已經(jīng)有15人回復(fù)
NEB的計(jì)算出現(xiàn)負(fù)值，測試的其他能量都比初始的00和終止的10低已經(jīng)有17人回復(fù)
VASP計(jì)算ZnO的性質(zhì)，為什么ENCUT的測試曲線是這個(gè)樣子的？已經(jīng)有10人回復(fù)
【重新開帖】KPOINTS的測試腳本問題已經(jīng)有2人回復(fù)
ENCUT測試已經(jīng)有12人回復(fù)
求助：帶電缺陷HSE06總能計(jì)算結(jié)果匪夷所思，請賜教！已經(jīng)有11人回復(fù)
[求助成功]VASP 計(jì)算 graphene 的錯(cuò)誤已經(jīng)有15人回復(fù)
緊急求助vasp5.2 計(jì)算氧化鈰表面吸附遇到的若干問題已經(jīng)有4人回復(fù)
利用GGA+U計(jì)算ZnO帶隙的問題已經(jīng)有3人回復(fù)
vasp4.6問題已經(jīng)有2人回復(fù)
為什么我進(jìn)行ENCUT測試時(shí)是這個(gè)樣子的？問題出在哪里？已經(jīng)有8人回復(fù)
【求助】關(guān)于測試SIGMA的問題已經(jīng)有9人回復(fù)

1樓 2012-11-04 10:18:29

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

回帖支持 ( 顯示支持度最高的前 50 名 )

wzbhit

至尊木蟲 (著名寫手)

應(yīng)助: 84 (初中生)
金幣: 11959.7
散金: 26
紅花: 36
帖子: 1162
在線: 256.7小時(shí)
蟲號(hào): 1128043
注冊: 2010-10-21
專業(yè): 理論和計(jì)算化學(xué)

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
liliangfang: 金幣+2, 謝謝提示 2012-11-04 19:47:13
Gina88: 金幣+5, ★★★★★最佳答案 2012-11-04 23:34:25

收斂標(biāo)準(zhǔn)是1meV/atom。sigma=0.2時(shí)，得到TS=-0.089262 除以70=-0.001275已經(jīng)在可接受的范圍了。cutoff和kpoint測試一樣的收斂收斂標(biāo)準(zhǔn)，沒必要超胞能量變化在1meV，精度太高，因?yàn)槟愕捏w系70atom，屬于比較大的體系。

贊一下(2人)

回復(fù)此樓

從明天開始堅(jiān)持吃早飯呀！

3樓2012-11-04 11:54:45

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

普通回帖

colinchen2

金蟲 (小有名氣)

應(yīng)助: 15 (小學(xué)生)
金幣: 1652
帖子: 59
在線: 33.7小時(shí)
蟲號(hào): 1536228
注冊: 2011-12-13
性別: GG
專業(yè): 半導(dǎo)體物理

【答案】應(yīng)助回帖

★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
liliangfang: 金幣+1, 謝謝指教 2012-11-04 19:47:06
Gina88: 金幣+3, ★★★很有幫助 2012-11-04 23:33:59

是應(yīng)該不同的吧.......

你可以看一下侯大師的教程啊，上面有提到
ISMEAR & SIGMA部分
選一個(gè)絕對值最小的......

贊一下

回復(fù)此樓

2樓2012-11-04 10:21:38

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

Gina88

木蟲 (正式寫手)

1ST強(qiáng)帖: 1
應(yīng)助: 9 (幼兒園)
金幣: 3328.8
散金: 124
紅花: 8
帖子: 395
在線: 113小時(shí)
蟲號(hào): 1047684
注冊: 2010-06-26
專業(yè): 理論和計(jì)算化學(xué)

引用回帖:

2樓: Originally posted by colinchen2 at 2012-11-04 10:21:38
是應(yīng)該不同的吧.......

你可以看一下侯大師的教程啊，上面有提到
ISMEAR & SIGMA部分
選一個(gè)絕對值最小的......

非常感謝您的回答！
選一個(gè)絕對值最小的是什么意思呢，能夠具體一點(diǎn)�。�

贊一下

回復(fù)此樓

4樓2012-11-04 23:33:41

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

Gina88

木蟲 (正式寫手)

1ST強(qiáng)帖: 1
應(yīng)助: 9 (幼兒園)
金幣: 3328.8
散金: 124
紅花: 8
帖子: 395
在線: 113小時(shí)
蟲號(hào): 1047684
注冊: 2010-06-26
專業(yè): 理論和計(jì)算化學(xué)

引用回帖:

3樓: Originally posted by wzbhit at 2012-11-04 11:54:45
收斂標(biāo)準(zhǔn)是1meV/atom。sigma=0.2時(shí)，得到TS=-0.089262 除以70=-0.001275已經(jīng)在可接受的范圍了。cutoff和kpoint測試一樣的收斂收斂標(biāo)準(zhǔn)，沒必要超胞能量變化在1meV，精度太高，因?yàn)槟愕捏w系70atom，屬于比較大的體系 ...

您說的意思是不一定需要TS小于1meV/atom，只需要量級(jí)在1meV/atom？
還有，如果按照測試，sigma要比0.2大好多，如果后面的計(jì)算sigma也設(shè)置得比0.2大好多，得到的TS卻比1meV/atom大好多啊，遇到這種情況該怎么辦呢？

還有，您說cutoff和kpoint測試一樣的收斂標(biāo)準(zhǔn)，指的是delta-total energy達(dá)到1meV/atom即可認(rèn)為收斂，不必一定達(dá)到手冊上說的delta-total energy達(dá)到1meV？比如：對于我的體系，delta-total energy只要小于70meV即可認(rèn)為收斂，而不必小于1meV？

贊一下

回復(fù)此樓

5樓2012-11-04 23:43:04

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

colinchen2

金蟲 (小有名氣)

應(yīng)助: 15 (小學(xué)生)
金幣: 1652
帖子: 59
在線: 33.7小時(shí)
蟲號(hào): 1536228
注冊: 2011-12-13
性別: GG
專業(yè): 半導(dǎo)體物理

引用回帖:

4樓: Originally posted by Gina88 at 2012-11-04 23:33:41
非常感謝您的回答！
選一個(gè)絕對值最小的是什么意思呢，能夠具體一點(diǎn)��？...

...我樓下那位大俠講的很詳細(xì)了嗯～

贊一下

回復(fù)此樓

6樓2012-11-05 07:26:43

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wzbhit

至尊木蟲 (著名寫手)

應(yīng)助: 84 (初中生)
金幣: 11959.7
散金: 26
紅花: 36
帖子: 1162
在線: 256.7小時(shí)
蟲號(hào): 1128043
注冊: 2010-10-21
專業(yè): 理論和計(jì)算化學(xué)

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★
Gina88: 金幣+7, ★★★★★最佳答案 2012-11-05 09:40:33

引用回帖:

5樓: Originally posted by Gina88 at 2012-11-04 23:43:04
您說的意思是不一定需要TS小于1meV/atom，只需要量級(jí)在1meV/atom？
還有，如果按照測試，sigma要比0.2大好多，如果后面的計(jì)算sigma也設(shè)置得比0.2大好多，得到的TS卻比1meV/atom大好多啊，遇到這種情況該怎么辦呢？ ...

問題1. 是的，量級(jí)相同即可。
問題2。sigma值在計(jì)算中要一樣，sigma一般經(jīng)驗(yàn)取0.2 0.1 0.05等值，太多了，就會(huì)可能就會(huì)出現(xiàn)你說的后面大好多的情況。建議盡可能小，注意太小不好收斂。
問題3. 你的理解正確

贊一下

回復(fù)此樓

從明天開始堅(jiān)持吃早飯呀！

7樓2012-11-05 09:20:31

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

Gina88

木蟲 (正式寫手)

1ST強(qiáng)帖: 1
應(yīng)助: 9 (幼兒園)
金幣: 3328.8
散金: 124
紅花: 8
帖子: 395
在線: 113小時(shí)
蟲號(hào): 1047684
注冊: 2010-06-26
專業(yè): 理論和計(jì)算化學(xué)

引用回帖:

7樓: Originally posted by wzbhit at 2012-11-05 09:20:31
問題1. 是的，量級(jí)相同即可。
問題2。sigma值在計(jì)算中要一樣，sigma一般經(jīng)驗(yàn)取0.2 0.1 0.05等值，太多了，就會(huì)可能就會(huì)出現(xiàn)你說的后面大好多的情況。建議盡可能小，注意太小不好收斂。
問題3. 你的理解正確...

非常感謝！

那這么說來，前面測試就沒什么用啦？
還是測試標(biāo)準(zhǔn)更高，以保證后面的計(jì)算滿足TS在1meV/atom量級(jí)？
大俠幫我看看我的sigma測試結(jié)果，怎么不是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

贊一下

回復(fù)此樓

8樓2012-11-05 09:48:31

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wzbhit

至尊木蟲 (著名寫手)

應(yīng)助: 84 (初中生)
金幣: 11959.7
散金: 26
紅花: 36
帖子: 1162
在線: 256.7小時(shí)
蟲號(hào): 1128043
注冊: 2010-10-21
專業(yè): 理論和計(jì)算化學(xué)

【答案】應(yīng)助回帖

★ ★ ★ ★ ★
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
非常感謝！

那這么說來，前面測試就沒什么用啦？
還是測試標(biāo)準(zhǔn)更高，以保證后面的計(jì)算滿足TS在1meV/atom量級(jí)？
大俠幫我看看我的sigma測試結(jié)果，怎么不是sigma越大，TS（EENTRO）越大�。�
0.10 -0.021321 ...

變化趨勢比較詭異，看不出來原因，把INCAR貼出來看看

贊一下

回復(fù)此樓

從明天開始堅(jiān)持吃早飯呀！

9樓2012-11-05 11:44:51

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

Gina88

木蟲 (正式寫手)

1ST強(qiáng)帖: 1
應(yīng)助: 9 (幼兒園)
金幣: 3328.8
散金: 124
紅花: 8
帖子: 395
在線: 113小時(shí)
蟲號(hào): 1047684
注冊: 2010-06-26
專業(yè): 理論和計(jì)算化學(xué)

引用回帖:

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.
!

贊一下

回復(fù)此樓

10樓2012-11-05 21:34:45

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 Gina88 的主題更新

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 326求調(diào)劑 +4	mlpqaz03 2026-03-15	4/200	2026-03-21 19:10 by ColorlessPI
[考研] 286求調(diào)劑 +8	Faune 2026-03-21	8/400	2026-03-21 17:16 by 學(xué)員8dgXkO
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境總分308求調(diào)劑 +7	墨墨漠 2026-03-20	7/350	2026-03-21 16:36 by barlinike
[考研] 0805材料320求調(diào)劑 +3	深海物語 2026-03-20	3/150	2026-03-21 15:46 by 無際的草原
[考研] 求調(diào)劑 +6	Mqqqqqq 2026-03-19	6/300	2026-03-21 08:04 by JourneyLucky
[考研] 346求調(diào)劑[0856] +4	WayneLim327 2026-03-16	7/350	2026-03-21 04:02 by JourneyLucky
[考研] 299求調(diào)劑 +6	△小透明* 2026-03-17	6/300	2026-03-21 02:42 by JourneyLucky
[考研] 271材料工程求調(diào)劑 +8	.6lL 2026-03-18	8/400	2026-03-21 00:58 by JourneyLucky
[考研] 311求調(diào)劑 +5	冬十三 2026-03-18	5/250	2026-03-21 00:16 by JourneyLucky
[考研] 294求調(diào)劑材料與化工專碩 +15	陌の森林 2026-03-18	15/750	2026-03-20 23:28 by JourneyLucky
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑，總分315（英一） +5	sbdksD 2026-03-19	5/250	2026-03-20 22:10 by luoyongfeng
[考研] 0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章 +22	rare12345 2026-03-18	22/1100	2026-03-20 20:39 by zhukairuo
[考研] 353求調(diào)劑 +3	拉鉤不許變 2026-03-20	3/150	2026-03-20 19:56 by JourneyLucky
[考研] 319求調(diào)劑 +3	小力氣珂珂 2026-03-20	3/150	2026-03-20 19:47 by JourneyLucky
[考研] 一志愿福大288有機(jī)化學(xué)，求調(diào)劑 +3	小木蟲200408204 2026-03-18	3/150	2026-03-19 13:31 by houyaoxu
[考研] 本科鄭州大學(xué)物理學(xué)院，一志愿華科070200學(xué)碩，346求調(diào)劑 +4	我不是一根蔥 2026-03-18	4/200	2026-03-19 09:11 by 浮云166
[考研] 312求調(diào)劑 +8	陌宸希 2026-03-16	9/450	2026-03-18 12:39 by Linda Hu
[考研] 308求調(diào)劑 +4	是Lupa啊 2026-03-16	4/200	2026-03-17 17:12 by ruiyingmiao
[考研] 一志愿，福州大學(xué)材料專碩339分求調(diào)劑 +3	木子momo青爭 2026-03-15	3/150	2026-03-17 07:52 by laoshidan
[考研] 318求調(diào)劑 +3	Yanyali 2026-03-15	3/150	2026-03-16 16:41 by houyaoxu

24小時(shí)熱門版塊排行榜

[求助] ISMEAR=1的測試問題

» 猜你喜歡

» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助: