计算化学公社

标题: 求助计算福井函数f0和绘制出来的等值面图不一致的问题 [打印本页]

作者
Author: ppy412 时间: 2024-11-28 16:51
标题: 求助计算福井函数f0和绘制出来的等值面图不一致的问题
各位老师好，我使用Gaussian程序按照博文（http://sobereva.com/484）对一个自由基进行了福井函数的计算，然后绘制了福井函数的等值面图，但是发现福井函数等值面图和数据对不上：两个O原子的f0都是0.15，但是在等值面上其中有一个氧的等值面没有显示，请老师指点。

作者
Author: sobereva 时间: 2024-11-28 17:35
福井函数明明是Multiwfn算的而不是Gaussian算的，描述时不要有误导性

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

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

作者
Author: sobereva 时间: 2024-11-28 18:59
N 1.wfn改名为N+1.wfn，和N.wfn、N-1.wfn放到Multiwfn目录下，启动Multiwfn载入N.wfn，然后输入
22
3
3
之后选3观看f0，用0.01等值面，如下所示没有任何问题
(, 下载次数 Times of downloads: 16)

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

作者
Author: ppy412 时间: 2024-11-28 19:28
谢谢卢老师，+符号问题我这边是正常的，可能是上传的时候出现了问题。

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

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