| 查看: 12084 | 回复: 22 | |||||||||||
| 【奖励】 本帖被评价16次，作者zdhlover增加金币 12.45 个 | |||||||||||
[资源]
【分享】用VASP进行Partial Charge分析实例
|
|||||||||||
|
[分享] 用VASP进行Partial Charge分析实例 VASP Version : 4.6 在这篇文章中，我将首先介绍Partial Charge的概念，以及如何用VASP具体的计算Partial Charge。首先，所谓的Partial Charge是针对与Total Charge来说的，指的是某个能量范围、某个K点或者某个特定的态所对应的电荷密度。在文献中最常见的是价带顶部，导带底部，表面态或者局域态所对应的Partial Charge。通过分析这些态所对应的Partial Charge，可以得到体系的一些性质，比如局域态具体的是局域在哪个原子上等。我将通过具体的例子说明如何用VASP进行Partial Charge Analysis。 进行Partial Charge Analysis的第一步是进行自洽的计算，得到体系的电子结构。这一步的计算采用通常的INCAR和KPOINTS文件。在自洽计算结束后，我们需要保存WAVECAR文件。（通过在INCAR文件中设置LWAVE=TRUE实现）在这个例子中，假设我们需要计算一个硅纳米线的导带和价带的Partial Charge。硅纳米线的结构如下： 第二步是画出能带结构，以决定你需要画哪条能带的那个K点的态所对应的Partial Charge。关于具体如何用VASP画能带，请参见用VASP4.6计算晶体硅能带实例一文。我们得到硅纳米线的能带结构如下： 画能带时有些小技巧。你可以用一些支持列模块的编辑器，如UltraEdit，将OUTCAR里的各个K点所对应的本征值粘贴到Origin中。这一步完成后，在Origin中做一个矩阵转置，然后将K点坐标贴到第一列，并将其设为X坐标。如此画出来的基本上就是能带图了。在Origin中可以通过设置纵轴范围来更加清楚的区分费米能级附近的各条能带。如上的硅纳米线所对应的能带结构图如下： 决定画哪条能带，或者那些感兴趣的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。这种模式最利于分析费米面附近的波函数的性质。 用第一种方法，我们可以得到硅纳米线价带顶部和导带底部的Partial Charge如下： [注意]此文的部分结果将发表于Phys. Rev. B |
回复此楼» 收录本帖的淘帖专辑推荐
 VASP | 第一性原理相关文档 | vasp | Photochemistry |
| VASP相关 | siesta & vasp | 概念教程等 | VASP计算方法 |
| VASP blabla |
» 猜你喜欢
人气不行了
已经有10人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有6人回复
-
我面上完蛋了
已经有10人回复
-
能否退出参与的面上项目解除限项
已经有24人回复
-
filecode，4个jtjc了
已经有17人回复
-
没有任何消息-是不是就凉了
已经有6人回复
-
2026国自然函评费到账
已经有19人回复
-
建议基金发布提前给出明确的时间点
已经有16人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有21人回复
-
估计是周四
已经有4人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 求vasp输出文件的分析
已经有5人回复
- 【求助】怎样用vasp计算charge transfer
已经有10人回复
- 关于用origin画Si能带图时特殊点纵轴的添加
已经有10人回复
- vasp 计算 partial charge 时k点选取的问题
已经有19人回复
- 三种charge density的区别？
已经有21人回复
- vasp画能带图的疑问
已经有12人回复
- charge的意义 高分求解
已经有5人回复
- vasp计算能带，刚开始就服务器就卡了
已经有20人回复
- vasp计算中，如何选择正确的K点
已经有7人回复
- 请问用vasp做charge density difference步骤是什么
已经有14人回复
- 能带计算 vasp
已经有11人回复
- 分享vasp的一些可用的免费资源（入门级）
已经有79人回复
- Bader charge analysis 资料
已经有684人回复
- VASP 算不准氧气键能
已经有16人回复
- 用vasp计算能带结构时报错
已经有5人回复
- 【求助】vasp计算Partial Charge，能带交叉的如何确定IBAND？
已经有13人回复
- 【求助】请教vasp如何使本征值按每条能带输出
已经有18人回复
3楼2009-06-30 13:56:29
4楼2009-06-30 16:10:37
5楼2009-06-30 16:33:11
6楼2009-07-01 20:12:32
7楼2009-09-01 05:10:29
8楼2009-09-01 07:48:58
10楼2009-09-05 11:25:05
11楼2010-03-06 15:19:24
12楼2010-03-22 08:59:45
|
请问楼主，我要计算边缘态的Partial Charge，能带交叉的怎么设置IBAND？这种情况是不是需要分段计算了！谢谢！ 详见： http://muchong.com/bbs/viewthread.php?tid=3022926&page=1#pid96405 |
13楼2011-03-30 19:27:48
14楼2011-03-31 13:53:13
15楼2011-03-31 15:12:22
16楼2012-06-23 15:12:19
17楼2014-02-21 10:41:06
18楼2014-08-29 16:34:31
19楼2015-03-10 15:59:17
★★★ 三星级,支持鼓励
|
你好 我在算partial charge时也遇到了问题：在自洽计算得到CHGCAR WAVECAR后 进行partial charge计算。输入 INCAR :INCAR: ISTART=1 ICHARG=11 ENCUT=400 NSW=0 LCHAG=.TRUE. ISPIN=2 LPARD=.TRUE. IBAND=80 81 NBMOD=2 KPUSE=1 LSEPB=.TRUE. LSEPK=.TRUE. KPOINTS: k-points along high symmetry lines 1 Reciprocal 0.000000 0.000000 0.000000 1.00 提交任务后错误提示： LDA part: xc-table for Ceperly-Alder, standard interpolation found WAVECAR, reading the header number of k-points has changed, file: 172 present: 1 trying to continue reading WAVECAR, but it might fail WARNING: stress and forces are not correct POSCAR, INCAR and KPOINTS ok, starting setup FFT: planning ...( 1 ) reading WAVECAR ERROR: while reading WAVECAR, plane wave coefficients changed 18395 |
20楼2015-03-21 13:38:10
21楼2017-10-20 10:53:10
23楼2018-12-12 20:01:05
简单回复
NKDMS2楼
2009-06-27 00:35
回复
2009-09-02 12:38
回复
周宾22楼
2017-10-23 08:47
回复
五星好评 顶 [ 发自手机版 http://muchong.com/3g ]
