应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 求助关于电荷密度计算的问题
81/1返回列表
查看: 2149  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

xixi1007

木虫 (小有名气)

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

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

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

» 猜你喜欢

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

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2015-03-26 15:14:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)

【答案】应助回帖

感谢参与，应助指数 +1
改变了键长就改变了原子坐标，二者不能作电子密度差。
赞一下 回复此楼
2楼2015-03-26 15:27:43
已阅   回复此楼   关注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的回帖

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的回帖

zhou2009

版主 (著名写手)

【答案】应助回帖

实在要作，而又不违反原则，可以这样做：
1、优化出CH3CONHCH3，作>N-H与>N和H的密度差。>N和H的坐标取>N-H中的。可以看到二者成键前后电子的变化。
2、改变N-H，增加键长dr，作此时>N-H与>N和H的密度差。>N和H的坐标取>N-H中的。可以看到二者键长增长前后电子的变化。
3、改变N-H，减少键长dr，作此时>N-H与>N和H的密度差。>N和H的坐标取>N-H中的。可以看到二者键长减少前后电子的变化。
比较这三者的图形，看NH键长增加dr和减小dr后电子密度差的变化。
赞一下 回复此楼
6楼2015-03-27 17:15:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xixi1007

木虫 (小有名气)
引用回帖:
6楼: Originally posted by zhou2009 at 2015-03-27 17:15:47
实在要作，而又不违反原则，可以这样做：
1、优化出CH3CONHCH3，作>N-H与>N和H的密度差。>N和H的坐标取>N-H中的。可以看到二者成键前后电子的变化。
2、改变N-H，增加键长dr，作此时>N-H与>N和 ...

我刚做这个电荷密度，不是很懂，对multiwfn程序也不太了解。怎样做N-H，与N和H的密度差，是从优化后的结果中，单独取N-H，N，H的坐标，然后分别做这三部分的密度，之后相减吗。
赞一下 回复此楼
7楼2015-03-30 09:58:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
xixi1007: 金币+10, 非常感谢您的帮助！ 2015-03-30 15:29:15
引用回帖:
7楼: Originally posted by xixi1007 at 2015-03-30 09:58:38
我刚做这个电荷密度，不是很懂，对multiwfn程序也不太了解。怎样做N-H，与N和H的密度差，是从优化后的结果中，单独取N-H，N，H的坐标，然后分别做这三部分的密度，之后相减吗。...

是的。

>N-H与>N和H的密度相减。>N是包含其上所连基团的。
这三者都取>N-H中的坐标，进行计算，不可优化，加关键词nosymm，不让它作对称性分析，否则坐标又改变了。

最好是在计算时加关键词：0utput=wfn，要求输出wfn文件，这样在multiwfn中作密度差更方便。
详见：《使用Multiwfn作电子密度差图》
http://muchong.com/bbs/viewthread.php?tid=3980314
赞一下 回复此楼
8楼2015-03-30 11:10:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xixi1007 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 中青基了要发朋友圈吗？ +3 349506619 2026-08-28 3/150 2026-08-29 10:49 by 木风9012
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 投票:  有多少人是今天查系统知道结果的？ +16 爱看书的可乐 2026-08-26 18/900 2026-08-29 10:18 by winsaint
[考博] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 ASdOkHsho7FD 2026-08-28 5/250 2026-08-29 09:49 by G6APbkg8SA6w
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 ASdOkHsho7FD 2026-08-28 4/200 2026-08-29 09:48 by G6APbkg8SA6w
[找工作] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 ASdOkHsho7FD 2026-08-28 5/250 2026-08-29 09:28 by G6APbkg8SA6w
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 为什么资助数各大高校都创新高，自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 基金不中，共勉 +11 eulota 2026-08-26 11/550 2026-08-28 14:22 by 火星超人xi
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜？ 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解！ +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 今天务委会开完了，明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 国合现在查不到了吗？ +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 在坚冰还盖着北海的时候，我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 某些机构，以效率低为荣，以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
信息提示
关闭
请填处理意见
关闭