求助:寻找DME与O2反应的过渡态
求助:寻找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. 返回小木虫查看更多
有时候用scf=qc可以收敛
你的反应物做过结构优化吗?是放在一起优化的还是?氧气的位置是如何确定的呢
9999 是opt没有收敛。 SCF=qc 还是对于 SCF不收敛比较有用。(L502 error)
我看命令已经加了很多来帮助收敛了
如果看能量是振荡了 所以不converge 也可以试试改变opt步长(maxstep=1,noTrust)
或者还可以加上opt=GDIIS 虽然其实我觉得加上了用处可能也有限。
不行的话可能还是得修改修改初猜结构 可能是太不准了。
另外 用HF 不行 没准用DFT就能converge了 都可以试试。
o2 is a triplet species
谢谢,问题已经解决
,
谢谢,问题已经解决。
你好,请问对于以O2作为RH分子H获取得到R基团和O2H基团,需要特殊考虑O2 is a triplet species吗,用B3LYP和M06-2x来进行过渡态优化可以吗?