计算化学公社

标题: 关于T2态优化及输出其波函数的问题 [打印本页]

作者
Author: Kalinite 时间: 2021-3-15 01:50
标题: 关于T2态优化及输出其波函数的问题
本帖最后由 Kalinite 于 2021-3-15 01:51 编辑

如题，想要在Gaussian16中优化一个分子的T2态，并通过Multiwfn求其自旋密度差，以确定单电子主要分布在什么位置。这两种方法何者正确呢？
(1) TD(nstates=5,root=2,triplet) density=current out=wfn ... 电荷0 自旋多重度1
(2) TD(nstates=5,root=2,triplet) density=current out=wfn ... 电荷0 自旋多重度3
看到这篇帖子(http://bbs.keinsci.com/thread-10672-1-1.html)中提到应该是用(1)，但是得到的.wfn和.fchk文件用Multiwfn加载均识别为闭壳层波函数，自旋密度处处相等。

问题：
(1) 电荷和自旋多重度指定的是参考态的电荷和自旋多重度吗？
(2) 应该采用哪一种设定进行优化？

作者
Author: sobereva 时间: 2021-3-15 04:26
(1)正确。因为参考态是单重态。处处相当是因为TDDFT根本没法给出自旋密度，应该用Multiwfn算odd electron function，这个事最近我在论坛说过多次，首页google框搜索

2021-Aug-14补充：已经专门写了相关的文章
使用Multiwfn计算odd electron density考察激发态单电子分布
http://sobereva.com/583（http://bbs.keinsci.com/thread-24645-1-1.html）

作者
Author: Kalinite 时间: 2021-3-15 12:26

sobereva 发表于 2021-3-15 04:26
(1)正确。因为参考态是单重态。处处相当是因为TDDFT根本没法给出自旋密度，应该用Multiwfn算odd electron f ...

已经参照手册成功解决，感谢sob老师。
另外还有一个问题，因为T1是UKS算的，因此用的自旋密度差分析的单电子分布，对T1和T2所用方法不一致会不会被审稿人质疑。
我试图也对T1做OED，这时候显示的等值面遍布整个分子，是不是应该用手册上自旋极化单重态的做法写pop=no算一次波函数呢？

作者
Author: sobereva 时间: 2021-3-16 04:35

Kalinite 发表于 2021-3-15 12:26
已经参照手册成功解决，感谢sob老师。
另外还有一个问题，因为T1是UKS算的，因此用的自旋密度差分析的单 ...

当前完全不涉及自旋极化单重态
你也用TDDFT算T1就完了，但结果与UKS有甚至不小的差异并不意外。原理上来说UKS更合理

作者
Author: Kalinite 时间: 2021-4-3 10:09

sobereva 发表于 2021-3-16 04:35
当前完全不涉及自旋极化单重态
你也用TDDFT算T1就完了，但结果与UKS有甚至不小的差异并不意外。原理上来 ...

请问sob老师有没有比较推荐可以引用的，能够佐证T1用UKS优于TDDFT的文献，想在文章中引用。

作者
Author: sobereva 时间: 2021-4-3 15:44

Kalinite 发表于 2021-4-3 10:09
请问sob老师有没有比较推荐可以引用的，能够佐证T1用UKS优于TDDFT的文献，想在文章中引用。

手头上没有现成的文献
我这是从原理上来说，一般懂理论的人都能明白。UKS优化一方面是快得多，另一方面UKS能考虑轨道弛豫（优化轨道）

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