24小时热门版块排行榜    

查看: 1356  |  回复: 6

默默七月

新虫 (小有名气)

[求助] K 点收敛性测试comment文件无能量输出已有1人参与

如题,进行K 点收敛性测试comment文件无能量输出,按的是侯老师教程写的脚本,提交命令也改成自己的。下面是comment文件
7 0.000000
9 0.000000
11 0.000000
13 0.000000
15 0.000000
17 0.000000
19 0.000000
脚本文件
cd $PBS_O_WORKDIR
NP=`cat $PBS_NODEFILE|wc -l`
echo $PBS_O_WORKDIR
#!/bin/sh
rm WAVECAR
for i in  7 9 11 13 15 17 19
do
cat >KPOINTS <<!
Automatic generation
0
Monhkorst-Pack
$i  $i  1
0.0 0.0 0.0
!
echo "k mesh= $i x $i x " ; mpirun -np $NP -machinefile $PBS_NODEFILE /public2/home/wuli/liupei5134/bin/vasp53
KP=`grep "irreducible" OUTCAR|tail -1 |awk '{printf"%12.6f\n",$5}'`
echo $i $KP $E >>comment
done
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

g_wcai

木虫 (正式写手)

开心

【答案】应助回帖


感谢参与,应助指数 +1
liliangfang: 金币+1, 谢谢交流 2015-12-01 21:12:58
这个问题可能遇到过,相当于提交了几个不同k点设置的计算,你的comment出来时,你去查看一下你提交的任务完成了吗?任务没完成出来的ENErGe当然为零!!!

发自小木虫Android客户端
哈哈,笑口常开!
2楼2015-11-24 15:59:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

默默七月

新虫 (小有名气)

引用回帖:
2楼: Originally posted by g_wcai at 2015-11-24 15:59:07
这个问题可能遇到过,相当于提交了几个不同k点设置的计算,你的comment出来时,你去查看一下你提交的任务完成了吗?任务没完成出来的ENErGe当然为零!!!

谢谢,但是我的任务计算完成了,请问还有其他原因吗
3楼2015-11-24 16:51:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

g_wcai

木虫 (正式写手)

开心

【答案】应助回帖


liliangfang: 金币+1, 谢谢交流 2015-12-01 21:13:04
那你查看一下Outcar中能量是多少,如果是o,就有可能你的程序有问题!

发自小木虫Android客户端
哈哈,笑口常开!
4楼2015-11-24 17:17:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

默默七月

新虫 (小有名气)

引用回帖:
4楼: Originally posted by g_wcai at 2015-11-24 17:17:32
那你查看一下Outcar中能量是多少,如果是o,就有可能你的程序有问题!

好的,谢谢
5楼2015-11-24 18:22:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

g_wcai

木虫 (正式写手)

开心

【答案】应助回帖


liliangfang: 金币+1, 谢谢交流 2015-12-01 21:13:12
还有,你程序中KP=`grep "irreducible" OUTCAR|tail -1 |awk '{printf"%12.6f\n",$5}'`这个语句是不是取能量值?!
如果是可能有问题,irreducible是什么能量?一般用grep "TOTEN" OUTCAR!
哈哈,笑口常开!
6楼2015-11-24 20:56:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

andyhox

铜虫 (初入文坛)

mpirun -np $NP -machinefile $PBS_NODEFILE /public2/home/wuli/liupei5134/bin/vasp53
请问这一行是什么意思   新手小白
7楼2019-03-14 20:03:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 默默七月 的主题更新
信息提示
请填处理意见