含自由H结构的几何优化不收敛问题

Iris_1

需要用ORCA计算2-fluoro-2,2-dinitroethanol分子失去－ＯＨ基团中的H原子的能垒，准备用NEBTS方法。在准备产物结构时，用Aｖｏｇａｄｒｏ软件断开H键，并用GAFF力场优化了一下，用力场优化后的结构作为初始结构进行OPT　FREQ计算，结果提示优化圈数已经达到上限，优化不收敛，而且最后一步的ｘｙｚ文件H原子又与O原子相连了。

ERROR !!!       The optimization did not converge but reached the maximum number of optimization cycles.
       As a subsequent Frequencies calculation has been requested ORCA will abort at this point of the run.
       Please restart the calculation with the lowest energy geometry and/or a larger maxiter for the geometry optimization.
猜测可能是初始的优化结构不合理，求问H原子的位置如何摆放使得初始结构更合理；
如果不是初始结构的问题，该如何解决呢？谢谢解答～

附件为结构优化的初始结构图及ｘｙｚ文件
ｉｎｐ文件为：
! B3LYP D4 def2-TZVP def2/J RIJCOSX OPT FREQ verytightSCF miniprint nopop
%maxcore 1500
%pal nprocs 4 end

*xyz 0 1
C          1.68175        0.39929       -0.28502
C          0.65916        1.36320        0.36890
N          0.86953        1.39450        1.90569
O          1.96280        1.73657        2.32801
O         -0.04559        0.98693        2.60152
N          0.86054        2.83351       -0.11093
O          0.07265        3.67243        0.29718
O          1.85110        3.09749       -0.77420
F         -0.56147        1.09710       -0.14435
O          1.43952       -0.95282        0.11553
H          2.69706        0.68335        0.00353
H          1.59805        0.46593       -1.37271
H          3.83933       -1.71013        2.81878
*

sobereva

拿GAFF优化断键的状态完全莫名其妙，GAFF根本不可能靠谱地优化自由基体系。不要乱用分子力场
而且本身这种过程就是没势垒的

Iris_1

sobereva 发表于 2023-6-11 21:39
拿GAFF优化断键的状态完全莫名其妙，GAFF根本不可能靠谱地优化自由基体系。不要乱用分子力场
而且本身这种 ...

谢谢sob老师的回复！冒昧再请教一下您，请问这种含自由基的结构如何优化才能收敛呢，以及您说的H键断裂没有势垒能否具体说一下（为什么是没有势垒的呢）。刚开始接触量化知识面还比较窄，感谢您的指导！

sobereva

Iris_1 发表于 2023-6-11 21:56
谢谢sob老师的回复！冒昧再请教一下您，请问这种含自由基的结构如何优化才能收敛呢，以及您说的H键断裂没 ...

一律用量子化学方法优化
只有一个过程中同时有能量升高和降低因素才可能有势垒。简单的O-H键断裂几乎没有能量降低因素

Iris_1

sobereva 发表于 2023-6-12 10:41
一律用量子化学方法优化
只有一个过程中同时有能量升高和降低因素才可能有势垒。简单的O-H键断裂几乎没 ...

明白了，谢谢sob老师指导！