计算化学公社

标题: 求助关于Gamess 优化结构时读取轨道出错的问题 [打印本页]

作者
Author: Finale 时间: 2019-11-27 10:34
标题: 求助关于Gamess 优化结构时读取轨道出错的问题
本帖最后由 Finale 于 2019-11-27 10:35 编辑

我在使用 Gamess 进行优化结构的时候不收敛，于是将使用 Gaussian 优化成功后生成的 fchk 文件通过 Multiwfn 转化成为 gamess-us 所能读取的轨道格式，但是总是出现 POSSIBLY A DAMAGED OR MANGLED ORBITAL INPUT GROUP? 的错误。不知道各位同学老师是否遇到过这样的问题？请不吝赐教，谢谢！

ps：
1、我读过相关帖子，说是 ”gamess 存在bug：对于存在对称性的分子，gamess可能会把分子输入坐标转动到标准方位（gamess与gaussian的标准方位往往不一致），但初始轨道不做转动，这种情况下采用gaussian的初始轨道无助于收敛，必须去掉对称性，或者禁止分子转动“ —— 为此我在用 gaussian 优化的时候加入了 nosymm 关键字（因为我这个分子是对称性分子），但是仍没有解决问题

2、后面我又尝试了使用 Chemcraft 来打开 fchk 文件，然后 tool -> orbitals -> convert orbitals into input format -> Gamess-us format 来转换成 gamess 的轨道格式，仍然出现相同问题

3、我的体系是一个 Sc 和联苯的配合物，基态是二重态，所以高斯优化的时候用的是非限制性方法，在 gamess 中所用方法仍为 uwb97xd，基组为 6-31G**，inp 和 out 文件见附件

作者
Author: zjxitcc 时间: 2019-11-27 12:53
gjf, fchk文件都没上传，信息太少，看不出两个程序的基函数数目和分子轨道数目。要传轨道，一定要保证两个程序的基函数和分子轨道数目分别相等。先排查这个因素：高斯6-31G(d,p)默认是反人类的6D 7F，为了后续将轨道传给GAMESS，你需要在高斯计算里关键词添加6D 10F，最好再添上int=nobasistransform关键词。如果已经算完6D 7F，只需使用guess=read续算即可，收敛飞速。

另外，GAMESS一些垃圾关键词如TIMLIM=9999999999就不用写了。ISPHER=1和diis=.t. dirscf=.t. soscf=.f. shift=.t. damp=.t.
fdiff=.f. NPUNCH=1这一串什么鬼其实都没用。如果你传过去能很快收敛，用不着这些词。

最后，就算上述问题都攻克了，也未必能保证GAMESS能很好收敛，只能说看运气，因为不同程序的积分格点、同一泛函的实现方式略有差异。如果传过去能很快收敛就万事大吉了。

作者
Author: Finale 时间: 2019-11-27 15:41

zjxitcc 发表于 2019-11-27 12:53
gjf, fchk文件都没上传，信息太少，看不出两个程序的基函数数目和分子轨道数目。要传轨道，一定要保证两个 ...

谢谢您！忘了说了之前算 Gaussian 的时候我已经加上了 6D 10F 但问题仍然没有解决
看了您的回复，我发现了问题所在：两个程序用的基组不同，这导致基函数肯定也不同，十分感谢！

欢迎光临计算化学公社 (http://bbs.keinsci.com/)