计算化学公社

标题: 过渡态碳氧成键老是有好几个负本征值，不收敛怎么办？ [打印本页]

作者
Author: dlj 时间: 2019-9-28 15:42
标题: 过渡态碳氧成键老是有好几个负本征值，不收敛怎么办？
过渡态中碳氧成键的结构在附件中，一直不收敛，催化剂用的是AgSbF6,一直有两个及多个负本征值，一直不收敛，不然就是碳氧键走开。用的方法基组如下；结构在附件中
%rwf=x.rwf
%nosave
%chk=ts2b-1-b3lyp-dzvp-opt-sol.chk
%mem=80MW
%nproc=1
#p b3lyp/dgdzvp scfcon=6 iop(1/8=6) opt=(ts,calcfc,noeigentest,nrscale,
modredundant) volume=tight scrf=(read,PCM,solvent=Ether)
freq=noraman int(grid=99590)

ts2b-1-b3lyp-dzvp-opt1-sol stepn=23

0 1
  6          0    -3.712321 -3.790445 0.692645
  6          0    -3.545296 -3.349894 2.153661
  6          0    -2.334562 -2.466642 2.399325
  6          0    -1.635439 -1.903874 1.324192
  6          0    -1.945761 -2.422012 -0.070152
  6          0    -1.986816 -2.135321 3.713492
  6          0    -0.677868 -0.892479 1.565673
  6          0    -0.324971 -0.538374 2.905752
  6          0    -0.971078 -1.211486 3.971442
  6          0    -1.515186 -1.580886 -1.291700
  1          0    -0.572870 -1.971279 -1.682795
  1          0    -4.432568 -2.778670 2.459566
  8          0    -0.576087 -0.949428 5.252085
  1          0    -1.066593 -1.532830 5.864990
  6          0       0.696154 0.510053 3.140180
  8          0       1.692748 0.715008 2.438346
  8          0       0.443694 1.278973 4.210724
  6          0       1.435120 2.283317 4.527882
  1          0       1.075634 2.760366 5.439340
  1          0       2.409105 1.818368 4.694879
  1          0       1.513377 3.011391 3.717157
  1          0    -2.508136 -2.602092 4.548536
  1          0    -3.511572 -4.235163 2.799428
  8          0    -3.342624 -2.720120 -0.196942
  1          0    -3.032783 -4.635676 0.494328
  8          0    -0.161274 -0.214363 0.537429
  1          0    -1.405306 -3.374731 -0.189637
  1          0    -2.272132 -1.741536 -2.071189
  6          0    -1.322733 -0.107243 -1.150507
  8          0    -0.412271 0.391342 -1.959071
  6          0    -2.371829 0.941844 -0.841452
  6          0    -0.503537 1.863226 -2.031282
  6          0    -1.569582 2.252975 -1.002904
  1          0       0.490464 2.256876 -1.829360
  1          0    -0.787633 2.092999 -3.062454
  1          0    -1.077246 2.469961 -0.049919
  6          0    -2.434314 3.450744 -1.404988
  1          0    -3.242158 3.596755 -0.680677
  1          0    -2.927427 3.267514 -2.369895
  6          0    -1.654371 4.754996 -1.506336
  6          0    -2.420367 6.092688 -1.714452
  8          0    -0.435033 4.762522 -1.422384
  6          0    -1.598157 7.012966 -2.633003
  1          0    -0.597884 7.182025 -2.232503
  1          0    -1.501508 6.559420 -3.624716
  1          0    -2.097854 7.978372 -2.739294
  6          0    -2.563142 6.754217 -0.328700
  6          0    -3.843167 5.898244 -2.295479
  1          0    -3.789526 5.437298 -3.287092
  1          0    -4.480550 5.282142 -1.657450
  1          0    -4.326724 6.872840 -2.412484
  8          0    -2.324484 7.919790 -0.083903
  8          0    -3.038286 5.881639 0.588100
  6          0    -3.245552 6.408920 1.922502
  1          0    -3.610226 5.568275 2.511965
  1          0    -2.304646 6.786146 2.328274
  1          0    -3.984494 7.212634 1.898892
  6          0    -5.143378 -4.212513 0.367641
  1          0    -5.470165 -4.920243 1.141121
  1          0    -5.791988 -3.330875 0.458150
  6          0    -5.310650 -4.850700 -1.018373
  1          0    -4.950608 -4.152778 -1.783682
  1          0    -4.669255 -5.740671 -1.084193
  6          0    -6.761003 -5.248511 -1.327860
  1          0    -7.125326 -5.935523 -0.550682
  1          0    -7.402658 -4.357640 -1.271347
  6          0    -6.932954 -5.907835 -2.704104
  1          0    -6.568507 -5.222823 -3.481546
  1          0    -6.293993 -6.799404 -2.759795
  6          0    -8.383924 -6.301569 -3.010836
  1          0    -9.043892 -5.425904 -3.001211
  1          0    -8.767479 -7.012599 -2.269415
  1          0    -8.470007 -6.771461 -3.996960
