24小时热门版块排行榜    

查看: 1120  |  回复: 11
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

kite2005

银虫 (小有名气)

[交流] 【求助】用TD做激发能似乎进入死循环了

命令行如下:
#p PBE1PBE/6-31+G(d) TD(singlet, NStates=6) scfcon=7 test
正常的情况下,类似分子只需约十几个小时就可计算完,但现在一个分子算了十几天了还处在不断的循环中,不知是什么原因,现将部分循环内容贴至下面,还望大侠们多多指点
Iteration    11 Dimension   144 NMult  2640
Cannot handle 2e integral symmetry, ISym2E=1.
CISAX:  IP=  1 NPass=   1 NMax=  12.
CISAX will form  12 AO SS matrices at one time.
NMat= 12 NSing= 12.
PrismC:  NFx=      2048 NFxT=        12 NFxU=        12.
Cannot handle 2e integral symmetry, ISym2E=1.
CISAX:  IP=  1 NPass=   1 NMax=  12.
CISAX will form  12 AO SS matrices at one time.
NMat= 12 NSing= 12.
PrismC:  NFx=      2048 NFxT=        12 NFxU=        12.
New state      1 was old state      3
Root      1 not converged, maximum delta is    0.516669365043637
New state      2 was old state      1
Root      2 not converged, maximum delta is    0.197406735008100
New state      3 was old state      6
Root      3 not converged, maximum delta is    0.136594616961202
New state      4 was old state      3
Root      4 not converged, maximum delta is    0.113178240838583
New state      5 was old state      4
Root      5 not converged, maximum delta is    0.131739810272319
New state      6 was old state      5
Root      6 not converged, maximum delta is    0.183663207407156
Excitation Energies [eV] at current iteration:
Root      1 :     1.520261072481056   Change is   -1.217243276404173
Root      2 :     2.360011931923294   Change is    0.435827506597912
Root      3 :     2.677578443751468   Change is   -0.826738085163761
Root      4 :     2.829977037376423   Change is    0.092472688491195
Root      5 :     2.950892467024729   Change is    0.037604877077249
Root      6 :     3.230262796012697   Change is    0.015591024266489
Iteration    12 Dimension   156 NMult  2652
Cannot handle 2e integral symmetry, ISym2E=1.
CISAX:  IP=  1 NPass=   1 NMax=  12.
CISAX will form  12 AO SS matrices at one time.
NMat= 12 NSing= 12.
PrismC:  NFx=      2048 NFxT=        12 NFxU=        12.
Cannot handle 2e integral symmetry, ISym2E=1.
CISAX:  IP=  1 NPass=   1 NMax=  12.
CISAX will form  12 AO SS matrices at one time.
NMat= 12 NSing= 12.
PrismC:  NFx=      2048 NFxT=        12 NFxU=        12.
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yuhuobuku

木虫 (正式写手)

★ ★ ★ ★ ★
xuefei06(金币+1,VIP+0):thanks! 3-3 23:00
kite2005(金币+2,VIP+0):谢谢参与 3-6 08:45
kite2005(金币+2,VIP+0):感谢参与 4-1 09:51
引用回帖:
Originally posted by loovfnd at 2009-3-3 17:28:
scfcon=7 这个关键词是什么意思,我算这个没用过这个scfcon=7。

收敛标准。scf关键字,conver=n选项。数字n代表10的-n次幂
薄酒可以忘忧,丑妻可以白头,徐行不必驷马,称身不必狐裘平淡才是生活的主旋律
5楼2009-03-03 19:01:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答

loovfnd

至尊木虫 (著名写手)

★ ★
kite2005(金币+2,VIP+0):感谢参与 4-1 09:51
scfcon=7 这个关键词是什么意思,我算这个没用过这个scfcon=7。
莫强求!
2楼2009-03-03 17:28:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

S07111072

铁杆木虫 (著名写手)

★ ★ ★ ★ ★
yuhuobuku(金币+1):欢迎参加讨论
kite2005(金币+2,VIP+0):谢谢参与 3-6 08:45
kite2005(金币+2,VIP+0):感谢参与 4-1 09:51
scfcon=7指的是收敛步数, 现在scfcon=7不收敛,可以试试scfcon=6
3楼2009-03-03 18:04:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

★ ★ ★ ★
yuhuobuku(金币+1):欢迎参加讨论
kite2005(金币+3,VIP+0):我是菜鸟,我理解这个TD是做垂直激发能,跟结构有关系么??未见输出文件中有构型变化的内容,请问是否能更详细说明 3-4 08:47
看一下是否是在某一个结构附件振荡……
4楼2009-03-03 18:06:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见