24小时热门版块排行榜    

查看: 1923  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见