Matrix Isolation and ab Initio Studies of 1:1 Hydrogen-Bonded Complexes HCN−H2O and HNC−H2O Produced by Photolysis of Formaldoxime

UV Laser-Induced Photodecomposition of Matrix-Isolated Salicylhydroxamic Acid: Identification of New Isocyanate Complexes

Photochemical reactions of salicylhydroxamic acid were induced using tunable UV laser radiation followed by FTIR spectroscopy. Four pairs of co-products were experimentally found to appear in the photolysis: C6H4(OH)NCO⋯H2O (1), C6H4(OH)C(O)N⋯H2O (2), C6H4(OH)2⋯HNCO (3), and C6H4(OH)NHOH⋯CO (4). The comparison of the theoretical spectra with the experimental ones allowed us to determine the structures of the complexes formed in the matrices. The mechanisms of the reaction channels leading to the formation of the photoproducts were proposed. It was concluded that the first step in the formation of the complexes (1), (2), and (3) was the scission of the N-O bond, whereas the creation of complex (4) was due to cleavage of the C-N bond.

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.

H2O-HCN complex: A new potential energy surface and intermolecular rovibrational states from rigorous quantum calculations

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.

Theoretical Investigation of Dissociation versus Intramolecular Rearrangements in Aminohydroxymethylene

Aminohydroxymethylene (H₂N-C̈-OH) is the simplest aminooxycarbene which is a heteroatom stabilized carbene. This highly reactive molecule was prepared in an Ar matrix in a recent experimental work. Unimolecular reactivity of this astrochemically important molecule was investigated and only fragmentations were identified contrary to the observations of both fragmentations and intramolecular rearrangements in other hydroxycarbenes. These rearrangement reactions form the corresponding imine and carbonyl compounds. In the present work, direct chemical dynamics simulations of unimolecular chemistry of aminohydroxymethylene were performed in the gas phase to study atomic level dissociation mechanisms. Classical trajectories were generated on-the-fly using potentials and gradients computed at the density functional B3LYP/6-31+G* level of electronic structure theory. Simulation results showed that intramolecular rearrangements accompany fragmentations during the unimolecular decay process of aminohydroxymethylene. However, the average lifetime of the intermediate isomers were found to be only few picoseconds which might not have been long enough for detection in the experiments.

Weakly Bound Complex Formation between HCN and CH3Cl: A Matrix-Isolation and Computational Study

The matrix-isolated infrared spectrum of a hydrogen cyanide-methyl chloride complex was investigated in a solid argon matrix. HCN and CH₃Cl were co-condensed onto a substrate held at 10 K with an excess of argon gas, and the infrared spectrum was measured using Fourier-transform infrared spectroscopy. Quantum chemical geometry optimization, harmonic frequency, and natural bonding orbital calculations indicate stabilized hydrogen- and halogen-bonded structures. The two resulting weakly bound complexes are both composed of one CH₃Cl molecule bound to a (HCN)₃ subunit, where the three HCN molecules are bound head-to-tail in a ring formation. Our study suggests that─in the presence of CH₃Cl─the formation of (HCN)₃ is promoted through complexation. Since HCN aggregates are an important precursor to prebiotic monomers (amino acids and nucleobases) and other life-bearing polymers, this study has astrophysical implications toward the search for life in space.

Direct evidence for a radiation-induced synthesis of acetonitrile and isoacetonitrile from a 1 : 1 CH4HCN complex at cryogenic temperatures: is it a missing link between inorganic and prebiotic astrochemistry?

