24小时热门版块排行榜    

查看: 1466  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 考研调剂 +3 呼呼?~+123456 2026-03-21 3/150 2026-03-21 20:04 by 无际的草原
[考研] 326求调剂 +4 mlpqaz03 2026-03-15 4/200 2026-03-21 19:10 by ColorlessPI
[考研] 311求调剂 +7 26研0 2026-03-15 7/350 2026-03-21 19:06 by ColorlessPI
[考研] 0805 316求调剂 +3 大雪深藏 2026-03-18 3/150 2026-03-21 18:55 by 学员8dgXkO
[考研] 297求调剂 +3 喜欢还是不甘心 2026-03-20 3/150 2026-03-21 18:33 by 学员8dgXkO
[考研] 求调剂 +3 13341 2026-03-20 3/150 2026-03-21 18:28 by 学员8dgXkO
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +6 嘉年新程 2026-03-15 6/300 2026-03-21 17:07 by Dream007008
[考研] 求调剂 +3 白QF 2026-03-21 3/150 2026-03-21 13:12 by zhukairuo
[考研] 306求调剂 +4 chuanzhu川烛 2026-03-18 4/200 2026-03-21 08:25 by laoshidan
[考研] 346求调剂[0856] +4 WayneLim327 2026-03-16 7/350 2026-03-21 04:02 by JourneyLucky
[考研] 307求调剂 +10 冷笙123 2026-03-17 10/500 2026-03-21 01:54 by JourneyLucky
[考研] 材料专硕英一数二306 +7 z1z2z3879 2026-03-18 7/350 2026-03-20 23:48 by JourneyLucky
[考研] 308求调剂 +3 阿姐阿姐家啊 2026-03-18 3/150 2026-03-20 23:24 by JourneyLucky
[考研] 考研调剂求学校推荐 +3 伯乐29 2026-03-18 5/250 2026-03-20 22:59 by JourneyLucky
[考研] 329求调剂 +9 想上学吖吖 2026-03-19 9/450 2026-03-20 22:01 by luoyongfeng
[考研] 260求调剂 +3 朱芷琳 2026-03-20 3/150 2026-03-20 20:35 by 学员8dgXkO
[考研] 一志愿吉林大学材料学硕321求调剂 +11 Ymlll 2026-03-18 15/750 2026-03-20 19:40 by 丁丁*
[考研] 求调剂 +3 @taotao 2026-03-20 3/150 2026-03-20 19:35 by JourneyLucky
[考研] 求调剂 +3 暗涌afhb 2026-03-16 3/150 2026-03-20 00:28 by 河南大学校友
[考研] 344求调剂 +6 knight344 2026-03-16 7/350 2026-03-18 20:13 by walc
信息提示
请填处理意见