计算化学公社

标题: Multiwfn输出的Fock矩阵与out文件相差较大 [打印本页]

作者
Author: 铭律凉 时间: 2021-11-3 17:03
标题: Multiwfn输出的Fock矩阵与out文件相差较大
本帖最后由铭律凉于 2021-11-4 14:39 编辑

我按照manual计算Fock矩阵FOCK.txt
但是它与gaussian输出的fock矩阵差别有点大
manual里虽然提到会有些许误差，但我感觉差的还是有点远
请问问题出在哪里？

关于附件：
  ethane.fchk
  FOCK.txt :multiwfn输出的fock矩阵
  echtane2.out :gaussian输出文件，从文本最末向上搜索fock matrix可以找到最后的fock矩阵

作者
Author: snljty 时间: 2021-11-3 17:49
首先注意一下Fock矩阵不是密度矩阵。你的最后一行写的有问题。
我用你的文件重跑了一下，得到结果和Gaussian输出没有明显差别，和你的FOCK.txt明显不一样。用Gaussian读你的波函数重算导出的NBO 47文件也没有明显差别。怀疑是你操作有问题，建议贴上你的操作步骤。

作者
Author: 铭律凉 时间: 2021-11-4 12:51
本帖最后由铭律凉于 2021-11-4 15:01 编辑

snljty 发表于 2021-11-3 17:49
首先注意一下Fock矩阵不是密度矩阵。你的最后一行写的有问题。
我用你的文件重跑了一下，得到结果和Gaussi ...

我最后一行刚才写错了，但是我确实查找的是fock matrix
我查找的是fock matrix，应该没问题吧而且得到的这个FOCK.txt也不是密度矩阵
我确实是100 17 然后输入输出路径，还是一样的结果，FOCK.txt
但是之后我把multiwfn重启了重新读了fchk然后算100 17 结果这次就对上了
我愿意保证我始终是这么操作的

我之后重现了bug，你可以尝试，我以此通过6号check&modify wafefunction中的功能将重叠矩阵打印在屏幕上之后重新回到100 17 后输出的fock矩阵就不对了
不知道是否属于bug还是属于我的违规操作导致的？

作者
Author: sobereva 时间: 2021-11-4 14:58

铭律凉发表于 2021-11-4 12:51
我最后一行刚才写错了，但是我确实索引的是fock matrix
我索引的是fock matrix，应该没问题吧而且得到的这 ...

是个bug，输出重叠矩阵后程序会对重叠矩阵对角化以输出其本征值，之后没复原重叠矩阵。
已经修正这个问题，重新下载Multiwfn即可

不要用“索引”这个词，语义不明

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