| 查看: 5989 | 回复: 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朵）
» 猜你喜欢
求助
已经有4人回复
-
0805材料320求调剂
已经有3人回复
-
297求调剂
已经有11人回复
-
268求调剂
已经有9人回复
-
一志愿西安交通大学材料工程专业 282分求调剂
已经有10人回复
-
278求调剂
已经有7人回复
-
26考研一志愿中国石油大学(华东)305分求调剂
已经有5人回复
-
280求调剂
已经有9人回复
-
279求调剂
已经有5人回复
-
286求调剂
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于vasp计算导带底或价带顶电荷密度的疑惑(partial charge density)
已经有4人回复
- vasp 计算 partial charge 时k点选取的问题
已经有19人回复
- 【求助】vasp计算Partial Charge，能带交叉的如何确定IBAND？
已经有13人回复
- 【求助】两种计算电荷密度的方法得到的结果相差很大
已经有9人回复
- 【分享】用VASP进行Partial Charge分析实例
已经有22人回复
送红花一朵
 |
2楼2014-06-12 22:02:01
yushi_xing
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 8.3
- 红花: 1
- 帖子: 95
- 在线: 16.8小时
- 虫号: 4822415
- 注册: 2016-07-06
- 性别: MM
- 专业: 海洋工程
3楼2016-07-18 16:29:47
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
送红花一朵
小木虫: 金币+0.5, 给个红包，谢谢回帖
送红花一朵 |
本帖内容被屏蔽 |
4楼2021-09-16 21:24:30
40