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

登錄注冊

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

24小時熱門版塊排行榜

返回列表

cuicui2397

新蟲 (小有名氣)

應(yīng)助: 5 (幼兒園)
金幣: 1611.8
紅花: 2
帖子: 211
在線: 130.4小時
蟲號: 2198831
注冊: 2012-12-21
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

[求助] vasp階段能測試求助，謝謝�。〖鼻�。

看的侯博士的教程，里面有測試Si的截斷能，但我自己去測試老是出現(xiàn)問題，求解��！謝謝。
POSCAR文件：
Si-Diamond
5.430
0.0 0.5 0.5
0.5 0.0 0.5
0.5 0.5 0.0
2
Direct
0.0 0.0 0.0
0.25 0.25 0.25

KPOINTS文件：
Automatic generation
0
Monhkorst-Pack
9 9 9
0.0 0.0 0.0

用來確定ENCUT的腳本程序為run_ecut，其內(nèi)容是：
#!/bin/sh
rm -rf WAVECAR
for i in 150 200 250 300 350 400
do
cat > INCAR <<!
SYSTEM = Si-Diamond
ENCUT = $i
ISTART = 0; ICHARG = 2
ISMEAR = -5
PREC = Accurate
!
echo "ENCUT = $i eV"; mpirun -np 2 vasp
E='grep "TOTEN" OUTCAR |tail-1|awk '{printf"%12.6f \n",$5}''
echo $i $E >>comment
done

同時選擇了贗勢文件POTCAR，把以上四個文件POSCAR、KPOINTS、run_ecut和POTCAR放入我的計算目錄下，在終端輸入sh.run_ecut命令，在生成的comment里是這樣的：
150 grep "TOTEN" OUTCAR |tail-1|awk {printf%12.6f\n,}
200 grep "TOTEN" OUTCAR |tail-1|awk {printf%12.6f\n,}
250 grep "TOTEN" OUTCAR |tail-1|awk {printf%12.6f\n,}
300 grep "TOTEN" OUTCAR |tail-1|awk {printf%12.6f\n,}
350 grep "TOTEN" OUTCAR |tail-1|awk {printf%12.6f\n,}
400 grep "TOTEN" OUTCAR |tail-1|awk {printf%12.6f\n,}
教程中第二列顯示的是能量值，請問這是什么問題��？謝謝�。�

回復(fù)此樓

» 猜你喜歡

一志愿西南交通專碩材料355 本科雙非求調(diào)劑已經(jīng)有3人回復(fù)
295復(fù)試調(diào)劑已經(jīng)有4人回復(fù)
298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑已經(jīng)有8人回復(fù)
材料學(xué)碩297已過四六級求調(diào)劑推薦已經(jīng)有4人回復(fù)
能源材料化學(xué)課題組招收碩士研究生8-10名已經(jīng)有14人回復(fù)
0856調(diào)劑，是學(xué)校就去已經(jīng)有8人回復(fù)
材料080500調(diào)劑求收留已經(jīng)有5人回復(fù)
271材料工程求調(diào)劑已經(jīng)有7人回復(fù)
265求調(diào)劑已經(jīng)有8人回復(fù)
085410人工智能專碩317求調(diào)劑（0854都可以）已經(jīng)有4人回復(fù)

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

vasp運行，第一離子步結(jié)束后自動停了已經(jīng)有7人回復(fù)
vasp中，做截斷能測試時，測試收斂的標(biāo)準已經(jīng)有10人回復(fù)
深圳超算中心vasp編譯及任務(wù)提交問題已經(jīng)有25人回復(fù)
關(guān)于VASP一些瑣事的求助（K點，吸附優(yōu)化等）已經(jīng)有12人回復(fù)
vasp并行測試時出錯。已經(jīng)有7人回復(fù)
vasp并行出錯，急于求助�。�！已經(jīng)有13人回復(fù)
vasp計算：優(yōu)化能收斂，但同樣的設(shè)置靜態(tài)計算不收斂已經(jīng)有3人回復(fù)
VASP計算時的幾個細節(jié)設(shè)置問題，希望有經(jīng)驗的人過來指點一下！已經(jīng)有23人回復(fù)
VASP 截斷能測試急求已經(jīng)有8人回復(fù)
【求助成功】vasp測試K點收斂性的腳本問題已經(jīng)有20人回復(fù)
【求助】VASP怎么做晶格參數(shù)的結(jié)構(gòu)優(yōu)化已經(jīng)有17人回復(fù)
【求助】求將vasp的輸出文件轉(zhuǎn)化成xyz坐標(biāo)的腳本已經(jīng)有11人回復(fù)
【求助】VASP計算電子態(tài)密度已經(jīng)有34人回復(fù)
【求助】vasp做表面計算時真空層的厚度已經(jīng)有13人回復(fù)
【求助】VASP能對體系加電場計算嗎？已經(jīng)有18人回復(fù)
【討論】vasp結(jié)構(gòu)優(yōu)化之前的測試已經(jīng)有11人回復(fù)

