10. 耦合簇、大小广延性与诊断
本课位置： 第 10/10 课。先完成前一课的推导；使用下方练习检查理解。
先修知识： 二次量子化与电子相关
1. 对连接激发算符取指数
耦合簇写 \(|\Psi\rangle=e^T|\Phi_0\rangle\)，其中 \(T=T_1+T_2+\cdots\) 是相对参考行列式的激发算符。CCSD 在 \(T\) 中保留连接单、双激发。指数展开为
因此 CCSD 的行列式展开含更高激发秩的非连接乘积，与 CISD 不同，也不等价于完整三、四激发。连接与非连接指代数贡献的组织方式，不是说真实粒子或片段是否在空间中连接。
2. 推导投影方程
代入薛定谔方程，左乘 \(e^{-T}\)，定义 \(\bar H=e^{-T}He^T\)。向参考及激发行列式投影得到
Baker-Campbell-Hausdorff 展开为 \(\bar H=H+[H,T]+\tfrac12[[H,T],T]+\cdots\)。标准有限基组双体哈密顿量及纯激发算符的连接展开在有限对易阶终止。所得振幅方程是非线性的，需要迭代；能量变化小不能代替振幅残差小。
3. 从因子分解理解大小广延性
非相互作用片段具有可分解参考，且 \(T=T_A+T_B\) 时，两个算符对易，因此
同时描述两个片段相关所需的非连接交叉乘积自动出现，满足相应条件时能量可加。这对不同大小分子的比较和解离极限重要，但不能保证单参考在各片段及整个断键路径中始终合理。
4. 完整例题：两个片段的展开
令各片段近似态为 \((1+tD_A)|\Phi_A\rangle\)、\((1+tD_B)|\Phi_B\rangle\)，其中 \(D\) 是双替换。乘积为 \(1+tD_A+tD_B+t^2D_AD_B\)。全局 CISD 缺少最后项；CC 双激发指数中的 \(\tfrac12(T_A+T_B)^2\) 自动产生它，因为两个对易交叉项抵消二分之一因子。这是结构性代数解释，不是实际执行的分子基准计算。
5. 精度、三激发与响应
CCSD(T) 在 CCSD 上增加特定微扰三激发修正，不是迭代 CCSDT。常规形式代价由 CCSD 的六次方增加到三激发修正的七次方。强近简并、大振幅或不稳定参考可使三激发微扰失效。应检查参考性质及诊断，不能把任意统一阈值当成准确性的证明。
相似变换后的哈密顿量通常非厄米，标准投影 CC 能量不具有变分上界。性质导数需要适当左态或拉格朗日响应，不能直接把参考 bra 当作相关 ket 的共轭。运动方程激发态方法也有自己的截断与状态性状限制。
6. 练习与解释
\(T_2^2/2\) 是连接四激发吗？ 不是，是连接双激发的乘积。独立的连接四激发自由度需要 \(T_4\)。
CC 能量比基准更低，能否据此判定更好？ 不能。投影 CC 无变分排序，需要系统截断、基组及参考态检验。
为什么报告残差范数？ 非线性迭代可能能量几乎不动，而振幅仍未满足投影方程。残差直接检验实际要求解的方程。
解析教学示意图，不是模拟数据。
进一步的概念检验
相关振幅依赖轨道选择，性质应结合完整近似解释。占据或虚空间旋转可改变数值组织，却不能补入遗漏的连接激发秩。局域相关阈值需要独立收敛检验，不能只写 CCSD(T) 标签就认为完全等价于正则结果。对于能量差，各结构的截断和局域设置应平衡，检查改变阈值后差值是否稳定。
7. 参考资料与学习联系
- DePrince group academic HF tutorial
- DePrince group academic tutorial collection
- Psi4 official coupled-cluster documentation
- Helgaker, Jørgensen and Olsen: Molecular Electronic-Structure Theory
本课推导及模型算术为原创教学构造；这里没有报告实际执行的分子模拟。
8. 关联理论与实践
量子力学 · 密度泛函理论 · 量子蒙特卡罗 · Gaussian