Search

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

We present the computational methodology, which for the first time allows rigorous twelve-dimensional (12D) quantum calculations of the coupled intramolecular and intermolecular vibrational states of hydrogen-bonded trimers of flexible diatomic molecules. Its starting point is the approach that we introduced recently for fully coupled 9D quantum calculations of the intermolecular vibrational states of noncovalently bound trimers comprised of diatomics treated as rigid. In this paper, it is extended to include the intramolecular stretching coordinates of the three diatomic monomers. The cornerstone of our 12D methodology is the partitioning of the full vibrational Hamiltonian of the trimer into two reduced-dimension Hamiltonians, one in 9D for the intermolecular degrees of freedom (DOFs) and another in 3D for the intramolecular vibrations of the trimer, and a remainder term. These two Hamiltonians are diagonalized separately, and a fraction of their respective 9D and 3D eigenstates is included in the 12D product contracted basis for both the intra- and intermolecular DOFs, in which the matrix of the full 12D vibrational Hamiltonian of the trimer is diagonalized. This methodology is implemented in the 12D quantum calculations of the coupled intra- and intermolecular vibrational states of the hydrogen-bonded HF trimer on an ab initio calculated potential energy surface (PES). The calculations encompass the one- and two-quanta intramolecular HF-stretch excited vibrational states of the trimer and low-energy intermolecular vibrational states in the intramolecular vibrational manifolds of interest. They reveal several interesting manifestations of significant coupling between the intra- and intermolecular vibrational modes of (HF)3. The 12D calculations also show that the frequencies of the v = 1, 2 HF stretching states of the HF trimer are strongly redshifted in comparison to those of the isolated HF monomer. Moreover, the magnitudes of these trimer redshifts are much larger than that of the redshift for the stretching fundamental of the donor-HF moiety in (HF)2, most likely due to the cooperative hydrogen bonding in (HF)3. The agreement between the 12D results and the limited spectroscopic data for the HF trimer, while satisfactory, leaves room for improvement and points to the need for a more accurate PES. [ABSTRACT FROM AUTHOR]

We present the computational methodology that allows rigorous and efficient nine-dimensional (9D) quantum calculations of the intermolecular vibrational states of noncovalently bound trimers of diatomic molecules, with the monomers treated as rigid. The full 9D vibrational Hamiltonian of the trimer is partitioned into a 3D "frame" (or stretching) Hamiltonian and a 6D "bend" Hamiltonian. These two Hamiltonians are diagonalized separately, and a certain number of their lowest-energy eigenstates is included in the final 9D product contracted basis in which the full 9D intermolecular vibrational Hamiltonian is diagonalized. This methodology is applied to the 9D calculations of the intermolecular vibrational levels of (HF)3, a prototypical hydrogen-bonded trimer, on the rigid-monomer version of an ab initio calculated potential energy surface (PES). They are the first to include fully the stretch-bend coupling present in the trimer. The frequencies of all bending fundamentals considered from the present 9D calculations are about 10% lower than those from the earlier quantum 6D calculations that considered only the bending modes of the HF trimer. This means that the stretch-bend coupling is strong, and it is imperative to include it in any accurate treatment of the (HF)3 vibrations aiming to assess the accuracy of the PES employed. Moreover, the 9D results are in better agreement with the limited available spectroscopic data that those from the 6D calculations. In addition, the 9D results show sensitivity to the value of the HF bond length, equilibrium or vibrationally averaged, used in the calculations. The implication is that full-dimensional 12D quantum calculations will be required to obtain definitive vibrational excitation energies for a given PES. Our study also demonstrates that the nonadditive three-body interactions are very significant in (HF)3 and have to be included in order to obtain accurate intermolecular vibrational energy levels of the trimer. [ABSTRACT FROM AUTHOR]

