过渡态优化有明显负值eigenvalue，最后频率计算无虚频是什么原因？

brad_ling

用Gaussian16 进行过渡态优化任务关键词
#p opt=(calcfc,ts,noeigen,nofreeze) freq 6-311+g(d,p) guess=read  um062x


任务正常结束
Normal termination of Gaussian 16 at Thu May 16 10:10:19 2024.
Normal termination of Gaussian 16 at Thu May 16 10:14:53 2024.


$ grep  'Eigenvalues ---   -' AA.out
     Eigenvalues ---   -0.02758   0.00022   0.00141   0.00385   0.00775
     Eigenvalues ---   -0.02550   0.00022   0.00139   0.00392   0.00737
     Eigenvalues ---   -0.02587   0.00022   0.00143   0.00392   0.00807
     Eigenvalues ---   -0.02598   0.00022   0.00149   0.00392   0.00876
     Eigenvalues ---   -0.02565   0.00022   0.00151   0.00392   0.00899
     Eigenvalues ---   -0.02534   0.00022   0.00155   0.00392   0.00899
     Eigenvalues ---   -0.02515   0.00022   0.00159   0.00392   0.00899
     Eigenvalues ---   -0.02498   0.00022   0.00159   0.00391   0.00898
     Eigenvalues ---   -0.02488   0.00021   0.00163   0.00390   0.00868
     Eigenvalues ---   -0.02495   0.00021   0.00166   0.00389   0.00819
     Eigenvalues ---   -0.02510   0.00022   0.00168   0.00388   0.00778
     Eigenvalues ---   -0.02522   0.00023   0.00168   0.00386   0.00746
     Eigenvalues ---   -0.02519   0.00025   0.00167   0.00378   0.00698
     Eigenvalues ---   -0.02503   0.00029   0.00167   0.00359   0.00632
     Eigenvalues ---   -0.02482   0.00036   0.00166   0.00326   0.00585
     Eigenvalues ---   -0.02479   0.00040   0.00166   0.00297   0.00569
     Eigenvalues ---   -0.02492   0.00036   0.00166   0.00287   0.00573
     Eigenvalues ---   -0.02496   0.00027   0.00166   0.00288   0.00580
     Eigenvalues ---   -0.02494   0.00025   0.00166   0.00298   0.00581
     Eigenvalues ---   -0.02491   0.00021   0.00166   0.00303   0.00578
     Eigenvalues ---   -0.02489   0.00015   0.00169   0.00304   0.00580
     Eigenvalues ---   -0.02489   0.00009   0.00178   0.00302   0.00577
     Eigenvalues ---   -0.02486  -0.00034   0.00216   0.00303   0.00582
     Eigenvalues ---   -0.02485  -0.00111   0.00279   0.00311   0.00581
     Eigenvalues ---   -0.02487  -0.00158   0.00284   0.00346   0.00591
     Eigenvalues ---   -0.02486  -0.00128   0.00284   0.00366   0.00594
     Eigenvalues ---   -0.02486  -0.00027   0.00281   0.00366   0.00594
     Eigenvalues ---   -0.02498   0.00038   0.00190   0.00313   0.00505
     Eigenvalues ---   -0.02465   0.00025   0.00092   0.00301   0.00485
     Eigenvalues ---   -0.02449   0.00035   0.00140   0.00287   0.00485
     Eigenvalues ---   -0.02414   0.00043   0.00152   0.00286   0.00481
     Eigenvalues ---   -0.02389   0.00050   0.00171   0.00282   0.00475


最后频率计算无虚频
grep 'Frequencies' AA.out
Frequencies --     12.3121                38.9422                52.7731
Frequencies --     56.1296                63.1460                93.5382
Frequencies --    101.7566               141.2316               151.0076
Frequencies --    172.2848               198.9081               214.4991
Frequencies --    230.8935               242.8249               317.4377
Frequencies --    363.3635               373.3726               413.1285
Frequencies --    433.1921               486.1094               522.3230
Frequencies --    570.9798               608.1146               620.0899
Frequencies --    658.2377               704.2252               723.2934
Frequencies --    744.7122               779.2446               788.0564
Frequencies --    866.7618               929.2262               951.6038
Frequencies --    961.2149              1021.7574              1024.4368
Frequencies --   1052.3150              1083.9827              1089.6606
Frequencies --   1120.7711              1135.4677              1253.1965
Frequencies --   1327.1779              1369.7424              1439.0630
Frequencies --   1479.0355              1492.7965              1554.3011
Frequencies --   1590.6363              1645.2896              1679.2185


结构上明显也是一个过渡态。请问这是什么原因？

Freeman

那些对应的是还没收敛时近似的hessian本征值

努力学计算

Freeman 发表于 2024-5-17 15:50
那些对应的是还没收敛时近似的hessian本征值

请问老师，我是固定了底层原子来搜的过渡态，成功搜出来了。然后用的这个固定底层原子的模型结构计算频率，在outmol文件里也显示finished successfully，但频率分析的时候，强度那一列都是0，这是因为固定了底层，所以强度都为0嘛，但outmol文件最后提示成功完成，所以我不太明白是不是与底层固定有关，计算频率是不是需要都不固定呀？就是只固定底层原子也不行

Freeman

努力学计算 发表于 2024-7-24 08:51
请问老师，我是固定了底层原子来搜的过渡态，成功搜出来了。然后用的这个固定底层原子的模型结构计算频率 ...

这是个新问题了，应该开个新贴。