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

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

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

返回列表

cuicui2397

新蟲 (小有名氣)

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

[求助] vasp階段能測(cè)試求助，謝謝啊！急求。

看的侯博士的教程，里面有測(cè)試Si的截?cái)嗄�，但我自己去測(cè)試?yán)鲜浅霈F(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的腳本程序?yàn)閞un_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

同時(shí)選擇了贗勢(shì)文件POTCAR，把以上四個(gè)文件POSCAR、KPOINTS、run_ecut和POTCAR放入我的計(jì)算目錄下，在終端輸入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,}
教程中第二列顯示的是能量值，請(qǐng)問這是什么問題�。恐x謝�。�

回復(fù)此樓

» 猜你喜歡

一志愿中南化學(xué)337求調(diào)劑已經(jīng)有4人回復(fù)
281求調(diào)劑（0805）已經(jīng)有23人回復(fù)
265求調(diào)劑已經(jīng)有9人回復(fù)
279分求調(diào)劑一志愿211 已經(jīng)有10人回復(fù)
招收調(diào)劑碩士已經(jīng)有12人回復(fù)
296求調(diào)劑已經(jīng)有8人回復(fù)
招收博士1-2人已經(jīng)有3人回復(fù)
085700資源與環(huán)境308求調(diào)劑已經(jīng)有11人回復(fù)
274求調(diào)劑已經(jīng)有8人回復(fù)
290求調(diào)劑已經(jīng)有4人回復(fù)

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

