| 3 | 1/1 | 返回列表 |
| 查看: 1043 | 回復: 2 | ||
GTF1902金蟲 (小有名氣)
|
[求助]
【重新開帖】KPOINTS的測試腳本問題
|
|
|
我用vasp優(yōu)化結構 先測試K點時出了點問題,請大家看完 下面是我的一個優(yō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)化結構前 先進行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ā)表 | |
|---|---|---|---|---|
|
[考研] 306求調(diào)劑 +7 | chuanzhu川燭 2026-03-18 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 289材料與化工(085600)B區(qū)求調(diào)劑 +3 | 這么名字咋樣 2026-03-22 | 4/200 |
|
|
[考研] 287求調(diào)劑 +8 | 晨昏線與星海 2026-03-19 | 9/450 |
|
|
[考研] 275求調(diào)劑 +6 | shansx 2026-03-22 | 8/400 |
|
|
[考研] 303求調(diào)劑 +5 | 安憶靈 2026-03-22 | 6/300 |
|
|
[考研] 生物學一志愿985,分數(shù)349求調(diào)劑 +4 | zxts12 2026-03-21 | 7/350 |
|
|
[考研] 初試 317 +7 | 半拉月丙 2026-03-20 | 7/350 |
|
|
[考研] 材料求調(diào)劑 +5 | @taotao 2026-03-21 | 5/250 |
|
|
[考研] 278求調(diào)劑 +9 | 煙火先于春 2026-03-17 | 9/450 |
|
|
[考研] 313求調(diào)劑 +4 | 肆叁貳壹22 2026-03-19 | 4/200 |
|
|
[考研] 一志愿山大07化學 332分 四六級已過 本科山東雙非 求調(diào)劑! +3 | 不想理你 2026-03-16 | 3/150 |
|
|
[考研] 301求調(diào)劑 +10 | yy要上岸呀 2026-03-17 | 10/500 |
|
|
[考研] 二本跨考鄭大材料306英一數(shù)二 +3 | z1z2z3879 2026-03-17 | 3/150 |
|
|
[考研] 265求調(diào)劑 +9 | 梁梁校校 2026-03-17 | 9/450 |
|
|
[考研] 材料專碩英一數(shù)二306 +7 | z1z2z3879 2026-03-18 | 7/350 |
|
|
[考研] 南京大學化學376求調(diào)劑 +3 | hisfailed 2026-03-19 | 6/300 |
|
|
[考研] 一志愿南京理工大學085701資源與環(huán)境302分求調(diào)劑 +4 | 葵梓衛(wèi)隊 2026-03-18 | 6/300 |
|
|
[考研] 中南大學化學學碩337求調(diào)劑 +3 | niko- 2026-03-19 | 6/300 |
|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|
[考研] 283求調(diào)劑 +3 | 聽風就是雨; 2026-03-16 | 3/150 |
|