计算化学公社

标题: 金属配合物配位键力场参数的拟合的问题 [打印本页]

作者
Author: 哈哈 时间: 2022-8-12 10:41
标题: 金属配合物配位键力场参数的拟合的问题
老师，想做金属配合物配位键力场参数的拟合，先用mseminario得到harmonic进行了拟合，k为58.431693，r为1.98，其拟合的曲线如图1。想采用morse对配位键进行拟合，拟合的曲线为图2，通过拟合QM该键的扫描得到键的morse表达式，其表达式为：y=75.3637*{1-exp[-1.4463*(r-1.98)]}^2。通过这个表达式求得的k=315.289117。morse的D和harmonic的k存在一个转化关系a=(k/2D )^(1/2),看由D转化的k和mseminario计算的k差别大，想问这其中是哪有什么问题吗？麻烦了

作者
Author: sobereva 时间: 2022-8-12 10:49
modified Seminario以及sobtop里我提出的m2Seminario都是直接基于Hessian计算k的，不是拟合出来的，http://sobereva.com/soft/Sobtop有例子

作者
Author: 哈哈 时间: 2022-8-17 16:23

sobereva 发表于 2022-8-12 10:49
modified Seminario以及sobtop里我提出的m2Seminario都是直接基于Hessian计算k的，不是拟合出来的，http:// ...

是的老师，mseminario是计算出的。我以为说拟合的morse的D推出的k会和mseminario是计算出的k接近，没想到差别那么大，所以有疑问。我那张拟合的harmonic图是为了为了看计算的k的合理性做的。

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