Morse势函数拟合的问题

lc123 · 发表于 Post on 2022-3-30 20:18:26

老师您好，我想重现一遍您之前上课时用Origin程序对乙烷C-C键势能曲线拟合的例子，对这里面的一个参数k(力的常数）不太理解，例子如下，请问这里的k初猜2000是从out文件里面读取的吗，如果是的话请问在out文件里面哪里啊，我没有找到这个力常数啊，谢谢老师

wzkchem5 · 发表于 Post on 2022-3-30 21:19:26

这个2000只是一个初始猜测值，需要根据经验提供，不是从文件里读的。如果不知道这个值应该选多少，就应该选不同的值多试几次，分别拟合，看看哪个猜测值拟合得最好。

lc123 · 发表于 Post on 2022-3-31 09:06:06

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

嗯嗯好的，谢谢老师~

lc123 · 发表于 Post on 2022-3-31 11:13:15

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

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

sobereva · 发表于 Post on 2022-3-31 15:22:17

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

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

lc123 · 发表于 Post on 2022-3-31 16:33:05

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

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

sobereva · 发表于 Post on 2022-4-1 00:55:48

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

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

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

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

[辅助/分析程序] Morse势函数拟合的问题