We present quantum five-dimensional bound-state calculations of the fully coupled intermolecular rovibrational states of H2O–CO2 and D2O–CO2 van der Waals (vdW) complexes in the rigid-monomer approximation for the total angular momentum J values of 0, 1, and 2. A rigid-monomer version of the recent ab initio full-dimensional (12D) potential energy surface of H2O–CO2 [Q. Wang and J. M. Bowman, J. Chem. Phys. 147, 161714 (2017)] is employed. This treatment provides for the first time a rigorous and comprehensive description of the intermolecular rovibrational level structure of the two isotopologues that includes the internal-rotation tunneling splittings and their considerable sensitivity to rotational and intermolecular vibrational excitations, as well as the rotational constants of the two vdW complexes. Two approaches are used in the calculations, which differ in the definition of the dimer-fixed (DF) frame and the coordinates associated with them. We demonstrate that with the approach introduced in this work, where the DF frame is fixed to the CO2 moiety, highly accurate results are obtained using significantly smaller basis sets in comparison to those for the alternative approach. In addition, the resulting wavefunctions tend to lend themselves better to physical interpretation and assignment. The H2O–CO2 ground-state internal-rotation tunneling splittings, the rotational transition frequencies, and the rotational constants of both vdW complexes are in excellent agreement with the experimental results. The calculated intermolecular vibrational fundamentals agree well with the scant terahertz spectroscopy data for these complexes in cryogenic neon matrices. [ABSTRACT FROM AUTHOR]

We present efficient yet rigorous, full-dimensional quantum bound-state calculations of the fully coupled J = 0 and one intra- and intermolecular rovibrational levels of H2O–CO and D2O–CO complexes. The new ab initio nine-dimensional (9D) potential energy surface (PES) [Y. Liu and J. Li, Phys. Chem. Chem. Phys. 21, 24101 (2019)] is employed. In the spirit of the recently introduced general procedure [P. M. Felker and Z. Bačić J. Chem. Phys. ´ 151, 024305 (2019)], the 9D rovibrational Hamiltonian is partitioned into a 5D (rigid-monomer) intermolecular Hamiltonian, two intramolecular vibrational Hamiltonians—one for the water monomer (3D) and another for the CO monomer (1D), and a 9D remainder term. The low-energy eigenstates of the three reduced-dimension Hamiltonians are used to build up the 9D product contracted basis, in which the matrix of the full rovibrational Hamiltonian is diagonalized. In line with the findings of our earlier study referenced above, the 5D intermolecular eigenstates included in the 9D bases extend up to at most 230 cm−1 above the lowestenergy state of the given parity, much less than the intramolecular fundamentals of the two complexes that span the range of energies from about 1200 cm−1 to 3800 cm−1 . The resulting Hamiltonian matrices are small for the 9D quantum problem considered, ≈ 10 000 for J = 0 and 13 500 for J = 1 calculations, allowing for direct diagonalization. The 9D calculations permit exploring a number of features of the rovibrational level structure of H2O–CO and D2O–CO that are beyond the quantum 5D rigid-monomer treatments reported to date. These include the differences in the magnitudes of the hydrogen-exchange tunneling splittings computed in 9D and 5D, the sensitivity of the tunneling splittings to the intramolecular vibrational excitation, the frequency shifts of the intramolecular vibrational modes, which, depending on the mode, can be either blue- or redshifts, and the effects of the excitation of the intramolecular fundamentals on the low-lying intermolecular eigenstates. Also examined is the extent of the eigenstate delocalization over the two minima on the PES. Whenever possible, a comparison is made with the experimental data in the literature. [ABSTRACT FROM AUTHOR]

