Weakly bound molecular dimers: Intramolecular vibrational fundamentals, overtones, and tunneling splittings from full-dimensional quantum calculations using compact contracted bases of intramolecular and low-energy rigid-monomer intermolecular eigenstates

HF Trimer: A New Full-Dimensional Potential Energy Surface and Rigorous 12D Quantum Calculations of Vibrational States

HF trimer, as the smallest and the lightest cyclic hydrogen-bonded (HB) cluster, has long been a favorite prototype system for spectroscopic and theoretical investigations of the structure, energetics, spectroscopy, and dynamics of hydrogen-bond networks. Recently, rigorous quantum 12D calculations of the coupled intra- and intermolecular vibrations of this fundamental HB trimer (J. Chem. Phys. 2023, 158, 234109) were performed, employing an older ab initio-based many-body potential energy surface (PES). While the theoretical results were found to be in reasonably good agreement with the available spectroscopic data, it was also evident that it is highly desirable to develop a more accurate 12D PES of HF trimer. Motivated by this, here we report a new, and the first fully ab initio 12D PES of this paradigmatic system. Approximately 42,540 geometries were sampled and calculated at the level of CCSD(T)-F12a/AVTZ. The permutationally invariant polynomial-neural network based Δ-machine learning approach (J. Phys. Chem. Lett. 2022, 13, 4729) was employed to perform cost-efficient calculations of the basis-set-superposition error (BSSE) correction. By strategically selecting data points, this approach facilitated the construction of a high-precision PES with BSSE correction, while requiring only a minimal number of BSSE value computations. The fitting error of the final PES is only 0.035 kcal/mol. To assess its performance, the 12D fully coupled quantum calculations of excited intra- and intermolecular vibrational states of HF trimer are carried out using the rigorous methodology developed by us earlier. The results are found to be in a significantly better agreement with the available spectroscopic data than those obtained with the previously existing semiempirical 12D PES.

Revisiting Numerical Solutions of Weakly Bound Noble Gases' Vibrational Energy Levels Modeled by the Improved Lennard-Jones Potential

We revisit the numerical solutions of vibrational eigenstates of weakly bound homonuclear and heteronuclear noble gas pairs by applying a Fortran program based on the Numerov method. The harmonic, Lennard-Jones (LJ), Morse, Tang-Toennies (TT), and Improved Lennard-Jones (ILJ) potential models have been implemented to represent the potential energy curves (PECs). The obtained vibrational energies spectrum was tested on the experimental data and accurate ab initio calculations at CCSD(T)/CBS level. The vibrational eigenvalues and eigenfunctions can be reproduced accurately within the ILJ potential model. Moreover, considering the calculated lifetime of van der Waals (vdW) complexes, the implementation of ILJ rather than standard LJ potential model has a significant impact on the systems dynamics by providing more representative atomic trajectories when the function is incorporated in force fields for molecular dynamics (MD) simulations. Overall, the ILJ function is the best suited potential model for the representation of vibrational motions and the determination of vibrational energy levels of weakly bound systems, both at equilibrium and non-equilibrium conditions.

Machine-learning to predict anharmonic frequencies: a study of models and transferability

With more and more accurate electronic structure methods at hand, the inclusion of anharmonic effects in the post-processing of such data towards thermochemical properties is the next step. In this context, the description of anharmonicity has been an important topic of physical chemistry and chemical physics for a long time. In this study, anharmonic frequencies of various hydrogen-halides and halogenated hydrocarbon molecular clusters are calculated using harmonic as well as explicitly anharmonic methods, i.e., normal mode analysis and vibrational self-consistent field. Simple harmonic model based descriptors were used to predict anharmonic frequencies via multilinear regression and gradient boosting regression. Gradient boosting regression is capable of predicting reliable anharmonic data and even the simple multilinear regression model yields reasonable predictions that can account for mode-to-mode couplings. Moreover, the transferability to unseen chemical systems is assessed and it is confirmed that the machine-learned models can be applied to larger, unseen molecules.

A two-step quadrature-based variational calculation of ro-vibrational levels and wavefunctions of CO2 using a bisector-x molecule-fixed frame

