| 查看: 2077 | 回复: 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比较有意义。 大家觉得呢。 当然，如果事先确定一条带是由某个原子周围电子态产生的，把整个能带的合电荷密度画出来以证明这点，这是可以的。 |
回复此楼» 猜你喜欢
284求调剂
已经有10人回复
-
一志愿山东大学药学学硕求调剂
已经有4人回复
-
07化学280分求调剂
已经有4人回复
-
298-一志愿中国农业大学-求调剂
已经有12人回复
-
求材料，环境专业调剂
已经有3人回复
-
335求调剂
已经有5人回复
-
求调剂
已经有7人回复
-
一志愿吉大化学322求调剂
已经有4人回复
-
环境学硕288求调剂
已经有8人回复
-
341求调剂(一志愿湖南大学070300)
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- vasp计算在进行静态计算时，出错。
已经有3人回复
栾清luanqing
铜虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 416.8
- 散金: 60
- 帖子: 23
- 在线: 22.7小时
- 虫号: 1241621
- 注册: 2011-03-22
- 性别: MM
- 专业: 原子和分子物理
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
那请问，K点怎样才能在K空间采样完全？我最近计算的占据数就出现了误差，麻烦请楼主给看看，请不吝赐教！http://muchong.com/bbs/viewthread.php?tid=10442146&fpage=1 |
6楼2016-06-28 20:44:53
蛋蛋小童鞋
版主 (著名写手)
- 应助: 71 (初中生)
- 贵宾: 0.058
- 金币: 13668.4
- 散金: 497
- 红花: 29
- 沙发: 1
- 帖子: 2066
- 在线: 753.4小时
- 虫号: 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
