请问可以让Multiwfn打印偶极积分到mwfn文件吗？

Freeman

社长，您好。

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

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

非常感谢！

zjxitcc

我提醒一个额外的问题，如果你用不同软件同一计算级别的molden文件->mwfn文件，然后做分析，你可能会看到结果各不相同，那不是数值误差。如果你真的一开始用的是molden，那你需要检查一下。

Freeman

zjxitcc 发表于 2024-6-27 15:24
我提醒一个额外的问题，如果你用不同软件同一计算级别的molden文件->mwfn文件，然后做分析，你可能会看到结 ...

谢谢提醒。molden文件，除了orca赝势导致核电荷数有问题之外，还有什么问题吗？我暂时只转换过fchk成mwfn。

sobereva

Multiwfn主功能6里专门有选项输出各种积分，只不过不是输出到mwfn文件里，而是输出到单独的文本文件。

Freeman

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

复制代码

高斯也喜欢用这种格式打印矩阵。我很久以前写过读取这种格式的代码，非常折磨，十分别扭。请问有什么简便方法让机器读取这些矩阵吗？

sobereva

Freeman 发表于 2024-7-2 18:57
谢谢。可是打印出来的矩阵是这种格式的，

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

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

1. 
   
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可以直接读这种格式的矩阵，可以借用