鸢算科技：Cu111 选能窗密度与py4vasp STM实操
bash preview.sh重画参考图；bash run.sh依次准备最小HDF读取载体、提取数值、运行原生图形工具、绘制三张PNG。results/data为数值，tools-output为8个原生HTML，plots为PNG，figures为随包参考图。使用新解压目录，输出目录受保护，不重新运行DFT。
实际LPARD使用VASP6.6.1 GPU/HDF5、NVHPC25.3/CUDA12.8，P100与V10016G单GPU，复用自己的Cu111五层电子态；550eV、Gamma9x9x1、33能带。两能窗-0.5到0和0到+0.5eV都相对EF。输入完整INCAR/POSCAR/KPOINTS在inputs，匹配的PAW/WAVECAR/CHGCAR需自行准备，本包不含。
LPARD用单独目录并复制输入，不把原始SCF的CHGCAR硬链接或软链接过来。本例原生程序会重写该任务目录中的CHGCAR，即使LCHARG=false；原SCF文件必须保留，绘图读取partialdensity数据。
两套40x40x420选定数值保留全分辨率，无插值或手工归一化。carriers是自有NPZ和字段配方，运行时构造明确标记的最小HDF读取载体，不是完整原始VASP输出。
后处理py4vasp0.11.3/Python3.11，静态图pymatgen/2026.9.24环境的NumPy2.3.5/Matplotlib3.11.2；已分配资源内1CPU即可。各目录受保护，重做请换新解压目录。
原生HDF部分电荷单位通过选定态k权重独立核对：mean(raw)=50/81和8/9电子；物理密度rho=raw/V，不再除网格点数。它们是能窗态密度的权重，不应当作实际额外载流子数。
此py4vasp版本to_stm的恒高显示值为gaussian(rho)/Ngrid*enhancement。包内保留原生HTML，同时静态图用正确rho单位，显示倍率不是实际实验电流。
恒流API接受带nA标签的参数，但这里明确作为软件密度阈值映射，未做针尖/势垒/电流标定。阈值取2A最近平面平滑密度的中位数，传给API的值=threshold/Ngrid*1e9。比较两能窗时阈值分别选择，不把色深直接比较为实验电流。
此版恒流graph.data是由网格索引缩放的相对坐标；转物理相对高度用native_data*c/Nz*interpolation_factor，不直接给其原值贴angstrom标签。脚本与上方真空中独立线性等密度面求根核对，默认10和加密50均满足两子网格步长误差界。
输入tip_height=2/3A对应实际最近网格高度2.026508/2.969177A，需与请求值区分。当前常高数据采用默认4网格点Gaussian平滑，sigma单位是网格点，不是A。
恒流图以最小高度为零，只显示起伏；不作为绝对针尖高度。插值因子50改善z采样，不会创造新的真实电子结构信息。静态恒高图共用色标，恒流两面板分别用色标，线剖面保留原生台阶，无平滑伪造。
本例为非磁Cu的简化LDOS-STM流程，不含真实针尖态、隧穿矩阵元或全实验条件，不将结果声称为已收敛实验STM电流。更换体系需收敛结构、空带、k点、真空与能窗。
