24小时热门版块排行榜    

查看: 1490  |  回复: 9

明8221

新虫 (初入文坛)

[求助] 求助:寻找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.
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sky诺烯

新虫 (初入文坛)

2楼2017-01-07 16:45:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

涵源亦可

银虫 (小有名气)

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

发自小木虫Android客户端
3楼2017-01-13 14:50:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paramecium86

版主 (著名写手)

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

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

另外 用HF 不行 没准用DFT就能converge了 都可以试试。
4楼2017-01-14 05:14:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (著名写手)

本帖仅楼主可见

» 本帖已获得的红花(最新10朵)

5楼2017-01-14 06:19:41
已阅   申请QC强帖   回复此楼   编辑   查看我的主页

明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)
或者还可以加上 ...

谢谢,问题已经解决。
6楼2017-01-15 14:10:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

明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)
或者还可以加上 ...

谢谢,问题已经解决。
7楼2017-01-15 14:11:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

明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来进行过渡态优化可以吗?
8楼2017-01-15 14:15:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (著名写手)

本帖仅楼主可见
9楼2017-01-15 22:58:15
已阅   申请QC强帖   回复此楼   编辑   查看我的主页

明8221

新虫 (初入文坛)

引用回帖:
9楼: Originally posted by chesterliu at 2017-01-15 22:58:15
M062X is fine. Multiplicity is tricky but necessarily needed to be considered.

Thank you.
10楼2017-01-17 16:20:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 明8221 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 gy1nBQXYQJqL 2026-08-29 8/400 2026-08-30 12:14 by l0VvVHGBGRLv
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 11:41 by l0VvVHGBGRLv
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +6 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 11:19 by l0VvVHGBGRLv
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ASdOkHsho7FD 2026-08-28 6/300 2026-08-30 11:10 by l0VvVHGBGRLv
[基金申请] 投票:  有多少人是今天查系统知道结果的? +16 爱看书的可乐 2026-08-26 18/900 2026-08-30 10:33 by winsaint
[基金申请] 面上意见出来了 +9 黄鸟于飞Chao 2026-08-29 17/850 2026-08-30 10:28 by 孤独的英雄6
[基金申请] 国社科又开始会评了,不知道这次命运如何 +3 雨打竹帘 2026-08-30 3/150 2026-08-30 09:43 by KNIGHT2011
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +6 gdd2018 2026-08-28 11/550 2026-08-30 08:58 by 超级无敌华子
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 G6APbkg8SA6w 2026-08-29 4/200 2026-08-30 07:28 by ZPa0EcMwuECS
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 G6APbkg8SA6w 2026-08-29 4/200 2026-08-30 07:28 by ZPa0EcMwuECS
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 ASdOkHsho7FD 2026-08-28 13/650 2026-08-30 05:38 by ZPa0EcMwuECS
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 02:09 by ZPa0EcMwuECS
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 01:34 by ZPa0EcMwuECS
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 ASdOkHsho7FD 2026-08-28 5/250 2026-08-30 01:03 by ZPa0EcMwuECS
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
信息提示
请填处理意见