24小时热门版块排行榜    

查看: 3026  |  回复: 19
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

前田庆次

木虫 (著名写手)

[求助] 求救:愁人的bader電荷分析……大侠救命

本来是想复现一下文献的情景,石墨烯吸附金属Sn原子,顺便做了一下bader电荷分析,这下可不得了,T,B,H位的吸附,结果都很诡异……实在是感觉很难解释,使用的赝势是Sn_d,PBE泛函,下面直接上bader电荷分析得到的ACF.dat
T位:
#         X           Y           Z        CHARGE     MIN DIST    ATOMIC VOL
--------------------------------------------------------------------------------
    1      2.8516      0.0000     14.9709      5.8719      0.6070     23.5489
    2      1.4260      0.0000     14.9714      5.7519      0.6067     23.8618
    3      4.9873     -1.2382     14.9815      6.0269      0.6111     22.6046
    4      3.5641     -1.2353     14.9659      6.1072      0.6078     24.4330
    5      7.1294     -2.4750     14.9815      6.0683      0.6505     22.8459
    6      5.7029     -2.4722     14.9671      5.9549      0.6081     22.3019
    7      9.2696     -3.7055     14.9709      5.9455      0.6070     23.9497
    8      7.8435     -3.7061     14.9659      6.0589      0.6606     23.9458
    9      4.9873      1.2382     14.9815      6.0269      0.6111     22.5968
   10      3.5641      1.2353     14.9659      6.1072      0.6078     24.4339
   11      5.6958      0.0000      0.0228      6.1554      0.6149     18.6998
   12      9.2741     -1.2368     14.9815      6.1445      0.6111     23.3988
   13      7.8475     -1.2423      0.0228      6.0355      0.6149     18.1314
   14     11.4088     -2.4694     14.9565      5.7609      0.6057     23.4011
   15      9.9832     -2.4707     14.9659      6.0254      0.6078     23.6145
   16      7.1294      2.4750     14.9815      6.0683      0.6505     22.8459
   17      5.7029      2.4722     14.9671      5.9549      0.6081     22.3019
   18      9.2741      1.2368     14.9815      6.1445      0.6111     23.3992
   19      7.8475      1.2423      0.0228      6.0355      0.6149     18.1314
   20     11.4087      0.0000     14.9520      6.0695      0.6065     25.7021
   21      9.9850      0.0000     14.9671      5.8643      0.6081     21.5529
   22     13.5474     -1.2348     14.9520      5.9383      0.6061     24.3916
   23     12.1210     -1.2357     14.9517      6.0805      0.6608     25.2588
   24      9.2696      3.7055     14.9709      5.9455      0.6070     23.9497
   25      7.8435      3.7061     14.9659      6.0589      0.6606     23.9458
   26     11.4088      2.4694     14.9565      5.7609      0.6057     23.3990
   27      9.9832      2.4707     14.9659      6.0254      0.6078     23.6143
   28     13.5474      1.2348     14.9520      5.9383      0.6061     24.3916
   29     12.1210      1.2357     14.9517      6.0805      0.6608     25.2590
   30     15.6859      0.0000     14.9565      5.7991      0.6057     23.9656
   31     14.2614      0.0000     14.9517      5.8728      0.6067     23.9852
   32      7.1303      0.0000      0.0507      6.2647      0.6604     15.8271
   33      7.1303      0.0000      2.7788     48.0338      1.3418    534.3591
--------------------------------------------------------------------------------
    VACUUM CHARGE:               0.0000
    VACUUM VOLUME:               0.0000
    NUMBER OF ELECTRONS:       239.9773

第33个原子是Sn原子,其余的都是C原子,怎么总电子数是239.9773呢?怎么出来了小数?另外Sn是第五十号元素,怎么出来了48个电子?绝对不可能是转移了2个电子到石墨烯上,再来看B位的:
    #         X           Y           Z        CHARGE     MIN DIST    ATOMIC VOL
