stevencong

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 8
帖子: 18
在線: 8.4小時
蟲號: 1732947
注冊: 2012-04-02
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

[求助] 有關(guān)磁性體系計算能帶時遇到的問題

把全部金幣貼進來了，求大牛們指導(dǎo)啊~~

現(xiàn)在在做磁性納米管的能帶計算，但是不太清楚計算方法。我的計算過程是這樣的：

1. 計算體系的Charge Density，輸入文件如下（其中MAGMOM是從弛豫后的OUTCAR文件中抄下來的）：

SYSTEM = Nanotube
ENCUT = 800
ISTART = 0; ICHARG = 2
ISMEAR = -5
ISPIN = 2
MAGMOM = 0.03 0.03 -0.01 0.007 0.03 0.03 0.038 0.038 0.006 -0.01 0.039 0.038 0.03 0.03 -0.01 0.006 0.03 0.03 0.039 0.039 0.007 -0.01 0.038 0.038 0.03 0.03 -0.01 0.006 0.03 0.03 0.038 0.038 0.006 -0.01 0.038 0.039
PREC = Accurate

2. 再做BS計算，輸入文件為：

SYSTEM = Nanotube
ENCUT = 800
ISTART = 1; ICHARG = 11
NBANDS = 120
ISPIN = 2
MAGMOM = 0.03 0.03 -0.01 0.007 0.03 0.03 0.038 0.038 0.006 -0.01 0.039 0.038 0.03 0.03 -0.01 0.006 0.03 0.03 0.039 0.039 0.007 -0.01 0.038 0.038 0.03 0.03 -0.01 0.006 0.03 0.03 0.038 0.038 0.006 -0.01 0.038 0.039
PREC = Accurate

但是每次在計算BS時VASP剛開始運行就掛掉了，out文件輸出：

running on 16 nodes
distr:  one band on 1 nodes, 16 groups
vasp.5.2.2 15Apr09 complex
POSCAR found type information on POSCAR  C
POSCAR found :  1 types and 36 ions

