计算化学公社

标题: fukui函数测亲电性问题 [打印本页]

作者
Author: miya 时间: 2016-5-20 17:04
标题: fukui函数测亲电性问题
sob老师，我想再请教下关于fukui函数测算亲电性，亲核性的具体。我已经可以测算出N,N+1，N-1分子对应的每个原子的Hirshfeld电荷，但是我不清楚您给的图中f-，f+底下对应的每个原子的Hirshfeld电荷是怎么来的（下图红色标出的），我按照浮亏函数直接先测算亲电性，让后返回操作界面后直接输入 7 // Population analysis，1 // Hirshfeld population，但是得出的结果显示的原子电荷值跟单独的N是一样的，请问老师到底应该怎么操作啊

作者
Author: sobereva 时间: 2016-5-20 23:04
看手册4.7.3节的例子

作者
Author: miya 时间: 2016-5-29 17:05

sobereva 发表于 2016-5-20 23:04
看手册4.7.3节的例子

又打搅了sob老师，我看了，但是没有写出怎么得到对应的具体原子N-1，N+1状态下的Hirshfeld电荷，我也翻了下4.5.4节，也没有发现额，求老师可以具体一点么，我卡在这块了

作者
Author: sobereva 时间: 2016-5-29 18:26

miya 发表于 2016-5-29 17:05
又打搅了sob老师，我看了，但是没有写出怎么得到对应的具体原子N-1，N+1状态下的Hirshfeld电荷，我也翻了 ...

你载入N-1和N+1状态的.wfn文件，按照手册4.7节的例子就算出电荷来了，超级容易
会算中性的不可能不会算离子态的，操作过程完全一样，仅输入文件不一样而已

作者
Author: miya 时间: 2016-5-30 15:22

sobereva 发表于 2016-5-29 18:26
你载入N-1和N+1状态的.wfn文件，按照手册4.7节的例子就算出电荷来了，超级容易
会算中性的不可能不会 ...

谢谢sob老师，搞定啦

作者
Author: QuantumHan 时间: 2018-4-24 12:06
借楼主问的Fukui函数预测反应亲电性的问题再请教下Sob老师。博文http://sobereva.com/402这段有些迷糊。。先在Gaussian中在B3LYP/6-31G*下优化中性状态的晕苯，将得到的chk文件转换为fch后命名为N.fch。然后基于优化的中性的晕苯的结构，将电荷和自旋多重度改为1 2（即+1阳离子态），算个单点任务（切勿做优化任务），将得到的chk转换并命名为N-1.fch
单点计算N-1的时候基组可以与输入文件为N时的一致么，只是不优化。比如输入文件为N时用的#p m062X/6-311+G(d,p) opt freq N-1时直接写 #p m062x/6-311+G(d,p)，还是得算单点能提高一个级别才可行呢？谢谢老师

作者
Author: sobereva 时间: 2018-4-24 14:19

QuantumHan 发表于 2018-4-24 12:06
借楼主问的Fukui函数预测反应亲电性的问题再请教下Sob老师。博文http://sobereva.com/402这段有些迷糊。。 ...

需要一致

作者
Author: QuantumHan 时间: 2018-4-24 15:16

sobereva 发表于 2018-4-24 14:19
需要一致

老师不好意思，我没表述清楚。我是想问，优化用的#p m062X/6-311+G(d,p) opt freq . 将优化后的N和N-1两个文件保存成两个单点输入文件，这两个计算单点能否还是仍用#p m062X/6-311+G(d,p)呢，还是需要提升一个级别比如用双杂化泛函？因为在算势能面能垒的时候用的后HF方法，而后HF方法计算又显得比较昂贵，只是想得到这个N-1的.chk文件作分析

作者
Author: sobereva 时间: 2018-4-24 17:09

QuantumHan 发表于 2018-4-24 15:16
老师不好意思，我没表述清楚。我是想问，优化用的#p m062X/6-311+G(d,p) opt freq . 将优化后的N和N-1两 ...

#p m062X/6-311+G(d,p) opt freq优化N，直接就得到了N.fch
然后再用N.fch保存成新的.gjf文件，设成N-1状态的单点任务，还是用m062X/6-311+G(d,p)，算完了就有了N-1.fch。通过这俩fch求差得到福井函数。
没必要用后HF。电子密度相关的量对计算级别敏感度远远低于能量。

作者
Author: QuantumHan 时间: 2018-4-27 20:44

sobereva 发表于 2018-4-24 17:09
#p m062X/6-311+G(d,p) opt freq优化N，直接就得到了N.fch
然后再用N.fch保存成新的.gjf文件，设成N-1 ...

谢谢老师的指导。

作者
Author: liuyjhx 时间: 2019-12-21 14:32

sobereva 发表于 2018-4-24 14:19
需要一致

老师为了更准确我可以用b3lyp/genecp优化出来的结构，然后把N，N-1状态下都在b2plypd3/def2tzvp 基组下计算单点能来考查亲电位点吗？

作者
Author: sobereva 时间: 2019-12-22 00:36

liuyjhx 发表于 2019-12-21 14:32
老师为了更准确我可以用b3lyp/genecp优化出来的结构，然后把N，N-1状态下都在b2plypd3/def2tzvp 基组下计 ...

没意义
算电子能量更准的不代表算的波函数也更准，何况波函数分析对波函数质量要求并不高，用普通泛函足够了。而且本来Gaussian算双杂化就很慢，产生双杂化的密度更耗时

欢迎光临计算化学公社 (http://bbs.keinsci.com/)