The confinement of water molecules is vital in fields from biology to nanotechnology. The conditions allowing confinement in small finite polycyclic aromatic hydrocarbons (PAHs) are unclear, yet are crucial for understanding confinement in larger systems. Here, we report a computational study of water cluster confinement within PAHs dimers. Our results serve as a model for larger carbon allotropes and for understanding molecular interactions in confined systems. We identified size and structural motifs allowing confinement and demonstrated the motifs in various PAHs systems. We show that optimal OH⋯π interactions between water clusters and the PAH dimer permit optimal confinement to occur. However, the lack of such interactions leads to the formation of CH⋯O interactions, resulting in less ideal confinement. Confinement of layered clusters is also possible, provided that the optimal OH⋯π interactions are conserved.

The growth mechanism of hydrocarbons in ionizing environments, such as the interstellar medium (ISM), and some combustion conditions remains incompletely understood. Ab initio molecular dynamics (AIMD) simulations and molecular beam vacuum-UV (VUV) photoionization mass spectrometry experiments were performed to understand the ion-molecule growth mechanism of small acetylene clusters (up to hexamers). A dramatic dependence of product distribution on the ionization conditions is demonstrated experimentally and understood from simulations. The products change from reactive fragmentation products in a higher temperature, higher density gas regime toward a very cold collision-free cluster regime that is dominated by products whose empirical formula is (C₂H₂) _n⁺, just like ionized acetylene clusters. The fragmentation products result from reactive ion-molecule collisions in a comparatively higher pressure and temperature regime followed by unimolecular decomposition. The isolated ionized clusters display rich dynamics that contain bonded C₄H₄⁺ and C₆H₆⁺ structures solvated with one or more neutral acetylene molecules. Such species contain large amounts (>2 eV) of excess internal energy. The role of the solvent acetylene molecules is to affect the barrier crossing dynamics in the potential energy surface (PES) between (C₂H₂)_n⁺ isomers and provide evaporative cooling to dissipate the excess internal energy and stabilize products including the aromatic ring of the benzene cation. Formation of the benzene cation is demonstrated in AIMD simulations of acetylene clusters with n > 3, as well as other metastable C₆H₆⁺ isomers. These results suggest a path for aromatic ring formation in cold acetylene-rich environments such as parts of the ISM.

Polycyclic aromatic hydrocarbons (PAHs) may comprise up to 20% of the carbon budget in our galaxy and most PAHs condense onto water-rich icy grain mantles. Benzene-water clusters have been invoked as model systems for studying the photo-processing of water ice mantles containing PAHs. However, there is a paucity of information on larger aromatics, where the extended π cloud could affect photo-processing. In this study, tunable vacuum ultraviolet (VUV) photoionization of naphthalene-water clusters Nx(H2O)y (N denotes naphthalene) is performed using synchrotron radiation and analyzed by reflectron time-of-flight mass spectrometry. Naphthalene clusters up to x = 4 are generated as are naphthalene-water clusters up to y = 25. At low photon energy (<11 eV), the naphthalene moiety is ionized and there is no proton transfer from N+ to the water sub-cluster, which is very different from the benzene-water system. Protonated products, N[(H2O)xH]+ and OH radical addition product (NOH)[(H2O)xH]+ are generated above 11 eV, suggesting that water sub-clusters dominate the dynamics at high photon energies. Ab initio calculations are performed to decipher the experimental results. Energetics of the neutral structures N(H2O)1-4 and their photoionized counterparts are calculated, including ionization on the N moiety as well as on the water sub-cluster. Energy decomposition analysis (EDA) is performed to understand trends in the binding between the naphthalene and the water sub-cluster in the ionized species.

