金属配合物配位键力场参数的拟合的问题

哈哈

老师，想做金属配合物配位键力场参数的拟合，先用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差别大，想问这其中是哪有什么问题吗？麻烦了

sobereva

modified Seminario以及sobtop里我提出的m2Seminario都是直接基于Hessian计算k的，不是拟合出来的，http://sobereva.com/soft/Sobtop有例子

哈哈

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

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