[交流] 一种简洁快速的HSE能带计算方法 已有1人参与

最近研究了下用PWmat计算HSE能带，发现速度还是很快的，而且输入参数设置也是很简单的，主要是分为三个步骤： 1.对体系进行结构优化，设置参数也是很简单，主要是设置参数输入文件etot.input中的参数 “JOB=relax”和“RELAX_DETAIL”，如果优化也想用杂化泛函HSE，将参数“XCFUNCTIONAL”设置为"hse"，如果不用杂化泛函HSE，可以不用设置该参数，直接用默认值（LDA）,具体的设置如下图 http://muchongimg.xmcimg.com/data/bcs/2017/0726/w343h6383842_1501034304_821.png 2.对优化后的结构进行自洽计算，将参数"JOB"设置为“scf”，并且一定要把参数“XCFUNCTIONAL”设置为"hse"，即进行HSE自洽，具体参数设置如下图 http://muchongimg.xmcimg.com/data/bcs/2017/0726/w414h6383842_1501035597_746.png 3.进行非自洽计算，将参数"JOB"设置为“nonscf”，并且需要读入上一步自洽得到的势函数文件即OUT.VR，具体就是将OUT.VR拷贝到IN.VR，并在参数输入文件etot.input中加入“IN.VR = T”，另外还需要在K点文件IN.KPT中设置高对称点，在etot.input中加入“IN.KPT = T”，具体参数设置如下图 http://muchongimg.xmcimg.com/data/bcs/2017/0726/w351h6383842_1501036330_798.png 最后用PWmat自带的数据处理工具plot_band_structure.x进行数据处理，得到能带数据文件bandstructure_1.txt用于画图。 具体测试了一个GaAs原胞的例子，对于HSE结构优化、HSE自洽和HSE非自洽三步总共用时4个小时，这个速度还是很快的，具体得到的能带结构图如下图 http://muchongimg.xmcimg.com/data/bcs/2017/0726/w149h6383842_1501037135_605.png relax.png scf.png nonscf.png band.png

回复此楼