LOBSTER 自定义STO基组输入

下载完整输入、原生输出与三张图的绘图程序，保存为example.zip：

python3 -m zipfile -e example.zip .
cd lobster-bands-custom-examples
bash results.sh
bash plot.sh

使用LOBSTER6.0.0官方g++9.3.0静态OpenMP程序，单线程。VASP6.6.1 GPU环境为NVHPC25.3、CUDA12.8、HPC-X2.22.1、MKL2023.2。POTCAR、WAVECAR和完整DFT输出不分发；完整安全输入与原生LOBSTER数值输出在资料包中，重算使用您自己生成的同源波函数。

1. 本例到底改变了什么#

我们从此前Si分析实际输出的基函数参数构造自定义STO文件，保留选用的3s/3p径向参数和六位小数系数，只演示自定义基组接口。没有重新拟合材料能量，也不把它称为性能或精度升级。

cat inputs/Si-custom-fullshell/lobsterin
head -n 16 inputs/Si-custom-fullshell/Si-fullshell-derived.sto

输入选择：

basisSet custom
basisFunctions Si 3s 3p
customSTOforAtom Si Si-fullshell-derived.sto
writeBasisFunctions

本版本也实际接受手册示例中的useBasisSet custom写法；两套输入都在资料包中，不通过忽略未知关键词的退出码来猜测是否真的选中了自定义基组。

2. 看懂STO文件头与列数#

本例采用完整壳层编号布局：

si, Z=14
FirstOrbitals: 1 2 0 0
Orbitals: 3 2 0 0
Primitives: 12 11 0 0

后面S行包含q、指数alpha和三个s壳层的系数列；P行包含q、alpha和两个p壳层的系数列。alpha单位为Bohr⁻¹，q和系数按原生格式无量纲。不能删去未使用的列后仍保留原文件头，也不能只改元素名称就认为文件适用于另一种元素。

本例实际只请求3s/3p；未选用的1s/2s/2p槽用明确占位参数补全文件布局，不用于物理核心态计算，也不是原子的全电子基组拟合。若要使用核心轨道、其他n/l或不同电子结构，需准备真正对应的完整基组数据，不能套用这些占位参数。

图展示输入参数，未用未说明的归一化公式拼出一条所谓真实原子径向波函数。

3. 提交并核对实际采用的基组#

sbatch job.slurm Si-custom-fullshell --source "$(pwd)/my-si-vasp"
sbatch job.slurm Si-custom-alias-fullshell --source "$(pwd)/my-si-vasp"
grep -A 5 'setting up local basis functions' reference-results/Si-custom-fullshell/lobsterout

原生输出明确标记：

Si (custom) 3s 3p_y 3p_z 3p_x

两个原子共8个局域函数；多余PAW带从第9条起按本例常规策略忽略。绝对电荷spilling约1.46%，仍需连同真实采用的基组、电子数和插值质量判断；不因为输入文件是自定义的就自动认为更准确。

换了基组就重新投影，不加载此前基组的projectionData.lobster。基组质量的一般判断见基组与投影质量。