计算化学公社

标题: 同一基态结构下的激发态优化b3lyp和cam-b3lyp的结果差别十分大 [打印本页]

作者
Author: 一起挖煤 时间: 2020-4-9 14:51
标题: 同一基态结构下的激发态优化b3lyp和cam-b3lyp的结果差别十分大
请问我在已经用b3lyp/6-31G*优化之后的基态结构之上分别使用b3lyp/6-31G*和cam-b3lyp/6-31G*进行激发态优化，结果发现优化出的结果不仅仅构型完全不同，而且out文件中的结果也是完全不同。我感觉b3lyp优化的激发态结果更符合我的实验结果但是对于我的分子受到激发后是电荷转移激发类型的这样的分子不是应该cam-b3lyp的结果不是应该更好吗？

作者
Author: liyuanhe211 时间: 2020-4-9 17:24
这很可能是追踪到不同的电子态去了，你CAM-B3LYP的S2是B3LYP的S1，CAM-B3LYP的S1是B3LYP的S2，电子态都不一样显然激发态极小点结构和激发态的发射光谱都有很大区别。应当在优化时适当的选择root。

作者
Author: 一起挖煤 时间: 2020-4-9 18:33

liyuanhe211 发表于 2020-4-9 17:24
这很可能是追踪到不同的电子态去了，你CAM-B3LYP的S2是B3LYP的S1，CAM-B3LYP的S1是B3LYP的S2，电子态都不一 ...

谢谢李老师，确实跟您说的一样，B3LYP优化出的激发态结构中两个苯环是垂直的，而在CAM-B3LYP的优化结构中两个苯环是基本共平面的。
请问老师
1 我能否使用B3LYP的这个激发态优化结果，尽管我的分子受激是一个电荷转移激发的类型，但是我觉得这个结果对我的实验结果来说比较合理。
2 您说的合适的选择root,这点我不明白您的意思，因为我是研究它的荧光发射所以我就是都默认的1，您能否解释一下我该怎么操作

作者
Author: liyuanhe211 时间: 2020-4-9 18:44

一起挖煤发表于 2020-4-9 18:33
谢谢李老师，确实跟您说的一样，B3LYP优化出的激发态结构中两个苯环是垂直的，而在CAM-B3LYP的优化结构中 ...

root设几在优化时就追踪哪个态。是荧光发射就都应该设1是错误的，首先应当定好你要研究哪个态，而不是默认的出来啥是啥。
不应当现在就考虑能不能按照实验结果“选泛函”，因为你其中一个计算可以说是错误的，先纠正了再说。

作者
Author: 一起挖煤 时间: 2020-4-9 19:32

liyuanhe211 发表于 2020-4-9 18:44
root设几在优化时就追踪哪个态。是荧光发射就都应该设1是错误的，首先应当定好你要研究哪个态，而不是默 ...

谢谢李老师可是一般荧光发射遵循卡莎规则，是由S0与S1态之间进行的跃迁所以不是应该就是研究S1态吗

作者
Author: chiweijie 时间: 2020-4-9 20:16
B3LYP优化的激发态结构很可能是有问题的。

作者
Author: liyuanhe211 时间: 2020-4-9 20:35

一起挖煤发表于 2020-4-9 19:32
谢谢李老师可是一般荧光发射遵循卡莎规则，是由S0与S1态之间进行的跃迁所以不是应该就是研究S1态吗

一方面其实它不一定对，第二是不同几何结构下的S1可不一定都是同一个电子态。

作者
Author: 一起挖煤 时间: 2020-4-9 20:38
本帖最后由一起挖煤于 2020-4-9 20:40 编辑

chiweijie 发表于 2020-4-9 20:16
B3LYP优化的激发态结构很可能是有问题的。

谢谢您的回复。我用B3LYP优化出的结构是一个TICT态而且波长与我的预想差不多，比较符合我的实验结果也确实是一个TICT态，但是CAM-B3LYP的波长差的太远了。我也知道一般对于这种电荷转移激发的激发态优化 B3LYP的结果确实是不太好 CAM-B3LYP和WB97XD的优化结果会比较好。所以我现在就很迷茫了

作者
Author: 一起挖煤 时间: 2020-4-9 20:48

liyuanhe211 发表于 2020-4-9 20:35
一方面其实它不一定对，第二是不同几何结构下的S1可不一定都是同一个电子态。

那您的意思是我应该对CAM-B3LYP 去ROOT S2吗不好意思李老师我刚接触这方面半个月问的问题很小白还望您谅解

作者
Author: liyuanhe211 时间: 2020-4-9 21:25
本帖最后由 liyuanhe211 于 2020-4-9 21:28 编辑

一起挖煤发表于 2020-4-9 20:48
那您的意思是我应该对CAM-B3LYP 去ROOT S2吗不好意思李老师我刚接触这方面半个月问的问题很小白 ...

我不知道你想要哪个态，你自己看激发涉及的轨道等情况来定你想要哪个态，电子态对应哪个态就root=几，有可能一开始root=2，优化完了就变成S1了也不一定；或者反过来基态极小点激发的时候是S1，结果在这个激发态电子态的势能面上的极小点处它成了S2了也有可能。

(, 下载次数 Times of downloads: 42)

作者
Author: chiweijie 时间: 2020-4-9 21:28

