skip to main content


Title: Automated Construction of Potential Energy Surfaces Suitable to Describe van der Waals Complexes With Highly-Excited Nascent Molecules: The Rotational Spectra of Ar–CS( v ) and Ar–SiS( v )
Some reactions produce extremely hot nascent-products which nevertheless can form sufficiently long-lived van der Waals (vdW) complexes—with atoms or molecules from a bath gas—as to be observed via microwave spectroscopy. Theoretical calculations of such unbound resonance-states can be much more challenging than ordinary bound-state calculations depending on the approach employed. One encounters not only the floppy, and perhaps multi-welled potential energy surface (PES) characteristic of vdWs complexes, but in addition must contend with excitation of the intramolecular modes and its corresponding influence on the PES. Straightforward computation of the (resonance) rovibrational levels of interest, involves the added complication of the unbound nature of the wavefunction, often treated with techniques such as introducing a complex absorbing potential. Here, we have demonstrated that a simplified approach of making a series of vibrationally effective PESs for the intermolecular coordinates—one for each reaction product vibrational quantum number of interest—can produce vdW levels for the complex with spectroscopic accuracy. This requires constructing a series of appropriately weighted lower-dimensional PESs for which we use our freely available PES-fitting code AUTOSURF. The applications of this study are the Ar–CS and Ar–SiS complexes, which are isovalent to Ar–CO and Ar–SiO, the latter of which we considered in a previously reported study. Using a series of vibrationally effective PESs, rovibrational levels and predicted microwave transition frequencies for both complexes were computed variationally. A series of shifting rotational transition frequencies were also computed as a function of the diatom vibrational quantum number. The predicted transitions were used to guide and inform an experimental effort to make complementary observations. Comparisons are given for the transitions that are within the range of the spectrometer and were successfully recorded. Calculations of the rovibrational level pattern agree to within 0.2 % with experimental measurements.  more » « less
Award ID(s):
1908576
NSF-PAR ID:
10149708
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
The Journal of Physical Chemistry A
ISSN:
1089-5639
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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.

     
    more » « less
  2. 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. 
    more » « less
  3. 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.

     
    more » « less
  4. Two-color, two-photon laser-induced fluorescence experiments were performed to probe the intermolecular interactions within the Ar + I2(E, vE = 0–3) potential energy surfaces. Spectra were recorded using the lowest-energy T-shaped level and an excited intermolecular vibrational level with bending excitation within the Ar + I2(B, vB = 23) potential as intermediate levels to guide the spectral assignments. Progressions of intermolecular stretching and bending levels bound within the Ar + I2(E, vE) potentials were identified, and their vibrational frequencies were determined. The harmonic frequency and anharmonic constant for the bending vibrational mode were determined to be ωe(b) ∼ 34.8 cm−1 and ωeχe(b) ∼ 0.3 cm−1. The frequency and anharmonic constant for the stretching mode were found to be the same as reported previously [V.V. Baturo, et al. Chem. Phys. Lett. 647 (2016) 161], ωe(s) = 37.2(1.1) cm−1 and ωeχe(s) = 1.8(2) cm−1. 
    more » « less
  5. 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.

     
    more » « less