计算化学公社

标题: 关于IRC相关问题的请教 [打印本页]

作者
Author: 978142355 时间: 2016-5-24 09:27
标题: 关于IRC相关问题的请教
依然是那本exploring chemistry with electronic structure methods (second edition)。再看到Si+与硅烷的过渡态那个例子的时候（200页），里面有一个IRC=(RCFC,StepSize=30,MaxPoints=15) scf=qc。如图1所示。我的问题是：1.StepSize=30是设定歩幅大小，那么其单位是什么，高斯手册上是一个样，而这本书又是另外一个样，图1是这本书的，图2是高斯手册的。（PS：其一，这个单位怎么推导过来；其二，在这本书中如果默认是0.1amu-1bohr，那么设定StepSize=30是不是应该是3amu-1bohr了？）
2.我自己通过这个例子，算了水过渡态的IRC曲线，其中图3是输入文件，图4为输出文件用GV显示的图像。我设定了MaxPoints=21，不是应该总共有42个点才对吗，为何31个点就计算结束了？
在此先谢谢大家。

作者
Author: 陌上芬菲 时间: 2016-5-24 12:09
如果一边已到达它认为的平衡态了就自动停了，（可以结合RMS-force看）。步长，不了解，坐等高手。

作者
Author: 978142355 时间: 2016-5-24 14:04
本帖最后由 978142355 于 2016-5-24 19:35 编辑

陌上芬菲发表于 2016-5-24 12:09
如果一边已到达它认为的平衡态了就自动停了，（可以结合RMS-force看）。步长，不了解，坐等高手。

恩，谢谢你的解答，谢谢讨论

。十分感谢你的回答，撒花。

作者
Author: sobereva 时间: 2016-5-24 15:49
stepsize写正数是0.01 bohr为单位，写负数是0.01 amu^1/2 Bohr为单位。以当前版本手册为准。

作者
Author: 978142355 时间: 2016-5-24 19:34
本帖最后由 978142355 于 2016-5-24 20:06 编辑

sobereva 发表于 2016-5-24 15:49
stepsize写正数是0.01 bohr为单位，写负数是0.01 amu^1/2 Bohr为单位。以当前版本手册为准。

有老师这句话我就放心了，同时我也理解为什么会31歩就结束了。谢谢sob老师，也谢谢@陌上芬菲，谢谢交流讨论。不过sob老师为什么stepsize=N设为负值，为何单位就变了呢，其单位是怎么过来的；amu1/2-Bohr，中间那个减号是什么意思，我查了一下，说乘号的有，说除号的也有；以及什么情况下设置为负值（PS：我看见过的都是设置成正值的情况）？问题有些多，基础比较薄弱，还请sob老师别见怪。

作者
Author: sobereva 时间: 2016-5-24 22:16

978142355 发表于 2016-5-24 19:34
有老师这句话我就放心了，同时我也理解为什么会31歩就结束了。谢谢sob老师，也谢谢@陌上芬菲，谢谢交流讨 ...

IRC是在质权坐标下定义的，所以如果你在笛卡尔坐标下定义步长，程序也会自动转化成质权坐标下的步长
那个“减号”是相乘

作者
Author: 978142355 时间: 2016-5-25 09:40

sobereva 发表于 2016-5-24 22:16
IRC是在质权坐标下定义的，所以如果你在笛卡尔坐标下定义步长，程序也会自动转化成质权坐标下的步长
那 ...

原来如此，明白了，谢谢sob老师的指点

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