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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » MS » 【求助】Dmol3中hirshfeld电荷布局分布的问题
61/1返回列表
查看: 1517  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿华中农业071010，320求调剂 +17 困困困困坤坤 2026-04-14 19/950 2026-04-17 20:08 by 关一盏灯cd
[考研] 294求调剂 +7 淡然654321 2026-04-17 8/400 2026-04-17 16:36 by wutongshun
[考研] 297，工科调剂? +4 河南农业大学-能 2026-04-14 4/200 2026-04-16 22:52 by wulijun2012
[考研] 材料相关专业344求调剂双非工科学校或课题组 +23 hualkop 2026-04-12 25/1250 2026-04-16 22:12 by SUSE_CL
[考研] 一志愿沪9，生物学326求调剂 +9 刘墨墨 2026-04-15 9/450 2026-04-16 17:14 by 崔崔崔cccc
[考研] 291求调剂 +11 关忆北. 2026-04-14 11/550 2026-04-16 15:18 by jiahl2024
[考研] 求调剂 +11 小聂爱学习 2026-04-11 15/750 2026-04-15 21:57 by noqvsozv
[考研] 通信工程求调剂！！！ +6 zlb770521 2026-04-14 6/300 2026-04-15 20:00 by 学员JpLReM
[考研] 各位老师好，求调剂，本科211，一志愿天津大学生物与医药学硕，差两名录取。 +11 路六六jjj 2026-04-13 11/550 2026-04-14 16:01 by zs92450
[考研] 271求调剂 +35 2261744733 2026-04-11 41/2050 2026-04-14 15:36 by zs92450
[教师之家] 转长聘了 +7 简单化xn 2026-04-13 7/350 2026-04-14 14:50 by xindong
[考研] 105500药学求调剂 +4 x_skys 2026-04-12 4/200 2026-04-14 13:37 by rndfc
[考研] B区0809 ，数一英一，290 求调剂 +3 泠潍1111 2026-04-12 4/200 2026-04-13 20:35 by 学员JpLReM
[考研] 302求调剂 +10 易！? 2026-04-13 10/500 2026-04-13 19:04 by lbsjt
[考研] 339求调剂 +4 hanwudada 2026-04-12 4/200 2026-04-13 12:03 by 蓝云思雨
[考研] 344 材料专业 求调剂211 无地域要求 +8 hualkop 2026-04-11 8/400 2026-04-12 22:24 by fqwang
[考研] 346分，工科0854求调剂，专硕 +6 moser233 2026-04-12 7/350 2026-04-12 22:11 by fqwang
[考研] 339求调剂 +8 hanwudada 2026-04-11 9/450 2026-04-12 15:36 by laoshidan
[考研] 331求调剂 +5 王国帅 2026-04-11 5/250 2026-04-11 22:56 by 溪涧流水
[考研] 270求调剂 +14 杨乐369 2026-04-11 14/700 2026-04-11 20:16 by 蓝云思雨
信息提示
关闭
请填处理意见
关闭