DETERMINING QUANTUM BOUND-STATE EIGENVALUES AND EIGENVECTORS AS FUNCTIONS OF PARAMETERS IN THE HAMILTONIAN - AN EFFICIENT EVOLUTIONARY APPROACH

This paper addresses the problem of finding the quantum bound-state en ergy eigenvalues and eigenvectors as functions of a set of continuous parameters characterizing a Hamiltonian. A recent paper introduced a p arametric equations of motion (PEM) method for this purpose, and the p resent work extends the method to allow for the analysis of a single e nergy level and its wavefunction. After solving the Schrodinger equati on for its nth eigenvalue and eigenvector, evaluated at a reference va lue of the Hamiltonian's parameters, the differential equations of the single-state PEM (ss-PEM) method are used to propagate the nth energy level and its eigenfunction through the entire parameter space of the Hamiltonian. The new ss-PEM method, which reduces the number of diffe rential equations to be solved, appears more efficient than diagonaliz ation when the energy is sought at a moderate number of values for the parameters in the Hamiltonian. The PEM methods are extended to treat non-orthogonal basis sets that facilitate more rapid convergence of th e solutions. The energy of the ss-PEM, which is always an upper bound to the true energy, is exact in the limit of a complete basis set. Con nections of the method are made to linear variational calculations, Da lgarno-Lewis perturbation theory and the original PEM methods. Sets of non-orthogonal Chebyshev polynomials are employed in illustrations of the ss-PEM method to determine (a) the ground-state energy as a funct ion of internuclear separation in the hydrogen molecular ion, and (b) the ground-state energy of two electron ions as a function of nuclear charge. The calculation with the two-electron ions involves two parame ters, the nuclear charge and a basis set parameter that influences the distribution of the nodes of the Chebyshev basis functions. Evolution of the basis set parameter to improve the energies of the ions sugges ts an additional application of the ss-PEM method in which quantum ene rgies are minimized with respect to nonlinear basis set parameters. Th e ss-PEM method offers an effective tool for mapping the solutions of the Schrodinger equation as a function of model parameters in the Hami ltonian.