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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » encut测试是否需要修改要提交的作业脚本
51/1返回列表
查看: 1333  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
本帖产生 2 个 1ST强帖 ，点击这里进行查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

dingfengbo

新虫 (小有名气)

[求助] encut测试是否需要修改要提交的作业脚本

我想做encut测试这是我的encut脚本：
#!/bin/sh
rm WAVECAR
for a in 400 450 500 550 600 650 700 750
do
cat >INCAR < SYSTEM=CGa
LPLANE=.TRUE.
NPAR=8
LREAL = .TRUE.
PREC    = Accurate
EDIFF   =1E-6
EDIFFG  = -0.03      
IALGO   = 48         
NELMIN  = 4         
ISYM    = 0         
ISPIN = 1
NBANDS = 240
LCHARG  = .TRUE.
LWAVE   = .TRUE.     
LVTOT   = .FALSE.
NBLOCK  = 1         
NSW     = 10
IBRION  = 2
POTIM=0.3
ISMEAR = 0   
SIGMA = 0.05
ENCUT=$a
ISTART=0;ICHARG=2
!
echo "ENCUT=$a eV"
vasp
E=`grep "TOTEN" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5 }'`
echo $a $E >>comment
done
测试的时候我想将之提交，下面是我的作业脚本：
#!/bin/sh
#PBS -l nodes=2:ppn=8
#PBS -u XX
#PBS -N YY

JOBINFOR=$PBS_JOBID
MASTERNODE=`hostname`
SCRATCHDIR=$PBS_JOBID
NCPU=`wc -l < $PBS_NODEFILE`

SERVER=$PBS_O_HOST
WORKDIR=$PBS_O_WORKDIR
MKDIR=/bin/mkdir
RSH=/usr/bin/rsh
CP=/bin/cp
cat $PBS_NODEFILE > $WORKDIR/mpihosts
LAUNCH="mpirun --hostfile  $WORKDIR/mpihosts -np $NCPU"
PROGRAMEXEC="/opt/vasp/vasp"

