Skip to content

Orbital Rotation Parameterization

A converged set of SCF molecular orbitals \( \{C_0\} \) can be rotated into a nearby set through a unitary transformation generated by an anti-Hermitian matrix \( \kappa \):

\[ C(\kappa) = C_0 \, e^{\kappa}, \qquad \kappa^{\dagger} = -\kappa \]

Only rotations that mix an occupied orbital \( i \) with a virtual orbital \( a \) change the electronic energy (occupied–occupied and virtual–virtual rotations merely relabel orbitals within the same subspace). The independent parameters are therefore the occupied–virtual block elements \( \kappa_{ai} \).

To second order, the energy as a function of these parameters is

\[ E(\kappa) = E_0 + \sum_{ai} \frac{\partial E}{\partial \kappa_{ai}} \kappa_{ai} + \frac{1}{2} \sum_{ai,bj} \frac{\partial^2 E}{\partial \kappa_{ai}\, \partial \kappa_{bj}} \kappa_{ai}\kappa_{bj} + \dots \]

At SCF convergence, the first derivatives vanish identically:

\[ \frac{\partial E}{\partial \kappa_{ai}} = F_{ai} = 0 \]

This is Brillouin's theorem — it guarantees a stationary point, but says nothing about the sign of the second-derivative term. That second-derivative matrix,

\[ H_{ai,bj} = \frac{\partial^2 E}{\partial \kappa_{ai}\, \partial \kappa_{bj}} \]

is exactly the electronic (orbital) Hessian, also called the stability matrix, examined by a stable calculation.