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

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

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

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

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

谢谢老师的指导。

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。电子密度相关的量对计算级别敏感度远远低于能量。

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文件作分析

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

需要一致

借楼主问的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)，还是得算单点能提高一个级别才可行呢？谢谢老师

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

谢谢sob老师，搞定啦

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

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

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

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

看手册4.7.3节的例子