24小时热门版块排行榜

返回列表

【奖励】本帖被评价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人回复
磺酰氟产物，毕不了业了！已经有4人回复
论文终于录用啦！满足毕业条件了已经有16人回复
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球已经有8人回复

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

能带计算的问题已经有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的回帖

gaoenlai110

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 428.6
帖子: 77
在线: 81.1小时
虫号: 2303907

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

http://muchong.com/bbs/viewthread.php?tid=3623104
根据这个官网指导，第二步似乎没必要啊。另外，HSE计算时需要拷贝上一步的WAVECAR吗？我没拷贝也计算出结果了

赞一下

回复此楼

45楼2016-01-13 14:45:02

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

gaoenlai110

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 428.6
帖子: 77
在线: 81.1小时
虫号: 2303907

引用回帖:

46楼: Originally posted by liliangfang at 2016-01-13 19:52:12
我现在已经不做这方面了，改行做土木了。这是我好几年之前做的了，当时这个方法还没完善。我当时也是根据各方面的资料验证和总结的方法，会有一些问题出现。当然现在方法成熟了你可以自己总结和改善。希望对你有帮 ...

谢谢。新行业加油！

发自小木虫Android客户端

回复此楼

47楼2016-01-13 23:20:03

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

相关版块跳转我要订阅楼主 liliangfang 的主题更新

返回列表

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

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