Mixed complexes of acetylene-ethylene are studied using vacuum-ultraviolet (VUV) photoionization mass spectrometry and theoretical calculations. These complexes are produced and ionized at different distances from the exit of a continuous nozzle followed by reflectron time-of-flight mass spectrometry detection. Acetylene, with a higher ionization energy (11.4 eV) than ethylene (10.6 eV), allows for tuning the VUV energy and initializing reactions either from a C2H2(+) or a C2H4(+) cation. Pure acetylene and ethylene expansions are separately carried out to compare, contrast, and hence identify products from the mixed expansion: these are C3H3(+) (m/z = 39), C4H5(+) (m/z = 53), and C5H5(+) (m/z = 65). Intensity distributions of C2H2, C2H4, their dimers and reactions products are plotted as a function of ionization distance. These distributions suggest that association mechanisms play a crucial role in product formation closer to the nozzle. Photoionization efficiency (PIE) curves of the mixed complexes demonstrate rising edges closer to both ethylene and acetylene ionization energies. We use density functional theory (ωB97X-V/aug-cc-pVTZ) to study the structures of the neutral and ionized dimers, calculate their adiabatic and vertical ionization energies, as well as the energetics of different isomers on the potential energy surface (PES). Upon ionization, vibrationally excited clusters can use the extra energy to access different isomers on the PES. At farther ionization distances from the nozzle, where the number densities are lower, unimolecular decay is expected to be the dominant mechanism. We discuss the possible decay pathways from the different isomers on the PES and examine the ones that are energetically accessible.

The growth mechanisms of organic molecules in an ionizing environment such as the interstellar medium are not completely understood. Here we examine by means of ab initio molecular dynamics (AIMD) simulations and density functional theory (DFT) computations the possibility of bond formation and molecular growth upon ionization of van der Waals clusters of pure HCN clusters, and mixed clusters of HCN and HCCH, both of which are widespread in the interstellar medium. Ionization of van der Waals clusters can potentially lead to growth in low temperature and low-density environments. Our results show, that upon ionization of the pure HCN clusters, strongly bound stable structures are formed that contain NH bonds, and growth beyond pairwise HCN molecules is seen only in a small percentage of cases. In contrast, mixed clusters, where HCCH is preferentially ionized over HCN, can grow up to 3 or 4 units long with new carbon-carbon and carbon-nitrogen covalent bonds. Moreover, cyclic molecules formed, such as the radical cation of pyridine, which is a prebiotic molecule. The results presented here are significant as they provide a feasible pathway for molecular growth of small organic molecules containing both carbon and nitrogen in cold and relatively denser environments such as in dense molecular clouds but closer to the photo-dissociation regions, and protoplanetary disks. In the mechanism we propose, first, a neutral van der Waals cluster is formed. Once the cluster is formed it can undergo photoionization which leads to chemical reactivity without any reaction barrier.

We investigated the formation mechanisms of the nucleobases adenine and guanine and the nucleobase analogues hypoxanthine, xanthine, isoguanine, and 2,6-diaminopurine in a UV-irradiated mixed 10:1 H₂O:NH₃ ice seeded with precursor purine by using ab initio and density functional theory computations. Our quantum chemical investigations suggest that a multistep reaction mechanism involving purine cation, hydroxyl and amino radicals, together with water and ammonia, explains the experimentally obtained products in an independent study. The relative abundances of these products appear to largely follow from relative thermodynamic stabilities. The key role of the purine cation is likely to be the reason why purine is not functionalized in pure ammonia ice, where cations are promptly neutralized by free electrons from NH₃ ionization. Amine group addition to purine is slightly favored over hydroxyl group attachment based on energetics, but hydroxyl is much more abundant due to higher abundance of H₂O. The amino group is preferentially attached to the 6 position, giving 6-aminopurine, that is, adenine, while the hydroxyl group is preferentially attached to the 2 position, leading to 2-hydroxypurine. A second substitution by hydroxyl or amino group occurs at either the 6 or the 2 position depending on the first substitution. Given that H₂O is far more abundant than NH₃ in the experimentally studied ices (as well as based on interstellar abundances), xanthine and isoguanine are expected to be the most abundant bi-substituted photoproducts. Key Words: Astrophysical ice-Abiotic organic synthesis-Nucleic acids-Origin of life-RNA world. Astrobiology 17, 771-785.

