TD opt计算出错

zhongyuabc

在计算tdopt时，输入文件如下：%chk=tdopt320.chk
%mem=8GB
%nprocshared=16
#M06/genecp TD(NSTATES=20,root=1) opt

Title Card Required

0 1
In               0    0.17307600   -0.05858500    0.07312900
O                0    2.23583000   -0.80600100    0.89519900
O                0    1.23702700   -1.14471100   -1.52027900
O                0    4.16164600   -1.76187600    0.24128100
O                0    3.14068500   -2.14271800   -2.16754100

.....
In              -1    5.09018900   -2.69396200   -1.44218000
O                0    5.28251600   -4.52300900   -1.06407000
H                0    5.94303700   -5.05926700   -1.49497900
In              -1   -4.16570900    3.21541500   -1.75448600
O                0   -2.18634700    2.55795000   -2.27541900
O                0   -3.51978900    2.18710400   -0.01856500
O                0   -0.48090900    1.34179000   -1.47872900
O                0   -1.80650200    0.97431300    0.77741000
C                0   -1.61142400    1.85313400   -1.40958200
C                0   -2.38554700    1.63885100   -0.08590000
In              -1   -4.10758200   -3.54892200    1.47652700


O N C H 0
6-311G(d,p)
****
In 0
LanL2DZ
****

In 0
LanL2DZ


输出文件的末尾显示如下
No map to state      1
You need to solve for more vectors in order to follow this state.
Error termination via Lnk1e in /home/g09/g09/l914.exe at Sat Jan 14 14:16:41 2017.


我看了之前其他人遇到的问题，公社论坛上回复加大Nstates试试，不过我还是有几个问题不明白，想请教一下各位前辈们：
1.我研究的是root=1，按理说Nstates=20已经是可以了，难道还需要加大Nstates？
2.我看了输出文件中最后部分荧光发射波长的计算，荧光发射波长竟为5310.95 nm，这是否可以知道homo与lumo之间能级间隔很小？
Excitation energies and oscillator strengths:

Excited State   1:      Singlet-A      0.2335 eV 5310.95 nm  f=0.0007  <S**2>=0.000
     141 ->142         0.70686
This state for optimization and/or second-order correction.
Copying the excited state density for this state as the 1-particle RhoCI density.

Excited State   2:      Singlet-A      2.5742 eV  481.65 nm  f=0.0000  <S**2>=0.000
     141 ->143        -0.69927

Excited State   3:      Singlet-A      2.7243 eV  455.11 nm  f=0.0002  <S**2>=0.000
     137 ->142        -0.18126
     138 ->142        -0.56218
     139 ->142         0.37374



3.在输出文件的前半部分，出现的是第一激发态的能量比第二激发态能量高，这是因为出现了能级交错吗？
Excitation energies and oscillator strengths:

Excited State   1:      Singlet-A      2.8343 eV  437.44 nm  f=0.0009  <S**2>=0.000
     138 ->142        -0.47790
     139 ->142         0.34560
     140 ->142         0.36395
This state for optimization and/or second-order correction.
Copying the excited state density for this state as the 1-particle RhoCI density.

Excited State   2:      Singlet-A      2.7638 eV  448.60 nm  f=0.0000  <S**2>=0.000
     141 ->142        -0.70215

Excited State   3:      Singlet-A      3.0116 eV  411.69 nm  f=0.0001  <S**2>=0.000
     138 ->142         0.35680
     139 ->142         0.59437

4.出现 这样情况的体系是否不能单靠增大NStates来解决？是否有一些其他的推荐办法？
附件中是完整的输入文件和输出文件

小范范1989

第一个问题，当时存了一个pdf，给你看看
第二个问题，可能是存在了势能面交叉。我记得论坛了问过好多次了。

zhongyuabc

小范范1989 发表于 2017-1-15 16:18
第一个问题，当时存了一个pdf，给你看看
第二个问题，可能是存在了势能面交叉。我记得论坛了问过好多次了 ...

谢谢您，我之前看过了高斯官方关于no map to state **的回答的，只是还是有以上四个问题不太明白。非常感谢您的解答

ggdh

参考这个帖子http://bbs.keinsci.com/forum.php?mod=viewthread&tid=4905
尝试以下解决方案：
1。 换gaussian09 D01 或者E01版
或者
2。把td里面的root=1 去掉。

sobereva

1、4 优化S1，当前都已经nstates=20了，继续增大不会有用。原理和bug在ggdh提到的帖子里已有讨论
另外，root=1本身就是默认的，写了完全多余

2 你实际看一下HOMO-LUMO gap是多少就完了。这和没有必然关系，用一些HF成分不太高的泛函算一些大共轭体系容易出现ghost态，就是激发能非常小，f也非常小。你可以试试BHandHLYP、M06-2X这些高HF成份的泛函或wB97X，如果不再有这个态说明这个太是artifical的。

3 是因为优化过程中能级顺序发生了变化所致

zhongyuabc

ggdh 发表于 2017-1-16 00:01
参考这个帖子http://bbs.keinsci.com/forum.php?mod=viewthread&tid=4905
尝试以下解决方案：
1。 换gaus ...

好的，谢谢您

zhongyuabc

sobereva 发表于 2017-1-16 06:46
1、4 优化S1，当前都已经nstates=20了，继续增大不会有用。原理和bug在ggdh提到的帖子里已有讨论
另外，ro ...

好的，谢谢您，我去试试

zhongyuabc

sobereva 发表于 2017-1-16 06:46
1、4 优化S1，当前都已经nstates=20了，继续增大不会有用。原理和bug在ggdh提到的帖子里已有讨论
另外，ro ...

老师，我用了M06-2X/3-21G优化激发态结构，能够优化成功，但是就是光谱峰相比实验值是蓝移，我现在准备用BhandHLYP/6-31G*试试看。