Si16 Fourthorder + FourPhonon 真实受力实例
先看参考图：bash preview.sh
重建力常数：bash prepare.sh → bash reconstruct.sh
运行15个CPU输运条件：sbatch cpu.sh（8核AVX512 CPU资源），随后查看results/data、plots/figures。
声子辅助工具：bash phonons.sh，读取同一份真实二阶常数，由Phonopy4.8.3画能带与DOS。
本包使用152组逻辑受力，来源是119个独立真实VASP6.6.1静态DFT、33个相同构型复用；不重算电子结构。
H=0.01/0.02/0.04 Angstrom，2x2x2 16原子，最近邻截断；位移工具内部H单位是nm，程序在work/programs私人副本修改H，共享程序不改。
三组各48逻辑位移；三阶8组。真实二阶由原Si16DFT数据复用，符号与原子顺序均已转换。
reconstruct.sh生成的XML只包含实际数值力，是给原生解析器的最小数据载体，不是完整VASP电子输出。
QE接口例是同一份真实VASP受力按晶体对称性/原子排列及Ry/bohr单位转换后的格式互检；没有运行QE电子计算，不能用它证明两套赝势一致。
H010的平移和规则偏差8.51e-6 eV/A^4超过预设1e-6，保留作位移噪声对照，不用于输运。H040有36独立构型力检查：RMSE1.8555e-4 eV/A。
3^3/4^3/5^3热导未收敛，只是终端操作例，非Si材料收敛预测。CPU为定量参考；GPU逐支散射率对照不通过预设0.5%，没有把热导接近等同于逐模一致。
CPU支持RTA、完整迭代、同位素、抽样。GPU当前单卡单MPI，只用确定性RTA，不支持这里的迭代/抽样例。
num_sample_process=10000只是一组抽样示例；原生用random_seed()而无固定种子入口，不能宣称独立复跑完全相同或统计收敛。
包内reference是实际计算参考，work是新重构/新输运结果，results是独立提取；不会写回reference。
安装：Fourthorder0.1-fix20260618 Python3.11/NumPy1.26.4/SciPy1.13.1/GCC13.3；Thirdorder1.1.4；FourPhonon1.3-fix20260924(固定5b814b3)，CPU oneAPI2024.2/IntelMPI2021.13/MKL顺序，GPU NVHPC25.3/CUDA12.8，P100和V10016。
仅自有输入、受力、IFC和数值工具；不含PAW、波函数、完整电子输出、第三方源码/二进制或私人执行凭据。
