计算化学公社

标题: Gaussian计算有机反应机理中三线态能量的问题 [打印本页]

作者
Author: BOBOCHI05 时间: 2017-10-24 20:21
标题: Gaussian计算有机反应机理中三线态能量的问题
本帖最后由 BOBOCHI05 于 2017-10-24 20:23 编辑

各位老师，我想请教一下，我要研究一个有机反应的机理，如下图所示，我想要知道catalyst在T1状态的能量能否将reactant从S0激发到T1，是不是应该计算catalyst从S0到T1的绝热吸收，再计算reactant从S0到T1的绝热吸收，然后比较catalyst和reactant的绝热吸收。
如果是计算绝热吸收的话，看过论坛之前的帖子是这样计算的：
绝热吸收：电子从基态被激发到激发态，结构也从基态极小点结构变化到激发态极小点结构。计算绝热吸收能包含以下三步
1.优化基态几何结构，在最后一步得到基态极小点能量
2.优化激发态几何结构，在最后一步得到激发态极小点能量
3.将第二步得到的能量与第一步得到的能量求差值

我想请教的是上述2中计算的设置应该是哪种：(1) opt td=(triplets,nstates=5,root=1) b3lyp/6-31+g(d) charge=0 Multiplicity=1
(2) opt td=(triplets,nstates=5,root=1) b3lyp/6-31+g(d) charge=0 Multiplicity=3

作者
Author: sobereva 时间: 2017-10-25 14:49
基态自旋多重度是几Multiplicity就是几，TDDFT是以基态作为参考态计算激发态的

作者
Author: BOBOCHI05 时间: 2017-10-25 19:33
老师好，我的体系基态的自旋多重度是1，我用opt td=(triplets,nstates=5,root=1) b3lyp/6-31g(d) charge=0 Multiplicity=1作为关键词计算上述绝热吸收的第2步，让我疑惑的是优化过程中能量是不断升高的，如图片1所示，优化出来的结构苯环都不在一个平面上了，如图片2所示，这步的输出文件为附件中的文件，我不知道这个计算哪里出了问题，还请老师指点一下，谢谢。

作者
Author: sobereva 时间: 2017-10-26 12:09

BOBOCHI05 发表于 2017-10-25 19:33
老师好，我的体系基态的自旋多重度是1，我用opt td=(triplets,nstates=5,root=1) b3lyp/6-31g(d) charge= ...

T1你直接用0 3优化多好，TDDFT白费很多时间，准确度还不如UKS
三重态时候体系结构发生显著扭曲、变形是很正常的。
你读的能量是基态参考态（S0）的能量，因为优化过程中逐渐偏离S0，S0能量升高是显然的。

作者
Author: BOBOCHI05 时间: 2017-10-26 16:18

sobereva 发表于 2017-10-26 12:09
T1你直接用0 3优化多好，TDDFT白费很多时间，准确度还不如UKS
三重态时候体系结构发生显著扭曲、变形 ...

老师，绝热吸收的具体步骤是不是这样的（体系基态电荷为0，自旋多重度为1）：
1.优化基态几何结构 #opt b3lyp/6-31g(d) charge=0 multiplicity=1
2.在基态平衡几何构型的基础上优化三重态几何结构 #opt b3lyp/6-31g(d) charge=0 multiplicity=3
3.绝热吸收是1中基态极小点能量与2中激发态极小点能量的差值

作者
Author: sobereva 时间: 2017-10-26 22:31

BOBOCHI05 发表于 2017-10-26 16:18
老师，绝热吸收的具体步骤是不是这样的（体系基态电荷为0，自旋多重度为1）：
1.优化基态几何结构 #op ...

是

作者
Author: BOBOCHI05 时间: 2017-10-27 10:00

sobereva 发表于 2017-10-26 22:31
是

多谢老师~

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