Investigation of the Proton-Bound Dimer of Dihydrogen Phosphate and Formate Using Infrared Spectroscopy in Helium Droplets.

Understanding the structural and dynamic properties of proton-bound complexes is crucial for elucidating fundamental aspects of chemical reactivity and molecular interactions. In this work, the proton-bound complex between dihydrogen phosphate and formate, and its deuterated counterparts, is investigated using IR action spectroscopy in helium droplets. Contrary to the initial expectation that the stronger phosphoric acid would donate a proton to formate, both experiment and theory show that all exchangeable protons are located in the phosphate moiety. The experimental spectra show good agreement with both scaled harmonic and VPT2 anharmonic calculations, indicating that anharmonic effects are small. Some H-bending modes of the nondeuterated complex are found to be sensitive to the helium environment. In the case of the partially deuterated complexes, the experiments indicate that internal dynamics leads to isomeric interconversion upon IR excitation.

Water confinement in small polycyclic aromatic hydrocarbons.

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 influence of adsorption geometry on the reduction affinity of nitroaromatics on Au(111).

Chemoselective reduction of nitro groups in multifunctional nitroaromatics is a challenging catalytic process with high interest due to the importance of the resulting anilines for the chemical industry. Molecular-level understanding of the ways by which adsorption geometry of nitroaromatics influence their affinity toward nitro reduction will enable the development of highly selective reactions. Herein, taking advantage of the well-ordered self-assembly of para- and ortho-nitrothiophenol (p-NTP and o-NTP, respectively) monolayers on Au(111), we examined the correlation between adsorption geometry and nitro reduction affinity. The anchoring geometry of NTPs and their nitro reduction affinity were determined by conducting polarized X-ray absorption spectroscopy while the influence of NTPs' adsorption geometry on the interaction with the Au surface was analyzed by density functional theory (DFT) calculations. Exposure of surface anchored p-NTPs to reducing conditions led to their reorientation from a tilt angle of 52° to 25°, which enabled strong interactions between the π system of the molecules and the Au surface. Direct correlation was identified between the surface proximity of the nitro group, its parallel position to the surface and the resulting reduction yield. The asymmetric structure of o-NTP led to a tilted adsorption geometry in which the nitro group was rotated away from the plane of the aromatic ring and therefore was positioned parallel and in high proximity to the Au surface. This positioning led to surface-bonding that involved the oxygen atoms of o-NTP. The higher surface proximity and stronger surface interactions of the nitro group in o-NTP enabled nitro reduction already at 180 °C, while in p-NTP nitro reduction was achieved only at 230 °C, due to the longer distance between the NO2 group and the Au surface that led to weaker adsorbate-surface interactions. Thus, parallel positioning of the nitro group and high surface proximity were found as essential descriptors for nitro reduction affinity in both p-NTP and o-NTP on the Au surface. These findings provide explicit guidelines for tuning the reactant and surface properties in order to control the reactant's adsorption geometry for selective nitro reduction in multifunctional nitroaromatics.

A combined theoretical and experimental study of small anthracene-water clusters.

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.

Accelerated Multiphosphorylated Peptide Synthesis.

Preparing phosphorylated peptides with multiple adjacent phosphorylations is synthetically difficult, leads to ß-elimination, results in low yields, and is extremely slow. We combined synthetic chemical methodologies with computational studies and engineering approaches to develop a strategy that takes advantage of fast stirring, high temperature, and a very low concentration of 1,8-diazabicyclo[5.4.0]undec-7-ene (DBU) to produce multiphosphorylated peptides at an extremely rapid time and high purity.

Isomerization of hydrogen cyanide and hydrogen isocyanide in a cluster environment: quantum chemical study.

