Wavelets in curvilinear coordinate quantum calculations: H2+ electronic states

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.

Rovibrational spectroscopy calculations of neon dimer using a phase space truncated Weyl-Heisenberg wavelet basis

In a series of earlier articles [B. Poirier J. Theor. Comput. Chem. 2, 65 (2003); B. Poirier and A. Salam J. Chem. Phys. 121, 1690 (2004); B. Poirier and A. Salam J. Chem. Phys. 121, 1740 (2004)], a new method was introduced for performing exact quantum dynamics calculations in a manner that formally defeats exponential scaling with system dimensionality. The method combines an optimally localized, orthogonal Weyl-Heisenberg wavelet basis set with a simple phase space truncation scheme, and has already been applied to model systems up to 17 degrees of freedom (DOF's). In this paper, the approach is applied for the first time to a real molecular system (neon dimer), necessitating the development of an efficient numerical scheme for representing arbitrary potential energy functions in the wavelet representation. All bound rovibrational energy levels of neon dimer are computed, using both one DOF radial coordinate calculations and a three DOF Cartesian coordinate calculation. Even at such low dimensionalities, the approach is found to be competitive with another state-of-the-art method applied to the same system [J. Montgomery and B. Poirier J. Chem. Phys. 119, 6609 (2003)].

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

Wavelets provide potentially useful quantum bases for coupled anharmonic vibrational modes in polyatomic molecules as well as many other problems. A single compact support wavelet family provides a flexible basis with properties of orthogonality, localization, customizable resolution, and systematic improvability for general types of one-dimensional and separable systems. While direct product wavelet bases can be used in coupled multidimensional problems, exponential scaling of basis size with dimensionality ultimately provides limits on the number of coupled modes that can be treated simultaneously in exact quantum calculations. The molecular self-consistent-field plus configuration-interaction method is used here in multimode wavelet calculations to reduce the basis size without sacrificing flexibility or the ability to systematically control errors. Both two-dimensional Cartesian coordinate and three-dimensional curvilinear coordinate systems are examined with wavelets serving as universal bases in each case. The first example uses standard Daubechies [Ten Lectures on Wavelets (SIAM, Philadelphia (1992)] wavelets for each mode and the second adapts symmlet wavelets to intervals for each of the curvilinear coordinates.

Additions to the class of symmetric-antisymmetric multiwavelets: Derivation and use as quantum basis functions

Multiwavelet bases have been shown recently to apply to a variety of quantum problems. There are, however, only a few multiwavelet families that have been defined to date. Chui-Lian-type symmetric and antisymmetric multiwavelets are derived here that equal and exceed the polynomial interpolating power of previously available examples. Adaptations to domain edges are made with a view to use in curvilinear coordinate molecular calculations. The new highest-order multiwavelet family is shown to provide uniformly better performance for (i) basis representation of terms such as 1r(2) in near approach to the singularity at r=0 and (ii) eigenvalue calculation of a bending Hamiltonian taken from a curvilinear model of the ground-state vibrations of nitrosyl chloride.

Multiresolution quantum chemistry: Basic theory and initial applications

We describe a multiresolution solver for the all-electron local density approximation Kohn-Sham equations for general polyatomic molecules. The resulting solutions are obtained to a user-specified precision and the computational cost of applying all operators scales linearly with the number of parameters. The construction and use of separated forms for operators (here, the Green's functions for the Poisson and bound-state Helmholtz equations) enable practical computation in three and higher dimensions. Initial applications include the alkali-earth atoms down to strontium and the water and benzene molecules.

Quantum dynamics calculations using symmetrized, orthogonal Weyl-Heisenberg wavelets with a phase space truncation scheme. III. Representations and calculations

In a previous paper [J. Theo. Comput. Chem. 2, 65 (2003)], one of the authors (B.P.) presented a method for solving the multidimensional Schrodinger equation, using modified Wilson-Daubechies wavelets, and a simple phase space truncation scheme. Unprecedented numerical efficiency was achieved, enabling a ten-dimensional calculation of nearly 600 eigenvalues to be performed using direct matrix diagonalization techniques. In a second paper [J. Chem. Phys. 121, 1690 (2004)], and in this paper, we extend and elaborate upon the previous work in several important ways. The second paper focuses on construction and optimization of the wavelength functions, from theoretical and numerical viewpoints, and also examines their localization. This paper deals with their use in representations and eigenproblem calculations, which are extended to 15-dimensional systems. Even higher dimensionalities are possible using more sophisticated linear algebra techniques. This approach is ideally suited to rovibrational spectroscopy applications, but can be used in any context where differential equations are involved.

Quantum dynamics calculations using symmetrized, orthogonal Weyl-Heisenberg wavelets with a phase space truncation scheme. II. Construction and optimization

In this paper, we extend and elaborate upon a wavelet method first presented in a previous publication [B. Poirier, J. Theo. Comput. Chem. 2, 65 (2003)]. In particular, we focus on construction and optimization of the wavelet functions, from theoretical and numerical viewpoints, and also examine their localization properties. The wavelets used are modified Wilson-Daubechies wavelets, which in conjunction with a simple phase space truncation scheme, enable one to solve the multidimensional Schrodinger equation. This approach is ideally suited to rovibrational spectroscopy applications, but can be used in any context where differential equations are involved.

Multiscale quantum propagation using compact-support wavelets in space and time

Orthogonal compact-support Daubechies wavelets are employed as bases for both space and time variables in the solution of the time-dependent Schrodinger equation. Initial value conditions are enforced using special early-time wavelets analogous to edge wavelets used in boundary-value problems. It is shown that the quantum equations may be solved directly and accurately in the discrete wavelet representation, an important finding for the eventual goal of highly adaptive multiresolution Schrodinger equation solvers. While the temporal part of the basis is not sharp in either time or frequency, the Chebyshev method used for pure time-domain propagations is adapted to use in the mixed domain and is able to take advantage of Hamiltonian matrix sparseness. The orthogonal separation into different time scales is determined theoretically to persist throughout the evolution and is demonstrated numerically in a partially adaptive treatment of scattering from an asymmetric Eckart barrier.