24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 2009  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 专硕304找调剂,一线城市最好 +3 李lsl李 2026-04-05 3/150 2026-04-06 12:16 by ffffjjjj
[考研] 专硕0854初试考材科基,求调剂 +6 3220548044 2026-04-06 9/450 2026-04-06 10:26 by barlinike
[考研] 0817化学工程与技术求调剂,一志愿中海洋319 +14 lv945 2026-04-04 14/700 2026-04-06 10:20 by 蓝云思雨
[考研] 332求调剂 +17 小小孟... 2026-04-05 18/900 2026-04-06 09:51 by 蓝云思雨
[考研] 085600,320分求调剂 +16 大馋小子 2026-04-04 17/850 2026-04-06 07:58 by MOF_Catal
[考研] 生物与医药求调剂 +5 heguanhua 2026-04-05 6/300 2026-04-05 22:58 by Hdyxbekcb
[考研] 322求调剂 +3 嗯哼哼恒 2026-04-05 3/150 2026-04-05 19:52 by nepu_uu
[考研] 304求调剂(085602,过四级,一志愿985) +15 化工人999 2026-04-04 15/750 2026-04-05 16:25 by 我是电风扇r
[考研] 一志愿上海海洋大学083200食品学硕,求调剂,接受其他专业083200 +4 what张 2026-04-04 5/250 2026-04-05 14:07 by chw1980_0
[考研] 353求调剂 +10 MayUxw1 2026-04-03 10/500 2026-04-05 09:23 by 无际的草原
[考研] 083200 333求调剂 +3 十二!! 2026-04-04 3/150 2026-04-05 08:28 by barlinike
[考研] 344材料与化工调剂 +9 调剂上岸玘 2026-04-03 9/450 2026-04-04 23:10 by happyddm
[考研] 一志愿中国石油大学化学工程323分求调剂 +4 化工专硕323分 2026-04-03 6/300 2026-04-03 22:12 by dongzh2009
[考研] 286求调剂 +8 lim0922 2026-04-02 8/400 2026-04-03 20:19 by rzh123456
[考研] 材料科学与工程339求调剂 +12 hyz0119 2026-03-31 13/650 2026-04-03 18:33 by ls刘帅
[考研] 301求调剂 +14 A_JiXing 2026-04-01 14/700 2026-04-03 18:31 by ls刘帅
[考研] 08工科275分求调剂 +14 AaAa7420 2026-03-31 14/700 2026-04-03 11:13 by cocolv
[考研] 生物学327,求调剂 +5 书上的梅子 2026-04-01 6/300 2026-04-02 06:47 by ilovexiaobin
[考研] 310分求调剂 +4 成功上岸wang 2026-04-01 4/200 2026-04-01 20:35 by liu823948201
[考研] 267求调剂 +13 uiybh 2026-03-31 13/650 2026-04-01 10:25 by 探123
信息提示
请填处理意见