ZnS Gamma LO–TO与非解析修正
先在新目录解压，bash preview.sh 重画参考图；bash run.sh 使用原生DFPT导出的完整自有数值运行真实Phonopy Wang Gamma NAC，再与独立Born/介电rank-one矩阵、LST关系核对。只需1核后处理，不重新计算DFT。复跑换新目录。
reference/包含完整2原子DFPT力常数/质量加权Hessian、Born张量、电子与离子介电张量、结构与实际PAW质量。没有POTCAR/波函数/完整XML/OUTCAR或厂商程序。
要计算自己的ZnS，先准备合法POTCAR，在inputs/ZnS-DFPT中sbatch run.slurm；结束后回包根目录，module load phonopy/4.8.3，python read_vasp_response.py inputs/ZnS-DFPT --output my-response，再python analyze_nac.py --data my-response --output my-results及python plot.py --data my-results --output my-plots。
本模型仅2原子原胞Gamma DFPT，不包含完整有限q色散。三个方向为笛卡尔100/110/111，程序转换为实际原胞倒格矢坐标。Phonopy4.8.3，VASP6.6.1 CPU oneAPI2024.2/IntelMPI2021.13；原生输入550eV，EDIFF1e-6，Gamma7^3，固定a5.45A。现有V10016分区模板调用CPU程序，GPU仅分配不用于该CPU计算。
TO259.59139cm-1，Wang NAC LO326.50019cm-1；LST响应326.41684cm-1，约0.0255%差别保留。未施加ASR/修正Born；原声学近零-0.498/-0.498/-0.396cm-1保留，不能凭此断言材料结构不稳定。图中箭头为共同缩放后的模式方向，不是真实热振幅。
