计算化学公社

标题: 求助Gaussian09 TS方法计算过渡态不收敛 [打印本页]

作者
Author: 忘不掉的回忆 时间: 2022-5-25 10:50
标题: 求助Gaussian09 TS方法计算过渡态不收敛
本帖最后由忘不掉的回忆于 2022-5-26 09:56 编辑

老师您好，我用TS方法求过渡态总是不收敛麻烦老师看下

%mem=10GB
%nprocshared=48
%chk=guodutai111.chk
#p opt=(calcfc,ts,noeigen,cartesian,gdiis) freq ub3lyp/genecp em=gd3bj

Title Card Required

0 2

Mg             -1 2.19973600 -1.82960800 0.86639300
Mg             -1 5.47384100 -0.40038700 0.45787400
Mg             0 -0.83136600 -0.08755200 -0.02829800
Mg             -1 -6.88157400 0.38111700 -0.91546700
Mg             -1 -3.97917800 -1.43860300 0.17943700
Cl             0 1.47761600 -0.52520800 -1.18696500
Cl             0 0.18067600 2.12558800 0.37978500
Cl             -1 0.03365000 -1.07092100 1.95945600
Cl             0 -5.15969700 -0.74038200 -2.03734100
Cl             -1 -1.81096000 -2.19082400 -0.91187800
Cl             0 3.42654900 0.19565200 1.81090300
Cl             0 4.71947500 1.25122600 -1.13006400
Cl             -1 7.64164100 -1.15630200 -0.63474300
Cl             -1 4.36646300 -2.58318500 -0.22570500
Cl             -1 -9.04795800 1.13632200 0.17654400
Cl             0 -2.96961500 0.67541700 0.67749900
Cl             -1 -6.14490700 -0.68090400 1.27241100
Ti             0 2.46971100 1.66107700 -0.28307500
C             0 3.16117400 3.89482700 0.46173200
H             0 4.19832600 4.25688600 0.47413200
H             0 3.21846300 2.96456800 1.11824400
H             0 2.51152400 4.55912200 1.04767700
C             0 2.63775000 3.60173500 -0.92518200
H             0 3.37294700 3.71224300 -1.73133000
H             0 1.67415700 4.06437900 -1.16796900
C             0 2.35661133 4.20305262 -2.31835919
C             0 1.11305154 3.67936848 -2.19233927
H             0 2.49284108 5.26434422 -2.31709421
H             0 3.20223582 3.55523656 -2.41912417
C             0 0.91698240 2.15190209 -2.19416132
C             0 1.98800143 1.33141579 -2.32178383
H             0 -0.06487214 1.73842660 -2.09466421
H             0 1.85177130 0.27012426 -2.32305244
H             0 2.96985616 1.74489126 -2.42127930
C             0 -0.10401544 4.61173925 -2.04731280
H             0 -0.96694619 4.13999398 -2.46886578
H             0 -0.27791827 4.81262211 -1.01082648
H             0 0.08719367 5.53041774 -2.56148097

Mg Cl C H
def2svp
****
Ti
SDD
****

Ti
SDD

>>>>>>>>>> Convergence criterion not met.

SCF Done:  E(UB3LYP) =  -6854.47252012    A.U. after  129 cycles

         NFock=128  Conv=0.37D-05    -V/T= 2.0078

<Sx>= 0.0000 <Sy>= 0.0000 <Sz>= 0.5000 <S**2>= 0.7624 S= 0.5062

<L.S>= 0.000000000000E+00

KE= 6.801193950916D+03 PE=-2.748206648241D+04 EE= 8.222181247925D+03

Annihilation of the first spin contaminant:

S**2 before annihilation    0.7624, after    0.7501

Convergence failure -- run terminated.

Error termination via Lnk1e in /home[img]file:///C:\Users\李康康\AppData\Roaming\Tencent\QQTempSys\{)_G62``_[8)(8ND1X0_M7L.gif[/img]mba/g09/l502.exe at Wed May 25 06:21:56 2022.

Job cpu time:    4 days 10 hours 28 minutes 36.4 seconds.

File lengths (MBytes):  RWF= 2167 Int=    0 D2E=    0 Chk=    28 Scr=    1

作者
Author: 一颗赛艇 时间: 2022-5-25 11:58
你算高斯不看结构合不合理就扔进去算吗？

作者
Author: snljty2 时间: 2022-5-25 12:20
你这建模...这画的都是个啥结构...各种二面角完全不对，按照你的有机化学和无机化学基本知识，先把模型建合理了，熟悉基本的优化到极小值结构计算等以后，再去上过渡态计算这种明显需要更多经验和chemical insight的任务。

作者
Author: RAL 时间: 2022-5-25 14:15
这看上去怕是想用簇模型研究MgCl2表面掺杂Ti催化的某种反应...先看看sob老师是怎么做的吧，至少你取的表面太小了
《使用量子化学程序基于簇模型计算金属表面吸附问题》http://sobereva.com/540

作者
Author: sobereva 时间: 2022-5-26 03:42
以后记得输入文件别直接贴，而应当作为附件上传
当前空行明显都不合理，赝势前头空两行的话赝势信息都根本读不进去。你不上传输入文件别人也不知道是你发帖时候是不注意导致出现的额外空行还是真就是空行数不对

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