求助几种利用Gaussian比较不同分子之间的亲电反应活性的做法

yzy1995123

Sob老师和各位老师好，我的研究物质（10个不同结构的分子）均发生亲电反应，反应速率不同，我想从不同分子的亲电反应活性解释差异，阅读相关资料和帖子后思考如下，有几种不同的实现方式：
一、先预测反应位点，然后对预测反应位点的下述值进行比较
1. 利用Hrishfeld原子电荷计算简缩福井函数，将预测反应位点的f -- (r)乘全局亲电性Electrophilicity Index ω（softness），得到预测反应位点的局部亲电指数，不同分子之间比较：
分子反应位点的局部亲电指数越高，分子反应活性越高。
2. 比较预测反应位点的静电势值：
表面静电势更小（更负）的反应活性更高。
二、 不预测反应位点，直接比较
1. 直接比较亲电指数ωcubic  ：
比较数值相对大小而非绝对值，亲电指数ωcubic越高，分子反应活性越高。
2.比较不同分子之间的HOMO轨道值 / LUMO-HOMO 能量差：
HOMO能量越高或者LUMO-HOMO 能量差越小，反应活性越高。
请问以上方法哪条可行呢？各个量的计算都是必要的吗，是否有更简便的比较不同分子间的反应活性的方法呢？
希望能得到各位老师的点拨，谢谢！

sobereva

如置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示，求助帖必须在帖子标题明确体现出此帖内容是求助或提问，并清楚、准确反映出帖子具体内容，避免有任何歧义和含糊性，仔细看http://bbs.keinsci.com/thread-9348-1-1.html。我已把你的不恰当标题 “利用Gaussian比较不同分子之间的亲电反应活性” 改了，以后务必注意

sobereva

记住Hirshfeld怎么拼

对于横向对比不同分子亲电反应活性，只需要考虑这些：

平均局部离子化能（ALIE），看下文
使用Multiwfn的定量分子表面分析功能预测反应位点、分析分子间相互作用
http://sobereva.com/159
如果在感兴趣的原子上没有分子表面上的ALIE极值点以对比考察，可以用Multiwfn计算原子局部表面上的ALIE平均值，参看（把下文里的分析对象从静电势改成ALIE即可）
谈谈怎么计算“原子的静电势”
http://sobereva.com/641（http://bbs.keinsci.com/thread-29998-1-1.html）

简缩局部软度s-，看下文
使用Multiwfn超级方便地计算出概念密度泛函理论中定义的各种量
http://sobereva.com/484（http://bbs.keinsci.com/thread-13328-1-1.html）

原子电荷。但是如《原子电荷计算方法的对比》（http://www.whxb.pku.edu.cn/CN/abstract/abstract27818.shtml）里所述，这并不总靠谱。

也可以计算原子局部范德华表面的静电势平均值，看上面的641博文。

yzy1995123

好的，以后注意！谢谢老师十分详尽的回复！