--------------------------------------------------------------------------------
    1      2.8514      0.0010     14.9563      5.8738      0.6059     24.0667
    2      1.4252      0.0015     14.9591      5.9500      0.6073     24.8776
    3      4.9870     -1.2358     14.9714      5.8895      0.6095     21.9642
    4      3.5640     -1.2346     14.9500      6.1247      0.6073     24.9959
    5      7.1258     -2.4781      0.0100      6.0652      0.6305     20.3430
    6      5.7015     -2.4732     14.9714      6.1060      0.6095     23.4929
    7      9.2695     -3.7052     14.9743      6.0564      0.6064     23.7231
    8      7.8423     -3.7076     14.9772      5.9370      0.6075     22.9285
    9      4.9884      1.2355     14.9772      5.8684      0.6075     22.4948
   10      3.5639      1.2346     14.9595      6.0923      0.6077     25.3325
   11      5.6950      0.0002      0.0100      6.1628      0.6154     20.9643
   12      9.2786     -1.2352      0.0100      6.1628      0.6154     20.9643
   13      7.8463     -1.2402      0.0654      6.0734      0.6570     15.3270
   14     11.4097     -2.4696     14.9595      5.9590      0.6077     24.2475
   15      9.9852     -2.4705     14.9772      5.8684      0.6075     22.4936
   16      7.1313      2.4726     14.9772      5.9370      0.6075     22.9291
   17      5.7041      2.4702     14.9743      6.0564      0.6064     23.7231
   18      9.2721      1.2382     14.9714      6.1060      0.6095     23.4929
   19      7.8478      1.2431      0.0100      6.0652      0.6305     20.3430
   20     11.4096     -0.0004     14.9500      6.1247      0.6073     24.9950
   21      9.9866      0.0008     14.9714      5.8895      0.6095     21.9640
   22     13.5484     -1.2365     14.9591      5.9500      0.6073     24.8778
   23     12.1221     -1.2360     14.9563      5.8738      0.6059     24.0667
   24      9.2689      3.7061     14.9563      5.9286      0.6059     24.5125
   25      7.8443      3.7059     14.9595      6.0882      0.6596     25.1661
   26     11.4089      2.4692     14.9424      5.9288      0.6065     24.4923
   27      9.9826      2.4711     14.9500      6.0396      0.6073     24.4588
   28     13.5473      1.2339     14.9500      6.0396      0.6073     24.4578
   29     12.1210      1.2358     14.9424      5.7575      0.6052     23.2667
   30     15.6856     -0.0009     14.9595      5.9036      0.6062     23.9125
   31     14.2610     -0.0011     14.9563      5.9286      0.6059     24.5125
   32      7.1273      0.0052      0.0654      6.2068      0.6570     16.0987
   33      7.4868     -0.6175      2.7575     77.8800      1.3584    532.5623
--------------------------------------------------------------------------------
    VACUUM CHARGE:               0.0000
    VACUUM VOLUME:               0.0000
    NUMBER OF ELECTRONS:       269.8937
电子总数居然变成了269.8937……,此时Sn的电子数理解为多少呢?78?那么就是失去了0.12个电子?
这时我已经抓狂了,但是更为难理解的是,H位的结果:
    #         X           Y           Z        CHARGE     MIN DIST    ATOMIC VOL
