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.
- NSF-PAR ID:
- 10412654
- Date Published:
- Journal Name:
- Physical Chemistry Chemical Physics
- Volume:
- 24
- Issue:
- 40
- ISSN:
- 1463-9076
- Page Range / eLocation ID:
- 24655 to 24676
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
-
more » « less
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.
-
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
-
Symmetry-adapted perturbation theory (SAPT) is a popular and versatile tool to compute and decompose noncovalent interaction energies between molecules. The intramolecular SAPT (ISAPT) variant provides a similar energy decomposition between two nonbonded fragments of the same molecule, covalently connected by a third fragment. In this work, we explore an alternative approach where the noncovalent interaction is singled out by a range separation of the Coulomb potential. We investigate two common splittings of the 1/r potential into long-range and short-range parts based on the Gaussian and error functions, and approximate either the entire intermolecular/interfragment interaction or only its attractive terms by the long-range contribution. These range separation schemes are tested for a number of intermolecular and intramolecular complexes. We find that the energy corrections from range-separated SAPT or ISAPT are in reasonable agreement with complete SAPT/ISAPT data. This result should be contrasted with the inability of the long-range multipole expansion to describe crucial short-range charge penetration and exchange effects; it shows that the long-range interaction potential does not just recover the asymptotic interaction energy but also provides a useful account of short-range terms. The best consistency is attained for the error-function separation applied to all interaction terms, both attractive and repulsive. This study is the first step toward a fragmentation-free decomposition of intramolecular nonbonded energy.more » « less
-
Spectroscopy is a key analytical tool that provides valuable insight into molecular structure and is widely used to identify chemical samples. Tagging spectroscopy is a form of action spectroscopy in which the absorption of a single photon by a molecular ion is detected via the loss of a weakly attached, inert “tag” particle (e.g. He, Ne, N2).1–3 The absorption spectrum is derived from the tag loss rate as a function of incident radiation frequency. To date, all spectroscopy of gas phase polyatomic molecules has been restricted to large molecular ensembles, complicating spectral interpretation by the presence of multiple chemical and isomeric species. Here we present a novel tagging spectroscopic scheme to analyze the purest possible sample: a single gas phase molecule. We demonstrate this technique with the measurement of the infrared spectrum of a single tropylium (C7H7+ ) molecular ion; to our knowledge the first recorded spectrum of a single gas phase polyatomic molecule. Our method’s high sensitivity revealed spectral features previously unobserved using traditional tagging methods.4 Our approach in principle enables analysis of multi-component mixtures by identifying constituent molecules one at a time. Single molecule sensitivity extends action spectroscopy to rare samples, such as those of extraterrestrial origin,5,6 or to reactive reaction intermediates formed at number densities too low for traditional action methods.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1039/d2cp04005k
