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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » MS » 【求助】Dmol3中hirshfeld电荷布局分布的问题
51/1返回列表
查看: 1441  |  回复: 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的回帖

topten-gg

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

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

qscao

木虫 (小有名气)
★ ★
小木虫(金币+0.5):恭喜抢沙发，给个红包
aylayl08(金币+1,VIP+0):谢谢参与 10-30 17:36
小数表示该原子带的电荷
赞一下(7人) 回复此楼
2楼2009-10-30 15:23:07
 已阅   回复此楼   关注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的回帖
查看全部 6 个回答
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