求助使用calcall后，优化和频率计算四项收敛标准仍然不同的原因

BLL130

我尝试使用，TPSSh-D3(BJ)结合Def2-TZVP基组，在SMD水溶剂模型下优化闭壳层的[EDTA-Co(III)]-离子（软件为G16 C.02版本），优化结束显示收敛，随后做freq显示最大位移一项未收敛（单独取出最后一帧的结构做的freq，并非自定义IOp未传到到后续任务的问题）。以为是由于优化过程中使用近似Hessian矩阵、而频率计算中使用精确Hessian矩阵所致，故使用"opt=calcall"优化最后一帧结构。
比较奇怪的是"opt=calcall"的结果显示输入文件的这个结构四项标准已经是全部收敛了的，得出了与频率计算任务不同的结果（四项标准的数值均不同，不仅仅是位移不同）。而检查二者给出的单点能、频率值都是相同的，因此不清楚为什么会两个任务给出不同的结果。我十分清楚这个优化出的结构已经是可以用的（即使频率计算，最大位移一项也仅仅是略超出标准），只是好奇会有什么因素使得二者结果产生差别。

附两个任务的输入、输出文件，欢迎大家赐教，谢谢！

hebrewsnabla

首先，两个任务的计算的Force应该是一样的，这个可以通过输出文件中 Forces (Hartree/Bohr)这一段来确认。为什么Max Force/RMS Force不一样呢？我猜是因为freq任务统计的是cartesian force而opt=calcall任务统计的是内坐标force。

至于为什么Displacement不同，应该是两个任务采用的优化算法不同（当然，坐标也不同），这个可以从输出文件第87行看出来。

LittlePupil

有点好奇，如果你把最后一帧结构拿出来，chk路径写频率分析的chk文件，然后opt选项写rcfc而不是calcall会是什么结果？

BLL130

hebrewsnabla 发表于 2022-9-11 21:17
首先，两个任务的计算的Force应该是一样的，这个可以通过输出文件中 Forces (Hartree/Bohr)这一段来确认。 ...

感谢回复。请问这种情况下，如果要求高的话（比如算出来差别还比较大，或者很高精度的计算），有什么解决方法呢？

BLL130

LittlePupil 发表于 2022-9-12 00:02
有点好奇，如果你把最后一帧结构拿出来，chk路径写频率分析的chk文件，然后opt选项写rcfc而不是calcall会是 ...

是怀疑收敛到了不同的波函数吗？SCF Done输出的能量是相同的，感觉应该不是这个问题。

hebrewsnabla

BLL130 发表于 2022-9-13 21:34
感谢回复。请问这种情况下，如果要求高的话（比如算出来差别还比较大，或者很高精度的计算），有什么解决 ...

也许两个任务都加上tight能让差别小一点。

但是原理上应该不会出现差别比较大的情况。本来displacement并不是一个有明确定义的东西，它是依赖于算法的，固然我们可以调整opt任务的iop来使其和freq任务的结果一样，但是我觉得没什么意思。

LittlePupil

BLL130 发表于 2022-9-13 21:34
是怀疑收敛到了不同的波函数吗？SCF Done输出的能量是相同的，感觉应该不是这个问题。

我不是这个意思，只是想通过读freq任务的Hessian矩阵证实2L的猜测。

我赞同楼上“原理上不会出现差别比较大的情况”这一看法，毕竟Gaussian优化算法的稳健性是业内闻名的，不太可能会对displacement的计算出现较大偏差。

BLL130

LittlePupil 发表于 2022-9-13 23:20
我不是这个意思，只是想通过读freq任务的Hessian矩阵证实2L的猜测。

我赞同楼上“原理上不会出现差别 ...

试了一下，是全部收敛的，但是计算结果和calcall的并不相同，神奇。