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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » MS » 【求助】Dmol3中hirshfeld电荷布局分布的问题
北京石油化工学院2026年研究生招生接收调剂公告
61/1返回列表
查看: 1513  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

topten-gg

银虫 (正式写手)

[交流] 【求助】Dmol3中hirshfeld电荷布局分布的问题

我计算了一个中性体系 想用hirshfeld方法看一看电荷的布局分布 并与常用的米利根比较一下
发现在结果文件中出现的总电荷数不是0
如下
Charge partitioning by Hirshfeld method: (   0.00056)
  N     1 charge                          -0.0529
  C     2 charge                           0.0496
  H     3 charge                           0.0670
  N     4 charge                          -0.2077
  C     5 charge                          -0.0119
  C     6 charge                           0.1048
  N     7 charge                          -0.1385
  H     8 charge                           0.1363
………………（这是input中设置Hirshfeld_analysis            charge的情况）

Charge partitioning by Hirshfeld method: (   0.00056)
  N     1 charge                          -0.0529
          dipoles x,y,z  [au]             -0.0274    0.0343    0.0002
  C     2 charge                           0.0496
          dipoles x,y,z  [au]             -0.0534    0.0507    0.0000
  H     3 charge                           0.0670
          dipoles x,y,z  [au]             -0.1619   -0.0671    0.0000
  N     4 charge                          -0.2077
          dipoles x,y,z  [au]              0.0135    0.2289    0.0000
  C     5 charge                          -0.0119
          dipoles x,y,z  [au]              0.0157   -0.0536    0.0000
  C     6 charge                           0.1048
          dipoles x,y,z  [au]              0.0191   -0.0381    0.0000
  N     7 charge                          -0.1385
          dipoles x,y,z  [au]              0.0019   -0.0174    0.0000
  H     8 charge                           0.1363
……………………（这是input中设置Hirshfeld_analysis            dipole的情况）

或者
Charge partitioning by Hirshfeld method: (   0.00056)
  N     1 charge                          -0.0529
          dipoles x,y,z  [au]             -0.0274    0.0343    0.0002
          quadrupoles xx yy zz yz xz xy    0.0405    0.0023   -0.0429    0.0003   -0.0003   -0.0358
  C     2 charge                           0.0496
          dipoles x,y,z  [au]             -0.0534    0.0507    0.0000
          quadrupoles xx yy zz yz xz xy   -0.1215   -0.1341    0.2556    0.0003   -0.0001    0.1300
  H     3 charge                           0.0670
          dipoles x,y,z  [au]             -0.1619   -0.0671    0.0000
          quadrupoles xx yy zz yz xz xy   -0.0390    0.0464   -0.0074    0.0000    0.0000   -0.0387
  N     4 charge                          -0.2077
          dipoles x,y,z  [au]              0.0135    0.2289    0.0000
          quadrupoles xx yy zz yz xz xy    0.2547   -0.3819    0.1272   -0.0001    0.0000    0.0816
  C     5 charge                          -0.0119
          dipoles x,y,z  [au]              0.0157   -0.0536    0.0000
          quadrupoles xx yy zz yz xz xy   -0.1625   -0.0630    0.2255    0.0000    0.0001   -0.1269
…………………………（这是input中设置Hirshfeld_analysis            quadrupole的情况）


这里出现的非0的小数是怎么回事
请教大侠们怎么解决


bow！！
回复此楼

» 猜你喜欢
1楼 2009-10-30 10:37:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qscao

木虫 (小有名气)
★ ★
小木虫(金币+0.5):恭喜抢沙发，给个红包
aylayl08(金币+1,VIP+0):谢谢参与 10-30 17:36
小数表示该原子带的电荷
赞一下(7人) 回复此楼
2楼2009-10-30 15:23:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

topten-gg

银虫 (正式写手)
引用回帖:
Originally posted by qscao at 2009-10-30 15:23:
小数表示该原子带的电荷

但是总的原子所带的电荷加起来 并不是0啊
赞一下 回复此楼
3楼2009-10-30 15:38:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

topten-gg

银虫 (正式写手)
高人 现身吧
回复此楼
4楼2009-10-31 09:49:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yzcluster

金虫 (著名写手)

