Wannier90 位点对称适配：GaAs的s与p轨道

下载完整输入、参考数据和绘图程序，保存为example.zip：

python3 -m zipfile -e example.zip .
cd wannier-utilities-examples
bash results.sh
bash plot.sh

环境使用wannier90/4.0.3-serial，Intel ifx2024.2、AVX2、MKL顺序库；MPI构建使用Intel MPI2021.13。本文小例子用单线程即可。job.slurm申请CPU分区资源，再由run.py在新的work/算例名下运行。查看参考结果不需要重新计算；重复计算时换用新的解压目录。

1. 需要什么文件#

本例使用官方4.0.3 GaAs As-s/p固定DFT矩阵，在集群实际执行对称适配Wannier化。原始DFT不是这次重新计算的；本例展示如何消费已有的dmn和重叠矩阵。

ls inputs/GaAs-site-symmetry
cat inputs/GaAs-site-symmetry/GaAs.win

除通常的win/amn/mmn/eig外，需要与同一k网格、同一能带顺序对应的GaAs.dmn。dmn包含对称表示信息，不能从别的材料随意复制。

关键输入为：

num_bands = 4
num_wann = 4
site_symmetry = true
symmetrize_eps = 1e-9
write_hr = true
write_u_matrices = true

完整输入中的As-s/p投影、晶胞和k点以随包文件为准。site_symmetry不同于普通改变投影标签；它还要求相应对称矩阵。

2. 提交与读结果#

sbatch job.slurm GaAs-site-symmetry
grep -A 9 'Final State' reference-results/GaAs-site-symmetry/GaAs.wout

实际得到一个s型函数和三个p型函数：s展宽1.63710680Ų，三个p展宽均为2.83312865Ų，总展宽10.13649274Ų。四个原始中心都在As位点附近，周期等价的笛卡尔表示约为(-1.412903,1.412903,1.412903)Å，与As位置的差小于2×10⁻⁷Å。

中心在原子位点、三个p展宽相同，适合用作本例对称适配的直观检查；是否可用于您的紧束缚模型，仍要对照原始能带、能窗与目标轨道。不要用最小展宽一个标准比较对称适配和不受对称约束的MLWF。

原理和输入要求见官方对称适配GaAs教程。