--------------------------------------------------------------------------------
    1      2.8521      0.0000     14.9956      5.8922      0.6061     24.1005
    2      1.4263      0.0000     14.9948      5.7997      0.6059     23.0273
    3      4.9906     -1.2353     14.9942      6.0374      0.6066     23.2391
    4      3.5650     -1.2348     14.9948      5.9258      0.6062     23.5879
    5      7.1290     -2.4721     14.9941      6.0428      0.6086     22.1198
    6      5.7037     -2.4703     14.9942      5.9587      0.6066     22.4406
    7      9.2694     -3.7056     14.9942      5.9322      0.6061     22.3256
    8      7.8433     -3.7056     14.9942      6.0655      0.6603     23.3402
    9      4.9906      1.2353     14.9942      6.0374      0.6066     23.2393
   10      3.5650      1.2348     14.9948      5.9258      0.6062     23.5875
   11      5.7018      0.0000     14.9941      6.1274      0.6086     22.5912
   12      9.2708     -1.2375      0.0012      6.1289      0.6075     19.5602
   13      7.8419     -1.2375      0.0012      6.2640      0.6616     20.4372
   14     11.4090     -2.4703     14.9942      5.9595      0.6066     22.4889
   15      9.9836     -2.4721     14.9941      6.0428      0.6086     22.1198
   16      7.1290      2.4721     14.9941      6.0428      0.6086     22.1198
   17      5.7037      2.4703     14.9942      5.9587      0.6066     22.3667
   18      9.2708      1.2375      0.0012      6.1289      0.6075     19.5860
   19      7.8419      1.2375      0.0012      6.2632      0.6616     20.4260
   20     11.4109      0.0000     14.9941      6.1274      0.6086     22.5912
   21      9.9853      0.0000      0.0012      5.8953      0.6075     18.2777
   22     13.5476     -1.2348     14.9948      5.9258      0.6062     23.5879
   23     12.1220     -1.2353     14.9942      6.0374      0.6066     23.2395
   24      9.2694      3.7056     14.9942      5.9322      0.6061     22.2453
   25      7.8433      3.7056     14.9942      6.0655      0.6603     23.3402
   26     11.4090      2.4703     14.9942      5.9595      0.6066     22.4156
   27      9.9836      2.4721     14.9941      6.0428      0.6086     22.1198
   28     13.5476      1.2348     14.9948      5.9258      0.6062     23.5877
   29     12.1220      1.2353     14.9942      6.0374      0.6066     23.2389
   30     15.6864      0.0000     14.9948      5.8005      0.6059     23.1192
   31     14.2605      0.0000     14.9956      5.8922      0.6061     24.1005
   32      7.1274      0.0000      0.0012      5.8953      0.6075     18.2777
   33      8.5563      0.0000      3.1380     57.0072      1.6683    555.2029
--------------------------------------------------------------------------------
    VACUUM CHARGE:               0.0000
    VACUUM VOLUME:               0.0000
    NUMBER OF ELECTRONS:       249.0782
电子总数又变了,居然成了249.0782……Sn的原子总数又变成了57.0072……

想知道,怎么会出现这样奇怪的结果……
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

长枪一横花飘零,松风追月伴我行。无双人间世如梦,倾奇万世永留名。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

前田庆次

木虫 (著名写手)

引用回帖:
16楼: Originally posted by lfhuang at 2013-12-04 20:18:38
Mulliken, or Löwdin population
这两者用的都是自然的原子轨道,不像Bader,太artificial了,很多结果太出乎物理之外了,虽然数学上看貌似挺漂亮的。...

不好意思,大俠,我是vasp,好像做不了Lowdin,能做Mullken嗎?應該如何做呢?小木蟲上看了一圈沒找到,請指教
长枪一横花飘零,松风追月伴我行。无双人间世如梦,倾奇万世永留名。
17楼2013-12-04 21:04:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 20 个回答

ultra1989

金虫 (小有名气)

【答案】应助回帖


感谢参与,应助指数 +1
liliangfang: 金币+1, 谢谢交流 2013-12-03 07:42:25
有小数是正常的,你看下bader的算法就知道了,暴力分割的结果
不畏浮云遮望眼
2楼2013-12-02 22:52:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ultra1989

金虫 (小有名气)

另外你确定这几个体系用的是同一个赝势?
不畏浮云遮望眼
3楼2013-12-02 22:56:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

前田庆次

木虫 (著名写手)

引用回帖:
3楼: Originally posted by ultra1989 at 2013-12-02 22:56:35
另外你确定这几个体系用的是同一个赝势?

绝对采用同一个,都是Sn_d,而且您可以发现这三个体系的总电子数都不相同
长枪一横花飘零,松风追月伴我行。无双人间世如梦,倾奇万世永留名。
4楼2013-12-02 22:58:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见