In this paper, we propose a new two-step strategy for computing ro-vibrational energy levels and wavefunctions of a triatomic molecule and apply it to CO2. A two-step method [J. Tennyson and B. T. Sutcliffe, Mol. Phys., 1986, 58, 1067] uses a basis whose functions are products of K-dependent "vibrational" functions and symmetric top functions. K is the quantum number for the molecule-fixed z component of the angular momentum. For a linear molecule, a two-step method is efficient because the Hamiltonian used to compute the basis functions includes the largest coupling term. The most important distinguishing feature of the two-step method we propose is that it uses an associated Legendre basis and quadrature rather than a K-dependent discrete variable representation. This reduces the cost of the calculation and simplifies the method. We have computed ro-vibrational energy levels with J up to 100 for CO2, on an accurate available potential energy surface which is known as the AMES-2 PES and present a subset of those levels. We have converged most levels up to 20 000 cm-1 to 0.0001 cm-1.

HCl trimer: HCl-stretch excited intramolecular and intermolecular vibrational states from 12D fully coupled quantum calculations employing contracted intra- and inter-molecular bases

We present fully coupled, full-dimensional quantum calculations of the inter- and intra-molecular vibrational states of HCl trimer, a paradigmatic hydrogen-bonded molecular trimer. They are performed utilizing the recently developed methodology for the rigorous 12D quantum treatment of the vibrations of the noncovalently bound trimers of flexible diatomic molecules [Felker and Bačić, J. Chem. Phys. 158, 234109 (2023)], which was previously applied to the HF trimer by us. In this work, the many-body 12D potential energy surface (PES) of (HCl)3 [Mancini and Bowman, J. Phys. Chem. A 118, 7367 (2014)] is employed. The calculations extend to the intramolecular HCl-stretch excited vibrational states of the trimer with one- and two-quanta, together with the low-energy intermolecular vibrational states in the two excited v = 1 intramolecular vibrational manifolds. They reveal significant coupling between the intra- and inter-molecular vibrational modes. The 12D calculations also show that the frequencies of the v = 1 HCl stretching states of the HCl trimer are significantly redshifted relative to those of the isolated HCl monomer. Detailed comparison is made between the results of the 12D calculations on the two-body PES, obtained by removing the three-body term from the original 2 + 3-body PES, and those computed on the 2 + 3-body PES. It demonstrates that the three-body interactions have a strong effect on the trimer binding energy as well as on its intra- and inter-molecular vibrational energy levels. Comparison with the available spectroscopic data shows that good agreement with the experiment is achieved only if the three-body interactions are included. Some low-energy vibrational states localized in a secondary minimum of the PES are characterized as well.

H2O-HCN complex: A new potential energy surface and intermolecular rovibrational states from rigorous quantum calculations

In this work the H2O-HCN complex is quantitatively characterized in two ways. First, we report a new rigid-monomer 5D intermolecular potential energy surface (PES) for this complex, calculated using the symmetry-adapted perturbation theory based on density functional theory method. The PES is based on 2833 ab initio points computed employing the aug-cc-pVQZ basis set, utilizing the autoPES code, which provides a site-site analytical fit with the long-range region given by perturbation theory. Next, we present the results of the quantum 5D calculations of the fully coupled intermolecular rovibrational states of the H2O-HCN complex for the total angular momentum J values of 0, 1, and 2, performed on the new PES. These calculations rely on the quantum bound-state methodology developed by us recently and applied to a variety of noncovalently bound binary molecular complexes. The vibrationally averaged ground-state geometry of H2O-HCN determined from the quantum 5D calculations agrees very well with that from the microwave spectroscopic measurements. In addition, the computed ground-state rotational transition frequencies, as well as the B and C rotational constants calculated for the ground state of the complex, are in excellent agreement with the experimental values. The assignment of the calculated intermolecular vibrational states of the H2O-HCN complex is surprisingly challenging. It turns out that only the excitations of the intermolecular stretch mode can be assigned with confidence. The coupling among the angular degrees of freedom (DOFs) of the complex is unusually strong, and as a result most of the excited intermolecular states are unassigned. On the other hand, the coupling of the radial, intermolecular stretch mode and the angular DOFs is weak, allowing straightforward assignment of the excitation of the former.

Computing vibrational spectra using a new collocation method with a pruned basis and more points than basis functions: Avoiding quadrature

We present a new collocation method for computing the vibrational spectrum of a polyatomic molecule. Some form of quadrature or collocation is necessary when the potential energy surface does not have a simple form that simplifies the calculation of the potential matrix elements required to do a variational calculation. With quadrature, better accuracy is obtained by using more points than basis functions. To achieve the same advantage with collocation, we introduce a collocation method with more points than basis functions. Critically important, the method can be used with a large basis because it is incorporated into an iterative eigensolver. Previous collocation methods with more points than functions were incompatible with iterative eigensolvers. We test the new ideas by computing energy levels of molecules with as many as six atoms. We use pruned bases but expect the new method to be advantageous whenever one uses a basis for which it is not possible to find an accurate quadrature with about as many points as there are basis functions. For our test molecules, accurate energy levels are obtained even using non-optimal, simple, equally spaced points.

