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.
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
- 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
-
more » « less
-
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
-
Ion time-of-flight velocity-map imaging was used to measure the kinetic-energy distributions of the I2 ion-pair fragments formed after photoexcitation of Ar⋯I2 complexes to intermolecular vibrational levels bound within the Ar + I2 (E, vE = 0–2) potential energy surfaces. The kinetic-energy distributions of the I2 products indicate that complexes in the Ar⋯I2 (E, vE) levels preferentially dissociate into I2 in the D and β ion-pair states with no change in I2 vibrational excitation. The energetics of the levels prepared suggest that there is a non-adiabatic coupling of the initially prepared levels with the continuum of states lying above the Ar + I2 (D, vD = vE) and Ar + I2 (β, vβ = vE) dissociation limits. The angular anisotropies of the I2 product signals collected for many of the Ar⋯I2 (E, vE) levels have maxima parallel to the laser polarization axis. This contradicts expectations for the prompt dissociation of complexes with T-shaped geometries, which would result in images with maxima perpendicular to the polarization axis. These anisotropies suggest that there is a perturbation of the transition moment in these clusters or there are additional intermolecular interactions, likely those sampled while traversing above the attractive wells of the lower-energy potentials during dissociation. I2 (D′, vD′) products are also identified when preparing several of the low-lying levels localized in the T-shaped well of the Ar + I2 (E, vE = 0–2) potentials, and they are formed in multiple νD′ vibrational levels spanning energy ranges up to 500 cm−1.more » « less
-
New computational and experimental studies have been carried out for the MgCCH radical in its X2Σ+ state. Coupled cluster theory with single, double, and perturbative triples, CCSD(T), was used in conjunction with post-CCSD(T) and scalar relativistic additive corrections to compute vibrational quartic force fields for this molecule. From the quartic force fields, higher-order spectroscopic properties, including rotational constants, were obtained. In tandem, the five lowest energy rotational transitions for MgCCH, N = 1→0 through N = 5→4, were measured for the first time using Fourier transform microwave/millimeter wave methods in the frequency range 9 -50 GHz. The radical was created in the Discharge Assisted Laser Ablation Source (DALAS) developed in the Ziurys group. A combined fit of these data with previous millimeter direct absorption measurements have yielded the most accurate rotational constants for MgCCH to date. The computed principle rotational constant lies within 1.51-1.65 MHz of the experimental one, validating the computational approach. High-level theory was then applied to produce accurate rovibrational spectroscopic constants for MgCCH+, including a rotational constant of B0 = 5354.5–5359.5 MHz.. These new predictions will further the experimental study of MgCCH+, and aid in the low-temperature characterization of MgCCH, detected towards the circumstellar shell of IRC+10216, a carbon-rich star.more » « less
-
more » « less
ABSTRACT This paper presents rate coefficients for transitions between rotational levels of the A-type and E-type nuclear spin modifications of methanol induced by collisions with molecular hydrogen. These rate coefficients are required for an accurate determination of methanol abundance in the interstellar medium, where local thermodynamic equilibrium conditions generally do not apply. Time-independent close-coupling quantum scattering calculations have been employed to calculate cross-sections and rate coefficients for the (de-)excitation of methanol in collisions with para- and ortho-H2. These calculations utilized a potential energy surface (PES) for the interaction of methanol with H2 recently computed by the explicitly correlated CCSD(T)-F12a coupled-cluster method that employed a correlation-consistent aug-cc-pVTZ basis. Rate coefficients for temperatures ranging from 3 to 250 K were calculated for all transitions among the first 76 rotational levels of both A-type and E-type methanol, whose energies are less than or equal to 170 K. These rate coefficients are compared with those by Rabli and Flower who carried out coupled-state calculations using a PES computed by second-order many-body perturbation theory. Simple radiative transfer calculations using the present set of rate coefficients are also reported and compared with such calculations using the rate coefficients previously computed by Rabli and Flower.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1021/acs.jpca.0c02685
