解决了一个Multiwfn产生mkl文件作为ORCA当初猜波函数的问题

sobereva

之前写过一篇文章
将Gaussian等程序收敛的波函数作为ORCA的初猜波函数的方法
http://sobereva.com/517（http://bbs.keinsci.com/thread-14905-1-1.html）

当时做测试的时候只考虑了带有d基函数的情况。但ORCA的基函数归一化条件比较疯狂，F(+3,-3) 、G(+3,-3,+4,-4) 、H(+3,-3,+4,-4)都是归一化为-1而不是1。虽然Multiwfn在载入ORCA的molden文件时考虑了，但当时导出mkl文件时没考虑这点，导致含有>=f基函数的时候，Gaussian里收敛的波函数按照上文给ORCA当初猜的时候还是需要做很多轮迭代，甚至根本不收敛。

今天在官网http://sobereva.com/multiwfn更新的Multiwfn已经解决了此问题。下面的附件里给了一个测试例子，是HF/cc-pVQZ下计算CH3NH2三重态（此时最高涉及g基函数），先用Gaussian做的计算，然后fch文件用Multiwfn转化为了gau.mkl并进一步用orca_2mkl转化为了gau.gbw，在ORCA计算时作为初猜波函数使用，可见一轮就收敛了，非常顺利，而且结果和Gaussian最后给出的精确一致。请使用Multiwfn这个功能的用户尽快更新程序，若遇到问题欢迎反馈。

--------------
SCF ITERATIONS
--------------
ITER       Energy         Delta-E        Max-DP      RMS-DP      [F,P]     Damp
               ***  Starting incremental Fock matrix formation  ***
  0    -95.0605877667   0.000000000000 0.00000177  0.00000001  0.0000032 0.7000
                 **** Energy Check signals convergence ****
               *****************************************************
               *                     SUCCESS                       *
               *           SCF CONVERGED AFTER   1 CYCLES          *
               *****************************************************


此例涉及的所有文件见此
Desktop.rar (1.06 MB, 下载次数 Times of downloads: 42)

2020-7-1 06:38 上传 Uploaded

点击下载
Click to download

neocc

请问社长自定义基组和ecp得到的gbw转成的mkl UO2-orca.mkl (304.98 KB, 下载次数 Times of downloads: 2)

2022-3-3 10:01 上传 Uploaded

点击下载
Click to download

，multiwfn分析ESP会闪退请问是什么原因呢
具体基组信息请见：ORCA 5.0.2如何手动添加基组和ECP信息 http://bbs.keinsci.com/thread-27897-1-1.html

sobereva

neocc 发表于 2022-3-3 10:02
请问社长自定义基组和ecp得到的gbw转成的mkl，multiwfn分析ESP会闪退请问是什么原因呢
具体基组信息请见： ...

Multiwfn根本不支持用mkl格式做波函数分析，支持什么下文和手册2.5节都明确说了
详谈Multiwfn支持的输入文件类型、产生方法以及相互转换
http://sobereva.com/379（http://bbs.keinsci.com/thread-6020-1-1.html）

neocc

sobereva 发表于 2022-3-3 21:26
Multiwfn根本不支持用mkl格式做波函数分析，支持什么下文和手册2.5节都明确说了
详谈Multiwfn支持的输入 ...

抱歉社长，我用错了orca_2mkl命令，误把mkl文件当成了xxx.molden.input
molden的读取文件非常顺利