计算化学公社

 找回密码 Forget password
 注册 Register

过渡态搜索失败,求指导

查看数: 11854 | 评论数: 15 | 收藏 Add to favorites 1
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2021-8-20 15:01

正文摘要:

cut-off :500, k-POINTS: 8*8*8, RMS: 0.01 outmol:QST cycle limit reached QST cycle:                        &nb ...

回复 Reply

计算小白请关照 发表于 Post on 2021-10-11 16:14:17
wzkchem5 发表于 2021-10-11 15:12
不是基态不基态的问题,而是是不是局部极小值点的问题。初态和末态必须都是局部极小值点,结果才是有意义 ...

好的,谢谢前辈耐心指导
wzkchem5 发表于 Post on 2021-10-11 15:12:28
计算小白请关照 发表于 2021-10-11 02:42
我明白你的意思了,如果初始和末态原本就不是基态,可能在这两者之间就没有过渡态了。,如果尝试降低RMS ...

不是基态不基态的问题,而是是不是局部极小值点的问题。初态和末态必须都是局部极小值点,结果才是有意义的。所以必须做频率计算,没有任何替代方法(如果不做频率计算确定没有虚频,文章根本不可能发表,早晚也要做)。你说的步长问题根本就没在点子上
计算小白请关照 发表于 Post on 2021-10-11 09:42:23
wzkchem5 发表于 2021-10-10 23:42
建议算一下频率看看有没有虚频。原则上两个没有虚频的结构之间的单分子反应必然是有势垒的(只不过不一定 ...

我明白你的意思了,如果初始和末态原本就不是基态,可能在这两者之间就没有过渡态了。,如果尝试降低RMS的值,让扫描步长快点,会不会就有可能找到过渡态,之前会不会有可能是扫描速度太慢,在局部最小值之间一直震荡?
wzkchem5 发表于 Post on 2021-10-10 23:42:11
计算小白请关照 发表于 2021-10-10 14:47
反应物和产物都优化过,但是没有计算频率。

建议算一下频率看看有没有虚频。原则上两个没有虚频的结构之间的单分子反应必然是有势垒的(只不过不一定只经过一个势垒),但是如果其中一个结构有虚频,那就有可能是无垒过程,这种情况下你就永远也找不到过渡态
计算小白请关照 发表于 Post on 2021-10-10 21:47:59
wzkchem5 发表于 2021-10-10 15:58
从图看,搞不好这个反应的过渡态本来就不存在。你的反应物和产物都优化过且算过频率吗?确定没有虚频吗?

反应物和产物都优化过,但是没有计算频率。
wzkchem5 发表于 Post on 2021-10-10 15:58:23
本帖最后由 wzkchem5 于 2021-10-10 09:00 编辑
计算小白请关照 发表于 2021-10-10 07:28
还请前辈们在指导一二,谢谢
从图看,搞不好这个反应的过渡态本来就不存在。你的反应物和产物都优化过且算过频率吗?确定没有虚频吗?
计算小白请关照 发表于 Post on 2021-10-10 14:28:58
wzkchem5 发表于 2021-8-22 02:53
我觉得可以先试试把Max. number of QST steps改大

还请前辈们在指导一二,谢谢
计算小白请关照 发表于 Post on 2021-8-22 13:40:07
wzkchem5 发表于 2021-8-22 02:53
我觉得可以先试试把Max. number of QST steps改大

谢谢,我先试试
wzkchem5 发表于 Post on 2021-8-22 02:53:08
计算小白请关照 发表于 2021-8-21 16:30
这是大致的计算参数设置

我觉得可以先试试把Max. number of QST steps改大
计算小白请关照 发表于 Post on 2021-8-21 21:36:42
wzkchem5 发表于 2021-8-21 02:08
检查一下反应过程中波函数有没有突变,比如自旋布居数突然变化。对于成键类型比较复杂的第一过渡周期元素 ...

在Outmol文件中如何看 波函数变化? AS diffuse calculation.castep (331.58 KB, 下载次数 Times of downloads: 2) 。您说的进行TS精修就是这里的TS confirmation吗?
wzkchem5 发表于 Post on 2021-8-21 02:08:33
计算小白请关照 发表于 2021-8-20 16:23
反应是这样子的,就是V原子(灰色的)从初始反应物的位置跳跃到生成物对应的空位上,只有这一条直接的扩 ...

检查一下反应过程中波函数有没有突变,比如自旋布居数突然变化。对于成键类型比较复杂的第一过渡周期元素,常常有不止一个稳定的波函数,如果几何优化中途在两个稳定波函数之间跳跃,就有可能导致不收敛。
这个不是导致不收敛的概率最大的原因,但是是最容易忽略的原因。检查了波函数没问题以后,可以尝试以现在这个没有完全收敛的TS为初始结构,直接做过渡态精修试试

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-19 00:17 , Processed in 0.193040 second(s), 25 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list