Nitriles are important constituents of extraterrestrial media. Nitriles are supposed to play a crucial role in prebiotic chemistry occurring in the interstellar medium. In this work, we have investigated the low-temperature radiation-induced transformations of a 1 : 1 CH₄HCN complex as a plausible precursor of the simplest nitriles using the matrix isolation approach with FTIR spectroscopic detection. The parent complexes isolated in a noble gas (Ng) matrix were obtained by deposition of the CH₄/HCN/Ng gaseous mixture and characterized by comparison of experimental complexation-induced shifts of the HCN fundamentals with the results of the ab initio calculations. It was found that the X-ray irradiation of low-temperature matrices containing the isolated 1 : 1 CH₄HCN complex resulted in the formation of acetonitrile (CH₃CN) and isoacetonitrile (CH₃NC) and it appears to be the first experimental evidence for the formation of C₂ nitriles (acetonitrile and isoacetonitrile) from such a "building block". Additionally, a 1 : 1 CH₄HNC complex was tentatively assigned to the irradiated Ar and Kr matrices. It is demonstrated that the matrix has a strong effect on the CH₃CN/CH₃NC yield ratio, which dramatically increases in the row Ar < Kr < Xe. Also, the efficiency of the radiation-induced formation of the CH₄HNC complex was shown to decrease from Ar to Kr. It is believed that the proposed pathway for acetonitrile formation may be a significant step in the radiation-induced evolution leading to complex organic molecules and biomolecules under astrochemical conditions. Furthermore, the obtained results provide a prominent example of the impact of very weak intermolecular interactions on the radiation-induced transformations in cold media.

Aminohydroxymethylene (H2N-C̈-OH), the Simplest Aminooxycarbene

We generated and isolated hitherto unreported aminohydroxymethylene (1, aminohydroxycarbene) in solid Ar via pyrolysis of oxalic acid monoamide (2). Astrochemically relevant carbene 1 is persistent under cryogenic conditions and only decomposes to HNCO + H₂ and NH₃ + CO upon irradiation of the matrix at 254 nm. This photoreactivity is contrary to other hydroxycarbenes and aminomethylene, which undergo [1,2]H shifts to the corresponding carbonyls or imine. The experimental data are well supported by the results of CCSD(T)/cc-pVTZ and B3LYP/6-311++G(3df,3pd) computations.

Formation of Formamide from HCN + H2O: A Computational Study on the Roles of a Second H2O as a Catalyst, as a Spectator, and as a Reactant

Formamide (NH₂CHO), being the smallest and fundamental building block of life (with a peptide linkage), has recently been able to attract much interests, in the field of astrochemistry, astrophysics, and astrobiology. In this work, using quantum mechanical computations, reactions between HCN and H₂O, leading to the formation of formamide, have been analyzed. For the first time, an alternative and competing reaction channel, which proceeds via a geminal diol intermediate, for the formation of formamide, has been proposed. In this alternative channel, an extra water molecule (second H₂O) was found to be acting as a reactant, in the second step of the reaction path. Effects of second H₂O molecule in the reaction paths, providing catalytic assistance to the reaction or behaving like a spectator (concept is introduced for the first time for this reaction), have also been analyzed. Usefulness of spectator behavior is highlighted for the reactions happening on the rigid water-ice surfaces, where the water-ice may not be getting involved for any catalytic assistance. In light of catalytic assistances provided by the second H₂O, prominent effects in reducing the barrier heights drastically (even for the second step of the reaction, the barrier height was found to be below the reactants), through a hydrogen relay transport mechanism, were observed. In addition to the mechanism studies, interstellar feasibilities of all the reaction channels and their significances are discussed in detail.

Characterization of the HCNCO complex and its radiation-induced transformation to HNCCO in cold media: an experimental and theoretical investigation

The HCNCO complex and its X-ray induced transformation to HNCCO in solid noble gas (Ng) matrices (Ng = Ne, Ar, Kr, Xe) was first characterized by matrix isolation FTIR spectroscopy at 5 K. The HCNCO complex was obtained by deposition of HCN/CO/Ng gaseous mixtures. The assignment was based on extensive quantum chemical calculations at the CCSD(T) level of theory. The calculations predicted two computationally stable structures for HCNCO and three stable structures for HNCCO. However, only the most energetically favorable linear structures corresponding to the co-ordination between the H atom of HCN (HNC) and the C atom of CO have been found experimentally. The HCNCO complex demonstrates a considerable red shift of the H-C stretching vibrations (-24 to -38 cm^-1, depending on the matrix) and a blue shift of the HCN bending vibrations (+29 to +32 cm^-1) with respect to that of the HCN monomer, while the C[double bond, length as m-dash]O stretching mode is blue-shifted by 15 to 20 cm^-1 as compared to the CO monomer. The HNCCO complex reveals a strong red shift of the H-N bending (-77 to -118 cm^-1) and a strong blue shift of the HNC bending mode (ca. +100 cm^-1) as compared to the HNC monomer, whereas the C[double bond, length as m-dash]O stretching is blue-shifted by 24 to 29 cm^-1 with respect to that of the CO monomer. The interaction energies were determined to be 1.01 and 1.87 kcal mol^-1 for HCNCO and HNCCO, respectively. It was found that the formation of complexes with CO had a remarkable effect on the radiation-induced transformations of HCN. While the dissociation of HCN to H and CN is suppressed in complexes, the isomerization of HCN to HNC is strongly catalyzed by the complexation with CO. The astrochemical implications of the results are discussed.

