24小时热门版块排行榜    

查看: 2262  |  回复: 8

wzjxq

银虫 (初入文坛)

[交流] 【求助】过渡态搜索问题

gaussian新手,目前照着文献作C+NO——CO+N的路径分析。文献中的路径如下:

TS1,M2都得到了,TS2结构一直搜索不到。
TS方法输入文件:
%Chk=cnoTS_2.chk
#T b3lyp/6-31G(d) opt=(TS,noeigen,calcfc) Freq

C+NO-CO+1/2N2

0,2
N
O,1,1.37
C,2,1.30,1,124.0

上面的结构是没辙之后直接用的文献中结果,还是出错。
>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(UB+HF-LYP) =  -167.774383685     A.U. after  129 cycles
             Convg  =    0.1807D-03             -V/T =  2.0093
             S**2   =   0.8331
Annihilation of the first spin contaminant:
S**2 before annihilation     0.8331,   after     0.7504
Convergence failure -- run terminated.
Error termination via Lnk1e in C:\G03W\l502.exe at Sat Oct 09 10:05:29 2010.
Job cpu time:  0 days  0 hours  3 minutes 28.0 seconds.
File lengths (MBytes):  RWF=     11 Int=      0 D2E=      0 Chk=     12 Scr=      1
循环次数加到500,还是无法收敛。

是输入文档中有什么因素没有考虑吗?请大家解惑!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

manson1998

金虫 (著名写手)

OLED达人

引用回帖:
Originally posted by wzjxq at 2010-10-09 10:16:31:
gaussian新手,目前照着文献作C+NO——CO+N的路径分析。文献中的路径如下:

TS1,M2都得到了,TS2结构一直搜索不到。
TS方法输入文件:
%Chk=c ...

用qst2/3吧
惺惺相惜情不自禁
2楼2010-10-09 10:36:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

erylingjet

铁杆木虫 (著名写手)

小分子的过渡态应该对结构很敏感,呵呵
3楼2010-10-09 10:49:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

314202528

铜虫 (小有名气)

建议用qst2方法试一下
被高斯弄的焦头烂额中....
4楼2010-10-09 13:36:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wzjxq

银虫 (初入文坛)

qst2用过了
qst2
%Chk=cnoQST2.chk
#T b3lyp/6-31G(d) opt=(QST2,Noeigen,calcfc,maxcyc=500)  scf=(maxcyc=500) Freq

C+NO-CO+1/2N2

0,2
O
C,1,1.405
N,1,1.427,2,56.5

Product

0,2
O
C,1,1.185
N,1,1.346,2,180.0

出现l9999.exe错误,看结构已经偏离很远了。
Error termination request processed by link 9999.
Error termination via Lnk1e in C:\G03W\l9999.exe at Sat Oct 09 12:43:50 2010.
Job cpu time:  0 days  0 hours 14 minutes  9.0 seconds.
File lengths (MBytes):  RWF=     11 Int=      0 D2E=      0 Chk=      4 Scr=      1

想不通啊,明明是很简单,只有三原子的结构。。。。
5楼2010-10-09 14:15:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

冬天里的骄阳

铁杆木虫 (正式写手)

Convergence failure -- run terminated.
Error termination via Lnk1e in C:\G03W\l502.exe at Sat Oct 09 10:05:29 2010.


还是用TS法,关键词加上iop(5/13=1)避免上述错误
6楼2010-10-09 14:40:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

manson1998

金虫 (著名写手)

OLED达人

引用回帖:
Originally posted by wzjxq at 2010-10-09 14:15:14:
qst2用过了
qst2
%Chk=cnoQST2.chk
#T b3lyp/6-31G(d) opt=(QST2,Noeigen,calcfc,maxcyc=500)  scf=(maxcyc=500) Freq

C+NO-CO+1/2N2

0,2
O
C,1,1.405
N,1,1.427,2,56.5

Product

0,2
O
C, ...

还有scan嘛 三个原子计算量很小 可以用scan试试
惺惺相惜情不自禁
7楼2010-10-09 17:14:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wzjxq

银虫 (初入文坛)

引用回帖:
Originally posted by 冬天里的骄阳 at 2010-10-09 14:40:18:
Convergence failure -- run terminated.
Error termination via Lnk1e in C:\G03W\l502.exe at Sat Oct 09 10:05:29 2010.


还是用TS法,关键词加上iop(5/13=1)避免上述错误

得到了一个稳定结构  N1-O2-C3 有一个200多的虚频,
! R1    R(1,2)                  2.027          -DE/DX =    0.0                 !
! R2    R(2,3)                  1.1445         -DE/DX =    0.0                 !
! A1    A(1,2,3)              122.7163         -DE/DX =    0.0                 !

跟文献不一样,这完全就是CO和N了

作IRC出现以下错误,没找到一个点。。
Harris functional with IExCor=  402 diagonalized for initial guess.
ExpMin= 1.69D-01 ExpMax= 5.48D+03 ExpMxC= 8.25D+02 IAcc=1 IRadAn=         1 AccDes= 1.00D-06
HarFok:  IExCor= 402 AccDes= 1.00D-06 IRadAn=         1 IDoV=1
ScaDFX=  1.000000  1.000000  1.000000  1.000000
Integral accuracy reduced to 1.0D-05 until final iterations.
A guess orbital has undefined symmetry, so N**3 symmetry is turned off.
Initial convergence to 1.0D-05 achieved.  Increase integral accuracy.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Density matrix is not changing but DIIS error= 3.35D-06 CofLast= 4.02D-01.
The SCF is confused.
Error termination via Lnk1e in C:\G03W\l502.exe at Sat Oct 09 17:24:22 2010.
Job cpu time:  0 days  0 hours  7 minutes  7.0 seconds.
File lengths (MBytes):  RWF=     16 Int=      0 D2E=      0 Chk=     12 Scr=      1

[ Last edited by wzjxq on 2010-10-9 at 17:50 ]
8楼2010-10-09 17:43:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fatpig8832

铁杆木虫 (著名写手)

我做过阴离子的[C, O, N],楼主图中的ts2,用b3lyp,b3pw91,b98,bhandh时O-N键超过2.0,但用mp2和ccsd(t)来做时却在1.4~1.5左右,比较正常,估计用ccsd来做也一样。

我的猜测是:dft方法不行...
9楼2010-10-13 12:48:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wzjxq 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +9 Ldrop2023 2026-08-13 9/450 2026-08-16 13:24 by wukun_007
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +32 医学老男孩 2026-08-13 72/3600 2026-08-16 13:15 by 医学老男孩
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +10 majunge000 2026-08-11 12/600 2026-08-16 08:18 by xli1984
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] filecode=后面第一个是大写字母 +6 wangze12014 2026-08-14 7/350 2026-08-15 20:12 by gltch
[基金申请] filecode +8 cratir 2026-08-14 12/600 2026-08-15 18:08 by zyfgau
[基金申请] 关于Filecode分析方法 +10 majunge000 2026-08-10 13/650 2026-08-15 12:26 by hanpeng972
[硕博家园] 请教兼职经验 +3 是阿文鸭 2026-08-09 3/150 2026-08-14 12:07 by HER12025
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
信息提示
请填处理意见