Doubling of Chebyshev correlation function for calculating narrow resonances using low-storage filter diagonalization

Full-dimensional quantum mechanical calculations for the tunneling behavior of HOCO dissociation to H + CO2

The tunneling behavior during HOCO dissociation to H + CO₂ was investigated by full-dimensional quantum mechanical calculations based on an accurate global potential energy surface. The tunneling lifetimes for the low-lying 1500 vibrational states were calculated using the low-storage filter diagonalization method after a 1 million-step Chebyshev propagation. In the calculated energy range, the lifetimes of different vibrational states with similar energy are found to differ by 3-4 orders of magnitude, and the lower limit for these tunneling lifetimes is consistent with the experimental results reported by Continetti et al. For the given vibrational progressions, the lifetime of the vibrational states decreases with the increasing energy level, which is consistent with the results of 1D simulation calculations reported by Bowman, but the declining curve obviously fluctuates, and the declining slope is significantly different from that obtained by 1D simulation. Due to a difference in the effective barrier width, the mode-specific behavior of the tunneling effect is manifested in that the O-C-O' and H-O-C bend modes lead to the largest enhancement and an inhibitory effect on the tunneling process, respectively.

FURTHER GENERALIZATION AND NUMERICAL IMPLEMENTATION OF PSEUDO-TIME SCHRÖDINGER EQUATIONS FOR QUANTUM SCATTERING CALCULATIONS

We review and further develop the recently introduced numerical approach [Phys. Rev. Lett. 86, 5031, (2001)] for scattering calculations based on a so called pseudo-time Schrödinger equation, which is in turn a modification of the damped Chebyshev polynomial expansion scheme [J. Chem. Phys. 103, 2903, (1995)]. The method utilizes a special energy-dependent form for the absorbing potential in the time-independent Schrödinger equation, in which the complex energy spectrum is mapped inside the unit disk Ek → uk, where uk are the eigenvalues of some explicitly known sparse matrix U. Most importantly for the numerical implementation, all the physical eigenvalues uk are the extreme eigenvalues of U (i.e. |uk| ≈ 1 for resonances and |uk| = 1 for the bound states), which allows one to extract these eigenvalues very efficiently by harmonic inversion of a pseudo-time autocorrelation function y(t) = ϕ T Ut ϕ using the filter diagonalization method. The computation of y(t) up to time t = 2T requires only T sparse real matrix-vector multiplications. We describe and compare different schemes, effectively corresponding to different choices of the energy-dependent absorbing potential, and test them numerically by calculating resonances of the HCO molecule. Our numerical tests suggest an optimal scheme that provide accurate estimates for most resonance states using a single autocorrelation function.

Revelation of non-statistical behavior in HO2 vibration by a newab initiopotential energy surface

The hydroperoxyl radical (HO2) has long been considered as a prototype for statistical vibrational dynamics. In this work, however, it is shown that the bound state energy levels (up to the dissociation threshold) and low-lying resonances of the HO2 system (J=0) obtained on a new ab initio potential energy surface exhibit surprisingly large regularity. The implications of the non-statistical behavior of the HO2 system in unimolecular and bimolecular reactions are discussed.

HO2 Ro-Vibrational Bound-State Calculations for Large Angular Momentum: J = 30, 40, and 50

The Lanczos homogeneous filter diagonalization method and the real Chebyshev filter diagonalization scheme incorporating doubling of the autocorrelation functions have been employed to compute the HO2 ro-vibrational states for high total angular momenta, J = 30, 40, and 50. For such computationally challenging calculations, we have adopted a parallel computing strategy to perform the matrix-vector multiplications. Low-lying bound states and high-lying bound states close to the dissociation threshold are reported. For low-lying bound states, a spectroscopic assignment has been attempted and the widely used approximate J-shifting method has been tested for this deep-well system. For high-lying bound states, the attempted spectroscopic assignments as well as the J-shifting approximation fail because of very strong Coriolis mixing, indicating that the Coriolis couplings are important for this system.

Unimolecular rovibrational bound and resonance states for large angular momentum: J=20 calculations for HO2

We explore the calculation of unimolecular bound states and resonances for deep-well species at large angular momentum using a Chebychev filter diagonalization scheme incorporating doubling of the autocorrelation function as presented recently by Neumaier and Mandelshtam [Phys. Rev. Lett. 86, 5031 (2001)]. The method has been employed to compute the challenging J=20 bound and resonance states for the HO2 system. The methodology has firstly been tested for J=2 in comparison with previous calculations, and then extended to J=20 using a parallel computing strategy. The quantum J-specific unimolecular dissociation rates for HO2-->H+O2 in the energy range from 2.114 to 2.596 eV have been reported for the first time, and comparisons with the results of Troe and co-workers [J. Chem. Phys. 113, 11019 (2000) Phys. Chem. Chem. Phys. 2, 631 (2000)] from statistical adiabatic channel method/classical trajectory calculations have been made. For most of the energies, the reported statistical adiabatic channel method/classical trajectory rate constants agree well with the average of the fluctuating quantum-mechanical rates. Near the dissociation threshold, quantum rates fluctuate more severely, but their average is still in agreement with the statistical adiabatic channel method/classical trajectory results.

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.

Resonances of CH2(ãA11) and their roles in unimolecular and bimolecular reactions

Low-lying resonances of the CH2(a 1A1) system (J=0) in an accurate ab initio potential energy surface are studied using a filter-diagonalization method. The width of these resonances fluctuates by more than two orders of magnitude and on average increases with the energy. Analysis of the resonance states concludes that the unimolecular decay of the excited molecular system near the dissociation threshold is neither mode specific nor statistical state specific. This is apparently due to remnant regularity embedded in the largely chaotic classical phase space, as evidenced by periodic orbit analysis. As a result, the Rice-Ramsperger-Kassel-Marcus and statistical adiabatic channel models overestimate the average unimolecular decay rate. The implications of the resonances for the bimolecular C(1D)+H2 reaction are also discussed.

Converged quantum calculations of HO2 bound states and resonances for J=6 and 10

Bound and resonance states of HO(2) are calculated quantum mechanically using both the Lanczos homogeneous filter diagonalization method and the real Chebyshev filter diagonalization method for nonzero total angular momentum J=6 and 10, using a parallel computing strategy. For bound states, agreement between the two methods is quite satisfactory; for resonances, while the energies are in good agreement, the widths are in general agreement. The quantum nonzero-J specific unimolecular dissociation rates for HO(2) are also calculated.

Pseudotime Schrödinger equation with absorbing potential for quantum scattering calculations

The Schrödinger equation (Hpsi) (r) = [E+u(E)W(r)]psi(r) with an energy-dependent complex absorbing potential -u(E)W(r), associated with a scattering system, can be reduced for a special choice of u(E) to a harmonic inversion problem of a discrete pseudotime correlation function y(t) = phi(T)U(t)phi. An efficient formula for Green's function matrix elements is also derived. Since the exact propagation up to time 2t can be done with only approximately t real matrix-vector products, this gives an unprecedently efficient scheme for accurate calculations of quantum spectra for possibly very large systems.