gaussian菜鸟常见问题分析

gaussian菜鸟常见问题分析

ID:15957129

大小:62.23 KB

页数:24页

时间:2018-08-06

gaussian菜鸟常见问题分析_第1页
gaussian菜鸟常见问题分析_第2页
gaussian菜鸟常见问题分析_第3页
gaussian菜鸟常见问题分析_第4页
gaussian菜鸟常见问题分析_第5页
资源描述:

《gaussian菜鸟常见问题分析》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库

1、Gaussian菜鸟常见问题分析1.检查是否有初始文件错误在命令行中加入%kJobL301or%kJobL302如果通过则一般初始文件ok。常见初级错误:a.自旋多重度错误b.变量赋值为整数c.变量没有赋值或多重赋值d.键角小于等于0度,大于等于180度e.分子描述后面没有空行f.二面角判断错误,造成两个原子距离过近g.分子描述一行内两次参考同一原子,或参考原子共线2.SCF(自洽场)不收敛则一般是L502错误省却情况做64个cycle迭代(G03缺省128cycles)a.修改坐标,使之合理b.改变初始猜Guess=Huckel或其他的,看Guess关键词。c.增加叠代次数SCFCY

2、C=N(对小分子作计算时最好不要增加,很可能结构不合理)d.iop(5/13=1)这样忽略不收敛,继续往下做。3.分子对称性改变a.修改坐标,强制高对称性或放松对称性b.给出精确的、对称性确定的角度和二面角。如CH4的角度给到109.47122c.放松对称性判据Symm=loosed.不做对称性检查iop(2/16=1)(最好加这个选项)iop(2/16=2)则保持新的对称性来计算4.Opt时收敛的问题a.修改坐标,使之合理b.增加叠代次数optcyc=N5.优化过渡态,若势能面太平缓,则不好找到。iop(1/8=10)默认30(下一个结构和该结构的差别0.3),可改成10。如果每一步

3、都要用到小的步长,应该加opt(notrustupdate)6.在CI(组态)方法中如QCISD(T),CCSD(T),CID方法中,省却最大循环50,若出错(L913错误)解决方法:#PQCISD(maxcyc=N)注:N≤5127.优化过渡态opt=TS(给出过渡态)opt=qst2(给出反应物和产物)opt=qst3(给出反应物和产物和过渡态)a.用G03时的出错opt=ts必须加FC(forceconstant)写法:opt=(TS,calcFc)oropt=(TS,calchffc)计算HF力常数,对QCISD,CCSD等方法用;oropt=(TS,modRedundant)

4、(最好写这个)b.如果计算采用QCISD计算(不好计算FC)则写为QCISDopt=(TS,calcHFFC)(用HF计算FC)8.无法写大的Scratch文件RWFa.劈裂RWF文件%rwf=loc1,size1,loc2,size2,….,locN,-1b.改变计算方法MP2=Direct可以少占硬盘空间c.限制最大硬盘maxdisk=NGB,*MB,有些系统写2GB会出错,可以写2000MB9.FOPT出错原因是变量数与分子自由度数不相等。可用POPT或直接用OPT10.优化过渡态只能做一个STEP原因是负本征数目不对添加iop(1/11)=1或者noeigentest开壳层和闭

5、壳层以及多重度的选择闭壳层计算就是对于多重度是1的体系,此时α和β的电子数目相同,可以把α和β配对,成对的α和β使用同一个轨道,一个轨道上填充2个电子。开壳层计算就是对α和β电子分别计算,一个轨道上只填充1个电子,一般来说,多重度是1时,开壳层计算和闭壳层计算会给出相同的结果。限制开壳层计算是对多重度大于1的体系,此时α和β的电子数目不同,设有m个α和n个β电子,mn,那么让前n个轨道上每个填充一个α和一个β,剩下的m-n个α电子再填充m-n个轨道。即前n个轨道是闭的(每个轨道2个电子),后m-n个轨道是开的(每个轨道1个电子)在高斯中,以HF为例,闭、开、限制开壳层计算分别是RHF,

6、UHF,ROHF。如果只写HF,则按下面的方式取默认方法:对多重度是1的体系,默认为RHF,对多重度大于1的体系,默认是UHF。总结如下图。其中M是多重度。此主题相关图片如下:关于自旋多重度定义:多重度=2S+1,S=n*1/2,n为单电子数。所以,关键是单电子的数目是多少。当有偶数个电子时,例如O2,共有16个电子,那么单电子数目可能是0,即8个alpha和8个beta电子配对,对应单重态,但是也可能是有9个α电子和7个β电子,那么能成对的是7对,还剩2个α没有配对,于是n=2,对应的是多重度3。同理还可以有多重度5,7,9,.一般而言,是多重度低的能量低,最稳定,所以,一般来说,偶

7、数电子的体系多重度就是1。但是也有例外,例如O2就是一个大家都知道的例子,它的基态是三重态,其单重态反而是激发态。所以对于未知的体系,还是算几个保险一点,看哪个能量更低。所以,总结一下,就是电子数目是偶数,未成对电子数目n=0,2,4,6,.自旋多重度是1,3,5,7,.电子数目是奇数,未成对电子数目n=1,3,5,7,.自旋多重度是2,4,6,8,.多数情况是多重度低的能量低,有时(特别是有"磁"性的时候,例如顺磁的O2,以及Fe啊什么的),

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。