计算化学公社

标题: 请问计算MECP的计算方法，应该与几何优化方法一致，还是与单点能计算的计算方法一致？ [打印本页]

作者
Author: HexSona 时间: 2025-9-14 22:28
标题: 请问计算MECP的计算方法，应该与几何优化方法一致，还是与单点能计算的计算方法一致？
各位老师好，本人正在用sobMECP计算最小能量交叉点，如图所示，Int3→MECP→Int1（几何优化用的是PBE0/def2SVP；单点能用的是M06-2X/def2TZVP）。但考虑到MECP点一般不算频率，也就无法进行热力学校正。那么如图所示，为了对比MECP与另外两个自旋多重态的能量差，计算MECP时应该选用什么计算方法？PBE0/def2SVP？还是M06-2X/def2TZVP

作者
Author: sobereva 时间: 2025-9-16 03:59
TPSSh/def-TZVP下做极小点优化和MECP优化，直接取电子能量对比

作者
Author: HexSona 时间: 2025-9-16 11:01

sobereva 发表于 2025-9-16 03:59
TPSSh/def-TZVP下做极小点优化和MECP优化，直接取电子能量对比

好的，谢谢sob老师。请问是单独对Int3、Int1这两个中间体用TPSSh/defTZVP做几何优化，其他结构仍用我之前的计算方法（PBE0/def2SVP）计算，是吗？
对比其他结构与Int3、Int1的自由能时，用PBE0/def2SVP与M06-2X/def2TZVP的热力学校正后的数据；而对比Int3、Int1与MECP的能量时，直接用TPSSh/defTZVP的电子能量，是这个意思吗？

作者
Author: sobereva 时间: 2025-9-16 21:04

HexSona 发表于 2025-9-16 11:01
好的，谢谢sob老师。请问是单独对Int3、Int1这两个中间体用TPSSh/defTZVP做几何优化，其他结构仍用我之前 ...

所有驻点和MECP的优化方法保持统一。用PBE0也行。但def2-SVP对于计算能量差来说太糙

是

作者
Author: HexSona 时间: 2025-9-16 22:34

sobereva 发表于 2025-9-16 21:04
所有驻点和MECP的优化方法保持统一。用PBE0也行。但def2-SVP对于计算能量差来说太糙

是

好的，明白了。非常感谢sob老师的耐心解答

作者
Author: HexSona 时间: 2025-9-19 16:38

sobereva 发表于 2025-9-16 21:04
所有驻点和MECP的优化方法保持统一。用PBE0也行。但def2-SVP对于计算能量差来说太糙

是

sob老师您好，我又来提问了。若MECP计算过程迟迟不收敛，能否适当放宽收敛标准？
我在上述结构上加入了阴离子（PF6-，六氟磷酸根离子）后寻找MECP，但是收敛难度显著增大，主要是Max Gradient El.一项迟迟不收敛，其他四项均已收敛。
我从中提取了最接近收敛标准的一帧结构，并将MECP.f里的STPMX = 0.1d0改为0.03d0后重新搜索MECP，但仍然迟迟无法收敛，甚至距离收敛标准越来越远。其中最接近收敛标准的为：

Step Number 1 -- MECP not yet converged

Max Gradient El.: 0.000908 (0.000700) NO
RMS Gradient El.: 0.000203 (0.000500)  YES
Max Change of X: 0.001362 (0.004000)  YES
RMS Change of X: 0.000353 (0.002500)  YES
Difference in E: 0.000004 (0.000050)  YES

可否将该帧结构用于发表？

以下是输出文件，记载了每一帧的收敛信息

作者
Author: sobereva 时间: 2025-9-20 02:44

HexSona 发表于 2025-9-19 16:38
sob老师您好，我又来提问了。若MECP计算过程迟迟不收敛，能否适当放宽收敛标准？
我在上述结构上加入了 ...

这收敛得就算可以了

作者
Author: HexSona 时间: 2025-9-20 14:41

sobereva 发表于 2025-9-20 02:44
这收敛得就算可以了

好嘞，谢谢sob老师的指点

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