Ghost levels and near-variational forms of the discrete variable representation: Application to H2O

A theoretical study on the infrared signatures of proton-bound rare gas dimers (Rg-H+-Rg), Rg = {Ne, Ar, Kr, and Xe}

The infrared spectrum of proton-bound rare gas dimers has been extensively studied via matrix isolation spectroscopy. However, little attention has been paid on their spectrum in the gas phase. Most of the Rg₂H⁺ has not been detected outside the matrix environment. Recently, Ar_nH⁺ (n = 3-7) has been first detected in the gas-phase [D. C. McDonald et al., J. Chem. Phys. 145, 231101 (2016)]. In that work, anharmonic theory can reproduce the observed vibrational structure. In this paper, we extend the existing theory to examine the vibrational signatures of Rg₂H⁺, Rg = {Ne, Ar, Kr, and Xe}. The successive binding of Rg to H⁺ was investigated through the calculation of stepwise formation energies. It was found that this binding is anti-cooperative. High-level full-dimensional potential energy surfaces at the CCSD(T)/aug-cc-pVQZ//MP2/aug-cc-pVQZ were constructed and used in the anharmonic calculation via discrete variable representation. We found that the potential coupling between the symmetric and asymmetric Rg-H⁺ stretch (ν₁ and ν₃ respectively) causes a series of bright n₁ν₁ + ν₃ progressions. From Ne₂H⁺ to Xe₂H⁺, an enhancement of intensities for these bands was observed.

Accuracy of a hybrid finite-element method for solving a scattering Schrödinger equation

This hybrid method [finite-element discrete variable representation (FE-DVR)], introduced by Resigno and McCurdy [Phys. Rev. A 62, 032706 (2000)], uses Lagrange polynomials in each partition, rather than "hat" functions or Gaussian functions. These polynomials are discrete variable representation functions, and they are orthogonal under the Gauss-Lobatto quadrature discretization approximation. Accuracy analyses of this method are performed for the case of a one-dimensional Schrödinger equation with various types of local and nonlocal potentials for scattering boundary conditions. The accuracy is ascertained by a comparison with a spectral Chebyshev integral equation method, accurate to 1:10⁻¹¹. For an accuracy of the phase shift of 1:10⁻⁸, the FE-DVR method is found to be 100 times faster than a sixth-order finite-difference method (Numerov), it is easy to program, and it can routinely achieve an accuracy of better than 1:10⁻⁸ for the numerical examples studied.

Variational properties of the discrete variable representation: discrete variable representation via effective operators

A variational finite basis representation/discrete variable representation (FBR/DVR) Hamiltonian operator has been introduced. By calculating its matrix elements exactly one obtains, depending on the choice of the basis set, either a variational FBR or a variational DVR. The domain of grid points on which the FBR/DVR is variational has been shown to consist of the subsets of the set of grid points one obtains by diagonalizing commuting variational basis representations of the coordinate operators. The variational property implies that the optimal of the subsets of a fixed number of points, i.e., the subset which gives the possible highest accuracy eigenpairs, gives the DVR of the smallest trace. The symmetry properties of the variational FBR/DVR Hamiltonian operator are analyzed and methods to incorporate symmetry into FBR/DVR calculations are discussed. It is shown how the Fourier-basis FBR/DVR suitable to solving periodic systems arise within the theory presented. Numerical examples are given to illustrate the theoretical results. The use of variational effective Hamiltonian and coordinate operators has been instrumental in this study. They have been introduced in a novel way by exploiting quasi-Hermiticity.

Spectral convergence of the quadrature discretization method in the solution of the Schrödinger and Fokker-Planck equations: Comparison with sinc methods

Spectral methods based on nonclassical polynomials and Fourier basis functions or sinc interpolation techniques are compared for several eigenvalue problems for the Fokker-Planck and Schrodinger equations. A very rapid spectral convergence of the eigenvalues versus the number of quadrature points is obtained with the quadrature discretization method (QDM) and the appropriate choice of the weight function. The QDM is a pseudospectral method and the rate of convergence is compared with the sinc method reported by Wei [J. Chem. Phys., 110, 8930 (1999)]. In general, sinc methods based on Fourier basis functions with a uniform grid provide a much slower convergence. The paper considers Fokker-Planck equations (and analogous Schrodinger equations) for the thermalization of electrons in atomic moderators and for a quartic potential employed to model chemical reactions. The solution of the Schrodinger equation for the vibrational states of I2 with a Morse potential is also considered.

Computing resonance energies, widths, and wave functions using a Lanczos method in real arithmetic

We introduce new ideas for calculating resonance energies and widths. It is shown that a non-Hermitian-Lanczos approach can be used to compute eigenvalues of H+W, where H is the Hamiltonian and W is a complex absorbing potential (CAP), without evaluating complex matrix-vector products. This is done by exploiting the link between a CAP-modified Hamiltonian matrix and a real but nonsymmetric matrix U suggested by Mandelshtam and Neumaier [J. Theor. Comput. Chem. 1, 1 (2002)] and using a coupled-two-term Lanczos procedure. We use approximate resonance eigenvectors obtained from the non-Hermitian-Lanczos algorithm and a very good CAP to obtain very accurate energies and widths without solving eigenvalue problems for many values of the CAP strength parameter and searching for cusps. The method is applied to the resonances of HCO. We compare properties of the method with those of established approaches.

Contracted basis Lanczos methods for computing numerically exact rovibrational levels of methane

