| 4 | 1/1 | 返回列表 |
| 查看: 1316 | 回復(fù): 3 | ||
cuicui2397新蟲(chóng) (小有名氣)
|
[求助]
vasp中,結(jié)構(gòu)參數(shù)優(yōu)化時(shí),出現(xiàn)錯(cuò)誤,求助!謝謝!
|
|
侯老師的vasp教程中,針對(duì)復(fù)雜情況的結(jié)構(gòu)參數(shù)優(yōu)化出現(xiàn)問(wèn)題,求助! 復(fù)雜的情況,以六角結(jié)構(gòu)Mg的晶格常數(shù)為例 1.先進(jìn)行一次體積保持不變的粒子馳豫的計(jì)算(ISIF=5) 2再把優(yōu)化得到的CONTCAR拷貝成POSCAR,進(jìn)行一次靜態(tài)的計(jì)算 3對(duì)a的值取10個(gè)左右的點(diǎn),每個(gè)點(diǎn)重復(fù)上面兩步,得到靜態(tài)計(jì)算下volume-Etot的關(guān)系。 這三部的腳步可以通過(guò)以下腳本run_cell來(lái)進(jìn)行,內(nèi)容如下: #!/bin/sh rm WAVECAR for i in 5.20 5.23 5.26 5.29 5.32 5.35 5.38 5.41 5.44 5.47 do cat > POSCAR << ! ureadopef $i 0 -1 0 0.8660254 0.5 0 0 0 1.6230530 2 Direct 0.666666667 0.333333333 0.75 0.333333333 0.666666667 0.25 ! cp INCAR.relax INCAR echo "a=$i angestrom"; mpiexec -np 8 vasp cp CONTCAR POSCAR cp INCAR.static INCAR echo "a= $i angstrom ";mpiexec -np 8 vasp E=`grep "TOTEN" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5}'` V=`grep "volume" OUTCAR | tail -1 | awk '{printf "%12.4f \n", $5}'` echo $V $E >>EtVo.dat done 4.采用狀態(tài)方程擬合得到的平衡態(tài)下的體積,體彈性模量。 5.在該體積下,重復(fù)上面第2步和第3步,得到平衡狀態(tài)下的其他晶胞參數(shù)。這一步也就是:在的到了E(V)曲線后,通過(guò)狀態(tài)方程擬合得到平衡狀態(tài)下的體積,計(jì)算出上面腳本中的變量$i的值,并改變$i的循環(huán)值,再運(yùn)行run_cell計(jì)算一次,得到其他的結(jié)構(gòu)參數(shù)c和位置u。 我的問(wèn)題是:在執(zhí)行sh run.cell命令時(shí),出現(xiàn)錯(cuò)誤提示是這樣的: run_cell: line 15: syntax error near unexpected token `newline' run_cell: line 15: `!' 請(qǐng)問(wèn)這是什么錯(cuò)誤。恐x謝!我改了改!附近的內(nèi)容,還是不行。 |
新蟲(chóng) (小有名氣)
銀蟲(chóng) (小有名氣)
新蟲(chóng) (著名寫(xiě)手)
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 281求調(diào)劑(0805) +12 | 煙汐憶海 2026-03-16 | 23/1150 |
|
|---|---|---|---|---|
|
[考研] 一志愿西南交通 專(zhuān)碩 材料355 本科雙非 求調(diào)劑 +4 | 西南交通專(zhuān)材355 2026-03-19 | 4/200 |
|
|
[考研] 298-一志愿中國(guó)農(nóng)業(yè)大學(xué)-求調(diào)劑 +8 | 手機(jī)用戶(hù) 2026-03-17 | 8/400 |
|
|
[考研] 0856調(diào)劑,是學(xué)校就去 +7 | sllhht 2026-03-19 | 8/400 |
|
|
[考研] 求調(diào)劑 +3 | 暗涌afhb 2026-03-16 | 3/150 |
|
|
[考研] 081700化工學(xué)碩調(diào)劑 +3 | 【1】 2026-03-16 | 3/150 |
|
|
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +3 | 葵梓衛(wèi)隊(duì) 2026-03-18 | 5/250 |
|
|
[考研] 286求調(diào)劑 +6 | lemonzzn 2026-03-16 | 10/500 |
|
|
[考研] 324分 085600材料化工求調(diào)劑 +3 | llllkkkhh 2026-03-18 | 3/150 |
|
|
[考研] 085600材料與化工求調(diào)劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|
|
[考研] 330求調(diào)劑 +3 | 小材化本科 2026-03-18 | 3/150 |
|
|
[考研] 311求調(diào)劑 +4 | 冬十三 2026-03-18 | 4/200 |
|
|
[考研] 一志愿西南交大,求調(diào)劑 +4 | 材化逐夢(mèng)人 2026-03-18 | 4/200 |
|
|
[考研]
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級(jí)食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
[考博] 26申博 +4 | 八6八68 2026-03-16 | 4/200 |
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
[考研] 297求調(diào)劑 +4 | 學(xué)海漂泊 2026-03-13 | 4/200 |
|