24小时热门版块排行榜    

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

洁如玉

铁虫 (小有名气)

[交流] 分享“Partial charge density计算-VASP” 已有3人参与

自己初学,受教很多。故抄过来,望能帮助到更多的新人。

文章出处:http://plianghust.blog.163.com/b ... 090822012892478779/

Partial charge density计算-VASP  

2012-09-09 14:47:08|  分类: 量子计算 |举报|字号 订阅
Partial charge density计算或称为Band decomposed charge density计算,即计算特定的某个(或某些)k点和本征值(这些k点和本征值是相互对应的)所对应的本征波函数的平方(也就是电荷密度)。特别是用在 STM的计算中,以及分析特定能量范围内或能量点的化学键特征(或atomic characteristic)。

LPARD:可赋予的值为.TRUE.或.FALSE.,它的默认值是.FALSE.,当为.TRUE.时,表示读入自洽收敛的CHGCAR和WAVECAR并进行Band decomposed charge density计算。
IBAND:它的赋值就是你想要计算的第几条能带或哪几条能带(比如要计算第4、5、6条能带,那么就设置IBAND = 4 5 6)。此时NBMOD的值就是所要计算的能带的条数(本例中,只计算一条能带,那么设置为1)。它和EINT不能一起用。
EINT:另一种指定所要计算能带的方式,它是指定计算某能量范围的Band decomposed charge density,一般是设置为两个实数,比如EINT= 4.00 5.00,此时NBOMD的值设置为-2。如果只设置了一个数,那么表示计算从EINT到费米能级这个范围内的,此时NBOMD的值为-3.
NBOMD: 当设置0时,表示计算全部的电荷密度(并包括了未占据的能带);如果设置-1,就是通常的电荷密度计算,这也是它的默认值。
KPUSE:指定所要计算的k点(哪个或哪几个),要计算一个,并相应的在KPOINTS设置所要计算的k点。如果要计算好几个,那么设置KPUSE的值为这些k点在KPOINTS中的序号。
LSEPB:指定是不是要把计算的partial charge density按每个带分别写到各自对应的文件PARCHG.nb.(设置为.TRUE.)中,还是把它们合并写到一个文件中(相当于把各个带对应的partial charge density加起来,设置为.FALSE.)。默认值为.FALSE.。
LSEPK:指定是不是把要计算的partial charge density按每个k分别写到各自对应的文件PARCHG.nk.(设置为.TRUE.)中,还是把它们合并写到一个文件中(相当于把各个带对应的partial charge density加起来,设置为.FALSE.)。默认值为.FALSE.。

决定画哪条能带,或者那些感兴趣的K点之后,有如下几种方法计算不同的Partial Charge。如果你希望计算价带顶端的Partial Charge,则需要首先通过能带结构图确定价带的能带标号。需要注意,进行Partial Charge分析必须要保留有自洽计算的WAVECAR才可以。

第一种Partial Charge分析的INCAR
ISTART =      1    job   : 0-new  1-cont  2-samecut
ICHARG =      1   charge: 1-file 2-atom 10-const

LPARD=.TRUE.
IBAND= 20 21 22 23
KPUSE= 1 2 3 4
LSEPB=.TRUE.
LSEPK=.TRUE.

这样的INCAR给出的是指定能带,指定K点所对应的Partial Charge。分析导带、价带等的Partial Charge特性,通常采用的都是这种模式。

第二种Partial Charge分析的INCAR
ISTART =      1    job   : 0-new  1-cont  2-samecut
ICHARG =      1   charge: 1-file 2-atom 10-const

LPARD=.TRUE.
EINT = -10.3 -5.1
LSEPB=.FALSE.
LSEPK=.FALSE.

这样的INCAR给出的是在[-10.3 -5.1]能量之间的Partial Charge。这种模式适合于分析某个能量区间内的波函数的性质。

第三种Partial Charge分析的INCAR
ISTART =      1    job   : 0-new  1-cont  2-samecut
ICHARG =      1   charge: 1-file 2-atom 10-const

LPARD=.TRUE.
NBMOD=-3
EINT = -1
LSEPB=.FALSE.
LSEPK=.FALSE.

这样的INCAR给出的是从[Ef-1.0 Ef]能量之间的Partial Charge。这种模式最利于分析费米面附近的波函数的性质。
回复此楼

» 收录本帖的淘帖专辑推荐

电荷密度分离

» 本帖已获得的红花(最新10朵)

» 猜你喜欢

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

加油少年!!!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小木虫: 金币+0.5, 给个红包,谢谢回帖
送红花一朵
本帖内容被屏蔽

4楼2021-09-16 21:24:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

yushi_xing

新虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
楼主求助,怎么赋予的电荷,譬如硫酸根离子,那个电荷怎么赋予的呢

发自小木虫Android客户端
3楼2016-07-18 16:29:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 什么时候开奖? +4 CrisMessi 2026-08-18 4/200 2026-08-19 00:53 by gatelove
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 27/1350 2026-08-18 21:44 by Noways
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +3 瞬息宇宙 2026-08-15 3/150 2026-08-18 21:40 by xiaohuixin
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 9/450 2026-08-18 20:09 by 奶牛小黑
[基金申请] 朋友圈看到的 +4 wangzilk 2026-08-18 5/250 2026-08-18 19:25 by wangzilk
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 时间戳变了,能看出什么问题? +8 基诺咪客 2026-08-17 11/550 2026-08-18 17:13 by 悠然熊猫
[基金申请] 重要消息,中午系统在维护 +9 yuleib84 2026-08-18 10/500 2026-08-18 16:52 by archvillain
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +4 Tide man 2026-08-14 5/250 2026-08-18 11:35 by Tide man
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
信息提示
请填处理意见