The hierarchical expansion of the kinetic energy operator in curvilinear coordinates extended to the vibrational configuration interaction method

VSCF/VCI theory based on the Podolsky Hamiltonian

While the vibrational spectra of semi-rigid molecules can be computed on approaches relying on the Watson Hamiltonian, floppy molecules or molecular clusters are better described by Hamiltonians, which are capable of dealing with any curvilinear coordinates. It is the kinetic energy operator (KEO) of these Hamiltonians, which render the correlated calculations relying on them rather costly. Novel implementation of vibrational self-consistent field theory and vibrational configuration interaction theory on the basis of the Podolsky Hamiltonian are reported, in which the inverse of the metric tensor, i.e., the G matrix, is represented by an n-mode expansion expressed in terms of polynomials. An analysis of the importance of the individual terms of the KEO with respect to the truncation orders of the n-mode expansion is provided. Benchmark calculations have been performed for the cis-HOPO and methanimine, H2CNH, molecules and are compared to experimental data and to calculations based on the Watson Hamiltonian and the internal coordinate path Hamiltonian.

Theoretical Insights into the Vibrational Structure of Carbon Dioxide Rare-Gas Complexes

Two new flexible-monomer two-body ab initio potential energy surfaces (PESs) for the neon and krypton van der Waals complexes with carbon dioxide were developed, extending our previous work on the Ar-CO2 molecule. The accuracy of the PESs was validated by their agreement with the vibrational spectrum of the rare-gas complexes. The intermolecular and intramolecular vibrational excitation energies were computed at the vibrational self-consistent field and vibrational configuration interaction levels of theory. Overall, the agreement between theory and experiment is excellent throughout the vibrational spectra. The observed slight splitting of the bending modes, resulting from their nondegeneracy in the complexes, is confirmed by our computations, and the results qualitatively agree with the experiment. The splitting increases with increasing polarizability of the rare-gas atom. Additionally, we explain a discrepancy in the mode assignment in the intermolecular region of the neon complex with our VCI character assignment.

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.

Determining internal coordinate sets for optimal representation of molecular vibration

Arising from the harmonic approximation in solving the vibrational Schrödinger equation, normal modes dissect molecular vibrations into distinct degrees of freedom. Normal modes are widely used as they give rise to descriptive vibrational notations and are convenient for expanding anharmonic potential energy surfaces as an alternative to higher-order Taylor series representations. Usually, normal modes are expressed in Cartesian coordinates, which bears drawbacks that can be overcome by switching to internal coordinates. Considering vibrational notations, normal modes with delocalized characters are difficult to denote, but internal coordinates offer a route to clearer notations. Based on the Hessian, normal mode decomposition schemes for a given set of internal coordinates can describe a normal mode by its contributions from internal coordinates. However, choosing a set of internal coordinates is not straightforward. While the Hessian provides unique sets of normal modes, various internal coordinate sets are possible for a given system. In the present work, we employ a normal mode decomposition scheme to choose an optimal set. Therefore, we screen reasonable sets based on topology and symmetry considerations and rely on a metric that minimizes coupling between internal coordinates. Ultimately, the Nomodeco toolkit presented here generates internal coordinate sets to find an optimal set for representing molecular vibrations. The resulting contribution tables can be used to clarify vibrational notations. We test our scheme on small to mid-sized molecules, showing how the space of definable internal coordinate sets can significantly be reduced.

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.

Exploring the anharmonic vibrational structure of carbon dioxide trimers

Our previously developed mbCO2 potential [O. Sode and J. N. Cherry, J. Comput. Chem. 38, 2763 (2017)] is used to describe the vibrational structure of the intermolecular motions of the CO₂ trimers: barrel-shaped and cyclic trimers. Anharmonic corrections are accounted for using the vibrational self-consistent field theory, vibrational second-order Møller-Plesset perturbation (VMP2) theory, and vibrational configuration interaction (VCI) methods and compared with experimental observations. For the cyclic structure, we revise the assignments of two previously observed experimental peaks based on our VCI and VMP2 results. We note that the experimental band observed near 13 cm^-1 is the out-of-phase out-of-plane degenerate motion with E″ symmetry, while the peak observed at 18 cm^-1 likely corresponds to the symmetric out-of-plane torsion A″ vibration. Since the VCI treatment of the vibrational motions accounts for vibrational mixing and delocalization, overtones and combination bands were also observed and quantified in the intermolecular regions of the two trimer isomers.

