计算化学公社

标题: 求助stable关键词与激发能过小 [打印本页]

作者
Author:
zhanfei    时间: 2014-10-30 13:39
标题: 求助stable关键词与激发能过小
前面提要如题进行TD计算时激发能过小
Excited State   1:  3.013-A      0.2810 eV 4411.54 nm  f=0.0000  <S**2>=2.020
     76B -> 79B       -0.25932
     78B -> 79B        0.96322
This state for optimization and/or second-order correction.
Total Energy, E(TD-HF/TD-KS) =  -976.549782768   
Copying the excited state density for this state as the 1-particle RhoCI density.

Excited State   2:  3.014-A      0.5147 eV 2408.76 nm  f=0.0000  <S**2>=2.021
     75B -> 79B        0.99794


sob sennsei:最好先确认结构无误,并用stabe检测一下基态波函数稳定性,然后如果没问题再从物理意义上讨论。gap很小主要是过渡金属d轨道能级接近简并造成的,这也是这类体系电子静态相关效应往往较强的原因。

结构应该是没问题,附件是stable的输出文件。
问题是
1此时即使前几个激发态激发能过小,但是我关注的是400-500nm的MLCT,这部分计算结果是否也会受影响?
2 stable 关键词 的输出怎么理解?
我看输出文件中最后不断重复(Enter /workfs/cac/dwang/test_g09c/g09-c01/g09/l914.exe)进入l914开始激发态计算的SCF稳定性检测
UHF ground state
Ignoring spin specification with unrestricted ground state.
MDV=  1610612736 DFT=T DoStab=F Mixed=T DoRPA=T DoScal=F NonHer=T
到Convergence achieved on expansion vectors.
是看这个Convergence on wavefunction:    0.001000000000000吗?
谢谢大家


作者
Author:
sobereva    时间: 2014-10-30 16:55
1 如果计算无误的话,应该不影响。但建议也试试其它泛函如cam-b3lyp、wb97X,看看是否都这样

2 使用姿势不对。
你就写b3lyp/genecp stable就行了,其它的都不写,这样先看看基态的波函数是否稳定。

PS:确信这个结构的基态是三重态么?并且确信这是带两个正电荷么?

作者
Author:
zhanfei    时间: 2014-10-30 23:17
sobereva 发表于 2014-10-30 16:55
1 如果计算无误的话,应该不影响。但建议也试试其它泛函如cam-b3lyp、wb97X,看看是否都这样

2 使用姿势 ...

分子带电荷没问题,弱弱的问一下带电分子,两个未成对电子会导致什么吗?
作者
Author:
sobereva    时间: 2014-10-30 23:22
也不会导致什么,只不过大多体系基态是单重态
作者
Author:
高雪皎    时间: 2014-10-31 15:29
我也算激发能,算出来第一激发能竟然是负的,我用了stable=Opt 检查了波函数稳定性,结果是稳定的:The wavefunction is stable under the perturbations considered.  不知道原因,现在正在试着调整分子多重度,不知道可不可行。




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3