计算化学公社

标题: Multiwfn保存的molden支持h函数，请问是如何排序的？ [打印本页]

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

作者
Author: zjxitcc 时间: 2020-6-1 18:04
我也对这个问题挺好奇的，是跟哪个软件保持一致么还是？molden的H问题偶尔会是个令人头疼的问题

作者
Author: niobium 时间: 2020-6-1 18:51
两位z老师可以看一下源代码啊

作者
Author: sobereva 时间: 2020-6-2 18:11
我今天刚更新了官网上的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格式。

作者
Author: beefly 时间: 2020-6-4 17:58
本帖最后由 beefly 于 2020-6-4 17:59 编辑

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以后，就不能通过检测了。

%chk=test
# hf/gen nosym
test calculation
-1 1
He -0.56122609 0.59932388 0.04945757
H -0.22790420 -0.34348921 0.04945757
H -0.22788699 1.17072405 0.86595430
H -0.21788699 1.07172405 -0.76703916
He 0
h 1 1.0
1.0 1.0
****
H 0
s 1 1.0
0.7 1.0
****

复制代码

作者
Author: sobereva 时间: 2020-6-6 11:50

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将之误当成笛卡尔型的考虑。
现在这两个问题都解决了。

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