分子动力学模拟涉及非标残基,考虑到Multiwfn简单易用,因此尝试用其进行了RESP电荷的计算。但我在进行电荷拟合的时候发现,尽管设置了电荷约束,但是被约束区域求出的总电荷与约束值之间在百分位存在误差。如设置电荷约束为“8-12,15-34,38-73 0.0259”,则被约束区域实际总电荷为0.011486,若使用Gaussian的out文件进行拟合则区域总电荷为0.022539。
Gaussian的fch、out文件及电荷约束文件:
(, 下载次数 Times of downloads: 18)