不要认为文献里就是合理的 拟合静电势电荷计算时候需要设定分布于范德华表面附近的拟合点,原则上来说,拟合点的密度越大,拟合结果越好。6/42=N设的是MK方式设置拟合点的时候每平方埃内的点为N个。N默认为1,明显偏小,即得到的拟合静电势电荷质量差。6/42=6刻意用来提高拟合点密度以提高拟合静电势电荷的数值层面的质量。要是加这个有害无益的话,下文的脚本里自然不会加上这个 计算RESP原子电荷的超级懒人脚本(一行命令就算出结果) http://sobereva.com/476(http://bbs.keinsci.com/thread-12858-1-1.html) 直接基于波函数文件用Multiwfn计算RESP电荷的时候也可以直接手动设拟合点分布情况,如下所示,可见默认设置下等同于IOp(6/42=6)时候Gaussian产生的拟合点分布(注:虽然拟合点密度相同,但Multiwfn产生的拟合点的坐标和Gaussian略有差异,都是合理的,这导致的结果差异可忽略不计) ![]() 原理上,结果会随着拟合点密度的增加趋于收敛,默认设6是在精度和耗时上的较好权衡。继续增大对结果没有明显改进,而减小则可能令结果精度打不可忽视的折扣。 另外,Multiwfn手册4.7.8节专门介绍了怎么可视化拟合点分布以及拟合误差,感兴趣可以看看,对理解原理有好处。 |
参与人数Participants 1 | eV +4 | 收起 理由Reason |
---|---|---|
| + 4 | 谢谢 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2025-8-17 11:20 , Processed in 1.262738 second(s), 26 queries , Gzip On.