计算化学公社

标题: 对于opt任务，fchk文件存的坐标是优化以后的吗？ [打印本页]

作者
Author: 让你变成回忆 时间: 2018-6-2 11:15
标题: 对于opt任务，fchk文件存的坐标是优化以后的吗？
如题，最近在写程序遇到的困惑。基于fch文件，然后从中读取关键词”Current cartesian coordinates”，然后把单位从Bohr转换成了angstrom，这样得到了一个坐标；
另一种方式是把优化结束以后的log文件用Gview打开，然后重新保存为输入文件。这样又得到一个坐标；
我对比两种方式得到的坐标，在小数点后面5位或者6位以后不一样，请问这样对计算结果有影响吗？我用Multiwfn也试了下，产生的坐标也是在小数点后5位与log文件有差别。

作者
Author: zjxitcc 时间: 2018-6-2 12:24
本帖最后由 zjxitcc 于 2018-6-2 12:26 编辑

“有无影响”是个简单粗暴的二分概念，并不适合这么问。应该说“影响有多大”。坐标小数点后5位有差别时，算个单点，能量可能会在小数点后第7、8位有区别，此时讨论这种问题没有意义，因为计算数值误差等非重现性因素已经达到同样误差水平，当你换算到kcal/mol单位时，发现是看不出区别的。
况且，当前最准的电子结构方法，预测键长也不过是在0.001A级别，更后面的小数已经不是有效数字了。
另外，你从Bohr转换为Angstrom，小数点后第6、7位也不一定准，其中还有误差呢

作者
Author: sobereva 时间: 2018-6-2 12:28
是优化以后的

作者
Author: 让你变成回忆 时间: 2018-6-2 14:00

zjxitcc 发表于 2018-6-2 12:24
“有无影响”是个简单粗暴的二分概念，并不适合这么问。应该说“影响有多大”。坐标小数点后5位有差别时， ...

谢谢~~

作者
Author: 让你变成回忆 时间: 2018-6-2 14:00

sobereva 发表于 2018-6-2 12:28
是优化以后的

谢谢~

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