1樓 2013-11-14 20:16:44

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

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

甲斐之虎

銀蟲 (正式寫手)

應(yīng)助: 60 (初中生)
金幣: 302.3
散金: 220
紅花: 7
帖子: 970
在線: 465小時
蟲號: 2273045
注冊: 2013-02-03
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

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

★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
cuicui2397(fzx2008代發(fā)): 金幣+3 2013-11-14 21:25:06
sunyang1988: 金幣+1, 謝謝交流 2013-11-15 18:40:10
cuicui2397: 金幣+5, ★★★★★最佳答案 2013-11-23 10:10:40

那個，有可能是符號有問題，那個E='grep "TOTEN" OUTCAR |tail-1|awk '{printf"%12.6f \n",$5}''
是`（ESC下面那個鍵，數(shù)字1左面那個鍵），而不是'(單引號)，你試試看

贊一下(1人)

回復(fù)此樓

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

cuicui2397

別人笑我太瘋顛，我笑他人看不穿。不見五陵豪杰墓，無花無酒鋤作田。

2樓2013-11-14 20:30:12

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

甲斐之虎

銀蟲 (正式寫手)

應(yīng)助: 60 (初中生)
金幣: 302.3
散金: 220
紅花: 7
帖子: 970
在線: 465小時
蟲號: 2273045
注冊: 2013-02-03
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
sunyang1988: 金幣+1, 謝謝交流 2013-11-15 18:40:17
cuicui2397: 金幣+12, ★★★很有幫助 2013-11-23 10:11:24

贗勢是不是不一樣，你有沒有選正確的贗勢？侯老師的那個教程，有沒有說選什么贗勢？

贊一下(1人)

回復(fù)此樓

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

cuicui2397

別人笑我太瘋顛，我笑他人看不穿。不見五陵豪杰墓，無花無酒鋤作田。

4樓2013-11-14 23:06:04

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

cuicui2397

新蟲 (小有名氣)

應(yīng)助: 5 (幼兒園)
金幣: 1611.8
紅花: 2
帖子: 211
在線: 130.4小時
蟲號: 2198831
注冊: 2012-12-21
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

7樓: Originally posted by 甲斐之虎 at 2013-11-15 12:48:07
謝謝你的花，但是請別送花了，浪費金幣，留著以后提問題吧，有的用呢！
我說的不是這個意思，而是你采用了何種方法？PAW，UUSP？何種泛函91，PBE，LDA？這些不一樣，得出結(jié)果也不一樣...

嗯，呵呵，我自己選用的是paw_pbe贗勢，而侯的教程選用的是uspp、LDA。之前沒看清楚，給過來之后，結(jié)果和教程中一致了。謝謝你啊！不過我還想請教一下，針對具體體系，選擇贗勢是paw或是uspp，交換關(guān)聯(lián)函數(shù)選擇是LDA或是GGA，有什么根據(jù)嗎？謝謝！

贊一下(1人)

回復(fù)此樓

8樓2013-11-15 16:03:10

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

微6涼

鐵蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 202.8
帖子: 89
在線: 35.6小時
蟲號: 2724445
注冊: 2013-10-14
專業(yè): 金屬非晶態(tài)、準晶和納米晶

請問，我按照您這樣的做法，總會出現(xiàn)這樣一些結(jié)果Try `rm --help' for more information.
'un_ecut: line 4: syntax error near unexpected token `do
'un_ecut: line 4: `do 這是怎么回事

贊一下(1人)

回復(fù)此樓

