| 查看: 2292 | 回复: 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比较有意义。 大家觉得呢。 当然，如果事先确定一条带是由某个原子周围电子态产生的，把整个能带的合电荷密度画出来以证明这点，这是可以的。 |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有15人回复
-
澳大利亚 Murdoch University 全奖博士招生（3个名额）地质化工冶金领域
已经有20人回复
-
静等基金结果
已经有21人回复
-
FileCode能看出啥？
已经有28人回复
-
阿姨
已经有3人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有6人回复
-
filecode
已经有16人回复
-
重要来源：本周末出结果
已经有10人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有6人回复
-
不应该看fileCode
已经有9人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- vasp计算在进行静态计算时，出错。
已经有3人回复
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
1.?????????band?????K????о??????????????????????????????????????????????????????????????????????????????3eV???????band??????????2eV ?????????????????????У????????????????????????????????Χ??????????????????????????? ?????????????????????? > ??????????????????f?? ???^?@?? NBMOD= -3 ????????????????? band structure ?? partial charge decomposed ?? ?? band structure ?? partial charge ?? NBMOD = -2 ???^??? ???P? NBMOD = -3 ?????????e? ????~ ?????????? STM simulation |
9楼2016-06-29 19:04:47
蛋蛋小童鞋
版主 (著名写手)
- 应助: 71 (初中生)
- 贵宾: 0.058
- 金币: 13671.6
- 散金: 497
- 红花: 29
- 沙发: 1
- 帖子: 2070
- 在线: 753.8小时
- 虫号: 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
jugengfans
金虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 773.9
- 散金: 151
- 红花: 3
- 帖子: 537
- 在线: 205.6小时
- 虫号: 291609
- 注册: 2006-10-29
- 专业: 凝聚态物性 II ：电子结构
5楼2016-04-21 06:46:33
栾清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
binarec
木虫 (小有名气)
- 应助: 5 (幼儿园)
- 金币: 3512.3
- 红花: 4
- 帖子: 167
- 在线: 171.7小时
- 虫号: 893268
- 注册: 2009-11-04
- 性别: GG
- 专业: 金属材料表面科学与工程
7楼2016-06-29 00:24:36
jugengfans
金虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 773.9
- 散金: 151
- 红花: 3
- 帖子: 537
- 在线: 205.6小时
- 虫号: 291609
- 注册: 2006-10-29
- 专业: 凝聚态物性 II ：电子结构
8楼2016-06-29 02:36:08
jugengfans
金虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 773.9
- 散金: 151
- 红花: 3
- 帖子: 537
- 在线: 205.6小时
- 虫号: 291609
- 注册: 2006-10-29
- 专业: 凝聚态物性 II ：电子结构
|
谢谢回复。有人说计算band的时候即使是line mode，也会进行一个看不见的均匀撒点的。line mode里面设的K点不占权重。所以对费米能的影响仅仅是smearing=0和smearing=-5的区别 发自小木虫IOS客户端 |
10楼2016-06-30 08:20:04
