24小时热门版块排行榜

返回列表

当前主题已经存档。

wqd198686

铁杆木虫 (正式写手)

应助: 1 (幼儿园)
金币: 6113.3
散金: 80
红花: 1
帖子: 505
在线: 185.1小时
虫号: 204319
注册: 2006-03-02
专业: 理论和计算化学

[交流] 【转载】Gaussian常见问题分析

来源：厦大理论化学研究中心
http://210.34.15.126/cgi-bin/topic.cgi?forum=3&topic=2061

1．检查是否有初始文件错误
在命令行中加入 %kJob L301 or %kJob L302
如果通过则一般初始文件ok。
常见初级错误：
a. 自旋多重度错误
b. 变量赋值为整数
c. 变量没有赋值或多重赋值
d. 键角小于等于0度，大于等于180度
e. 分子描述后面没有空行
f. 二面角判断错误，造成两个原子距离过近
g. 分子描述一行内两次参考同一原子,或参考原子共线

2． SCF(自洽场)不收敛则一般是L502错误
省却情况做64个cycle迭代(G03缺省128 cycles)
a. 修改坐标，使之合理
b. 改变初始猜  Guess=Huckel 或其他的，看Guess关键词。
c. 增加叠代次数 SCFCYC=N （对小分子作计算时最好不要增加，很可能结构不合理）
d. iop（5/13=1）这样忽略不收敛，继续往下做。

3．分子对称性改变
a. 修改坐标，强制高对称性或放松对称性
b. 给出精确的、对称性确定的角度和二面角。如CH4的角度给到109.47122
c. 放松对称性判据  Symm=loose
d. 不做对称性检查 iop(2/16=1) （最好加这个选项） iop(2/16=2) 则保持新的对称性来计算

4． Opt时收敛的问题
a. 修改坐标，使之合理
b. 增加叠代次数optcyc=N

5．优化过渡态，若势能面太平缓，则不好找到。
iop(1/8=10) 默认30（下一个结构和该结构的差别0.3Å），可改成10。如果每一步都要用到小的步长，应该加opt(notrustupdate)

6．在CI（组态）方法中
如QCISD(T),CCSD(T),CID方法中，省却最大循环50，若出错（L913错误）
解决方法：#P QCISD(maxcyc=N) 注：N≤512

7．优化过渡态
opt=TS (给出过渡态)
opt=qst2 (给出反应物和产物)
opt=qst3 (给出反应物和产物和过渡态)
a.  用G03时的出错 opt=ts 必须加FC (force constant)写法：opt=(TS, calcFc)
or opt=(TS,calchffc)计算HF力常数，对QCISD，CCSD等方法用；
or opt=(TS,modRedundant)  (最好写这个)
b. 如果计算采用QCISD计算（不好计算FC）
则写为QCISD opt=(TS, calcHFFC) (用HF计算FC)

8. 无法写大的Scratch文件RWF
a. 劈裂RWF文件 %rwf=loc1,size1,loc2,size2,……..,locN,-1
b. 改变计算方法 MP2=Direct可以少占硬盘空间
c. 限制最大硬盘 maxdisk=N GB，****MB,有些系统写2GB会出错，可以写2000MB

9. FOPT出错原因是变量数与分子自由度数不相等。可用POPT 或直接用OPT

10. 优化过渡态只能做一个STEP 原因是负本征数目不对添加 iop（1/11）=1或者noeigentest

回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

1楼2007-10-14 21:46:44

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 wqd198686 的主题更新

返回列表