审稿人问有机自由基自旋污染S2值

Zheng123

有一条审稿意见提到了有机自由基的S2值在使用CCSD（t）方法时，有机自由基中的s2值会与0.75相差很大，我看了我的pvtz文件在0.75到0.99的都有，这数值正常吗，（DFT计算输出文件的S2值大概都在0.75左右）审稿人说让制成表格。顺便想请教一下这条意见该怎么回复

zjxitcc

注意，如果你是用的Gaussian做UCCSD(T)计算，输出文件中<S**2>是UHF的，不是UCCSD(T)的，切勿错误指认（给审稿人进一步授之以柄）。如果是二重态自由基，<S**2>为0.99的话，自旋污染不是很大（意思就是偏离理想值0.75不是很多，不是很要紧）。如果结构只有几十个（计算量不大），建议全部检验波函数稳定性（如果不懂这个意思，请在该论坛上搜索，往期有无数帖子提到），然后回复的时候跟审稿人说全部检验过波函数稳定性，并且本文给出的结果 对应的波函数都是稳定的。当然，如果发现某一个不稳定，则要重新计算/继续优化该结构。

正常回复，列表就列表。

Zheng123

zjxitcc 发表于 2024-7-13 21:06
注意，如果你是用的Gaussian做UCCSD(T)计算，输出文件中是UHF的，不是UCCSD(T)的，切勿错误指认（给审稿人 ...

啊我CCSD（T）计算实在gaussian里做的，那请问这个S2值从哪查找呢，我看了别的帖子误以为<S**2>就是CCSD（t）的s2值

zjxitcc

Zheng123 发表于 2024-7-13 21:37
啊我CCSD（T）计算实在gaussian里做的，那请问这个S2值从哪查找呢，我看了别的帖子误以为就是CCSD（t）的 ...

据我所知没有程序可以计算UCCSD(T)的<S**2>。你所列的表里要有一行是UHF的<S**2>。

Zheng123

zjxitcc 发表于 2024-7-13 21:40
据我所知没有程序可以计算UCCSD(T)的。你所列的表里要有一行是UHF的。

好的谢谢 我回复意见里可以这么说吗 列表里放DFT和UHF的

zjxitcc

Zheng123 发表于 2024-7-13 21:55
好的谢谢 我回复意见里可以这么说吗 列表里放DFT和UHF的

可以，写的时候要写明具体泛函。例如用的是B3LYP泛函，如果是二重态体系，应写明UB3LYP。

Zheng123

zjxitcc 发表于 2024-7-13 22:36
可以，写的时候要写明具体泛函。例如用的是B3LYP泛函，如果是二重态体系，应写明UB3LYP。

好的 十分感谢！