计算化学公社

标题: TD opt计算出错 [打印本页]

作者
Author: zhongyuabc 时间: 2017-1-15 15:59
标题: TD opt计算出错
在计算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来解决？是否有一些其他的推荐办法？
附件中是完整的输入文件和输出文件

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

作者
Author: zhongyuabc 时间: 2017-1-15 17:17

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

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

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

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

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

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

作者
Author: zhongyuabc 时间: 2017-1-16 09:43

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

好的，谢谢您

作者
Author: zhongyuabc 时间: 2017-1-16 09:44

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

好的，谢谢您，我去试试

作者
Author: zhongyuabc 时间: 2017-1-16 19:37

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

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)