10. 微扰理论与响应
本课位置： 第 10/10 课。先完成前一课的推导；使用下方练习检查理解。
先修知识： 本征态与矩阵元
1. 明确小参数与能隙
写 \(H=H_0+\lambda V\)，其中 \(H_0|n\rangle=E_n^{(0)}|n\rangle\)。假设目标为孤立非简并离散能级，扰动矩阵元相对能隙足够小。同时展开能量与波函数，采用中间归一化 \(\langle n|\psi_n\rangle=1\)。比较一阶项得
投影到原态得到 \(E_n^{(1)}=V_{nn}\)；投影到其他态得到 \(\langle m|\psi_n^{(1)}\rangle=V_{mn}/(E_n^{(0)}-E_m^{(0)})\)。再把二阶方程投影回原态，有
小分母直接暴露近似风险。基态的耦合态在上方，故二阶修正非正；激发态可能有正负两类贡献。这是幂级数系数，并非每一阶都具有变分上界。
2. 用可解模型检验展开
取 \(H=\begin{pmatrix}0&v\\v&\Delta\end{pmatrix}\)，\(\Delta>0\)。精确低根为 \((\Delta-\sqrt{\Delta^2+4v^2})/2\)。在 \(|v|/\Delta\ll1\) 时展开根号，得到
比值为零点一时二阶结果接近精确值；能隙趋零时，二阶表达式发散，而精确能量保持有限并趋于 \(-|v|\)。此时应在耦合的低能子空间内直接对角化，不能继续把小分母代入非简并公式。这也是理解多参考电子相关的起点。
3. 简并微扰的正确步骤
若若干态共享同一未扰动能量，先取投影 \(P\)，对角化 \(PVP\)。所得本征向量是一级分裂后适当的组合，然后才处理与补空间的耦合。原简并空间中的任意基底，其各对角矩阵元一般不代表独立的物理能级位移。对称性常能帮助缩小矩阵并检验禁止耦合。
4. 完整例题：电场响应
沿 \(z\) 的电场与偶极相互作用为 \(V=-F\hat\mu_z\)。无永久偶极、宇称为偶的基态一阶位移为零；二阶能量写成 \(-\alpha F^2/2\)，其中
两能级教学模型的间隔为 \(\Delta\)、跃迁偶极为 \(d\)，便有 \(\alpha=2d^2/\Delta\)。能隙减半使弱场极化率加倍，同时也使混合参数 \(Fd/\Delta\) 更快变大，线性响应适用的场强范围缩小。很大的响应可能意味着接近微扰失效边界。
5. 练习、解释与局限
若扰动在本征基中完全对角，二阶混合修正为何为零？ 因为所有非对角矩阵元为零。能级仍可有一阶位移，但本征向量不变。
若存在永久偶极会怎样？ 能量中先出现 \(-F\langle\mu_z\rangle\) 线性项，应将永久极化与诱导极化分开。
连续谱要用积分替代部分求和，共振和强场需要不同处理。实际响应计算中，基组与电子相关会同时影响能隙及跃迁矩阵元；总能量收敛不能保证极化率准确。应比较不同场强的有限差分结果，并检查线性或二次拟合是否稳定。
解析教学示意图，不是模拟数据。
进一步的概念检验
用恒定力扰动谐振子，基态经位置算符只耦合到第一激发态。代入其矩阵元及间隔得到 \(\Delta E=-F^2/(2m\omega^2)\)。直接对 \(m\omega^2x^2/2-Fx\) 配方也得到同一常数位移，并把平衡位置移到 \(F/(m\omega^2)\)。两条独立代数路径检验二阶修正的符号与因子。该模型特殊地由二阶能量给出精确答案，不代表一般微扰都在二阶终止。
6. 参考资料与学习联系
本课推导及模型算术为原创教学构造；这里没有报告实际执行的分子模拟。
7. 关联理论与实践
分子轨道与波函数方法 · 量子蒙特卡罗 · Gaussian · VASP