| 3 | 1/1 | 返回列表 |
| 查看: 1045 | 回復: 2 | ||
GTF1902金蟲 (小有名氣)
|
[求助]
【重新開帖】KPOINTS的測試腳本問題
|
|
我用vasp優(yōu)化結(jié)構(gòu) 先測試K點時出了點問題,請大家看完 下面是我的一個優(yōu)化結(jié)構(gòu)vasp文件 INCAR : System = PCCOF ISIF = 2 ISMEAR = 0 ISTART = 0 ICHARG = 2 SIGMA = 0.1 NSW = 0 PREC = Medium ENCUT = 500 GGA = 91 VOSKOWN = 1 IBRION = 2 EDIFF = 1E-4 EDIFFG = -0.03 LREAL = Auto IALGO = 48 LALANE = .TRUE NPAR = 6 LSCALU = .FALSE. NSIM = 4 RUN_K: #!/bin/sh rm -f WAVECAR for i in 1 2 3 4 5 6 7 do cat>KPOINTS<<! Automatic generation 0 Monkhorst-Pack $i $i $i 0.0 0.0 0.0 ! echo"k mesh = $i x $i x $i"; time /share/home/sci_028/sourcecode/vasp.5.2/vasp>out.log 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 還沒算完 但是 comment文件非常怪異:大俠請看 1 1 -569.842023 1 1 -569.842023 1 1 -569.842023 1 1 -569.842023 1 1 -569.842023 1 1 -569.842023 2 1 -600.814240 2 1 -600.814240 2 1 -600.814240 2 1 -600.814240 2 1 -600.814240 2 1 -600.814240 3 8 -599.488993 3 8 -599.488993 3 8 -599.488993 3 8 -599.488993 3 8 -599.488993 3 8 -599.488993 4 8 -599.575041 4 8 -599.575041 4 8 -599.575041 4 8 -599.575041 4 8 -599.575041 4 8 -599.575041 5 27 -599.519665 5 27 -599.519665 5 27 -599.519665 5 27 -599.519665 5 27 -599.519665 5 27 -599.519665 剛測試完第五個點 問題是為啥同樣的數(shù)據(jù)要出現(xiàn)6次啊 是不是我的測試腳本文件出問題了 這是我第二遍算了 第一遍測試4個K點也是同樣的數(shù)據(jù)在comment文件中顯示六次 是不是電腦重復計算了六次啊,那就虧大了 我是剛學習 也不知哪里出問題了 [ Last edited by ben_ladeng on 2013-5-15 at 18:21 ] |
木蟲 (正式寫手)

金蟲 (小有名氣)
|
我用vasp優(yōu)化結(jié)構(gòu)前 先進行K測試,下面是我的incar 和K_RUN文件 System = PCCOF ISIF = 2 ISMEAR = 0 ISTART = 0 ICHARG = 2 SIGMA = 0.1 NSW = 0 PREC = Medium ENCUT = 500 GGA = 91 VOSKOWN = 1 IBRION = 2 EDIFF = 1E-4 EDIFFG = -0.03 LREAL = Auto IALGO = 48 LALANE = .TRUE NPAR = 6 LSCALU = .FALSE. NSIM = 4 我是用先用下面的txt提交任務: APP_NAME=scivip OMP_NUM_THREADS=1 NP_PER_NODE=6 MY_MPI_TYPE=openmpi MY_MPI_HOME=/share/home/sci_028/sourcecode/openmpi-1.4.2 NP=6 RUN="/share/home/sci_028/gtf/vasp/shiyansan/run_k.txt >out.log" 下面是的 run_k.txt: #!/bin/sh rm -f WAVECAR for i in 1 2 3 4 5 6 7 do cat > KPOINTS <<! Automatic generation 0 Monkhorst-Pack $i $i $i 0.0 0.0 0.0 ! echo "k mesh = $i x $i x $i"; time /share/home/sci_028/sourcecode/vasp.5.2/vasp >out.log 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 但是comment輸出卻為: 1 1 -569.842023 1 1 -569.842023 1 1 -569.842023 1 1 -569.842023 1 1 -569.842023 1 1 -569.842023 2 1 -600.814240 2 1 -600.814240 2 1 -600.814240 2 1 -600.814240 2 1 -600.814240 2 1 -600.814240 3 8 -599.488993 3 8 -599.488993 3 8 -599.488993 3 8 -599.488993 3 8 -599.488993 3 8 -599.488993 4 8 -599.575041 4 8 -599.575041 4 8 -599.575041 4 8 -599.575041 4 8 -599.575041 4 8 -599.575041 5 27 -599.519665 5 27 -599.519665 5 27 -599.519665 5 27 -599.519665 5 27 -599.519665 5 27 -599.519665 6 27 -599.540837 6 27 -599.540837 6 27 -599.540837 6 27 -599.540837 6 27 -599.540837 6 27 -599.540837 7 64 -599.557505 7 64 -599.557505 7 64 -599.557505 7 64 -599.557505 7 64 -599.557505 7 64 -599.557505 也就說沒個k點輸出了六次 是不是計算了六次啊 真的好郁悶 那就虧大了 問題出在哪里? |
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿武理材料工程348求調(diào)劑 +6 |  ̄^ ̄゜汗 2026-03-19 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 接收2026碩士調(diào)劑(學碩+專碩) +4 | allen-yin 2026-03-23 | 6/300 |
|
|
[考研] 284求調(diào)劑 +6 | Zhao anqi 2026-03-22 | 6/300 |
|
|
[考研] 307求調(diào)劑 +11 | 冷笙123 2026-03-17 | 11/550 |
|
|
[考研] 298求調(diào)劑一志愿211 +3 | 上岸6666@ 2026-03-20 | 3/150 |
|
|
[考研] 303求調(diào)劑 +5 | 安憶靈 2026-03-22 | 6/300 |
|
|
[考研] 求調(diào)劑 +7 | Auroracx 2026-03-22 | 7/350 |
|
|
[考研] 考研調(diào)劑 +4 | 來好運來來來 2026-03-21 | 4/200 |
|
|
[考研] 297求調(diào)劑 +3 | 喜歡還是不甘心 2026-03-20 | 3/150 |
|
|
[考研] 一志愿深大,0703化學,總分302,求調(diào)劑 +4 | 七月-七七 2026-03-21 | 4/200 |
|
|
[考研] 311求調(diào)劑 +3 | 勇敢的小吳 2026-03-20 | 3/150 |
|
|
[考研] 材料 271求調(diào)劑 +5 | 展信悅_ 2026-03-21 | 5/250 |
|
|
[考研] 265求調(diào)劑 +12 | 梁梁校校 2026-03-19 | 14/700 |
|
|
[考研] 346求調(diào)劑[0856] +4 | WayneLim327 2026-03-16 | 7/350 |
|
|
[考研] 一志愿華中科技大學,080502,354分求調(diào)劑 +5 | 守候夕陽CF 2026-03-18 | 5/250 |
|
|
[考研] 一志愿武理材料305分求調(diào)劑 +6 | 想上岸的鯉魚 2026-03-18 | 7/350 |
|
|
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3 | 葵梓衛(wèi)隊 2026-03-20 | 3/150 |
|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|
[考研] 一志愿福大288有機化學,求調(diào)劑 +3 | 小木蟲200408204 2026-03-18 | 3/150 |
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|