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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » HSE（GW类似）计算能带总结
51/1返回列表
查看: 20858  |  回复: 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的回帖

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的回帖
查看全部 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的回帖

lishen3314859

银虫 (初入文坛)

★ 一星级,一般

这个不是GW , 这个是HＳＥ，如果用ＧＷ方法是不可能计算出能带结构的。
赞一下 回复此楼
4楼2011-09-23 22:17:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liliangfang

荣誉版主 (著名写手)
引用回帖:
3楼: Originally posted by njuxyh at 2011-09-23 18:17:07:
0.50000000000000    0.00000000000000    0.50000000000000             1
0.500000  0.000000  0.500000  0
这两个点在KPOINTS 中，楼主难道你不报错吗？

下一个权重是0 计算能带
赞一下(1人) 回复此楼
5楼2011-09-24 09:19:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 52 个回答
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿沪9，326生物学求相关专业调剂 +4 刘墨墨 2026-04-09 4/200 2026-04-10 12:07 by pengliang8036
[考研] 301求调剂 +4 XYPLR 2026-04-05 6/300 2026-04-10 09:17 by Delta2012
[考研] 08600生物与医药-327 +10 18755400796 2026-04-05 10/500 2026-04-10 08:14 by kangsm
[考研] 已调剂 +18 柴郡猫_ 2026-04-09 19/950 2026-04-09 22:10 by 柴郡猫_
[考研] 284求调剂 +7 让我上岸吧阿西 2026-04-09 7/350 2026-04-09 18:59 by haironglove
[考研] 求机械专硕297第二批调剂 +5 拾柒12。 2026-04-08 5/250 2026-04-09 16:43 by 允当适度
[考研] 085400电子信息类（川大控制工程）求调剂可跨专业 求老师联系 +3 626776879 2026-04-08 3/150 2026-04-09 16:05 by 猪会飞
[考研] 265求调剂 +4 风说她早忘了 2026-04-07 4/200 2026-04-09 13:59 by only周
[考研] 一志愿985初试354分生物调剂 +3 031001 2026-04-06 3/150 2026-04-09 00:30 by Evan_Liu
[考研] 二次调剂求老师收留 +3 笑笑袁 2026-04-08 3/150 2026-04-08 23:50 by 醉在风里
[考研] 生物学308分求调剂（一志愿华东师大）做过分子实验 +6 相信必会光芒万� 2026-04-07 7/350 2026-04-08 16:49 by tjzhao
[考研] 326分，一志愿沪9，求生物学调剂 +4 刘墨墨 2026-04-05 4/200 2026-04-08 06:22 by lijunpoly
[考研] 本科生物信息学，总分362 求07 08调剂 +6 q小倩1210 2026-04-06 6/300 2026-04-07 19:40 by macy2011
[考研] 306求调剂 +3 15287505595 2026-04-03 3/150 2026-04-07 18:08 by 蓝云思雨
[考研] 305分求调剂 +3 哈_哈_哈_哈_哈 2026-04-04 5/250 2026-04-07 14:49 by 哈_哈_哈_哈_哈
[考研] 286求调剂 +20 Faune 2026-04-06 20/1000 2026-04-07 11:33 by 诗与自由
[考研] 308求调剂 +3 终不似从前 2026-04-05 3/150 2026-04-05 20:07 by 啵啵啵0119
[考研] 348求调剂 +6 wukira 2026-04-04 6/300 2026-04-05 18:11 by 猪会飞
[考研] 一志愿沪9，求生物学调剂，326分 +6 刘墨墨 2026-04-04 6/300 2026-04-04 19:44 by 唐沐儿
[考研] 化工求调剂 +11 荔香芝士椰奶 2026-04-03 11/550 2026-04-03 22:06 by 啵啵啵0119
信息提示
关闭
请填处理意见
关闭