vasp運(yùn)行，第一離子步結(jié)束后自動(dòng)停了已經(jīng)有7人回復(fù)
vasp中，做截?cái)嗄軠y(cè)試時(shí)，測(cè)試收斂的標(biāo)準(zhǔn) 已經(jīng)有10人回復(fù)
深圳超算中心vasp編譯及任務(wù)提交問題已經(jīng)有25人回復(fù)
關(guān)于VASP一些瑣事的求助（K點(diǎn)，吸附優(yōu)化等）已經(jīng)有12人回復(fù)
vasp并行測(cè)試時(shí)出錯(cuò)。已經(jīng)有7人回復(fù)
vasp并行出錯(cuò)，急于求助�。�！已經(jīng)有13人回復(fù)
vasp計(jì)算：優(yōu)化能收斂，但同樣的設(shè)置靜態(tài)計(jì)算不收斂已經(jīng)有3人回復(fù)
VASP計(jì)算時(shí)的幾個(gè)細(xì)節(jié)設(shè)置問題，希望有經(jīng)驗(yàn)的人過來指點(diǎn)一下！已經(jīng)有23人回復(fù)
VASP 截?cái)嗄軠y(cè)試急求已經(jīng)有8人回復(fù)
【求助成功】vasp測(cè)試K點(diǎn)收斂性的腳本問題已經(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計(jì)算電子態(tài)密度已經(jīng)有34人回復(fù)
【求助】vasp做表面計(jì)算時(shí)真空層的厚度已經(jīng)有13人回復(fù)
【求助】VASP能對(duì)體系加電場計(jì)算嗎？已經(jīng)有18人回復(fù)
【討論】vasp結(jié)構(gòu)優(yōu)化之前的測(cè)試已經(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小時(shí)
蟲號(hào): 2273045
注冊(cè): 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

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

贊一下(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小時(shí)
蟲號(hào): 2273045
注冊(cè): 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

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

贊一下(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小時(shí)
蟲號(hào): 2198831
注冊(cè): 2012-12-21
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

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

嗯，呵呵，我自己選用的是paw_pbe贗勢(shì)，而侯的教程選用的是uspp、LDA。之前沒看清楚，給過來之后，結(jié)果和教程中一致了。謝謝你啊！不過我還想請(qǐng)教一下，針對(duì)具體體系，選擇贗勢(shì)是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小時(shí)
蟲號(hào): 2724445
注冊(cè): 2013-10-14
專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶

請(qǐng)問，我按照您這樣的做法，總會(huì)出現(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小時(shí)
蟲號(hào): 2198831
注冊(cè): 2012-12-21
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

送紅花一朵

引用回帖:

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

你好，首先謝謝你，我把符號(hào)改過之后，第二列出現(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左右就足夠了，在這個(gè)例子中，因此，選擇ENCUT=250eV。

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

贊一下

回復(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小時(shí)
蟲號(hào): 2198831
注冊(cè): 2012-12-21
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

送紅花一朵

引用回帖:

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

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

贊一下

回復(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小時(shí)
蟲號(hào): 2198831
注冊(cè): 2012-12-21
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

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

贊一下

回復(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小時(shí)
蟲號(hào): 2273045
注冊(cè): 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的贗勢(shì)，貌似更不對(duì)了，能量值是這樣的：
150  -11.407850
200  -10.938620
250  -10.825592
300  -10.813574
350  -10.830185
400  -10.849576
請(qǐng)問你一般是根據(jù)什么選擇贗勢(shì)��？謝謝！...

謝謝你的花，但是請(qǐng)別送花了，浪費(fèi)金幣，留著以后提問題吧，有的用呢！
我說的不是這個(gè)意思，而是你采用了何種方法？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小時(shí)
蟲號(hào): 2669270
注冊(cè): 2013-09-21
性別: GG
專業(yè): 地球內(nèi)部物理學(xué)

引用回帖:

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

請(qǐng)問是個(gè)什么情況？

贊一下

回復(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小時(shí)
蟲號(hào): 2669270
注冊(cè): 2013-09-21
性別: GG
專業(yè): 地球內(nèi)部物理學(xué)

請(qǐng)問這個(gè)例子中你的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ā)表

[考研] 工科材料085601 279求調(diào)劑 +6	困于星晨 2026-03-17	8/400	2026-03-20 11:24 by kkcoco25
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +7	Doleres 2026-03-19	7/350	2026-03-20 10:39 by EBSD
[考研] 一志愿南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩 +4	@taotao 2026-03-20	4/200	2026-03-20 10:37 by EBSD
[考研] 271材料工程求調(diào)劑 +7	.6lL 2026-03-18	7/350	2026-03-20 09:10 by xingguangj
[考研] 288求調(diào)劑 +15	于海海海海 2026-03-19	15/750	2026-03-19 22:41 by 學(xué)員8dgXkO
[考研] 0703化學(xué)調(diào)劑，六級(jí)已過，有科研經(jīng)歷 +12	曦熙兮 2026-03-15	12/600	2026-03-19 19:42 by maocaozhuxi
[考研] 一志愿福大288有機(jī)化學(xué)，求調(diào)劑 +3	小木蟲200408204 2026-03-18	3/150	2026-03-19 13:31 by houyaoxu
[考研] 328求調(diào)劑，英語六級(jí)551，有科研經(jīng)歷 +4	生物工程調(diào)劑 2026-03-16	12/600	2026-03-19 11:10 by 生物工程調(diào)劑
[考研] 一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑 +10	Liwangman 2026-03-15	10/500	2026-03-19 10:25 by 無際的草原
[考研] 304求調(diào)劑 +6	司空. 2026-03-18	6/300	2026-03-18 23:03 by 星空星月
[考研] 材料專業(yè)求調(diào)劑 +5	hanamiko 2026-03-18	5/250	2026-03-18 20:19 by 星空星月
[考研] 【同濟(jì)軟件】軟件（085405）考研求調(diào)劑 +3	2026eternal 2026-03-18	3/150	2026-03-18 19:09 by 搏擊518
[考研] 一志愿武理材料305分求調(diào)劑 +5	想上岸的鯉魚 2026-03-18	6/300	2026-03-18 17:53 by 無際的草原
[考研] 297求調(diào)劑 +8	戲精丹丹丹 2026-03-17	8/400	2026-03-18 14:30 by laoshidan
[考研] 070300化學(xué)319求調(diào)劑 +6	錦鯉0909 2026-03-17	6/300	2026-03-18 13:22 by Iveryant
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6	薛云鵬 2026-03-15	6/300	2026-03-17 11:05 by 學(xué)員h26Tkc
[考研] 機(jī)械專碩325，尋找調(diào)劑院校 +3	y9999 2026-03-15	5/250	2026-03-16 19:58 by y9999
[考研] 326求調(diào)劑 +3	mlpqaz03 2026-03-15	3/150	2026-03-16 07:33 by Iveryant
[考研] 327求調(diào)劑 +6	拾光任染 2026-03-15	11/550	2026-03-15 22:47 by 拾光任染
[考研] 復(fù)試調(diào)劑 +3	呼呼？~+123456 2026-03-14	3/150	2026-03-14 16:53 by WTUChen

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

[求助] vasp階段能測(cè)試求助，謝謝啊！急求。

» 猜你喜歡

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

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

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

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

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

[求助] vasp階段能測(cè)試求助，謝謝啊！急求。

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