Gaussian中稳定波函数后的fchk文件转gbw文件求助

hdhxx123

老师们好，我今天在Linux系统下用Gaussian对一个分子进行初步的波函数稳定，但稳定后却无法将fchk文件转变为gbw，一直卡在orca_2mkl的那一步，不知道有什么办法可以解决，以及不知道出了什么问题，附件里是对应的gjf、out与fchk文件，老师们能不能帮我看一下是什么情况，谢谢老师们！

hdhxx123

zjxitcc 发表于 2020-12-12 14:07
你输入文件中关键词
#p b3lyp/sdd stable=opt 5d em=gd3bj int=(nobasistransform,acc2e=12) nosymm scf(q ...

谢谢老师！我一开始是用的下面那些def2svp,ma-def2sv(p)，用了一堆办法不收敛，qc都不收敛，没办法只能转SDD，我再想想用其他方法先得到收敛初猜再读看看

zjxitcc

你输入文件中关键词
#p b3lyp/sdd stable=opt 5d em=gd3bj int=(nobasistransform,acc2e=12) nosymm scf(qc,conver=6,vshift=500) scrf=smd
坐标底下却写了自定义基组和赝势数据。

需仔细阅读sob老师博文《详解Gaussian中混合基组、自定义基组和赝势基组的输入》http://sobereva.com/60

一、这关键词写的十分不合理，不建议传给ORCA。你在#行写了SDD，导致（1）坐标的底下的基组和赝势全没用；（2）所有原子的基组变成了SDD。然而SDD默认是(6D,10F)，你写了5D，但你忘记写7F，导致算的仍然是10F。另外，非金属原子应照常用合适基组，不要用SDD。
二、哪有一上来就写qc,conver=6的，这收敛的波函数精度低且不一定可靠（尽管经过了stable=opt检验）。正常高斯算完不能收敛的话，也有很多收敛技巧选择，十分不建议直接就写qc,conver=6。
三、不是错误，但需知道：这是单点计算，em=gd3bj不影响波函数（即不影响轨道系数），写没写对传轨道没有任何影响。

待你改正完之后还有问题的话，再上来发帖询问。若文件正确，但还是慢，那是orca_2mkl的问题，这个程序本身写的不好，比较慢。