Mathematics for Nerds: Wavefunction Stability¶
This section derives the theory behind the stability test, starting from the orbital-rotation parameterization of the SCF energy and ending with the full iterative algorithm Gaussian uses to check it.
- Orbital Rotation Parameterization, how the wavefunction is expanded as a function of occupied–virtual orbital mixing.
- The Electronic Hessian (A and B Matrices), the second-derivative matrix that defines stability, and its link to RPA/TDHF theory.
- Types of Wavefunction Instabilities, real vs. complex, singlet vs. triplet.
- The Stability Eigenvalue Problem, why the static stability test reduces to a real, symmetric eigenvalue problem.
- Connection to RPA/TDHF Equations, why stability analysis and excited-state theory share the same machinery.
- Davidson Iterative Diagonalization, how Gaussian finds the lowest eigenvalues without ever forming the full Hessian.
- Complete Wavefunction Stability Algorithm, the full procedure, start to finish.