各位老师好，我想计算一系列结构相似的药物分子的Hirshfeld电荷，根据电荷大小来判断氧原子等氢键供受体形成氢键的潜力。就是说，如果氧原子A的Hirshfeld电荷为-0.27；氧原子B的Hirshfeld为-0.25，是否可以判断A作为氢键受体能形成更稳定的氢键？这种判断是否正确呢，能否用在正式的文章中？
对于计算细节，我使用Gaussian在B3LYP/6-31G(d,p)下，使用SMD模型，加GD3BJ矫正对小分子进行几何优化，再用Multiwfn中的“Hirshfeld atomic charge”-->" 1 Use build-in sphericalized atomic densities in free-states (more convenient)" 进行计算，精度是否合适？
如果使用“2 Provide wavefunction file of involved elements by yourself or invoke Gaussian to automatically calculate them”调用Gaussian进行计算，用b3lyp/6-311g(d,p) scrf=smd em=gd3bj 这个级别精度够吗？会不会被审稿人质疑呢？
谢谢各位老师！