24小时热门版块排行榜    

查看: 1658  |  回复: 6

低调de往前走

新虫 (初入文坛)

[求助] 求助:Vasp计算Bader电荷问题,Ben-NH2中氨基的净Bader电荷计算值与理论相反

用Vasp计算氨基苯(Ben-NH2)中氨基的Bader电荷,结果如下:
ACF.dat文件:      
6、7是氨基上的两个H原子,14是N原子,所以氨基现在的电荷是0.4784+0.4644+6.4067=7.3495
减去氨基的7个价电子,也就是氨基多了0.3495电子,而理论上苯环上的氨基应该是电子给体
为什么会多了电子呢,请大神帮忙看一看。
    #         X                Y              Z          CHARGE     MIN DIST    ATOMIC VOL
--------------------------------------------------------------------------------
    1      7.8081       9.4480      9.5463      0.8880      0.2912           764.7257
    2     12.8020      9.4859      9.5165      0.9014      0.2992          1048.0155
    3     11.5513     11.6320     9.6024      0.8547      0.1684          985.8867
    4     11.5729      7.2986      9.4594      0.8499      0.1749         1326.9372
    5      9.0634       7.3036      9.4604      0.8422      0.2288         1212.0266
    6      7.9297     11.8769      9.3592      0.4784      0.0708        1169.3560
    7      9.3898     12.7242      9.3840      0.4644      0.0291        823.1664
    8      8.9071       9.4561      9.5406      4.0407      0.4981         72.2097
    9      9.5924     10.6888      9.5860      3.8599      0.3285         40.0318
   10     11.0037   10.6795      9.5710      3.9833      0.5308         69.6469
   11     11.7047     9.4694      9.5267      4.0632      0.4956        75.2214
   12     11.0200     8.2446      9.4925      4.2242      0.5326        86.0362
   13      9.6165      8.2507      9.4962      4.1429      0.5510        69.8123
   14      8.8857     11.9011      9.7247     6.4067      0.6562        256.9276
--------------------------------------------------------------------------------
    VACUUM CHARGE:               0.0000
    VACUUM VOLUME:               0.0000
    NUMBER OF ELECTRONS:        36.0000

INCAR文件:
SYSTEM = NH2-BEN
ISTART = 0
ICHARG = 2
NELM = 100         
NELMIN = 2  
#LCHARG = F
#LWAVE= F
PREC = Medium
GGA = PE
ENCUT = 300
EDIFF = 1E-4
EDIFFG = -0.02
NSW = 1000
ISYM = 2
NPAR = 1
ISMEAR = -5
IBRION = 2
#ISPIN = 2
#MAGMOM = 48*0 46*0 10  10
#LDAU = .TURE.
#LDAUTYPE = 2
#LDAUL = -1 -1 3
#LDAUU = 0 0 7.4
#LDAUJ = 0 0 0.5
#LDAUPRINT = 2      
LORBIT = 10
LREAL = Auto
ISIF = 2
#NGX = 58  
#NGY = 54  
#NGZ = 90

结构是在MS中构建的,Ben-NH2分子的晶格是20埃的立方体。
回复此楼
苦心人天不负
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

低调de往前走

新虫 (初入文坛)

自己顶,不要沉啊
苦心人天不负
2楼2017-08-06 22:45:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

南墙2014

木虫 (正式写手)

顶一下,我前天的算的结果是N多课0.31e。同样困惑。
3楼2017-08-07 15:32:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

低调de往前走

新虫 (初入文坛)

引用回帖:
3楼: Originally posted by 南墙2014 at 2017-08-07 15:32:43
顶一下,我前天的算的结果是N多课0.31e。同样困惑。

搞不明白什么原因

发自小木虫IOS客户端
苦心人天不负
4楼2017-08-09 08:42:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

低调de往前走

新虫 (初入文坛)

苦心人天不负
5楼2017-08-10 22:30:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

低调de往前走

新虫 (初入文坛)

苦心人天不负
6楼2017-08-14 08:55:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

鱼儿Pamela

铁虫 (初入文坛)

我遇到了跟你同样的问题,请问你现在解决了这个问题吗
7楼2020-12-27 14:09:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 低调de往前走 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 Qj98d49dHFb9 2026-08-25 4/200 2026-08-26 03:08 by cNXvBfCpiZOM
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 7K1CJE38xLG4 2026-08-25 4/200 2026-08-26 02:04 by cNXvBfCpiZOM
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 19/950 2026-08-25 17:57 by chick875
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 22/1100 2026-08-25 15:36 by 医学老男孩
[基金申请] 没有任何消息-是不是就凉了 +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
[基金申请] 今天基金会出结果吗?20260819 +17 kkkl_v 2026-08-19 18/900 2026-08-25 09:41 by windflowerwy
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 范进中举一文的中心思想 +7 炎黄贵胄 2026-08-22 8/400 2026-08-25 08:48 by ZJTJZ
[基金申请] 只有每年这种时候来逛逛小木虫 +25 yaoyewhu2008 2026-08-20 27/1350 2026-08-25 08:01 by Equinoxhua
[基金申请] 能否退出参与的面上项目解除限项 +21 koalala 2026-08-24 24/1200 2026-08-24 19:25 by 家与远方
[基金申请] filecode,4个jtjc了 +14 ziyangfang 2026-08-19 17/850 2026-08-24 18:37 by 哈哈蛤?
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
信息提示
请填处理意见