24小时热门版块排行榜    

查看: 4788  |  回复: 8

小樱123

新虫 (初入文坛)

[求助] 如何在MS6.0的DMol3模块下,分析差分电荷密度,具体的操作怎样,新手求助,不胜感激

还有,差分电荷密度和形变电荷密度有什么区别?
另外,有木有大神会DMol3模块的布局分析啊?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

★ ★ ★
gmy1990: 金币+3, 感谢应助 2015-01-23 23:07:37
差分电荷指的是片段结合为体系前后电荷密度差,形变电荷指的是原子结合成体系前后的电荷密度差。

DMol3不能直接进行差分电荷的计算,需要手工处理。

如果一定要做,你可以将片段A和片段B结合的体系AB进行结构优化并计算电荷密度,然后对优化后的结构分别删去A(剩下B)进行保存以及删去B(剩下A)进行保存,之后分别按照完全一样精度以及一样的密度生成格点计算剩余的AB单独的片段的单点能及电荷密度,然后分别将grd文件提取出来,对三个grd文件对应的数据进行做差。所以这个过程建议采用PBC模式使得计算密度格点的边界以及网格一致,否则不能够做差。一般建议格点的生成精度不要小于fine。

至于三个格点文件都能够通过记事本打开,也可以改名字为txt通过origin打开,或者自己写程序使三个格点文件一一对应进行做差。

至于布居分析,DMol对于周期体系只能进行Mulliken和Hirshfeld布居,孤立体系还可以做overlap相关以及ESP的分析。这些在计算的对话框能够非常简单地实现。
不一定挂在论坛,计算问题问题欢迎留言。
2楼2015-01-23 18:51:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ustc

金虫 (著名写手)

好东西!

[ 发自小木虫客户端 ]
3楼2015-01-24 08:14:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小樱123

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by 卡开发发 at 2015-01-23 18:51:26
差分电荷指的是片段结合为体系前后电荷密度差,形变电荷指的是原子结合成体系前后的电荷密度差。

DMol3不能直接进行差分电荷的计算,需要手工处理。

如果一定要做,你可以将片段A和片段B结合的体系AB进行结构 ...

形变电荷密度和差分电荷密度有啥区别呢?
4楼2015-01-30 16:57:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

★ ★
gmy1990: 金币+2 2015-01-31 01:55:48
引用回帖:
4楼: Originally posted by 小樱123 at 2015-01-30 16:57:29
形变电荷密度和差分电荷密度有啥区别呢?...

形变电荷指的是形成整个体系原子间的电荷得失,差分电荷指的是两个参考片段之间形成整个体系这些片段电荷的得失。

举个例子,比如CH3OH,如果形变电荷指的就是形成CH3OH之后各个原子电荷密度的变化;而差分电荷指的是:以CH3与OH形成CH3OH,那么形成前后电荷密度的变化,当然,参考的片段不同结果也是不同的,具体要看想考察那两个部分的结合,如H和CH3O得到的差分电荷与前面就不同。对于双原子分子,差分电荷等于形变电荷。
不一定挂在论坛,计算问题问题欢迎留言。
5楼2015-01-30 21:01:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goudiaocao

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by 卡开发发 at 2015-01-23 18:51:26
差分电荷指的是片段结合为体系前后电荷密度差,形变电荷指的是原子结合成体系前后的电荷密度差。

DMol3不能直接进行差分电荷的计算,需要手工处理。

如果一定要做,你可以将片段A和片段B结合的体系AB进行结构 ...

前辈你好,应该如何对三个格点不一致的grd文件做差?您说的采用PBC模式是什么啊?
6楼2017-01-16 10:17:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

引用回帖:
6楼: Originally posted by goudiaocao at 2017-01-16 10:17:18
前辈你好,应该如何对三个格点不一致的grd文件做差?您说的采用PBC模式是什么啊?...

格点不一致的话只能通过插值等方法(这就降低了可靠性)处理到格点一致,否则没有程序能够处理。PBC指的是周期性边界条件,只有这样你才能够保证同一原子能够在分数坐标重合以及格点边界尺寸一致。

实际上你要清楚电荷密度差在数值处理上的具体过程,实际就是密度的格点在实空间一一对应进行做差,要满足这个条件的话,你可以考虑一下应该怎么处理才不会引进其他的数值误差,看看是否和我前面提到的一致。
不一定挂在论坛,计算问题问题欢迎留言。
7楼2017-01-16 12:33:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goudiaocao

新虫 (初入文坛)

引用回帖:
7楼: Originally posted by 卡开发发 at 2017-01-16 12:33:04
格点不一致的话只能通过插值等方法(这就降低了可靠性)处理到格点一致,否则没有程序能够处理。PBC指的是周期性边界条件,只有这样你才能够保证同一原子能够在分数坐标重合以及格点边界尺寸一致。

实际上你要清 ...

受教  非常感谢!

发自小木虫IOS客户端
8楼2017-01-16 14:42:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

记忆中的红罂

新虫 (小有名气)

引用回帖:
5楼: Originally posted by 卡开发发 at 2015-01-30 21:01:46
形变电荷指的是形成整个体系原子间的电荷得失,差分电荷指的是两个参考片段之间形成整个体系这些片段电荷的得失。

举个例子,比如CH3OH,如果形变电荷指的就是形成CH3OH之后各个原子电荷密度的变化;而差分电荷 ...

前辈你好,我想问下在MS中的DMOL3中的deformation density 是不是差分电荷呢?那对于团簇体系来说,是不是只有差分电荷而不存在形变电荷密度呢?
9楼2017-10-19 08:26:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 小樱123 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 ASdOkHsho7FD 2026-08-28 5/250 2026-08-29 11:49 by jCd0dEvKHShX
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ASdOkHsho7FD 2026-08-28 6/300 2026-08-29 11:40 by jCd0dEvKHShX
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 ASdOkHsho7FD 2026-08-28 6/300 2026-08-29 11:29 by jCd0dEvKHShX
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 ASdOkHsho7FD 2026-08-28 9/450 2026-08-29 11:18 by jCd0dEvKHShX
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +14 晴天加油 2026-08-26 15/750 2026-08-29 11:06 by biily
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 9/450 2026-08-29 09:59 by jklily
[基金申请] 有没有仍没收到信息的 +3 德尚中行 2026-08-27 3/150 2026-08-28 15:24 by rongshuxia
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
信息提示
请填处理意见