Engineering topics

Carrington, Tucker

Publications and source records attributed to Carrington, Tucker.

4 records

DOE OSTI2016

Using a pruned, nondirect product basis in conjunction with the multi-configuration time-dependent Hartree (MCTDH) method

In this paper, we propose a pruned, nondirect product multi-configuration time dependent Hartree (MCTDH) method for solving the Schrödinger equation. MCTDH uses optimized 1D basis functions, called single particle functions, but the size of the standard direct product MCTDH basis scales exponentially with D, the number of coordinates. We compare the pruned approach to standard MCTDH calculations for basis sizes small enough that the latter are possible and demonstrate that pruning the basis reduces the CPU cost of computing vibrational energy levels of acetonitrile (D = 12) by more than two orders of magnitude. Using the pruned method, it is possible to do calculations with larger bases, for which the cost of standard MCTDH calculations is prohibitive. Pruning the basis complicates the evaluation of matrix-vector products. In this paper, they are done term by term for a sum-of-products Hamiltonian. When no attempt is made to exploit the fact that matrices representing some of the factors of a term are identity matrices, one needs only to carefully constrain indices. In this paper, we develop new ideas that make it possible to further reduce the CPU time by exploiting identity matrices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2016

Assessing the utility of phase-space-localized basis functions: Exploiting direct product structure and a new basis function selection procedure

In this paper we show that it is possible to use an iterative eigensolver in conjunction with Halverson and Poirier’s symmetrized Gaussian (SG) basis [T. Halverson and B. Poirier, J. Chem. Phys. 137, 224101 (2012)] to compute accurate vibrational energy levels of molecules with as many as five atoms. This is done, without storing and manipulating large matrices, by solving a regular eigenvalue problem that makes it possible to exploit direct-product structure. These ideas are combined with a new procedure for selecting which basis functions to use. The SG basis we work with is orders of magnitude smaller than the basis made by using a classical energy criterion. We find significant convergence errors in previous calculations with SG bases. For sum-of-product Hamiltonians, SG bases large enough to compute accurate levels are orders of magnitude larger than even simple pruned bases composed of products of harmonic oscillator functions.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2016

Calculated vibrational states of ozone up to dissociation

A new accurate global potential energy surface for the ground electronic state of ozone [R. Dawes et al., J. Chem. Phys. 139, 201103 (2013)] was published fairly recently. The topography near dissociation differs significantly from previous surfaces, without spurious submerged reefs and corresponding van der Waals wells. This has enabled significantly improved descriptions of scattering processes, capturing the negative temperature dependence and large kinetic isotope effects in exchange reaction rates. The exchange reactivity was found to depend on the character of near-threshold resonances and their overlap with reactant and product wavefunctions, which in turn are sensitive to the potential. We present global “three-well” calculations of all bound vibrational states of three isotopic combinations of ozone ( 48 O 3 , 16 O 2 18 O, 16 O 2 17 O) for J = 0 and J = 1 with a focus on the character and density of highly excited states and discuss their impact on the ozone isotopic anomaly. The calculations were done using a parallel symmetry-adapted Lanczos method with the RV3 code. Some comparisons were made with results obtained with the improved relaxation method implemented in the Heidelberg multi-configuration time-dependent Hartree code.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2015

A multi-dimensional Smolyak collocation method in curvilinear coordinates for computing vibrational spectra

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.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