24小时热门版块排行榜    

查看: 1986  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 289求调剂 +5 步川酷紫123 2026-03-11 5/250 2026-03-15 00:45 by kruisytel
[考研] 本科南京大学一志愿川大药学327 +3 麦田耕者 2026-03-14 3/150 2026-03-14 20:04 by 外星文明
[基金申请] 现在如何回避去年的某一个专家,不知道名字 +3 zk200107 2026-03-12 6/300 2026-03-14 17:13 by zk200107
[考研] 267一志愿南京工业大学0817化工求调剂 +5 SUICHILD 2026-03-12 5/250 2026-03-14 14:53 by jean5056
[考研] 求调剂 +3 清风问长安 2026-03-09 3/150 2026-03-14 02:15 by JourneyLucky
[考研] 2026考研调剂+本科延边大学+山东大学+生物化学与分子生物学+有项目经验 +3 ccdsscjy 2026-03-09 6/300 2026-03-14 02:14 by JourneyLucky
[考研] 308求调剂 +4 是Lupa啊 2026-03-09 4/200 2026-03-14 02:06 by tranquil_ya
[考研] 085600求调剂 +3 a邵星池 2026-03-09 3/150 2026-03-14 01:32 by JourneyLucky
[考研] 一志愿湖师大化学289求调剂 +6 XMCMM3.14159 2026-03-10 6/300 2026-03-14 00:28 by JourneyLucky
[考研] 0703化学调剂 +4 快乐的香蕉 2026-03-11 4/200 2026-03-13 22:41 by JourneyLucky
[考研] 085600调剂 +5 漾漾123sun 2026-03-12 5/250 2026-03-13 22:06 by 星空星月
[考研] 26调剂/材料科学与工程/总分295/求收留 +9 2026调剂侠 2026-03-12 9/450 2026-03-13 20:46 by 18595523086
[考研] 材料与化工085600调剂求老师收留 +9 jiaanl 2026-03-11 9/450 2026-03-13 20:22 by JourneyLucky
[考研] 【0856】化学工程(085602)313 分,本科学科评估A类院校化学工程与工艺,诚求调剂 +7 小刘快快上岸 2026-03-11 7/350 2026-03-13 16:06 by ruiyingmiao
[考研] 308求调剂 +3 是Lupa啊 2026-03-12 3/150 2026-03-13 14:30 by 求调剂zz
[考研] 0703一志愿211 285分求调剂 +4 ly3471z 2026-03-13 4/200 2026-03-13 13:00 by JourneyLucky
[考博] 26读博 +4 Rui135246 2026-03-12 10/500 2026-03-13 07:15 by gaobiao
[考研] 081200-11408-276学硕求调剂 +3 崔wj 2026-03-12 4/200 2026-03-12 19:33 by 求调剂zz
[考研] 0857环境调剂 +5 熠熠_11 2026-03-10 5/250 2026-03-11 10:59 by wang_dand
[考研] 收调剂 +7 调剂的考研学生 2026-03-10 7/350 2026-03-10 17:57 by 麦茶汤圆
信息提示
请填处理意见