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

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.

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.

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.

A COMPARATIVE STUDY ON INTERMOLECULAR HYDROGEN BOND INTERACTIONS IN MOLECULAR DIMERS USING DIFFERENT LEVELS OF COMPUTATIONAL METHODS

Hydrogen bond interactions in biological systems are important scientific issues but are challenging for their theoretical determinations at quantum-mechanical level of theory. Due to the different approximations, the available theoretical approaches often predict diverse hydrogen bond lengths and strengths. In this work, we evaluated the reliabilities of a number of widely used theoretical approaches including HF, SVWN, BLYP, PW91, B3LYP, BH and HLYP, B97D, M06L, MP2, and DFTB-D in studying hydrogen bonding, by calculating the hydrogen bond lengths and binding energies of 23 dimers formed by HCOOH , NH3 and Glycine. We also compared the effects of STO-3G, 6-31+G**, 6-311++G** and 6-311++G(2df,2p) basis sets on the results. Our result shows that, M06L, B3LYP and BHandHLYP methods can predict accurate dimer structures with a moderate basis set. Moreover, DFTB-D also gives reasonably reliable results with high efficiency and satisfactory precision, being a good choice for studying complex structures which contain hydrogen bonds.

THE COOPERATIVITY BETWEEN HYDROGEN AND HALOGEN BONDS

The cooperativity between hydrogen bonds and halogen bonds in X–HCN–Y ( X: C2H2, H2O, NH3, HCI, HCN, HF; Y: HF, BrF, Br2 is analyzed with MP2/6-311++G(d, p) and DFT/6-311++G(d, p) calculations using the B3LYP and mPW1PW91 hybrid functionals. The results from the quantum chemical calculations are typically clustered in groups according to the Y-ligand. By choosing the X–HCN–HF group as reference it is possible to describe the interaction between the hydrogen and the halogen bond with a two-parameter model. The value of the first parameter of the model describes the contribution of the X -ligand to the interbond cooperativity in the reference cluster. The second parameter of our model quantifies the changes in interbond cooperativity upon varying the Y -ligand. This simple model can be used to predict the cooperativity in X–HCN–Y trimers with reasonable accuracy and thereby to organize the results systematically. It is further shown that the conclusions drawn from this ordering scheme are independent from the computational method and thereby generally applicable.

Nitrile bonds as infrared probes of electrostatics in ribonuclease S

Three different nitrile-containing amino acids, p-cyanophenylalanine, m-cyanophenylalanine, and S-cyanohomocysteine, have been introduced near the active site of the semisynthetic enzyme ribonuclease S (RNase S) to serve as probes of electrostatic fields. Vibrational Stark spectra, measured directly on the probe-modified proteins, confirm the predominance of the linear Stark tuning rate in describing the sensitivity of the nitrile stretch to external electric fields, a necessary property for interpreting observed frequency shifts as a quantitative measure of local electric fields that can be compared with simulations. The X-ray structures of these nitrile-modified RNase variants and enzymatic assays demonstrate minimal perturbation to the structure and function, respectively, by the probes and provide a context for understanding the influence of the environment on the nitrile stretching frequency. We examine the ability of simulation techniques to recapitulate the spectroscopic properties of these nitriles as a means to directly test a computational electrostatic model for proteins, specifically that in the ubiquitous Amber-99 force field. Although qualitative agreement between theory and experiment is observed for the largest shifts, substantial discrepancies are observed in some cases, highlighting the ongoing need for experimental metrics to inform the development of theoretical models of electrostatic fields in proteins.

IR spectroscopic testing of surfaces in water ice and in icy mixtures with prussic acid or ammonia

An experimental technique is developed for spectral studies of ice films deposited at 77 K onto the ZnSe inner windows of a cell, so that the spectra could be registered in the presence of gaseous adsorbate. Pure H2O or D2O ice as well as HCN/D2O and ND3/D2O mixed icy films with different dopant/water ratios (1:10, 1:5, and 1:1) were investigated. The surface area of ice deposited at 77 K H2O estimated from adsorption measurements was about 160 m2/g. Bands of dangling hydroxyl groups disappear on raising the temperature up to 130-160 K when the changes of bulk absorption provide evidence for a phase transition from amorphous to polycrystalline ice. Surface properties of icy films were characterized by low-temperature adsorption of CO and CHF3. Insertion of small doses of HCN or ammonia does not change the acid/base strength of dangling hydroxyl groups or coordinately unsaturated surface oxygen atoms, but it changes the proportion between the concentrations of these sites as compared with pure water ice. For high dopant concentrations, the dangling hydroxyls were not observed, and the dominant adsorption sites for CO are likely to be the unsaturated oxygen atoms, while serious structural changes occur in the bulk of ices.

Lewis acid-base interactions in weakly bound formaldehyde complexes with CO2, HCN, and FCN: considerations on the cooperative H-bonding effects

Ab initio quantum chemistry calculations reveal that HCN and mainly FCN can form Lewis acid-base complexes with formaldehyde associated with cooperative H bonds, as first noticed by Wallen et al. (Blatchford, M. A.; Raveendran, P.; Wallen, S. L. J. Am. Chem. Soc. 2002, 124, 14818-14819) for CO2-philic materials under supercritical conditions. The present results, obtained with MP2(Full)/aug-cc-pVDZ calculations, show that the degeneracy of the nu(2) mode in free HCN or FCN is removed upon complexation in the same fashion as that of CO2. The splitting of these bands along with the electron structure analysis provides substantial evidence of the interaction of electron lone pairs of the carbonyl oxygen with the electron-deficient carbon atom of the cyanides. Also, this work investigates the role of H bonds acting as additional stabilizing interactions in the complexes by performing the energetic and geometric characterization.

Gas-phase basicities of polyfunctional molecules. Part 1: Theory and methods

The experimental and theoretical methods of determination of gas-phase basicities, proton affinities and protonation entropies are presented in a tutorial form. Particularities and limitations of these methods when applied to polyfunctional molecules are emphasized. Structural effects during the protonation process in the gas-phase and their consequences on the corresponding thermochemistry are reviewed and classified. The role of the nature of the basic site (protonation on non-bonded electron pairs or on pi-electron systems) and of substituent effects (electrostatic and resonance) are first examined. Then, linear correlations observed between gas-phase basicities and ionization energies or substituent constants are recalled. Hydrogen bonding plays a special part in proton transfer reactions and in the protonation characteristics of polyfunctional molecules. A survey of the main properties of intermolecular and intramolecular hydrogen bonding in both neutral and protonated species is proposed. Consequences on the protonation thermochemistry, particularly of polyfunctional molecules are discussed. Finally, chemical reactions which may potentially occur inside protonated clusters during the measurement of gas-phase basicities or inside a protonated polyfunctional molecule is examined. Examples of bond dissociations with hydride or alkyl migrations, proton transport catalysis, tautomerization, cyclization, ring opening and nucleophilic substitution are presented to illustrate the potentially complex chemistry that may accompany the protonation of polyfunctional molecules.

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.