求助stable关键词与激发能过小

zhanfei

前面提要如题进行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吗？
谢谢大家

sobereva

1 如果计算无误的话，应该不影响。但建议也试试其它泛函如cam-b3lyp、wb97X，看看是否都这样

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

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

zhanfei

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

2 使用姿势 ...

分子带电荷没问题，弱弱的问一下带电分子，两个未成对电子会导致什么吗？

sobereva

也不会导致什么，只不过大多体系基态是单重态

高雪皎

我也算激发能，算出来第一激发能竟然是负的，我用了stable=Opt 检查了波函数稳定性，结果是稳定的：The wavefunction is stable under the perturbations considered.  不知道原因，现在正在试着调整分子多重度，不知道可不可行。