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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » MS » 关于dmol3计算表面吸附分子电荷转移的问题请教
91/1返回列表
查看: 2964  |  回复: 8
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【悬赏金币】回答本帖问题，作者青涩的小学弟将赠送您 30 个金币

青涩的小学弟

新虫 (初入文坛)

[求助] 关于dmol3计算表面吸附分子电荷转移的问题请教

请问高手，在用dmol3模拟小分子表面吸附的过程时，电荷转移具体是怎么得到的？
我重复了MoS2吸附气体小分子的论文，用论文里的方法计算的结构和吸附能都符合，只有电荷转移与论文中的结果不符。论文也是用的Hirshfeld计算的。
我在计算时，是勾选了population analysis——Hirshfeld analysis——  charge ——run 然后在outmol中直接将Hirshfeld 下面charge 的读数与文献中作对比，每次都符合不上。想请问一下 ，实际是的电荷转移是怎么得到的，是有什么式子吗？？

就比如我只算一个小分子的电荷转移，都与文献的结果差很多，下面是我算用GGA-PBE算的一个H2S在格子里的电荷转移，直接找到outmol中的Hirshfeld分析

Charge partitioning by Hirshfeld method: (  0.0005   -0.0000)
  H     1 charge   spin                    0.0495   -0.0000
  H     2 charge   spin                    0.0495    0.0000
  S     3 charge   spin                   -0.0985    0.0000
请问，我能直接说S原子带-0.0985的负电荷吗？

基础比较薄弱，这个问题困扰好几天了，请高手指教
回复此楼

» 收录本帖的淘贴专辑推荐

MS计算模拟

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)
1楼2019-09-06 11:01:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

青涩的小学弟

新虫 (初入文坛)
自己顶一下
回复此楼
2楼2019-09-06 11:18:18
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

青涩的小学弟

新虫 (初入文坛)
请问，吸附分子的带电量怎么计算啊，是放在格子里吗，还是直接几何优化得到原子带电量
赞一下 回复此楼
3楼2019-09-06 18:48:11
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur
引用回帖:
3楼: Originally posted by 青涩的小学弟 at 2019-09-06 18:48:11
请问，吸附分子的带电量怎么计算啊，是放在格子里吗，还是直接几何优化得到原子带电量

1、分子不用放到格子里面。
2、MS版本不同结果可能会有些许差异，对于Hirshfeld应该不会很敏感。
赞一下(1人) 回复此楼

» 本帖已获得的红花（最新10朵）

青涩的小学弟
不一定挂在论坛，计算问题问题欢迎发帖或留言。
4楼2019-09-07 04:59:14
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

青涩的小学弟

新虫 (初入文坛)
送红花一朵
引用回帖:
4楼: Originally posted by 卡开发发 at 2019-09-07 04:59:14
1、分子不用放到格子里面。
2、MS版本不同结果可能会有些许差异，对于Hirshfeld应该不会很敏感。...

谢谢老师回复。
老师您好，那我计算分子表面吸附电荷转移通过population analysis——Hirshfeld analysis——  charge ——run 过程得到输出文件里读取分子的charge值就是分子转移的电荷量吗？
还是应该与单独分子下的charge量作对比才行？  谢谢老师指教


麻烦您点个应助呀
赞一下(1人) 回复此楼
5楼2019-09-07 11:16:14
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur
引用回帖:
5楼: Originally posted by 青涩的小学弟 at 2019-09-07 11:16:14
谢谢老师回复。
老师您好，那我计算分子表面吸附电荷转移通过population analysis——Hirshfeld analysis——  charge ——run 过程得到输出文件里读取分子的charge值就是分子转移的电荷量吗？
还是应该与单独分 ...

那个是分子带的净电荷量（也就是分子各个原子布居分析得到的电子数目减去原子核电荷数），然后需要对相互作用前后的这个量进行做差得到的才是电荷转移的量。

不过既然是布居分析，那么划分的方法不同结果总是会有些差异，其实不妨看看电荷密度差的情况结合来说明。
赞一下(1人) 回复此楼

» 本帖已获得的红花（最新10朵）

青涩的小学弟
不一定挂在论坛，计算问题问题欢迎发帖或留言。
6楼2019-09-07 13:41:45
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

青涩的小学弟

新虫 (初入文坛)
送红花一朵
引用回帖:
6楼: Originally posted by 卡开发发 at 2019-09-07 13:41:45
那个是分子带的净电荷量（也就是分子各个原子布居分析得到的电子数目减去原子核电荷数），然后需要对相互作用前后的这个量进行做差得到的才是电荷转移的量。

不过既然是布居分析，那么划分的方法不同结果总是会 ...