Theoretical study of HCN-water interaction: five dimensional potential energy surfaces

A new five-dimensional potential energy surface is calculated at the coupled-cluster CCSD(T) level of theory for the HCN-water system, treating both monomers as rigid rotors. The associated methodology, which combines extensive ab initio calculations of moderate accuracy (CCSD(T)/AVDZ) and a fitting procedure involving a much lower angular coverage with more accurate ab initio calculations (CCSD(T)/CBS), is described in detail. This methodology provides a time-saving approach to compute quantitatively accurate potential energy surfaces with reasonable computational effort. Our potential reproduces the main features reported in the literature, and will allow us to perform the first quantum and semi-classical simulations of the collisional dynamic on this system.

Matrix-isolation studies of noncovalent interactions: more sophisticated approaches

Noncovalent interactions are crucial for many physical, chemical, and biological phenomena. Matrix isolation is a powerful method to study noncovalent interactions, including hydrogen-bonded species, and it has been extensively used in this field. However, there are difficult situations, such as in the case of species that are impossible to prepare in the gas phase. In this article, we describe some advanced approaches allowing studies of complexes that are problematic for the traditional methods. Photolysis of a suitable precursor in a matrix can lead to a large concentration of 1:1 complexes, which are otherwise very difficult to prepare (e.g., the H2O···O complex). Photolysis of species combined with annealing can lead to complexes of molecules with mobile atoms (e.g., the same H2O···O complex). Simultaneous photolysis of two species combined with annealing can produce complexes of radicals via reactions of the photogenerated complexes with mobile atoms (e.g., the H2O···HCO complex). Interaction of noble-gas (Ng) hydrides with other species is another topic (e.g., the N2···HArF complex) and very large blue shifts of the H-Ng stretching modes are normally observed for these systems. Complexes and dimers of the higher-energy conformer of formic acid have been prepared by using selective vibrational excitation of the ground-state conformer. The higher-energy conformer of formic acid can be efficiently stabilized in the complexes with strong hydrogen bonding. We also consider some problematic cases when the changes in the vibrational frequencies of the 1:1 complexes are very small (e.g., the phenol···Xe complex) and when the complex formation is prevented by strong solvation in the matrix (e.g., species in solid xenon).

Formaldoxime hydrogen bonded complexes with ammonia and hydrogen chloride

An infrared spectroscopic and MP2/6-311++G(2d,2p) study of hydrogen bonded complexes of formaldoxime with ammonia and hydrogen chloride trapped in solid argon matrices is reported. Both 1:1 and 1:2 complexes between formaldoxime and ammonia, hydrogen chloride have been identified in the CH2NOH/NH3/Ar, CH2NOH/HCl/Ar matrices, respectively, their structures were determined by comparison of the spectra with the results of calculations. In the 1:1 complexes present in the argon matrices the OH group of formaldoxime acts as a proton donor for ammonia and the nitrogen atom acts as a proton acceptor for hydrogen chloride. In the 1:2 complexes ammonia or hydrogen chloride dimers interact both with the OH group and the nitrogen atom of CH2NOH to form seven membered cyclic structures stabilized by three hydrogen bonds. The theoretical spectra generally agree well with the experimental ones, but they seriously underestimate the shift of the OH stretch for the 1:1 CH2NOH⋯NH3 complex.

Infrared spectroscopic characterization of hydrogen-bonded propylene oxide − ethanol and propylene oxide − 2-fluoroethanol complexes isolated in solid neon matrices

