请问如何用multiwfn画orca两个轨道间的密度差图

michaelm

有ORCA的gbw文件，需要制作其中两个轨道产生的密度差图。
几点尝试，
1）因为使用了ecp，用orca_2aim没有得到wfn或者wfx文件。
2）不知道为什么，用orca_2mkl得到的molden文件不包含virtual orbital。但是我需要绘制的恰恰是一个occupied和一个virtual轨道之间的密度差。
3）用orca_plot产生了需要的两个轨道的cube文件。

所以我想用这两个轨道的cube文件，用multiwfn来产生相应密度之间的密度差文件，存为新的cube文件。能否请大家指点一下，在multiwfn里怎么做？

谢谢！

snljty

Multiwfn的主功能13有格点运算功能，先加载一个，按照提示加载另一个就行。
不过“不知道为什么，用orca_2mkl得到的molden文件不包含virtual orbital”这句话不是很理解，我是没遇到过这种问题，能否贴一下输入输出文件，至少是关键词？
照理说用molden文件，改一下里面的原子序数即核电荷数为实际考虑ECP以后的有效核电荷数，然后用Multiwfn读取molden的波函数信息也能做密度差分析。

sobereva

看Multiwfn手册4.13.2节关于Multiwfn的格点数据运算功能的例子
“用orca_2mkl得到的molden文件不包含virtual orbital”我也从没遇到过。应争取把这个解决，之后结合Multiwfn做各种分析省事得多。

michaelm

谢谢两位。这两天多琢磨了一阵multiwfn在cube文件上的用法，搞明白了不少东东。谢谢指点了。