跳转至
课程笔记 · 08

8. 参数拟合与可迁移性

本课位置: 第 8/10 课。先完成前一课的推导;使用下方练习检查理解。

1. 目标与假设

紧束缚参数代表一个训练适用域。需要明确参考包含哪些结构、应变、成分、磁态和能量窗口。只拟合一个平衡结构能带的模型,可能在压缩或断键时失败。可迁移性要通过目标范围内的独立测试证明,而不能由训练残差小推断。

目标函数应使用与研究目的相关的观测量。能带可以重点加权费米面附近或规定窗口;若目标是力和总能,还需完整且一致的能量模型,不能只依赖单电子谱。对称相关参数应共同约束。正则化可限制不稳定或过大的参数,但它有需要选择和记录的强度。若同时拟合能量和力,各项单位决定权重含义,不能直接把不同量纲平方相加。

仅拟合能量可能无法唯一确定参数。共同现场能平移使所有能带平移;轨道规范变化能改变矩阵元而保持本征值。能带交叉还使按序号配对变得危险,某条排序能带可能交换轨道性质。必要时应使用对称标签、子空间或波函数重叠进行对应。可辨识性分析帮助区分真正受到参考约束的参数组合与任意选择。

验证集应按结构划分,而不是把同一个结构相邻波矢点随机拆分。高度相似样本共享信息,容易产生过度乐观测试误差。应保留独立应变或化学环境,同时报告平均误差与最差情况,并检查物理约束。键长或配位超出训练范围时要明确属于外推。即使离散能量点拟合很好,径向函数和截断不平滑仍会使力不可靠。

2. 逐步推导

逐行阅读等式及其假设。DFT 方程默认采用原子单位,另有说明除外;紧束缚参数保持明确能量和长度单位。匕首表示共轭转置,坐标上的撇号表示独立积分变量。

\[ \begin{aligned} \mathcal L(\boldsymbol\theta)&=\sum_{s,n,\mathbf k}w_{sn\mathbf k}[E^{TB}_{sn\mathbf k}(\boldsymbol\theta)-E^{ref}_{sn\mathbf k}]^2+\lambda\|\boldsymbol\theta\|^2,\\ t(R)&=t_0e^{-\alpha(R-R_0)},\quad\frac{dt}{dR}=-\alpha t(R),\\ \mathrm{RMSE}&=\sqrt{\frac1M\sum_i(y_i-\hat y_i)^2},\\ H_{ij}&=H_{ji}^*. \end{aligned} \]

2.1. 总能拟合不仅需要本征值

模型复现占据能带仍可能因其他能量项缺失而预测错误键长。力模型应分清电子、排斥、电荷和相互作用贡献。参考混合不同自旋或泛函时,要判断是目标差异还是不兼容标签。交叉验证应保留完整环境并报告适用域。不同参数向量若给出相同损失,就不能对各参数作唯一物理解读;域内预测仍可能稳定,但外推尤其不确定。

3. 计算示例

已知相同现场能且无重叠的双站点能级分裂为二倍跳跃绝对值。分裂四电子伏特确定绝对值为二,却不能在未规定轨道相位时确定符号。若两个能量都知道,平均值确定现场能,半差确定跳跃绝对值。

4. 习题与解释性答案

习题. 推导指数跳跃在小键长变化下的响应。

解释性答案. 小变化下指数约为 1−αδR,所以跳跃变化约为 −αt0δR。α 的单位是长度倒数。忽略导数的拟合可能复现能量,却给出不可靠电子晶格耦合。

进一步检查. 写出算例各参数的单位和适用范围。改变一个假设,指出推导哪一步必须重做。正确解释应说明受到影响的约束、算符或边界条件,而不只是说结果会变化。

5. 误区与局限

平均误差小不证明近简并相的排序正确。应将误差尺度与具体准备解释的能量差比较。

最近邻链及余弦色散

图为原创教学示意,不是数值材料模拟输出。

6. 联系与来源

相关:密度泛函基础 · 多体蒙特卡罗方法

DFTB+

解释、推导和假设算例为原创教学综合。所引文献建立相关理论,不是这些教学数值的来源。

 依照本课假设绘制的原创解析教学图,不是模拟结果。

依照本课假设绘制的原创解析教学图,不是模拟结果。

7. 关联理论与实践

量子力学 · 密度泛函理论 · ABACUS · DFTB+


上一课 · 课程目录 · 下一课