Mixed quantum/classical calculations of rotationally inelastic scattering in the CO + CO system: a comparison with fully quantum results.

An updated version of the CO + CO potential energy surface from [R. Dawes, X. G. Wang and T. Carrington, J. Phys. Chem. A 2013, 117, 7612] is presented, that incorporates an improved treatment of the asymptotic behavior. It is found that this new surface is only slightly different from the other popular PES available for this system in the literature [G. W. M. Vissers, P. E. S. Wormer and A. Van Der Avoird, Phys. Chem. Chem. Phys. 2003, 5, 4767]. The differences are quantified by expanding both surfaces over a set of analytic functions and comparing the behavior of expansion coefficients along the molecule-molecule distance R. It is shown that all expansion coefficients behave similarly, except in the very high energy range at small R where the PES is repulsive. That difference has no effect on low collision-energy dynamics, which is explored via inelastic scattering calculations carried out using the MQCT program which implements the mixed quantum/classical theory for molecular energy exchange processes. The validity of MQCT predictions of state-to-state transition cross sections for CO + CO is also tested by comparison against full-quantum coupled-states calculations. In all cases MQCT gives reliable results, except at very low collision energy where the full-quantum calculations predict strong oscillations of state-to-state transition cross sections due to resonances. For strong transitions with large cross sections, the results of MQCT are reliable, especially at higher collision energy. For weaker transitions, and lower collision energies, the cross sections predicted by MQCT may be up to a factor of 2-3 different from those obtained by full-quantum calculations.

Rovibrational states calculations of the H2O-HCN heterodimer with the multiconfiguration time dependent Hartree method.

Water and hydrogen cyanide are two of the most common species in space and the atmosphere with the ability of binding to form dimers such as H2O-HCN. In the literature, while calculations characterizing various properties of the H2O-HCN cluster (equilibrium distance, vibrational frequencies and rotational constants) have been done in the past, extensive calculations of the rovibrational states of this system using a reliable quantum dynamical approach have yet to be reported. In this work, we intend to mend that by performing the first calculation of the rovibrational states of the H2O-HCN van der Waals complex on a recently developed potential energy surface. We use the block improved relaxation procedure implemented in the Heidelberg MultiConfiguration Time-Dependent Hartree (MCTDH) package to compute the states of the H2O-HCN isomer, from which we extract the transition frequencies and rotational constants of the complex. We further adapt an approach first suggested by Wang and Carrington-and supported here by analysis routines of the Heidelberg MCTDH package-to properly characterize the computed rovibrational states. The subsequent assignment of rovibrational states was done by theoretical analysis and visual inspection of the wavefunctions. Our simulations provide a Zero Point Energy (ZPE) and intermolecular vibrational frequencies in good agreement with past ab initio calculations. The transition frequencies and rotational constants obtained from our simulations match well with the available experimental data. This work has the broad aim to propose the MCTDH approach as a reliable option to compute and characterize rovibrational states of van der Waals complexes such as the current one.

Ab initio quantum scattering calculations and a new potential energy surface for the HCl(X1Σ+)-O2(X3Σg-) system: Collision-induced line shape parameters for O2-perturbed R(0) 0-0 line in H35Cl.

The remote sensing of abundance and properties of HCl-the main atmospheric reservoir of Cl atoms that directly participate in ozone depletion-is important for monitoring the partitioning of chlorine between "ozone-depleting" and "reservoir" species. Such remote studies require knowledge of the shapes of molecular resonances of HCl, which are perturbed by collisions with the molecules of the surrounding air. In this work, we report the first fully quantum calculations of collisional perturbations of the shape of a pure rotational line in H35Cl perturbed by an air-relevant molecule [as the first model system we choose the R(0) line in HCl perturbed by O2]. The calculations are performed on our new highly accurate HCl(X1Σ+)-O2(X3Σg-) potential energy surface. In addition to pressure broadening and shift, we also determine their speed dependencies and the complex Dicke parameter. This gives important input to the community discussion on the physical meaning of the complex Dicke parameter and its relevance for atmospheric spectra (previously, the complex Dicke parameter for such systems was mainly determined from phenomenological fits to experimental spectra and the physical meaning of its value in that context is questionable). We also calculate the temperature dependence of the line shape parameters and obtain agreement with the available experimental data. We estimate the total combined uncertainties of our calculations at 2% relative root-mean-square error in the simulated line shape at 296 K. This result constitutes an important step toward computational population of spectroscopic databases with accurate ab initio line shape parameters for molecular systems of terrestrial atmospheric importance.

