Growth mechanism of aromatic prebiotic molecules: insights from different processes of ion-molecule reactions in benzonitrile-ammonia and benzonitrile-methylamine clusters

Benzonitrile (BN, C6H5CN) has been detected in the cold molecular cloud Taurus molecular cloud-1 (TMC-1) in 2018, which is suggested to be a precursor in the formation of more complex nitrogen-containing aromatic interstellar compounds. In this study, we utilized mass-selected infrared (IR) photodissociation spectroscopy and quantum chemical calculations to investigate the structures and gaseous ion-molecule reactions of benzonitrile-ammonia (BN-NH3) and benzonitrile-methylamine (BN-MA) clusters. The spectral observations indicate that the cyclic hydrogen bonding structure predominates in both neutral clusters. After VUV (118 nm) single-photon ionization, a new C-N covalent bond formed between BN and NH3 in the (BN-NH3)+ cluster. However, proton sharing constitutes the primary structure of the (BN-MA)+ cluster. The two nitrogen-containing interstellar molecules react with BN to yield distinct products due to difference in charge distribution and molecular polarity in the ionized clusters. The reactions of BN with other molecules contribute to our understanding of the growth mechanisms of complex interstellar molecules.

Hydrogen bond network structures of protonated dimethylamine clusters H+(DMA)n (n = 3-7)

Infrared spectroscopy of protonated dimethylamine clusters, H+(DMA)n, (n = 3-7), and their Ar-tagged clusters was performed in the NH and CH stretching vibrational region to explore their hydrogen bond network structures. A stable isomer search and vibrational spectral simulations of the clusters were also carried out to support the interpretations of the observed spectra. Weakly hydrogen-bonded NH stretching vibrational bands, which are characteristic of cyclic structures of small-sized protonated clusters, are observed in the spectra of the Ar-tagged clusters of n ≥ 5, while only linear chain type structures are suggested for the Ar-tagged clusters of n = 3-4 and the bare clusters of all the sizes. These results demonstrate that the size and temperature dependence of the hydrogen bond network structures of the protonated dimethylamine clusters is analogous to that of protonated monohydric alcohol clusters.

Infrared spectra of SinH4n-1+ ions (n = 2-8): inorganic H-(Si-H)n-1 hydride wires of penta-coordinated Si in 3c-2e and charge-inverted hydrogen bonds

SinHm+ cations are important constituents in silane plasmas and astrochemical environments. Protonated disilane (Si2H7+) was shown to have a symmetric three-centre two-electron (3c-2e) Si-H-Si bond that can also be considered as a strong ionic charge-inverted hydrogen bond with polarity Siδ+-Hδ--Siδ+. Herein, we extend our previous work to larger SinH4n-1+ cations, formally resulting from adding SiH4 molecules to a SiH3+ core. Infrared spectra of size-selected SinH4n-1+ ions (n = 2-8) produced in a cold SiH4/H2/He plasma expansion are analysed in the SiH stretch range by complementary dispersion-corrected density functional theory calculations (B3LYP-D3/aug-cc-pVTZ) to reveal their bonding characteristics and cluster growth. The ions with n = 2-4 form a linear inorganic H-(Si-H)n hydride wire with adjacent Si-H-Si 3c-2e bridges, whose strength decreases with n, as evident from their characteristic and strongly IR active SiH stretch fundamentals in the range 1850-2100 cm-1. These 3c-2e bonds result from the lowest-energy valence orbitals, and their high stability arises from their delocalization along the whole hydride wire. For SinH4n-1+ with n ≥ 5, the added SiH4 ligands form weak van der Waals bonds to the Si4H19+ chain. Significantly, because the SinH4n-1+ hydride wires are based on penta-coordinated Si atoms leading to supersaturated hydrosilane ions, analogous wires cannot be formed by isovalent carbon.

A Tug of War between the Self- and Cross-associating Hydrogen Bonds in Neutral Ammonia-Water Clusters: Energetic Insights by Molecular Tailoring Approach