14樓2015-08-03 11:50:15

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

普通回帖

cuicui2397

新蟲 (小有名氣)

應(yīng)助: 5 (幼兒園)
金幣: 1611.8
紅花: 2
帖子: 211
在線: 130.4小時
蟲號: 2198831
注冊: 2012-12-21
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

送紅花一朵

引用回帖:

2樓: Originally posted by 甲斐之虎 at 2013-11-14 20:30:12
那個，有可能是符號有問題，那個E='grep "TOTEN" OUTCAR |tail-1|awk '{printf"%12.6f \n",$5}''
是`（ESC下面那個鍵，數(shù)字1左面那個鍵），而不是'(單引號)，你試試看

你好，首先謝謝你，我把符號改過之后，第二列出現(xiàn)了能量值。不過和教程上差別很大，你覺得是什么原因造成的啊？我的結(jié)果是這樣的：
150 -10.744031
200 -10.795431
250 -10.822115
300 -10.836639
350 -10.840896
400 -10.842717

而教程中是這樣寫的：
150  -11.900655
200  -11.938864
250  -11.944599
300  -11.945248
350  -11.945503
400  -11.945622 總能變化在0.001eV左右就足夠了，在這個例子中，因此，選擇ENCUT=250eV。

我這樣算出的沒什么參考價值啊，求解，謝謝！

贊一下

回復(fù)此樓

3樓2013-11-14 21:18:24

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

cuicui2397

新蟲 (小有名氣)

應(yīng)助: 5 (幼兒園)
金幣: 1611.8
紅花: 2
帖子: 211
在線: 130.4小時
蟲號: 2198831
注冊: 2012-12-21
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

送紅花一朵

引用回帖:

4樓: Originally posted by 甲斐之虎 at 2013-11-14 23:06:04
贗勢是不是不一樣，你有沒有選正確的贗勢？侯老師的那個教程，有沒有說選什么贗勢？

針對這個例子來說，侯的教程里沒有說具體選什么贗勢，而Si有兩個贗勢，分別是Si 245和Si_h 380。之前侯的教程里說一般優(yōu)先選擇Si 245，所以我選擇了245那個，不知道對不？那我再試試380那個贗勢。

贊一下

回復(fù)此樓

5樓2013-11-15 09:54:46

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

cuicui2397

新蟲 (小有名氣)

應(yīng)助: 5 (幼兒園)
金幣: 1611.8
紅花: 2
帖子: 211
在線: 130.4小時
蟲號: 2198831
注冊: 2012-12-21
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

4樓: Originally posted by 甲斐之虎 at 2013-11-14 23:06:04
贗勢是不是不一樣，你有沒有選正確的贗勢？侯老師的那個教程，有沒有說選什么贗勢？

改成Si 380的贗勢，貌似更不對了，能量值是這樣的：
150  -11.407850
200  -10.938620
250  -10.825592
300  -10.813574
350  -10.830185
400  -10.849576
請問你一般是根據(jù)什么選擇贗勢�。恐x謝！

贊一下

回復(fù)此樓

6樓2013-11-15 10:29:02

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

甲斐之虎

銀蟲 (正式寫手)

應(yīng)助: 60 (初中生)
金幣: 302.3
散金: 220
紅花: 7
帖子: 970
在線: 465小時
蟲號: 2273045
注冊: 2013-02-03
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

★
sunyang1988: 金幣+1, 謝謝交流 2013-11-15 18:40:33

引用回帖:

6樓: Originally posted by cuicui2397 at 2013-11-15 10:29:02
改成Si 380的贗勢，貌似更不對了，能量值是這樣的：
150  -11.407850
200  -10.938620
250  -10.825592
300  -10.813574
350  -10.830185
400  -10.849576
請問你一般是根據(jù)什么選擇贗勢�。恐x謝！...

謝謝你的花，但是請別送花了，浪費金幣，留著以后提問題吧，有的用呢！
我說的不是這個意思，而是你采用了何種方法？PAW，UUSP？何種泛函91，PBE，LDA？這些不一樣，得出結(jié)果也不一樣

贊一下

回復(fù)此樓

別人笑我太瘋顛，我笑他人看不穿。不見五陵豪杰墓，無花無酒鋤作田。

7樓2013-11-15 12:48:07

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

忠義二郎

銅蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 88.5
帖子: 18
在線: 12.5小時
蟲號: 2669270
注冊: 2013-09-21
性別: GG
專業(yè): 地球內(nèi)部物理學(xué)

引用回帖:

請教一個問題，同樣是侯博士的這個測試Si的一個例子，我沒有UUSP的贗勢文件，所有就選了LDA、PBE的幾個Si的贗勢文件試一下，可是不論用哪個POTCAR，在comment中總是
150 0.0000
200 0.0000

請問是個什么情況？

贊一下

回復(fù)此樓

http://sdustzk.lofter.com

9樓2013-12-02 10:30:58

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

忠義二郎

銅蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 88.5
帖子: 18
在線: 12.5小時
蟲號: 2669270
注冊: 2013-09-21
性別: GG
專業(yè): 地球內(nèi)部物理學(xué)

請問這個例子中你的INCAR是怎么寫的？

贊一下

回復(fù)此樓

http://sdustzk.lofter.com

10樓2013-12-02 10:34:18

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

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

返回列表

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

[考研] 材料學(xué)碩297已過四六級求調(diào)劑推薦 +4	adaie 2026-03-19	4/200	2026-03-20 09:31 by yuncha
[考研] 271材料工程求調(diào)劑 +7	.6lL 2026-03-18	7/350	2026-03-20 09:10 by xingguangj
[考研] 317求調(diào)劑 +3	申子申申 2026-03-19	7/350	2026-03-20 08:53 by 申子申申
[考研] 0703化學(xué)調(diào)劑，六級已過，有科研經(jīng)歷 +12	曦熙兮 2026-03-15	12/600	2026-03-19 19:42 by maocaozhuxi
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +5	楓橋ZL 2026-03-18	7/350	2026-03-19 14:52 by 功夫瘋狂
[考研] 0703化學(xué)調(diào)劑，求各位老師收留 +10	秋有木北 2026-03-14	10/500	2026-03-19 05:52 by anny19840123
[考研] 材料工程專碩調(diào)劑 +5	204818@lcx 2026-03-17	6/300	2026-03-18 22:55 by 204818@lcx
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +3	生物工程調(diào)劑 2026-03-17	7/350	2026-03-18 20:41 by Wangjingyue
[考研] 085600材料與化工 +5	安全上岸！ 2026-03-16	5/250	2026-03-18 15:33 by cmz0325
[考研] 070300化學(xué)319求調(diào)劑 +6	錦鯉0909 2026-03-17	6/300	2026-03-18 13:22 by Iveryant
[考研] 0854，計算機類招收調(diào)劑 +3	胡辣湯放糖 2026-03-15	6/300	2026-03-18 12:09 by 上岸上岸……..
[考研] 278求調(diào)劑 +5	煙火先于春 2026-03-17	5/250	2026-03-18 08:43 by 星空星月
[考研] 考研求調(diào)劑 +3	橘頌. 2026-03-17	4/200	2026-03-17 21:43 by 有只貍奴
[考研] 268求調(diào)劑 +8	一定有學(xué)上- 2026-03-14	9/450	2026-03-17 17:47 by laoshidan
[考研] 326求調(diào)劑 +5	上岸的小葡 2026-03-15	6/300	2026-03-17 17:26 by ruiyingmiao
[考研] 材料與化工專碩調(diào)劑 +5	heming3743 2026-03-16	5/250	2026-03-17 14:03 by 勇敢太監(jiān)王公公
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3	尖角小荷 2026-03-16	6/300	2026-03-16 23:21 by 尖角小荷
[考研] 333求調(diào)劑 +3	文思客 2026-03-16	7/350	2026-03-16 18:21 by 文思客
[考研] 318求調(diào)劑 +3	Yanyali 2026-03-15	3/150	2026-03-16 16:41 by houyaoxu
[考研] 26考研一志愿中國石油大學(xué)(華東)305分求調(diào)劑 +3	嘉年新程 2026-03-15	3/150	2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿

24小時熱門版塊排行榜

[求助] vasp階段能測試求助，謝謝�。〖鼻�。

» 猜你喜歡

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

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

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

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

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

[求助] vasp階段能測試求助，謝謝�。〖鼻�。

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