24小时热门版块排行榜    

查看: 1000  |  回复: 7
当前主题已经存档。

雪狼790

木虫 (正式写手)

[交流] 【其他】关于不能收敛,再次请教

上次的PCM问题解决了,但又出现了新的问题。

在计算过渡态或负离子时,用了PCM方法,无论是CPCM还是IEFPCM,都总是不能收敛。我查了坛子里关于收敛的几个贴子,试着用了SCF=tight和nosymm,在我这儿都没效果。贼郁闷

敬请各位大侠指教。

拜谢了先。
回复此楼
人生,郁闷是常态,我只是在每天力求变态而已
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

poppyly

铁虫 (小有名气)


lei0736(金币+1,VIP+0):谢谢 7-13 14:50
你把你的输入文件和输入文件贴出来呀,我们可以帮你试试
2楼2009-07-13 10:42:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

雪狼790

木虫 (正式写手)

输入文件是:

%chk=anion1.chk
%mem=180MW
%nproc=1
# opt=(calcfc,tight,noeigentest) freq=noraman b3lyp/6-31g
scrf=(iefpcm,solvent=water) nosymm geom=connectivity

Title Card Required

-1 1
C                  2.46886100    2.79624500    0.35610300



输出文件的结尾部分:

Sum of Mulliken charges=  -1.00000
Electronic spatial extent (au):  =  5583.4784
Charge=    -1.0000 electrons
Dipole moment (field-independent basis, Debye):
    X=     5.4592    Y=    10.2715    Z=    -3.7174  Tot=    12.2117
Quadrupole moment (field-independent basis, Debye-Ang):
   XX=  -125.4455   YY=  -129.1495   ZZ=  -125.4620
   XY=   -10.9941   XZ=   -13.3242   YZ=     4.5798
Traceless Quadrupole moment (field-independent basis, Debye-Ang):
   XX=     1.2402   YY=    -2.4639   ZZ=     1.2237
   XY=   -10.9941   XZ=   -13.3242   YZ=     4.5798
Octapole moment (field-independent basis, Debye-Ang**2):
  XXX=    26.6982  YYY=    65.8760  ZZZ=   -22.7019  XYY=    64.2020
  XXY=    42.7981  XXZ=   -29.7102  XZZ=    19.1459  YZZ=    -6.2078
  YYZ=   -19.2994  XYZ=    31.8921
Hexadecapole moment (field-independent basis, Debye-Ang**3):
XXXX= -4631.5026 YYYY= -1968.2059 ZZZZ=  -721.6718 XXXY=   -99.5782
XXXZ=  -183.2045 YYYX=  -155.8201 YYYZ=    63.2964 ZZZX=   -49.7033
ZZZY=    46.9189 XXYY= -1128.1435 XXZZ=  -947.0148 YYZZ=  -412.0342
XXYZ=     6.6839 YYXZ=   -92.4755 ZZXY=   -25.0067
N-N= 1.452349779807D+03 E-N=-5.744468407309D+03  KE= 1.193312712187D+03


SCIENCE AT ITS BEST PROVIDES US WITH BETTER QUESTIONS,
NOT ABSOLUTE ANSWERS

                                    -- NORMAN COUSINS, 1976
Error termination request processed by link 9999.
Error termination via Lnk1e in C:\G03W\l9999.exe at Wed Jul 08 21:07:59 2009.
Job cpu time:  2 days 10 hours 58 minutes 41.0 seconds.
File lengths (MBytes):  RWF=    412 Int=      0 D2E=      0 Chk=     12 Scr=      1
人生,郁闷是常态,我只是在每天力求变态而已
3楼2009-07-13 11:26:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

heather889

木虫 (小有名气)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
lei0736(金币+2,VIP+0):谢谢 7-13 20:47
你可以增加循环次数
例如命令:
opt=(maxcycle=200) freq optcyc=300 scfcyc=900
此外,建议你检查初始结构是不是合理呀
4楼2009-07-13 17:46:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

雪狼790

木虫 (正式写手)

初始结构是不用PCM时优化的结构,按说应该没啥问题。

我用maxcycle试试吧。

多谢指点
人生,郁闷是常态,我只是在每天力求变态而已
5楼2009-07-13 17:50:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

3867826

木虫 (著名写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
lei0736(金币+1,VIP+0):谢谢 7-13 20:47
增加循环次数 或者降低收敛标准
拒稿,终于体会到了!乱写接受,认真写拒稿,这就是2和4的区别吗?
6楼2009-07-13 19:12:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

neweroica

木虫 (著名写手)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
gwdavid(金币+2,VIP+0):谢谢! 7-14 10:43
阴离子应该加弥散吧。

另外,用Guess=Huckel试试
Simpler, stronger, more tolerant
7楼2009-07-14 06:27:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

雪狼790

木虫 (正式写手)

多谢各位大侠
人生,郁闷是常态,我只是在每天力求变态而已
8楼2009-07-14 09:00:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 雪狼790 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见