气相中过渡态振荡不收敛

夜千

方法是M062x，液相中已经找到此过渡态，气相中要么直接是产物，要么非常规律的振荡不收敛。已经按照教程加入了不同的关键词，包括calcall,gdiis,maxstep等，仍然不收敛。请问现在该怎么修改命令呢
这是我的输入文件
# M062x/6-31G(d,p) opt(calcall,maxcycle=200,TS,noeigen) freq=noraman integral=ultrafine temperature=313.15

Title Card Required

0 1
C                  3.61576900   -1.07496700   -0.47556600
C                  3.07262300    0.22416900   -0.29557000
C                  3.93548900    1.31743600   -0.37023500
C                  5.28982900    1.10763200   -0.64393300
C                  5.79382400   -0.17596300   -0.83268700
C                  4.95608800   -1.29151800   -0.74364600
C                  1.38742900   -1.40281900   -0.07171000
H                  3.58620300    2.32948300   -0.20630800
H                  5.95647800    1.96137700   -0.70553100
H                  6.84855600   -0.31685200   -1.04463700
H                  5.34388900   -2.29680200   -0.87472600
N                  2.60797400   -2.01921400   -0.32822400
C                  2.75126900   -3.43623900   -0.53904800
H                  3.39094100   -3.88892800    0.22617200
H                  1.75267900   -3.87088000   -0.48053000
H                  3.18256900   -3.63823100   -1.52509400
O                  0.32366900   -2.01400100    0.06169900
C                 -0.82854000   -1.17041200    2.62301500
C                 -0.09489100   -0.12876400    3.01716900
C                  0.73778100    2.20416200    0.04622400
C                 -0.50731800    2.99518000    0.17787900
F                 -0.24258100    4.27814700    0.43778900
F                 -1.28859800    2.97580700   -0.93366600
F                 -1.31060100    2.53103500    1.16761000
C                  1.66352900    0.02813000   -0.02646700
N                  0.65601900    0.90216800    0.11052500
C                 -1.81669300   -1.08791700    1.54542000
C                 -2.93872100   -1.92081900    1.49746500
C                 -2.50504900    0.00773600   -0.41033700
C                 -3.86149000   -1.77483900    0.47654500
H                 -3.07363000   -2.67199800    2.26633300
C                 -3.63382400   -0.78709500   -0.47372600
H                 -2.28278600    0.77160000   -1.14885500
H                 -4.74686900   -2.39511500    0.39915300
N                 -1.64056200   -0.15997200    0.58614300
N                 -4.58381900   -0.58967100   -1.57050700
O                 -5.54978500   -1.32716700   -1.60303800
O                 -4.33751500    0.29711600   -2.36381000
H                 -0.20587600    0.85436000    2.56944300
H                 -0.74005300   -2.13831600    3.10647000
H                  1.63617500    2.74210800   -0.22224700
H                  0.62691800   -0.23292500    3.81898800
H                 -0.57003200    0.34690000    0.42665900

Shaoqz

震荡不收敛的有优化出一个yes吗？（RMS Force项）感觉可能是气相TS和液相差的比较多所以初结构不合理了

夜千

Shaoqz 发表于 2019-5-11 20:44
震荡不收敛的有优化出一个yes吗？（RMS Force项）感觉可能是气相TS和液相差的比较多所以初结构不合理了

我所上传的输入文件中没有一个yes，但是如果输入命令改为# M062x/6-31G(d,p) opt(calcfc,maxcycle=200,TS,noeigen,gdiis,maxstep=3,notrust) freq=noraman integral=ultrafine temperature=313.15
就只有RMS有一个yes

Shaoqz

夜千 发表于 2019-5-11 21:49
我所上传的输入文件中没有一个yes，但是如果输入命令改为# M062x/6-31G(d,p) opt(calcfc,maxcycle=200,TS ...

有试过把rms force比较小且结构合理的作为新的初结构重投吗

sobereva

气相下和溶剂模型下势能面可能存在明显不同，未必溶剂下存在的TS在气相下就有，或者可能都有但结构偏离较大。不确定的话可以做势能面扫描辅助判断

写=noraman完全多余
常见的多余的和被滥用的Gaussian关键词
http://sobereva.com/331（http://bbs.keinsci.com/thread-3460-1-1.html）

夜千

Shaoqz 发表于 2019-5-12 01:05
有试过把rms force比较小且结构合理的作为新的初结构重投吗

有，仍然振荡

夜千

sobereva 发表于 2019-5-12 10:05
气相下和溶剂模型下势能面可能存在明显不同，未必溶剂下存在的TS在气相下就有，或者可能都有但结构偏离较大 ...

谢谢老师，我用B3LYP找到了此过渡态的气相，能否说明他在气相中是存在的呢，只是可能由于M062找过渡态弱一点

sobereva

夜千 发表于 2019-5-13 20:56
谢谢老师，我用B3LYP找到了此过渡态的气相，能否说明他在气相中是存在的呢，只是可能由于M062找过渡态弱 ...

对有机体系，如果M06-2X找不到，底气不足。至少用其它的靠谱的比如wB97XD也能找到才比较容易说得过去