24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1349  |  回复: 14

didadida321

金虫 (正式写手)

[求助] encut测试时,为什么得到的comment中没有对应的能量?

在对encut测试时,为什么得到的comment中没有对应的能量?刚学习vasp,望各位帮帮忙,谢了。
回复此楼

» 猜你喜欢

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

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

liliangfang

荣誉版主 (著名写手)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
didadida321: 金币+2, ★★★很有帮助 2012-05-04 09:48:48
这个可能是你的脚本有点问题,请检查一下
2楼2012-05-04 09:00:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liliangfang

荣誉版主 (著名写手)

【答案】应助回帖


sunyang1988: 金币+1, 谢谢帮助 2012-05-04 12:08:10
#!/bin/sh
rm WAVECAR
for i in 250 300 350 400 450 500 550 600
do
cat > INCAR < System = Si
ENCUT = $i
ISTART = 0
ICHARG = 2
ISMEAR = 0
SIGMA = 0.1
PREC = Accurate
!
echo "ENCUT = $i eV" ;  mpdboot ; mpirun -np 4 vasp
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
这个你最好不是用复制的,自己打进去,这个mpdboot ; mpirun -np 4 vasp
根据你需要修改

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : encut
  • 2012-05-04 09:05:15, 391 bytes
3楼2012-05-04 09:05:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

didadida321

金虫 (正式写手)

引用回帖:
3楼: Originally posted by liliangfang at 2012-05-04 09:05:16:
#!/bin/sh
rm WAVECAR
for i in 250 300 350 400 450 500 550 600
do
cat > INCAR <<!
System = Si
ENCUT = $i
ISTART = 0
ICHARG = 2
ISMEAR = 0
SIGMA = 0.1
PREC = Accurate
!
echo  ...

谢谢了
4楼2012-05-04 09:48:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

techfan

铁虫 (小有名气)

请问,我按照这个脚本修改了一下,想得到SIGMA数值,可结果却如下,
0.10 0.000000
0.12 0.000000
0.14 0.000000
0.16 0.000000
0.18 0.000000
0.20 0.000000
0.22 0.000000
0.24 0.000000
0.26 0.000000
0.28 0.000000
0.30 0.000000
后来查看了一下OUTCAR,发现
entropy T*S    EENTRO =         0.00000000
而且,循环多次,EENTRO一直都算0,这是什么情况,如何修改?谢谢
5楼2012-05-04 10:15:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liliangfang

荣誉版主 (著名写手)

【答案】应助回帖

★ ★
sunyang1988: 金币+2, 谢谢帮助 2012-05-04 12:08:00
引用回帖:
5楼: Originally posted by techfan at 2012-05-04 10:15:02:
请问,我按照这个脚本修改了一下,想得到SIGMA数值,可结果却如下,
0.10 0.000000
0.12 0.000000
0.14 0.000000
0.16 0.000000
0.18 0.000000
0.20 0.000000
0.22 0.000000
0.24 0.000000
0.26 0.00000 ...

这是SIGMA,KPOINTS,晶格常数的脚本,你自己修改一下就可以了

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : k
  • 2012-05-04 10:37:09, 343 bytes
  • 附件 2 : rvasp
  • 2012-05-04 10:37:09, 450 bytes
  • 附件 3 : sigma
  • 2012-05-04 10:37:09, 419 bytes
6楼2012-05-04 10:37:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

techfan

铁虫 (小有名气)

引用回帖:
6楼: Originally posted by liliangfang at 2012-05-04 10:37:11:
这是SIGMA,KPOINTS,晶格常数的脚本,你自己修改一下就可以了

谢谢了,我用你提供的脚本试试看
7楼2012-05-04 12:42:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

techfan

铁虫 (小有名气)

实在不好意思打扰了,我用了你提供的脚本,结果依然没有变化,
entropy T*S    EENTRO =         0.00000000
而且,循环多次,EENTRO一直都算0,
到底怎么改呢?555
8楼2012-05-04 13:49:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liliangfang

荣誉版主 (著名写手)

【答案】应助回帖


sunyang1988: 金币+1, 谢谢交流 2012-05-04 18:31:34
引用回帖:
8楼: Originally posted by techfan at 2012-05-04 13:49:06:
实在不好意思打扰了,我用了你提供的脚本,结果依然没有变化,
entropy T*S    EENTRO =         0.00000000
而且,循环多次,EENTRO一直都算0,
到底怎么改呢?555

#!/bin/sh
rm WAVECAR
for i in 0.01 0.02 0.04 0.05 0.06 0.08 0.10 0.12 0.14 0.16 0.18 0.20 do
cat > INCAR < System = Si
ENCUT = 350
ISTART = 0
ICHARG = 2
ISMEAR = 0
SIGMA = $i
PREC = Accurate
!
echo "SIGMA = $i eV" ;  vasp
TS=`grep "EENTRO" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5 }'`
KP=`grep "irreducible" OUTCAR | tail -1 | awk '{printf "%5i \n", $2 }'`
echo $i $TS >>comment
done

红色的地方需要你自己修改,如果不行那就不是脚本的问题啦
9楼2012-05-04 15:42:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

techfan

铁虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
是的,我已经安装自己的情况修改了,能顺利运行,问题就是
所有循环得到的EENTRO全部为0,导致
0.10 0.000000
0.12 0.000000
0.14 0.000000
0.16 0.000000
0.18 0.000000
0.20 0.000000
0.22 0.000000
0.24 0.000000
0.26 0.000000
0.28 0.000000
0.30 0.000000
不知道为啥?
10楼2012-05-04 15:51:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 didadida321 的主题更新
信息提示
请填处理意见