24小时热门版块排行榜    

查看: 1135  |  回复: 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的回帖
回帖置顶 ( 共有2个 )

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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 什么时候开奖? +4 CrisMessi 2026-08-18 4/200 2026-08-19 00:53 by gatelove
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 27/1350 2026-08-18 21:44 by Noways
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +3 瞬息宇宙 2026-08-15 3/150 2026-08-18 21:40 by xiaohuixin
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 9/450 2026-08-18 20:09 by 奶牛小黑
[基金申请] 朋友圈看到的 +4 wangzilk 2026-08-18 5/250 2026-08-18 19:25 by wangzilk
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +51 医学老男孩 2026-08-13 110/5500 2026-08-18 19:07 by icerianxia
[基金申请] 时间戳变了,能看出什么问题? +8 基诺咪客 2026-08-17 11/550 2026-08-18 17:13 by 悠然熊猫
[基金申请] 重要消息,中午系统在维护 +9 yuleib84 2026-08-18 10/500 2026-08-18 16:52 by archvillain
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +4 Tide man 2026-08-14 5/250 2026-08-18 11:35 by Tide man
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
信息提示
请填处理意见