Structural sensitivity of CH vibrational band in methyl benzoate

The CH vibrational bands of methyl benzoate are studied to understand its coupling pattern with other vibrational bands of the biological molecule. This will facilitate to understand the biological structure and dynamics in spectroscopic as well as in microscopic study. Due to the congested spectroscopic pattern, near degeneracy, and strong anharmonicity of the CH stretch vibrations, assignment of the CH vibrational frequencies are often misleading. Anharmonic vibrational frequency calculation with multidimensional potential energy surface interprets the CH vibrational spectra more accurately. In this article we have presented the importance of multidimensional potential energy surface in anharmonic vibrational frequency calculation and discuss the unexpected red shift of asymmetric CH stretch vibration of methyl group. The CD stretch vibrational band which is splitted to double peaks due to the Fermi resonance is also discussed here.

Pushing the limits in accurate vibrational structure calculations: anharmonic frequencies of lithium fluoride clusters (LiF)n, n = 2-10

The vibrational spectra of a series of small lithium fluoride clusters, i.e. (LiF)n, n = 2-10, were studied by vibrational configuration interaction (VCI) calculations relying on potential energy surfaces including three-mode coupling terms and being obtained from explicitly correlated local coupled cluster calculations. Due to the account for anharmonicity effects, the simulated spectra allow for a direct comparison with experimental data and may thus help to identify clusters in the experiments. Even structurally closely related clusters can clearly be distinguished by infrared spectroscopy. The largest system in this study required more than 1000 basis functions in the electronic structure calculations and more than 10(7) configurations in the vibrational structure calculations and became computationally feasible only due to a combination of different approximations and highly parallelized algorithms.

Efficient anharmonic vibrational spectroscopy for large molecules using local-mode coordinates

This article presents a general computational approach for efficient simulations of anharmonic vibrational spectra in chemical systems. An automated local-mode vibrational approach is presented, which borrows techniques from localized molecular orbitals in electronic structure theory. This approach generates spatially localized vibrational modes, in contrast to the delocalization exhibited by canonical normal modes. The method is rigorously tested across a series of chemical systems, ranging from small molecules to large water clusters and a protonated dipeptide. It is interfaced with exact, grid-based approaches, as well as vibrational self-consistent field methods. Most significantly, this new set of reference coordinates exhibits a well-behaved spatial decay of mode couplings, which allows for a systematic, a priori truncation of mode couplings and increased computational efficiency. Convergence can typically be reached by including modes within only about 4 Å. The local nature of this truncation suggests particular promise for the ab initio simulation of anharmonic vibrational motion in large systems, where connection to experimental spectra is currently most challenging.

Vibrational spectroscopy of Methyl benzoate

Methyl benzoate is studied as a model compound for the development of new IR pulse schemes with possible applicability to biomolecules.

Numerically constructed internal-coordinate Hamiltonian with Eckart embedding and its application for the inversion tunneling of ammonia

It is shown that the use of an Eckart-frame embedding with a kinetic energy operator expressed in curvilinear internal coordinates becomes feasible and straightforward to implement for arbitrary molecular compositions and internal coordinates if the operator is defined numerically over a (discrete variable representation) grid. The algorithm proposed utilizes the transformation method of Dymarsky and Kudin to maintain the rotational Eckart condition. In order to demonstrate the applicability and flexibility of our approach the non-rigid ammonia molecule is considered and the corresponding rotational-vibrational energy levels and wave functions are computed using kinetic energy operators with three different embeddings. Two of them fulfill the Eckart conditions corresponding to a trigonal pyramidal (C3v) and a trigonal planar (D3h) reference structure and the third one is a non-Eckart frame. The computed energy levels are, of course, identical, and the structure of the three different wave-function representations are analyzed in terms of the rigid rotor functions for a symmetric top. The possible advantages of one frame representation over another are discussed concerning the interpretation of the rovibrational states in terms of the traditional rigid rotor labels.