The infrared absorption spectra of hydrogen-bonded complexes of propylene oxide with either ethanol or 2-fluoroethanol have been recorded in neon matrices. Mixtures of propylene oxide and ethanol or propylene oxide and 2-fluoroethanol vapors were mixed with an excess of neon gas and deposited onto a KBr substrate at 4.2 K. The results indicate that hydrogen-bonded complexes were formed with propylene oxide as the hydrogen bond acceptor and either ethanol or 2-fluoroethanol as the hydrogen bond donors. The features assigned to the O−H stretch were red-shifted by 175 and 193 cm−1 for the ethanol- and 2-fluoroethanol-containing complexes, respectively. The difference in red shifts can be accounted for due to the greater acidity of 2-fluroethanol. Deuterium isotope experiments were conducted to help confirm the assignment of the O–H stretch for the complexes. As well, structures and infrared spectra were calculated using B3LYP/6-311++G(2d,2p) calculations and were used to compare with the experimental spectra. A “scaling equation” rather than a scaling factor was used and is shown to greatly increase the utility of the calculations when comparing with experimental spectra. An examination of the O–H stretching red shifts for many hydrogen-bound complexes reveals a relationship between the shift and the difference between the acidity of the hydrogen bond donor and the basicity of the hydrogen bond acceptor (the enthalpy of proton transfer). Both hydrogen-bonded complexes and proton-bound complexes appear to have a maximum in the reduced frequency value that corresponds to complexes where the hydrogen/proton are equally shared between the two bases.

Complexation of formaldoxime with water. Infrared matrix isolation and theoretical studies

The 1:1, 1:2 and 2:1 formaldoxime-water complexes isolated in the argon matrices have been studied by help of FTIR spectroscopy and MP2/6-311++G(2d,2p) method. The calculations predicted the stability of the three CH(2)NOH···H(2)O isomeric complexes, three CH(2)NOH···(H(2)O)(2) ones and one (CH(2)NOH)(2)···H(2)O complex. The analysis of the experimental spectra and their comparison with theoretical ones indicated that both the 1:1 and 1:2 complexes trapped in solid argon have the most stable cyclic structures stabilized by the O-H···O and O-H···N bonds between the formaldoxime and water molecules. In the 1:2 complex formaldoxime interacts with the water dimer, one H(2)O molecule acts as a proton acceptor for the OH group of formaldoxime whereas the second H(2)O molecule acts as a proton donor toward the nitrogen atom of the formaldoxime molecule. In the (CH(2)NOH)(2)···H(2)O complex the OH group of the water molecule acts as a proton donor toward one of the oxygen atoms of the formaldoxime cyclic dimer.

Water/cyanobutadiyne complexes: an infrared matrix isolation and theoretical study

The structures and energies of the 1:1 HC5N:H2O complexes in solid argon matrices have been investigated using FTIR spectroscopy and ab initio calculations, at the B3LYP/6-31G** and MP2/6-31G** levels of theory. Two types of 1:1 complexes are observed. The first one corresponds to the NH structure characterized by a hydrogen bond between H2O and the nitrogen of HC5N. The second corresponds to the OH form that involves a van der Waals interaction between the hydrogen of HC5N and the oxygen of water. HC5N can thus act either as an electrophile or as a nucleophile in complexes with water.

Infrared−Infrared Double Resonance Spectroscopy of the Isomers of Acetylene−HCN and Cyanoacetylene−HCN in Helium Nanodroplets

Infrared-infrared double resonance spectroscopy is used to probe the vibrational dynamics of molecular complexes solvated in helium nanodroplets. We report results for the acetylene-HCN and cyanoacetylene-HCN binary complexes, each having two stable isomers. We find that vibrational excitation of an acetylene-HCN complex results in a population transfer to the other isomer. Photoinduced isomerization is found to be dependent on both the initially excited vibrational mode and the identity of the acetylene-HCN isomer. However, population transfer is not observed for the cyanoacetylene-HCN complexes. The results are rationalized in terms of the ab initio intermolecular potential energy surfaces for the two systems with particular emphasis on the long-range barriers to rearrangement.

