CCSD不收敛怎么解决

mfdsrax2

看了教程，只有处理SCF不收敛，优化不收敛，IRC不收敛的解决方法，没有CCSD不收敛的解决方法。CCSD部分的输出文件在附件里，这个应该怎么解决？

zjxitcc

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

1. grep 'E(CORR)=' ccsd.txt

复制代码

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

1. #p CCSD(SaveAmplitudes)/cc-pVDZ
   
2. 
   
3. --Link1--
   
4. #p CCSD(ReadAmplitudes)/cc-pVTZ guess=read geom=allcheck
   

复制代码

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

mfdsrax2

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

能量忽大忽小，居然还有1.23D+05，那就只能从从cc-pvdz读振幅了？

scf

molpro里有shift和diis的一些参数可以调整，如果Gaussian有些iop能调整可能还有点办法，问问客服？

biogon

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

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

mfdsrax2

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

请问什么是多参考特征

zjxitcc

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

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

mfdsrax2

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

我搜了一下，多参考态特征就一个帖子，多参考特征就3个，还是不知道是怎么一回事

snljty2

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

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

mfdsrax2

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

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

zjxitcc

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

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

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