小白第一次做计算,gaussian 结构优化求助
输出文件结果与文献结果相差巨大。
scf done: e(rb3lyp) = -3744.02980095 a.u. after 129 cycles
nfock=128 conv=0.24d+01 -v/t= 3.2022
convergence failure -- run terminated.
error termination via lnk1e in e:\gs09\g09w\l502.exe at wed jul 01 15:38:09 2020.
job cpu time: 0 days 0 hours 1 minutes 45.0 seconds.
file lengths (mbytes): rwf= 5 int= 0 d2e= 0 chk= 1 scr= 1
文献原文:通过构型优化,频率验证,单点能计算得出合理的几何构型鞍点结构和基于该构型的最小能量值。通过计算吸附或键合前后的能量差值,考察反应物分子在催化剂表面的吸附特性。
(dft-b3lyp参数方法,au采用lanl2dz基组,cl采用6-31g(d,p)基组计算。)
在用gaussian09重复时是否还需要设置一些什么?(输入文件、输出文献已附)
今日热帖
这个错误本身的意思是 scf不收敛。遇到这种情况常规的应对方法可以看 我这个教程里的 http://muchong.com/t-14242735-1 第2.2节,如图所示
如果你检查到教程中c部分 ,就可以发现是因为你的基组写的不对。你对au原子用了lanl2dz基组,但是基组的赝势部分你得单独定义。你可以看同样的教程中的1.1节 就是讲怎么写混合基组和赝势的。 截图也附上了
502错误_page_1.jpg
502错误_page_2.jpg
赝势写法.png,
非常感谢您的回答,但是我写入赝势基组后(如图)仍不收敛,并且在提前进行了uff优化后,结果也是不收敛。自旋多重度2、4、6、8·····都尝试了一下,均不收敛。
还烦请再次抽时间指教,谢谢。
批注 2020-07-02 112744.png
你输入文件下面的赝势格式写的并不对, 你不能分开写赝势。要把它们写在一起。比如这样
Au Cu 0
lanl2dz
改好之后如果运行还不收敛 就用截图里后面那些关键词加上试一试。关键词那么多呢 你都可以试着用用