[交流] 【求助】Gaussian09作业中断，可否接着断点继续计算？ 已有6人参与

#!/bin/sh

if [ $# -ne 1 ]
then
        echo "Usage: $0 Gaussian output file"
        exit 1
fi

if [ ! -e $1 ]
then
        echo "$1 not existed!"
        exit 1
fi

scf_iters=`grep "SCF Done" $1 | wc -l`
natoms=`grep "NAtoms" $1 | awk '{print $2}'| head -1`
left_iters=$(($natoms*6+1-$scf_iters))
echo ""
echo "The SCF has runned for $scf_iters cycles"
echo -e "Still need \033[45m\033[5m $left_iters \033[0m cycles"
echo ""

Originally posted by xmdxn at 2010-10-16 14:27:36:


我现在又换了6-31+g*和6-31g*继续优化和频率分析，激发态的优化很快，一般两个小时就可以做出来，但是频率分析还是要跑很久，目前跑了两千五百分钟还是没有跑出来，不停地出现
                             ...

Originally posted by xmdxn at 2010-10-16 18:18:47:


原来是数值梯度哦。Gaussian09说明书上写的应该是说一阶导为积分的。
谢谢这位虫友，应该可以接着算，利用rwf文件可以的，似乎chk文件也可以。
用ultraedit搜索字符串也可以看出来scf done了几次O(∩_∩)O~ ...

Originally posted by xmdxn at 2010-10-18 15:00:11:
ls这位虫友，我用6-31g*重新做频率分析，终于跑出来了。我的体系是C20H12，32个原子，但是scf迭代了44次。照6N+1来算，不止44次啊。。。又困惑了，G09的TD频率分析的梯度到底是什么状况哦。。。