计算化学公社

标题: IRC＝Cartesian关键字的问题 [打印本页]

作者
Author: Will 时间: 2019-10-31 17:59
标题: IRC＝Cartesian关键字的问题
看到一篇文献说高斯的IRC的默认单位是“amu −1/2 Bohr”，如果加上irc＝cartesian关键词之后就变成Bohr了，那么这两个单位的换算关系是怎样的呢？还有就是我加了这个关键词以后跑的一个IRC感觉没跑完就正常结束了，没加这个关键词的时候能够完整的跑完，这是为什么？

作者
Author: Mikasa 时间: 2019-11-1 05:23
IRC默认是在mass-weighted coordinate下计算的q_mw=sqrt(M)·q_cart，单位自然是amu^-1/2·Bohr，cartesian coordinate下直接就是Bohr。
坐标转换可以参考Chen, Z. Rotation procedure in intrinsic reaction coordinate calculations. Theoretica chimica acta 1989, 75, 481-484.

作者
Author: sobereva 时间: 2019-11-1 06:54
不要写irc＝cartesian，否则根本不是一般意义的IRC

作者
Author: Will 时间: 2019-11-1 09:20

Mikasa 发表于 2019-11-1 05:23
IRC默认是在mass-weighted coordinate下计算的q_mw=sqrt(M)·q_cart，单位自然是amu^-1/2·Bohr，cartesian ...

好的，多谢！

作者
Author: Will 时间: 2019-11-1 09:24

sobereva 发表于 2019-11-1 06:54
不要写irc＝cartesian，否则根本不是一般意义的IRC

社长，我这并不是要算irc，而是要读出他的宽度代一个算tunneling几率的JWKB公式。另外请问ORCA里的NEB和这个irc＝cartesian有什么区别呢？

作者
Author: Will 时间: 2019-11-1 09:30
关于没跑完就自己结束的问题，后来换了一个基组解决了

作者
Author: sobereva 时间: 2019-11-2 08:44

Will 发表于 2019-11-1 09:24
社长，我这并不是要算irc，而是要读出他的宽度代一个算tunneling几率的JWKB公式。另外请问ORCA里的NEB和 ...

算法的差别、NEB的介绍看
过渡态、反应路径的计算方法及相关问题
http://sobereva.com/44

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