应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】用TD做激发能似乎进入死循环了
121/2返回列表12下一页
查看: 1115  |  回复: 11
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

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.
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)
1楼 2009-03-03 09:29:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

loovfnd

至尊木虫 (著名写手)
★ ★
kite2005(金币+2,VIP+0):感谢参与 4-1 09:51
scfcon=7 这个关键词是什么意思，我算这个没用过这个scfcon=7。
赞一下(2人) 回复此楼
莫强求！
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
赞一下(5人) 回复此楼
3楼2009-03-03 18:04:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)
★ ★ ★ ★
yuhuobuku(金币+1):欢迎参加讨论
kite2005(金币+3,VIP+0):我是菜鸟，我理解这个TD是做垂直激发能，跟结构有关系么？？未见输出文件中有构型变化的内容，请问是否能更详细说明 3-4 08:47
看一下是否是在某一个结构附件振荡……
赞一下(5人) 回复此楼
4楼2009-03-03 18:06:23
已阅   回复此楼   关注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次幂
赞一下(14人) 回复此楼
薄酒可以忘忧，丑妻可以白头，徐行不必驷马，称身不必狐裘平淡才是生活的主旋律
5楼2009-03-03 19:01:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kite2005

银虫 (小有名气)

回禀各位大侠

以前是此一步骤不收敛，加scfcon＝7后收敛了，未知我上述的死循环跟scfcon有关么？？
Cycle  25  Pass 1  IDiag  1:
RMSU=  3.90D-06    CP:  1.00D+00  7.73D-01  1.27D+00  1.02D+00  9.26D-01
PrismC:  NFx=      2048 NFxT=        12 NFxU=        12.
E= -1604.57320979898     Delta-E=       -0.000000000407 Rises=F Damp=F
DIIS: error= 1.18D-07 at cycle  25 NSaved=  20.
NSaved=20 IEnMin=16 EnMin= -1604.57320980467     IErMin=14 ErrMin= 1.17D-07
ErrMax= 1.18D-07 EMaxC= 1.00D-01 BMatC= 2.19D-12 BMatP= 1.83D-12
IDIUse=1 WtCom= 1.00D+00 WtEn= 0.00D+00
Coeff-Com: -0.269D-05-0.161D-04-0.188D-04-0.269D-03-0.268D-03-0.396D-03
Coeff-Com:  0.249D-02 0.129D-01 0.314D-01 0.471D-01 0.394D+00-0.422D+00
Coeff-Com:  0.282D+00 0.131D+00-0.315D+00 0.393D+00 0.464D+00-0.155D+00
Coeff-Com:  0.655D-01 0.707D-01
Coeff:     -0.269D-05-0.161D-04-0.188D-04-0.269D-03-0.268D-03-0.396D-03
Coeff:      0.249D-02 0.129D-01 0.314D-01 0.471D-01 0.394D+00-0.422D+00
Coeff:      0.282D+00 0.131D+00-0.315D+00 0.393D+00 0.464D+00-0.155D+00
Coeff:      0.655D-01 0.707D-01
Gap=     0.117 Goal=   None    Shift=    0.000
RMSDP=5.45D-08 MaxDP=2.78D-06 DE=-4.07D-10 OVMax= 5.03D-09

SCF Done:  E(RPBE+HF-PBE) =  -1604.57320980     A.U. after   25 cycles
             Convg  =    0.5448D-07             -V/T =  2.0093
             S**2   =   0.0000
KE= 1.589755105497D+03 PE=-1.030287435286D+04 EE= 3.823560461715D+03
Leave Link  502 at Mon Feb 16 17:26:03 2009, MaxMem=  196608000 cpu:   13445.0
(Enter C:\G03W\l801.exe)
ExpMin= 4.38D-02 ExpMax= 5.48D+03 ExpMxC= 8.25D+02 IAcc=3 IRadAn=         5 AccDes= 0.00D+00
HarFok:  IExCor= 205 AccDes= 0.00D+00 IRadAn=         5 IDoV=1
ScaDFX=  1.000000  1.000000  1.000000  1.000000
Largest valence mixing into a core orbital is  6.63D-05
Largest core mixing into a valence orbital is  1.16D-05
Range of M.O.s used for correlation:    37   726
NBasis=   736 NAE=   127 NBE=   127 NFC=    36 NFV=     0
NROrb=    690 NOA=    91 NOB=    91 NVA=   599 NVB=   599
赞一下 回复此楼
6楼2009-03-04 08:45:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kite2005

银虫 (小有名气)

re

自己顶一下，希望大侠们多多参与
赞一下 回复此楼
7楼2009-03-06 08:46:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kite2005

银虫 (小有名气)

大虾们都在哪里

求助，哈哈，自顶
回复此楼
8楼2009-03-08 17:27:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bigcontinent

金虫 (正式写手)
★ ★ ★ ★ ★ ★ ★ ★
yuhuobuku(金币+1,VIP+0):欢迎参加讨论 3-9 13:07
kite2005(金币+5,VIP+0):谢谢参与 3-24 14:21
kite2005(金币+2,VIP+0):感谢 4-1 09:52
应该是激发态的“自洽”没有收敛吧，不知道“自洽”这个词是否准确，我以前也遇到过类似的问题不过当时是我的态选的过多，你的问题可能出在初始构型，换一个构型再试试。
赞一下(9人) 回复此楼
9楼2009-03-09 13:02:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hjy_lxq

金虫 (小有名气)
P我觉得可以改一下你的构型，同时改变收敛步数！！！
赞一下 回复此楼
10楼2009-11-05 11:58:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 kite2005 的主题更新
121/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