能否计算S1→S2激发能？

redfoxs

冰释之川 发表于 2017-4-30 06:50
首先这A-B-C是基于S1-Minima来计算的。我们分别将B、A、C三态定义为E0(S1), E1(S1)和E2(S1), 其中括号中 ...

非常感谢冰老师和小范范老师，这个问题我清楚了。
还有几个问题，在算S0的td时我们设置td(nstate=5)是表示算在S0的opt结构下往上最低的五个激发态，即S0→S1，S0→S2，。。S0→S5，这里的root默认是1，那么当我们设置root=4时，即S0时，td(nstate=5，root=4)是什么意思呢，sobereva老师说如果"只是计算激发态信息则不需要管此设定"。这里的激发态信息是指什么，跃迁吗？ 设置root值是不是只有和opt联用才有意义，比如opt td(nstate=5，root=4) 表示优化第四激发态？此外，对opt td(nstate=5，root=2)，是先做td再opt吗，opt 和td的先后顺序（opt td& td opt）有差别吗？

redfoxs

冰释之川 发表于 2017-5-1 08:22
（1）即S0时，td(nstate=5，root=4)是什么意思呢？
和默认的root=1一样，只不过如下信息会变成你设置的第 ...

非常详细，感谢冰老师

redfoxs

冰释之川 发表于 2017-5-1 08:22
（1）即S0时，td(nstate=5，root=4)是什么意思呢？
和默认的root=1一样，只不过如下信息会变成你设置的第 ...

非常感谢冰老师详细的回答，对第四个问题我还是不太明白，
1.对opt td(nstate=5，root=2)，他是在对S2做迭代优化，然后每一步优化都做nstate=5的TDDFT吗？还是在基态（因为我们输入的点群结构是基态的）的结构下做TDDFT，然后对每一个激发态（这里nstate=5, 所有是5个激发态）做优化，但是只给出S2的优化后的结构呢？或者是他对每一个激发态做优化，然后对S2做nstate=5的TDDFT？
2.在opt td(nstate=5，root=2)完成后，HOMO，LUMO轨道能级和电子云分布也是通过fchk文件读出了？

另一个单独的问题是：
3.在opt td(nstate=5，root=1)时，我们得到的Excited state 1是S0-S1的激发能也就是S1-S0的发射能，那么opt td(nstate=5，root=2)中Excited state 1得到的是哪个跃迁的能量了？

redfoxs

冰释之川 发表于 2017-5-1 08:22
（1）即S0时，td(nstate=5，root=4)是什么意思呢？
和默认的root=1一样，只不过如下信息会变成你设置的第 ...

冰老师，还有一个问题，见图片，现在我知道
E(AB)=td
E(CD)=opt td(root=1)
E(DE)=opt td(root=1) 后E2-E1
那么E(FH) 怎么得到了？是否是E(FH)=opt td(root=2)?
能否得到E(HG) ？

谢谢啦

冰释之川

redfoxs 发表于 2017-5-4 22:41
冰老师，还有一个问题，见图片，现在我知道
E(AB)=td
E(CD)=opt td(root=1)

非常感谢冰老师详细的回答，对第四个问题我还是不太明白，
1.对opt td(nstate=5，root=2)，他是在对S2做迭代优化，然后每一步优化都做nstate=5的TDDFT吗？还是在基态（因为我们输入的点群结构是基态的）的结构下做TDDFT，然后对每一个激发态（这里nstate=5, 所有是5个激发态）做优化，但是只给出S2的优化后的结构呢？或者是他对每一个激发态做优化，然后对S2做nstate=5的TDDFT？
每一步优化都做TDDFT（nstates =5），追踪的态是root指定的S2, 其它态不管的，所以优化完以后得到的结构就是S2-Minima


2.在opt td(nstate=5，root=2)完成后，HOMO，LUMO轨道能级和电子云分布也是通过fchk文件读出了？

是的，不过注意，此时的电子密度是 S2-Minima结构下的基态密度。

另一个单独的问题是：
3.在opt td(nstate=5，root=1)时，我们得到的Excited state 1是S0-S1的激发能也就是S1-S0的发射能，那么opt td(nstate=5，root=2)中Excited state 1得到的是哪个跃迁的能量了？



这里的Excited state 1是指 S2-Minima结构下 S0 -> S1的跃迁，不过你不用理他，关注的应该是 Excited state 2

冰释之川

redfoxs 发表于 2017-5-4 22:41
冰老师，还有一个问题，见图片，现在我知道
E(AB)=td
E(CD)=opt td(root=1)

最后那个问题：

