应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » HSE（GW类似）计算能带总结
51/1返回列表
查看: 20128  |  回复: 51
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【奖励】 本帖被评价17次，作者liliangfang增加金币 11.6 个
本帖产生 1 个 1ST强帖 ，点击这里进行查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

liliangfang

荣誉版主 (著名写手)

[资源] HSE（GW类似）计算能带总结

1.DFT
INCAR
SYSTEM = Si-Diamond
ENCUT = 400
ISTART = 0
ISYM = 0
ICHARG = 2
ISMEAR = 0
SIGMA = 0.05
LORBIT = 11
NSW = 0
IBRION = -1
ISIF = 2
EDIFF = 0.00001
EDIFFG = -0.001

KPOINTS
Automatic mesh
0
G
4   4   4
0   0   0

POSCAR
Si-Diamond:                             
   5.48000000000000     
     0.0000000000000000    0.5000000000000000    0.5000000000000000
     0.5000000000000000    0.0000000000000000    0.5000000000000000
     0.5000000000000000    0.5000000000000000    0.0000000000000000
   Si
   2
Direct
  0.8750000000000000  0.8750000000000000  0.8750000000000000
  0.1250000000000000  0.1250000000000000  0.1250000000000000

  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
vasp计算
2.HSE06
INCAR
设置ICHARG =2 ，IBRION =1 ,NELMIN =5 or 其他值
SYSTEM = Si-Diamond
ISTART = 1
ICHARG = 2
EDIFF = 0.00001
EDIFFG = -0.001
ENCUT = 400
ENAUG = 800
LREAL = .FALSE.
LWAVE = .TRUE.
LCHARG = .TRUE.
NELM = 200
NSW = 0
IBRION = -1
LMAXMIX = 4
ISMEAR = 0
SIGMA = 0.1
NSIM = 4
ISYM = 0
LHFCALC = .TRUE.
HFSCREEN = 0.2
ALGO = D
TIME = 0.4
ENCUTFOCK = 0
AEXX = 0.25
其他文件从DFT文件夹拷贝过来，vasp计算
3.band计算
INCAR
SYSTEM = Si-Diamond
ISTART = 1
ICHARG = 2
EDIFF = 0.00001
EDIFFG = -0.001
ENCUT = 400
ENAUG = 800
LREAL = .FALSE.
LWAVE = .TRUE.
LCHARG = .TRUE.
NELM = 200
NSW = 0
IBRION = 1
LMAXMIX = 4
ISMEAR = 0
SIGMA = 0.1
NSIM = 4
ISYM = 0
LHFCALC = .TRUE.
HFSCREEN = 0.2
ALGO = N
TIME = 0.4
ENCUTFOCK = 0
AEXX = 0.25
NELMIN = 5

KPOINTS
从HSE06文件夹的IBZKPT拷贝到KPOINTS,然后把高对称点拷贝在后面，权重为0
更新k点总数
Automatically generated mesh
      127
