24小时热门版块排行榜    

查看: 1456  |  回复: 14
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

didadida321

金虫 (正式写手)

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

在对encut测试时,为什么得到的comment中没有对应的能量?刚学习vasp,望各位帮帮忙,谢了。
回复此楼
已阅   回复此楼   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] (出分啦!!!)26 考研调剂 20w+调剂信息火速更新中!【捡漏】qq 群:1087823172 +3 db58 2026-02-28 79/3950 2026-02-28 17:50 by nbf68
[考研] 085600材料工程一志愿中科大总分312求调剂 +4 吃宵夜1 2026-02-28 6/300 2026-02-28 17:38 by 吃宵夜1
[考研] 311求调剂 +6 南迦720 2026-02-28 6/300 2026-02-28 17:17 by 刘茉莉
[考研] 高分子化学与物理调剂 +3 好好好1233 2026-02-28 4/200 2026-02-28 17:13 by sandychj
[考研] 化工专硕348,一志愿985求调剂 +3 弗格个 2026-02-28 5/250 2026-02-28 17:04 by sandychj
[考博] 博士推荐 +4 花儿笑? 2026-02-21 5/250 2026-02-28 16:55 by mumin1990
[考研] 0856材料求调剂 +8 hyf hyf hyf 2026-02-28 9/450 2026-02-28 16:27 by etapple
[考研] 285求调剂 +3 满头大汗的学生 2026-02-28 3/150 2026-02-28 16:22 by 无际的草原
[考研] 265分求调剂不调专业和学校有行学上就 +4 礼堂丁真258 2026-02-28 6/300 2026-02-28 16:18 by 求调剂zz
[考博] 26申博 +3 想申博! 2026-02-26 3/150 2026-02-28 16:07 by nxgogo
[考研] 材料类求调剂 +3 wana_kiko 2026-02-28 3/150 2026-02-28 15:03 by lature00
[考研] 295求调剂 +4 19171856320 2026-02-28 4/200 2026-02-28 13:39 by ms629
[考研] 0856调剂 +3 刘梦微 2026-02-28 3/150 2026-02-28 13:22 by houyaoxu
[考研] 304求调剂 +5 曼殊2266 2026-02-28 6/300 2026-02-28 12:44 by 迷糊CCPs
[硕博家园] 博士自荐 +6 科研狗111 2026-02-26 9/450 2026-02-28 12:32 by seaskyy
[考研] 272求调剂 +3 田智友 2026-02-28 3/150 2026-02-28 12:31 by 王加浩to
[考研] 298求调剂 +4 axyz3 2026-02-28 4/200 2026-02-28 11:21 by wang_dand
[基金申请] 面上可以超过30页吧? +12 阿拉贡aragon 2026-02-22 13/650 2026-02-26 22:09 by Hahaxia
[教师之家] 版面费该交吗 +14 苹果在哪里 2026-02-22 17/850 2026-02-26 11:55 by AGI智创机器人
[硕博家园] 【博士招生】太原理工大学2026化工博士 +4 N1ce_try 2026-02-24 8/400 2026-02-26 08:40 by N1ce_try
信息提示
请填处理意见