计算化学公社

标题: 不同泛函的波函数信息 [打印本页]

作者
Author: 赵云跳槽 时间: 2018-1-10 17:13
标题: 不同泛函的波函数信息
本帖最后由赵云跳槽于 2018-1-10 17:15 编辑

背景：使用DFT或MP2的fchk（实际使用SCF density或MP2 density，以前文献都是使用MP2的）作为波函数文件导入到DMA（分布式多极分析）中，得到原子的电荷，偶极，四极等信息，通过计算后发现，DFT或MP2的fchk对这些多极信息的影响不敏感，从而计算得到的分子间相互作用变化不大

解决：虽然不敏感，但还是得解释下原因
      对此，本人的解决方案是通过波函数的对应性分析来说明
      选择了两个泛函M11和B3LYP（与金标准比较的一好一坏的两泛函）
      （选取比较的是HOMO-4 ~ LUMO+4）
(, 下载次数 Times of downloads: 53)

问题：1. 我使用的是SCF density或MP2 density，但分析的时候是通过自然轨道或正则轨道，是否合适？
         2. 前面说过 DFT或MP2的fchk对结果不敏感，通过我图片中的数据，能否说明 “在波函数层次上就比较小，使用这些波函数产生的结果自然也小” 这样的结论？

作者
Author: sobereva 时间: 2018-1-10 17:36
1 不理解你的意思
2 如果你是指DFT和MP2的DMA结果相差不大，这容易理解，因为二者都考虑了电子相关作用。在我的比较原子电荷的文章中也体现出DFT和MP2结果差异不大，但是和HF差异非常大。不建议从轨道角度上讨论，本身不好分析，也不容易直接说明问题。（对不同级别的密度做密度差图，或者再取模做个积分之类，更容易考察差异程度）

作者
Author: 赵云跳槽 时间: 2018-1-10 19:27

sobereva 发表于 2018-1-10 17:36
1 不理解你的意思
2 如果你是指DFT和MP2的DMA结果相差不大，这容易理解，因为二者都考虑了电子相关作用。 ...

我再问下，如果做密度差图，对于MP2方法来说，它使用的是fchk中的SCF density还是MP2 density ？

作者
Author: sobereva 时间: 2018-1-10 19:32

赵云跳槽发表于 2018-1-10 19:27
我再问下，如果做密度差图，对于MP2方法来说，它使用的是fchk中的SCF density还是MP2 density ？

如果你是指Multiwfn，载入高斯直接产生的MP2任务的fch，用的是HF密度，因为此时fch里的轨道是HF的

作者
Author: 赵云跳槽 时间: 2018-1-10 21:36

sobereva 发表于 2018-1-10 19:32
如果你是指Multiwfn，载入高斯直接产生的MP2任务的fch，用的是HF密度，因为此时fch里的轨道是HF的

有点迷惑，我加了density关键词，算出来有MP2 density，我使用的就是这块的密度，怎么就扯到HF密度上了

作者
Author: sobereva 时间: 2018-1-11 05:45

赵云跳槽发表于 2018-1-10 21:36
有点迷惑，我加了density关键词，算出来有MP2 density，我使用的就是这块的密度，怎么就扯到HF密度上了

详谈Multiwfn支持的输入文件类型、产生方法以及相互转换
http://sobereva.com/379
这里详细说了，直接载入输入文件，Multiwfn读取的只是轨道，而不是密度矩阵。轨道对应什么波函数，Multiwfn就分析的什么波函数

如果要利用fch里的不同的密度矩阵，需要按本文做法把密度矩阵专门转化成对应的轨道
在Multiwfn中基于fch产生自然轨道的方法与激发态波函数、自旋自然轨道分析实例
http://sobereva.com/403

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