Temperature Dependence of the Electronic Absorption Spectrum of NO2.

The nitrogen dioxide (NO2) radical is composed of the two most abundant elements in the atmosphere, where it can be formed in a variety of ways including combustion, detonation of energetic materials, and lightning. Relevant also to smog and ozone cycles, together these processes span a wide range of temperatures. Remarkably, high-resolution NO2 electronic absorption spectra have only been reported in a narrow range below about 300 K. Previously, we reported [ J. Phys. Chem. A 2021, 125, 5519-5533] the construction of quasi-diabatic potential energy surfaces (PESs) for the lowest four electronic states (X̃, Ã, B̃, and C̃) of NO2. In addition to three-dimensional PESs based on explicitly correlated MRCI(Q)-F12/VTZ-F12 ab initio data, the geometry dependence of each component of the dipoles and transition dipoles was also mapped into fitted surfaces. The multiconfigurational time-dependent Hartree (MCTDH) method was then used to compute the 0 K electronic absorption spectrum (from the ground rovibrational initial state) employing those energy and transition dipole surfaces. Here, in an extension of that work, we report an investigation into the effects of elevated temperature on the spectrum, considering the effects of the population of rotationally and vibrationally excited initial states. The calculations are complemented by new experimental measurements. Spectral contributions from hundreds of rotational states up to N = 20 and from 200 individually-characterized vibrational states were computed. A spectral simulation tool was developed that enables modeling the spectrum at various temperatures─by weighting individual spectral contributions via the partition function, or for pure excited initial states, which can be probed via transient absorption spectroscopy. We validate these results against experimental absorption spectroscopy data at high temperatures, as well as via a new measurement from the (1,0,1) initial vibrational state.

Collisional excitation of C2H- by H2: New interaction potential and scattering calculations.

Interstellar anions play an important role in astrochemistry as being tracers of the physical and chemical conditions in cold molecular clouds and circumstellar gas. The local thermodynamic equilibrium is generally not fulfilled in media where anions are detected and radiative and collisional data are required to model the observed lines. The C2H- anion has not yet been detected in the interstellar medium; however, collisional data could be used for non-LTE models that would help in identifying the most intense lines. For this purpose, we have computed the first 4D potential energy surface (PES) of the C2H--H2 complex using an explicitly correlated coupled-cluster approach. The PES is characterized by a single deep minimum with a well-depth of 924.96 cm-1. From this interaction potential, we derived excitation cross sections and rate coefficients of C2H- induced by collisions with para- and ortho-H2. The results obtained for collisions with para-H2 are compared to previous calculations performed using a 2D-PES obtained from an average over H2 rotations.

Ab initio quantum scattering calculations for the CO-O2 system and a new CO-O2 potential energy surface: O2 and air broadening of the R(0) line in CO.

We present ab initio calculations of the collisional broadening of the R(0) pure rotational line in CO (at 115 GHz) perturbed by O2. Our calculations are done in a fully quantum way by solving close-coupling quantum-scattering equations without any approximations. We also report a new, highly accurate CO-O2 potential energy surface on which we did the quantum-scattering calculations. The calculated collisional broadening agrees with the available experimental data in a wide temperature range. The calculated collisional shift is negligible compared to the broadening, which is also consistent with the experimental data. We combine this result with our previous calculations for the same line in CO perturbed by N2 [Józwiak et al., J. Chem. Phys. 154, 054314 (2021)]; the obtained air-perturbed broadening of the R(0) pure rotational line in CO and its temperature dependence perfectly agree with the HITRAN database. This result constitutes an important step toward developing a methodology for providing accurate ab initio reference data on spectroscopic collisional line-shape parameters for molecular systems relevant to the Earth's atmosphere and for populating spectroscopic line-by-line databases.

Theoretical study of the CO2-O2 van der Waals complex: potential energy surface and applications.

