版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

>有奖起名 (21735)
>公派出国 (10138)
>找工作 (9709)
>博后之家 (9490)
>考研 (8901)
>海外博后 (7056)
>留学生活 (6240)
>考博 (5913)
>留学DIY (5461)
>教师之家 (5457)
>职场人生 (4612)
>访问学者 (4252)
>论坛更新日志 (3590)
>文学芳草园 (3510)
>虫友互识 (3468)
>育儿交流 (3106)

返回列表

【奖励】本帖被评价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 ]

回复此楼

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

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

能带计算的问题已经有5人回复
到底什么是HSE06计算已经有3人回复
HSE能带计算出错已经有6人回复
关于能带自洽计算是否完成的判断已经有4人回复
HSE06计算的势文件已经有4人回复
hse计算中ICHARG=0，不收敛已经有10人回复
vasp画能带图的疑问已经有12人回复
HSE能带计算中ISYM参数的问题已经有7人回复
HSE06计算突然停止？已经有6人回复
百草枯原药有没有，95%含量的原粉？！已经有4人回复
【重点讨论】计算能带时如何确定高对称点和k-points path? 已经有40人回复
【原创】ICHARG=11和+U的计算已经有9人回复

1楼 2011-09-19 16:02:14

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

feifeikeke

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 2959.9
帖子: 159
在线: 194.2小时
虫号: 964759

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

楼主，算能带时为何把参数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

木虫 (著名写手)

应助: 1 (幼儿园)
金币: 1802
帖子: 1320
在线: 771.5小时
虫号: 455673

★★★ 三星级,支持鼓励

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

银虫 (初入文坛)

应助: 0 (幼儿园)
金币: 136.5
帖子: 9
在线: 84.8小时
虫号: 503613

★ 一星级,一般

这个不是GW , 这个是HＳＥ，如果用ＧＷ方法是不可能计算出能带结构的。

赞一下

回复此楼

4楼2011-09-23 22:17:04

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 52 个回答

☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级

普通表情龙兔虎猫高级回复 (可上传附件)

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[基金申请] 时间戳变了，能看出什么问题？ +4	基诺咪客 2026-08-17	5/250	2026-08-18 08:03 by 基诺咪客
[基金申请] 93BebMhtakh前后11位开头都是大写 +4	且听虎啸 2026-08-17	5/250	2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6	angus9576 2026-08-17	11/550	2026-08-17 23:57 by angus9576
[基金申请] filecode=后面第一个是大写字母 +8	wangze12014 2026-08-14	10/500	2026-08-17 17:05 by xter9665
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜？ +12	Ldrop2023 2026-08-13	15/750	2026-08-17 15:02 by 小豌豆_发芽
[基金申请] 今天系统多次维护，明天很可能放榜！ +8	zju2000 2026-08-16	9/450	2026-08-17 12:20 by lmz0216
[基金申请] 时间戳又变了8-15 +13	archvillain 2026-08-15	25/1250	2026-08-16 20:13 by zhaosm1982
[基金申请] 2027广东省杰青 +3	奶牛小黑 2026-08-15	6/300	2026-08-16 20:07 by 奶牛小黑
[基金申请] 快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。 +3	Tide man 2026-08-14	3/150	2026-08-16 17:47 by jurkat.1640
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7	启萌科技 2026-08-12	26/1300	2026-08-16 12:35 by 启萌科技
[精细化工] 招聘金属平磨液，抛光液研发工程师 +3	小天0311 2026-08-14	3/150	2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友，我要去昆明玩几天，回来见。 +7	Tide man 2026-08-14	8/400	2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果？ +4	yuleib84 2026-08-11	4/200	2026-08-14 23:05 by lfy8008
[硕博家园] 读博的好处 +4	lnee 2026-08-11	4/200	2026-08-14 10:20 by ahsoarli
[基金申请] 重要来源：本周末出结果 +10	瞬息宇宙 2026-08-12	10/500	2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7	且听虎啸 2026-08-12	9/450	2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了，巨变 +3	WH3796 2026-08-12	4/200	2026-08-13 14:13 by 小木虫6752397
[基金申请] 结合人工智能，周易传统文化，filecode打分制来了，3分以上希望很大。 +3	Tide man 2026-08-12	4/200	2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见，大家看看几个A，几个B? +11	Tide man 2026-08-11	11/550	2026-08-13 07:35 by 撸猫猫
[基金申请] 什么时候出结果，有咨询渠道？？？ +3	Tide man 2026-08-11	3/150	2026-08-11 17:54 by kudofaye