当前位置: 首页 > 量子化学 >求QST2方法输入文件

求QST2方法输入文件

作者 默默地别动
来源: 小木虫 250 5 举报帖子
+关注

求QST2方法输入文件完整模板,最好使用赝势

 返回小木虫查看更多

今日热帖
  • 精华评论
  • paramecium86

    首先GAUSSIAN的QST方法其实算法很差的。稍微复杂一点儿的过渡态根本就根据反应物和产物猜不出来过渡态初始结构。导致各种会报错。
    建议还是使用 OPT=TS方法直接优化过渡态。而过渡态初始结构 根据反应物结构 用柔性扫描等方法去确定成功率要高很多。
    不过既然你要例子我就给你发一个算例。但是我个人很不推荐用QST 方法。或者说至少不是首选方法。

    算例 La系金属有机物。 对La用SDD 而对其它原子使用6-311+g(d,p), 输入文件中的反应物和产物结构是已经优化过的. 如果你用的g09 就把关键词里的 recalc=3 去掉就可以。

    %nprocshared=20
    %mem=16GB
    %chk=QST2_example.chk
    # opt(qst2,recalc=3) ub3lyp/genecp

    reactant

    0 2
    La                 0.28580300   -0.00811400   -0.00021600
    N                 -2.26630400    0.06638300    0.00171300
    H                 -2.63322500   -0.40365100   -0.83140800
    H                 -2.63760600    1.02030300    0.01176100
    H                 -2.63194400   -0.42064000    0.82558500

    product

    0 2
    La                 0.23245900    0.00480200   -0.00001500
    N                 -1.98036400   -0.17390000   -0.00087500
    H                 -2.58317900   -0.14448700   -0.82076700
    H                  0.64272200    2.06348700    0.01531700
    H                 -2.58298400   -0.16460400    0.81966700

    La 0
    SDD
    ****
    N H 0
    6-311+G(d,p)
    ****

    La 0
    SDD,

  • 默默地别动

    引用回帖:
    2楼: Originally posted by paramecium86 at 2020-05-07 16:04:16
    首先GAUSSIAN的QST方法其实算法很差的。稍微复杂一点儿的过渡态根本就根据反应物和产物猜不出来过渡态初始结构。导致各种会报错。
    建议还是使用 OPT=TS方法直接优化过渡态。而过渡态初始结构 根据反应物结构 用柔性 ...

    老师,我就按照这个输入,但是出现了错误。关键词里需要加freq吗?
    Old curvilinear step not converged, using linear step:
    SCX= 3.94D+01 DXMaxT= 6.18-317 SCLim= 0.00D+00 Fact= 0.00D+00
    RedCar/ORedCr failed for GTrans.
    Error termination via Lnk1e in /home/gauss/local/g09_D01_em64t/g09/l101.exe at Fri May  8 20:05:23 2020.
    Job cpu time:       0 days  0 hours  0 minutes 14.4 seconds.
    File lengths (MBytes):  RWF=      5 Int=      0 D2E=      0 Chk=      1 Scr=      1

  • paramecium86

    引用回帖:
    3楼: Originally posted by 默默地别动 at 2020-05-07 23:23:14
    老师,我就按照这个输入,但是出现了错误。关键词里需要加freq吗?
    Old curvilinear step not converged, using linear step:
    SCX= 3.94D+01 DXMaxT= 6.18-317 SCLim= 0.00D+00 Fact= 0.00D+00
    RedCar/ORedCr ...

    这就是我上一个帖子提到的。QST2的算法很傻。不是你给什么样的初始结构 和 末态结构都能猜得出来的。这个提示就是高斯没猜出来TS结构所以出错。所以一般你自己都猜不出来大概过渡态结构的体系,QST算法也大概率猜不出来。

  • paramecium86

    引用回帖:
    3楼: Originally posted by 默默地别动 at 2020-05-07 23:23:14
    老师,我就按照这个输入,但是出现了错误。关键词里需要加freq吗?
    Old curvilinear step not converged, using linear step:
    SCX= 3.94D+01 DXMaxT= 6.18-317 SCLim= 0.00D+00 Fact= 0.00D+00
    RedCar/ORedCr ...

    计算方法就是这样的,也不需要改什么。根据你给的结构 QST能猜出来,就能开始算。猜不出来非常快就会报错。所以我一直都建议用OPT=TS方法的

  • 默默地别动

    引用回帖:
    5楼: Originally posted by paramecium86 at 2020-05-08 13:43:46
    计算方法就是这样的,也不需要改什么。根据你给的结构 QST能猜出来,就能开始算。猜不出来非常快就会报错。所以我一直都建议用OPT=TS方法的...

    老师,在找过渡态时,除了调节可能是过渡态的两原子距离,还有什么会影响找过渡态?0.01埃的单位长度也不影响吗?高斯小白,求老师解答

猜你喜欢
下载小木虫APP
与700万科研达人随时交流
  • 二维码
  • IOS
  • 安卓