计算化学公社

标题: 求助想通过能量分解提取具有LJ物理意义的项做范德华参数拟合 [打印本页]

作者
Author: hfc 时间: 2023-10-13 22:27
标题: 求助想通过能量分解提取具有LJ物理意义的项做范德华参数拟合
目前需要拟合某几个原子类型的LJ参数，想通过EDA能量分解进行提取具有LJ物理意义的项，然后拟合该势能曲线，体系为Mg2+ 和P2O74-之间的相互作用，我的问题是：
1. 使用sobEDA时， C12是应该拟合Ex+Erep, 还是Erep
2. C6项是应该拟合Ec, 还是EDFTc+Edc
3. Eorb 所考虑的电子转移相关效应是否应该可以隐式的考虑在LJ势中如果考虑系数如何确定？这整一项Eorb应该怎么加入到拟合当中？
4. 该体系使用sobEDA分解能量是否可靠，还是需要通过SAPT2

作者
Author: sobereva 时间: 2023-10-14 00:16
如置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示，求助帖必须在帖子标题明确体现出此帖内容是求助或提问，并清楚、准确反映出帖子具体内容，避免有任何歧义和含糊性，仔细看http://bbs.keinsci.com/thread-9348-1-1.html。我已把你的不恰当标题 “范德华参数拟合” 改了，以后务必注意

作者
Author: sobereva 时间: 2023-10-14 00:23
我不知道这么拟合的意义何在，比起直接对总作用能曲线拟合参数并没什么额外优势

作者
Author: hfc 时间: 2023-10-14 09:58
本帖最后由 hfc 于 2023-10-14 10:09 编辑

sobereva 发表于 2023-10-14 00:23
我不知道这么拟合的意义何在，比起直接对总作用能曲线拟合参数并没什么额外优势

好的 sob老师下次一定注意标题的问题！是因为该拟合总能量曲线两者相关性很差 C6参数甚至会出现负值，即并不能用简单的12-6势+点电荷模型就能拟合好。所以我想到了对照LJ势的物理意义（虽然是猜的），进行拟合，希望Sob老师能解决下我上面提出的问题，十分感谢！
附上拟合的图其中绿色为QM计算的势能扫描（拉远Mg2+与PO43-的距离的刚性扫描势能面），蓝色为原始力场参数结果，橙色为调整LJ来拟合QM的结果。

作者
Author: sobereva 时间: 2023-10-14 17:38

hfc 发表于 2023-10-14 09:58
好的 sob老师下次一定注意标题的问题！是因为该拟合总能量曲线两者相关性很差 C6参数甚至会出现负值， ...

可以看看此文的讨论：J. Chem. Theory Comput. 2015, 11, 2609−2618
J. Phys. Chem. Lett. 2014, 5, 2670−2674

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