We investigated the formation mechanisms of the nucleobases adenine and guanine and the nucleobase analogues hypoxanthine, xanthine, isoguanine, and 2,6-diaminopurine in a UV-irradiated mixed 10:1 H₂O:NH₃ ice seeded with precursor purine by using ab initio and density functional theory computations. Our quantum chemical investigations suggest that a multistep reaction mechanism involving purine cation, hydroxyl and amino radicals, together with water and ammonia, explains the experimentally obtained products in an independent study. The relative abundances of these products appear to largely follow from relative thermodynamic stabilities. The key role of the purine cation is likely to be the reason why purine is not functionalized in pure ammonia ice, where cations are promptly neutralized by free electrons from NH₃ ionization. Amine group addition to purine is slightly favored over hydroxyl group attachment based on energetics, but hydroxyl is much more abundant due to higher abundance of H₂O. The amino group is preferentially attached to the 6 position, giving 6-aminopurine, that is, adenine, while the hydroxyl group is preferentially attached to the 2 position, leading to 2-hydroxypurine. A second substitution by hydroxyl or amino group occurs at either the 6 or the 2 position depending on the first substitution. Given that H₂O is far more abundant than NH₃ in the experimentally studied ices (as well as based on interstellar abundances), xanthine and isoguanine are expected to be the most abundant bi-substituted photoproducts. Key Words: Astrophysical ice-Abiotic organic synthesis-Nucleic acids-Origin of life-RNA world. Astrobiology 17, 771-785.

Water-cluster interactions with polycyclic aromatic hydrocarbons (PAHs) are of paramount interest in many chemical and biological processes. We report a study of anthracene monomers and dimers with water (up to four)-cluster systems utilizing molecular beam vacuum-UV photoionization mass spectrometry and density functional calculations. Structural loss in photoionization efficiency curves when adding water indicates that various isomers are generated, while theory indicates only a slight shift in energy in photoionization states of different isomers. Calculations reveal that the energetic tendency of water is to remain clustered and not to disperse around the PAH. Theoretically, we observe water confinement exclusively in the case of four water clusters and only when the anthracenes are in a cross configuration due to optimal OH⋯π interactions, indicating dependence on the size and structure of the PAH. Furthermore theory sheds light on the structural changes that occur in water upon ionization of anthracene, due to the optimal interactions of the resulting hole and water hydrogen atoms.

N-heterocyclic carbenes (NHCs) have been widely utilized for the formation of self-assembled monolayers (SAMs) on various surfaces. The main methodologies for preparation of NHCs-based SAMs either requires inert atmosphere and strong base for deprotonation of imidazolium precursors or the use of specifically-synthesized precursors such as NHC(H)[HCO₃] salts or NHC-CO₂ adducts. Herein, we demonstrate an electrochemical approach for surface-anchoring of NHCs which overcomes the need for dry environment, addition of exogenous strong base or restricting synthetic steps. In the electrochemical deposition, water reduction reaction is used to generate high concentration of hydroxide ions in proximity to a metal electrode. Imidazolium cations were deprotonated by hydroxide ions, leading to carbenes formation that self-assembled on the electrode's surface. SAMs of NO₂-functionalized NHCs and dimethyl-benzimidazole were electrochemically deposited on Au films. SAMs of NHCs were also electrochemically deposited on Pt, Pd and Ag films, demonstrating the wide metal scope of this deposition technique.

A summary of the technical advances that are incorporated in the fourth major release of the Q-Chem quantum chemistry program is provided, covering approximately the last seven years. These include developments in density functional theory methods and algorithms, nuclear magnetic resonance (NMR) property evaluation, coupled cluster and perturbation theories, methods for electronically excited and open-shell species, tools for treating extended environments, algorithms for walking on potential surfaces, analysis tools, energy and electron transfer modelling, parallel computing capabilities, and graphical user interfaces. In addition, a selection of example case studies that illustrate these capabilities is given. These include extensive benchmarks of the comparative accuracy of modern density functionals for bonded and non-bonded interactions, tests of attenuated second order Møller-Plesset (MP2) methods for intermolecular interactions, a variety of parallel performance benchmarks, and tests of the accuracy of implicit solvation models. Some specific chemical examples include calculations on the strongly correlated Cr2 dimer, exploring zeolite-catalysed ethane dehydrogenation, energy decomposition analysis of a charged ter-molecular complex arising from glycerol photoionisation, and natural transition orbitals for a Frenkel exciton state in a nine-unit model of a self-assembling nanotube.