金属有机化合物催化酯化反应过渡态搜索，结构一直震荡不收敛

liangcx

金属有机化合物催化反应过渡态搜索，发现结构一直震荡不收敛
输入文件相关参数如下：
# opt=(ts,modredundant,gdiis,noeigentest,maxstep=3,notrust) freq b3lyp/gen
em=gd3 optcyc=200 pseudo=read

Title Card Required

0 1
C                 -0.96112310    1.32829369    0.00000000
C                 -2.33098110    1.55971169   -0.03189300
C                 -2.87201710    2.37490369   -1.03563400
C                 -2.03780510    2.95512969   -2.00074400
C                 -0.66623910    2.72119369   -1.96577100
C                 -0.12363710    1.90476669   -0.96517100
H                 -0.51631410    0.70492769    0.76792200
H                 -2.97971210    1.13063269    0.72108400
H                 -2.47614610    3.58645869   -2.76567300
H                 -0.01545310    3.16894069   -2.70767900
C                 -4.32497310    2.66316369   -1.16809800
C                  1.33675690    1.61785569   -0.87715700
O                 -4.86603210    3.47566769   -1.84454700
O                  1.84883790    0.91400569   -0.03047400
O                 -5.37935561    0.19495474   -1.56926196
H                 -4.70078661   -0.49194726   -1.61106296
O                  2.05328090    2.22867269   -1.85085300
H                  2.98047690    1.98000369   -1.70306100
C                 -5.45672810    3.05822242    2.40994638
H                 -6.51475251    3.29958659    2.54420201
H                 -5.11398574    2.38184158    3.19933077
C                 -4.63904173    4.34415947    2.43944698
H                 -3.56556385    4.14506129    2.35168665
H                 -4.90375647    5.03400526    1.63695276
C                 -4.91458119    5.06058749    3.76167714
H                 -4.73431527    4.38120435    4.61363634
H                 -5.97556965    5.35225343    3.79076942
O                 -5.39647935    2.25045400    1.18354697
H                 -6.03442497    0.73897955    0.48210948
O                 -4.05217392    6.19043515    3.81599509
H                 -4.38484751    6.78807235    4.49616233
C                 -8.44232353    3.45842844   -4.20777227
C                 -9.57724078    4.37753989   -3.73285136
C                 -9.12518085    5.37637146   -2.65460538
H                 -8.83412511    2.75619572   -4.96169828
H                 -9.97053779    4.94065364   -4.59184585
H                 -8.74546539    4.80788001   -1.78584974
O                 -8.15664473    6.28627110   -3.12934288
O                 -7.34251384    4.15382303   -4.77137305
H                 -8.10305917    2.84403549   -3.35637283
H                -10.39724924    3.75806637   -3.34293322
H                 -9.99766131    5.94347116   -2.29657242
C                 -6.33609359    8.12595987   -5.91986532
C                 -4.42934331    7.11667493   -4.53781800
C                 -4.88820022    7.65921192   -5.89924189
H                 -6.58868235    8.52826123   -6.90338741
H                 -4.41891321    7.96088965   -3.82492338
H                 -4.85781853    6.91198456   -6.70011500
Al                -6.91092597    5.77534128   -4.24167759
O                 -5.21319085    6.06645999   -4.04000179
O                 -7.17067698    6.94811673   -5.72363781
H                 -3.38743969    6.77616644   -4.63230977
H                 -4.23466963    8.48937994   -6.19590253
H                 -6.57647219    8.86346576   -5.14429243
H                 -8.07730983    7.04462117   -6.14007411

B 11 15 D
B 15 29 D
B 29 28 D
B 28 11 D

Al 0
Lanl2Dz
****
H C O
6-31g*
****

Al 0
Lanl2Dz

输出文件报错如下：
Error termination request processed by link 9999.
Optimization stopped.
    -- Number of steps exceeded,  NStep= 200

sobereva

甭用modredundant试图不去计算精确的初始Hessian，这会导致优化过渡态成功几率大降
Al用赝势基组毫无意义

并且监控结构变化，别光看曲线图。结构收敛到了离谱的结构再怎么跑都是白搭

liangcx

sobereva 发表于 2024-10-20 14:34
甭用modredundant试图不去计算精确的初始Hessian，这会导致优化过渡态成功几率大降
Al用赝势基组毫无意义
...

好的，谢谢sob老师解答