Ab Initio Study of the Isomeric Equilibrium of the HCN···H2O and H2O···HCN Hydrogen-Bonded Clusters

An ab initio study of the stability, spectroscopic properties, and isomeric equilibrium of the hydrogen-bonded HCN...H2O and H2O...HCN isomers is presented. Density functional theory and perturbative second-order MP2 and coupled-cluster CCSD(T) calculations were carried out and binding energies obtained with correlation-consistent basis sets including extrapolation to the infinity basis set level. At the best theoretical level, CCSD(T), the H2O...HCN complex is more stable than the HCN...H2O complex by ca. 6.3 kJ mol(-1). Rotational and vibrational spectra, including anharmonic corrections, are calculated. These calculated spectroscopic data are used to obtain thermochemical contributions to the thermodynamic functions and hence the Gibbs free energy. The relative free energies are used to estimate the equilibrium constant for isomerism. We find that under typical conditions of supersonic expansion experiments (T < 150 K) H2O...HCN is essentially the only isomer present. Furthermore, our calculations indicate that the hydrogen-bonded cluster becomes favorable over the separated moieties at temperatures below 200 K.

Matrix isolation Fourier transform infrared study of photodecomposition of formimidic acid

The UV isomerization of formamide (HCONH2) trapped in xenon, nitrogen, argon, and neon cryogenic matrices has been monitored by Fourier transform infrared (FT-IR) spectroscopy. Formamide monomer is the only species present in the matrices after deposition; when UV-selective irradiation was carried out at 240 nm, the n --> pi transition allowed us to observe the formation of several isomers of formimidic acid [H(OH)C=NH]. On these latter species, we carried out selective IR irradiation of their OH stretching mode and compared the experimental and theoretical (B3LYP/6-311+G(2d,2p)) sets of bands. This study allowed us to characterize for the first time all the isomers of formimidic acid. We have then studied the vacuum UV photodecomposition (lambda > 160 nm) of this molecule at 10 K in argon and xenon matrices. Several primary photoproducts such as HCN.H2O, HNC.H2O, and HNCO.H2 complexes, yielded by dehydration and dehydrogenation processes, were characterized.

Theoretical study of photoacidity of HCN: the effect of complexation with water

The character of the hydrogen bonding and the excited state proton transfer (ESPT) in the model system HCN...H(2)O is investigated. The PES of the two lowest excited states of the H(2)O...HCN complex was calculated using the CASPT2 method. The nonadiabatic coupling of the two states of the (pi-->pi*) and (pi-->sigma*) character is responsible for the excited state proton/hydrogen transfer. Compared to the ground state, the barrier for this process is significantly smaller. An increased number of water molecules in the complex with cyclic hydrogen-bonded network causes a large blue shift of the state of the (pi-->sigma*) character. The question of the dissociation of the complex in its excited state is also addressed.

Vacuum Ultraviolet (VUV) Photodecomposition of Urea Isolated in Cryogenic Matrix: First Detection of Isourea

Vacuum ultraviolet (VUV) irradiation at wavelengths of lambda > 160 nm of urea-h4 (NH2CONH2) and urea-d4 (ND2COND2) has been monitored by Fourier transform infrared spectroscopy in argon and xenon matrixes. Several primary photoproducts, such as HNCO:NH3 (isocyanic acid:ammonia), CO:N2H4 (carbon monoxide:hydrazine) molecular complexes, and isourea (H2N(OH)C=NH), which is reported for the first time, were characterized. The assignment of complexes was achieved by co-depositing the pairs of respective species, whereas the isourea identification was based on the comparison between the experimental and theoretical (B3LYP) infrared spectra. Isourea is found in the argon matrix in its most stable (s-Z)-(E) configuration. It is an intermediate in the VUV decomposition process; its dehydration leads to the NH2CN:H2O complex. In the xenon matrix, the photochemistry of urea yields the HNCO:NH3 complex as a major product, whereas the CO:N2H4 complex is observed in trace amounts. The observed differences between the argon and xenon matrixes suggest the crossing between S1 and T1 potential surfaces of urea to be responsible for the formation of the HNCO:NH3 complex. A comparison is also performed with other carboxamides, such as formamide (HCONH2) or acetamide (CH3CONH2).