Multimode wavelet basis calculations via the molecular self-consistent-field plus configuration-interaction method

Optimization of numerical orbitals using the Helmholtz kernel

We present an integration scheme for optimizing the orbitals in numerical electronic structure calculations on general molecules. The orbital optimization is performed by integrating the Helmholtz kernel in the double bubble and cube basis, where bubbles represent the steep part of the functions in the vicinity of the nuclei, whereas the remaining cube part is expanded on an equidistant three-dimensional grid. The bubbles' part is treated by using one-center expansions of the Helmholtz kernel in spherical harmonics multiplied with modified spherical Bessel functions of the first and second kinds. The angular part of the bubble functions can be integrated analytically, whereas the radial part is integrated numerically. The cube part is integrated using a similar method as we previously implemented for numerically integrating two-electron potentials. The behavior of the integrand of the auxiliary dimension introduced by the integral transformation of the Helmholtz kernel has also been investigated. The correctness of the implementation has been checked by performing Hartree-Fock self-consistent-field calculations on H₂, H₂O, and CO. The obtained energies are compared with reference values in the literature showing that an accuracy of 10^-4 to 10^-7 E_h can be obtained with our approach.

Matrix-free application of Hamiltonian operators in Coifman wavelet bases

A means of evaluating the action of Hamiltonian operators on functions expanded in orthogonal compact support wavelet bases is developed, avoiding the direct construction and storage of operator matrices that complicate extension to coupled multidimensional quantum applications. Application of a potential energy operator is accomplished by simple multiplication of the two sets of expansion coefficients without any convolution. The errors of this coefficient product approximation are quantified and lead to use of particular generalized coiflet bases, derived here, that maximize the number of moment conditions satisfied by the scaling function. This is at the expense of the number of vanishing moments of the wavelet function (approximation order), which appears to be a disadvantage but is shown surmountable. In particular, application of the kinetic energy operator, which is accomplished through the use of one-dimensional (1D) [or at most two-dimensional (2D)] differentiation filters, then degrades in accuracy if the standard choice is made. However, it is determined that use of high-order finite-difference filters yields strongly reduced absolute errors. Eigensolvers that ordinarily use only matrix-vector multiplications, such as the Lanczos algorithm, can then be used with this more efficient procedure. Applications are made to anharmonic vibrational problems: a 1D Morse oscillator, a 2D model of proton transfer, and three-dimensional vibrations of nitrosyl chloride on a global potential energy surface.

Two-dimensional quantum propagation using wavelets in space and time

A recent method for solving the time-dependent Schrodinger equation has been developed using expansions in compact-support wavelet bases in both space and time [H. Wang et al., J. Chem. Phys. 121, 7647 (2004)]. This method represents an exact quantum mixed time-frequency approach, with special initial temporal wavelets used to solve the initial value problem. The present work is a first extension of the method to multiple spatial dimensions applied to a simple two-dimensional (2D) coupled anharmonic oscillator problem. A wavelet-discretized version of norm preservation for time-independent Hamiltonians discovered in the earlier one-dimensional investigation is verified to hold as well in 2D and, by implication, in higher numbers of spatial dimensions. The wavelet bases are not restricted to rectangular domains, a fact which is exploited here in a 2D adaptive version of the algorithm.