计算化学公社

标题: 求助:dimer方法计算得到的过渡态比猜想的末态能量低 [打印本页]

作者
Author: zdbikobc 时间: 2025-12-8 19:16
标题: 求助:dimer方法计算得到的过渡态比猜想的末态能量低
sob老师及各位老师好，我在cp2k中用dimer方法计算过渡态，在dimer之前先计算了cl-neb，但没达到收敛，用100+轮scf以后能量最高点作为过渡态初猜结构，用初始和末态生成了DIMER_VECTOR，dimer计算达到四个收敛，虚频也只有一个，但过渡态能量比末态结构低，后续又单独高精度计算了一次单点能依旧如此，现在疑惑是否根据初末态给出了矢量后，dimer寻找到的反应路径依旧可能不对应猜测路径？接下来应该如何找到正确的初末态结构呢？@sobereva

作者
Author: sobereva 时间: 2025-12-8 19:40
按照虚频模式前、后各微调一次结构，然后做几何优化，获得与这个过渡态直接相连的极小点
你当前情况说明得到的过渡态与之前的极小点并不直接相连（除非是之前的极小点计算偶然收敛到了不稳定波函数导致能量被严重高估）

作者
Author: zdbikobc 时间: 2025-12-9 10:56

sobereva 发表于 2025-12-8 19:40
按照虚频模式前、后各微调一次结构，然后做几何优化，获得与这个过渡态直接相连的极小点
你当前情况说明得 ...

多谢sob老师，我对之前优化后的末态结构进行频率计算发现存在两个较大的虚频，可能并没有达到极小点；我在ts结构基础上按照虚频对应的矢量微调可以优化得到极小点

但目前问题是直接相连的极小点并不能描述我想计算的过程（我计算的是脱水过程，微调得到的初态就已经是脱水之后的构型）。这种情况下是因为反应路径不止一个TS吗？想要计算脱水过程应该如何进行TS结构初猜呢，请老师指教

作者
Author: sobereva 时间: 2025-12-9 14:42

zdbikobc 发表于 2025-12-9 10:56
多谢sob老师，我对之前优化后的末态结构进行频率计算发现存在两个较大的虚频，可能并没有达到极小点；我 ...

根据对机理的认识，重新摆过渡态初猜做dimer，或者根据无虚频的正确的极小点结构跑CI-NEB

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