求助计算福井函数f0和绘制出来的等值面图不一致的问题

ppy412

各位老师好，我使用Gaussian程序按照博文（http://sobereva.com/484）对一个自由基进行了福井函数的计算，然后绘制了福井函数的等值面图，但是发现福井函数等值面图和数据对不上：两个O原子的f0都是0.15，但是在等值面上其中有一个氧的等值面没有显示，请老师指点。

sobereva

福井函数明明是Multiwfn算的而不是Gaussian算的，描述时不要有误导性

计算方式不当，等值面都被截断了。严格照着博文操作，从头把博文里的例子完整重复一遍，确保用的Multiwfn最新版本、没有手动做多余的事

ppy412

谢谢老师指点，我刚刚重复例子做了一次，首先输入22/1/回车/输入0 2/输入-1 1/输入1 1得到了三个gjf文件，然后调用Gaussian进行计算，得到了三个wfn文件，然后输入2，得到简缩福井函数的数据（没有做关于ωcubic和亲电描述符（ε）的计算）然后输入3，计算福井函数等值面，输入3选择高质量格点，输入3得到f0的等值面图，但依然是有截断的，还是没有解决这个问题。

sobereva

N 1.wfn改名为N+1.wfn，和N.wfn、N-1.wfn放到Multiwfn目录下，启动Multiwfn载入N.wfn，然后输入
22
3
3
之后选3观看f0，用0.01等值面，如下所示没有任何问题


选7得到f0.cub，用VMD按下文显示，也没任何问题
在VMD里将cube文件瞬间绘制成效果极佳的等值面图的方法
http://sobereva.com/483（http://bbs.keinsci.com/thread-13329-1-1.html）

ppy412

谢谢卢老师，+符号问题我这边是正常的，可能是上传的时候出现了问题。

问题已经解决了，按照卢老师的步骤成功了，归结原因应该是我没有拖拽N.wfn进Multiwfn导致，而是直接让其读取。