| 5 | 1/1 | 返回列表 |
| 查看: 2605 | 回復(fù): 14 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
[求助]
VASP中腳本程序問題
|
|||
|
我初學(xué)VASP,看侯老師的教程中有時(shí)候需要編寫腳本程序,比如優(yōu)化晶格常數(shù)有如下腳本:#!/bin/sh rm wAVECAR for i in 3.80 3.85 3.90 3.95 4.00 4.05 4.10 4.15 4.20 4.25 4.30 do cat>POSCAR< Al-fcc $i 0.0 0.5 0.5 0.5 0.0 0.5 0.5 0.5 0.0 1 Direct 0.0 0.0 0.0 ! echo "a=$i angstrom";time vasp E='grep "TOTEN" OUTCAR |tail -1|awk '{printf "%12.6f\n",$5}'' V='grep "volume" OUTCAR |tail -1|awk '{prinf "%12.4f\n",$5}'' echo $V $E>>EtVo.dat done 我想問的是這個(gè)是寫在POSCAR文件里嗎?如果不是的話就是另寫一個(gè)腳本吧?我們學(xué)校的高性能計(jì)算中心利用PBS提交作業(yè),請(qǐng)問如果另寫一個(gè)腳本的話我的PBS文件中應(yīng)該加入什么命令呢?謝謝了 |
第一原理資料貼 |
金蟲 (正式寫手)
|
這種小程序是隨便建個(gè)文檔寫入,寫完保存后變成執(zhí)行文件就好 最后,關(guān)于pbs上面如何提交作業(yè)的話我是這樣寫的,你可以參考下 #!/bin/sh #PBS -N vasp #PBS -d ./ #PBS -o ./ #PBS -e ./ #PBS -l nodes=1:ppn=4 #PBS -l walltime=10:00:00 #PBS -l mem= #PBS -q serique #PBS -V ulimit -s unlimited source /opt/intel/Compiler/11.1/046/bin/intel64/iccvars_intel64.sh source /opt/intel/Compiler/11.1/046/bin/intel64/ifortvars_intel64.sh source /opt/intel/mkl/10.2.2.025/tools/environment/mklvarsem64t.sh source /opt/intel/impi/3.2.2/bin64/mpivars.sh EXEC=/***/***/vasp_5.2 ################################################################ rm WAVECAR for i in 250 300 350 400 450 500 550 600 do cat > INCAR <<! System = Si ENCUT = $i ISTART = 0 ICHARG = 2 ISMEAR = 0 SIGMA = 0.1 PREC = Accurate ! echo "ENCUT = $i eV" ; /***/***/ vasp_5.2(這里是寫vasp的路徑,你自己是什么寫什么) E=`grep "TOTEN" OUTCAR | tail -1 | awk '{printf "%12.6f\n",$5}'` KP=`grep "irreducible" OUTCAR | tail -1 | awk '{printf "%5i \n", $2 }'` echo $i $KP $E >> comment done |
金蟲 (正式寫手)
木蟲 (小有名氣)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 275求調(diào)劑 +11 | jjjjjjjjjjl 2026-03-27 | 11/550 |
|
|---|---|---|---|---|
|
[考研] 一志愿北京科技大學(xué)085601材料工程英一數(shù)二初試總分335求調(diào)劑 +7 | 雙馬尾痞老板2 2026-04-01 | 7/350 |
|
|
[考研] 271求調(diào)劑 +14 | 勒布朗@ 2026-03-31 | 19/950 |
|
|
[考博] 26年申博 +3 | staryer 2026-03-30 | 4/200 |
|
|
[考研] 材料專碩322分 +10 | 哈哈哈吼吼吼哈 2026-04-01 | 10/500 |
|
|
[考研] 310分求調(diào)劑 +4 | 成功上岸wang 2026-04-01 | 4/200 |
|
|
[考研] 生物學(xué)學(xué)碩,一志愿湖南大學(xué),初試成績338 +8 | YYYYYNNNNN 2026-03-26 | 10/500 |
|
|
[碩博家園] 考研調(diào)劑 +5 | 駱駝男人 2026-04-01 | 5/250 |
|
|
[考研] 08工科,295,接受跨專業(yè)調(diào)劑 +6 | lmnlzy 2026-03-31 | 6/300 |
|
|
[考研] 考研材料工程351分調(diào)劑 +5 | 整個(gè)好的 2026-03-31 | 5/250 |
|
|
[考研] 254材料與化工求調(diào)劑 +3 | 翰冬林楠 2026-03-30 | 4/200 |
|
|
[考研] 江蘇蘇北高校誠邀調(diào)劑同學(xué) +3 | zzll406 2026-03-31 | 3/150 |
|
|
[考研] 0856 335分 +9 | cccchenso 2026-03-29 | 9/450 |
|
|
[考研] 353求調(diào)劑 +3 | 江上楓_26 2026-03-28 | 3/150 |
|
|
[考研] 本科211生物醫(yī)學(xué)工程085409求調(diào)劑339分 +7 | 里子木yy 2026-03-29 | 7/350 |
|
|
[考研] 考研調(diào)劑 +5 | Sanmu-124 2026-03-26 | 5/250 |
|
|
[考研] 085701求調(diào)劑初試286分 +5 | secret0328 2026-03-28 | 5/250 |
|
|
[考研] 283求調(diào)劑(080500) +14 | A child 2026-03-27 | 14/700 |
|
|
[考研] 求調(diào)劑,一志愿 南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩 +6 | @taotao 2026-03-26 | 7/350 |
|
|
[考研] 一志愿南京航空航天大學(xué)材料學(xué)碩求調(diào)劑 +3 | @taotao 2026-03-28 | 3/150 |
|