| 4 | 1/1 | 返回列表 |
| 查看: 1315 | 回復(fù): 3 | ||
cuicui2397新蟲 (小有名氣)
|
[求助]
vasp中,結(jié)構(gòu)參數(shù)優(yōu)化時(shí),出現(xiàn)錯(cuò)誤,求助!謝謝!
|
|
侯老師的vasp教程中,針對復(fù)雜情況的結(jié)構(gòu)參數(shù)優(yōu)化出現(xiàn)問題,求助! 復(fù)雜的情況,以六角結(jié)構(gòu)Mg的晶格常數(shù)為例 1.先進(jìn)行一次體積保持不變的粒子馳豫的計(jì)算(ISIF=5) 2再把優(yōu)化得到的CONTCAR拷貝成POSCAR,進(jìn)行一次靜態(tài)的計(jì)算 3對a的值取10個(gè)左右的點(diǎn),每個(gè)點(diǎn)重復(fù)上面兩步,得到靜態(tài)計(jì)算下volume-Etot的關(guān)系。 這三部的腳步可以通過以下腳本run_cell來進(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)曲線后,通過狀態(tài)方程擬合得到平衡狀態(tài)下的體積,計(jì)算出上面腳本中的變量$i的值,并改變$i的循環(huán)值,再運(yùn)行run_cell計(jì)算一次,得到其他的結(jié)構(gòu)參數(shù)c和位置u。 我的問題是:在執(zhí)行sh run.cell命令時(shí),出現(xiàn)錯(cuò)誤提示是這樣的: run_cell: line 15: syntax error near unexpected token `newline' run_cell: line 15: `!' 請問這是什么錯(cuò)誤啊?謝謝!我改了改!附近的內(nèi)容,還是不行。 |
新蟲 (小有名氣)
銀蟲 (小有名氣)
新蟲 (著名寫手)
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿西南交通 專碩 材料355 本科雙非 求調(diào)劑 +4 | 西南交通專材355 2026-03-19 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 317求調(diào)劑 +4 | 申子申申 2026-03-19 | 8/400 |
|
|
[考研] 317求調(diào)劑 +4 | 申子申申 2026-03-19 | 9/450 |
|
|
[考研] 招收調(diào)劑碩士 +4 | lidianxing 2026-03-19 | 11/550 |
|
|
[考研] 296求調(diào)劑 +3 | www_q 2026-03-18 | 7/350 |
|
|
[考研] 08工學(xué)調(diào)劑 +4 | 用戶573181 2026-03-20 | 4/200 |
|
|
[考研] 294求調(diào)劑材料與化工專碩 +14 | 陌の森林 2026-03-18 | 14/700 |
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 18889395102 2026-03-18 | 4/200 |
|
|
[考研] 085601材料工程專碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
[考研] 材料與化工求調(diào)劑 +7 | 為學(xué)666 2026-03-16 | 7/350 |
|
|
[考研] 328求調(diào)劑,英語六級551,有科研經(jīng)歷 +4 | 生物工程調(diào)劑 2026-03-16 | 12/600 |
|
|
[考研] 化學(xué)工程321分求調(diào)劑 +15 | 大米飯! 2026-03-15 | 18/900 |
|
|
[考研] 268求調(diào)劑 +6 | 簡單點(diǎn)0 2026-03-17 | 6/300 |
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
[考研] 一志愿蘇州大學(xué)材料工程(085601)專碩有科研經(jīng)歷三項(xiàng)國獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級立項(xiàng) +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
[考研] 326求調(diào)劑 +4 | 諾貝爾化學(xué)獎(jiǎng)覬?/a> 2026-03-15 | 7/350 |
|
|
[考研] 321求調(diào)劑 +5 | 大米飯! 2026-03-15 | 5/250 |
|
|
[考研] 085601材料工程315分求調(diào)劑 +3 | yang_0104 2026-03-15 | 3/150 |
|
|
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3 | 麥田耕者 2026-03-14 | 3/150 |
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|