24小时热门版块排行榜    

查看: 2025  |  回复: 2

北漂2012

新虫 (初入文坛)

[交流] VASP parital charge 分波电荷密度计算中的一个问题 已有2人参与

VASP parital charge 分波电荷密度计算中的一个问题:
体系:一维的石墨烯纳米带。自洽计算选择的K点:1*20*1;因为体系在Y方向是周期的。
画能带如下选择K点61个,KPOINTS如下:
K-points
61
Line-mode
rec
0.0       0.0      0.0  !gamma
0.0      -0.5      0 .0  !Y
#############
通过观察能带图:
(1)我现在要画某个能量区间的电荷密度
按照官方说明书以及参考网上观点。我把自洽的WAVECAR和CHGCAR拷贝到新的目录下,当然还有POSCAR,POTCAR
(注意这时候我用的KPOINTS书,画能带的KPOINTS而不是自洽的KPOINTS !!,这一点说明书上没有特别说明,我也一直很困惑,到底用哪一个KPOINTS)
INCAR我设置如下
ISTART=1
ICHARG=1
NBMOD=-3
EINT=-0.1  0.1
LPARD=.TRUE.
###########
其他设置同自洽一致。
按理说,程序应该判断有多少k点,对应多少条能带上的能量落在我选择的能量区间上,然后把这些满足条件k点对应的电荷密度累加起来。
问题是:做partial charge计算时候,我的KPOINTS文件,应该用自洽的,还是用画能带的?如果用自洽的,我只有20个K点,太少了,可能没有一个k点对应的能量落在我选择的能量区间中。
(2)能带上对应的某条能带上某个K点,的态对应的电荷密度,用自洽的WAVECAR和CHGCAR,KPOINTS我用的是画能带的KPOINTS.
INCAR的设置如下:
ISTART=1
ICHARG=1
NBMOD=2
LPARD=.TRUE.
IBAND=95
KPUSE= 21  
LSEPB=.TRUE.
LSEPK=.TRUE.  
########
按照理解,我这里应该是用画能带时候的KPOINTS,或者自己把要计算的K点坐标直接写进KPOINTS。
但是我发现程序好像并没有按照我给的KPOINTS去计算,它依然是按照自洽的K点去计算的。
问题:在计算partial charge 我想得到某条能带,某个K点处的电荷密度(这个K点,刚好也不在我自洽的K点之中),程序会按照我此时自己提供的KPOINTS文件中的K点去计算吗?
回复此楼

» 猜你喜欢

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

» 本主题相关价值贴推荐,对您同样有帮助:

ibelieveicanfly,ibelieveeicantouchthesky.Ificanseeit,thenicandoit,ifijustbelieveit,thethereisnothingtoit.
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

菲顶顶

铜虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
楼主解决这个问题了吗?就是计算某一条能带上的某一个K点的电荷密度KPOINTS文件到底是如何设置的
2楼2015-03-25 19:42:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zh11111

银虫 (正式写手)

---


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
2楼: Originally posted by 菲顶顶 at 2015-03-25 19:42:51
楼主解决这个问题了吗?就是计算某一条能带上的某一个K点的电荷密度KPOINTS文件到底是如何设置的

我在计算parital charge时遇到了问题,请见我的求助帖http://muchong.com/bbs/viewthread.php?tid=8684252
谢谢!
---
3楼2015-05-26 08:53:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 北漂2012 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见