Hydrogen cyanide (HCN) and its isomer hydrogen isocyanide (HNC) are omnipresent in the interstellar medium (ISM). The ratio between the two isomers serves as an indicator of the physical conditions in different areas of the ISM. As such, the isomerization process between the two isomers has been extensively studied on the neutral potential energy surface. Moreover, HCN and HNC are thought to be precursors of important organic molecules, such as adenine. Here, we use quantum chemistry calculations and ab initio molecular dynamics simulations to focus on the chemistry that occurs upon ionization of pure HNC clusters. We demonstrate that upon ionization of HNC clusters, a distonic ion CN⋯HCNH+ is formed, and this formation is accompanied by HNC-to-HCN isomerization. Moreover, we show that the cluster environment and the network of hydrogen bonds are crucial for the isomerization process to occur and for the stabilization of the clusters. We demonstrate that, in contrast to HNC clusters, isomerization of ionized HCN clusters can occur only for the larger clusters. In addition, we discuss the formation of aminonitrile cation in the clusters and propose a barrierless route for diaminonitrile, a known precursor of amino acids and nucleobases, to form.

The Effect of Cluster Size on the Intra-Cluster Ionic Polymerization Process.

Polyaromatic hydrocarbons (PAHs) are widespread in the interstellar medium (ISM). The abundance and relevance of PAHs call for a clear understanding of their formation mechanisms, which, to date, have not been completely deciphered. Of particular interest is the formation of benzene, the basic building block of PAHs. It has been shown that the ionization of neutral clusters can lead to an intra-cluster ionic polymerization process that results in molecular growth. Ab-initio molecular dynamics (AIMD) studies in clusters consisting of 3-6 units of acetylene modeling ionization events under ISM conditions have shown maximum aggregation of three acetylene molecules forming bonded C6H6+ species; the larger the number of acetylene molecules, the higher the production of C6H6+. These results lead to the question of whether clusters larger than those studied thus far promote aggregation beyond three acetylene units and whether larger clusters can result in higher C6H6+ production. In this study, we report results from AIMD simulations modeling the ionization of 10 and 20 acetylene clusters. The simulations show aggregation of up to four acetylene units producing bonded C8H8+. Interestingly, C8H8+ bicyclic species were identified, setting a precedent for their astrochemical identification. Comparable reactivity rates were shown with 10 and 20 acetylene clusters.

Molecular dynamics reveals formation path of benzonitrile and other molecules in conditions relevant to the interstellar medium.

Polycyclic aromatic hydrocarbons and polycyclic aromatic nitrogen heterocycles are believed to be widespread in different areas of the interstellar medium. However, the astronomical detection of specific aromatic molecules is extremely challenging. As a result, only a few aromatic molecules have been identified, and very little is known about how they are formed in different areas of the interstellar medium. Recently, McGuire et al. [Science 359, 202-205 (2018)] detected the simple aromatic molecule benzonitrile in Taurus Molecular Cloud-1. Although benzonitrile has been observed, the molecular mechanism for its formation is still unknown. In this study, we use quantum chemistry and ab initio molecular dynamics to model ionization processes of van der Waals clusters containing cyanoacetylene and acetylene molecules. We demonstrate computationally that the clusters' ionization leads to molecular formation. For pure cyanoacetylene clusters, we observe bond formation among two and three monomer units, whereas in mixed clusters, bond formation can also occur in up to four units. We show that the large amount of energy available to the system after ionization can lead to barrier crossing and the formation of complex molecules. Our study reveals the rich chemistry that is observed upon ionization of the clusters, with a wide variety of molecules being formed. Benzonitrile is among the observed molecules, and we study the potential energy path for its formation. These results also offer insights that can guide astronomers in their search for aromatic molecules in the interstellar medium.

Stereoselective Diels-Alder Reactions of gem-Diborylalkenes: Toward the Synthesis of gem-Diboron-Based Polymers.