In the present work, the energies of various types of individual HBs observed in neutral (NH₃ )_m (H₂ O)_n , (m+n=2 to 7) clusters were estimated using the molecular tailoring approach (MTA)-based method. The calculated individual HB energies suggest that the O-H…N HBs are the strongest (1.21 to 12.49 kcal mol^-1 ). The next ones are the O-H…O (3.97 to 9.30 kcal mol^-1 ) HBs. The strengths of N-H…N (1.09 to 5.29 kcal mol^-1 ) and N-H…O (2.85 to 5.56 kcal mol^-1 ) HBs are the weakest. The HB energies in dimers also follow this rank ordering. However, the HB energies in dimers are much smaller than those obtained by the MTA-based method due to the loss in cooperativity contribution in the dimers. Thus, the calculated cooperativity contributions, for different types of HBs, fall in the range 0.64 to 5.73 kcal mol^-1 . We wish to emphasize based on the energetic rank ordering obtained by the MTA-based method that the O-H of water is a better HB donor than the N-H of ammonia. The reasons for the observed energetic rank ordering are two folds: (i) intrinsically stronger O-H…N HBs than the O-H…O ones as revealed by dimer energies and (ii) the higher cooperativity contribution in the former than the later ones. Indeed, the MTA-based method is useful in providing the missing energetic rank ordering of various type of HBs in neutral (NH₃ )_m (H₂ O)_n clusters, in the literature.

Identifying UiO-67 Metal-Organic Framework Defects and Binding Sites through Ammonia Adsorption

Ammonia is a widely used toxic industrial chemical that can cause severe respiratory ailments. Therefore, understanding and developing materials for its efficient capture and controlled release is necessary. One such class of materials is 3D porous metal-organic frameworks (MOFs) with exceptional surface areas and robust structures, ideal for gas storage/transport applications. Herein, interactions between ammonia and UiO-67-X (X: H, NH₂ , CH₃ ) zirconium MOFs were studied under cryogenic, ultrahigh vacuum (UHV) conditions using temperature-programmed desorption mass spectrometry (TPD-MS) and in-situ temperature-programmed infrared (TP-IR) spectroscopy. Ammonia was observed to interact with μ₃ -OH groups present on the secondary building unit of UiO-67-X MOFs via hydrogen bonding. TP-IR studies revealed that under cryogenic UHV conditions, UiO-67-X MOFs are stable towards ammonia sorption. Interestingly, an increase in the intensity of the C-H stretching mode of the MOF linkers was detected upon ammonia exposure, attributed to NH-π interactions with linkers. These same binding interactions were observed in grand canonical Monte Carlo simulations. Based on TPD-MS, binding strength of ammonia to three MOFs was determined to be approximately 60 kJ mol^-1 , suggesting physisorption of ammonia to UiO-67-X. In addition, missing linker defect sites, consisting of H₂ O coordinated to Zr⁴⁺ sites, were detected through the formation of nNH₃ ⋅H₂ O clusters, characterized through in-situ IR spectroscopy. Structures consistent with these assignments were identified through density functional theory calculations. Tracking these bands through adsorption on thermally activated MOFs gave insight into the dehydroxylation process of UiO-67 MOFs. This highlights an advantage of using NH₃ for the structural analysis of MOFs and developing an understanding of interactions between ammonia and UiO-67-X zirconium MOFs, while also providing directions for the development of stable materials for efficient toxic gas sorption.

Appraisal of individual hydrogen bond strengths and cooperativity in ammonia clusters via a molecular tailoring approach

In this work, we propose and test a method, based on the molecular tailoring approach (MTA), for the evaluation of individual hydrogen bond (HB) energies in ammonia (NH₃)_n clusters. This methodology was tested, in our earlier work, on water clusters. Liquid ammonia being a universal, non-aqueous ionizing solvent, such information of individual HB strength is indispensable in many studies. The estimated HB energies by an MTA-based method, in (NH₃)_n for n = 3-8, were calculated to be in the range of 0.65 to 5.54 kcal mol^-1 with the cooperativity contribution falling between -0.54 and 1.88 kcal mol^-1 both calculated at the MP2(full)/aug-cc-pVTZ level of theory. It is seen that the strong HBs in (NH₃)_n clusters were additionally strengthened by the large contribution of HB cooperativity. The accuracy of these estimated HB energies was validated by approximately estimating the molecular energy of a given cluster by adding the sum of HB energies to the sum of monomer energies. This approximately estimated molecular energy of a given cluster was found to be in excellent agreement with the actual calculated values. The negligibly small difference (less than 5.6 kcal mol^-1) in these two values suggests that the estimated individual HB energies in ammonia clusters are quite reliable. Furthermore, these estimated HB energies by MTA are in excellent qualitative agreement with the other indirect measures of HB strength, such as HB bond distances and angles, N-H stretching frequency and the electron density values at the (3,-1) bond critical points.

Spectroscopic evidence of the C-N covalent bond formed between two interstellar molecules (ISM): acrylonitrile and ammonia

