symfc 截断半径与超胞比较

截断半径决定力常数模型保留多远的原子相互作用；超胞决定周期模型能区分的空间范围。两者要一起检查。

下载完整Si实例包，解压到自己的工作目录。包内保留123个真实静态力条件、独立训练/验证/测试划分、14个拟合条件、原生导出文件与绘图脚本。

1. 同一超胞改变截断#

已有 labels/ 后：

bash fit.sh Si16-fc2-c3p0
bash fit.sh Si16-fc2-c4p2
bash fit.sh Si16-fc2-c5p0

这三组均为二阶模型，截断分别为3.0、4.2、5.0 Å。本例的 cutoff 字典显式列出各阶，避免二阶和三阶不小心使用不同的默认范围。

2. 改变超胞#

bash fit.sh Si54-fc2-c5p0
bash fit.sh Si128-fc2-c5p0

原胞2原子分别重复3和4倍，得到54、128原子。各体系的随机训练与测试结构独立生成，因此跨体系RMSE比较并非同一份测试集的逐构型误差差值。

3→4.2 Å通常改善本例测试误差；4.2→5 Å的收益较小，但Γ光学频率仍有变化。这说明需要同时看误差和目标性质，不能只观察某一条误差曲线。

小超胞中的长程相互作用可能通过周期映像混合。不同超胞相同截断下的独立参数数目也可能不同；这不是把短超胞参数直接复制到长超胞就能解决的问题。

3. 看声子而不只看力误差#

bash analyze.sh
head analysis-reference/fits.csv

CSV包含训练/验证/测试RMSE、截断、原子数和Γ光学平均频率。实际Phonopy色散另见：

4. 如何选参数#

先用验证集判断模型范围与阶数，再用未参与选择的测试集估计泛化误差；最后检查自己关心的声子、热导或自由能。更大的超胞、更高的阶数与更低的训练误差，都不能单独证明目标物性收敛。

本组只有三个有限超胞、几个截断和一个NEP势，属于实际操作与数值敏感性演示，不把图中的某个点称为收敛Si材料参数。