Although gem-diborylalkenes are known to be among the most valuable reagents in modern organic synthesis, providing a rapid access to a wide array of transformations, including the construction of C-C and C-heteroatom bonds, their use as dienophile-reactive groups has been rare. Herein we report the Diels-Alder (DA) reaction of (unsymmetrical) gem-diborylalkenes. These reactions provide a general and efficient method for the stereoselective conversion of gem-diborylalkenes to rapidly access 1,1-bisborylcyclohexenes. Using the same DA reaction manifold with borylated-dienes and gem-diborylalkenes, we also developed a concise, highly regioselective synthesis of 1,1,2-tris- and 1,1,3,4-tetrakis(boronates)cyclohexenes, a family of compounds that currently lack efficient synthetic access. Furthermore, DFT calculations provided insight into the underlying factors that control the chemo-, regio-, and stereoselectivity of these DA reactions. This method also provides stereodivergent syntheses of gem-diborylnorbornenes. The utility of the gem-diborylnorbornene building blocks was demonstrated by ring-opening metathesis polymerization (ROMP), providing a highly modular approach to the first synthesis of the gem-diboron-based polymers. Additionally, these polymers have been successfully submitted to postpolymerization modification reactions. Given its simplicity and versatility, we believe that this novel DA and ROMP approach holds great promise for organoboron synthesis as well as organoboron-based polymers and that it will result in more novel transformations in both academic and industrial research.

Electrochemical deposition of N-heterocyclic carbene monolayers on metal surfaces.

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)[HCO3] salts or NHC-CO2 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 NO2-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.

Molecular growth upon ionization of van der Waals clusters containing HCCH and HCN is a pathway to prebiotic molecules.

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.

Photoredox-Mediated Reaction of gem-Diborylalkenes: Reactivity Toward Diverse 1,1-Bisborylalkanes.

The use of gem-diborylalkenes as radical-reactive groups is explored for the first time. These reactions provide an efficient and general method for the photochemical conversion of gem-diborylalkenes to rapidly access 1,1-bisborylalkanes. This method exploits a novel photoredox decarboxylative radical addition to gem-diborylalkenes to afford α-gem-diboryl carbon-centered radicals, which benefit from additional stability by virtue of an interaction with the empty p-orbitals on borons. The reaction offers a highly modular and regioselective approach to γ-amino gem-diborylalkanes. Furthermore, EPR spectroscopy and DFT calculations have provided insight into the radical mechanism underlying the photochemistry reaction and the stability of the bis-metalated radicals, respectively.

Probing solvation and reactivity in ionized polycyclic aromatic hydrocarbon-water clusters with photoionization mass spectrometry and electronic structure calculations.

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.

Mechanisms of the Formation of Adenine, Guanine, and Their Analogues in UV-Irradiated Mixed NH3:H2O Molecular Ices Containing Purine.

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 H2O:NH3 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 NH3 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 H2O. 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 H2O is far more abundant than NH3 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.

Ab initio dynamics and photoionization mass spectrometry reveal ion-molecule pathways from ionized acetylene clusters to benzene cation.

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 (C2H2) 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 C4H4+ and C6H6+ 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 (C2H2)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 C6H6+ isomers. These results suggest a path for aromatic ring formation in cold acetylene-rich environments such as parts of the ISM.

Probing Ionic Complexes of Ethylene and Acetylene with Vacuum-Ultraviolet Radiation.

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.

Seniority zero pair coupled cluster doubles theory.

Coupled cluster theory with single and double excitations accurately describes weak electron correlation but is known to fail in cases of strong static correlation. Fascinatingly, however, pair coupled cluster doubles (p-CCD), a simplified version of the theory limited to pair excitations that preserve the seniority of the reference determinant (i.e., the number of unpaired electrons), has mean field computational cost and is an excellent approximation to the full configuration interaction (FCI) of the paired space provided that the orbital basis defining the pairing scheme is adequately optimized. In previous work, we have shown that optimization of the pairing scheme in the seniority zero FCI leads to a very accurate description of static correlation. The same conclusion extends to p-CCD if the orbitals are optimized to make the p-CCD energy stationary. We here demonstrate these results with numerous examples. We also explore the contributions of different seniority sectors to the coupled cluster doubles (CCD) correlation energy using different orbital bases. We consider both Hartree-Fock and Brueckner orbitals, and the role of orbital localization. We show how one can pair the orbitals so that the role of the Brueckner orbitals at the CCD level is retained at the p-CCD level. Moreover, we explore ways of extending CCD to accurately describe strongly correlated systems.