51          0       5.309273 -1.337984 -1.024451
  9          0       5.604979 0.426279 -0.256438
  9          0       6.276991 -0.837875 -2.628800
  9          0       3.619026 -0.617564 -1.769230
  9          0       4.894079 -3.078784 -1.773916
  9          0       6.883589 -2.037331 -0.134986
  9          0       4.196801 -1.761341 0.551510
47          0       2.168892 -0.231240 0.255274
  1          0    -3.050879 0.841329 -1.708250
  6          0    -3.220250 0.834580 0.424906
  1          0    -3.727562 -0.130527 0.486445
  1          0    -2.616447 0.972743 1.323648
  1          0    -3.989465 1.613137 0.407577

  29 26

作者
Author: sobereva 时间: 2019-9-29 02:15
你这输入文件全是错误，病入膏肓，算出来的结果根本没法用

垃圾关键词太多，仔细看
常见的多余的和被滥用的Gaussian关键词
http://sobereva.com/331（http://bbs.keinsci.com/thread-3460-1-1.html）

拿dgdzvp描述Ag、Sb还不用相对论哈密顿纯属瞎算，就算算出来想要的结果，最后碰上内行审稿人肯定要挨批，让你重算一遍

对Ag、Sb这样的重元素该用什么仔细看
谈谈赝势基组的选用
http://sobereva.com/373（http://bbs.keinsci.com/thread-5625-1-1.html）

写int(grid=99590)毫无意义，纯粹瞎浪费时间

写volume=tight干什么，要算体积用Multiwfn算严格准确得多
谈谈分子体积的计算
http://sobereva.com/102

scrf里写read还不写要自定义的溶剂设置完全莫名其妙

PCM本来就是默认的关键词，完全多余。怎么正确使用溶剂模型仔细看
谈谈隐式溶剂模型下溶解自由能和体系自由能的计算
http://sobereva.com/327（http://bbs.keinsci.com/thread-3345-1-1.html）

写scfcon=6等于对数据不负责任，除非你最后再默认收敛限下重新优化。这还容易带来其它问题。而且做振动分析更是不能接受这样的SCF收敛限

你这体系存在显著弱相互作用还不加D3，很容易挨批，仔细看
谈谈“计算时是否需要加DFT-D3色散校正？”
http://sobereva.com/413（http://bbs.keinsci.com/thread-9772-1-1.html）

没事别冻结

关键词改成B3LYP/genecp em=GD3BJ opt(TS,calcfc,noeigen) scrf=solvent=ether freq

末尾空一行写
C H O F
6-31G*
****
Ag Sb
SDD
****

Ag Sb
SDD

如果优化不到期望的结果到时候再说。另参考
简谈Gaussian里找过渡态的关键词opt=TS和QST2、3
http://sobereva.com/460（http://bbs.keinsci.com/thread-12093-1-1.html）

作者
Author: dlj 时间: 2019-9-29 14:44
谢谢您的解答，我会改正后算的。

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