求助用高斯09在PM6D3下算电荷报错

yjb

输入文件：
%chk=input
%nprocshared=36
%mem=20GB
#p PM6D3 NoSymm Force Punch=Derivatives Charge Prop=(Field,Read)

报错：
Missing nuc terms in AM1 ESP.
Error termination via Lnk1e in /share/home/yjb/soft/g09/l602.exe at Mon Mar  7 23:48:21 2022.

本人目前在使用gromacs+Gaussian 跑QM/MM MD的蛋白体系模拟。用DFT或者HF并没有这个问题，上网查了这个问题发现很少提到。
故请教各位老师

Best regard！
Thank you very much！

wzkchem5

首先，为什么要用PM6算电荷？PM6本身就不是为算电荷设计的，即使能算也不能保证结果可靠

sobereva

帖子标题别随便用叹号，极富攻击性，这在置顶的新社员必读贴里明确说了。给你把叹号去掉了

程序不支持半经验方法算静电势而已。半经验方法产生的静电势质量很差，就算程序能算出来结果没什么实际意义，好在程序报错了，免得你用很烂的电荷做QM/MM。做QM/MM的人一定要对原子电荷的特点、原子电荷对静电势重现性、波函数的档次等问题有基本的了解，否则QM和MM部分静电作用根本没法合理描述。原子电荷对静电势重现性问题看《原子电荷计算方法的对比》（http://www.whxb.pku.edu.cn/CN/abstract/abstract27818.shtml）里的专门讨论

yjb

sobereva 发表于 2022-3-8 22:05
帖子标题别随便用叹号，极富攻击性，这在置顶的新社员必读贴里明确说了。给你把叹号去掉了

程序不支持半 ...

谢谢老师， 我去好好看看。