qe计算电子结构

结合案例图分析，第一张案例图横轴是高对称路径Γ→X→M→Γ，纵轴是能量。红色虚线标出费米能级E_F≈2.0 eV。能带在费米面附近有多条带交叉，说明体系具有金属性。Γ-X方向色散较陡，X-M方向较平，暗示电子在Γ-X方向迁移率更高。图中标注”scf → bands.x”，说明这是先自洽得到电荷密度，再用bands.x做非自洽能带计算的标准流程。Bands.x读取的是pw.x输出的电荷密度文件，不会重新优化波函数，因此能带更平滑。

第二张案例图横轴是波函数截断能ecutwfc（Ry），纵轴是总能量（Ry）。可以看到截断能从15 Ry增加到约50 Ry时，总能量快速下降；50 Ry之后曲线趋于平坦，图中标注”ΔE<1 meV/atom”对应45 Ry以后。这意味着对于该体系，45 Ry的ecutwfc已经足够收敛，再增加截断能只会浪费机时而不会显著改变结果。

从能带图和截断能收敛测试的组合来看，这是一个用Quantum ESPRESSO（QE）计算电子结构的项目——客户想看看某个立方晶系金属氧化物的能带特征，以及截断能收敛情况。能带图回答了体系的电子结构性质（金属性/半导体性），收敛测试回答了”用多大截断能才够”这个最基础的问题。

项目计算参数上，我用QE的pw.x做自洽计算，pseudo选择GBRV赝势库。K点网格取8×8×8 Monkhorst-Pack，因为体系是立方而且含金属，需要较高密度。smearing用Methfessel-Paxton，展宽0.02 Ry。bands.x计算时沿高对称路径取了100个点，平滑效果很好。我最终选用了50 Ry作为生产计算的截断能，留了一点安全余量。

QE与VASP相比有几个不同点。第一，QE的输入文件是文本形式，需要手动写&SYSTEM、&ELECTRONS等namelist，对初学着门槛稍高。第二，ecutrho（电荷密度截断能）默认是ecutwfc的4倍，对于含氧体系我建议保持这个比例或更高，因为氧的2p轨道局域。第三，QE的赝势库很多，我用的是ONCV或GBRV，换pseudopotential时要重新做收敛测试，不能直接把VASP的截断能搬过来。

我踩过的两个坑值得一提。一个是k点路径写法：QE的k点路径用crystal_b格式，如果晶胞不是标准原胞，高对称点坐标会错，导致能带形状与文献不符。我先用xcrysden确认了倒空间路径，再写bands.in。第二个是非自洽网格：bands.x不需要自洽网格和能带网格一致，但要保证电荷密度文件存在且K点路径正确。如果忘记设置restart_mode=’from_scratch’或prefix一致，会报找不到文件的错误。

这张图对客户的意义是：50 Ry截断能+8×8×8 K点是可靠且经济的计算设置，能带显示金属性且费米面处态密度不低，适合做输运或催化研究。关于QE计算电子结构的更多技巧，我整理在qe计算电子结构栏目里。