计算福井函数预测反应活性位点过程求助

Luffy-D-M

向大神请教一下三个问题，我是用LUX系统下Gaussian做的N,N-,N+1的优化和单点任务，得出了.chk文件后转化为.fchk文件，用Multiwfn计算的Hrishfeld原子电荷，再后续计算福井函数时遇到一些问题，故向各位大佬请教，不胜感激！！1.附件中的OPT-G,OPT-G1,OPT-G2.fchk文件分别表示q（N+1）净电荷数为-1的原子电荷，q（N）净电荷数为0的原子电荷，q（N-1）净电荷数为+1的原子电荷，为什么我用福井函数公式f+=q(N+1)-q(N)，f0=(f++f-)/2，f-=q(N)-q(N-1)算的f0,f-，f+都为负值？2.用福井函数怎么表述自由基反应的反应活性位点，是f0越趋于0越容易发生自由基反应呢？3.用Multiwfn算福井函数时遇到错误，当我出入完一个中性分子的.fchk文件后-,examples\phenol_N-1.wfn  // “-“ is subtraction sign. Property of the firstly loaded file (namely phenol.wfn) will be subtracted by corresponding property of phenol_N-1.wfn到这里是要输入一个N-1的.fchk文件吗？这里不太懂。

sobereva

如置顶的新社员必读贴和论坛首页的公告栏所示，求助帖必须在帖子标题明确体现出“提问”、“求助”要素及具体内容（仔细看http://bbs.keinsci.com/thread-9348-1-1.html），我已把你的不恰当标题“福井函数预测反应活性位点”改了，以后务必注意。

1 光是第一步就错了。计算福井函数根本就不能对N+1和N-1态单独优化，这在手册4.5.4节，以及此文里都明确强调了
使用Multiwfn作电子密度差图
http://sobereva.com/113

缺乏基本知识的话直接用此文的做法，又方便还省得搞错
使用Multiwfn超级方便地计算出概念密度泛函理论中定义的各种量
http://sobereva.com/484（http://bbs.keinsci.com/thread-13328-1-1.html）

2 f0越大越容易自由基反应

3 有什么错误？
手册4.5.4节里已经说明了N-1、phenol_N-1.wfn是什么含义。手册里说的是wfn，为什么你要提fchk文件？虽然此处用fchk文件也没问题
手册4.5.4节写得超级清楚，搞不懂的反复一个字一个字多阅读几遍

Luffy-D-M

sobereva 发表于 2019-12-14 11:19
如置顶的新社员必读贴和论坛首页的公告栏所示，求助帖必须在帖子标题明确体现出“提问”、“求助”要素及具 ...

谢谢老师指点！还是得认真看文章。看了那篇文章后，照着重新用Multiwfn算了一遍，都算出来了。可是最后的f0的值有点差强人意，和我的结果不是很一致，所以现在准备用ALIE再算一下，看看结果如何。谢谢老师！！