E(FH)是这么得到的，取opt td(root=2)最后一次TDDFT输出信息里的Excited state 2激发能。而E(HG)=E(FH) - E(FG)

其中E(FG)取opt td(root=2)最后一次TDDFT输出信息里的Excited state 1激发能

redfoxs

冰释之川 发表于 2017-5-4 23:07
最后那个问题：

E(FH)是这么得到的，取opt td(root=2)最后一次TDDFT输出信息里的Excited state 2激发 ...

谢谢冰老师，这个问题我懂了，就是说不管怎么设置root值，opt td后的结构都是在S0的势能曲线上找到那个对应的结构做的跃迁计算，是吗？

redfoxs

冰释之川 发表于 2017-5-4 23:07
非常感谢冰老师详细的回答，对第四个问题我还是不太明白，
1.对opt td(nstate=5，root=2)， ...

1.对opt td(nstate=5，root=2)，他是在对S2做迭代优化，然后每一步优化都做nstate=5的TDDFT吗？还是在基态（因为我们输入的点群结构是基态的）的结构下做TDDFT，然后对每一个激发态（这里nstate=5, 所有是5个激发态）做优化，但是只给出S2的优化后的结构呢？或者是他对每一个激发态做优化，然后对S2做nstate=5的TDDFT？
每一步优化都做TDDFT（nstates =5），追踪的态是root指定的S2, 其它态不管的，所以优化完以后得到的结构就是S2-Minima

1.就是说其实这个opt是在对S2进行，然后每一步都做TDDFT，一直到opt找到那个能量最低点，计算才结束是吗，和对基态的优化是一样的？
2.如果这样的话，是不是说nstate的值对最后的优化结果没有影响，因为对S2的opt是一直找到最小值才结束

冰释之川

redfoxs 发表于 2017-5-4 23:47
谢谢冰老师，这个问题我懂了，就是说不管怎么设置root值，opt td后的结构都是在S0的势能曲线上找到那个对 ...

你给的结构确定了，那么做TDDFT显然都是垂直激发的，比如你在S2结构下做了TDDFT，输出信息里的各个激发态其实都是S2结构下的垂直激发， root为多少都无所谓了。只有当你要对某个激发态感兴趣（例如优化某个激发态，或者读取某个激发态下的偶极矩，原子电荷分布），你才需要设置root。

冰释之川

redfoxs 发表于 2017-5-5 00:03
1.对opt td(nstate=5，root=2)，他是在对S2做迭代优化，然后每一步优化都做nstate=5的TDDFT吗？还是在基 ...

参考《Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法》中“3. Gaussian中的TDDFT计算”节：http://sobereva.com/314

redfoxs

冰释之川 发表于 2017-5-5 00:07
你给的结构确定了，那么做TDDFT显然都是垂直激发的，比如你在S2结构下做了TDDFT，输出信息里的各个激发态 ...

好的，非常感谢冰老师

redfoxs

冰释之川 发表于 2017-5-4 23:07
非常感谢冰老师详细的回答，对第四个问题我还是不太明白，
1.对opt td(nstate=5，root=2)， ...

2.在opt td(nstate=5，root=2)完成后，HOMO，LUMO轨道能级和电子云分布也是通过fchk文件读出了？

是的，不过注意，此时的电子密度是 S2-Minima结构下的基态密度。


不好意思，冰老师，查了好久也没看到怎么在gaussview下怎么做出S2激发态的HOMO,LUMO电子分布图？

冰释之川

redfoxs 发表于 2017-5-5 03:51
2.在opt td(nstate=5，root=2)完成后，HOMO，LUMO轨道能级和电子云分布也是通过fchk文件读出了？

是的 ...

redfoxs

冰释之川 发表于 2017-5-5 08:37

谢谢冰老师，是这样的，我有一个新星分子，上面有四个发色团，但两个发色团与另两个的连接方式不一样，我一直想能不能模拟是哪一组在发色荧光，刚好看到有篇文章计算了激发态的电子分布，所以我也想试试看在激发态时的电子分布。http://pubs.acs.org/doi/abs/10.1021/jp1109042

http://pubs.acs.org/doi/full/10.1021/jacs.5b04075

http://pubs.rsc.org/en/content/articlehtml/2013/dt/c3dt51571k

冰释之川

要鉴别哪些基团参与了发射过程，对于垂直发射做个电子密度差就可以直观的考差了（例如：http://pubs.acs.org/doi/abs/10.1021/acs.jpcc.7b01360，还有就是【第三幅图】）。当然，做自然跃迁密度（NTO）【第二幅图】，或者直接查看基态和激发态结构下的激发/发射归属以及轨道对贡献【第一幅图】 都是可以表征的