Reciprocal lattice
    0.00000000000000    0.00000000000000    0.00000000000000             1
    0.25000000000000    0.00000000000000    0.00000000000000             2
    0.50000000000000    0.00000000000000    0.00000000000000             1
    0.00000000000000    0.25000000000000    0.00000000000000             2
    0.25000000000000    0.25000000000000    0.00000000000000             2
    0.50000000000000    0.25000000000000    0.00000000000000             2
   -0.25000000000000    0.25000000000000    0.00000000000000             2
    0.00000000000000    0.50000000000000    0.00000000000000             1
    0.25000000000000    0.50000000000000    0.00000000000000             2
    0.50000000000000    0.50000000000000    0.00000000000000             1
    0.00000000000000    0.00000000000000    0.25000000000000             2
    0.25000000000000    0.00000000000000    0.25000000000000             2
    0.50000000000000    0.00000000000000    0.25000000000000             2
   -0.25000000000000    0.00000000000000    0.25000000000000             2
    0.00000000000000    0.25000000000000    0.25000000000000             2
    0.25000000000000    0.25000000000000    0.25000000000000             2
    0.50000000000000    0.25000000000000    0.25000000000000             2
   -0.25000000000000    0.25000000000000    0.25000000000000             2
    0.00000000000000    0.50000000000000    0.25000000000000             2
    0.25000000000000    0.50000000000000    0.25000000000000             2
    0.50000000000000    0.50000000000000    0.25000000000000             2
   -0.25000000000000    0.50000000000000    0.25000000000000             2
    0.00000000000000   -0.25000000000000    0.25000000000000             2
    0.25000000000000   -0.25000000000000    0.25000000000000             2
    0.50000000000000   -0.25000000000000    0.25000000000000             2
   -0.25000000000000   -0.25000000000000    0.25000000000000             2
    0.00000000000000    0.00000000000000    0.50000000000000             1
    0.25000000000000    0.00000000000000    0.50000000000000             2
    0.50000000000000    0.00000000000000    0.50000000000000             1
    0.00000000000000    0.25000000000000    0.50000000000000             2
    0.25000000000000    0.25000000000000    0.50000000000000             2
    0.50000000000000    0.25000000000000    0.50000000000000             2
   -0.25000000000000    0.25000000000000    0.50000000000000             2
    0.00000000000000    0.50000000000000    0.50000000000000             1
    0.25000000000000    0.50000000000000    0.50000000000000             2
    0.50000000000000    0.50000000000000    0.50000000000000             1
  0.500000  0.000000  0.500000  0
  0.475000  0.000000  0.475000  0
  0.450000  0.000000  0.450000  0
  0.425000  0.000000  0.425000  0
  0.400000  0.000000  0.400000  0
  0.375000  0.000000  0.375000  0
  0.350000  0.000000  0.350000  0
  0.325000  0.000000  0.325000  0
  0.300000  0.000000  0.300000  0
  0.275000  0.000000  0.275000  0
  0.250000  0.000000  0.250000  0
  0.225000  0.000000  0.225000  0
  0.200000  0.000000  0.200000  0
  0.175000  0.000000  0.175000  0
  0.150000  0.000000  0.150000  0
  0.125000  0.000000  0.125000  0
  0.100000  0.000000  0.100000  0
  0.075000  0.000000  0.075000  0
  0.050000  0.000000  0.050000  0
  0.025000  0.000000  0.025000  0
  0.000000  0.000000  0.000000  0
  0.025000  0.025000  0.025000  0
  0.050000  0.050000  0.050000  0
  0.075000  0.075000  0.075000  0
  0.100000  0.100000  0.100000  0
  0.125000  0.125000  0.125000  0
  0.150000  0.150000  0.150000  0
  0.175000  0.175000  0.175000  0
  0.200000  0.200000  0.200000  0
  0.225000  0.225000  0.225000  0
  0.250000  0.250000  0.250000  0
  0.275000  0.275000  0.275000  0
  0.300000  0.300000  0.300000  0
  0.325000  0.325000  0.325000  0
  0.350000  0.350000  0.350000  0
  0.375000  0.375000  0.375000  0
  0.400000  0.400000  0.400000  0
  0.425000  0.425000  0.425000  0
  0.450000  0.450000  0.450000  0
  0.475000  0.475000  0.475000  0
  0.500000  0.500000  0.500000  0
  0.500000  0.487500  0.512500  0
  0.500000  0.475000  0.525000  0
  0.500000  0.462500  0.537500  0
  0.500000  0.450000  0.550000  0
  0.500000  0.437500  0.562500  0
  0.500000  0.425000  0.575000  0
  0.500000  0.412500  0.587500  0
  0.500000  0.400000  0.600000  0
  0.500000  0.387500  0.612500  0
  0.500000  0.375000  0.625000  0
  0.500000  0.362500  0.637500  0
  0.500000  0.350000  0.650000  0
  0.500000  0.337500  0.662500  0
  0.500000  0.325000  0.675000  0
  0.500000  0.312500  0.687500  0
  0.500000  0.300000  0.700000  0
  0.500000  0.287500  0.712500  0
  0.500000  0.275000  0.725000  0
  0.500000  0.262500  0.737500  0
  0.500000  0.250000  0.750000  0
  0.487500  0.262500  0.750000  0
  0.475000  0.275000  0.750000  0
  0.462500  0.287500  0.750000  0
  0.450000  0.300000  0.750000  0
  0.437500  0.312500  0.750000  0
  0.425000  0.325000  0.750000  0
  0.412500  0.337500  0.750000  0
  0.400000  0.350000  0.750000  0
  0.387500  0.362500  0.750000  0
  0.375000  0.375000  0.750000  0
  0.356250  0.356250  0.712500  0
  0.337500  0.337500  0.675000  0
  0.318750  0.318750  0.637500  0
  0.300000  0.300000  0.600000  0
  0.281250  0.281250  0.562500  0
  0.262500  0.262500  0.525000  0
  0.243750  0.243750  0.487500  0
  0.225000  0.225000  0.450000  0
  0.206250  0.206250  0.412500  0
  0.187500  0.187500  0.375000  0
  0.168750  0.168750  0.337500  0
  0.150000  0.150000  0.300000  0
  0.131250  0.131250  0.262500  0
  0.112500  0.112500  0.225000  0
  0.093750  0.093750  0.187500  0
  0.075000  0.075000  0.150000  0
  0.056250  0.056250  0.112500  0
  0.037500  0.037500  0.075000  0
  0.018750  0.018750  0.037500  0
  0.000000  0.000000  0.000000  0
