鸢算科技：py4vasp电子—声子带隙与恒定弛豫时间输运实操
先看真实参考图：bash preview.sh。重新提取、原生工具、绘图：bash run.sh。
结果在results/data、tools-output和plots/figures。重复操作请使用新解压目录，输出目录受保护。
py4vasp0.11.3/Python3.11；静态图NumPy2.3.5/Matplotlib3.11.2。已分配资源内1CPU即可后处理，不新运行VASP。
数据来自真实VASP6.6.1GPU/HDF5、NVHPC25.3/CUDA12.8单GPU计算。EP使用已有真实声子势导数，新的第一份带隙计算保存电子态，随后任务复用收敛WAVECAR/CHGCAR，ICHARG11与ALGO Normal做固定电荷电子态更新及新的目标性质。
输入完整保存INCAR/POSCAR/KPOINTS/KPOINTS_ELPH。若要重新做电子结构计算，必须自己准备匹配PAW及声子势导数、需要时的前置波函数/电荷；这些不在本包中。
carriers是自有真实数值NPZ与字段配方，运行时创建显式标记的最小HDF读取载体，不是原始完整VASP输出。选定数组完整保留，未抽样或插值。
bandgap与transport默认read只列累加器数量；用[0]、[1]选择实例，再to_dict()/to_plotly()。Python索引0对应OUTCAR累加器1。
带隙：Si原胞2、550eV、SCF4^3/NB24，插值6^3/8^3、ELPH_NBANDS16，delta.05/.10eV，T0/300/600K。fundamental是实际采样网格的间接最低间隔，direct是同k最低直接间隔；不是将Gamma选中两态之差当全局带隙。数组renorm表示修正后的带隙，修正量应减bare。
源分析由实际稀疏bks_idx、band_start、KS能量与Fan/DW逐态重建全部温度带隙。Fan真实部与DW显式按态/温度相加，不能套用不匹配的广播形状。6/8网格和两种展宽仅作方法对照，未宣称q/k/空带收敛。
输运：真实原生CRTA、6^3插值、101能量节点、T300/600K、目标掺杂-1e19/+1e19cm^-3。tau10/20fs是假定值，不是本次由声子散射算出的寿命。负掺杂为空穴、正掺杂为电子；不能用这些值声称实际声子限制迁移率。
原生Seebeck单位V/K、Peltier单位V，OUTCAR以microV/K和microV打印；电导率S/m、迁移率cm^2/(V.s)、电子热导W/(m.K)。Pi=T*S转置，tau加倍时sigma/mobility/kappa翻倍、S/Pi不变。
输入目标掺杂和积分得到的实际多数载流子数要区分；n字段约8是电子数/晶胞，不是cm^-3。用实际ne/nh和晶胞体积换密度。小能量网格有离散积分误差。
为与该VASP版本逐位核对，脚本保留其EVTOJ常数及原代码单精度单位字面量提升规则；使用现代SI常数的约1e-7相对差远小于本小网格模型的不确定度。
