计算化学公社

标题: 高斯计算键能问题 [打印本页]

作者
Author: haibin369 时间: 2017-7-27 12:42
标题: 高斯计算键能问题
对于重金属与材料的配位键键能如何计算呢？
在gaussian中尝试计算O2的键能，输入文件如下，但是依然报错（信息如下）。对于简单例子都报错，不知道大体系该如何处理？
输入文件：
%mem=1GB
%chk=BSSE of O2
%nproc=2
#p b3lyp/6-31g* counterpoise=2

BSSE of O2

0,3,0,3,0,3
O  0.0    0.0  0.1  1
O  1.2145  0.0  0.0  2

输出：
>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(UB3LYP) =  -75.0606223086    A.U. after  129 cycles
         NFock=128  Conv=0.10D-05    -V/T= 2.0081
<Sx>= 0.0000 <Sy>= 0.0000 <Sz>= 1.0000 <S**2>= 2.0023 S= 1.0008
<L.S>= 0.000000000000E+00
KE= 7.446006017188D+01 PE=-1.776976220941D+02 EE= 2.817693961366D+01
Annihilation of the first spin contaminant:
S**2 before annihilation    2.0023, after    2.0000
Convergence failure -- run terminated.
Error termination via Lnk1e in D:\G09W\l502.exe at Thu Jul 27 11:23:36 2017.
Job cpu time:    0 days  0 hours  0 minutes 10.0 seconds.
File lengths (MBytes):  RWF=    5 Int=    0 D2E=    0 Chk=    1 Scr=    1

作者
Author: gauss98 时间: 2017-7-27 13:47
sob大神的文章

计算化学键键能时考虑BSSE不仅是多余的甚至是有害的
http://bbs.keinsci.com/forum.php ... &highlight=bsse

作者
Author: sobereva 时间: 2017-7-27 14:37
直接计算三重态氧原子和三重态O2分子，不可能报错

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