-----------------------------------------------------------------------------
|                                                                            |
|  ADVICE TO THIS USER RUNNING 'VASP/VAMP' (HEAR YOUR MASTER'S VOICE ...):  |
|                                                                            |
|    You have a (more or less) 'large supercell' and for larger cells    |
|    it might be more efficient to use real space projection operators    |
|    So try LREAL= Auto  in the INCAR file.                            |
|    Mind: If you want to do a very accurate calculations keep the       |
|    reciprocal projection scheme       (i.e. LREAL=.FALSE.)          |
|                                                                            |
-----------------------------------------------------------------------------

LDA part: xc-table for Pade appr. of Perdew
WARNING: stress and forces are not correct
POSCAR, INCAR and KPOINTS ok, starting setup
FFT: planning ...(          1 )

然后，然后就什么都沒有了。。。不知道哪位能解決這個問題��？？？在此先謝過了~

回復(fù)此樓

» 猜你喜歡

310求調(diào)劑已經(jīng)有5人回復(fù)
材料與化工272求調(diào)劑已經(jīng)有17人回復(fù)
材料工程專碩求調(diào)劑已經(jīng)有3人回復(fù)
337求調(diào)劑已經(jīng)有5人回復(fù)
材料與化工328分調(diào)劑已經(jīng)有9人回復(fù)
318一志愿吉林大學(xué)生物與醫(yī)藥求調(diào)劑已經(jīng)有5人回復(fù)
一志愿南開大學(xué)0710生物學(xué)359求調(diào)劑已經(jīng)有4人回復(fù)
求調(diào)劑已經(jīng)有6人回復(fù)
085600材料與化工調(diào)劑已經(jīng)有4人回復(fù)
327求調(diào)劑已經(jīng)有4人回復(fù)

» 本主題相關(guān)商家推薦: (我也要在這里推廣)

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

計算晶體的能帶譜圖發(fā)現(xiàn)了一個很奇怪的現(xiàn)象已經(jīng)有9人回復(fù)
【專家在線答疑貼之三】雜化泛函計算FAQ 已經(jīng)有28人回復(fù)
【求助】各位大俠！小女子求助鈣鈦礦結(jié)構(gòu)計算磁性的問題？已經(jīng)有36人回復(fù)
【求助】VASP計算體系性能沒有f能帶電子已經(jīng)有24人回復(fù)
【求助】castep 優(yōu)化鐵時 initial spin該如何設(shè)置已經(jīng)有26人回復(fù)
【求助】關(guān)于擬合得到體彈性模量&自旋極化&LDA+U等問題，望大家指點指點已經(jīng)有4人回復(fù)
【原創(chuàng)】DFT計算體會已經(jīng)有20人回復(fù)
【其他】在線答疑:VASP的問題已經(jīng)有314人回復(fù)
【求助】vasp noncollinear計算能帶的時候只能算spin-unpolarized嗎已經(jīng)有12人回復(fù)

1樓 2012-07-11 14:51:31

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

stevencong

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 8
帖子: 18
在線: 8.4小時
蟲號: 1732947
注冊: 2012-04-02
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

2樓: Originally posted by mousekingadv at 2012-07-11 15:48:06
這不是掛掉，這是條件太苛刻以致于一步都跑不出來……
ENCUT=800...手冊上第一個例子就是講的測試ENCUT吧...無論如何800太高了。
NBANDS并不需要手動提高，如果不關(guān)心高能空帶的結(jié)構(gòu)、也不需要計算光學(xué)性質(zhì)的話。 ...

結(jié)果又變成這樣了：

running on 8 nodes
distr:  one band on 1 nodes, 8 groups
vasp.5.2.2 15Apr09 complex
POSCAR found type information on POSCAR  C
POSCAR found :  1 types and 36 ions
LDA part: xc-table for Pade appr. of Perdew
WARNING: stress and forces are not correct
POSCAR, INCAR and KPOINTS ok, starting setup
WARNING: small aliasing (wrap around) errors must be expected
FFT: planning ...(          1 )
reading WAVECAR
WARNING: dimensions on CHGCAR file are different
ERROR: charge density could not be read from file CHGCAR for ICHARG>10
ERROR: charge density could not be read from file CHGCAR for ICHARG>10
ERROR: charge density could not be read from file CHGCAR for ICHARG>10
ERROR: charge density could not be read from file CHGCAR for ICHARG>10
ERROR: charge density could not be read from file CHGCAR for ICHARG>10
ERROR: charge density could not be read from file CHGCAR for ICHARG>10
ERROR: charge density could not be read from file CHGCAR for ICHARG>10
ERROR: charge density could not be read from file CHGCAR for ICHARG>10

我兩次計算所用的POSCAR應(yīng)該一樣啊，為什么Dimension會不一樣呢？

贊一下

回復(fù)此樓

4樓2012-07-12 00:19:55

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

查看全部 8 個回答

mousekingadv

金蟲 (小有名氣)

應(yīng)助: 51 (初中生)
金幣: 1141.8
紅花: 1
帖子: 140
在線: 28.7小時
蟲號: 1721002
注冊: 2012-03-27
專業(yè): 理論和計算化學(xué)

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
fzx2008: 金幣+3, 謝謝指教！ 2012-07-11 16:26:13
stevencong: 金幣+20, ★★★很有幫助 2012-07-12 14:34:45

這不是掛掉，這是條件太苛刻以致于一步都跑不出來……
ENCUT=800...手冊上第一個例子就是講的測試ENCUT吧...無論如何800太高了。
NBANDS并不需要手動提高，如果不關(guān)心高能空帶的結(jié)構(gòu)、也不需要計算光學(xué)性質(zhì)的話。
另外KPOINTS不能過多。樓主是不是沿著高對稱路徑“老板，每條路徑給我來50個K點”？先用較低的標(biāo)準(zhǔn)測試成功了，再在可接受的世界范圍內(nèi)提高精度。
LREAL=AUTO可以考慮寫上。

贊一下

回復(fù)此樓

2樓2012-07-11 15:48:06

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

stevencong

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 8
帖子: 18
在線: 8.4小時
蟲號: 1732947
注冊: 2012-04-02
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

這樣子，我先試試。謝謝你的回復(fù)！

贊一下

回復(fù)此樓

3樓2012-07-11 16:23:42

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

wutianhuan

鐵桿木蟲 (著名寫手)

應(yīng)助: 3 (幼兒園)
金幣: 6697.4
散金: 12
紅花: 3
帖子: 1723
在線: 216.2小時
蟲號: 284387
注冊: 2006-10-09
性別: MM
專業(yè): 金屬材料的合金相、相變及

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

★
感謝參與，應(yīng)助指數(shù) +1
xueht987: 金幣+1, 鼓勵交流！ 2012-07-13 09:17:30

ERROR: charge density could not be read from file CHGCAR for ICHARG>10
------------------------手冊里ICHARG有1、2和10一共3種選擇，你寫了11，所以程序不認識了，呵呵

贊一下

回復(fù)此樓

5樓2012-07-12 11:25:28

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

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

[考研] 085600，材料與化工321分求調(diào)劑 +10	大饞小子 2026-03-28	10/500	2026-03-29 23:35 by 飛行日記西
[考研] 307求調(diào)劑 +9	超級伊昂大王 2026-03-24	10/500	2026-03-29 20:26 by 永字號
[考研] 277跪求調(diào)劑 +6	1915668 2026-03-27	10/500	2026-03-29 16:03 by 王亮_大連醫(yī)科大
[考研] 297求調(diào)劑 +11	田洪有 2026-03-26	11/550	2026-03-29 13:14 by peike
[考研] 求調(diào)劑，一志愿南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +7	@taotao 2026-03-29	7/350	2026-03-29 12:03 by longlotian
[考研] 本科雙非材料，跨考一志愿華電085801電氣，283求調(diào)劑，任何專業(yè)都可以 +6	芝士雪baoo 2026-03-28	8/400	2026-03-29 08:16 by 松花缸1201
[考研] 一志愿太原理工安全工程300分，求調(diào)劑 +5	0857求調(diào)劑. 2026-03-24	6/300	2026-03-28 22:04 by zhq0425
[考研] 求調(diào)劑 +7	爭取九點睡 2026-03-28	8/400	2026-03-28 21:07 by 爭取九點睡
[考研] 一志愿廈門大學(xué)化學(xué)學(xué)碩307求調(diào)劑 +10	y7czhao 2026-03-26	10/500	2026-03-28 14:23 by 唐沐兒
[考研] 339求調(diào)劑，想調(diào)回江蘇 +6	烤麥芽 2026-03-27	8/400	2026-03-28 10:40 by 烤麥芽
[考研] 085405 考的11408求各位老師帶走 +3	Qiu學(xué)ing 2026-03-28	3/150	2026-03-28 09:19 by 樂呵呵的追夢人
[考研] 292求調(diào)劑 +14	鵝鵝鵝額額額額?/a> 2026-03-25	15/750	2026-03-28 08:45 by WYUMater
[考研] 0703化學(xué)求調(diào)劑，各位老師看看我！��！ +5	祁祺祺 2026-03-25	5/250	2026-03-27 21:44 by 東方豬豬
[考研] 266分求材料化工冶金礦業(yè)等專業(yè)的調(diào)劑 +4	哇呼哼呼哼 2026-03-26	4/200	2026-03-27 17:02 by zhyzzh
[考研] 284求調(diào)劑 +11	junqihahaha 2026-03-26	12/600	2026-03-27 04:37 by wxiongid
[考研] 化學(xué)調(diào)劑一志愿上海交通大學(xué)336分-本科上海211 +4	小魚愛有機 2026-03-25	4/200	2026-03-26 10:19 by aa331100
[考研] 0854AI CV方向招收調(diào)劑 +4	章小魚567 2026-03-23	4/200	2026-03-25 17:04 by CoderLoser
[考研] 各位老師您好：本人初試372分 +5	jj涌77 2026-03-25	6/300	2026-03-25 14:15 by mapenggao
[考研] 282求調(diào)劑 +3	wcq131415 2026-03-24	3/150	2026-03-25 12:16 by userper
[考研] 一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑 +3	開開心心沒煩惱 2026-03-23	4/200	2026-03-24 00:06 by 開開心心沒煩惱