We present a method for efficient calculation of intramolecular vibrational excitations of H2O inside C60, together with the low-energy intermolecular translation-rotation states within each intramolecular vibrational manifold. Apart from assuming rigid C60, this nine-dimensional (9D) quantum treatment is fully coupled. Following the recently introduced approach [P. M. Felker and Z. Bačić, J. Chem. Phys. 151, 024305 (2019)], the full 9D vibrational Hamiltonian of H2O@C60 is partitioned into two reduced-dimension Hamiltonians, a 6D one for the intermolecular vibrations and another in 3D for the intramolecular degrees of freedom, and a 9D remainder term. The two reduced-dimension Hamiltonians are diagonalized, and their eigenvectors are used to build up a product contracted basis in which the full vibrational Hamiltonian is diagonalized. The efficiency of this methodology derives from the insight of our earlier study referenced above that converged high-energy intramolecular vibrational excitations of weakly bound molecular complexes can be obtained from fully coupled quantum calculations where the full-dimensional product contracted basis includes only a small number of intermolecular vibrational eigenstates spanning the range of energies much below those of the intramolecular vibrational states of interest. In this study, the eigenstates included in the 6D intermolecular contacted basis extend to only 410 cm−1 above the ground state, which is much less than the H2O stretch and bend fundamentals, at ≈3700 and ≈1600 cm−1, respectively. The 9D calculations predict that the fundamentals of all three intramolecular modes, as well as the bend overtone, of the caged H2O are blueshifted relative to those of the gas-phase H2O, the two stretch modes much more so than the bend. Excitation of the bend mode affects the energies of the low-lying H2O rotational states significantly more than exciting either of the stretching modes. The center-of-mass translational fundamental is virtually unaffected by the excitation of any of the intramolecular vibrational modes. Further progress hinges on the experimental measurement of the vibrational frequency shifts in H2O@C60 and ab initio calculation of a high-quality 9D potential energy surface for this endohedral complex, neither of which is presently available. [ABSTRACT FROM AUTHOR]

We report the results of calculations pertaining to the HH intramolecular stretching fundamentals of (p - H 2 ) 2 encapsulated in the large cage of structure II clathrate hydrate. The eight-dimensional (8D) quantum treatment assumes rotationless (j = 0) H2 moieties and a rigid clathrate structure but is otherwise fully coupled. The ( H 2 ) 2 -clathrate interaction is constructed in a pairwise-additive fashion, by combining the ab initio H2–H2O pair potential for flexible H2 and rigid H2O [D. Lauvergnat et al., J. Chem. Phys. 150, 154303 (2019)] and the six-dimensional (6D) H2–H2 potential energy surface [R. J. Hinde, J. Chem. Phys. 128, 154308 (2008)]. The calculations are performed by first solving for the eigenstates of a reduced-dimension 6D "intermolecular" Hamiltonian extracted from the full 8D Hamiltonian by taking the H2 moieties to be rigid. An 8D contracted product basis for the solution of the full problem is then constructed from a small number of the lowest-energy 6D intermolecular eigenstates and two discrete variable representations covering the H2-monomer internuclear distances. Converged results are obtained already by including just the two lowest intermolecular eigenstates in the final 8D basis of dimension 128. The two HH vibrational stretching fundamentals are computed for three hydrate domains having an increasing number of H2O molecules. For the largest domain, the two fundamentals are found to be site-split by ∼0.5 cm−1 and to be redshifted by about 24 cm−1 from the free-H2 monomer stretch frequency, in excellent agreement with the experimental value of 26 cm−1. A first-order perturbation theory treatment gives results that are nearly identical to those of the 8D quantum calculations. [ABSTRACT FROM AUTHOR]

We report full-dimensional and fully coupled quantum bound-state calculations of the J=1 intra- and intermolecular rovibrational states of two isotopologues of the hydrogen chloride-water dimer, HCl-H2O (HH) and DCl-H2O (DH). The present study complements our recent theoretical investigations of the J=0 nine-dimensional (9D) vibrational level structure of these and two other H/D isotopologues of this noncovalently bound molecular complex, and employs the same accurate 9D permutation invariant polynomial-neural network potential energy surface. The calculations yield all intramolecular vibrational fundamentals of the HH and DH dimers and the low-energy intermolecular rovibrational states in these intramolecular vibrational manifolds. The results are compared with those of the 9D J=0 calculations of the same dimers. The energy differences between the K=1 and K=0 eigenstates exhibit pronounced variations with the intermolecular rovibrational states, for which a qualitative explanation is provided. [ABSTRACT FROM AUTHOR]

