24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 1491  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 305求调剂 +8 玛卡巴卡boom 2026-04-11 8/400 2026-04-14 09:04 by pengliang8036
[考研] 机械工程313分找工科调剂 +4 双一流本科机械 2026-04-08 4/200 2026-04-14 07:32 by Abskk
[考研] 一志愿双非085400电子信息344 求调剂,对材料和化学方向也感兴趣 +12 无情的小羊 2026-04-09 13/650 2026-04-13 14:17 by 张zhihao
[硕博家园] 新一代电子信息294求调剂 不挑学校 +7 Ytyt11 2026-04-09 8/400 2026-04-12 16:57 by ajpv风雷
[找工作] 山东高校教师考核超级无底线,员工过不下去啦 +4 qut2026 2026-04-09 9/450 2026-04-12 00:54 by qut2026
[考研] 277 数一104,学硕,求调剂 +21 瓶子PZ 2026-04-09 23/1150 2026-04-11 23:12 by labixiaoqiao
[考研] 材料工程日语考生求调剂 +7 0856?调剂 2026-04-10 7/350 2026-04-11 21:33 by 蓝云思雨
[考研] 303求调剂 +14 SereinQ 2026-04-10 15/750 2026-04-11 20:43 by 蓝云思雨
[考研] 086003调剂求助 +21 苏弋万 2026-04-09 22/1100 2026-04-11 20:25 by dongdian1
[考研] 284求调剂 +11 archer.. 2026-04-09 12/600 2026-04-11 20:23 by 蓝云思雨
[考研] 085400 328分 求调剂 +10 喂你一个大橙子 2026-04-09 14/700 2026-04-11 19:53 by lqspecial
[考研] 283求调剂 +22 那个噜子 2026-04-09 22/1100 2026-04-11 10:41 by 逆水乘风
[考研] 22408 327分求调剂 +4 韵风kon 2026-04-10 4/200 2026-04-11 09:51 by 猪会飞
[考研] 311求调剂 +13 xyp想读书 2026-04-10 14/700 2026-04-11 09:41 by 猪会飞
[考研] 化学工程与技术324调剂 +23 孙常华 2026-04-09 25/1250 2026-04-11 00:07 by 骑牛渡寒江
[考研] 一志愿矿大,材料工程专硕314分,0856可调都可以 +15 无懈可击的巨人 2026-04-09 15/750 2026-04-10 18:10 by hmn_wj
[考研] 344求调剂 +7 丶风雪夜归人丶 2026-04-09 7/350 2026-04-10 12:05 by pengliang8036
[考研] 化学工程与技术专业一志愿哈工程 291分B区 国家级大创负责人 有一作论文 +13 Emmy~ 2026-04-09 13/650 2026-04-09 14:47 by only周
[考研] 334求调剂 +16 Riot2025 2026-04-08 17/850 2026-04-09 09:28 by wdyheheeh
[考研] 085404,334分,求调剂 +5 sunjie8888 2026-04-08 8/400 2026-04-09 07:26 by sunjie8888
信息提示
请填处理意见