Stability of hemi-bonded vs proton-transferred structures of (H2O)2(+), (H2S)2(+), and (H2Se)2(+) studied with projected Hartree-Fock methods.

Hartree-Fock (HF) is known to suffer from drawbacks in the description of the relative stabilities between the hemi-bonded (HB) and proton-transferred (PT) isomers of the water dimer cation, (H2O)2(+). The energy difference predicted by HF is too large, approximately 27 kcal/mol, which is lowered to 7 kcal/mol when correlation effects are added. The error in HF has been previously attributed to the large dynamic correlation effects in the HB structure as well to the large symmetry breaking this structure exhibits. In this study we use the recently developed projected Hartree-Fock (PHF) methods to study the relative stability of the two isomers of (H2O)2(+) as well as its second and third row analogs, namely, (H2S)2(+) and (H2Se)2(+). In PHF, symmetries are broken and restored in a variation-after-projection approach and thus can deal easily with systems for which HF itself spontaneously breaks symmetry. We use different flavors of PHF (SUHF, KSUHF, SGHF, and KSGHF) to explore their ability in capturing dynamic correlation effects and to compare their performance to different wave function based methods. We study the role of the symmetry-breaking in the above systems, using wave function based methods with unrestricted and restricted wave functions as well as performing a single-shot symmetry restoration (a projection-after-variation scheme).

Seniority-based coupled cluster theory.

Doubly occupied configuration interaction (DOCI) with optimized orbitals often accurately describes strong correlations while working in a Hilbert space much smaller than that needed for full configuration interaction. However, the scaling of such calculations remains combinatorial with system size. Pair coupled cluster doubles (pCCD) is very successful in reproducing DOCI energetically, but can do so with low polynomial scaling (N(3), disregarding the two-electron integral transformation from atomic to molecular orbitals). We show here several examples illustrating the success of pCCD in reproducing both the DOCI energy and wave function and show how this success frequently comes about. What DOCI and pCCD lack are an effective treatment of dynamic correlations, which we here add by including higher-seniority cluster amplitudes which are excluded from pCCD. This frozen pair coupled cluster approach is comparable in cost to traditional closed-shell coupled cluster methods with results that are competitive for weakly correlated systems and often superior for the description of strongly correlated systems.

Calculation of transition dipole moment in fluorescent proteins--towards efficient energy transfer.

Förster Resonance Energy Transfer (FRET) between fluorescent proteins (FPs) is widely used to construct fluorescent sensor proteins, to study intracellular protein-protein interactions and to monitor conformational changes in multidomain proteins. Although FRET depends strongly on the orientation of the transition dipole moments (TDMs) of the donor and acceptor fluorophores, this orientation dependence is currently not taken into account in FRET sensor design. Similarly, studies that use FRET to derive structural constrains typically assume a κ(2) of 2/3 or use the TDM of green fluorescent protein, as this is the only FP for which the TDM has been determined experimentally. Here we used time-dependent density functional theory (TD-DFT) methods to calculate the TDM for a comprehensive list of commonly used fluorescent proteins. The method was validated against higher levels of calculation. Validation with model compounds and the experimentally determined TDM of GFP shows that the TDM is mostly determined by the structure of the π-conjugated fluorophore and is insensitive to non-conjugated side chains or the protein surrounding. Our calculations not only provide TDM for most of the currently used FPs, but also suggest an empirical rule that can be used to obtain the TDMs for newly developed fluorescent proteins in the future.