关于简缩福井函数与Hirshfeld原子电荷

傲娇的小螃蟹

现在需要计算简缩福井函数来预测分子上亲核取代的位点，看了相关内容后，还是不清楚怎样做，想请教一下：1，计算Hirshfeld电荷时优化N,N+1,N-1分子时，其构型需要保持一致吗，优化的方法和基组是根据具体的分子进行选择确定的吗？2，怎样用计算出 N,N+1，N-1分子对应的每个原子的Hirshfeld电荷的，是在Gaussian中进行计算吗？3，简缩福井函数是怎样基于Hirshfeld电荷进行预测亲电(亲核)性的。还请各位老师，大神们指教。

sobereva

1 常规方式优化中性状态(N)，基于其结构，计算N、N+1、N-1状态的原子电荷，代入公式。算N+1和N-1时不要优化，直接用N的结构计算单点。简缩福井函数计算过程在手册4.7.3节写明了。

2 Hirshfeld电荷显然用Multwfn计算，过程看Multiwfn手册4.7.1节的例子。

基本知识看
《亲电取代反应中活性位点预测方法的比较》http://www.whxb.pku.edu.cn/CN/abstract/abstract28694.shtml
《预测亲核反应位点方法的比较》https://link.springer.com/article/10.1007%2Fs11426-015-5494-7
概念密度泛函综述和重要文献合集：http://bbs.keinsci.com/forum.php?mod=viewthread&tid=384
幻灯片：Predicting reactive sites（http://sobereva.com/234）

傲娇的小螃蟹

sobereva 发表于 2018-5-15 20:25
1 常规方式优化中性状态(N)，基于其结构，计算N、N+1、N-1状态的原子电荷，代入公式。算N+1和N-1时不要优化 ...

谢谢sob大神的指点