Atomic spectral methods for molecular electronic structure calculations.

Abstract

Theoretical methods are reported for ab initio calculations of the adiabatic (Born-Oppenheimer) electronic wave functions and potential energy surfaces of molecules and other atomic aggregates. An outer product of complete sets of atomic eigenstates familiar from perturbation-theoretical treatments of long-range interactions is employed as a representational basis without prior enforcement of aggregate wave function antisymmetry. The nature and attributes of this atomic spectral-product basis are indicated, completeness proofs for representation of antisymmetric states provided, convergence of Schrodinger eigenstates in the basis established, and strategies for computational implemention of the theory described. A diabaticlike Hamiltonian matrix representative is obtained, which is additive in atomic-energy and pairwise-atomic interaction-energy matrices, providing a basis for molecular calculations in terms of the (Coulombic) interactions of the atomic constituents. The spectral-product basis is shown to contain the totally antisymmetric irreducible representation of the symmetric group of aggregate electron coordinate permutations once and only once, but to also span other (non-Pauli) symmetric group representations known to contain unphysical discrete states and associated continua in which the physically significant Schrodinger eigenstates are generally embedded. These unphysical representations are avoided by isolating the physical block of the Hamiltonian matrix with a unitary transformation obtained from the metric matrix of the explicitly antisymmetrized spectral-product basis. A formal proof of convergence is given in the limit of spectral closure to wave functions and energy surfaces obtained employing conventional prior antisymmetrization, but determined without repeated calculations of Hamiltonian matrix elements as integrals over explicitly antisymmetric aggregate basis states. Computational implementations of the theory employ efficient recursive methods which avoid explicit construction the metric matrix and do not require storage of the full Hamiltonian matrix to isolate the antisymmetric subspace of the spectral-product representation. Calculations of the lowest-lying singlet and triplet electronic states of the covalent electron pair bond (H(2)) illustrate the various theorems devised and demonstrate the degree of convergence achieved to values obtained employing conventional prior antisymmetrization. Concluding remarks place the atomic spectral-product development in the context of currently employed approaches for ab initio construction of adiabatic electronic eigenfunctions and potential energy surfaces, provide comparisons with earlier related approaches, and indicate prospects for more general applications of the method.