Computing excited OH stretch states of water dimer in 12D using contracted intermolecular and intramolecular basis functions

Due to the ubiquity and importance of water, water dimer has been intensively studied. Computing the (ro-)vibrational spectrum of water dimer is challenging. The potential has eight wells separated by low barriers, which makes harmonic approximations of limited utility. A variational approach is imperative, but difficult because there are 12 coupled vibrational coordinates. In this paper, we use a product contracted basis whose functions are products of intramolecular and intermolecular functions computed using an iterative eigensolver. An intermediate matrix F facilitates calculating matrix elements. Using F, it is possible to do calculations on a general potential without storing the potential on the full quadrature grid. We find that surprisingly many intermolecular functions are required. This is due to the importance of coupling between inter- and intra-molecular coordinates. The full G₁₆ symmetry of water dimer is exploited. We calculate, for the first time, monomer excited stretch states and compare P(1) transition frequencies with their experimental counterparts. We also compare with experimental vibrational shifts and tunneling splittings. Surprisingly, we find that the largest tunneling splitting, which does not involve the interchange of the two monomers, is smaller in the asymmetric stretch excited state than in the ground state. Differences between levels we compute and those obtained with a [6+6]D adiabatic approximation [Leforestier et al. J. Chem. Phys. 137 014305 (2012)] are ∼0.6 cm^-1 for states without monomer excitation, ∼4 cm^-1 for monomer excited bend states, and as large as ∼10 cm^-1 for monomer excited stretch states.

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.

A computational inspection of the dissociation energy of mid-sized organic dimers

The gas-phase value of the dissociation energy ( D0) is a key parameter employed in both experimental and theoretical descriptions of noncovalent complexes. The D0 data were obtained for a set of mid-sized organic dimers in their global minima which was located using geometry optimizations that applied ample basis sets together with either the conventional second-order Møller–Plesset (MP2) method or several dispersion-corrected density-functional theory (DFT-D) schemes. The harmonic vibrational zero-point (VZP) and deformation energies from the MP2 calculations were combined with electronic energies from the coupled cluster theory with singles, doubles, and iterative triples [CCSD(T)] extrapolated to the complete basis set (CBS) limit to estimate D0 with the aim of inspecting values that were most recently measured, and an analogous comparison was performed using the DFT-D data. In at least one case (namely, for the aniline⋯methane cluster), the D0 estimate that employed the CCSD(T)/CBS energies differed from experiment in the way that could not be explained by a possible deficiency in the VZP contribution. Curiously, one of the DFT-D schemes (namely, the B3LYP-D3/def2-QZVPPD) was able to reproduce all measured D0 values to within 1.0 kJ/mol from experimental error bars. These findings show the need for further measurements and computations of some of the complexes. In order to facilitate such studies, the physical nature of intermolecular interactions in the investigated dimers was analyzed by means of the DFT-based symmetry-adapted perturbation theory.

Vibrational Tunneling Spectra of Molecules with Asymmetric Wells: A Combined Vibrational Configuration Interaction and Instanton Approach

A combined approach that uses the vibrational configuration interaction (VCI) and semiclassical instanton theory was developed to study vibrational tunneling spectra of molecules with multiple wells in full dimensionality. The method can be applied to calculate low-lying vibrational states in the systems with an arbitrary number of minima, which are not necessarily equal in energy or shape. It was tested on a two-dimensional double-well model system and on malonaldehyde, and the calculations reproduced the exact quantum mechanical (QM) results with high accuracy. The method was subsequently applied to calculate the vibrational spectrum of the asymmetrically deuterated malonaldehyde with nondegenerate vibrational frequencies in the two wells. The spectrum is obtained at a cost of single-well VCI calculations used to calculate the local energies. The interactions between states of different wells are computed semiclassically using the instanton theory at a comparatively negligible computational cost. The method is particularly suited to systems in which the wells are separated by large potential barriers and tunneling splittings are small, for example, in some water clusters, when the exact QM methods come at a prohibitive computational cost.

Highly accurate HF dimer ab initio potential energy surface