Acrylonitrile (AN) and ammonia (NH3) are two important nitrogen-containing interstellar molecules in outer space, especially on Titan. Herein, we measured infrared (IR) spectra of neutral and cationic AN-NH3 complexes by VUV single-photon ionization combined with time-of-flight mass spectrometry. On combining IR spectra with the theoretical calculations, we found that the molecules prefer to form a single-ring cyclic H-bonded structure in the neutral AN-NH3 and (AN)2-NH3 clusters. However, after ionization of AN-NH3 and (AN)2-NH3 clusters, a new C-N-covalent bond is confirmed to form directly between AN and NH3, without any energy barrier in the cationic complexes. Moreover, in the ionized (AN)2-NH3 cluster, the covalent C-N bond prefers to form between AN and NH3 rather than the two AN groups. These results provide spectroscopic evidence of AN forming a new molecule with NH3, induced by VUV radiation. The formation of the new C-N bond broadens our knowledge on the evolution of the prebiotic nitrogen-containing molecules in space.

Understanding the Hydrogen-Bonded Clusters of Ammonia (NH3) n (n = 3-6): Insights from the Electronic Structure Theory

Although it is well known that hydrogen bonds commonly exist in ammonia clusters and play an important role, there are still many challenges in understanding the electronic structure properties of hydrogen bonds. In this paper, the geometric and electronic structure properties of cyclic ammonia clusters are investigated by using first-principles density functional theory (DFT) and the Møller-Plesset perturbation theory (MP2). The calculation results show that the pentamer and hexamer have deviated from the perfect plane, while the trimer and tetramer present planarization that has been confirmed by infrared (IR) spectra. The electronic structure analysis further shows that the covalent properties play a non-negligible role in hydrogen bonding. The results also indicate that the electronic structure facilitates structure planarization. Our work not only provides insight into the role and nature of hydrogen bonds in ammonia clusters but also provides a theoretical basis for frontier science in fields such as atmospheric haze and biomolecular functions.

IR-VUV spectroscopy of pyridine dimers, trimers and pyridine-ammonia complexes in a supersonic jet

The infrared spectra of the C-H stretching vibrations of (pyridine)m, m = 1-3, and the N-H stretching vibrations of (pyridine)m-(NH3)n, m = 1, 2; n = 1-4, complexes were investigated by infrared (IR)-vacuum ultraviolet (VUV) spectroscopy under jet-cooled conditions. The ionization potential (IP0) of the pyridine monomer was determined to be 74 546 cm-1 (9.242 eV), while its complexes showed only smooth curves of the ionization thresholds at ∼9 eV, indicating large structural changes in the ionic form. The pyridine monomer exhibits five main features with several satellite bands in the C-H stretching region at 3000-3200 cm-1. Anharmonic calculations including Fermi-resonance were carried out to analyze the candidates of the overtone and combination bands which can couple to the C-H stretching fundamentals. For (pyridine)2 and (pyridine)3, most C-H bands are blue-shifted by 3-5 cm-1 from those of the monomer. The structures revealed by random searching algorithms with density functional methods indicate that the π-stacked structure is most stable for (pyridine)2, while (pyridine)3 prefers the structures stabilized by dipole-dipole and C-Hπ interactions. For the (pyridine)m-(NH3)n complexes, the mass spectrum exhibited a wide range distribution of the complexes. The observed IR spectra in the N-H stretching vibrations of the complexes showed four main bands in the 3200-3450 cm-1 region. These features are very similar to those of (NH3)n complexes, and the bands are assigned to the anti-symmetric N-H stretching band (ν3), the symmetric N-H stretching (ν1) band, and the first overtone bands of the N-H bending vibrations (2ν4). The anharmonic calculations including the Fermi-resonance between ν1 and 2ν4 well reproduced the observed spectra.

Intramolecular Reactions in Ionized Ammonia Clusters: A Direct Ab Initio Molecular Dynamics Study

