| 5 | 1/1 | 返回列表 |
| 查看: 2604 | 回復(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ā)表 | |
|---|---|---|---|---|
|
[考研] 材料考研調(diào)劑 +4 | Gs大王 2026-04-02 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 366求調(diào)劑一志愿東北大學(xué) +3 | 運(yùn)氣來得若有似?/a> 2026-04-02 | 3/150 |
|
|
[考研] 085801 總分275 本科新能源 求調(diào)劑 +10 | bradoner 2026-04-01 | 11/550 |
|
|
[考研] 284求調(diào)劑 +12 | 小熊~~ 2026-03-31 | 12/600 |
|
|
[考研] 286求調(diào)劑 +5 | lim0922 2026-03-26 | 5/250 |
|
|
[考研] 省雙一流重點(diǎn)一本大學(xué)招收調(diào)劑 +4 | wwwwffffff 2026-03-31 | 7/350 |
|
|
[考研] 材料專碩322分 +9 | 哈哈哈吼吼吼哈 2026-04-01 | 9/450 |
|
|
[考研] 一志愿同濟(jì)大學(xué)323分(080500)求調(diào)劑 +4 | yikeniu 2026-04-01 | 4/200 |
|
|
[考研] 262求調(diào)劑 +9 | 勵(lì)志一定發(fā)文章 2026-03-31 | 10/500 |
|
|
[考研] 材料調(diào)劑 +11 | 一樣YWY 2026-03-31 | 11/550 |
|
|
[考研] 考研調(diào)劑 +11 | Amber00 2026-03-31 | 11/550 |
|
|
[考研] 085900土木水利336分求調(diào)劑 +3 | Zhangjiangj 2026-03-31 | 5/250 |
|
|
[考研] 材料工程專碩求調(diào)劑 +10 | hyl3153942 2026-03-29 | 10/500 |
|
|
[考研] 276求調(diào)劑 +3 | 趙久華 2026-03-29 | 3/150 |
|
|
[考研] 293分求調(diào)劑,外語為俄語 +5 | 加一一九 2026-03-31 | 5/250 |
|
|
[考研] 哈爾濱工業(yè)大學(xué)材料與化工專碩378求調(diào)劑 +3 | 塔比烏斯 2026-03-30 | 3/150 |
|
|
[考研] 085701求調(diào)劑初試286分 +5 | secret0328 2026-03-28 | 5/250 |
|
|
[考研] 086000生物與醫(yī)藥調(diào)劑 +5 | Feisty。 2026-03-28 | 9/450 |
|
|
[考研] 求調(diào)劑 +7 | 爭(zhēng)取九點(diǎn)睡 2026-03-28 | 8/400 |
|
|
[考研] 285求調(diào)劑 +4 | AZMK 2026-03-27 | 7/350 |
|