| 2 | 1/1 | 返回列表 |
| 查看: 956 | 回復(fù): 7 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
[求助]
ENCUT測(cè)試 能量不變
|
||
|
我用的測(cè)試ENCUT腳本為: 1 #!/bin/sh 2 rm -rf WAVECAR 3 for i in 0.1 1 5 4 #1 5 10 15 100 150 200 250 300 350 400 450 500 550 600 5 do 6 rm -rf WAVECAR 7 cat < INCAR << ! 8 SYSTEM = bcc Fe 9 ENCUT = $i 10 ISTART = 0 ; ICHARG = 2 11 ISMEAR = -5 12 PREC = Accurate 13 ! 14 echo "ENCUT = $i eV" 15 mpirun -hostfile ~/host2 -np 8 vasp 16 E=$(grep "TOTEN" OUTCAR | tail -1| awk '{printf "%15.9f \n", $5}') 17 echo $i $E >> comment 18 done 出來(lái)的comment 1 0.1 -15.435962000 2 1 -15.435962000 3 5 -15.435962000 4 1 -15.435962000 5 5 -15.435962000 6 10 -15.435962000 7 15 -15.435962000 8 100 -15.435962000 9 150 -15.435962000 10 200 -15.435962000 11 250 -15.435962000 12 300 -15.435962000 13 350 -15.435962000 14 400 -15.435962000 15 450 -15.435962000 16 500 -15.435962000 17 550 -15.435962000 18 600 -15.435962000 第一列是行號(hào),第二列是ENCUT,第三列是輸出的結(jié)果。 但是都一樣,不知道哪里錯(cuò)了,求助呀。 |

金蟲(chóng) (小有名氣)
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|