计算化学公社

标题: 计算福井函数预测反应活性位点过程求助 [打印本页]

作者
Author:
Luffy-D-M    时间: 2019-12-12 17:01
标题: 计算福井函数预测反应活性位点过程求助
向大神请教一下三个问题,我是用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文件吗?这里不太懂。


作者
Author:
sobereva    时间: 2019-12-14 11:19
如置顶的新社员必读贴和论坛首页的公告栏所示,求助帖必须在帖子标题明确体现出“提问”、“求助”要素及具体内容(仔细看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/484http://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节写得超级清楚,搞不懂的反复一个字一个字多阅读几遍
作者
Author:
Luffy-D-M    时间: 2019-12-16 10:29
sobereva 发表于 2019-12-14 11:19
如置顶的新社员必读贴和论坛首页的公告栏所示,求助帖必须在帖子标题明确体现出“提问”、“求助”要素及具 ...

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




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3