SeeK-path生成与检查高对称路径

本文使用鸢算科技集群实际运行的小算例。先下载输入、完整脚本与参考数据并解压，以下命令从解压目录开始。reference 保存原始参考结果，自己的运行放在 work，避免覆盖。

1 从结构生成高对称路径#

使用 pymatgen/2026.9.24 环境内的SeeK-path。读取晶胞后，软件会识别对称性，并基于标准化原胞给出路径。

module purge
module load pymatgen/2026.9.24
mkdir -p work/symmetry
cp scripts/symmetry_examples.py work/symmetry/
cd work/symmetry
python symmetry_examples.py seekpath-si

import seekpath
from ase.build import bulk
a = bulk('Si', 'diamond', a=5.43, cubic=True)
cell = (a.cell.array, a.get_scaled_positions(), a.numbers)
path = seekpath.get_explicit_k_path(cell, reference_distance=0.05)
print(path['path'])
print(path['point_coords'])
print(path['primitive_lattice'])

reference_distance=0.05控制倒空间采样间隔，单位Å⁻¹；它不是“每段20点”。实际点数随路径长度确定。

2 本例结果#

分别输入8原子Si常规胞及16原子超胞，均识别为227号空间群，输出标准化路径的6段、98个显式采样点。文件 si-path.json 保存端点和路径顺序，si-kpoints.csv 保存沿程距离与绝对k坐标。

cat seekpath-si/si-path.json
head seekpath-si/si-kpoints.csv

3 最重要的坐标基底检查#

默认输出对应标准化原胞的倒格矢。不能直接把输出的分数k坐标贴到任意原始超胞的KPOINTS里。 先确认计算所用POSCAR与输出原胞一致，或采用SeeK-path针对原始晶胞的接口并理解超胞能带折叠。本例只生成路径，没有运行电子能带计算。

查看和复现配图#

下载包 figures 中是高清图片，可点击正文图片放大。reference 中保存数据与结果表；重新绘图可执行：

module purge
module load ase/3.29.0
mkdir -p work/replot
cp -r reference work/replot/
cp plot_reference.py work/replot/
cd work/replot
python plot_reference.py

不同教程共用此资料包；同一个脚本不需要重复下载。重新运行算例时使用新的工作目录。

方法与接口参考：SeeK-path官方接口。