| 查看: 1928 | 回复: 10 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
jugengfans金虫 (正式写手)
|
[交流]
vasp中画Partial Charge的问题！感觉比较严重 已有6人参与
|
||
|
partial charge有三种画法，其中一种是设置一个能量范围，画出这个能量范围内的电荷密度。这个能量是从你设置的值到费米能级的，如： 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. 上面的init参数就表示是画从费米能下-1eV开始，到费米能的全部电子态的电荷密度。 问题来了！ 1.因为我们算band的时候，K点并没有均匀撒在全空间，所以算能带时得到的费米能级，和自洽计算的时候是不同的！ 比如费米能在自洽计算的时候是3eV，在计算band的时候就变成了2eV 那么，在这种画能带的方法中，程序肯定是按照错误的费米能级来给能量范围的，这样画出来不就没意义了吗！！ 那这个方法有什么存在价值呢。 2.另外，突然觉得某条能带所对应的电荷密度也是有问题的。因为它就是把你在Kpoints里面列的那几个点所对应的电荷密度加起来了！ 那如果你的line mode路径偏重于某一个方向，重复走了好几次，而某些方向走了一次，那这个时候总电荷密度就会偏向体现重复走的那些路径！ 这样一来，感觉只有分析能带N上的单个K点所对应的电荷密度才有意义。 也就是说，LSEPB=.FALSE, LSEPK=.FALSE 让程序把每个K点单独存放成CHGCAR比较有意义。 大家觉得呢。 当然，如果事先确定一条带是由某个原子周围电子态产生的，把整个能带的合电荷密度画出来以证明这点，这是可以的。 |
回复此楼» 猜你喜欢
参与限项
已经有3人回复
-
假如你的研究生提出不合理要求
已经有7人回复
-
实验室接单子
已经有4人回复
-
全日制（定向）博士
已经有4人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
不自信的我
已经有12人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
-
北核录用
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- vasp计算在进行静态计算时，出错。
已经有3人回复
binarec
木虫 (小有名气)
- 应助: 5 (幼儿园)
- 金币: 3512.3
- 红花: 4
- 帖子: 167
- 在线: 171.7小时
- 虫号: 893268
- 注册: 2009-11-04
- 性别: GG
- 专业: 金属材料表面科学与工程
7楼2016-06-29 00:24:36
蛋蛋小童鞋
版主 (著名写手)
- 应助: 71 (初中生)
- 贵宾: 0.058
- 金币: 13663.4
- 散金: 497
- 红花: 29
- 沙发: 1
- 帖子: 2041
- 在线: 753.1小时
- 虫号: 4083004
- 注册: 2015-09-18
- 专业: 理论和计算化学
- 管辖: 第一性原理
2楼2016-04-15 08:24:28
jugengfans
金虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 773.9
- 散金: 151
- 红花: 3
- 帖子: 537
- 在线: 205.6小时
- 虫号: 291609
- 注册: 2006-10-29
- 专业: 凝聚态物性 II ：电子结构
3楼2016-04-15 18:33:17
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
本帖内容被屏蔽 |
4楼2016-04-16 10:28:26