A highly accurate, (HF)₂ potential energy surface (PES) is constructed based on ab initio calculations performed at the coupled-cluster single double triple level of theory with an aug-cc-pVQZ-F12 basis set at about 152 000 points. A higher correlation correction is computed at coupled-cluster single double triple quadruple level for 2000 points and is considered alongside other more minor corrections due to relativity, core-valence correlation, and Born-Oppenheimer failure. The analytical surface constructed uses 500 constants to reproduce the ab initio points with a standard deviation of 0.3 cm^-1. Vibration-rotation-inversion energy levels of the HF dimer are computed for this PES by variational solution of the nuclear-motion Schrödinger equation using the program WAVR4. Calculations over an extended range of rotationally excited states show very good agreement with the experimental data. In particular, the known empirical rotational constants B for the ground vibrational states are predicted to better than about 2 MHz. B constants for excited vibrational states are reproduced several times more accurately than by previous calculations. This level of accuracy is shown to extend to higher excited inter-molecular vibrational states v and higher excited rotational quantum numbers (J, K_a).

Reduced-dimensional vibrational models of the water dimer

A model based on the finite-basis representation of a vibrational Hamiltonian expressed in internal coordinates is developed. The model relies on a many-mode, low-order expansion of both the kinetic energy operator and the potential energy surface (PES). Polyad truncations and energy ceilings are used to control the size of the vibrational basis to facilitate accurate computations of the OH stretch and HOH bend intramolecular transitions of the water dimer (H₂ ¹⁶O)₂. Advantages and potential pitfalls of the applied approximations are highlighted. The importance of choices related to the treatment of the kinetic energy operator in reduced-dimensional calculations and the accuracy of different water dimer PESs are discussed. A range of different reduced-dimensional computations are performed to investigate the wavenumber shifts in the intramolecular transitions caused by the coupling between the intra- and intermolecular modes. With the use of symmetry, full 12-dimensional vibrational energy levels of the water dimer are calculated, predicting accurately the experimentally observed intramolecular fundamentals. It is found that one can also predict accurate intramolecular transition wavenumbers for the water dimer by combining a set of computationally inexpensive reduced-dimensional calculations, thereby guiding future effective-Hamiltonian treatments.

Noncovalently bound molecular complexes beyond diatom–diatom systems: full-dimensional, fully coupled quantum calculations of rovibrational states

The methodological advances made in recent years have significantly extended the range and dimensionality of noncovalently bound molecular complexes for which full-dimensional quantum calculations of their rovibrational states are feasible.

DCl-H2O, HCl-D2O, and DCl-D2O Dimers: Inter- and Intramolecular Vibrational States and Frequency Shifts from Fully Coupled Quantum Calculations on a Full-Dimensional Neural Network Potential Energy Surface

We report full-dimensional and fully coupled quantum calculations of the inter- and intramolecular vibrational states of three isotopologues of the hydrogen chloride-water dimer: DCl-H₂O (DH), HCl-D₂O (HD), and DCl-D₂O (DD). The present study extends our recent theoretical investigation of the nine-dimensional (9D) vibrational level structure of the HCl-H₂O (HH) dimer [Liu, Y.; Li, J.; Felker, P. M.; Bačić, Z. Phys. Chem. Chem. Phys. 2021, 23, 7101-7114]. It employs the same accurate 9D permutation invariant polynomial-neural network potential energy surface and the highly efficient bound-state methodology. The objective of this work is to elucidate the isotopologue variations of a range of bound-state properties of the hydrogen chloride-water dimer and compare them to those of the HH dimer. In order to achieve this, for the isotopologues considered, the rigorous 9D quantum calculations performed encompass all intramolecular vibrational fundamentals, and their frequency shifts relative to the isolated monomer values, together with the low-lying intermolecular vibrational states in each of the intramolecular vibrational manifolds of interest. Moreover, for the ground state of each isotopologue, several informative vibrationally averaged intermolecular geometric properties of the dimer are computed, as well as the three rotational constants. The energies of the intermolecular inversion and rock modes, which mainly involve the motions of the water moiety, differ greatly for H₂O and D₂O, but are much less sensitive to whether the hydrogen chloride isotopologue is HCl or DCl. On the other hand, the excitation of the HCl/DCl stretch changes significantly the energies of the water inversion and rock modes. The DCl stretch frequency shift computed in 9D for the DD dimer, -114.91 cm^-1, agrees extremely well with the corresponding experimental value of -115.20 cm^-1 measured by Saykally and co-workers.

Fingerprint region of the formic acid dimer: variational vibrational computations in curvilinear coordinates

