对于opt任务，fchk文件存的坐标是优化以后的吗？

让你变成回忆 · 发表于 Post on 2018-6-2 11:15:29

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

zjxitcc · 发表于 Post on 2018-6-2 12:24:15

本帖最后由 zjxitcc 于 2018-6-2 12:26 编辑

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

sobereva · 发表于 Post on 2018-6-2 12:28:56

是优化以后的

让你变成回忆 · 发表于 Post on 2018-6-2 14:00:07

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

谢谢~~

让你变成回忆 · 发表于 Post on 2018-6-2 14:00:38

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

谢谢~

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

[Gaussian/gview] 对于opt任务，fchk文件存的坐标是优化以后的吗？

浏览过的版块