计算化学公社

标题: Morse势函数拟合的问题 [打印本页]

作者
Author: lc123 时间: 2022-3-30 20:18
标题: Morse势函数拟合的问题
老师您好，我想重现一遍您之前上课时用Origin程序对乙烷C-C键势能曲线拟合的例子，对这里面的一个参数k(力的常数）不太理解，例子如下，请问这里的k初猜2000是从out文件里面读取的吗，如果是的话请问在out文件里面哪里啊，我没有找到这个力常数啊，谢谢老师 (, 下载次数 Times of downloads: 23) (, 下载次数 Times of downloads: 26)

作者
Author: wzkchem5 时间: 2022-3-30 21:19
这个2000只是一个初始猜测值，需要根据经验提供，不是从文件里读的。如果不知道这个值应该选多少，就应该选不同的值多试几次，分别拟合，看看哪个猜测值拟合得最好。

作者
Author: lc123 时间: 2022-3-31 09:06

wzkchem5 发表于 2022-3-30 21:19
这个2000只是一个初始猜测值，需要根据经验提供，不是从文件里读的。如果不知道这个值应该选多少，就应该选 ...

嗯嗯好的，谢谢老师~

作者
Author: lc123 时间: 2022-3-31 11:13

wzkchem5 发表于 2022-3-30 21:19
这个2000只是一个初始猜测值，需要根据经验提供，不是从文件里读的。如果不知道这个值应该选多少，就应该选 ...

再次打扰老师一下了，我想问一下拟合出来的势函数标准误差在多少以内比较合适啊，谢谢老师

作者
Author: sobereva 时间: 2022-3-31 15:22

lc123 发表于 2022-3-31 11:13
再次打扰老师一下了，我想问一下拟合出来的势函数标准误差在多少以内比较合适啊，谢谢老师

把拟合的曲线和原始量子化学扫描出的数据点同时显示出来，肉眼一看就能判断拟合合理程度。基本吻合上就行了，不要教条

作者
Author: lc123 时间: 2022-3-31 16:33

sobereva 发表于 2022-3-31 15:22
把拟合的曲线和原始量子化学扫描出的数据点同时显示出来，肉眼一看就能判断拟合合理程度。基本吻合上就行 ...

嗯嗯好的，谢谢老师的解答。但我在拟合的过程中发现通过扫描然后拟合出来的键解离能，和通过计算A分子断裂后的B+C的焓相减得出的键解离能不一致，请问老师这种现象是正常的吗？如果正常的话，那这两种方法得出的键解离能哪个更准确一点呢？

作者
Author: sobereva 时间: 2022-4-1 00:55

lc123 发表于 2022-3-31 16:33
嗯嗯好的，谢谢老师的解答。但我在拟合的过程中发现通过扫描然后拟合出来的键解离能，和通过计算A分子断 ...

如果你用的理论方法有大小一致性（比如一般用的DFT泛函、微扰、耦合簇等），并且扫描方式得当（比如没有做出化学键异裂的扫描却用RKS来算这种事），两种方式得到的电子能量变化完全相同

扫描只能得到解离过程的电子能量变化，和解离焓没有可比性

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