Search

The present work intends to join and respond to the excellent and thoroughly documented rovibrational study of X. G. Wang and T. Carrington, Jr. [J. Chem. Phys. 154, 124112 (2021)] that used an approach tailored for floppy dimers with an analytic dimer Hamiltonian and a non-product basis set including Wigner D functions. It is shown in the present work that the GENIUSH black-box-type rovibrational method can approach the performance of the tailor-made computation for the example of the floppy methane–water dimer. Rovibrational transition energies and intensities are obtained in the black-box-type computation with a twice as large basis set and in excellent numerical agreement in comparison with the more efficient tailor-made approach. [ABSTRACT FROM AUTHOR]

The complex of a methane molecule and a fluoride anion represents a 12-dimensional (12D), four-well vibrational problem with multiple large-amplitude motions, which has challenged the quantum dynamics community for years. The present work reports vibrational band origins and tunneling splittings obtained in a full-dimensional variational vibrational computation using the GENIUSH program and the Smolyak quadrature scheme. The converged 12D vibrational band origins and tunneling splittings confirm complementary aspects of the earlier full- and reduced-dimensionality studies: (1) the tunneling splittings are smaller than 0.02 cm−1; (2) a single-well treatment is not sufficient (except perhaps the zero-point vibration) due to a significant anharmonicity over the wells; and thus, (3) a full-dimensional treatment appears to be necessary. The present computations extend to a higher energy range than earlier work, show that the tunneling splittings increase upon vibrational excitation of the complex, and indicate non-negligible "heavy-atom" tunneling. [ABSTRACT FROM AUTHOR]

Methodological progress is reported in the challenging direction of a black-box-type variational solution of the (ro)vibrational Schrödinger equation applicable to floppy, polyatomic systems with multiple large-amplitude motions. This progress is achieved through the combination of (i) the numerical kinetic-energy operator (KEO) approach of Mátyus et al. [J. Chem. Phys. 130, 134112 (2009)] and (ii) the Smolyak nonproduct grid method of Avila and Carrington, Jr. [J. Chem. Phys. 131, 174103 (2009)]. The numerical representation of the KEO makes it possible to choose internal coordinates and a body-fixed frame best suited for the molecular system. The Smolyak scheme reduces the size of the direct-product grid representation by orders of magnitude, while retaining some of the useful features of it. As a result, multidimensional (ro)vibrational states are computed with system-adapted coordinates, a compact basis- and grid-representation, and an iterative eigensolver. Details of the methodological developments and the first numerical applications are presented for the CH4·Ar complex treated in full (12D) vibrational dimensionality. [ABSTRACT FROM AUTHOR]

In this paper, we demonstrate that it is possible to apply collocation to compute vibrational energy levels of a five-atom molecule using an exact kinetic energy operator (with cross terms and coordinatedependent coefficients). This is made possible by using (1) a pruned basis of products of univariate functions; (2) a Smolyak grid made from nested sequences of grids for each coordinate; (3) a collocation method that obviates the need to solve a generalized eigenvalue problem; (4) an efficient sequential transformation between the (nondirect product) grid and the (nondirect product) basis representations; and (5) hierarchical univariate functions that make it possible to avoid storing large intermediate vectors. The accuracy of the method is confirmed by computing 500 vibrational energy levels of methane. [ABSTRACT FROM AUTHOR]

In this paper, we improve the collocation method for computing vibrational spectra that was presented in the work of Avila and Carrington, Jr. [J. Chem. Phys. 143, 214108 (2015)]. Known quadrature and collocation methods using a Smolyak grid require storing intermediate vectors with more elements than points on the Smolyak grid. This is due to the fact that grid labels are constrained among themselves and basis labels are constrained among themselves. We show that by using the so-called hierarchical basis functions, one can significantly reduce the memory required. In this paper, the intermediate vectors have only as many elements as the Smolyak grid. The ideas are tested by computing energy levels of CH2NH. [ABSTRACT FROM AUTHOR]

In this paper, we improve the collocation method for computing vibrational spectra that was presented in Avila and Carrington, Jr. [J. Chem. Phys. 139, 134114 (2013)]. Using an iterative eigensolver, energy levels and wavefunctions are determined from values of the potential on a Smolyak grid. The kinetic energy matrix-vector product is evaluated by transforming a vector labelled with (nondirect product) grid indices to a vector labelled by (nondirect product) basis indices. Both the transformation and application of the kinetic energy operator (KEO) scale favorably. Collocation facilitates dealing with complicated KEOs because it obviates the need to calculate integrals of coordinate dependent coefficients of differential operators. The ideas are tested by computing energy levels of HONO using a KEO in bond coordinates. [ABSTRACT FROM AUTHOR]

In this paper we propose and test a method for computing numerically exact vibrational energy levels of a molecule with six atoms. We use a pruned product basis, a non-product quadrature, the Lanczos algorithm, and the exact normal-coordinate kinetic energy operator (KEO) with the πtμπ term. The Lanczos algorithm is applied to a Hamiltonian with a KEO for which μ is evaluated at equilibrium. Eigenvalues and eigenvectors obtained from this calculation are used as a basis to obtain the final energy levels. The quadrature scheme is designed, so that integrals for the most important terms in the potential will be exact. The procedure is tested on C2H4. All 12 coordinates are treated explicitly. We need only ∼1.52 × 108 quadrature points. A product Gauss grid with which one could calculate the same energy levels has at least 5.67 × 1013 points. [ABSTRACT FROM AUTHOR]

