计算化学公社

标题: 是否可以更改IRC正常结束预设能量变化数值 [打印本页]

作者
Author: sun0615 时间: 2020-5-15 17:34
标题: 是否可以更改IRC正常结束预设能量变化数值
我要计算的反应是 A…C →B…C 的型态，其中A…C是键长约3.5的Van der Waals complex，我找到的TS是 A…C→B…C的TS，C在这个TS结构中只有非常微小的震动，B…C→BC是一个没有barrier的反应，由上述IRC的结果也可看出A…C→B…C的过程中，complex的键长由3.7慢慢缩短为3.0，但是因为IRC的结尾附近，B…C complex能量变化太小就正常结束了，加大stepsize也没用，加太大会导致IRC一开始就出错。我直接拿IRC最后一点B…C的结构跑opt，是会直接变为BC没错。我的问题是

1. IRC结束的预设能量变化定义数值是可以更改的吗? 我看了说明找不到这个关键词。

2. 如果无法用IRC让B…C跑到BC，我在B…C的opt时将maxstep设小一点，再加上NoTrust，如果opt过程中没有出现能量跳动，是可以convince reviewer我找到的TS就是A…C→BC的吗?

谢谢老师

作者
Author: hebrewsnabla 时间: 2020-5-15 19:04
本帖最后由 hebrewsnabla 于 2020-5-15 19:07 编辑

1. 有啊 IRC=Tight 或者VeryTight
2. 个人觉得现在的结果也凑合了，如果B...C到BC的变化很小的话

作者
Author: sobereva 时间: 2020-5-15 23:04

hebrewsnabla 发表于 2020-5-15 19:04
1. 有啊 IRC=Tight 或者VeryTight
2. 个人觉得现在的结果也凑合了，如果B...C到BC的变化很小的话

IRC=Tight、VeryTight设的不是IRC任务结束的判断阈值，而是IRC的校正步的收敛阈值，不是一码事

Gaussian的IRC的结束标准是通过Newton-Raphson法计算位移，看是否小于当前IRC任务设的步长，如果小于就结束，而非看能量变化。所以原理上来说，步长设得越小，越有可能走得越长

作者
Author: elvisng 时间: 2020-5-16 23:20

sobereva 发表于 2020-5-15 23:04
IRC=Tight、VeryTight设的不是IRC任务结束的判断阈值，而是IRC的校正步的收敛阈值，不是一码事

Gaussi ...

「Gaussian的IRC的结束标准是通过Newton-Raphson法计算位移，看是否小于当前IRC任务设的步长，如果小于就结束」

想請問一下, 那麼 stepsize 會同時影響 Newton-Raphson法计算的位移嗎? 如果會, 那麼把 stepsize 調小了也一樣會結束。請問除了這個之外, 還有別的方法嗎? 我有一個 IRC, 中間會經過一段很平坦的過程, 然後能量會突降繼續反應下去。IRC 計算很容易就停在平坦的中段那裡就結束, 沒法到產物。一直在煩惱該怎麼強迫 IRC 繼續跑下去....

作者
Author: sobereva 时间: 2020-5-17 03:38

elvisng 发表于 2020-5-16 23:20
「Gaussian的IRC的结束标准是通过Newton-Raphson法计算位移，看是否小于当前IRC任务设的步长，如果小于就 ...

不会

没有别的办法，这点只要你看过源代码便知。如果想加入别的判断标准得自己该源码重新编译。

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