24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1062  |  回复: 2

GTF1902

金虫 (小有名气)

[求助] 【重新开帖】KPOINTS的测试脚本问题

我用vasp优化结构 先测试K点时出了点问题,请大家看完 下面是我的一个优化结构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
刚测试完第五个点
问题是为啥同样的数据要出现6次啊 是不是我的测试脚本文件出问题了 这是我第二遍算了 第一遍测试4个K点也是同样的数据在comment文件中显示六次 是不是电脑重复计算了六次啊,那就亏大了 我是刚学习 也不知哪里出问题了

[ Last edited by ben_ladeng on 2013-5-15 at 18:21 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   关注TA 给TA发消息 送TA红花 TA的回帖

171713294

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★
感谢参与,应助指数 +1
GTF1902: 金币+2, 有帮助 2013-05-13 09:16:41
GTF1902: 回帖置顶 2013-05-13 10:01:05
franch: 金币+2, 谢谢回帖交流,, 2013-05-13 21:19:50
楼住,给你个建议,你吧你的问题一列列弄好,别这样一大片,没人会给你看的,重新弄好 我们给你看看
坚持,不放弃
2楼2013-05-13 08:45:57
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

GTF1902

金虫 (小有名气)

GTF1902: 回帖置顶 2013-05-13 10:01:55
引用回帖:
2楼: Originally posted by 171713294 at 2013-05-13 08:45:57
楼住,给你个建议,你吧你的问题一列列弄好,别这样一大片,没人会给你看的,重新弄好 我们给你看看

我用vasp优化结构前 先进行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楼2013-05-13 10:01:00
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 GTF1902 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求计算机方向调剂 +3 Toffee2 2026-04-16 6/300 2026-04-19 22:37 by ll叶
[考研] 通信工程求调剂!!! +7 zlb770521 2026-04-14 7/350 2026-04-19 20:56 by Equinoxhua
[考研] 26药学专硕105500求调剂 +7 喽哈加油 2026-04-13 8/400 2026-04-19 20:21 by Equinoxhua
[考研] 304求调剂 +8 castLight 2026-04-16 8/400 2026-04-19 17:14 by 中豫男
[考研] 327求调剂 +27 Xxjc1107. 2026-04-13 30/1500 2026-04-19 08:22 by cuisz
[考研] 0854求调剂 +23 门路摸摸 2026-04-15 27/1350 2026-04-19 01:59 by 烟雨流涯
[考研] 320求调剂 +5 深郊akm 2026-04-17 5/250 2026-04-18 19:52 by 王珺璞
[考研] 接受任何调剂 +6 也就是栗子 2026-04-17 7/350 2026-04-18 17:20 by 涵竹刘
[考研] 22408 312求调剂 +24 门路摸摸 2026-04-14 26/1300 2026-04-18 13:04 by wunaiy88
[考研] 收到复试调剂但是去不了 +8 小蜗牛* 2026-04-16 8/400 2026-04-18 11:15 by zixin2025
[考研] 260求调剂 +4 Zyt1314520.. 2026-04-17 5/250 2026-04-18 08:28 by babysonlkd
[考研] 急需调剂 +9 绝不放弃22 2026-04-15 10/500 2026-04-18 08:09 by chixmc
[考研] 297,工科调剂?河南农业大学本科 +14 河南农业大学-能 2026-04-14 14/700 2026-04-16 14:41 by dingyanbo1
[基金申请] RY:中国产出的科学垃圾论文,绝对数量和比例都世界第一 +7 zju2000 2026-04-14 18/900 2026-04-16 11:36 by 欢乐颂叶蓁
[考研] 一志愿A区211,22408 321求调剂 +6 随心所欲☆ 2026-04-15 7/350 2026-04-15 21:45 by lbsjt
[考研] 297工科调剂? +14 河南农业大学-能 2026-04-13 15/750 2026-04-15 13:25 by 黑科技矿业
[考研] 考研调剂 +13 长弓傲 2026-04-13 14/700 2026-04-14 14:44 by zs92450
[考研] 085408光电信息工程专硕355一志愿长春光机所调剂 +6 王ymaa 2026-04-13 13/650 2026-04-14 11:33 by 王ymaa
[考研] 245求调剂 +6 冰糖橘?汽水 2026-04-13 10/500 2026-04-14 10:49 by jyl0317
[考研] 085600材料与化工329分求调剂 +24 叶zilin 2026-04-13 25/1250 2026-04-14 09:20 by 试管破裂
信息提示
请填处理意见