24小时热门版块排行榜    

查看: 2171  |  回复: 7
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

xixi1007

木虫 (小有名气)

[求助] 求助关于电荷密度计算的问题 已有1人参与

以NMA分子CH3CONHCH3为例,对这个分子进行几何优化和频率计算后,得到其平衡态结构坐标。将NH键长增加dr和减小dr后,得到两个新的结构坐标,分别得到两个结构下的电子密度ρ(dr)和ρ(-dr),我想得到ρ(dr)-ρ(-dr)/2dr的图,具体应该怎么做呢。这种情况下两个结构的电子密度是用cubegen产生呢,还是用cube。得到两个cube文件后,将两个cube文件中的数据相减再除以2dr吗,如果是这样,之后应该怎样作图呢,得到的应该是类似以下图吧。感谢各位不吝赐教!

求助关于电荷密度计算的问题
2.JPG
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xixi1007

木虫 (小有名气)

引用回帖:
2楼: Originally posted by zhou2009 at 2015-03-26 15:27:43
改变了键长就改变了原子坐标,二者不能作电子密度差。

那如果我想做ρ(dr)-ρ(-dr)/2dr的图,应该怎么办呢
3楼2015-03-26 15:52:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

zhou2009

版主 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
改变了键长就改变了原子坐标,二者不能作电子密度差。
2楼2015-03-26 15:27:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
xixi1007: 金币+10, ★★★很有帮助, 非常感谢! 2015-03-27 08:47:59
引用回帖:
3楼: Originally posted by xixi1007 at 2015-03-26 15:52:15
那如果我想做ρ(dr)-ρ(-dr)/2dr的图,应该怎么办呢...

改变了原子坐标,二者作电子密度差是没有意义的,会得出错误的信息和结果。
实在想做,可以用分别求得的ρ(dr)和ρ(-dr)的两个cube文件,在multiwfn程序中进行cube运算:{ρ(dr)-ρ(-dr)}/2dr。
结果得到的还是一个cube文件,可以在GV中作图。
4楼2015-03-26 21:23:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)

【答案】应助回帖

为什么说改变了原子坐标,二者作电子密度差是没有意义的,这是因为原子坐标改变后,电子本身在密度差上就有很大的改变,再加上原子坐标改变引起了相互作用的电子密度改变,两种密度差叠加在一起,无法区分
5楼2015-03-27 15:44:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见