Ammonia cluster cations are a chemical species that has recently attracted considerable research attention as an ion-molecule reaction species in the planetary atmosphere, surface reaction species in materials chemistry, and super-alkali species. Reactions of the radical cation of an ammonia cluster, [(NH₃)_n]⁺ (n = 2-6), following the ionization of the parent neutral cluster, were investigated using direct ab initio molecular dynamics to elucidate the reactions of the ammonia cluster cation under astrochemical conditions. The calculations showed that two competing reaction channels-proton transfer (PT) channel and complex formation channel-operate after the ionization of neutral clusters. In the PT channel, a proton of NH₃⁺ was transferred to a neighboring ammonia molecule. The PT channel was found in all clusters (n = 2-6). Reaction via the PT channel became faster with increasing cluster size and saturated around n = 5-6. In the complex formation channel, a face-to-face complex having a H₃N-NH₃⁺ structure (with a N-N bond) was formed. This channel was found only in larger clusters (n = 5-6). Time scales of PT and complex formation channels were calculated to be 20-30 and 40-50 fs, respectively. The reaction mechanism was discussed based on the results of theoretical calculations.

Large-Sized Ammonia Clusters and Solvation Energies of the Proton in Ammonia

The absolute solvation energies (free energies and enthalpies) of the proton in ammonia are used to compute the pK_a of species embedded in ammonia. They are also used to compute the solvation energies of other ions in ammonia. Despite their importance, it is not possible to determine experimentally the solvation energies of the proton in a given solvent. We propose in this work a direct approach to compute the solvation energies of the proton in ammonia from large-sized neutral and protonated ammonia clusters. To undertake this investigation, we performed a geometry optimization of neutral and protonated ammonia 30-mer, 40-mer, and 50 mer to locate stable structures. These structures have been fully optimized at both APFD/6-31++g(d,p) and M06-2X/6-31++g(d,p) levels of theory. An infrared spectroscopic study of these structures has been provided to assess the reliability of our investigation. Using these structures, we have computed the absolute solvation free energy and the absolute solvation enthalpy of the proton in ammonia. It comes out that the absolute solvation free energy of the proton in ammonia is calculated to be -1192 kJ mol^-1 , whereas the absolute solvation enthalpy is evaluated to be -1214 kJ mol^-1 . © 2019 Wiley Periodicals, Inc.

Global and local minima of protonated acetonitrile clusters

Potential energy surfaces of protonated acetonitrile clusters have been explored to locate global and local minima energy structures. The structures are stabilized by strong hydrogen bonds, anti-parallel dimers, dipole–dipole and CH⋯N interactions.

Hydrogen bond network structures of protonated short-chain alcohol clusters

Protonated alcohol clusters enable extraction of the physical essence of the nature of hydrogen bond networks.

Aqueous Solvation of Ammonia and Ammonium: Probing Hydrogen Bond Motifs with FT-IR and Soft X-ray Spectroscopy

In a multifaceted investigation combining local soft X-ray and vibrational spectroscopic probes with ab initio molecular dynamics simulations, hydrogen-bonding interactions of two key principal amine compounds in aqueous solution, ammonia (NH₃) and ammonium ion (NH₄⁺), are quantitatively assessed in terms of electronic structure, solvation structure, and dynamics. From the X-ray measurements and complementary determination of the IR-active hydrogen stretching and bending modes of NH₃ and NH₄⁺ in aqueous solution, the picture emerges of a comparatively strongly hydrogen-bonded NH₄⁺ ion via N-H donating interactions, whereas NH₃ has a strongly accepting hydrogen bond with one water molecule at the nitrogen lone pair but only weakly N-H donating hydrogen bonds. In contrast to the case of hydrogen bonding among solvent water molecules, we find that energy mismatch between occupied orbitals of both the solutes NH₃ and NH₄⁺ and the surrounding water prevents strong mixing between orbitals upon hydrogen bonding and, thus, inhibits substantial charge transfer between solute and solvent. A close inspection of the calculated unoccupied molecular orbitals, in conjunction with experimentally measured N K-edge absorption spectra, reveals the different nature of the electronic structural effects of these two key principal amine compounds imposed by hydrogen bonding to water, where a pH-dependent excitation energy appears to be an intrinsic property. These results provide a benchmark for hydrogen bonding of other nitrogen-containing acids and bases.

An ab initio anharmonic approach to study vibrational spectra of small ammonia clusters

Fermi resonance between the N–H stretching (ν₁ and ν₃) and the overtone of N–H bending (2ν₄) in ammonia has hindered the interpretation and assignments of experimental spectra of small ammonia clusters.

Structures and spectroscopy of protonated ammonia clusters at different temperatures

Protonated ammonia clusters are all Eigen structures and the first solvation shell of the related ammonium ion core is saturated by four ammonia molecules.

Structures and relative stabilities of ammonia clusters at different temperatures: DFT vs. ab initio

The global minimum energy structures of (NH₃)_n=2–10are pointed out for the first time at a given temperature.