Using simultaneous diagonalization and trace minimization to make an efficient and simple multidimensional basis for solving the vibrational Schrodinger equation.
In this paper we improve the product simultaneous diagonalization (SD) basis method we previously proposed [J. Chem. Phys. 122, 134101 (2005)] and applied to solve the Schrodinger equation for the motion of nuclei on a potential surface. The improved method is tested using coupled complicated Hamiltonians with as many as 16 coordinates for which we can easily find numerically exact solutions. In a basis of sorted products of one-dimensional (1D) SD functions the Hamiltonian matrix is nearly diagonal. The localization of the 1D SD functions for coordinate qc depends on a parameter we denote alphac. In this paper we present a trace minimization scheme for choosing alphac to nearly block diagonalize the Hamiltonian matrix. Near-block diagonality makes it possible to truncate the matrix without degrading the accuracy of the lowest energy levels. We show that in the sorted product SD basis perturbation theory works extremely well. The trace minimization scheme is general and easy to implement.
Département de chimie, Université de Montréal, C.P. 6128, succursale Centre-ville, Montréal (Québec) H3C 3J7, Canada. r.dawes@umontreal.ca
Richard
R
Dawes
Tucker
T
Carrington