24小时热门版块排行榜    

查看: 1467  |  回复: 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿北京化工大学070300 学硕336求调剂 +3 vv迷 2026-03-21 4/200 2026-03-21 23:05 by f19980501
[考研] 初试 317 +7 半拉月丙 2026-03-20 7/350 2026-03-21 22:26 by peike
[考研] 化学调剂 +5 yzysaa 2026-03-21 5/250 2026-03-21 22:12 by peike
[考研] 333求调剂 +5 87639 2026-03-21 7/350 2026-03-21 19:31 by ColorlessPI
[考研] 0703化学调剂 +11 妮妮ninicgb 2026-03-15 15/750 2026-03-21 19:15 by ColorlessPI
[考研] 307求调剂 +3 余意卿 2026-03-18 3/150 2026-03-21 17:31 by ColorlessPI
[考研] 296求调剂 +4 www_q 2026-03-20 4/200 2026-03-21 17:26 by 学员8dgXkO
[考研] 298求调剂 +4 上岸6666@ 2026-03-20 4/200 2026-03-21 17:14 by 学员8dgXkO
[考研] 279求调剂 +5 红衣隐官 2026-03-21 5/250 2026-03-21 14:59 by lature00
[考研] 279分求调剂 一志愿211 +14 chaojifeixia 2026-03-19 15/750 2026-03-21 13:24 by zhukairuo
[考研] 304求调剂 +6 曼殊2266 2026-03-18 6/300 2026-03-21 00:32 by JourneyLucky
[考研] 288求调剂 +16 于海海海海 2026-03-19 16/800 2026-03-20 22:28 by JourneyLucky
[考研] 材料学硕297已过四六级求调剂推荐 +11 adaie 2026-03-19 11/550 2026-03-20 21:30 by laoshidan
[考研] 一志愿西安交通大学 学硕 354求调剂211或者双一流 +3 我想要读研究生 2026-03-20 3/150 2026-03-20 20:13 by JourneyLucky
[考研] 一志愿吉林大学材料学硕321求调剂 +11 Ymlll 2026-03-18 15/750 2026-03-20 19:40 by 丁丁*
[考研] 344求调剂 +6 knight344 2026-03-16 7/350 2026-03-18 20:13 by walc
[考博] 26博士申请 +3 1042136743 2026-03-17 3/150 2026-03-17 23:30 by 轻松不少随
[考博] 26申博 +4 八6八68 2026-03-16 4/200 2026-03-17 13:00 by 轻松不少随
[考研] 机械专硕325,寻找调剂院校 +3 y9999 2026-03-15 5/250 2026-03-16 19:58 by y9999
[考研] 321求调剂 +5 大米饭! 2026-03-15 5/250 2026-03-16 16:33 by houyaoxu
信息提示
请填处理意见