计算化学公社

标题: 含自由H结构的几何优化不收敛问题 [打印本页]

作者
Author: Iris_1 时间: 2023-6-11 12:04
标题: 含自由H结构的几何优化不收敛问题
本帖最后由 Iris_1 于 2023-6-11 12:08 编辑

需要用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
*

作者
Author: sobereva 时间: 2023-6-11 21:39
拿GAFF优化断键的状态完全莫名其妙，GAFF根本不可能靠谱地优化自由基体系。不要乱用分子力场
而且本身这种过程就是没势垒的

作者
Author: Iris_1 时间: 2023-6-11 21:56

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

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

作者
Author: sobereva 时间: 2023-6-12 10:41

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

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

作者
Author: Iris_1 时间: 2023-6-12 19:21

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

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

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