小木虫扫盲人
★ ★ ★ ★ ★
topten-gg(金币+3,VIP+0):很感谢 10-31 12:38
qasd(金币+2,VIP+0):xiexie 11-2 19:35
你给的HIRSHFELD应该不是input中给出的电荷分布，而是在SCF converged后电子性质计算的一部分。
我记得在SCF计算中，有一项所谓的charge density convergence，如果你的计算结果低于这个值，程序就认为已经收敛了，而此时不是真正数学上的0，所以总电荷可能是一个不为0的很小的数值，可以忽略不计。
这个问题不仅仅是HIRSHFELD方法，MULLIKEN，BADER，NBO等各种布局分析都会出现类似的情况（而且相同的布局分析在不同的程序中结果也不完全一样）。不管怎么说，一般的数值计算很小的误差往往可以忽略。
赞一下(11人) 回复此楼
5楼2009-10-31 10:40:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

topten-gg

银虫 (正式写手)
引用回帖:
Originally posted by yzcluster at 2009-10-31 10:40:
你给的HIRSHFELD应该不是input中给出的电荷分布，而是在SCF converged后电子性质计算的一部分。
我记得在SCF计算中，有一项所谓的charge density convergence，如果你的计算结果低于这个值，程序就认为已经收敛了 ...

那么这个不为0的数值怎么在电荷分布上消除呢？

谢谢你了
赞一下 回复此楼
6楼2009-10-31 12:38:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 topten-gg 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 334分机械专硕求调剂 +3 蛋花紫菜汤 2026-04-03 3/150 2026-04-07 14:49 by 逍遥cocoa
[考研] 化工学硕 285求调剂 +15 Wisjxn 2026-04-07 15/750 2026-04-07 14:45 by 上岸上岸，，
[考研] 285求调剂 +9 AZMK 2026-04-05 14/700 2026-04-07 13:01 by Sammy2
[考研] 302分求调剂 一志愿安徽大学085601 +12 zyx上岸！ 2026-04-04 12/600 2026-04-07 02:09 by BruceLiu320
[考研] 305求调剂 +3 77Qi 2026-04-06 3/150 2026-04-06 15:36 by going home
[考研] 材料工程310专硕调剂 +14 捞捞我…. 2026-04-04 15/750 2026-04-06 14:18 by lqwchd
[考研] 081700化学工程与技术 一志愿中海洋 323 求调剂学校 +18 披星河 2026-04-03 18/900 2026-04-06 13:55 by BruceLiu320
[考研] 372分，材料与化工，一志愿湖南大学，求调剂 +3 蓝笺片 2026-04-01 3/150 2026-04-06 09:04 by 无际的草原
[考研] 322求调剂 +3 嗯哼哼恒 2026-04-05 3/150 2026-04-05 19:52 by nepu_uu
[考研] 0854求调剂 +4 assdll 2026-04-04 4/200 2026-04-05 09:44 by zhq0425
[考研] 359求调剂 +7 hhhhaaaa$ 2026-04-04 7/350 2026-04-04 18:49 by imissbao
[考研] 268求调剂 +8 你好tg 2026-04-03 9/450 2026-04-04 05:08 by gswylq
[考研] 282求调剂 +20 ycy1201 2026-04-01 22/1100 2026-04-04 00:42 by userper
[考研] 350一志愿北京航空航天大学08500材料科学与工程求调剂 +5 kjnasfss 2026-04-03 5/250 2026-04-03 22:29 by 无际的草原
[考研] 重庆大学材料与化工085600，初试370+，求求调剂建议 +8 shzhou_ 2026-04-01 9/450 2026-04-03 09:31 by 蓝云思雨
[考研] 土木304求调剂 +4 兔突突突， 2026-04-02 5/250 2026-04-02 21:16 by 兔突突突，
[考研] 318求调剂 +3 笃行致远. 2026-03-31 4/200 2026-04-02 15:56 by Jaylen.
[考研] 324求调剂 +5 想上学求调 2026-04-01 6/300 2026-04-02 10:16 by sanrepian
[考研] 377求调剂 +3 RASKIN 2026-04-02 3/150 2026-04-02 09:45 by zzchen2000
[考研] 求调剂，一志愿北林食品与营养095500，301分，已过六级，有科研经历 +4 快乐储蓄罐 2026-03-31 4/200 2026-04-01 09:26 by JourneyLucky
信息提示
关闭
请填处理意见
关闭