We present a numerically exact calculation of rovibrational levels of a five-atom molecule. Two contracted basis Lanczos strategies are proposed. The first and preferred strategy is a two-stage contraction. Products of eigenfunctions of a four-dimensional (4D) stretch problem and eigenfunctions of 5D bend-rotation problems, one for each K, are used as basis functions for computing eigenfunctions and eigenvalues (for each K) of the Hamiltonian without the Coriolis coupling term, denoted H0. Finally, energy levels of the full Hamiltonian are calculated in a basis of the eigenfunctions of H0. The second strategy is a one-stage contraction in which energy levels of the full Hamiltonian are computed in the product contracted basis (without first computing eigenfunctions of H0). The two-stage contraction strategy, albeit more complicated, has the crucial advantage that it is trivial to parallelize the calculation so that the CPU and memory costs are independent of J. For the one-stage contraction strategy the CPU and memory costs of the difficult part of the calculation scale linearly with J. We use the polar coordinates associated with orthogonal Radau vectors and spherical harmonic type rovibrational basis functions. A parity-adapted rovibrational basis suitable for a five-atom molecule is proposed and employed to obtain bend-rotation eigenfunctions in the first step of both contraction methods. The effectiveness of the two methods is demonstrated by calculating a large number of converged J = 1 rovibrational levels of methane using a global potential energy surface.

New results for the OH (ν=0,j=0)+CO (ν=0,j=0)→H+CO2 reaction: Five- and full-dimensional quantum dynamical study on several potential energy surfaces

Full- [six-dimensional (6-D)] and reduced-dimensional [five-dimensional (5-D)] quantum wave packet calculations have been performed for the title reaction to obtain reaction probabilities deriving from the ground rovibrational states of OH and CO with total angular momentum J = 0. Three potential energy surfaces (PES) are studied, namely, those of Bradley and Schatz (BS), Yu, Muckerman, and Sears (YMS), and Lakin, Troya, Schatz, and Harding (LTSH). 6-D calculations are performed only for the BS PES, while 5-D results are reported for all three PES'. The 6-D results obtained in the present work improve on those previously reported, since a larger vibrational basis and a better representation of the OH and CO bonds has been introduced. In particular, we now employ a generalized Lanczos-Morse discrete variable representation for both the OH and CO vibrations. In a further improvement, the generalized discrete variable representation of the CO vibration is based on different CO intramolecular potentials for the asymptotic and product grids employed in our projection formalism. This new treatment of the vibrational bases allows for a large reduction in computation time with respect to our previous implementation of the wave packet method, for a given level of accuracy. As a result, we have been able to extend the range of collision energies for which we can obtain converged 6-D results to a higher energy (0.8 eV) than was possible before (0.5 eV). The comparison of the new 6-D and previous 5-D results for the BS PES shows good agreement of the general trend in the reaction probabilities over all collision energies considered (0.1-0.8 eV), while our previous 6-D calculation showed reaction probabilities that differed from the 5-D results by up to 10% between 0.5 and 0.8 eV. The 5-D reaction probabilities reveal interesting trends for the different PES'. In particular, at low energies (< 0.2 eV) the LTSH PES gives rise to much larger reactivity than the other PES', while at high energies (> 0.3 eV) its reaction probability decreases with respect to the BS and YMS PES', being more than a factor of 2 smaller at 0.8 eV. A 5-D calculation on a modified version of the LTSH surface shows that the van der Waals interaction in the entrance channel, which is not correctly described in the other PES' is largely responsible for its larger reactivity at low energies. The large difference between the 5-D reaction probabilities for the YMS and LTSH PES' serves to emphasize the importance of the van der Waals interaction for the reactivity at low energies, because most of the stationary point energies on the YMS and LTSH PES are rather similar, being in line with high-level ab initio information.

Mapped grid methods for long-range molecules and cold collisions

The paper discusses ways of improving the accuracy of numerical calculations for vibrational levels of diatomic molecules close to the dissociation limit or for ultracold collisions, in the framework of a grid representation. In order to avoid the implementation of very large grids, Kokoouline et al. [J. Chem. Phys. 110, 9865 (1999)] have proposed a mapping procedure through introduction of an adaptive coordinate x subjected to the variation of the local de Broglie wavelength as a function of the internuclear distance R. Some unphysical levels ("ghosts") then appear in the vibrational series computed via a mapped Fourier grid representation. In the present work the choice of the basis set is reexamined, and two alternative expansions are discussed: Sine functions and Hardy functions. It is shown that use of a basis set with fixed nodes at both grid ends is efficient to eliminate "ghost" solutions. It is further shown that the Hamiltonian matrix in the sine basis can be calculated very accurately by using an auxiliary basis of cosine functions, overcoming the problems arising from numerical calculation of the Jacobian J(x) of the R-->x coordinate transformation.

Solving the bound-state Schrodinger equation by reproducing kernel interpolation

Based on reproducing kernel Hilbert space theory and radial basis approximation theory, a grid method is developed for numerically solving the N-dimensional bound-state Schrodinger equation. Central to the method is the construction of an appropriate bounded reproducing kernel (RK) Lambda(alpha)( ||r ||) from the linear operator -nabla(2)(r)+lambda(2) where nabla(2)(r) is the N-dimensional Laplacian, lambda>0 is a parameter related to the binding energy of the system under study, and the real number alpha>N. The proposed (Sobolev) RK Lambda(alpha)(r,r(')) is shown to be a positive-definite radial basis function, and it matches the asymptotic solutions of the bound-state Schrodinger equation. Numerical tests for the one-dimensional (1D) Morse potential and 2D Henon-Heiles potential reveal that the method can accurately and efficiently yield all the energy levels up to the dissociation limit. Comparisons are also made with the results based on the distributed Gaussian basis method in the 1D case as well as the distributed approximating functional method in both 1D and 2D cases.