A four-dimensional-potential energy surface (4D-PES) of the atmospherically relevant carbon dioxide-oxygen molecule (CO2-O2) van der Waals complex is mapped using the ab initio explicitly correlated coupled cluster method with single, double, and perturbative triple excitations (UCCSD(T)-F12b), and extrapolation to the complete basis set (CBS) limit using the cc-pVTZ-F12/cc-pVQZ-F12 bases and the l-3 formula. An analytic representation of the 4D-PES was fitted using the method of interpolating moving least squares (IMLS). These calculations predict that the most stable configuration of CO2-O2 complex corresponds to a planar slipped-parallel structure with a binding energy of V ∼ -243 cm-1. Another isomer is found on the PES, corresponding to a non-planar cross-shaped structure, with V ∼ -218 cm-1. The transition structure connecting the two minima is found at V ∼ -211 cm-1. We also performed comparisons with some CO2-X van der Waals complexes. Moreover, we provide a SAPT analysis of this molecular system. Then, we discuss the complexation induced shifts of CO2 and O2. Afterwards, this new 4D-PES is employed to compute the second virial coefficient including temperature dependence. A comparison between quantities obtained in our calculations and those from experiments found close agreement attesting to the high quality of the PES and to the importance of considering a full description of the anisotropic potential for the derivation of thermophysical properties of CO2-O2 mixtures.

Unimolecular dissociation dynamics of electronically excited HCO(Ã2A''): rotational control of nonadiabatic decay.

The photoinduced unimolecular decay of the electronically excited HCO(Ã2A'') is investigated in a combined experimental-theoretical study. The molecule is excited to the (1, n2, 0) combination bands, which decay via Renner-Teller coupling to the ground electronic state. The rovibrational state distribution of the CO fragment was measured via the high-n Rydberg H-atom time-of-flight method and calculated using a wave packet method on an accurate set of potential energy surfaces. It is shown that the non-adiabatic decay rate is strongly modulated by the HCO rotational angular momentum, which leaves unique signatures in the product state distribution. The experimentally observed bimodal rotational distribution of the dominant CO(v = 0) fragment is likely due to decay of different vibronic states populated by the excitation and modulated by the excited state lifetime, which is in turn controlled by the parent rotational quantum number.

Fully quantum calculations of O2-N2 scattering using a new potential energy surface: Collisional perturbations of the oxygen 118 GHz fine structure line.

A proper description of the collisional perturbation of the shapes of molecular resonances is important for remote spectroscopic studies of the terrestrial atmosphere. Of particular relevance are the collisions between the O2 and N2 molecules-the two most abundant atmospheric species. In this work, we report a new highly accurate O2(X3Σg -)-N2(X1Σg +) potential energy surface and use it for performing the first quantum scattering calculations addressing line shapes for this system. We use it to model the shape of the 118 GHz fine structure line in O2 perturbed by collisions with N2 molecules, a benchmark system for testing our methodology in the case of an active molecule in a spin triplet state. The calculated collisional broadening of the line agrees well with the available experimental data over a wide temperature range relevant for the terrestrial atmosphere. This work constitutes a step toward populating the spectroscopic databases with ab initio line shape parameters for atmospherically relevant systems.

Rotational Modulation of Ã2A″-State Photodissociation of HCO via Renner-Teller Nonadiabatic Transitions.