A study of the intermolecular potential-energy surface (IPS) and the intermolecular states of the perprotonated and perdeuterated benzene–He complex is reported. From a fit to ab initio data computed within the coupled cluster singles and doubles including connected triples model for 280 interaction geometries, an analytic IPS including two- to four-body atom–atom terms is obtained. This IPS, and two other Lennard-Jones atom–atom surfaces from the literature, are each employed in dynamically exact (within the rigid-monomer approximation) calculations of J=0 intermolecular states of the isotopomers. Rotational constants and Raman-scattering coefficients for intermolecular vibrational transitions are also calculated for each of the three surfaces. The calculated results are compared with experimental results reported herein pertaining to intermolecular Raman spectra of benzene–He. The calculated rotational constants are compared with experimental values from the literature. The fitted IPS of this work leads to calculated observables that match the experimental results very well. The IPSs from the literature are not as successful, specifically in regard to the intermolecular Raman spectra. © 2003 American Institute of Physics. [ABSTRACT FROM AUTHOR]

The intermolecular potential-energy surface pertaining to the interaction between benzene and N[sub 2] is investigated theoretically and experimentally. Accurate intermolecular interaction energies are evaluated for the benzene-N[sub 2] van der Waals complex using the coupled cluster singles and doubles including connected triples [CCSD(T)] method and the aug-cc-pVDZ basis set extended with a set of 3s3p2d1f1g midbond functions. After fitting the energies to an analytic function, the intermolecular Schrödinger equation is solved to yield energies, rotational constants, and Raman-scattering coefficients for the lowest intermolecular levels of several benzene-N[sub 2] isotopomers. Experimentally, intermolecular Raman spectra of jet-cooled h[sub 6]- and d[sub 6]-benzene-N[sub 2] measured at 0.03 cm[sup -1] resolution by mass-selective, ionization-detected stimulated Raman spectroscopies are reported. Seven intermolecular bands are assigned for each isotopomer, including transitions involving intermolecular bending and stretching vibrations and internal rotation about the benzene C[sub 6] axis. These Raman data, together with measured rotational constants and binding energies obtained by other groups on benzene-N[sub 2], agree well with the theoretical results. Such agreement points to the promise of the quantum chemical methodology employed herein in future investigations of larger van der Waals complexes. [ABSTRACT FROM AUTHOR]

We evaluate the phenylacetylene-argon intermolecular potential energy surface by fitting a representative number of ab initio interaction energies to an analytic function. These energies are calculated at a grid of intermolecular geometries, using the CCSD(T) method and the aug-cc-pVDZ basis set extended with a series of 3s3p2d1f1g midbond functions. The potential is characterized by two equivalent global minima where the Ar atom is located above and below the phenylacetylene plane at a distance of 3.5781 Å from the molecular center of mass and at an angle of 9.08° with respect to the axis perpendicular to the phenylacetylene plane and containing the center of mass. The calculated interaction energy is -418.9 cm-1. To check further the potential, we obtain the rovibrational spectrum of the complex and the results are compared to the available experimental data. [ABSTRACT FROM AUTHOR]

Usingthe CCSD(T) model, we evaluated the intermolecular potentialenergy surfaces of the He–, Ne–, and Ar–phosgenecomplexes. We considered a representative number of intermoleculargeometries for which we calculated the corresponding interaction energieswith the augmented (He complex) and double augmented (Ne and Ar complexes)correlation-consistent polarized valence triple-ζ basis setsextended with a set of 3s3p2d1f1g midbond functions. These basis setswere selected after systematic basis set studies carried out at geometriesclose to those of the surface minima. The He–, Ne–,and Ar–phosgene surfaces were found to have absolute minimaof −72.1, −140.4, and −326.6 cm–1at distances between the rare-gas atom and the phosgene center ofmass of 3.184, 3.254, and 3.516 Å, respectively. The potentialswere further used in the evaluation of rovibrational states and therotational constants of the complexes, providing valuable resultsfor future experimental investigations. Comparing our results to thosepreviously available for other phosgene complexes, we suggest thatthe results for Cl2–phosgene should be revised. [ABSTRACT FROM AUTHOR]