24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

牧雪长空

铜虫 (正式写手)

应助: 1 (幼儿园)
金币: 1003.8
散金: 59
帖子: 373
在线: 212.7小时
虫号: 1465005
注册: 2011-10-28
专业: 实验地球化学和计算地球化

[求助] 用侯博士的脚本SIGMA测试出的TS都为0！求大侠指点

为什么我用侯博士的脚本测试的SIGMA值都是零呢，华丽丽的一串零。。如果进入OUTCAR看，会发现，第一个TS其实不是零的。也就是说，应该是我的脚本有错误，可是我又不知道错在哪里。我是在超算上运行的。所以，跟侯博士的原脚本有些差异。下面是我的脚本。
#!/bin/bash
rm WAVECAR
for i in 0.1 0.12 0.14
do
cat > INCAR < SYSTEM = CaCO3
ENCUT = 500
ISTART = 0;ICHARG = 2
NPAR = 4; LREAL = Auto
ISMEAR = 0; SIGMA = $i
PREC = Accurate
LREAL = Auto ; NPAR = 4
!
echo " SIGMA = $i eV "; time /public/software/vasp/vasp-openmpi-ifort/5.3.2/vasp5.3.k
TS=`grep "EENTRO" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5}'`
echo $i $TS >> comment
done

我在运行的时候，是用qsub提交的作业任务，在这个任务脚本里，路径则是这个程序run_sigma的当下路径了。

回复此楼