版塊導(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í)熱門版塊排行榜

返回列表

dfjk-123

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 29.9
散金: 757
帖子: 119
在線: 24.5小時(shí)
蟲號(hào): 665485
注冊: 2008-12-01
性別: MM
專業(yè): 理論和計(jì)算化學(xué)

[求助] 相同的vasp版本，不同機(jī)器上計(jì)算同一個(gè)體系，能量相差近10ev

各位：
我在兩個(gè)不同的LINUX集群下用vasp計(jì)算同一個(gè)體系，輸入文件完全相同。兩個(gè)集群安裝的都是vasp5.2.2版本，結(jié)構(gòu)優(yōu)化完成后能量居然相差了近10ev，幾何結(jié)構(gòu)兩個(gè)基本相同。大家說說這是怎么回事？

回復(fù)此樓

» 猜你喜歡

303求調(diào)劑已經(jīng)有3人回復(fù)
有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～已經(jīng)有7人回復(fù)
梁成偉老師課題組歡迎你的加入已經(jīng)有10人回復(fù)
一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項(xiàng)國獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級立項(xiàng) 已經(jīng)有8人回復(fù)
070300化學(xué)319求調(diào)劑已經(jīng)有3人回復(fù)
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑已經(jīng)有4人回復(fù)
268求調(diào)劑已經(jīng)有7人回復(fù)
290求調(diào)劑已經(jīng)有11人回復(fù)
材料與化工求調(diào)劑已經(jīng)有4人回復(fù)
293求調(diào)劑已經(jīng)有11人回復(fù)

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

不同版本vasp 計(jì)算相同的超胞的能量差了幾十個(gè)meV 這算是合理或者是有問題？已經(jīng)有10人回復(fù)
VASP靜態(tài)自洽計(jì)算，能量升高7eV，然后局域磁矩消失，請問是什么問題？已經(jīng)有8人回復(fù)
vasp計(jì)算單個(gè)氧原子和分子的能量已經(jīng)有20人回復(fù)
求助：vasp 計(jì)算固體相變，能量。已經(jīng)有3人回復(fù)
VASP計(jì)算反鐵磁和鐵磁態(tài)能量差高達(dá)1eV以上，是否可信？已經(jīng)有10人回復(fù)
求助：vasp 計(jì)算氫原子能量是-1.16ev 而不是 -13.6ev，why ? 已經(jīng)有18人回復(fù)
【求助】vasp靜態(tài)計(jì)算能量值錯(cuò)誤已經(jīng)有21人回復(fù)
【其他】vasp結(jié)合phonon計(jì)算材料不同溫度下能量的問題已經(jīng)有9人回復(fù)
【轉(zhuǎn)帖】用VASP計(jì)算氫分子的能量已經(jīng)有8人回復(fù)

1樓 2012-11-24 19:06:11

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

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

賀儀

鐵桿木蟲 (著名寫手)

1ST強(qiáng)帖: 5
應(yīng)助: 93 (初中生)
貴賓: 0.636
金幣: 5937
散金: 1063
紅花: 38
沙發(fā): 2
帖子: 2499
在線: 588.7小時(shí)
蟲號(hào): 995125
注冊: 2010-04-12
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

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

★
WDD880227: 金幣+1, 多謝指教 2012-11-28 14:31:47

引用回帖:

7樓: Originally posted by dfjk-123 at 2012-11-27 15:58:54
我后來試了NBANDS數(shù)不變，加大了K點(diǎn)，能量上升了。減少K點(diǎn)的話，warning就沒有了。
另外一個(gè)測試表明：增加了NBANDS數(shù)，K點(diǎn)數(shù)不變，能量下降了。在此基礎(chǔ)上又增加了K點(diǎn)數(shù)，能量在此基礎(chǔ)上又有下降。搞不清楚K點(diǎn)和 ...

