计算化学公社

标题: 解决了一个Multiwfn产生mkl文件作为ORCA当初猜波函数的问题 [打印本页]

作者
Author: sobereva 时间: 2020-7-1 06:44
标题: 解决了一个Multiwfn产生mkl文件作为ORCA当初猜波函数的问题
之前写过一篇文章
将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       *
            *****************************************************

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

作者
Author: neocc 时间: 2022-3-3 10:02
请问社长自定义基组和ecp得到的gbw转成的mkl (, 下载次数 Times of downloads: 2) ，multiwfn分析ESP会闪退请问是什么原因呢
具体基组信息请见：ORCA 5.0.2如何手动添加基组和ECP信息 http://bbs.keinsci.com/thread-27897-1-1.html

作者
Author: sobereva 时间: 2022-3-3 21:26

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）

作者
Author: neocc 时间: 2022-3-3 23:46

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

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

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