鸢算科技：缺陷结合能、费米面截面与液体动力学
各目录准备合法PBE PAW64 POTCAR，元素顺序与POSCAR一致。Si/Cu使用vasp_std，Na使用vasp_gam，run.sh经Slurm申请P100计算资源。
1. 缺陷：分别计算si-defect-bulk、si-vacancy、si-divacancy-near；固定共同晶胞，优化原子。完成后运行python extract_defect_binding.py。
2. 费米面：分别计算cu-primitive-mesh13与21，然后python extract_fermi.py。需要ASE、NumPy、spglib；module load ase/3.29.0提供本例环境。不使用磁性/SOC体系套用时间反演重建。
3. 液体：先na-melt，完成后在包根目录python prepare_md.py，选择1/2准备两条热化或3准备高温生产。热化完成后再选择4/5准备相应生产。每次进入目标目录sbatch run.sh。脚本沿用前一步CONTCAR的结构和所含速度；温度以生产轨迹实测为准。三条生产均完成1000步后python extract_liquid.py，读取XDATCAR、OSZICAR、OUTCAR。
附件POSCAR为本次实际初始状态；重新运行完整MD时应通过prepare_md.py接续自己计算的末态。参考RESULTS.json和figures不会自动随新任务更新。
Na仅16原子、固定盒和2ps生产，用于RDF/MSD/VACF工作流，不给出收敛的扩散活化能、黏度或热导。NVE不是结构优化，瞬时原子力不要求小于0.02eV/Å。
不含POTCAR、WAVECAR、CHGCAR、原始运行日志或用户私有数据。
