TDDFT计算激发态，S1结果不符的问题

yuyao_yang · 发表于 Post on 2022-7-22 13:23:24

用g16 TDDFT方法算了一个分子的激发态，看最后结果的时候，感觉怪怪的，S1的分子内部的键并没有被拉伸的感觉？求问是哪出来问题吗？

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

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

yuyao_yang · 发表于 Post on 2022-7-22 13:24:40

S1的结构大概长下面这样

pwzhou · 发表于 Post on 2022-7-22 13:53:10

你说的不符是不符合你的预期还是不符合实验结果？

Freeman · 发表于 Post on 2022-7-22 14:23:21

算激发态，没有设定泛函，大概默认用tdhf了，结果肯定不准。

yuyao_yang · 发表于 Post on 2022-7-22 14:56:33

本帖最后由 yuyao_yang 于 2022-7-22 15:03 编辑

原来这样，感谢！

yuyao_yang · 发表于 Post on 2022-7-23 09:28:16

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

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

yuyao_yang · 发表于 Post on 2022-7-23 09:33:14

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

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

pwzhou · 发表于 Post on 2022-7-23 13:26:03

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

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

yuyao_yang · 发表于 Post on 2022-7-23 14:20:58

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

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

sobereva · 发表于 Post on 2022-7-23 15:03:54

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优化

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Gaussian/gview] TDDFT计算激发态，S1结果不符的问题

浏览过的版块