Rotational study of the NH3–CO complex: Millimeter-wave measurements and ab initio calculations

The Radiation Chemistry of NH3···CO Complex in Cryogenic Media as Studied by Matrix Isolation

The NH₃···CO complex can be considered an important building block for cold synthetic astrochemistry leading to the formation of complex organic molecules, including key prebiotic species. In this work, we have studied the radiation-induced transformations of this complex in Ar, Kr, and Xe matrices using FTIR spectroscopy. On the basis of comparison with the quantum chemical calculations at the CCSD(T)/L2a_3 level of theory, it was found that the initial complex had the configuration with hydrogen bonding through the carbon atom of CO. Irradiation of the matrix isolated complex with X-rays at 6 K leads to the formation of a number of synthetic products, namely, HNCO (in all matrices), formamide NH₂CHO, NH₂CO, and HNCO-H₂ (in argon and krypton). The matrix effect on the product distribution was explained by the involvement of different excited states of the complex in their formation. It was suggested that formamide results from the singlet excited states while other species mainly originate from triplet excited states. The latter states are efficiently populated through ion-electron recombination (in all matrices) and through intersystem crossing (particularly, in xenon). High yield of the recombination triplet states is a feature of the processes induced by high-energy radiation (in contrast to direct photolysis). NCO, CN, and NO were found as minor secondary products at high adsorbed doses. The astrochemical implications of the obtained results are discussed.

Ab initio potential energy surface and microwave spectrum of the NH3-N2 van der Waals complex

We present a five-dimensional intermolecular potential energy surface (PES) of the NH₃-N₂ complex, bound state calculations, and new microwave (MW) measurements that provide information on the structure of this complex and a critical test of the potential. Ab initio calculations were carried out using the explicitly correlated coupled cluster [CCSD(T)-F12a] approach with the augmented correlation-consistent aug-cc-pVTZ basis set. The global minimum of the PES corresponds to a configuration in which the angle between the NH₃ symmetry axis and the intermolecular axis is 58.7° with the N atom of the NH₃ unit closest to the N₂ unit, which is nearly parallel to the NH₃ symmetry axis. The intermolecular distance is 7.01 a₀, and the binding energy D_e is 250.6 cm^-1. The bound rovibrational levels of the four nuclear spin isomers of the complex, which are formed when ortho/para (o/p)-NH₃ combines with (o/p)-N₂, were calculated on this intermolecular potential surface. The computed dissociation energies D₀ are 144.91 cm^-1, 146.50 cm^-1, 152.29 cm^-1, and 154.64 cm^-1 for (o)-NH₃-(o)-N₂, (o)-NH₃-(p)-N₂, (p)-NH₃-(o)-N₂, and (p)-NH₃-(p)-N₂, respectively. Guided by these calculations, the pure rotational transitions of the NH₃-N₂ van der Waals complex were observed in the frequency range of 13-27 GHz using the chirped-pulse Fourier-transform MW technique. A complicated hyperfine structure due to three quadrupole ¹⁴N nuclei was partly resolved and examined for all four nuclear spin isomers of the complex. Newly obtained data definitively established the K values (the projection of the angular momentum J on the intermolecular axis) for the lowest states of the different NH₃-N₂ nuclear spin isomers.

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.

Microwave spectra and nuclear quadrupole structure of the NH3-N2 van der Waals complex and its deuterated isotopologues

The microwave spectrum of the NH₃-N₂ van der Waals complex has been observed in a supersonic molecular jet expansion via broadband (2-8 GHz) chirped-pulse Fourier-transform microwave spectroscopy. Two pure rotational R(0) transitions (J = 1 - 0) with different hyperfine structure patterns were detected. One transition belongs to the (ortho)-NH₃-(ortho)-N₂ nuclear spin isomer in the ground K = 0 state reported earlier [G. T. Fraser et al., J. Chem. Phys. 84, 2472 (1986)], while another one is assigned to the (para)-NH₃-(para)-N₂ spin isomer in the K = 0 state not reported before (K is the projection of the total angular momentum J on the intermolecular axis). The complicated hyperfine structure arising from three quadrupole ¹⁴N nuclei of NH₃-N₂ was resolved for both transitions, and the quadrupole coupling constants associated with the NH₃ and N₂ subunits were precisely determined for the first time. These constants provided the dynamical information about the angular orientation of ammonia and nitrogen indicating that the average angle between the C ₃ axis of NH₃ and the N₂ axis is about 66°. The average van der Waals bond lengths are slightly different for (ortho)-NH₃-(ortho)-N₂ and (para)-NH₃-(para)-N₂ and amount to 3.678 Å and 3.732 Å, respectively. Similar results for the deuterated isotopologues, ND₃-N₂, NHD₂-N₂, and NH₂D-N₂, and their nuclear spin isomers were also obtained thus confirming and extending the analysis for the parent NH₃-N₂ complex.

Application of classical simulations for the computation of vibrational properties of free molecules

In this study, we investigate the ability of classical molecular dynamics (MD) and Monte-Carlo (MC) simulations for modeling the intramolecular vibrational motion. These simulations were used to compute thermally-averaged geometrical structures and infrared vibrational intensities for a benchmark set previously studied by gas electron diffraction (GED): CS₂, benzene, chloromethylthiocyanate, pyrazinamide and 9,12-I₂-1,2-closo-C₂B₁₀H₁₀. The MD sampling of NVT ensembles was performed using chains of Nose-Hoover thermostats (NH) as well as the generalized Langevin equation thermostat (GLE). The performance of the theoretical models based on the classical MD and MC simulations was compared with the experimental data and also with the alternative computational techniques: a conventional approach based on the Taylor expansion of potential energy surface, path-integral MD and MD with quantum-thermal bath (QTB) based on the generalized Langevin equation (GLE). A straightforward application of the classical simulations resulted, as expected, in poor accuracy of the calculated observables due to the complete neglect of quantum effects. However, the introduction of a posteriori quantum corrections significantly improved the situation. The application of these corrections for MD simulations of the systems with large-amplitude motions was demonstrated for chloromethylthiocyanate. The comparison of the theoretical vibrational spectra has revealed that the GLE thermostat used in this work is not applicable for this purpose. On the other hand, the NH chains yielded reasonably good results.