计算化学公社

标题: Multiwfn和Gaussian的Hirshfeld電荷計算 [打印本页]

作者
Author: loveworld9966 时间: 2018-7-7 01:48
标题: Multiwfn和Gaussian的Hirshfeld電荷計算
請教各大神,
我自己使用Gaussian設指令pop=Hirshfeld和Multiwfn計算Hirshfeld得到的電荷值不同......
請問是指令或是兩個軟件計算的方式不同嗎？
為什麼會造成差異呢？

感謝各位大神撥空解答

作者
Author: wbn 时间: 2018-7-7 03:03
本帖最后由 wbn 于 2018-7-7 03:06 编辑

如果是post-HF方法或者激发态有没有记得加density关键词

作者
Author: sobereva 时间: 2018-7-7 04:26
Hirshfeld电荷计算需要有原子自由状态的密度，Gaussian里内置的不知什么途径来的，肯定没Multiwfn自带的合理，Multiwfn自带的原子自由状态的球对称化的密度的构建细节看手册附录3。另外，Multiwfn也允许直接调用Gaussian产生指定级别下的原子波函数文件，不同级别下产生的，最终得到的Hirshfeld电荷也会有所不同。

简而言之，甭管Gaussian输出的，就按Multiwfn手册里的示例的做法计算Hirshfeld电荷就行了。

作者
Author: loveworld9966 时间: 2018-7-9 14:27
有無使用density=current的指令得到Gaussian的fch檔後，再丟入Multiwfn好像得到的Hirshfeld電荷也會不同？

作者
Author: sobereva 时间: 2018-7-9 18:21

loveworld9966 发表于 2018-7-9 14:27
有無使用density=current的指令得到Gaussian的fch檔後，再丟入Multiwfn好像得到的Hirshfeld電荷也會不同？

得说明具体用的关键词
用不用这个不会直接影响fch里记录的轨道，原理上不会影响Multiwfn的分析结果

详见
详谈Multiwfn支持的输入文件类型、产生方法以及相互转换
http://sobereva.com/379

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