计算化学公社

标题: 使用Multiwfn和GaussView绘制电子密度差图的区别在哪里？ [打印本页]

作者
Author: loovfnd 时间: 2018-8-30 09:20
标题: 使用Multiwfn和GaussView绘制电子密度差图的区别在哪里？
本帖最后由 loovfnd 于 2018-8-30 09:23 编辑

1.使用Multiwfn
首先优化了一个体系的基态结构。分别用（电荷和自旋多重度都是1, 1）%chk=gs.chk # b3lyp/6-311g(2d,p)得到gs.chk%chk=es.chk # td=(root=1,singlet) b3lyp/6-311g(2d,p) density得到es.chk

使用Multiwfn计算电子密度差，参照http://sobereva.com/113中方法计算过程很顺利，但在观察三维图的时候出现了问题，
在Multiwfn中调节iso值始终观察不到，故改用导出cube文件使用GaussView，也有错误提示，是设置问题吗?。如下图：
(, 下载次数 Times of downloads: 43) (, 下载次数 Times of downloads: 61)

2.使用GaussView
GaussView得到电子密度差只使用了es.chk，分别使用CI和SCF得到电子密度差，最后生成图形,过程如下图。
这里有个问题，如果Multiwfn能正常得到电子密度差图，GaussView仅使用es.chk和Multiwfn使用两个es.chk和gs.chk得到的差是一致的吗?如果不一致哪个更准确。
(, 下载次数 Times of downloads: 47) (, 下载次数 Times of downloads: 50)

作者
Author: sobereva 时间: 2018-8-30 15:09
1
Multiwfn中若看不到密度差等值面，而且不管isovalue调到多小也无济于事，说明计算过程错误，用任何其它程序也一样看不到

明显不能直接用fch文件求差，看：
在Multiwfn中基于fch产生自然轨道的方法与激发态波函数、自旋自然轨道分析实例
http://sobereva.com/403
否则相当于都是基态密度在求差，结果显然处处为0
Multiwfn用的是.fch里的轨道而不是密度矩阵
若非要用.fch，也得先按照上文的做法产生激发态自然轨道并写入激发态的.fch文件里

2 一样
Multiwfn绝对是最准确、最方便的，完全没有丝毫必要用gview做这事情

作者
Author: loovfnd 时间: 2018-8-30 15:36

sobereva 发表于 2018-8-30 15:09
1
Multiwfn中若看不到密度差等值面，而且不管isovalue调到多小也无济于事，说明计算过程错误，用任何其它 ...

明白了，最便利的还是使用wfn文件。

作者
Author: sobereva 时间: 2022-10-13 11:18

忘不掉的回忆发表于 2022-10-13 11:05
老师好，我也遇到了本贴的第一个问题。用的Multiwfn是2022-Oct-10更新的。
在Multiwfn中若看不到密度差 ...

0.0005 a.u.下的等值面，毫无问题

(, 下载次数 Times of downloads: 18)

作者
Author: 忘不掉的回忆 时间: 2022-10-13 11:54

sobereva 发表于 2022-10-13 11:18
0.0005 a.u.下的等值面，毫无问题

谢谢老师，是我调节的问题。
再次感谢社长。

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