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

登录注册

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

24小时热门版块排行榜

返回列表

HZ_Tian

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 49.3
散金: 4
帖子: 18
在线: 12.5小时
虫号: 6383842
注册: 2017-04-24
专业: 理论和计算化学

[交流] 一种简洁快速的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

回复此楼

» 猜你喜欢

1楼 2017-07-26 10:47:00

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

surelv

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 18.2
散金: 5
帖子: 53
在线: 82.8小时
虫号: 6035178
注册: 2017-03-17
专业: 理论和计算化学

我也用的pwmat