高斯优化激发态出错

小范范1989

flyskyzff 发表于 2014-12-11 09:37
请问优化t1,你这里说的“直接0 3计算就行了，不需要TDDFT”的意思是不是route部分只要# b3lyp/6-31* opt ...

对，也就是使用DFt的方式，0 3
http://emuch.net/bbs/viewthread.php?tid=8072030
   

小范范1989

sobereva 发表于 2014-12-7 15:43
又是三重态又是root=1，直接用0 3计算就行了，不需要TDDFT，又快又不会出现这些问题

sob好，之前发的帖子，是优化激发态的，按照你的方法，我优化出来了最低三重态，使用的0  3这样的方法，
但是，我有个问题，就是我优化出来后，然后td（50-50），结果如下，我不知道那个是激发态的三重态，哪个是单重态。谢谢sob
同时，我想的方法是，要是这样不能解决，我们可以在做完td（50-50）之后，在重新做td。这样，两个不同的，就是t1激发态s1或者是s2 s3....的。行得通吧？
谢谢sob

sobereva

你看<S**2>就能判断，见http://bbs.keinsci.com/forum.php?mod=viewthread&tid=111
貌似你是在0 3下做的TDDFT，因此造成<S**2>不为整数。不建议这样，TDDFT应当在0 1下做。

小范范1989

sobereva 发表于 2014-12-16 17:25
你看就能判断，见http://bbs.keinsci.com/forum.php?mod=viewthread&tid=111
貌似你是在0 3下做的TDDFT， ...

谢谢sob。我还是有点迷惑，我是在 03的情况下，优化出来了t1态，而不是使用0 1的方式优化出来t1，因为用0 1的方式，优化一直出错，而且也不准确，所以是在0 3的状态下优化出来的t1。
现在我在 03 优化出来的基础上，做的td（50-50），出现上面的结果，
按照您刚刚说的，TDDFT应该在 01下面做，那我应该怎么进行下一步？
是不是在 03 优化出来的基础上，然后采用0 3优化出来的坐标，然后 设置电荷为0，自旋为1，在此基础上做td？
谢谢sob，

sobereva

是的

小范范1989

sobereva 发表于 2014-12-16 20:03
是的

谢谢sob。我发现只要有不会的问题，来问sob，一定会得到解答。
sob不但知识面很广，而且为人和热情细心，谢谢sob。

小范范1989

sobereva 发表于 2014-12-16 20:03
是的

sob好，还是这个优化激发态的问题，我这不是也发过帖子，问过好多，也修改了好多，但是还是出问题，还请sob再次指点。
我的输入文件是：
%chk=test.chk
%mem=2GB
%nprocs=2
#p b3lyp/6-31g* opt td(nstates=10,singlets,root=1,direct)

test

0 1
N                  1.55801000    0.19783500   -0.79778700
C                  1.09085500   -0.95453600   -1.48926800
C                  2.14146600   -1.89758300   -1.54718200
C                 -0.13547400   -1.20956500   -2.10922100
C                  1.96005300   -3.11900500   -2.20189300
C                 -0.29383200   -2.43297000   -2.76144400
H                 -0.93851000   -0.48827500   -2.08290200
C                  0.73480500   -3.38468600   -2.80530000
H                  2.76816500   -3.84409000   -2.24412800
.....................（还有很多坐标）

错误时：
NMat=    20 NSing=    20.
New state      1 was old state      7
New state      2 was old state      8
New state      3 was old state      1
New state      4 was old state      6
New state      5 was old state      4
New state      6 was old state      3
New state      7 was old state      9
New state      8 was old state      5
New state      9 was old state     10
New state     10 was old state      2
Excitation Energies [eV] at current iteration:
Root      1 :     5.982060062346076
Root      2 :     6.383931700600978
Root      3 :     8.901901716469364
Root      4 :     9.700349738000840
Root      5 :    11.258633326878070
Root      6 :    11.308304626971730
Root      7 :    11.713305846440200
Root      8 :    12.248164103617310
Root      9 :    13.290151060085080
Root     10 :    13.512389675544020
Iteration     2 Dimension    40 NMult    40
Cannot handle 2e integral symmetry, ISym2E=1.
CISAX:  IP=  2 NPass=   1 NMax=    20.
CISAX will form    20 AO SS matrices at one time.
NMat=    20 NSing=    20.
Cannot handle 2e integral symmetry, ISym2E=1.
CISAX:  IP=  2 NPass=   1 NMax=    20.
CISAX will form    20 AO SS matrices at one time.
NMat=    20 NSing=    20.
New state      1 was old state      2
Root      1 not converged, maximum delta is    0.433400373447380
New state      2 was old state      1
Root      2 not converged, maximum delta is    0.528876387304794
No map to state      3
You need to solve for more vectors in order to follow this state.
Error termination via Lnk1e in /home/software/g09/l914.exe at Thu Dec 21 13:17:32 2006.
Job cpu time:  3 days  8 hours 11 minutes 34.4 seconds.
File lengths (MBytes):  RWF=   3198 Int=      0 D2E=      0 Chk=    442 Scr=      1

当时按照您说的，这样的错误 You need to solve for more vectors in order to follow this state，就是一个增加nstates，但是我都到20了，还在增加吗？而且，我看还有错误是不是不收敛？
Root      1 not converged, maximum delta is    0.433400373447380。
你说，是不是我应该重新调节一下我的分子构型？（但是，我的s1的初始构型，是优化的s0的构型）。或者是使用什么方法使其能计算出来？
谢谢sob，这个优化激发态，我已经弄了快2个月了，还是弄不出来， 而且我还没有计算激发态的频率。这样耗时间太长了。所以一次一次的来麻烦sob。望谅解。谢谢sob