By examining the product-state distribution of a prototypical nonadiabatic predissociation system, HCO(Ã2A″-X̃2A'), we demonstrate that the dissociation dynamics is strongly modulated by parent rotational quantum numbers. The predissociation of the nominal (νC-H = 0, νbend, νC-O = 1) vibronic levels of the Ã2A″ state surprisingly gives rise to both vibrational ground and excited states of the CO product, despite the assumed spectator nature of the CO moiety. This anomaly is attributed to the dependence of the lifetime of the vibronic resonance facilitated by the Renner-Teller interaction on the parent rotational angular momentum quantum numbers coupled with transient intensity borrowing from nearby vibronic resonances with νC-O = 0. This unique phenomenon is a purely quantum mechanical behavior that has no classical analogue.

The Low-Lying Electronic States of NO2: Potential Energy and Dipole Surfaces, Bound States, and Electronic Absorption Spectrum.

Nitrogen dioxide, NO2, is a free radical composed of the two most abundant elements in Earth's atmosphere, nitrogen and oxygen, and is relevant to atmospheric and combustion chemistry. The electronic structure of even its lowest-lying states is remarkably complex, with various conical intersections and Renner-Teller pairings, giving rise to complex and perturbed vibronic states. Here we report some analysis of the 18 molecular states of doublet spin-multiplicity formed by combining ground-state N(4Su) and O(3Pg) atoms. Three-dimensional potential energy surfaces were fit at the MRCI(Q)-F12/VTZ-F12 level, describing the lowest four (X̃, Ã, B̃, and C̃) electronic states. A properties-based diabatization procedure was applied to accommodate the intersections, producing energies in a quasidiabatic representation and yielding couplings that were also fit into surfaces. The low-lying vibrational levels on the ground X̃ state were computed and compared with experimental measurements. Compared to experiment, the lowest 125 calculated vibrational levels (up to 8500 cm-1 above the zero-point energy) have a root-mean-squared error of 16.5 cm-1. In addition, dipole moments for each of the lowest four electronic states-and the transition dipoles between them-were also computed and fit. With the coupled energy and dipole surfaces, the electronic spectrum was calculated in absolute intensity and compared with experimental measurements. Detailed structure in the experimental spectrum was successfully reproduced, and the total integrated intensity matches experiment to an accuracy of ∼1.5% with no empirical adjustments.

Inelastic scattering in isotopologues of O2-Ar: the effects of mass, symmetry, and density of states.

The two species considered here, O2 (oxygen molecule) and Ar (argon-atom), are both abundant components of Earth's atmosphere and hence familiar collision partners in this medium. O2 is quite reactive and extensively involved in atmospheric chemistry, including Chapman's cycle of the formation and destruction of ozone; while Ar, like N2, typically plays the nevertheless crucial role of inert collider. Inert species can provide stabilization to metastable encounter-complexes through the energy transfer associated with inelastic collisions. The interplay of collision frequency and energy transfer efficiency, with state lifetimes and species concentrations, contributes to the rich and varied chemistry and dynamics found in diverse environments ranging from planetary atmospheres to the interstellar and circumstellar media. The nature and density of bound and resonance states, coupled electronic states, symmetry, and nuclear spin-statistics can all play a role. Here, we systematically investigate some of those factors by looking at the O2-Ar system, comparing rigorous quantum-scattering calculations for the 16O16O-40Ar, 18O16O-40Ar, and 18O18O-40Ar isotope combinations. A new accurate potential energy surface was constructed for this purpose holding the O2 bond distance at its vibrationally averaged distance.

Collisional excitation of interstellar PN by H2: New interaction potential and scattering calculations.

Rotational excitation of interstellar PN molecules induced by collisions with H2 is investigated. We present the first ab initio four-dimensional potential energy surface (PES) for the PN-H2 van der Waals system. The PES was obtained using an explicitly correlated coupled cluster approach with single, double, and perturbative triple excitations [CCSD(T)-F12b]. The method of interpolating moving least squares was used to construct an analytical PES from these data. The equilibrium structure of the complex was found to be linear, with H2 aligned at the N end of the PN molecule, at an intermolecular separation of 4.2 Å. The corresponding well-depth is 224.3 cm-1. The dissociation energies were found to be 40.19 cm-1 and 75.05 cm-1 for complexes of PN with ortho-H2 and para-H2, respectively. Integral cross sections for rotational excitation in PN-H2 collisions were calculated using the new PES and were found to be strongly dependent on the rotational level of the H2 molecule. These new collisional data will be crucial to improve the estimation of PN abundance in the interstellar medium from observational spectra.

Spectroscopy and Scattering Studies Using Interpolated Ab Initio Potentials.

The Born-Oppenheimer potential energy surface (PES) has come a long way since its introduction in the 1920s, both conceptually and in predictive power for practical applications. Nevertheless, nearly 100 years later-despite astonishing advances in computational power-the state-of-the-art first-principles prediction of observables related to spectroscopy and scattering dynamics is surprisingly limited. For example, the water dimer, (H2O)2, with only six nuclei and 20 electrons, still presents a formidable challenge for full-dimensional variational calculations of bound states and is considered out of reach for rigorous scattering calculations. The extremely poor scaling of the most rigorous quantum methods is fundamental; however, recent progress in development of approximate methodologies has opened the door to fairly routine high-quality predictions, unthinkable 20 years ago. In this review, in relation to the workflow of spectroscopy and/or scattering studies, we summarize progress and challenges in the component areas of electronic structure calculations, PES fitting, and quantum dynamical calculations.

Computational study of the rovibrational spectrum of CO2-N2.

The CO2-N2 complex is formed from two key components of Earth's atmosphere, and as such, has received some attention from both experimental and theoretical studies. On the theory side, a potential energy surface (PES) based on high level ab initio data was reported [Nasri et al., J. Chem. Phys., 2015, 142, 174301] and then used in more recently reported rovibrational calculations [Lara-Moreno et al., Phys. Chem. Chem. Phys., 2019, 21, 3550]. Accuracy of about 1 percent was achieved for calculated rotational transitions of the ground vibrational state of the complex, compared with previously reported microwave spectra. However, a very recent measurement of the geared bending mode frequency [Barclay et al., J. Chem. Phys., 2020, 153, 014303] recorded a value of 21.4 cm-1, which is wildly different from the corresponding calculated value of 45.9 cm-1. To provide some insight into this discrepancy, we have constructed a new more accurate PES, and used it to perform highly converged variational rovibrational calculations. Our new results yield a value of 21.1 cm-1 for that bending frequency, in close agreement with the experiment. We also obtain significantly improved predicted rotational transitions. Finally, we note that a very shallow well, previously reported as a distinct second isomer, is not found on our new PES, but rather a transition structure is seen in that location.

Origin of the "odd" behavior in the ultraviolet photochemistry of ozone.

The origin of the even-odd rotational state population alternation in the 16O2(a1Δg) fragments resulting from the ultraviolet (UV) photodissociation of 16O3, a phenomenon first observed over 30 years ago, has been elucidated using full quantum theory. The calculated 16O2(a1Δg) rotational state distribution following the 266-nm photolysis of 60 K ozone shows a strong even-odd propensity, in excellent agreement with the new experimental rotational state distribution measured under the same conditions. Theory indicates that the even rotational states are significantly more populated than the adjacent odd rotational states because of a preference for the formation of the A' Λ-doublet, which can only occupy even rotational states due to the exchange symmetry of the two bosonic 16O nuclei, and thus not as a result of parity-selective curve crossing as previously proposed. For nonrotating ozone, its dissociation on the excited B1A' state dictates that only A' Λ-doublets are populated, due to symmetry conservation. This selection rule is relaxed for rotating parent molecules, but a preference still persists for A' Λ-doublets. The A''/A' ratio increases with increasing ozone rotational quantum number, and thus with increasing temperature, explaining the previously observed temperature dependence of the even-odd population alternation. In light of these results, it is concluded that the previously proposed parity-selective curve-crossing mechanism cannot be a source of heavy isotopic enrichment in the atmosphere.

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 just the floppy, and perhaps multiwelled potential energy surface (PES) characteristic of vdWs complexes, but in addition, one 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 wave function, 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.

A Global Full-Dimensional Potential Energy Surface for the K2Rb2 Complex and Its Lifetime.

A full-dimensional global potential energy surface for the KRb + KRb → K2 + Rb2 reaction is developed from 20 759 ab initio points calculated using a coupled cluster singles, doubles, and perturbative triples (CCSD(T)) method with effective core potentials, extrapolated to the complete basis set limit. The ab initio points are represented with high fidelity (root-mean-square error of 1.86 cm-1) using the permutation-invariant polynomial-neural network method, which enforces the permutation invariance of the potential with respect to exchange of identical nuclei. The potential energy surface features two D2h minima and one Cs minimum connected by the isomerization saddle points. The Rice-Ramsperger-Kassel-Marcus lifetime of the K2Rb2 reaction intermediate estimated using the potential energy surface is 227 ns, in reasonable agreement with the latest experimental measurement.

Rotationally inelastic scattering of O3-Ar: state-to-state rates with the multiconfigurational time dependent Hartree method.

The Chapman cycle, proposed in 1930, describes the various steps in the ongoing formation and destruction of stratospheric ozone. A key step in the formation process is the stabilization of metastable ozone molecules through collisions with a third body, usually an inert collider such as N2. The "ozone isotopic anomaly" refers to the observation of larger-than-expected atmospheric concentrations for certain ozone isotopologues. Previous studies point to the formation steps as the origin of this effect. A possibly key aspect of the ozone formation dynamics is that of the relative efficiencies of the collisional cooling of different isotopologues. Although the substitution of low-abundance 18O for 16O in O3 molecules corresponds to a relatively small net change in mass, related to this are some subtleties due to symmetry-breaking and a resulting more than doubling of the density of allowed states governed by nuclear-spin statistics for bosons. Recently, a highly accurate 3D potential energy surface (PES) describing O3-Ar interactions has been constructed and used to benchmark the low-lying rovibrational states of the complex. Here, using this new PES, we have studied the collisional energy-transfer dynamics using the MultiConfiguration Time Dependent Hartree method. A study of the rotationally inelastic scattering was performed for the parent 16O16O16O-Ar system and compared with that of the 16O16O18O-Ar isotopologue. The state-to-state cross-sections and rates from the 00,0 initial state to low lying excited states are reported. Analysis of these results yields insight into the interplay between small changes in the rotational constants of O3 and the reduced mass of the O3-Ar collision system, combined with that of the symmetry-breaking and introduction of a new denser manifold of allowed states.