Numerical on-the-fly implementation of the action of the kinetic energy operator on a vibrational wave function: application to methanol

Variational Vibrational States of Methanol (12D)

Full-dimensional (12D) vibrational states of the methanol molecule (CH3OH) have been computed using the GENIUSH-Smolyak approach and the potential energy surface from Qu and Bowman (2013). All vibrational energies are converged better than 0.5 cm-1 with respect to the basis and grid size up to the first overtone of the CO stretch, ca. 2000 cm-1 beyond the zero-point vibrational energy. About 70 torsion-vibration states are reported and assigned. The computed vibrational energies agree with the available experimental data within less than a few cm-1 in most cases, which confirms the good accuracy of the potential energy surface. The computations are carried out using curvilinear normal coordinates with the option of path-following coefficients, which minimize the coupling of the small- and large-amplitude motions. It is important to ensure tight numerical fulfillment of the C3v(M) molecular symmetry for every geometry and coefficient set used to define the curvilinear normal coordinates along the torsional coordinate to obtain a faithful description of degeneracy in this floppy system. The reported values may provide a computational reference for fundamental spectroscopy, astrochemistry, and for the search of the proton-to-electron mass ratio variation using the methanol molecule.

Efficient vibrationally correlated calculations using n-mode expansion-based kinetic energy operators

Due to its efficiency and flexibility, the n-mode expansion is a frequently used tool for representing molecular potential energy surfaces in quantum chemical simulations. In this work, we investigate the performance of n-mode expansion-based models of kinetic energy operators in general polyspherical coordinate systems. In particular, we assess the operators with respect to accuracy in vibrationally correlated calculations and their effect on potential energy surface construction with the adaptive density guided approach. Our results show that the n-mode expansion-based operator variants are reliable and systematically improvable approximations of the full kinetic energy operator. Moreover, we introduce a workflow to generate the n-mode expanded kinetic energy operators on-the-fly within the adaptive density guided approach. This scheme can be applied in studies of species and coordinate systems, for which an analytical form of the kinetic energy operator is not available.

Vibrationally correlated calculations in polyspherical coordinates: Taylor expansion-based kinetic energy operators

The efficiency of quantum chemical simulations of nuclear motion can in many cases greatly benefit from the application of curvilinear coordinate systems. This is rooted in the fact that a set of smartly selected curvilinear coordinates may represent the motion naturally well, thus decreasing the couplings between motions in these coordinates. In this study, we assess the validity of different Taylor expansion-based approximations of kinetic energy operators in a (curvilinear) polyspherical parametrization. To this end, we investigate the accuracy as well as the numerical performance of the approximations in time-independent vibrational coupled cluster and full vibrational interaction calculations for several test cases ranging from tri- to penta-atomic molecules. We find that several of the proposed schemes reproduce the vibrational ground state and excitation energies to a decent accuracy, justifying their application in future investigations. Furthermore, due to the restricted mode coupling and their inherent sum-of-products form, the new approximations open up the possibility of treating large molecular systems with efficient vibrational coupled cluster schemes in general coordinates.

Smolyak Scheme for solving the Schrödinger equation: Application to Malonaldehyde in Full Dimensionality

In 1963 Smolyak introduced an approach to overcome the exponential scaling with respect to the number of variables of the direct product size [S. A. Smolyak Soviet Mathematics Doklady, 4, 240 (1963)]. The main idea is to replace a single large direct product by a sum of selected small direct products. It was first used in quantum dynamics in 2009 by Avila and Carrington [G. Avila and T. Carrington, J. Chem. Phys., 131, 174103 (2009)]. Since then, several calculations have been published by Avila and Carrington and by other groups. In the present study, and to push the limit to larger and more complex systems, this scheme is combined with the use of an on-the-fly calculation of the kinetic energy operator and a Block-Davidson procedure to obtain eigenstates in our home-made Fortran codes, ElVibRot and Tnum-Tana. This was applied to compute the tunneling splitting of malonaldehyde in full dimensionality (21D) using the potential of Mizukami et al. [W. Mizukami, S. Habershon, and D.P. Tew, J. Chem. Phys. 141, 1443-10 (2014)]. Our tunneling splitting calculations, 21.7±0.3 cm-1 and 2.9±0.1 cm-1 , show excellent agreement with the experimental values, 21.6 cm-1 and 2.9 cm-1 for the normal isotopologue and the mono-deuterated one, respectively.

A numerical-tensorial "hybrid" nuclear motion Hamiltonian and dipole moment operator for spectra calculation of polyatomic nonrigid molecules

