24小时热门版块排行榜    

查看: 584  |  回复: 0

沐水公子

木虫 (初入文坛)

[求助] vasp计算电荷密度差成键处为负值?

计算电荷密度差,原子部分的根据大家的指示,将ICHARG改成12,将NELM改为0
然后用chgdiff.pl脚本操作./chgdiff CHGCAR1(原子) CHGCAR2(物质)
得到的CHGCAR_diff有的是对的如Si,但有的结果成键处是负值,与文献相反
是否有大神知道是怎么回事,检查过基本排除是CHGCAR相减正犯错误的问题,不然也不会有的得到是正确的
求各位指教!!!!!
现附上得到的结果与文献对比
1、Si ,1-1图为计算,1-2图为文献,成键处是正值,黄色
2、GeTe,2-1图为计算,2-2图为文献,我计算的成键处是负值,蓝色;文献中是黄色,正值
请问为什么会是反的啊啊啊啊!!!!
vasp计算电荷密度差成键处为负值?
1-1、Si-计算.jpg


vasp计算电荷密度差成键处为负值?-1
1-2、Si-文献.jpg


vasp计算电荷密度差成键处为负值?-2
2-1、GeTe-计算.jpg


vasp计算电荷密度差成键处为负值?-3
2-2、GeTe-文献.jpg
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

相关版块跳转 我要订阅楼主 沐水公子 的主题更新
信息提示
请填处理意见