计算化学公社

标题: 为什么用ORCA输出的molden文件计算CDFT中Multiwfn读不对能量 [打印本页]

作者
Author: david25 时间: 2022-9-2 14:50
标题: 为什么用ORCA输出的molden文件计算CDFT中Multiwfn读不对能量
本帖最后由 david25 于 2022-9-2 15:00 编辑

利用ORCA计算钴原子和两个甲基咪，控制参数如下
"BP RI def2-TZVP D3  KEEPDENS KDIIS SOSCF TightSCF "

电荷和multiplicity 设为  0 和 4，计算获得结果。
对于N+1的相关计算，设置为 -1 和 5；对于 N-1，设置为 1 和 5。而其余输入不变，得到结果。

输入orca_2mkl N -molden，产生molden文件，其余的产生对应wfn文件。
再利用Multiwfn的“100-2-5”产生对应的 wfn。
在做CDFT的计算时，我在终端上输入"Multiwfn N.molden.input"，进入Multiwfn后，输入“22-2”，自动获取目录下N.wfn, N-1. wfn, N+1.wfn,  获得了CDFT.txt文件。
但CDFT.txt中Electrophilicity为Nan值，还有其他的量要么是Nan要么是infinity, 详情见附件。注意到CDFT.txt对E(N) E(N+1) E(N-1) 以及对应的HOMO的值的记录如下图：
(, 下载次数 Times of downloads: 7)
检查三个计算的 final one-point energy 和 HOMO的能量，都是正常的值。在N.molden.input中才搜索到5.123575这个值。
为了解决这个问题，我做了如下尝试：
1、在作CDFT计算时，分别用了“Multiwfn N.wfn”和 “Multiwfn N.xyz”做启动，E_HOMO(N+1)会变对，其他的能量都为0，该为Nan值的继续为Nan值。
2、直接利用molden文件而不是wfn，此时E_HOMO都对，但是E(N),E(N+1), E(N-1)都为0，其余的该为Nan值的为Nan值。

对于ORCA，Multiwfn和 CDFT，我是新手中的新手。对于上面的问题，还请各位老师帮忙指导一下

作者
Author: sobereva 时间: 2022-9-2 15:07
怎么可能N+1和N-1都是5重态
molden文件里不记录总能量，molden转成wfn再用这功能没意义
直接让ORCA产生.wfn文件，或者让Multiwfn直接调用本机里的ORCA产生.wfn文件就完了。

此文说了
使用Multiwfn超级方便地计算出概念密度泛函理论中定义的各种量
http://sobereva.com/484（http://bbs.keinsci.com/thread-13328-1-1.html）

另外，如果你没有Gaussian而只有ORCA，也可以选择-2 Choose the quantum chemistry program used in option 1之后再选2。之后再选选项1产生输入文件时，Multiwfn产生的就是三个态的ORCA的输入文件N.inp、N+1.inp和N-1.inp，手动算完了之后就会在当前目录下产生N.wfn、N+1.wfn和N-1.wfn（其它产生的文件可删掉）。如果你把settings.ini里的orcapath设为了实际ORCA可执行文件路径，那么Multiwfn在产生输入文件后会问你是否自动调用ORCA去计算它们并在计算结束后删掉其它文件。

作者
Author: david25 时间: 2022-9-12 19:24

sobereva 发表于 2022-9-2 15:07
怎么可能N+1和N-1都是5重态
molden文件里不记录总能量，molden转成wfn再用这功能没意义
直接让ORCA产生.w ...

谢谢社长回复。根据社长回复我发现是我的多重态设置弄错了。

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