计算化学公社

标题: 请问可以让Multiwfn打印偶极积分到mwfn文件吗？ [打印本页]
作者
Author: Freeman    时间: 2024-6-27 13:36
标题: 请问可以让Multiwfn打印偶极积分到mwfn文件吗？
社长，您好。

我最近写了个波函数分析的程序，用来实现组内的相关idea。因为mwfn文件被Multiwfn支持，且有详细的文档，所以该程序从设计之初就是完全是依靠Multiwfn转换出的mwfn文件做分析的。mwfn文件里不仅有轨道系数和密度矩阵，更重要的是还有波函数分析不可或缺的电子积分，用起来非常方便。

不过现在Multiwfn好像只支持打印重叠积分到mwfn文件。目前我有一项工作，需要手搓一个Foster-Boys局域化。这就需要偶极积分了，可是Multiwfn并不打印。我知道Multiwfn是支持将偶极积分以易阅读的格式输出到屏幕上的，那么可不可以冒昧地请社长增加一项功能，使得Multiwfn可以打印偶极积分到mwfn文件？

非常感谢！


作者
Author: zjxitcc    时间: 2024-6-27 15:24
我提醒一个额外的问题，如果你用不同软件同一计算级别的molden文件->mwfn文件，然后做分析，你可能会看到结果各不相同，那不是数值误差。如果你真的一开始用的是molden，那你需要检查一下。
作者
Author: Freeman    时间: 2024-6-27 15:38
zjxitcc 发表于 2024-6-27 15:24
我提醒一个额外的问题，如果你用不同软件同一计算级别的molden文件->mwfn文件，然后做分析，你可能会看到结 ...

谢谢提醒。molden文件，除了orca赝势导致核电荷数有问题之外，还有什么问题吗？我暂时只转换过fchk成mwfn。
作者
Author: sobereva    时间: 2024-6-28 05:11
(, 下载次数 Times of downloads: 20)

(, 下载次数 Times of downloads: 18)


Multiwfn主功能6里专门有选项输出各种积分，只不过不是输出到mwfn文件里，而是输出到单独的文本文件。
作者
Author: Freeman    时间: 2024-7-2 18:57
sobereva 发表于 2024-6-28 05:11
Multiwfn主功能6里专门有选项输出各种积分，只不过不是输出到mwfn文件里，而是输出到单独的文本 ...

谢谢。可是打印出来的矩阵是这种格式的，
  1. 1 2 3 4 5
  2. 1
  3. 2
  4. 3
  5. 4
  6. 5
  7. 6
  8. 7
  9. 8
  10. 9
  11. 10
  12. 6 7 8 9 10
  13. 6
  14. 7
  15. 8
  16. 9
  17. 10
复制代码

高斯也喜欢用这种格式打印矩阵。我很久以前写过读取这种格式的代码，非常折磨，十分别扭。请问有什么简便方法让机器读取这些矩阵吗？
作者
Author: sobereva    时间: 2024-7-2 19:25
Freeman 发表于 2024-7-2 18:57
谢谢。可是打印出来的矩阵是这种格式的，

高斯也喜欢用这种格式打印矩阵。我很久以前写过读取这种格式 ...

看不懂你的意思
这是Multiwfn直接输出的

  2. **************** Electric dipole moment matrix (X component) ****************
  3.              1             2             3             4             5      
  4.      1  3.75851E-018
  5.      2  8.78327E-019  3.75851E-018
  6.      3 -7.27078E-002 -4.96075E-001  0.00000E+000
  7.      4  3.80983E-051  0.00000E+000  0.00000E+000  3.75851E-018
  8.      5  1.81568E-035  2.27653E-034 -1.56739E-017  0.00000E+000  3.75851E-018
  9.      6  6.28723E-019  2.87015E-018 -4.82583E-001  6.66664E-051  0.00000E+000
  10.      7 -1.00190E-002 -3.29090E-001  0.00000E+000 -6.75710E-034  0.00000E+000
  11.      8 -9.15262E-054  4.91970E-051 -6.75710E-034  1.88497E-018  0.00000E+000
  12.      9  1.67732E-035  0.00000E+000  0.00000E+000  0.00000E+000  1.88497E-018
  13.     10  1.36567E-022  0.00000E+000 -6.73666E-001  0.00000E+000  0.00000E+000
  14.     11  1.26029E-019  2.05615E-018 -2.24555E-001  0.00000E+000  9.34335E-037
  15.     12  1.26029E-019  2.05615E-018 -2.24555E-001 -1.06198E-051 -4.58754E-035
  16.     13 -1.75418E-035 -5.61898E-035  0.00000E+000 -3.88942E-001  0.00000E+000
  17. ...略
复制代码


Multiwfn里的util.f90中有subroutine readmatgau可以直接读这种格式的矩阵，可以借用



欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3