当前位置: 首页 > 量子化学 >求助:寻找DME与O2反应的过渡态

求助:寻找DME与O2反应的过渡态

作者 明8221
来源: 小木虫 450 9 举报帖子
+关注

求助:寻找DME与O2反应的过渡态,通过TS,QST2和QST3都无法收敛,求各位高手帮我找找解决办法,谢谢!
%chk=C:\Users\ming\Desktop\DMMnCal\B3LYP\CH3OCH3O2f3.chk
# opt=(calcall,qst3,maxcycle=1000,noeigentest) hf/3-21g geom=connectivity
scf=maxcycle=1000

Title Card Required

0 1
C                 -2.10732978   -0.07853403    0.00000000
H                 -1.10221849   -0.47366697    0.00461120
H                 -2.26586245    0.51304858    0.89748995
H                 -2.23629236    0.56768813   -0.86393050
O                 -2.99778413   -1.20073030   -0.04976905
C                 -4.38424039   -0.84240075   -0.06153636
H                 -4.63135615   -0.23696566   -0.92899566
H                 -4.94620999   -1.76314268   -0.09915401
H                 -4.66033960   -0.29129183    0.83297283
O                  0.26481359   -0.92623725    0.00646600
O                  1.19372436   -0.10842438   -0.09283296

1 2 1.0 3 1.0 4 1.0 5 1.0
2
3
4
5 6 1.0
6 7 1.0 8 1.0 9 1.0
7
8
9
10 11 2.0
11

Title Card Required

0 1
C                 -1.74459848   -0.15056216    0.51674170
H                 -0.62780816   -0.58959876    0.52186526
H                 -1.90313115    0.44102045    1.41423165
H                 -1.87356106    0.49566000   -0.34718880
O                 -2.63505283   -1.27275843    0.46697265
C                 -4.02150909   -0.91442888    0.45520534
H                 -4.26862485   -0.30899379   -0.41225396
H                 -4.58347869   -1.83517081    0.41758769
H                 -4.29760830   -0.36331996    1.34971453
O                  0.39931522   -0.92338973    0.52293355
O                  1.32822599   -0.10557686    0.42363459

1 2 0.5 3 1.0 4 1.0 5 1.0
2 10 0.5
3
4
5 6 1.0
6 7 1.0 8 1.0 9 1.0
7
8
9
10 11 1.5
11

Title Card Required

0 1
C                 -0.52025999    0.82290527    0.93457395
H                  0.90364766    0.26313361    0.94110648
H                 -0.67879266    1.41448788    1.83206390
H                 -0.64922257    1.46912743    0.07064345
O                 -1.41071434   -0.29929100    0.88480490
C                 -2.79717060    0.05903855    0.87303759
H                 -3.04428636    0.66447364    0.00557829
H                 -3.35914020   -0.86170338    0.83541994
H                 -3.07326981    0.61014747    1.76754678
O                  1.92857101   -0.04416437    0.94083050
O                  2.85748178    0.77364850    0.84153154

1 3 1.0 4 1.0 5 1.0
2 10 1.0
3
4
5 6 1.0
6 7 1.0 8 1.0 9 1.0
7
8
9
10 11 2.0
11

不收敛的提示是link 9999. 返回小木虫查看更多

今日热帖
  • 精华评论
  • sky诺烯

    有时候用scf=qc可以收敛

  • 涵源亦可

    你的反应物做过结构优化吗?是放在一起优化的还是?氧气的位置是如何确定的呢

  • paramecium86

    9999 是opt没有收敛。 SCF=qc 还是对于 SCF不收敛比较有用。(L502 error)
    我看命令已经加了很多来帮助收敛了
    如果看能量是振荡了 所以不converge 也可以试试改变opt步长(maxstep=1,noTrust)
    或者还可以加上opt=GDIIS  虽然其实我觉得加上了用处可能也有限。

    不行的话可能还是得修改修改初猜结构 可能是太不准了。

    另外 用HF 不行 没准用DFT就能converge了 都可以试试。

  • chesterliu

    o2 is a triplet species

  • 明8221

    引用回帖:
    4楼: Originally posted by paramecium86 at 2017-01-14 05:14:00
    9999 是opt没有收敛。 SCF=qc 还是对于 SCF不收敛比较有用。(L502 error)
    我看命令已经加了很多来帮助收敛了
    如果看能量是振荡了 所以不converge 也可以试试改变opt步长(maxstep=1,noTrust)
    或者还可以加上 ...

    谢谢,问题已经解决

  • 明8221

    引用回帖:
    4楼: Originally posted by paramecium86 at 2017-01-14 05:14:00
    9999 是opt没有收敛。 SCF=qc 还是对于 SCF不收敛比较有用。(L502 error)
    我看命令已经加了很多来帮助收敛了
    如果看能量是振荡了 所以不converge 也可以试试改变opt步长(maxstep=1,noTrust)
    或者还可以加上 ...

    谢谢,问题已经解决。

  • 明8221

    引用回帖:
    5楼: Originally posted by chesterliu at 2017-01-14 06:19:41
    o2 is a triplet species

    你好,请问对于以O2作为RH分子H获取得到R基团和O2H基团,需要特殊考虑O2 is a triplet species吗,用B3LYP和M06-2x来进行过渡态优化可以吗?

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