其他文件从HSE06中拷贝
vasp计算
4.从EIGENVAL中取权重为0的高对称点做能带图

欢迎大家指正
关于KPOINTS的问题，这里我说明一下，从自洽计算的IBZKPT拷贝到KPOINTS中，然后我用的是生成高对称k点的脚本syml生成k点，把这些点的权重改为0，放在KPOINTS的最后，syml在侯柱锋的vasp指南里有介绍
另外说一下，我自己亲自验证次法可行，以上是我对计算的总结

[ Last edited by liliangfang on 2012-2-22 at 16:13 ]
回复此楼

» 收录本帖的淘帖专辑推荐

精华网帖收集 研究生资料集 VASP and MS 第一性原理
材料计算模拟实用技巧 VASP VASP 第一性原理相关文档
第一性原理计算经验 计算-vasp 第一性原理和电化学 分子动力学
@个人收集@ 杂化泛函方面的研究 VASP相关 纳米及第一性原理方向
VASP计算 Graphene 模拟，理论 vasp
vasp initio 第一性 VASP学习
VASP报错建模问题集 VASP-study 科研工具

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-09-19 16:02:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feifeikeke

木虫 (小有名气)

★★★★★ 五星级,优秀推荐

楼主，算能带时为何把参数ALGO 从D 改成了N 呢，有什么影响么？
赞一下(2人) 回复此楼
40楼2014-02-13 21:44:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 52 个回答

liliangfang

荣誉版主 (著名写手)
计算过程中汇出行下面提示，可以无视
|                                                                             |
|           W    W    AA    RRRRR   N    N  II  N    N   GGGG   !!!           |
|           W    W   A  A   R    R  NN   N  II  NN   N  G    G  !!!           |
|           W    W  A    A  R    R  N N  N  II  N N  N  G       !!!           |
|           W WW W  AAAAAA  RRRRR   N  N N  II  N  N N  G  GGG   !            |
|           WW  WW  A    A  R   R   N   NN  II  N   NN  G    G                |
|           W    W  A    A  R    R  N    N  II  N    N   GGGG   !!!           |
|                                                                             |
|      Your generating k-point grid is not commensurate to the symmetry       |
|      of the lattice.  This can cause   slow convergence with respect        |
|      to k-points for HF type calculations                                   |
|      suggested SOLUTIONS:                                                   |
|       ) if not already the case, use automatic k-point generation           |
|       ) shift your grid to Gamma (G) (e.g. required for hex or fcc lattice) |
|                                                                             |
赞一下(1人) 回复此楼
2楼2011-09-19 16:04:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

njuxyh

木虫 (著名写手)

★★★ 三星级,支持鼓励

0.50000000000000    0.00000000000000    0.50000000000000             1
0.500000  0.000000  0.500000  0
这两个点在KPOINTS 中，楼主难道你不报错吗？
赞一下 回复此楼
3楼2011-09-23 18:17:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lishen3314859

银虫 (初入文坛)

★ 一星级,一般

这个不是GW , 这个是HＳＥ，如果用ＧＷ方法是不可能计算出能带结构的。
赞一下 回复此楼
4楼2011-09-23 22:17:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 52 个回答
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