In this paper we propose a new quadrature scheme for computing vibrational spectra and apply it, using a Lanczos algorithm, to CH3CN. All 12 coordinates are treated explicitly. We need only 157'419'523 quadrature points. It would not be possible to use a product Gauss grid because 33 853 318 889 472 product Gauss points would be required. The nonproduct quadrature we use is based on ideas of Smolyak, but they are extended so that they can be applied when one retains basis functions θn1(r1)...θnD(rD) that satisfy the condition α1n1 + ··· + αDnD <= b, where the αk are integers. We demonstrate that it is possible to exploit the structure of the grid to efficiently evaluate the matrix-vector products required to use the Lanczos algorithm. [ABSTRACT FROM AUTHOR]

The size of the quadrature grid required to compute potential matrix elements impedes solution of the vibrational Schrödinger equation if the potential does not have a simple form. This quadrature grid-size problem can make computing (ro)vibrational spectra impossible even if the size of the basis used to construct the Hamiltonian matrix is itself manageable. Potential matrix elements are typically computed with a direct product Gauss quadrature whose grid size scales as ND, where N is the number of points per coordinate and D is the number of dimensions. In this article we demonstrate that this problem can be mitigated by using a pruned basis set and a nonproduct Smolyak grid. The constituent 1D quadratures are designed for the weight functions important for vibrational calculations. For the SF6 stretch problem (D=6) we obtain accurate results with a grid that is more than two orders of magnitude smaller than the direct product Gauss grid. If D>6 we expect an even bigger reduction. [ABSTRACT FROM AUTHOR]

In this paper, we present a new collocation method for solving the Schroedinger equation. Collocation has the advantage that it obviates integrals. All previous collocation methods have, however, the crucial disadvantage that they require solving a generalized eigenvalue problem. By combining Lagrange-like functions with a Smolyak interpolant, we device a collocation method that does not require solving a generalized eigenvalue problem. We exploit the structure of the grid to develop an efficient algorithm for evaluating the matrix-vector products required to compute energy levels and wavefunctions. Energies systematically converge as the number of points and basis functions are increased. [ABSTRACT FROM AUTHOR]

In this paper, we present new basis pruning schemes and compatible quadrature grids for solving the vibrational Schrödinger equation. The new basis is designed to include the product basis functions coupled by the largest terms in the potential and important for computing low-lying vibrational levels. To solve the vibrational Schrödinger equation without approximating the potential, one must use quadrature to compute potential matrix elements. For a molecule with more than five atoms, the use of iterative methods is imperative, due to the size of the basis and the quadrature grid. When using iterative methods in conjunction with quadrature, it is important to evaluate matrix-vector products by doing sums sequentially. This is only possible if both the basis and the grid have structure. Although it is designed to include only functions coupled by the largest terms in the potential, the new basis and also the quadrature for doing integrals with the basis have enough structure to make efficient matrix-vector products possible. When results obtained with a multimode approximation to the potential are accurate enough, full-dimensional quadrature is not necessary. Using the quadrature methods of this paper, we evaluate the accuracy of calculations made by making multimode approximations. [ABSTRACT FROM AUTHOR]

The complex of a methane molecule and a fluoride anion represents a 12-dimensional (12D), four-well vibrational problem with multiple large-amplitude motions, which has challenged the quantum dynamics community for years. The present work reports vibrational band origins and tunneling splittings obtained in a full-dimensional variational vibrational computation using the GENIUSH program and the Smolyak quadrature scheme. The converged 12D vibrational band origins and tunneling splittings confirm complementary aspects of the earlier full- and reduced-dimensionality studies: (1) the tunneling splittings are smaller than 0.02 cm -1 ; (2) a single-well treatment is not sufficient (except perhaps the zero-point vibration) due to a significant anharmonicity over the wells; and thus, (3) a full-dimensional treatment appears to be necessary. The present computations extend to a higher energy range than earlier work, show that the tunneling splittings increase upon vibrational excitation of the complex, and indicate non-negligible "heavy-atom" tunneling.

In this paper, we demonstrate that it is possible to apply collocation to compute vibrational energy levels of a five-atom molecule using an exact kinetic energy operator (with cross terms and coordinate-dependent coefficients). This is made possible by using (1) a pruned basis of products of univariate functions; (2) a Smolyak grid made from nested sequences of grids for each coordinate; (3) a collocation method that obviates the need to solve a generalized eigenvalue problem; (4) an efficient sequential transformation between the (nondirect product) grid and the (nondirect product) basis representations; and (5) hierarchical univariate functions that make it possible to avoid storing large intermediate vectors. The accuracy of the method is confirmed by computing 500 vibrational energy levels of methane.

In this paper, we improve the collocation method for computing vibrational spectra that was presented in the work of Avila and Carrington, Jr. [J. Chem. Phys. 143, 214108 (2015)]. Known quadrature and collocation methods using a Smolyak grid require storing intermediate vectors with more elements than points on the Smolyak grid. This is due to the fact that grid labels are constrained among themselves and basis labels are constrained among themselves. We show that by using the so-called hierarchical basis functions, one can significantly reduce the memory required. In this paper, the intermediate vectors have only as many elements as the Smolyak grid. The ideas are tested by computing energy levels of CH 2 NH.

We propose a new method for obtaining potential energy surfaces in sum-of-products (SOP) form. If the number of terms is small enough, a SOP potential surface significantly reduces the cost of quantum dynamics calculations by obviating the need to do multidimensional integrals by quadrature. The method is based on a Smolyak interpolation technique and uses polynomial-like or spectral basis functions and 1D Lagrange-type functions. When written in terms of the basis functions from which the Lagrange-type functions are built, the Smolyak interpolant has only a modest number of terms. The ideas are tested for HONO (nitrous acid).