应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » encut测试时，为什么得到的comment中没有对应的能量？
51/1返回列表
查看: 1350  |  回复: 14
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

didadida321

金虫 (正式写手)

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

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

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2012-05-04 08:22:35
已阅   回复此楼   关注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的回帖
查看全部 15 个回答

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的回帖
查看全部 15 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