计算化学公社

标题: TDDFT计算激发态，S1结果不符的问题 [打印本页]

作者
Author: yuyao_yang 时间: 2022-7-22 13:23
标题: TDDFT计算激发态，S1结果不符的问题
用g16 TDDFT方法算了一个分子的激发态，看最后结果的时候，感觉怪怪的，S1的分子内部的键并没有被拉伸的感觉？求问是哪出来问题吗？

基态优化关键词：# B3LYP/def2TZVP opt
S1优化关键词：# opt TD(nstates=5,root=1)/def2tzvp freq

附上优化前基态结构1.mol2 和优化后结构2.mol2

作者
Author: yuyao_yang 时间: 2022-7-22 13:24
S1的结构大概长下面这样

作者
Author: pwzhou 时间: 2022-7-22 13:53
你说的不符是不符合你的预期还是不符合实验结果？

作者
Author: Freeman 时间: 2022-7-22 14:23
算激发态，没有设定泛函，大概默认用tdhf了，结果肯定不准。

作者
Author: yuyao_yang 时间: 2022-7-22 14:56
本帖最后由 yuyao_yang 于 2022-7-22 15:03 编辑

原来这样，感谢！

作者
Author: yuyao_yang 时间: 2022-7-23 09:28

Freeman 发表于 2022-7-22 14:23
算激发态，没有设定泛函，大概默认用tdhf了，结果肯定不准。

重现设定了泛函，用了pbe0/def2tzvp opt td=(singlet, nstate=3, root=1) freq=noraman,感觉结果还是这样，而且freq计算，只有一个分子，没看到多个，请问这可能是我哪里设置错了吗？

作者
Author: yuyao_yang 时间: 2022-7-23 09:33

pwzhou 发表于 2022-7-22 13:53
你说的不符是不符合你的预期还是不符合实验结果？

我在文献里看到了一个4CN-indole的S1结构，我的结构就在3位多了一个甲基，但是它文献里面的S1结构苯环上碳原子之间的键被拉得特别长，所以我就怀疑是不是我的计算有问题？

作者
Author: pwzhou 时间: 2022-7-23 13:26

yuyao_yang 发表于 2022-7-23 09:33
我在文献里看到了一个4CN-indole的S1结构，我的结构就在3位多了一个甲基，但是它文献里面的S1结构苯环上 ...

文献上是用什么方法做的？文献上的S1和你的计算的S1性质上是一样的吗？有时候不同方法计算的激发态顺序会不一致，最好给出更多的详细信息。
另外，gaussian里面的PBE0的写法应该是PBE1PBE，你直接写PBE0应该会报错的吧？

作者
Author: yuyao_yang 时间: 2022-7-23 14:20

pwzhou 发表于 2022-7-23 13:26
文献上是用什么方法做的？文献上的S1和你的计算的S1性质上是一样的吗？有时候不同方法计算的激发态顺序会 ...

文献上用的是CC2/cc-pVTZ方法，是的我在关键词用的是# pbe1pbe/def2tzvp opt td=(singlet, nstate=3, root=1) freq=noraman。

作者
Author: sobereva 时间: 2022-7-23 15:03

yuyao_yang 发表于 2022-7-23 09:28
重现设定了泛函，用了pbe0/def2tzvp opt td=(singlet, nstate=3, root=1) freq=noraman,感觉结果还是这样 ...

甭写毫无意义的noraman
仔细看
常见的多余的和被滥用的Gaussian关键词
http://sobereva.com/331（http://bbs.keinsci.com/thread-3460-1-1.html）

根本就不该写PBE0关键词，仔细看
http://bbs.keinsci.com/thread-13660-1-1.html

你若要重现文献就先用和文献相同的体系计算，重现成功再说自己的体系
如果证实TDDFT不理想，考虑用EOM-CCSD优化

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