明白了。谢谢老师指教， 老师您说的电荷密度差用doml算的话，对于表面吸附体系，Total density 、Deformation density、Spin density 这三者优先以哪个为主呢
赞一下 回复此楼
7楼2019-09-08 10:16:36
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur
引用回帖:
7楼: Originally posted by 青涩的小学弟 at 2019-09-08 10:16:36
明白了。谢谢老师指教， 老师您说的电荷密度差用doml算的话，对于表面吸附体系，Total density 、Deformation density、Spin density 这三者优先以哪个为主呢...

就是电荷密度差，即charge density difference，不是Total density 、Deformation density、Spin density，这个你得想办法把相互作用的片段用相同精度、相同基组、相同截断半径、相同k点、相同周期性边界、相同的格点密度分别进行计算后得到的grd文件进行做差得到。
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

青涩的小学弟
不一定挂在论坛，计算问题问题欢迎发帖或留言。
8楼2019-09-08 11:25:23
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

青涩的小学弟

新虫 (初入文坛)
送红花一朵
引用回帖:
8楼: Originally posted by 卡开发发 at 2019-09-08 11:25:23
就是电荷密度差，即charge density difference，不是Total density 、Deformation density、Spin density，这个你得想办法把相互作用的片段用相同精度、相同基组、相同截断半径、相同k点、相同周期性边界、相同的格 ...

谢谢您老师，明白了，您提的这个做法我在其他帖子的应助下看到您回复过，我去学习一下，非常感谢您
赞一下 回复此楼
9楼2019-09-09 16:30:06
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 青涩的小学弟 的主题更新
91/1返回列表
不应助 确定回帖应助 (注意：应助才可能被奖励，但不允许灌水，必须填写15个字符以上)
普通表情 龙 兔 虎 猫 论文翻译-EnPapers美助理教授团队-SCI英文润色
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] B口人才项目 +9 WOWO159357 2024-05-29 18/900 2024-06-01 13:33 by hahaqiamg
[硕博家园] 实验室太吵闹，无法安静学习，怎么办？ +6 utahh 2024-05-31 11/550 2024-06-01 13:27 by utahh
[硕博家园] 讨厌酒桌文化 +6 hahamyid 2024-05-31 6/300 2024-06-01 09:03 by lizhiweiouc
[考研] 研0二导师分到新来的博士后靠谱吗 +6 sone9 2024-05-31 6/300 2024-06-01 06:35 by sjianglei
[找工作] 找工作如此之难 +7 探123 2024-05-25 7/350 2024-05-31 15:58 by 平南小刀
[教师之家] 在大地上我们只过一生---看完我的阿勒泰上头了好几天，完结那天晚上几乎失眠 +10 瞬息宇宙 2024-05-27 11/550 2024-05-31 15:38 by 烟 火
[高分子] MMA预聚体光固化发雾问题求助 +3 惠亚金总 2024-05-29 10/500 2024-05-31 14:59 by 惠亚金总
[文学芳草园] 对对子啊 +5 天若孤独 2024-05-29 7/350 2024-05-31 09:00 by wjykycg
[论文投稿] 审稿专家比较坚定的让补充实验，但实在没法补充实验，修回还有希望吗？ (EPI+1) 3+3 qweasd12345 2024-05-29 6/300 2024-05-30 08:11 by qweasd12345
[硕博家园] 又想换工作 +15 brightmj 2024-05-27 20/1000 2024-05-29 23:25 by zyqchem
[考研] 研0 +3 徐小七七 2024-05-25 5/250 2024-05-29 13:49 by 1158057902
[基金申请] 信息学部函评结束了吗？ +6 ducan21 2024-05-28 7/350 2024-05-29 12:10 by WORLD0256
[论文投稿] 核心初审被拒，理由是“选题的意义不明确，文章写得不像是科技论文”，怎么改 5+3 工藤雷花樱 2024-05-27 8/400 2024-05-29 10:09 by topedit
[基金申请] 如果您是国自然评审专家 +3 丁香园账户 2024-05-28 3/150 2024-05-29 06:44 by gaohui8888
[基金申请] E05青基有几个评审 +4 KYXY123 2024-05-28 4/200 2024-05-28 19:25 by popt2t
[材料综合] 求助，刚入门MOF合成，实验总是失败 20+3 张月半 2024-05-25 8/400 2024-05-28 09:20 by 我想更懂你29
[基金申请] 面上基金会评专家，有回避机制吗？ +4 huang1991js 2024-05-27 4/200 2024-05-27 19:08 by 星火12
[基金申请] 感觉自然基金限制通过比例就是有点扯，学学B口，化学学部，不限制比例。 +10 wsjing 2024-05-26 14/700 2024-05-27 11:57 by kanmiaolucky
[硕博家园] 周日 +6 1加油哦棒 2024-05-26 9/450 2024-05-27 10:30 by hahamyid
[硕博家园] 2024博士招生 +3 big 混子 2024-05-26 3/150 2024-05-26 20:47 by 宁多缺毋滥
信息提示
关闭
请填处理意见
关闭