Curvilinear kinetic energy models are developed for variational nuclear motion computations including the inter- and the low-frequency intra-molecular degrees of freedom of the formic acid dimer. The coupling of the inter- and intra-molecular modes is studied by solving the vibrational Schrödinger equation for a series of vibrational models, from two up to ten active vibrational degrees of freedom by selecting various combinations of active modes and constrained coordinate values. Vibrational states, nodal assignment, and infrared vibrational intensity information is computed using the full-dimensional potential energy surface (PES) and electric dipole moment surface developed by Qu and Bowman [Phys. Chem. Chem. Phys., 2016, 18, 24835; J. Chem. Phys., 2018, 148, 241713]. Good results are obtained for several fundamental and combination bands in comparison with jet-cooled vibrational spectroscopy experiments, but the description of the ν₈ and ν₉ fundamental vibrations, which are close in energy and have the same symmetry, appears to be problematic. For further progress in comparison with experiment, the potential energy surface, and in particular, its multi-dimensional couplings representation, requires further improvement.

HCl–H2O dimer: an accurate full-dimensional potential energy surface and fully coupled quantum calculations of intra- and intermolecular vibrational states and frequency shifts

The present work reports a new full-dimensional potential energy surface (PES) of the HCl–H₂O dimer, and the first fully coupled 9D quantum calculations of the intra- and intermolecular vibrational states of the complex, utilizing this PES.

Tunneling splittings of vibrationally excited states using general instanton paths

A multidimensional semiclassical method for calculating tunneling splittings in vibrationally excited states of molecules using Cartesian coordinates is developed. It is an extension of the theory by Mil'nikov and Nakamura [J. Chem. Phys. 122, 124311 (2005)] to asymmetric paths that are necessary for calculating tunneling splitting patterns in multi-well systems, such as water clusters. Additionally, new terms are introduced in the description of the semiclassical wavefunction that drastically improves the splitting estimates for certain systems. The method is based on the instanton theory and builds the semiclassical wavefunction of the vibrationally excited states from the ground-state instanton wavefunction along the minimum action path and its harmonic neighborhood. The splittings of excited states are thus obtained at a negligible added numerical effort. The cost is concentrated, as for the ground-state splittings, in the instanton path optimization and the hessian evaluation along the path. The method can thus be applied without modification to many mid-sized molecules in full dimensionality and in combination with on-the-fly evaluation of electronic potentials. The tests were performed on several model potentials and on the water dimer.

A variational calculation of vibrational levels of vinyl radical

We report the vibrational energy levels of vinyl radical (VR) that are computed with a Lanczos eigensolver and a contracted basis. Many of the levels of the two previous VR variational calculations differ significantly and differ also from those reported in this paper. We identify the source of and correct symmetry errors on the potential energy surfaces used in the previous calculations. VR has two equivalent equilibrium structures. By plotting wavefunction cuts, we show that two tunneling paths play an important role. Using the computed wavefunctions, it is possible to assign many states and thereby to determine tunneling splittings that are compared with their experimental counterparts. Our computed red shift of the hot band at 2897.23 cm^-1, observed by Dong et al. [J Chem. Phys. 128, 044305 (2008)], is 4.47 cm^-1, which is close to the experimental value of 4.63 cm^-1.

Exact quantum dynamics background of dispersion interactions: case study for CH4·Ar in full (12) dimensions

A full-dimensional ab initio potential energy surface of spectroscopic quality is developed for the van-der-Waals complex of a methane molecule and an argon atom. Variational vibrational states are computed on this surface including all twelve (12) vibrational degrees of freedom of the methane-argon complex using the GENIUSH computer program and the Smolyak sparse grid method. The full-dimensional computations make it possible to study the fine details of the interaction and distortion effects and to make a direct assessment of the reduced-dimensionality models often used in the quantum dynamics study of weakly-bound complexes. A 12-dimensional (12D) vibrational computation including only a single harmonic oscillator basis function (9D) to describe the methane fragment (for which we use the ground-state effective structure as the reference structure) has a 0.40 cm^-1 root-mean-square error (rms) with respect to the converged 12D bound-state excitation energies, which is less than half of the rms of the 3D model set up with the 〈r〉₀ methane structure. Allowing 10 basis functions for the methane fragment in a 12D computation performs much better than the 3D models by reducing the rms of the bound state vibrational energies to 0.07 cm^-1. The full-dimensional potential energy surface correctly describes the dissociation of the system, which together with further development of the variational (ro)vibrational methodology opens a route to the study of the role of dispersion forces in the excited methane vibrations and the energy transfer from the intra- to the intermolecular vibrational modes.