Multiwfn保存的molden支持h函数，请问是如何排序的？

beefly

molden是常见的波函数据格式，但是一直不支持h和更高角动量的基函数（很多年前有人贡献过支持h函数的代码，但是未被作者采用）。如果仅仅为了数据传递而不用于显示，其实可以自己定义一套高角动量的基函数排序方式。最近，发现Multiwfn保存的molden文件可以使用h函数，但不知道直角h函数和球h函数如何排序。

zjxitcc

我也对这个问题挺好奇的，是跟哪个软件保持一致么还是？molden的H问题偶尔会是个令人头疼的问题

niobium

两位z老师可以看一下源代码啊

sobereva

我今天刚更新了官网上的Multiwfn最新版本，下面我针对当前的最新版本来说。

Molden格式没有定义标准的h的顺序，所以Multiwfn在导出和载入molden文件的时候都是用的Multiwfn的内部顺序，对于h的笛卡尔型是
"ZZZZZ","YZZZZ","YYZZZ","YYYZZ","YYYYZ","YYYYY","XZZZZ","XYZZZ","XYYZZ","XYYYZ","XYYYY","XXZZZ","XXYZZ","XXYYZ","XXYYY","XXXZZ","XXXYZ","XXXYY","XXXXZ","XXXXY","XXXXX"
这和fch文件、wfx文件完全一样

对于球谐型是
"H 0  ","H+1  ","H-1  ","H+2  ","H-2  ","H+3  ","H-3  ","H+4  ","H-4  ","H+5  ","H-5  "
这和fch文件完全一样

顺带一提，ORCA也能产生带h的molden文件，但是不知道是什么顺序，经测试带h时无法通过Multiwfn的测试。

由于molden格式不合理、不理想、不严格的地方很多，所以我之前提出了个mwfn格式（https://doi.org/10.26434/chemrxiv.11872524），希望能被业内支持，最终取代molden格式。

beefly

sobereva 发表于 2020-6-2 18:11
我今天刚更新了官网上的Multiwfn最新版本，下面我针对当前的最新版本来说。

Molden格式没有定义标准的h ...

谢谢！我已经试出来了，就是这个顺序。
通过与gaussian打印的MO因子对比，我发现Orca输出的molden文件h+/-3、+/-4对应的MO因子要乘以-1。但是即便这样也无法通过Multiwfn波函数检测，可能是检测程序对h函数的支持还有问题。例如，用gaussian做下面的计算，产生fch文件，可以通过检测。但是用100/2/6功能转成molden以后，就不能通过检测了。

1. %chk=test
   
2. # hf/gen nosym
   
3. 
   
4. test calculation
   
5. 
   
6. -1 1
   
7. He                -0.56122609    0.59932388    0.04945757
   
8. H                 -0.22790420   -0.34348921    0.04945757
   
9. H                 -0.22788699    1.17072405    0.86595430
   
10. H                 -0.21788699    1.07172405   -0.76703916
    
11. 
    
12. He 0
    
13. h 1 1.0
    
14. 1.0 1.0
    
15. ****
    
16. H 0
    
17. s 1 1.0
    
18. 0.7 1.0
    
19. ****
    

复制代码

sobereva

beefly 发表于 2020-6-4 17:58
谢谢！我已经试出来了，就是这个顺序。
通过与gaussian打印的MO因子对比，我发现Orca输出的molden文件h+ ...

谢谢提示关于h的系数符号需要反转的问题。

我今天更新了官网上的Multiwfn 3.7(dev)，已经完美支持了ORCA产生的含h角动量的.molden文件。
原先不能正常支持原因有两点：
(1) 系数没有修正
(2) ORCA导出含h的molden的时候没有在此文件里写上[11H]，导致Multiwfn将之误当成笛卡尔型的考虑。
现在这两个问题都解决了。