计算化学公社

标题: CCSD不收敛怎么解决 [打印本页]

作者
Author: mfdsrax2 时间: 2022-6-16 18:20
标题: CCSD不收敛怎么解决
看了教程，只有处理SCF不收敛，优化不收敛，IRC不收敛的解决方法，没有CCSD不收敛的解决方法。CCSD部分的输出文件在附件里，这个应该怎么解决？

作者
Author: zjxitcc 时间: 2022-6-16 18:42
本帖最后由 zjxitcc 于 2022-6-16 18:46 编辑

分析振幅迭代过程中能量变化，如果最后10圈能量一直在下降，且即将收敛，可以适当增加圈数，例如CCSD(maxcyc=100)。如何分析？运行以下命令

grep 'E(CORR)=' ccsd.txt

复制代码

如果最后十几圈一直在振荡，那增大maxcyc基本没有用处，此时可以尝试读取中小基组的振幅，举个例子

#p CCSD(SaveAmplitudes)/cc-pVDZ
--Link1--
#p CCSD(ReadAmplitudes)/cc-pVTZ guess=read geom=allcheck

复制代码

如果要做CCSD(T)计算，可以将最后一个(ReadAmplitudes)改为(ReadAmplitudes,T)，注意不要去改第一个(SaveAmplitudes)。如果各种技巧反复尝试都不行（此处尝试不含无脑maxcyc=1000），大概率当前体系有多组态/多参考态特征，常规CCSD无法准确描述。

作者
Author: mfdsrax2 时间: 2022-6-16 18:51
本帖最后由 mfdsrax2 于 2022-6-16 18:52 编辑

zjxitcc 发表于 2022-6-16 18:42
分析振幅迭代过程中能量变化，如果最后10圈能量一直在下降，且即将收敛，可以适当增加圈数，例如CCSD(maxcy ...

能量忽大忽小，居然还有1.23D+05，那就只能从从cc-pvdz读振幅了？
(, 下载次数 Times of downloads: 16)

作者
Author: scf 时间: 2022-6-17 00:02
molpro里有shift和diis的一些参数可以调整，如果Gaussian有些iop能调整可能还有点办法，问问客服？

作者
Author: biogon 时间: 2022-6-17 11:30

mfdsrax2 发表于 2022-6-16 18:51
能量忽大忽小，居然还有1.23D+05，那就只能从从cc-pvdz读振幅了？

这种还很有可能是体系有多参考特征了

作者
Author: mfdsrax2 时间: 2022-6-17 15:13

biogon 发表于 2022-6-17 11:30
这种还很有可能是体系有多参考特征了

请问什么是多参考特征

作者
Author: zjxitcc 时间: 2022-6-17 15:16

mfdsrax2 发表于 2022-6-17 15:13
请问什么是多参考特征

论坛上搜索即可，以往有大量帖子。

作者
Author: mfdsrax2 时间: 2022-6-17 17:57

zjxitcc 发表于 2022-6-17 15:16
论坛上搜索即可，以往有大量帖子。

我搜了一下，多参考态特征就一个帖子，多参考特征就3个，还是不知道是怎么一回事
(, 下载次数 Times of downloads: 26)
(, 下载次数 Times of downloads: 34)

作者
Author: snljty2 时间: 2022-6-17 18:27

mfdsrax2 发表于 2022-6-17 17:57
我搜了一下，多参考态特征就一个帖子，多参考特征就3个，还是不知道是怎么一回事

用个靠谱的搜索引擎。比如这样搜出来至少几百个贴。

作者
Author: mfdsrax2 时间: 2022-6-17 20:31

zjxitcc 发表于 2022-6-17 15:16
论坛上搜索即可，以往有大量帖子。

就算确定是多参考特征，然后怎么办？我在算反应机理，其他体系都是用CCSD（T）算的，就这个算不了

作者
Author: zjxitcc 时间: 2022-6-17 20:46
本帖最后由 zjxitcc 于 2022-6-17 20:50 编辑

mfdsrax2 发表于 2022-6-17 20:31
就算确定是多参考特征，然后怎么办？我在算反应机理，其他体系都是用CCSD（T）算的，就这个算不了

计算多参考态特征/诊断指标，有很多种指标，那些帖子里都有说，选一两种指标算即可；如果多参考态特征较为明显，需要做多组态/多参考态计算。
对你算过的那些例子，也可以挑一两个算一下多参考态特征/诊断指标。你当前能算不等于算得合理，收敛不等于结果可信。

另外，如果多参考态特征较弱，或者是简单的有机双自由基（当然，这句话就假设你是单重态了），图方便省事的话，可以用基于UHF的UCCSD(T)，注意UHF要得到自旋破缺单重态。你当前给的信息十分少，无法判断你用的是基于RHF的RCCSD(T)亦或是基于UHF的UCCSD(T)。

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