求助: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埃的立方体。 返回小木虫查看更多
自己顶,不要沉啊
顶一下,我前天的算的结果是N多课0.31e。同样困惑。
不要沉
我遇到了跟你同样的问题,请问你现在解决了这个问题吗,