The analysis and modeling of high-resolution spectra of nonrigid molecules require a specific Hamiltonian and group-theoretical formulation that differs significantly from that of more familiar rigid systems. Within the framework of Hougen-Bunker-Johns (HBJ) theory, this paper is devoted to the construction of a nonrigid Hamiltonian based on a suitable combination of numerical calculations for the nonrigid part in conjunction with the irreducible tensor operator method for the rigid part. For the first time, a variational calculation from ab initio potential energy surfaces is performed using the HBJ kinetic energy operator built from vibrational, large-amplitude motion, and rotational tensor operators expressed in terms of curvilinear and normal coordinates. Group theory for nonrigid molecules plays a central role in the characterization of the overall tunneling splittings and is discussed in the present approach. The construction of the dipole moment operator is also examined. Validation tests consisting of a careful convergence study of the energy levels as well as a comparison of results obtained from independent computer codes are given for the nonrigid molecules CH2, CH3, NH3, and H2O2. This work paves the way for the modeling of high-resolution spectra of larger nonrigid systems.

Exact quantum dynamics developments for floppy molecular systems and complexes

Molecular rotation, vibration, internal rotation, isomerization, tunneling, intermolecular dynamics of weakly and strongly interacting systems, intra-to-inter-molecular energy transfer, hindered rotation and hindered translation over surfaces are important types of molecular motions. Their fundamentally correct and detailed description can be obtained by solving the nuclear Schrödinger equation on a potential energy surface. Many of the chemically interesting processes involve quantum nuclear motions which are 'delocalized' over multiple potential energy wells. These 'large-amplitude' motions in addition to the high dimensionality of the vibrational problem represent challenges to the current (ro)vibrational methodology. A review of the quantum nuclear motion methodology is provided, current bottlenecks of solving the nuclear Schrödinger equation are identified, and solution strategies are reviewed. Technical details, computational results, and analysis of these results in terms of limiting models and spectroscopically relevant concepts are highlighted for selected numerical examples.

Quantum and Quantum-Classical Studies of the Photoisomerization of a Retinal Chromophore Model

We report an in-depth analysis of the photo-induced isomerization of the 2-cis-penta-2,4-dieniminium cation: a minimal model of the 11-cis retinal protonated Schiff base chromophore of the dim-light photoreceptor rhodopsin. Based on recently developed three-dimensional potentials parametrized on ab initio multi-state multi-configurational second-order perturbation theory data, we perform quantum-dynamical studies. In addition, simulations based on various quantum-classical methods, among which Tully surface hopping and the coupled-trajectory approach derived from the exact factorization, allow us to validate their performance against vibronic wavepacket propagation and, therefore, a purely quantum treatment. Quantum-dynamics results uncover qualitative differences with respect to the two-dimensional Hahn-Stock potentials, widely used as model potentials for the isomerization of the same chromophore, due to the increased dimensionality and three-mode correlation. Quantum-classical simulations show, instead, that three-dimensional model potentials are capable of capturing a number of features revealed by atomistic simulations and experimental observations. In particular, a recently reported vibrational phase relationship between double-bond torsion and hydrogen-out-of-plane modes critical for rhodopsin isomerization efficiency is correctly reproduced.

Derivation of ρ-dependent coordinate transformations for nonrigid molecules in the Hougen-Bunker-Johns formalism

In this paper, we report a series of transformations for the construction of a Hamiltonian model for nonrigid polyatomic molecules in the framework of the Hougen-Bunker-Johns formalism (HBJ). This model is expressed in normal mode coordinates for small vibrations and in a specific coordinate ρ to describe the large amplitude motion. For the first time, a general procedure linking the "true" curvilinear coordinates to ρ is proposed, allowing the expression of the potential energy part in the same coordinate representation as the kinetic energy operator, whatever the number of atoms. A Lie group-based method is also proposed for the derivation of the reference configuration in the internal axis system. This work opens new perspectives for future high-resolution spectroscopy studies of nonrigid, medium-sized molecules using HBJ-type Hamiltonians. Illustrative examples and computation of vibrational energy levels on semirigid and nonrigid molecules are given to validate this method.

Full-dimensional (12D) variational vibrational states of CH4·F-: Interplay of anharmonicity and tunneling

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.

Toward breaking the curse of dimensionality in (ro)vibrational computations of molecular systems with multiple large-amplitude motions

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 CH₄·Ar complex treated in full (12D) vibrational dimensionality.