‘我后來試了NBANDS數(shù)不變，加大了K點(diǎn)，能量上升了。減少K點(diǎn)的話，warning就沒有了�！�
這明顯是K點(diǎn)不夠，warming沒有了是因?yàn)樵谀承㎏點(diǎn)的最高帶上已經(jīng)沒有占據(jù)了，但這個(gè)沒有是錯(cuò)誤的提示。

對于K點(diǎn)的夠不夠你測試的時(shí)候不知道的參數(shù)就用默認(rèn)的，比如說NBANDS。測試的時(shí)候增加K點(diǎn)的密度，直到能量達(dá)到一個(gè)穩(wěn)定的值（誤差在0.幾個(gè)eV內(nèi)，越小越好）。

總的來講我認(rèn)為你的NBANDS數(shù)量不夠，NELECT/2 + NIONS/2 (non-spinpolarized)這個(gè)沒錯(cuò)，但你的電子數(shù)是不是數(shù)對了呢？

贊一下(1人)

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

dfjk-123

嘟啦啦嘟啦啦嘟

8樓2012-11-27 19:14:36

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

souledge

專家顧問 (著名寫手)

專家經(jīng)驗(yàn): +174
1ST強(qiáng)帖: 19
應(yīng)助: 303 (大學(xué)生)
金幣: 11306.4
散金: 1123
紅花: 108
帖子: 2814
在線: 916.1小時(shí)
蟲號(hào): 941486
注冊: 2010-01-12
專業(yè): 結(jié)構(gòu)陶瓷
管轄: 第一性原理

★
WDD880227: 金幣+1, 謝謝提示 2012-11-28 14:31:59

如果計(jì)算過程沒有警告問題，總能（別管你取哪個(gè)）有差別都無所謂，只要生成能之類都沒變化就可了。

贊一下(1人)

回復(fù)此樓

思想重于技巧，內(nèi)涵重于表象

9樓2012-11-27 19:35:22

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

longwen36

鐵桿木蟲 (正式寫手)

應(yīng)助: 11 (小學(xué)生)
金幣: 6165.9
紅花: 16
帖子: 344
在線: 208.4小時(shí)
蟲號(hào): 766090
注冊: 2009-05-08
性別: GG
專業(yè): 功能陶瓷

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

引用回帖:

2樓: Originally posted by dfjk-123 at 2012-11-24 19:13:05
剛檢查了兩個(gè)地方生成的OUTCAR文件，發(fā)現(xiàn)NBANDS數(shù)是不一樣的，這可能是導(dǎo)致能量差的原因。但是怎么會(huì)差這么大呢？我的INCAR中并沒有給定NBANDS數(shù)量啊。同樣版本的vasp，給出的能帶數(shù)應(yīng)該一樣才對啊。這是為什么呢？

貌似也我也有這樣的問題，
不同的機(jī)器NBANDS數(shù)不一樣，所以算出來的結(jié)果也不一樣。

INCAR里面沒有給NBANDS，它好像會(huì)根據(jù)你算的CPU數(shù)來設(shè)定。貌似與并行的CPU數(shù)有關(guān)。
所以最好在INCAR里面最好固定NBANDS數(shù)，就好像不會(huì)有這樣的問題。

贊一下(1人)

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

dfjk-123

16樓2012-12-04 16:53:22

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

普通回帖

dfjk-123

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 29.9
散金: 757
帖子: 119
在線: 24.5小時(shí)
蟲號(hào): 665485
注冊: 2008-12-01
性別: MM
專業(yè): 理論和計(jì)算化學(xué)

剛檢查了兩個(gè)地方生成的OUTCAR文件，發(fā)現(xiàn)NBANDS數(shù)是不一樣的，這可能是導(dǎo)致能量差的原因。但是怎么會(huì)差這么大呢？我的INCAR中并沒有給定NBANDS數(shù)量啊。同樣版本的vasp，給出的能帶數(shù)應(yīng)該一樣才對啊。這是為什么呢？

贊一下

回復(fù)此樓

2樓2012-11-24 19:13:05

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

賀儀

鐵桿木蟲 (著名寫手)

