计算化学公社

标题: 跃迁计算问题 [打印本页]

作者
Author: Theor_Comp 时间: 2017-12-4 22:46
标题: 跃迁计算问题
请问，优化激发态之后，比如S1态，得到输出结果中：Excited state 1的激发能代表的是垂直发光的能量。那么问题一：Excited state 2 以及Excited state 3等等激发能是否有意义，代表的是什么？问题二：在优化S1构型之后，在该结构上直接进行TD计算，得到的是从S1极小点向它的上能级（S2、S3等）跃迁的激发能以及振子强度吗？问题三：优化T1的结构后，在该结构直接做TD计算呢？

作者
Author: sobereva 时间: 2017-12-4 23:52
1 如果是指S1优化最后一步的，相当于S1极小点结构下基态到S2、S3的垂直激发能。没太大意义
2 同上
3 T1极小点结构下基态到各个激发态的激发能等信息

作者
Author: Theor_Comp 时间: 2017-12-4 23:57

sobereva 发表于 2017-12-4 23:52
1 如果是指S1优化最后一步的，相当于S1极小点结构下基态到S2、S3的垂直激发能。没太大意义
2 同上
3 T1极 ...

Sob老师您好，那能否计算S1极小点结构下S1态到其他各个激发态的激发能（S1-->S2、S1-->S3）或者是T1极小点结构下T1态到其他各个激发态的激发能（T1-->T2、T1-->T3）？？

作者
Author: sobereva 时间: 2017-12-5 00:00

Theor_Comp 发表于 2017-12-4 23:57
Sob老师您好，那能否计算S1极小点结构下S1态到其他各个激发态的激发能（S1-->S2、S1-->S3）或者是T1极小 ...

可以。你把那几个激发态的激发能求差即可。

作者
Author: Theor_Comp 时间: 2017-12-5 00:02

sobereva 发表于 2017-12-5 00:00
可以。你把那几个激发态的激发能求差即可。

Sob老师，主要是想得到S1-->S2、S1-->S3这些跃迁的振子强度f，这样才能得到光谱，不知gaussian能否做到？

作者
Author: sobereva 时间: 2017-12-5 00:03

Theor_Comp 发表于 2017-12-5 00:02
Sob老师，主要是想得到S1-->S2、S1-->S3这些跃迁的振子强度f，这样才能得到光谱，不知gaussian能否做到？

我并不清楚这么做有什么意义，激发不会从激发态激发。如果真是有极特殊情况，根据激发能还有激发态之间的跃迁偶极矩去计算f。激发态之间跃迁偶极矩的计算看
利用Multiwfn计算Gaussian输出的激发态之间的跃迁偶极矩
http://sobereva.com/227

作者
Author: Theor_Comp 时间: 2017-12-5 00:09

sobereva 发表于 2017-12-5 00:03
我并不清楚这么做有什么意义，激发不会从激发态激发。如果真是有极特殊情况，根据激发能还有激发态之间的 ...

谢谢Sob老师。
另外，优化S0结构后，在该结构上进行TD计算：
1. #p td=(root=1，nstates=5) b3lyp/6-31g*
2. #p td=(root=3，nstates=5,) b3lyp/6-31g*
这两种方式有什么不同呀？

作者
Author: sobereva 时间: 2017-12-5 01:22

Theor_Comp 发表于 2017-12-5 00:09
谢谢Sob老师。
另外，优化S0结构后，在该结构上进行TD计算：
1. #p td=(root=1，nstates=5) b3lyp/6-31 ...

结果一样

仔细看
Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法
http://sobereva.com/314

作者
Author: Theor_Comp 时间: 2017-12-5 09:57

sobereva 发表于 2017-12-5 01:22
结果一样

仔细看

（一）如果是想得到T1在S1极小点结构下的能量，是需要先优化S1，然后再算单点能吗？
即 #p opt td b3lyp/6-31g* scrf=()
0 1
#p b3lyp/6-31g* scrf=()
0 3
是这样吗？
（二）如果想得到T2态在S1极小点结构下的能量呢？

谢谢Sob老师及各位。

作者
Author: sobereva 时间: 2017-12-5 16:57

Theor_Comp 发表于 2017-12-5 09:57
（一）如果是想得到T1在S1极小点结构下的能量，是需要先优化S1，然后再算单点能吗？
即 #p opt td b3lyp ...

1 是
2 S1极小点下0 1 td(triplet)，把T2激发能与基态能量加和

作者
Author: Theor_Comp 时间: 2017-12-6 09:08
本帖最后由 Theor_Comp 于 2017-12-6 09:27 编辑

sobereva 发表于 2017-12-5 16:57
1 是
2 S1极小点下0 1 td(triplet)，把T2激发能与基态能量加和

Sob老师，另外：一、那关于电荷-自旋那的输入，什么情况下输入1，什么情况下是3呢？感觉除了使用UDFT优化T1时是输入3外，其他情况都是1吗？二、基于S1的极小点结构直接做单点计算（#p b3lyp/6-31g* scrf=()）或者是TD计算（#p td=nstates=5 b3lyp/6-31g* scrf=()）时，scrf里加不加eqsolv呢？
多谢Sob老师和各位了。

作者
Author: sobereva 时间: 2017-12-6 13:29

Theor_Comp 发表于 2017-12-6 09:08
Sob老师，另外：一、那关于电荷-自旋那的输入，什么情况下输入1，什么情况下是3呢？感觉除了使用UDFT优化 ...

1 基态（参考态）是什么情况就输入什么情况。看我的博文
2 弄清楚什么是平衡溶剂、非平衡溶剂再考虑加不加，博文里说了非平衡溶剂的问题，不懂就不要写。

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