Nitrous oxide dimer: A new potential energy surface and rovibrational spectrum of the nonpolar isomer

Theory cracks old data: Rovibrational energy levels of orthoH2-CO derived from experiment

Measurements of rovibrational spectra of clusters provide physical insight only if spectral lines can be assigned to pairs of quantum states, and further insight is obtained if one can deduce the quantitative energy-level pattern. Both steps can be so difficult that some measured spectra remain unassigned, one example is orthoH2-CO. To extend the scope of spectroscopic insights, we propose to use theoretical information in interpretation of spectra. We first performed high accuracy, full-dimensional calculations of the orthoH2-CO spectrum, at the highest practically achievable levels of electronic structure theory and quantum nuclear dynamics. Then, an iterative, theory-guided method developed here allowed us to fully interpret the spectrum of orthoH2-CO, extending the range of van der Waals clusters for which spectroscopy can provide physical insights.

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.

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.

Ab Initio Calculations of the Interaction Potential of the N2O-N2O Dimer: Strength of the Intermolecular Interactions and Physical Insights

N₂O, or nitrous oxide, is an important greenhouse gas with a significant impact on global warming and climate change. To accurately model the behavior of N₂O in the atmosphere, precise representations of its intermolecular force fields are required. First principles quantum mechanical calculations followed by appropriate fitting are commonly used to establish such force fields. However, fitting such force fields is challenging due to the complex mathematical functions that describe the molecular interactions of N₂O. As such, ongoing research is focused on improving our understanding of N₂O and developing more accurate models for use in climate modeling and other applications. In this study, we investigated the strength of the intermolecular interactions in the N₂O-N₂O dimer using the coupled-cluster theory with single, double, and perturbative triple excitation [CCSD(T)] method with the def2-QZVPP basis set. Our calculations provided a detailed understanding of the intermolecular forces that govern the stability and structure of the N₂O dimer. We found that the N₂O-N₂O dimer is stabilized by a combination of van der Waals forces and dipole-dipole interactions. The calculated interaction energy between the two N₂O molecules in the dimer was found to be -5.09 kcal/mol, which is in good agreement with previous theoretical and experimental results. Additionally, we analyzed the molecular properties of the N₂O-N₂O dimer, including its geometry and charge distribution. Our calculations provide a comprehensive understanding of the intermolecular interactions in the N₂O-N₂O dimer using the CCSD(T) method with the def2-QZVPP basis set by using the improved Lennard-Jones interaction potential method. These results can be used to improve our understanding of atmospheric chemistry and climate modeling, as well as to aid in the interpretation of experimental data.

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 O₂ and N₂ molecules-the two most abundant atmospheric species. In this work, we report a new highly accurate O₂(X³Σ_g ^-)-N₂(X¹Σ_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 O₂ perturbed by collisions with N₂ 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.

Understanding the structure of complex multidimensional wave functions. A case study of excited vibrational states of ammonia

Generalization of an earlier reduced-density-matrix-based vibrational assignment algorithm is given, applicable for systems exhibiting both large-amplitude motions, including tunneling, and degenerate vibrational modes. The algorithm developed is used to study the structure of the excited vibrational wave functions of the ammonia molecule, ¹⁴NH₃. Characterization of the complex dynamics of systems with several degenerate vibrations requires reconsidering the traditional degenerate-mode description given by vibrational angular momentum quantum numbers and switching to a symmetry-based approach that directly predicts state degeneracy and uncovers relations between degenerate modes. Out of the 600 distinct vibrational eigenstates of ammonia obtained by a full-dimensional variational computation, the developed methodology allows for the assignment of about 500 with meaningful labels. This study confirms that vibrationally excited states truly have modal character recognizable up to very high energies even for the non-trivial case of ammonia, a molecule which exhibits a tunneling motion and has two two-dimensional normal modes. The modal characteristics of the excited states and the interplay of the vibrational modes can be easily visualized by the reduced-density matrices, giving an insight into the complex modal behavior directed by symmetry.

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.

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, (H₂O)₂, 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.

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

Rotational excitation of interstellar PN molecules induced by collisions with H₂ is investigated. We present the first ab initio four-dimensional potential energy surface (PES) for the PN-H₂ 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 H₂ 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-H₂ and para-H₂, respectively. Integral cross sections for rotational excitation in PN-H₂ collisions were calculated using the new PES and were found to be strongly dependent on the rotational level of the H₂ molecule. These new collisional data will be crucial to improve the estimation of PN abundance in the interstellar medium from observational spectra.

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.

Observation of the elusive "oxygen-in" OCS dimer

The carbonyl sulfide (OCS) dimer serves as a prototype system for studying intermolecular forces between nonsymmetrical linear polyatomic molecules. Here, we performed a laser spectroscopic investigation of OCS dimers embedded in helium nanodroplets and found rovibrational bands corresponding to the non-polar "sulfur-in" and parallel polar dimers that have been extensively characterized in the gas phase, as well as a new non-polar "oxygen-in" dimer that has long been predicted by theory. Frequency alternations in the rotational branches along with the absence of a Stark effect provided important clues as to its assignment.

A combined on-the-fly/interpolation procedure for evaluating energy values needed in molecular simulations

We propose an algorithm for molecular dynamics or Monte Carlo simulations that uses an interpolation procedure to estimate potential energy values from energies and gradients evaluated previously at points of a simplicial mesh. We chose an interpolation procedure that is exact for harmonic systems and considered two possible mesh types: Delaunay triangulation and an alternative anisotropic triangulation designed to improve performance in anharmonic systems. The mesh is generated and updated on the fly during the simulation. The procedure is tested on two-dimensional quartic oscillators and on the path integral Monte Carlo evaluation of the HCN/DCN equilibrium isotope effect.

Computational study of the ro-vibrational spectrum of CO-CO2

An accurate ab initio ground-state intermolecular potential energy surface (PES) was determined for the CO-CO₂ van der Waals dimer. The Lanczos algorithm was used to compute rovibrational energies on this PES. For both the C-in and O-in T-shaped isomers, the fundamental transition frequencies agree well with previous experimental results. We confirm that the in-plane states previously observed are geared states. In addition, we have computed and assigned many other vibrational states. The rotational constants we determine from J = 1 energy levels agree well with their experimental counterparts. Planar and out-of-plane cuts of some of the wavefunctions we compute are quite different, indicating strong coupling between the bend and torsional modes. Because the stable isomers are T-shaped, vibration along the out-of-plane coordinates is very floppy. In CO-CO₂, when the molecule is out-of-plane, interconversion of the isomers is possible, but the barrier height is higher than the in-plane geared barrier height.

Full quantum calculation of the rovibrational states and intensities for a symmetric top-linear molecule dimer: Hamiltonian, basis set, and matrix elements

The rovibrational energy levels and intensities of the CH₃F-H₂ dimer have been obtained using our recent global intermolecular potential energy surface [X.-L. Zhang et al., J. Chem. Phys. 148, 124302 (2018)]. The Hamiltonian, basis set, and matrix elements are derived and given for a symmetric top-linear molecule complex. This approach to the generation of energy levels and wavefunctions can readily be utilized for studying the rovibrational spectra of other van der Waals complexes composed of a symmetric top molecule and a linear molecule, and may readily be extended to other complexes of nonlinear molecules and linear molecules. To confirm our method, the rovibrational levels of the H₂O-H₂ dimer have been computed and shown to be in good agreement with experiment and with previous theoretical results. The rovibrational Schrödinger equation has been solved using a Lanczos algorithm together with an uncoupled product basis set. As expected, dimers containing ortho-H₂ are more strongly bound than dimers containing para-H₂. Energies and wavefunctions of the discrete rovibrational levels of CH₃F-paraH₂ complexes obtained from the direct vibrationally averaged 5-dimensional potentials are in good agreement with the results of the reduced 3-dimensional adiabatic-hindered-rotor (AHR) approximation. Accurate calculations of the transition line strengths for the orthoCH₃F-paraH₂ complex are also carried out, and are consistent with results obtained using the AHR approximation. The microwave spectrum associated with the orthoCH₃F-orthoH₂ dimer has been predicted for the first time.

Infrared bands of CS2 dimer and trimer at 4.5 μm

We report observation of new infrared bands of (CS₂)₂ and (CS₂)₃ in the region of the CS₂ ν₁ + ν₃ combination band (at 4.5 µm) using a quantum cascade laser. The complexes are formed in a pulsed supersonic slit-jet expansion of a gas mixture of carbon disulfide in helium. We have previously shown that the most stable isomer of (CS₂)₂ is a cross-shaped structure with D_2d symmetry and that for (CS₂)₃ is a barrel-shaped structure with D₃ symmetry. The dimer has one doubly degenerate infrared-active band in the ν₁ + ν₃ region of the CS₂ monomer. This band is observed to have a rather small vibrational shift of -0.844 cm^-1. We expect one parallel and one perpendicular infrared-active band for the trimer but observe two parallel bands and one perpendicular band. Much larger vibrational shifts of -8.953 cm^-1 for the perpendicular band and -8.845 cm^-1 and +16.681 cm^-1 for the parallel bands are observed. Vibrational shifts and possible vibrational assignments, in the case of the parallel bands of the trimer, are discussed using group theoretical arguments.

Development of a potential energy surface for the O3–Ar system: rovibrational states of the complex

Collisional stabilization is an important step in the process of atmospheric formation of ozone.

Dynamical interference in the vibronic bond breaking reaction of HCO

First-principles treatments of quantum molecular reaction dynamics have reached the level of quantitative accuracy even in cases with strong non-Born-Oppenheimer effects. This achievement permits the interpretation of puzzling experimental phenomena related to dynamics governed by multiple coupled potential energy surfaces. We present a combined experimental and theoretical study of the photodissociation of formyl radical (HCO). Oscillations observed in the distribution of product states are found to arise from the interference of matter waves-a manifestation analogous to Young's double-slit experiment.

Ab initio study of the CO–N2 complex: a new highly accurate intermolecular potential energy surface and rovibrational spectrum

We present a highly accurate ab initio intermolecular potential-energy surface and rovibrational spectrum for the CO–N₂ complex.

Infrared spectrum and intermolecular potential energy surface of the CO–O2 dimer

The spectrum of the weakly-bound radical complex CO–O₂ is studied for the first time.

The vibration-rotation-tunneling levels of N2-H2O and N2-D2O

In this paper, we report vibration-rotation-tunneling levels of the van der Waals clusters N2-H2O and N2-D2O computed from an ab initio potential energy surface. The only dynamical approximation is that the monomers are rigid. We use a symmetry adapted Lanczos algorithm and an uncoupled product basis set. The pattern of the cluster's levels is complicated by splittings caused by H-H exchange tunneling (larger splitting) and N-N exchange tunneling (smaller splitting). An interesting result that emerges from our calculation is that whereas in N2-H2O, the symmetric H-H tunnelling state is below the anti-symmetric H-H tunnelling state for both K = 0 and K = 1, the order is reversed in N2-D2O for K = 1. The only experimental splitting measurements are the D-D exchange tunneling splittings reported by Zhu et al. [J. Chem. Phys. 139, 214309 (2013)] for N2-D2O in the v2 = 1 region of D2O. Due to the inverted order of the split levels, they measure the sum of the K = 0 and K = 1 tunneling splittings, which is in excellent agreement with our calculated result. Other splittings we predict, in particular those of N2-H2O, may guide future experiments.

Ab initio intermolecular potential energy surface and thermophysical properties of nitrous oxide

We present an analytical intermolecular potential energy surface (PES) for two rigid nitrous oxide (N2O) molecules derived from high-level quantum-chemical ab initio calculations. Interaction energies for 2018 N2O-N2O configurations were computed utilizing the counterpoise-corrected supermolecular approach at the CCSD(T) level of theory using basis sets up to aug-cc-pVQZ supplemented with bond functions. A site-site potential function with seven sites per N2O molecule was fitted to the pair interaction energies. We validated our PES by computing the second virial coefficient as well as shear viscosity and thermal conductivity in the dilute-gas limit. The values of these properties are substantiated by the best experimental data.

Two new methods for computing vibrational energy levels

I review two new ideas for coping with the size of large product basis sets and large product grids when one computes vibrational energy levels. The first is based on a tensor reduction scheme. It exploits advantages of a sum-of-products potential. The key idea is to use a basis each of whose function is a sum of optimized products and to compress the number of terms in each basis function. When the potential does not have the sum-of-products form, it is usually necessary to use quadrature. The second idea uses a nondirect product grid that has structure and is therefore compatible with efficient matrix–vector products.

Can a single water molecule really affect the HO2 + NO2 hydrogen abstraction reaction under tropospheric conditions?

During the HO₂ + NO₂ reaction, hydrogen abstraction by a single water molecule not only changes the branching ratio of HONO and HNO₂ formation, but also introduces different features with respect to the naked reaction, acting as a reactant that leads to the production of HNO₃.