1ST強(qiáng)帖: 5
應(yīng)助: 93 (初中生)
貴賓: 0.636
金幣: 5937
散金: 1063
紅花: 38
沙發(fā): 2
帖子: 2499
在線: 588.7小時(shí)
蟲號(hào): 995125
注冊: 2010-04-12
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

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

★ ★
感謝參與，應(yīng)助指數(shù) +1
franch: 金幣+2, 謝謝回帖交流， 2012-11-26 17:25:44

你比較的是什么能量？絕對的嗎？vasp屬于周期性的結(jié)構(gòu)，能量的零點(diǎn)沒法定義，在計(jì)算能量中有些數(shù)值是根據(jù)你所用的數(shù)據(jù)庫決定的，但是如果你取相對的能量的話不花有什么變化的

贊一下

回復(fù)此樓

嘟啦啦嘟啦啦嘟

3樓2012-11-25 10:04:22

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

dfjk-123

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 29.9
散金: 757
帖子: 119
在線: 24.5小時(shí)
蟲號(hào): 665485
注冊: 2008-12-01
性別: MM
專業(yè): 理論和計(jì)算化學(xué)

引用回帖:

3樓: Originally posted by 賀儀 at 2012-11-25 10:04:22
你比較的是什么能量？絕對的嗎？vasp屬于周期性的結(jié)構(gòu)，能量的零點(diǎn)沒法定義，在計(jì)算能量中有些數(shù)值是根據(jù)你所用的數(shù)據(jù)庫決定的，但是如果你取相對的能量的話不花有什么變化的

我是取的energy without entropy后面的能量。我又測試了下，將兩個(gè)測試的NBANDS數(shù)設(shè)置為一樣，計(jì)算出的結(jié)果基本上差不多了。我認(rèn)為還是NBANDS的原因。但是，INCAR文件沒有設(shè)置NBANDS數(shù)時(shí)，同樣的VASP版本，在不同機(jī)器上卻給的不一樣，這是為什么？
另外，資源有限，我的作業(yè)需要使用兩個(gè)機(jī)器做計(jì)算，如果得出的能量差別如此大，那算出來的相對能量還能準(zhǔn)確嗎？應(yīng)該相信哪一個(gè)？

贊一下

回復(fù)此樓

4樓2012-11-25 12:08:56

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

賀儀

鐵桿木蟲 (著名寫手)

1ST強(qiáng)帖: 5
應(yīng)助: 93 (初中生)
貴賓: 0.636
金幣: 5937
散金: 1063
紅花: 38
沙發(fā): 2
帖子: 2499
在線: 588.7小時(shí)
蟲號(hào): 995125
注冊: 2010-04-12
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

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

★ ★
franch: 金幣+2, 謝謝回帖交流， 2012-11-26 17:25:56

引用回帖:

4樓: Originally posted by dfjk-123 at 2012-11-25 12:08:56
我是取的energy without entropy后面的能量。我又測試了下，將兩個(gè)測試的NBANDS數(shù)設(shè)置為一樣，計(jì)算出的結(jié)果基本上差不多了。我認(rèn)為還是NBANDS的原因。但是，INCAR文件沒有設(shè)置NBANDS數(shù)時(shí)，同樣的VASP版本，在不同機(jī) ...

你的兩個(gè)電腦的設(shè)置等不一樣計(jì)算出來的都會(huì)有差別，影響很大的就是你編譯時(shí)以及計(jì)算時(shí)所調(diào)用的數(shù)據(jù)庫。對于絕對能量的比較最好是能夠在同一臺(tái)機(jī)械上計(jì)算。相對能量到無所謂。

贊一下

回復(fù)此樓

嘟啦啦嘟啦啦嘟

5樓2012-11-25 20:37:33

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

dfjk-123

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 29.9
散金: 757
帖子: 119
在線: 24.5小時(shí)
蟲號(hào): 665485
注冊: 2008-12-01
性別: MM
專業(yè): 理論和計(jì)算化學(xué)

引用回帖:

