We demonstrate that one can exhaustively determine the n‐bound eigenstates of a Hamiltonian H by constructing a sequence of supersymmetric (SUSY) partner Hamiltonians and invoking a time‐dependent quantum adiabatic switching algorithm for passage from the ground state of one to the other. The ground states of the initial pair H(0) and H(1) are constructed by solving the Riccati equation for the superpotential ϕ(0) for H(0) and adiabatically switching from the ground state Ψ 0(0) of H(0) to the ground state Ψ 0(1) of H(1). The charge operator Q 0+ is then used to recover the first excited state Ψ 1(0) of H(0). The procedure is repeated for the ground states of SUSY pairs H(n + 1) and H(n + 2), and appropriate charge operators lead to the excited states Ψ italicn + 2(0) of H(0) with
$n = 0,1,2, \cdots.$, thereby exhausting the full eigenspectrum of H(0). The workability of the proposed method is shown with several well‐known examples. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011