问题解决了。SCF的一般流程是：C->F->C->F......计算梯度应该在F->C之间，而不是C->F之间，因为在后一阶段，现有的C和现有的F已经是是恰好符合HFR方程的，Fmo的非对角元肯定为0，而前一阶段则不然。

当你的SCF收敛了的时候，轨道梯度就是0啊，没什么问题。没收敛的时候，也不满足HFR方程，此时(C^T)*Fao*C的OV block应该不是0，如果是0应该是你写错了。