5樓: Originally posted by 賀儀 at 2012-11-25 20:37:33
你的兩個(gè)電腦的設(shè)置等不一樣計(jì)算出來的都會(huì)有差別，影響很大的就是你編譯時(shí)以及計(jì)算時(shí)所調(diào)用的數(shù)據(jù)庫。對于絕對能量的比較最好是能夠在同一臺(tái)機(jī)械上計(jì)算。相對能量到無所謂。...

好的，謝謝。還有個(gè)問題想請教下，我在OUTCAR中發(fā)現(xiàn)了下面的警告：
|                                                                            |
|          W W AA RRRRR N N  II  N N GGGG !!!          |
|          W W A  A R R  NN N  II  NN N  G G  !!!          |
|          W W  A A  R R  N N  N  II  N N  N  G    !!!          |
|          W WW W  AAAAAA  RRRRR N  N N  II  N  N N  G  GGG !          |
|          WW  WW  A A  R R N NN  II  N NN  G G             |
|          W W  A A  R R  N N  II  N N GGGG !!!          |
|                                                                            |
|    Your highest band is occupied at some k-points! Unless you are       |
|    performing a calculation for an insulator or semiconductor, without |
|    unoccupied bands, you have included TOO FEW BANDS!! Please increase |
|    the parameter NBANDS in file 'INCAR' to ensure that the highest band |
|    is unoccupied at all k-points. It is always recommended to          |
|    include a few unoccupied bands to accelerate the convergence of       |
|    molecular dynamics runs (even for insulators or semiconductors).    |
|    Because the presence of unoccupied bands improves wavefunction       |
|    prediction, and helps to suppress 'band-crossings.'                   |
|    Following all k-points will be listed (with the Fermi weights of    |
|    the highest band given in paranthesis) ... :                         |
|                                                                            |
|                         1    ( 0.12000)                               |
|                         2    ( 0.05967)                               |
|                         4    ( 1.07820)                               |
|                         5    ( 0.22570)                               |
|                         8    ( 0.19163)                               |
|                         9    ( 1.07400)                               |
|                      10    ( 1.19262)                               |
|                      11    ( 0.11518)                               |
|                      12    ( 0.13723)                               |
|                      13    ( 0.37150)                               |
|                                                                            |
|    The total occupancy of band no. 448 is 0.72092 electrons ...    |
|                                                                            |
這個(gè)需要增大NBANDS到多少啊？按照VASP的說明手冊，NBANDS = NELECT/2 + NIONS/2 (non-spinpolarized) 這個(gè)值是足夠的啊，為什么會(huì)有這樣的警告？

贊一下

回復(fù)此樓

6樓2012-11-27 15:45:53

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

dfjk-123

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 29.9
散金: 757
帖子: 119
在線: 24.5小時(shí)
蟲號(hào): 665485
注冊: 2008-12-01
性別: MM
專業(yè): 理論和計(jì)算化學(xué)

引用回帖:

我后來試了NBANDS數(shù)不變，加大了K點(diǎn)，能量上升了。減少K點(diǎn)的話，warning就沒有了。
另外一個(gè)測試表明：增加了NBANDS數(shù)，K點(diǎn)數(shù)不變，能量下降了。在此基礎(chǔ)上又增加了K點(diǎn)數(shù)，能量在此基礎(chǔ)上又有下降。搞不清楚K點(diǎn)和NBANDS的關(guān)系到底是怎么一回事�。�
多謝回復(fù)！

贊一下

回復(fù)此樓

7樓2012-11-27 15:58:54

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

dfjk-123

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 29.9
散金: 757
帖子: 119
在線: 24.5小時(shí)
蟲號(hào): 665485
注冊: 2008-12-01
性別: MM
專業(yè): 理論和計(jì)算化學(xué)

送鮮花一朵

引用回帖:

8樓: Originally posted by 賀儀 at 2012-11-27 19:14:36
‘我后來試了NBANDS數(shù)不變，加大了K點(diǎn)，能量上升了。減少K點(diǎn)的話，warning就沒有了�！�
這明顯是K點(diǎn)不夠，warming沒有了是因?yàn)樵谀承㎏點(diǎn)的最高帶上已經(jīng)沒有占據(jù)了，但這個(gè)沒有是錯(cuò)誤的提示。