sobereva

nstates再加到40还不行就不是这个原因了。
可以换换泛函、基组再试。
也可以先用比如CIS优化，如果能顺利收敛，然后再切换回TDDFT继续优化。
也可以减小步长上限，或者用别的几何优化方法试试。
总之没有万能办法，各种路子都试试。

小范范1989

sobereva 发表于 2014-12-18 18:56
nstates再加到40还不行就不是这个原因了。
可以换换泛函、基组再试。
也可以先用比如CIS优化，如果能顺利 ...

嗯，好的，晓得了，谢谢sob:lol

小范范1989

回复错了sob，在我的另外一个帖子，sorry。看着好几个您的回帖，一下弄乱了，sorry

yopoxia

您好，最近我也在学习高斯。想请教您几个问题。
我看了一些帖子，都说步骤如下，您看我理解的对不
1、先在高斯里画出结构，然后比如用这个计算“#p opt freq B3LYP/6-31G”，得到优化后的原子坐标。
2、再用第一步优化后的坐标，写计算命令“#p B3LYP/6-31G td=(nstates=10)”
后面的我就不懂了。。。。还有root什么的，s1  t1就晕了。。。。第一激发态的结构是怎么得到的。。。
您能帮我解答一下么。。。具体的步骤是什么。。。

小范范1989

yopoxia 发表于 2014-12-25 13:11
您好，最近我也在学习高斯。想请教您几个问题。
我看了一些帖子，都说步骤如下，您看我理解的对不
1、先 ...

你好，很高兴和你讨论问题。
1：这个root=n，我的感觉，就是对那个态感兴趣，比如说，我要优化第一单重激发态，那么就要写（singlets，root=1）。如果说对第二个三重态感兴趣，（triplets，root=2）
2：s1代表第一单重激发态，是分子的一种状态，和基态差不多意思。同样t1代表第一三重激发态。
3：做法：#p  b3lyp/6-31g opt  这样优化出来基态坐标。  然后 使用基态坐标，该表头为#p b3lyp/6-31g opt td（singlets,root=1）这样就是优化s1。
#p b3lyp/6-31g opt td（triplets,root=1）这样是t1。
4：如果说计算出错，可以参考我的那几个帖子。有解决的方式，比如说加入direct等等

yopoxia

小范范1989 发表于 2014-12-25 15:16
你好，很高兴和你讨论问题。
1：这个root=n，我的感觉，就是对那个态感兴趣，比如说，我要优化第一单重 ...

嗯  非常感谢您的热心回复。
不过还是有地方不明白，希望您能继续解答，谢谢！
1、如果我是要算荧光的话，是选择哪个激发态呢？是（singlets，root=1）还是（singlets，root=2）还是其他或者（triplets，root=1）（triplets，root=2）
2、您说的第三条我好像明白了。我下面说例子，您看我理解的对不
比如我想算水的。
高斯里面画好，得到gif如下%chk=H2O.chk
%mem=6GB
%nprocshared=8
#p opt b3lyp/6-31G

H2O

0 1
O                  X1 Y1 Z1
H                  X2 Y2 Z2
H                  X3 Y3 Z3

然后再频率计算计算，取优化好的坐标，再计算，计算所用gjf内容如下
%chk=H2O.chk
%mem=6GB
%nprocshared=8
#p freq b3lyp/6-31G

H2O

0 1
O                  X1a Y1a Z1a
H                  X2a Y2a Z2a
H                  X3a Y3a Z3a

这样一来，优化opt和频率freq计算就做好了。
然后再用这个最新的坐标，再计算，gjf内容如下写
%chk=H2O.chk
%mem=6GB
%nprocshared=8
#p b3lyp/6-31g opt td（singlets,root=1）

H2O

0 1
O                  X1b Y1b Z1b
H                  X2b Y2b Z2b
H                  X3b Y3b Z3b

我的疑问是，上面计算的结果，就是最后一轮最大f对应的就是发射峰的位置了么？不是激发能？这个（singlets，root=1）还是（singlets，root=2）还是其他或者（triplets，root=1）（triplets，root=2）是怎么确定写哪个呢？

sobereva

yopoxia 发表于 2014-12-25 21:42
嗯  非常感谢您的热心回复。
不过还是有地方不明白，希望您能继续解答，谢谢！
1、如果我是要算荧光的 ...

荧光算的是singlets，root=1。root>1不符合kasha规则，triplet则对应的是磷光
最后一轮激发态优化产生的第一激发能就是荧光最大峰位置的近似。f对应的是强度而非峰的位置。

小范范1989

yopoxia 发表于 2014-12-25 21:42
嗯  非常感谢您的热心回复。
不过还是有地方不明白，希望您能继续解答，谢谢！
1、如果我是要算荧光的 ...

对，sob写的完全对，就是sob的意思。
荧光是从s1发出的，磷光是从t1发出的。
但是也不排除从高激发态发出荧光。
但是，一般都是从s1发错的。
所以，我们算荧光发射的时候，都是td（singlets，root=1），
singlets代表单重态，root=1代表第一个，所以加起来就是s1。
那么，最后一步，你在做#p b3lyp/6-31g opt td（singlets,root=1）的时候，计算结果会有3个td（默认的）
找到第一个td（一定是文件的最后的那个，可以从文件最低往上找），对应的能量，转化为波长，就是荧光波长，看看下面的图