计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 量子化学 (Quantum Chemistry) › IRC后优化的极小点似不在反应坐标上

IRC后优化的极小点似不在反应坐标上

查看数: 8696 | 评论数: 8 | 收藏 Add to favorites 1

关灯 | 提示：支持键盘翻页<-左右->

帖子模式

axuanyan

发布时间: 2020-12-4 05:20

正文摘要:

大家好，我计算完一个过渡金属Lewis酸催化的过渡态并跑了IRC，得到了较平滑的曲线共121点。基于反应坐标两端的结构进一步优化，发现其中一个结构的优化结果和反应坐标中的结构相差很大，不像是沿着反应坐标得到的极 ...

回复 Reply

sobereva 发表于 Post on 2020-12-4 16:14:43

axuanyan 发表于 2020-12-4 15:17
老师，对于150原子左右的IRC，分配28核的资源，几步重算一次力常数比较好？
我现在recalc=20，IRC的情况 ...

跟recalc没有直接关系

axuanyan 发表于 Post on 2020-12-4 14:57:54

本帖最后由 axuanyan 于 2020-12-4 15:19 编辑

FETEYA 发表于 2020-12-4 11:31
IRC的步长可以用StepSize关键词来调
步数越多就越接近极小点,用maxpoints=999后端点结构就和极小点很接 ...

150+原子，没这么大算力。也有可能我recalc的间隔设得太大，我现在是20步。

sobereva 发表于 Post on 2020-12-4 12:31:30

要想对IRC两边的结构优化后，得到的极小点结构相对于IRC端点变化不显著，前提是IRC已经跑得足够充分，比如IRC两侧的曲线已经接近水平、受力也比较小。

另外，有的反应路径本身就很长，还有肩峰，即便满足上述条件，亦有可能继续优化下去还会经历很长一段优化轨迹、对应结构大幅变化。如果怕优化到不是与当前过渡态直接连通的极小点，应当把优化的步长上限设小，最后优化成什么就是什么。

FETEYA 发表于 Post on 2020-12-4 11:31:53

snljty 发表于 2020-12-4 10:55
IRC步长是固定的，步数再多也不会正好跑到极小值点。

IRC的步长可以用StepSize关键词来调
步数越多就越接近极小点,用maxpoints=999后端点结构就和极小点很接近了. 所以我才怀疑题主是用了maxpoints=60步数不够多,不够接近极小点,所以出现两端的结构做优化后变化很大

snljty 发表于 Post on 2020-12-4 10:55:06

FETEYA 发表于 2020-12-4 10:07
IRC不一定会跑到极小点附近,只有步数和步长够大才会跑到极小点.
如果一定要让IRC跑到极小点,可以在IRC关键 ...

IRC步长是固定的，步数再多也不会正好跑到极小值点。

FETEYA 发表于 Post on 2020-12-4 10:07:11

IRC不一定会跑到极小点附近,只有步数和步长够大才会跑到极小点.
如果一定要让IRC跑到极小点,可以在IRC关键词中加maxpoints=999,和recalc=10保证力常数的准确性

CocuLiu 发表于 Post on 2020-12-4 07:59:46

能不能给出具体图示

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register