對于K點(diǎn)的夠不夠你 ...

NELECT和NIONS都在OUTCAR文件中可以查到，我沒有自己去數(shù)。另外，K點(diǎn)收斂性測試，我在1*1單胞的時(shí)候做的。變成3*3超胞后，我就將K點(diǎn)除以了3，有必要在超胞時(shí)還做K點(diǎn)測試嗎？謝謝回復(fù)！

贊一下

回復(fù)此樓

10樓2012-11-28 08:39:37

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

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

返回列表

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

[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3	TqlXswl 2026-03-16	7/350	2026-03-17 15:23 by TqlXswl
[考研] 283求調(diào)劑 +3	聽風(fēng)就是雨； 2026-03-16	3/150	2026-03-17 07:41 by 熱情沙漠
[考研] 一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑 +5	Liwangman 2026-03-15	5/250	2026-03-16 17:10 by 我的船我的海
[考研] 一志愿211 0703方向310分求調(diào)劑 +3	努力奮斗112 2026-03-15	3/150	2026-03-16 16:44 by houyaoxu
[考研] 318求調(diào)劑 +3	Yanyali 2026-03-15	3/150	2026-03-16 16:41 by houyaoxu
[考研] 0703化學(xué)調(diào)劑，求各位老師收留 +8	秋有木北 2026-03-14	8/400	2026-03-16 15:21 by 哦哦123
[考研] 070303 總分349求調(diào)劑 +3	LJY9966 2026-03-15	5/250	2026-03-16 14:24 by xwxstudy
[考研] 308求調(diào)劑 +3	是Lupa啊 2026-03-16	3/150	2026-03-16 10:07 by 求調(diào)劑zz
[考研] 326求調(diào)劑 +3	mlpqaz03 2026-03-15	3/150	2026-03-16 07:33 by Iveryant
[考研] 機(jī)械專碩調(diào)劑 +3	笨笨兔子 2026-03-12	3/150	2026-03-15 20:02 by 栗子粥?
[考研] 22408總分284求調(diào)劑 +3	InAspic 2026-03-13	3/150	2026-03-15 11:10 by zhq0425
[考研] 265求調(diào)劑 +4	威化餅07 2026-03-12	4/200	2026-03-14 17:23 by userper
[考研] 297求調(diào)劑 +4	學(xué)海漂泊 2026-03-13	4/200	2026-03-14 11:51 by 熱情沙漠
[考研] 331求調(diào)劑（0703有機(jī)化學(xué) +5	ZY-05 2026-03-13	6/300	2026-03-14 10:51 by Jy?
[考研] 266求調(diào)劑 +4	學(xué)員97LZgn 2026-03-13	4/200	2026-03-14 08:37 by zhukairuo
[考研] 281求調(diào)劑 +9	Koxui 2026-03-12	11/550	2026-03-13 20:50 by Koxui
[考研] 求調(diào)劑 +7	18880831720 2026-03-11	7/350	2026-03-13 16:10 by JourneyLucky
[考研] 工科材料085601 279求調(diào)劑 +8	困于星晨 2026-03-12	10/500	2026-03-13 15:42 by ms629
[考研] 308求調(diào)劑 +3	是Lupa啊 2026-03-12	3/150	2026-03-13 14:30 by 求調(diào)劑zz
[考研] 290求調(diào)劑 +3	ADT 2026-03-13	3/150	2026-03-13 10:19 by peike

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

[求助] 相同的vasp版本，不同機(jī)器上計(jì)算同一個(gè)體系，能量相差近10ev

» 猜你喜歡

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

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

» 本帖已獲得的紅花（最新10朵）

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

» 本帖已獲得的紅花（最新10朵）

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

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

[求助] 相同的vasp版本，不同機(jī)器上計(jì)算同一個(gè)體系，能量相差近10ev

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