计算化学公社

 找回密码 Forget password
 注册 Register

计算RESP电荷时Gaussian输入文件的关键字IOp求助

查看数: 3182 | 评论数: 2 | 收藏 Add to favorites 0
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2022-3-2 15:01

正文摘要:

本帖最后由 chemzhang 于 2022-3-2 15:10 编辑 (这个帖子半是求助半是经验交流) 我最近在根据一篇文献的报道使用Amber模拟计算氯化胆碱,准备输入文件时需要知道胆碱一价正离子Ch+的RESP电荷以下是参考文献中 ...

回复 Reply

sobereva 发表于 Post on 2022-3-2 22:43:53
不要认为文献里就是合理的

拟合静电势电荷计算时候需要设定分布于范德华表面附近的拟合点,原则上来说,拟合点的密度越大,拟合结果越好。6/42=N设的是MK方式设置拟合点的时候每平方埃内的点为N个。N默认为1,明显偏小,即得到的拟合静电势电荷质量差。6/42=6刻意用来提高拟合点密度以提高拟合静电势电荷的数值层面的质量。要是加这个有害无益的话,下文的脚本里自然不会加上这个
计算RESP原子电荷的超级懒人脚本(一行命令就算出结果)
http://sobereva.com/476http://bbs.keinsci.com/thread-12858-1-1.html

直接基于波函数文件用Multiwfn计算RESP电荷的时候也可以直接手动设拟合点分布情况,如下所示,可见默认设置下等同于IOp(6/42=6)时候Gaussian产生的拟合点分布(注:虽然拟合点密度相同,但Multiwfn产生的拟合点的坐标和Gaussian略有差异,都是合理的,这导致的结果差异可忽略不计)




原理上,结果会随着拟合点密度的增加趋于收敛,默认设6是在精度和耗时上的较好权衡。继续增大对结果没有明显改进,而减小则可能令结果精度打不可忽视的折扣。

另外,Multiwfn手册4.7.8节专门介绍了怎么可视化拟合点分布以及拟合误差,感兴趣可以看看,对理解原理有好处。

评分 Rate

参与人数
Participants 1
eV +4 收起 理由
Reason
chemzhang + 4 谢谢

查看全部评分 View all ratings

手机版 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.

快速回复 返回顶部 返回列表 Return to list