Skip to content

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.