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.

Ask a Question

Have a question about quantum chemistry, or anything discussed on this site? Send it our way and we will reply by email.

Prefer to fill it out directly? Open the form.