24小时热门版块排行榜    

查看: 1447  |  回复: 9

流浪者之辉

铁虫 (初入文坛)


[交流] 求问siesta中mulliken的问题

本人对单独的一个水分子进行优化,得到的mulliken分布如下:
mulliken: Atomic and Orbital Populations:

Species: H                  
Atom  Qatom  Qorb
               1s      1s      1Ppy    1Ppz    1Ppx   
   2  1.155   0.446   0.455   0.108   0.106   0.040
   3  1.155   0.446   0.455   0.108   0.104   0.042

Species: O                  
Atom  Qatom  Qorb
               2s      2s      2py     2pz     2px     2py     2pz     2px     
               2Pdxy   2Pdyz   2Pdz2   2Pdxz   2Pdx2-y2
   1  5.690   0.781   0.575   1.749   1.339   0.602   0.033   0.135   0.456
              0.000   0.003   0.003   0.012   0.002

mulliken: Qtot =        8.000
我的理解是以上的理解是mulliken的总电荷是价电子,根据上面氢原子得到电子而氧原子失去电子,跟实际不符合啊,但又找不到问题所在,求教大家帮忙。
下面文件是输入文件。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

流浪者之辉(金币+1): 谢谢参与
比较可能的原因我想可能是Mulliken Pop对基组依赖性比较强,对于赝轨道分析可能不太合适,结果有可能意义比较弱。Hirshfeld Pop一般不会出现这样的情况。
6楼2019-05-25 06:11:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

流浪者之辉

铁虫 (初入文坛)


引用回帖:
6楼: Originally posted by 卡开发发 at 2019-05-25 06:11:49
比较可能的原因我想可能是Mulliken Pop对基组依赖性比较强,对于赝轨道分析可能不太合适,结果有可能意义比较弱。Hirshfeld Pop一般不会出现这样的情况。

谢谢回答,再请教一下:如果我要反映必要有物理意义的结果,选择哪种电荷分析比较好呢,Hirshfeld还是bader。(按照上述情况siesta中mulliken是不是就完全不能合理实现)呢)
7楼2019-05-27 17:04:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
7楼: Originally posted by 流浪者之辉 at 2019-05-27 17:04:40
谢谢回答,再请教一下:如果我要反映必要有物理意义的结果,选择哪种电荷分析比较好呢,Hirshfeld还是bader。(按照上述情况siesta中mulliken是不是就完全不能合理实现)呢)...

我个人觉得Hirshfeld好点,但通常Hirshfeld的值可能不能还原静电势,Bader的那个zero-flux虽然数学意义有,但不一定符合化学经验(你可以参考http://sobereva.com/161)。基于Mulliken的划分有时候确实比较坑人,比如PDOS可能做出来是负的。
8楼2019-05-27 19:41:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

流浪者之辉

铁虫 (初入文坛)


引用回帖:
8楼: Originally posted by 卡开发发 at 2019-05-27 19:41:28
我个人觉得Hirshfeld好点,但通常Hirshfeld的值可能不能还原静电势,Bader的那个zero-flux虽然数学意义有,但不一定符合化学经验(你可以参考http://sobereva.com/161)。基于Mulliken的划分有时候确实比较坑人,比 ...

谢谢您的答复!!!
9楼2019-05-27 21:03:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
longwave2楼
2019-05-22 19:27   回复  
流浪者之辉(金币+1): 谢谢参与
祝福 发自小木虫Android客户端
tzynew3楼
2019-05-23 05:43   回复  
流浪者之辉(金币+1): 谢谢参与
2019-05-23 12:58   回复  
流浪者之辉(金币+1): 谢谢参与
nono20095楼
2019-05-23 22:18   回复  
流浪者之辉(金币+1): 谢谢参与
`
2019-05-30 02:11   回复  
相关版块跳转 我要订阅楼主 流浪者之辉 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +15 爱看书的可乐 2026-08-26 17/850 2026-08-27 10:35 by xxniao123
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +4 a089 2026-08-26 5/250 2026-08-27 09:24 by Howard28
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 17/850 2026-08-27 08:20 by 大萍1987
[基金申请] 国自然面上复盘~欢迎讨论 +7 晴天加油 2026-08-26 8/400 2026-08-27 07:53 by longfie172
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 04:33 by Ie9AyIAvGbvs
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 02:01 by Ie9AyIAvGbvs
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 怎么查啊 +5 huang1991js 2026-08-26 5/250 2026-08-26 16:09 by Equinoxhua
[基金申请] 范进中举一文的中心思想 +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 宝贝虫子
[基金申请] 科研孤儿太难了 +21 我4大白菜 2026-08-20 22/1100 2026-08-26 11:29 by nlgza
[基金申请] 国际合作可查了,中了面上 +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
信息提示
请填处理意见