${RSH} ${MASTERNODE}
cd ${WORKDIR}
${MKDIR} ${SCRATCHDIR}
cd ${SCRATCHDIR}
${CP} ${WORKDIR}/* .
${LAUNCH}  ${PROGRAMEXEC} > out 2>&1
请问我怎样修改作业脚本才能进行encut测试的计算呢？谢谢
回复此楼

» 收录本帖的淘帖专辑推荐

Photochemistry

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-04-22 16:37:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

haowenping

木虫 (著名写手)
引用回帖:
Originally posted by dingfengbo at 2011-04-22 16:37:10:
我想做encut测试这是我的encut脚本：
#!/bin/sh
rm WAVECAR
for a in 400 450 500 550 600 650 700 750
do
cat >INCAR <<!
SYSTEM=CGa
LPLANE=.TRUE.
NPAR=8
LREAL = .TRUE.
PREC    = Accura ...

你好，请问你有没有实现在pbs下优化你的那个截断能脚本？我现在也被这个问题困惑，如果你会了的话麻烦指教下我，谢谢了！
赞一下 回复此楼
人有多大胆，地有多大产！
3楼2011-05-07 19:07:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

bingmou

金虫 (著名写手)

【答案】应助回帖

★
youzhizhe(金币+1, 1ST强帖+1): 谢谢指导。 2011-04-23 08:13:11
如果知识测试小体系的话，建议不要提交，直接单机运行。如果是大体系的话，估计一个ecut就可以了，如果全部是金属，400eV足够；如果有O，N可能要700eV。其他的都可以采用赝势中默认值。
如果一定要提交，可以考虑把这句
PROGRAMEXEC="/opt/vasp/vasp"
改为
PROGRAMEXEC='脚本位置”
赞一下(1人) 回复此楼
2楼2011-04-22 22:58:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingfengbo

新虫 (小有名气)
引用回帖:
Originally posted by haowenping at 2011-05-07 19:07:03:
你好，请问你有没有实现在pbs下优化你的那个截断能脚本？我现在也被这个问题困惑，如果你会了的话麻烦指教下我，谢谢了！

你好！按照2楼说的就可以实现了，不知道你想问什么？
赞一下 回复此楼
4楼2011-05-08 12:52:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiao72379

木虫 (正式写手)

宅男甲

【答案】应助回帖

★
bingmou(金币+1, 1ST强帖+1): 感谢交流！ 2011-05-10 16:52:58
引用回帖:
Originally posted by bingmou at 2011-04-22 22:58:54:
如果知识测试小体系的话，建议不要提交，直接单机运行。如果是大体系的话，估计一个ecut就可以了，如果全部是金属，400eV足够；如果有O，N可能要700eV。其他的都可以采用赝势中默认值。
如果一定要提交，可以考虑 ...

个人意见：PROGRAMEXEC="/opt/vasp/vasp" 改为 PROGRAMEXEC='脚本位置”是不行的，最后一行其实就是一个运行并行的vasp程序的指令。如果改成脚本位置，那运行的程序也就是vasp都不能调用。只是楼主注意下变量PROGRAMEXEC="/opt/vasp/vasp"设置在你的机器上是正确的就行。还有就是我觉得这个脚本写得罗里罗嗦， 我整合个简单点的 看着清楚的 供楼主参考下：
#!/bin/bash
#PBS -N 5-zmos2nrbig
#PBS -j oe
#PBS -q vasp
#PBS -l walltime=9999:00:00
#PBS -l nodes=2:ppn=8

rm WAVECAR
for a in 400 450 500 550 600 650 700 750
do
cat >INCAR < SYSTEM=CGa
LPLANE=.TRUE.
NPAR=8
LREAL = .TRUE.
PREC    = Accurate
EDIFF   =1E-6
EDIFFG  = -0.03      
IALGO   = 48         
NELMIN  = 4         
ISYM    = 0         
ISPIN = 1
NBANDS = 240
LCHARG  = .TRUE.
LWAVE   = .TRUE.     
LVTOT   = .FALSE.
NBLOCK  = 1         
NSW     = 10
IBRION  = 2
POTIM=0.3
ISMEAR = 0   
SIGMA = 0.05
ENCUT=$a
ISTART=0;ICHARG=2
!
echo "ENCUT=$a eV"
mpirun XXX(此处是mpi的设置参数，根据自己的机器设置） ./vasp > ./out（此处将vasp程序 vasp必须的文件和脚本放同一目录，要嫌麻烦就在脚本里加几个变量把路径写进去）
E=`grep "TOTEN" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5 }'`
echo $a $E >>comment
done
赞一下(2人) 回复此楼
天行健君子以自强不息
5楼2011-05-10 11:04:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 招收调剂硕士 +4 lidianxing 2026-03-19 11/550 2026-03-20 11:05 by lidianxing
[考研] 材料080500调剂求收留 +6 一颗meteor 2026-03-13 6/300 2026-03-20 10:41 by EBSD
[考研] 0856调剂，是学校就去 +7 sllhht 2026-03-19 8/400 2026-03-20 09:30 by 每天只摆一小会
[考研] 0703化学调剂 +10 妮妮ninicgb 2026-03-15 14/700 2026-03-19 22:59 by 学员8dgXkO
[考研] 266求调剂 +5 阳阳哇塞 2026-03-14 10/500 2026-03-19 15:08 by 阳阳哇塞
[考研] 求调剂 +3 Mqqqqqq 2026-03-19 3/150 2026-03-19 14:11 by peike
[考研] 0703化学 305求调剂 +4 FY_yy 2026-03-14 4/200 2026-03-19 05:54 by anny19840123
[考研] 0703化学调剂，求各位老师收留 +10 秋有木北 2026-03-14 10/500 2026-03-19 05:52 by anny19840123
[考研] 331求调剂（0703有机化学 +7 ZY-05 2026-03-13 8/400 2026-03-18 14:13 by 007_lilei
[考研] 070300化学319求调剂 +6 锦鲤0909 2026-03-17 6/300 2026-03-18 13:22 by Iveryant
[考研] 303求调剂 +4 睿08 2026-03-17 6/300 2026-03-18 11:01 by Iveryant
[考研] 278求调剂 +5 烟火先于春 2026-03-17 5/250 2026-03-18 08:43 by 星空星月
[考研] 一志愿苏州大学材料工程（085601）专硕有科研经历三项国奖两个实用型专利一项省级立项 +6 大火山小火山 2026-03-16 8/400 2026-03-17 15:05 by 无懈可击111
[考研] 278求调剂 +3 Yy7400 2026-03-13 3/150 2026-03-17 08:24 by laoshidan
[考研] 一志愿，福州大学材料专硕339分求调剂 +3 木子momo青争 2026-03-15 3/150 2026-03-17 07:52 by laoshidan
[考研] 321求调剂 +5 大米饭！ 2026-03-15 5/250 2026-03-16 16:33 by houyaoxu
[考研] 326求调剂 +3 mlpqaz03 2026-03-15 3/150 2026-03-16 07:33 by Iveryant
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
[考研] 中科大材料专硕319求调剂 +3 孟鑫材料 2026-03-13 3/150 2026-03-14 18:10 by houyaoxu
[考研] 266求调剂 +4 学员97LZgn 2026-03-13 4/200 2026-03-14 08:37 by zhukairuo
信息提示
关闭
请填处理意见
关闭