求解，cfour计算是否可以用guess=mix

commanderwilson · 发表于 Post on 2025-6-21 05:28:29

是这样的，我试了一下用C4的CCSD(T)解析梯度算臭氧，用psi4的接口直接用的，计算加了关键词
reference uhf
cfour_hfstability=follow
进行计算，结果在检测到波函数稳定性以后，再次计算SCF就无法收敛了
于是把参考态换成了rhf，结果计算一下子就收敛了：

1.inp (309 Bytes, 下载次数 Times of downloads: 0)

1.out (439.74 KB, 下载次数 Times of downloads: 1)

2.inp (309 Bytes, 下载次数 Times of downloads: 0)

2.out (115.7 KB, 下载次数 Times of downloads: 0)
我知道用rhf的结果可能不行，所以我想问一下如何在cfour一边用UHF一边混合HOMO和LUMO，做自旋极化单重态的SCF

sobereva · 发表于 Post on 2025-6-22 08:10:52

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

commanderwilson · 发表于 Post on 2025-6-22 14:15:59

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

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

wzkchem5 · 发表于 Post on 2025-6-22 16:33:04

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

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

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

[MRCC/CFOUR] 求解，cfour计算是否可以用guess=mix

浏览过的版块