Hierarchical expansion of the kinetic energy operator in curvilinear coordinates for the vibrational self-consistent field 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.

Two-dimensional Infrared Spectroscopy Reveals Better Insights of Structure and Dynamics of Protein

Proteins play an important role in biological and biochemical processes taking place in the living system. To uncover these fundamental processes of the living system, it is an absolutely necessary task to understand the structure and dynamics of the protein. Vibrational spectroscopy is an established tool to explore protein structure and dynamics. In particular, two-dimensional infrared (2DIR) spectroscopy has already proven its versatility to explore the protein structure and its ultrafast dynamics, and it has essentially unprecedented time resolutions to observe the vibrational dynamics of the protein. Providing several examples from our theoretical and experimental efforts, it is established here that two-dimensional vibrational spectroscopy provides exceptionally more information than one-dimensional vibrational spectroscopy. The structural information of the protein is encoded in the position, shape, and strength of the peak in 2DIR spectra. The time evolution of the 2DIR spectra allows for the visualisation of molecular motions.

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.

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.

Largely reduced grid densities in a vibrational self-consistent field treatment do not significantly impact the resultingwavenumbers

Especially for larger molecules relevant to life sciences, vibrational self-consistent field (VSCF) calculations can become unmanageably demanding even when only first and second order potential coupling terms are considered. This paper investigates to what extent the grid density of the VSCF’s underlying potential energy surface can be reduced without sacrificing accuracy of the resulting wavenumbers. Including single-mode and pair contributions, a reduction to eight points per mode did not introduce a significant deviation but improved the computational efficiency by a factor of four. A mean unsigned deviation of 1.3% from the experiment could be maintained for the fifteen molecules under investigation and the approach was found to be applicable to rigid, semi-rigid and soft vibrational problems likewise. Deprotonated phosphoserine, stabilized by two intramolecular hydrogen bonds, was investigated as an exemplary application.

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.