24小时热门版块排行榜    

查看: 1987  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 308求调剂 +3 是Lupa啊 2026-03-16 3/150 2026-03-16 10:07 by 求调剂zz
[考研] 本人考085602 化学工程 专硕 +9 不知道叫什么! 2026-03-15 10/500 2026-03-16 10:04 by houyaoxu
[考博] 欢迎申博同学联系 +3 天道酬勤2026686 2026-03-10 7/350 2026-03-15 19:03 by 天道酬勤2026686
[考研] 268求调剂 +5 一定有学上- 2026-03-14 6/300 2026-03-14 22:20 by 运气yunqi
[考研] 080500,材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
[考研] 308求调剂 +4 是Lupa啊 2026-03-09 4/200 2026-03-14 02:06 by tranquil_ya
[考研] 307求调剂 +7 超级伊昂大王 2026-03-10 7/350 2026-03-14 00:49 by JourneyLucky
[考研] 341求调剂 +3 番茄头--- 2026-03-10 3/150 2026-03-13 23:07 by JourneyLucky
[考研] 085600调剂 +5 漾漾123sun 2026-03-12 5/250 2026-03-13 22:06 by 星空星月
[考研] 0856材料与化工301求调剂 +5 奕束光 2026-03-13 5/250 2026-03-13 22:00 by 星空星月
[考研] 四川大学085601材料工程专硕 初试294求调剂 +4 祝我们好在冬天 2026-03-11 4/200 2026-03-13 21:39 by peike
[考研] 301求调剂 +6 Liyouyumairs 2026-03-11 6/300 2026-03-13 20:11 by JourneyLucky
[考研] 材料工程调剂 +4 咪咪空空 2026-03-11 4/200 2026-03-13 19:57 by JourneyLucky
[考研] 274求调剂 +3 S.H1 2026-03-12 3/150 2026-03-13 15:15 by JourneyLucky
[考研] 求调剂 +3 程雨杭 2026-03-12 3/150 2026-03-13 15:06 by JourneyLucky
[考研] 081200-11408-276学硕求调剂 +3 崔wj 2026-03-12 4/200 2026-03-12 19:33 by 求调剂zz
[考博] 2026年博士申请 +3 QwQwQW10 2026-03-11 3/150 2026-03-12 17:58 by gxch43
[考研] 341求调剂 +4 捣蛋猪猪 2026-03-11 4/200 2026-03-12 14:47 by ruiyingmiao
[考博] 读博申请 +5 感dd 2026-03-10 7/350 2026-03-11 17:02 by QGZDSYS
[考研] 收调剂 +7 调剂的考研学生 2026-03-10 7/350 2026-03-10 17:57 by 麦茶汤圆
信息提示
请填处理意见