计算化学公社

标题: 气相中过渡态振荡不收敛 [打印本页]

作者
Author: 夜千 时间: 2019-5-11 17:39
标题: 气相中过渡态振荡不收敛
本帖最后由夜千于 2019-5-11 22:02 编辑

方法是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

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

作者
Author: 夜千 时间: 2019-5-11 21:49

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

作者
Author: Shaoqz 时间: 2019-5-12 01:05

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

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

作者
Author: sobereva 时间: 2019-5-12 10:05
气相下和溶剂模型下势能面可能存在明显不同，未必溶剂下存在的TS在气相下就有，或者可能都有但结构偏离较大。不确定的话可以做势能面扫描辅助判断

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

作者
Author: 夜千 时间: 2019-5-13 14:43

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

有，仍然振荡

作者
Author: 夜千 时间: 2019-5-13 20:56

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

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

作者
Author: sobereva 时间: 2019-5-14 03:26

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

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

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