一起挖煤发表于 2020-4-9 20:38
谢谢您的回复。我用B3LYP优化出的结构是一个TICT态而且波长与我的预想差不多，比较符合我的实验结果也确 ...

把结构摆出来，看结构也猜的八九不离十。如果你用B3LYP算这类分子，很多没有TICT的分子也会计算出来TICT.

作者
Author: 一起挖煤 时间: 2020-4-9 21:36

chiweijie 发表于 2020-4-9 21:28
把结构摆出来，看结构也猜的八九不离十。如果你用B3LYP算这类分子，很多没有TICT的分子也会计算出来TICT.

C:\Users\Lenovo\Desktop\QQ截图20200409213423.png

作者
Author: chiweijie 时间: 2020-4-9 21:41

一起挖煤发表于 2020-4-9 21:36

图片没贴上吧

作者
Author: 一起挖煤 时间: 2020-4-9 21:43
本帖最后由一起挖煤于 2020-4-9 22:11 编辑

结构如图

作者
Author: 一起挖煤 时间: 2020-4-9 21:43

chiweijie 发表于 2020-4-9 21:41
图片没贴上吧

不好意思放在下面了

作者
Author: 一起挖煤 时间: 2020-4-9 21:45

liyuanhe211 发表于 2020-4-9 21:25
我不知道你想要哪个态，你自己看激发涉及的轨道等情况来定你想要哪个态，电子态对应哪个态就root=几，有 ...

谢谢老师我按照您说的试一试

作者
Author: chiweijie 时间: 2020-4-9 21:53

一起挖煤发表于 2020-4-9 21:43
结构如图

很有可能TICT

作者
Author: 一起挖煤 时间: 2020-4-9 21:56

chiweijie 发表于 2020-4-9 21:53
很有可能TICT

它的实验现象确实是一个暗态的，您觉得B3LYP的结果可信度大吗

作者
Author: chiweijie 时间: 2020-4-9 22:00

一起挖煤发表于 2020-4-9 21:56
它的实验现象确实是一个暗态的，您觉得B3LYP的结果可信度大吗

不建议B3LYP, M06-2X结合CLR或SS溶剂模型基本就可以拿到合理的结果了.

作者
Author: 一起挖煤 时间: 2020-4-9 22:09

chiweijie 发表于 2020-4-9 22:00
不建议B3LYP, M06-2X结合CLR或SS溶剂模型基本就可以拿到合理的结果了.

谢谢您的指导，我会用您说的方法重新计算一遍。

作者
Author: sobereva 时间: 2020-4-10 19:19
此文没看过的话好好看看
谈谈势能面交叉对激发态优化的影响
http://sobereva.com/468（http://bbs.keinsci.com/thread-12318-1-1.html）

作者
Author: 一起挖煤 时间: 2020-4-11 17:29

sobereva 发表于 2020-4-10 19:19
此文没看过的话好好看看
谈谈势能面交叉对激发态优化的影响
http://sobereva.com/468（http://bbs.keinsc ...

谢谢 sob老师
我还有一个问题想请教您，我后来再在使用cam-b3lyp做TD opt的时候加了一个弥散，然后发现它优化出的结果与b3lyp做TD opt的结果类似了请问我这个做法可取吗然后它为什么会这样

作者
Author: sobereva 时间: 2020-4-11 23:07

一起挖煤发表于 2020-4-11 17:29
谢谢 sob老师
我还有一个问题想请教您，我后来再在使用cam-b3lyp做TD opt的时候加了一个弥散，然后发现 ...

弥散函数一般不会显著影响低阶、非里德堡激发的势能面的形状，但没准会令激发能顺序发生变化，先搞清楚原因是什么
没有什么可取不可取，本来B3LYP的结果又不是金标准。如果想先把基组层面的影响彻底撇清，就直接上个3-zeta档次的质量较好的基组，诸如def-TZVP。

作者
Author: 一起挖煤 时间: 2020-4-15 11:44

chiweijie 发表于 2020-4-9 22:00
不建议B3LYP, M06-2X结合CLR或SS溶剂模型基本就可以拿到合理的结果了.

您好 m062x与cam-b3lyp以及wb97xd得到的结果基本差不多，任然是很不符合实验数据的

作者
Author: chiweijie 时间: 2020-4-15 13:10

一起挖煤发表于 2020-4-15 11:44
您好 m062x与cam-b3lyp以及wb97xd得到的结果基本差不多，任然是很不符合实验数据的

采用了哪个溶剂模型？SS, CLR 还是LR?

作者
Author: 一起挖煤 时间: 2020-4-16 16:45

sobereva 发表于 2020-4-11 23:07
弥散函数一般不会显著影响低阶、非里德堡激发的势能面的形状，但没准会令激发能顺序发生变化，先搞清楚原 ...

请问sob老师，对于同一个结构的势能面扫描输入文件，仅仅改变其中的函数，所做的势能面扫描结果有很大的差异这是正常的吗？

作者
Author: sobereva 时间: 2020-4-19 00:48

一起挖煤发表于 2020-4-16 16:45
请问sob老师，对于同一个结构的势能面扫描输入文件，仅仅改变其中的函数，所做的势能面扫描结果有很大的 ...

那叫泛函不叫函数
不同DFT泛函下势能面可能有着定性的差异，一点也不奇怪

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