计算化学公社

标题: 求解，cfour计算是否可以用guess=mix [打印本页]

作者
Author: commanderwilson 时间: 2025-6-21 05:28
标题: 求解，cfour计算是否可以用guess=mix
是这样的，我试了一下用C4的CCSD(T)解析梯度算臭氧，用psi4的接口直接用的，计算加了关键词
reference uhf
cfour_hfstability=follow
进行计算，结果在检测到波函数稳定性以后，再次计算SCF就无法收敛了
于是把参考态换成了rhf，结果计算一下子就收敛了：
(, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 1) (, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 0)
我知道用rhf的结果可能不行，所以我想问一下如何在cfour一边用UHF一边混合HOMO和LUMO，做自旋极化单重态的SCF

作者
Author: sobereva 时间: 2025-6-22 08:10
臭氧的静态相关很强，用CCSD(T)优化可能还不如CCSD，若非要用耦合簇优化建议用能明显更好描述静态相关的CCSDT
若必须用CCSD(T)解析梯度优化不可，更建议用同样支持这个特征的ORCA做，用着明显更舒服

作者
Author: commanderwilson 时间: 2025-6-22 14:15

sobereva 发表于 2025-6-22 08:10
臭氧的静态相关很强，用CCSD(T)优化可能还不如CCSD，若非要用耦合簇优化建议用能明显更好描述静态相关的CCS ...

谢谢社长指导。
算臭氧这个例子更多是为了调试，因为最近Nature发了一个常温下合成链状氮原子簇N6的例子，我用ICE-NEVPT2跑了一下CAS(10,10)的活性空间，结果发现这个体系静态相关可能很强，基态的参考态成分只有0.6多一点，而Nature的原论文用的是耦合簇方法优化几何结构，所以我打算做计算来说明它们用的方法可能不大合适。
如果要用DFT处理这种静态相关很强的体系，是优先考虑纯泛函，还是为了获得对称破缺态使用高HF成分的杂化泛函呢

作者
Author: wzkchem5 时间: 2025-6-22 16:33

commanderwilson 发表于 2025-6-22 14:15
谢谢社长指导。
算臭氧这个例子更多是为了调试，因为最近Nature发了一个常温下合成链状氮原子簇N6的例子 ...

获得对称破缺态不是最终目的，没有必要为了获得对称破缺态而选择泛函。
这种体系无论怎么选泛函，即使结果是准的，也有很大碰巧因素。况且体系很小，所以我觉得可以直接用NEVPT2数值梯度优化试试看。
此外这么小的活性空间加ICE加速不大，反倒可能让CASSCF收敛变困难，导致总时间增加

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