24小时热门版块排行榜

>论坛更新日志 (4546)
>文献求助 (475)
>导师招生 (389)
>虫友互识 (350)
>考博 (167)
>论文投稿 (159)
>博后之家 (139)
>休闲灌水 (131)
>硕博家园 (116)
>考研 (96)
>公派出国 (92)
>基金申请 (79)
>论文道贺祈福 (71)
>教师之家 (71)
>找工作 (54)
>招聘信息布告栏 (45)

返回列表

低调de往前走

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 283.5
帖子: 6
在线: 3.8小时
虫号: 3203175
注册: 2014-05-13
性别: GG
专业: 凝聚态物性 II ：电子结构

[求助] 求助：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埃的立方体。

回复此楼