计算化学公社

标题: 求助：使用multiwfn进行简缩福井函数计算，请大家帮忙看看步骤和输入文件是否有错 [打印本页]

作者
Author: 嘟嘟嘟嘟 时间: 2021-5-24 21:31
标题: 求助：使用multiwfn进行简缩福井函数计算，请大家帮忙看看步骤和输入文件是否有错
本帖最后由嘟嘟嘟嘟于 2021-5-24 21:31 编辑

本人计算化学小白，摸索了一下计算简缩福井函数的步骤，此前已经参照社长的博文【使用Multiwfn超级方便地计算出概念密度泛函理论中定义的各种量http://bbs.keinsci.com/thread-13328-1-1.html】计算出了苯酚的各种量，现在参考文献计算了新的分子，发现结果差距较大，如F原子，N-1态分子，我计算出来的hirshfeld电荷为-0.0514，而文献中的值为-0.0857，请大家帮忙看下我计算的步骤是否有错，还有就是我选取的参数是否漏掉了文献中的一些东西，文献中的参数及计算结果见图。
计算步骤如下：
1.chemical draw 画出环丙沙星分子结构式，导出.gjf ,用xtb优化结构
2.g16优化结构，并计算频率，关键词：# opt freq b3lyp/6-311g(d,p)，结尾收敛判断标准四个YES
3.将第2步得到的.out文件拉入Gaussianview中，得到新的.gjf文件，导入multiwfn中得到不同状态的计算文件，关键词：#P B3LYP/6-311G** out=wfn nosymm
4.按照博文使用Multiwfn超级方便地计算出概念密度泛函理论中定义的各种量http://bbs.keinsci.com/thread-13328-1-1.html中的方法提交到multiwfn中计算简缩福井函数5.所有文件放在压缩包中，请大家帮忙看下，谢谢。

作者
Author: sobereva 时间: 2021-5-25 00:24
你的原子序号顺序和别人文章里都不一样，没法对比。另外注意构象也会影响计算结果，要和别人文章对比需要保证构象相同。

光从你的描述上看，你的计算过程并没有错误。如果你能精确重复出http://sobereva.com/484博文里的例子，就用完全相同的过程算你自己的体系就完了，算出来是什么结果就是什么结果。

只有GaussView没有Gaussianview

作者
Author: 嘟嘟嘟嘟 时间: 2021-5-25 09:18
感谢社长

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