Experimental and Theoretical Studies of the Reactivity and Thermochemistry of Dicyanamide: N(CN)2–

Acyl Phosphates as Chemically Fueled Building Blocks for Self-Sustaining Protocells

Lipids spontaneously assemble into vesicle-forming membranes. Such vesicles serve as compartments for even the simplest living systems. Vesicles have been extensively studied for constructing synthetic cells or as models for protocells-the cells hypothesized to have existed before life. These compartments exist almost always close to equilibrium. Life, however, exists out of equilibrium. In this work, we studied vesicle-based compartments regulated by a non-equilibrium chemical reaction network that converts activating agents. In this way, the compartments require a constant or periodic supply of activating agents to sustain themselves. Specifically, we use activating agents to condense carboxylates and phosphate esters into acyl phosphate-based lipids that form vesicles. These vesicles can only be sustained when condensing agents are present; without them, they decay. We demonstrate that the chemical reaction network can operate on prebiotic activating agents, opening the door to prebiotically plausible, self-sustainable protocells that compete for resources. In future work, such protocells should be endowed with a genotype, e.g., self-replicating RNA structures, to alter the protocell's behavior. Such protocells could enable Darwinian evolution in a prebiotically plausible chemical system.

Picking up Good Vibrations through Quartic Force Fields and Vibrational Perturbation Theory

Quartic force fields (QFFs) define sparse potential energy surfaces (compared to semiglobal surfaces) that are the cheapest and easiest means of computing anharmonic vibrational frequencies, especially when utilized with second-order vibrational perturbation theory (VPT2). However, flat and shallow potential surfaces are exceedingly difficult for QFFs to treat through a combination of numerical noise in the often numerically computed derivatives and in competing energy factors in the composite energies often utilized to provide high-level spectroscopic predictions. While some of these issues can be alleviated with analytic derivatives, hybrid QFFs, and intelligent choices in coordinate systems, the best practice is for predicting good molecular vibrations via QFFs is to understand what they cannot do, and this manuscript documents such cases where QFFs may fail.

Ion-molecule reactions of mass-selected ions

Gas-phase reactions of mass-selected ions with neutrals covers a very broad area of fundamental and applied mass spectrometry (MS). Oftentimes, ion-molecule reactions (IMR) can serve as a viable alternative to collision-induced dissociation and other ion dissociation techniques when using tandem MS. This review focuses on the literature pertaining applications of IMR since 2013. During the past decade considerable efforts have been made in analytical applications of IMR, including advances in one of the major techniques for characterization of unsaturated fatty acids and lipids, ozone-induced dissociation, and the development of a new technique for sequencing of large ions, hydrogen atom attachment/abstraction dissociation. Many advances have also been made in identifying gas-phase chemistry specific to a functional group in organic and biological compounds, which are useful in structure elucidation of analytes and differentiation of isomers/isobars. With "soft" ionization techniques like electrospray ionization having become mainstream for quite some time now, the efforts in the area of metal ion catalysis have firmly moved into exploring chemistry of ligated metal complexes in their "natural" oxidation states allowing to model individual steps of mechanisms in homogeneous catalysis, especially in combination with high-level DFT calculations. Finally, IMR continue to contribute to the body of knowledge in the area of chemistry of interstellar processes.

Unsupervised Reaction Pathways Search for the Oxidation of Hypergolic Ionic Liquids: 1-Ethyl-3-methylimidazolium Cyanoborohydride (EMIM+/CBH-) as a Case Study

Hypergolic ionic liquids have come under increased study for having several desirable properties as a fuel source. One particular ionic liquid, 1-ethyl-3-methylimidazolium/cyanoborohydride (EMIM⁺/CBH^-), and oxidant, nitric acid (HNO₃), has been reported to be hypergolic experimentally, but its mechanism is not well-understood at a mechanistic level. In this computational study, the reaction is first probed with ab initio molecular dynamics simulations to confirm that anion-oxidant interactions likely are the first step in the mechanism. Second, the potential energy surface of the anion-oxidant system is studied with an in-depth search over possible isomerizations, and a network of possible intermediates are found. The critical point search is unsupervised and thus has the potential of identifying structures that deviate from chemical intuition. Molecular graphs are employed for analyzing 3000+ intermediates found, and nudged elastic band calculations are employed to identify transition states between them. Finally, the reactivity of the system is discussed through examination of minimal energy paths connecting the reactant to various common products from hypergolic ionic liquid oxidation. Eight products are reported for this system: NO, N₂O, NO₂, HNO, HONO, HNO₂, HCN, and H₂O. All reaction paths leading to these exothermic products have overall reaction barriers of 6-7 kcal/mol.

Kinetics for the Reactions of H3O+(H2O)n=0-3 with Isoprene (2-Methyl-1,3-butadiene) as a Function of Temperature (300-500 K)

We report kinetics studies of H₃O⁺(H₂O)_n=0-3 with isoprene (2-methyl-1,3-butadiene, C₅H₈) as a function of temperature (300-500 K) measured using a flowing afterglow-selected ion flow tube. Results are supported by density functional (DFT) calculations at the B3LYP/def2-TZVP level. H₃O⁺ (n = 0) reacts with isoprene near the collision limit exclusively via proton transfer to form C₅H₉⁺. The first hydrate (n = 1) also reacts at the collision limit and only the proton transfer product is observed, although hydrated protonated isoprene may have been produced and dissociated thermally. Addition of a second water (n = 2) lowers the rate constant by about a factor of 10. The proton transfer of H₃O⁺(H₂O)₂ to isoprene is endothermic, but transfer of the water ligands lowers the thermicity and the likely process occurring is H₃O⁺(H₂O)₂ + C₅H₈ → C₅H₉⁺(H₂O)₂ + H₂O, followed by thermal dissociation of C₅H₉⁺(H₂O)₂. Statistical modeling indicates the amount of reactivity is consistent with the process being slightly endothermic, as is indicated by the DFT calculations. This reactivity was obscured in past experiments due to the presence of water in the reaction zone. The third hydrate is observed not to react and helps explain the past results for n = 2, as n = 2 and 3 were in equilibrium in that flow tube experiment. Very little dependence on temperature was found for the three species that did react. Finally, the C₅H₉⁺ proton transfer product further reacted with isoprene to produce mainly C₆H₉⁺ along with a small amount of clustering.

Generation and Identification of the Linear OCBNO and OBNCO Molecules with 24 Valence Electrons

Two structural isomers containing five second-row element atoms with 24 valence electrons were generated and identified by matrix-isolation IR spectroscopy and quantum chemical calculations. The OCBNO complex, which is produced by the reaction of boron atoms with mixtures of carbon monoxide and nitric oxide in solid neon, rearranges to the more stable OBNCO isomer on UV excitation. Bonding analysis indicates that the OCBNO complex is best described by the bonding interactions between a triplet-state boron cation with an electron configuration of (2s)0 (2pσ )0 (2pπ )2 and the CO/NO- ligands in the triplet state forming two degenerate electron-sharing π bonds and two ligand-to-boron dative σ bonds.

Molecular Dynamics Simulations, Reaction Pathway and Mechanism Dissection, and Kinetics Modeling of the Nitric Acid Oxidation of Dicyanamide and Dicyanoborohydride Anions

Direct dynamics simulations of HNO₃ with dicyanamide anion DCA^- (i.e., N(CN)₂^-) and dicyanoborohydride anion DCBH^- (i.e., BH₂(CN)₂^-) were performed at the B3LYP/6-31+G(d) level of theory in an attempt to elucidate the primary and secondary reactions in the two reaction systems. Guided by trajectory results, reaction coordinates and potential energy diagrams were mapped out for the oxidation of DCA^- and DCBH^- by one and two HNO₃ molecules, respectively, in the gas-phase and in the condensed-phase ionic liquids using the B3LYP/6-311++G(d,p) method. The oxidation of DCA^- by HNO₃ is initiated by proton transfer. The most important pathway leads to the formation of O₂N-NHC(O)NCN^-, and the latter reacts with a second HNO₃ to produce O₂N-NHC(O)NC(O)NH-NO₂^-(DNB^-). The oxidation of DCBH^- by HNO₃ may follow a similar mechanism as that of DCA^-, producing two analogue products: O₂N-NHC(O)BH₂CN^- and O₂N-NHC(O)BH₂C(O)NH-NO₂^-. Moreover, two new, unique reaction pathways were discovered for DCBH^- because of its boron-hydride group: (1) isomerization of DCBH^- to CNBH₂CN^- and CNBH₂NC^- and (2) H₂ elimination in which the proton in HNO₃ combines with a hydride-H in DCBH^-. The Rice-Ramsperger-Kassel-Marcus (RRKM) theory was utilized to calculate reaction kinetics and product branching ratios. The RRKM results indicate that the formation of DNB^- is exclusively important in the oxidation of DCA^-, whereas the same type of reaction is a minor channel in the oxidation of DCBH^-. In the latter case, H₂ elimination becomes dominating. The RRKM modeling also indicates that the oxidation rate constant of DCBH^- is higher than that of DCA^- by an order of magnitude. This rationalizes the enhanced preignition performance of DCBH^- over DCA^- with HNO₃.

Designing high-performance hypergolic propellants based on materials genome

A new generation of rocket propellants for deep space exploration, ionic liquid propellants, with long endurance and high stability, is attracting more and more attention. However, a major defect of ionic liquid propellants that restricts their application is the inadequate hypergolic reactivity between the fuel and the oxidant, and this defect results in local burnout and accidental explosions during the launch process. We propose a visualization model to show the features of structure, density, thermal stability, and hypergolic activity for estimating propellant performances and their application abilities. This propellant materials genome and visualization model greatly improves the efficiency and quality of developing high-performance propellants, which benefits the discovery of new advanced functional molecules in the field of energetic materials.

Electrospray ionization enters the final frontier: Mass spectrometry's role in understanding electrospray thrusters and their plumes

Electrospray thrusters using ionic liquid (IL)-based propellants are quickly gaining popularity in spacecraft design. Mass spectrometry is especially well-suited to provide important knowledge on the fundamentals of how these systems work and on evaluating their efficiencies and impacts, given that the operating principles of electrospray thrusters closely mimics the mass spectrometry experiment - in both ions are generated by electrospray and then enter a vacuum. Here, electrospray thruster technology and IL-based propellants are briefly introduced. This introduction is then followed by a discussion of mass spectrometry's current contribution to the study of IL-based electrospray thrusters - with a focus on electrospray, dissociation, and spectroscopy studies - and a brief discussion of areas ripe for immediate contributions from the mass spectrometry community.

Gas-Phase Reactions of Deprotonated Nucleobases with H, N, and O Atoms

Complex organic molecules, the hallmark of terrestrial life, are increasingly detected in exotic environments throughout the universe. Our studies probe the ion chemistry of these biomolecules. We report gas-phase reaction rate constants for five deprotonated nucleobases (adenine, cytosine, guanine, thymine, and uracil) reacting with the atomic species H, N, and O. Hydrogen atoms react at moderate rates via associative electron detachment. Oxygen atom reactions occur more rapidly, generating complex product distributions; reaction pathways include associative electron detachment, substitution of the hydrogen atom by an oxygen atom, and generation of OCN^-. Nitrogen atoms do not react with the nucleobase anions. The reaction thermodynamics were investigated computationally, and reported product channels are exothermic. Many of the proposed products have been observed in various astrochemical environments. These reactions provide insight into chemical processes that may occur at the boundaries between diffuse and dense interstellar clouds and in complex extraterrestrial ionospheres.

Computational Study of the Reaction of 1-Methyl-4-amino-1,2,4-triazolium Dicyanamide with NO2: From Reaction Dynamics to Potential Surfaces, Kinetics and Spectroscopy

Direct dynamics trajectories were calculated at the B3LYP/6-31G(d) level of theory in an attempt to understand the reaction of 1-methyl-4-amino-1,2,4-triazolium dicyanamide (MAT⁺DCA^-) with NO₂. The trajectories revealed an extensive intra-ion-pair proton transfer in MAT⁺DCA^-. The reaction pathways of the ensuing HDCA (i.e., HNCNCN) and [MAT⁺ - H_C5⁺] (i.e., deprotonated at C5-H of MAT⁺) molecules as well as DCA^- with NO₂ were identified. The reaction of NO₂ with HDCA and DCA^- produces HNC(-ONO)NCN and NCNC(-ONO)N^- or NCNCN-NO₂^-, respectively, whereas that with [MAT⁺ - H_C5⁺] results in the formation of 5-O-MAT (i.e., 4-amino-2-methyl-2,4-dihydro-3 H-1,2,4-triazo-3-one) + NO and [MAT⁺ - H₂⁺] + HNO₂. Using trajectories for guidance, structures of intermediates, transition states and products, and the corresponding reaction potential surfaces were elucidated at B3LYP/6-311++ G(d,p). Rice-Ramsperger-Kassel-Marcus (RRKM) theory was utilized to calculate the reaction rates and statistical product branching ratios. A comparison of direct dynamics simulations with RRKM modeling results indicate that the reactions of NO₂ with HDCA and DCA^- are nonstatistical. To validate our computational results, infrared and Raman spectra of MAT⁺DCA^- and its reaction products with NO₂ were calculated using an ionic liquid solvation model. The calculated spectra reproduced the vibrational frequencies detected in an earlier spectroscopic study of MAT⁺DCA^- droplets with NO₂ [ Brotton , S. J. ; J. Phys. Chem. Lett. 2017 , 8 , 6053 ].

Elucidating the Reactivity of O2 (a 1Δg): A Study with Amino Acid Anions and Related Sulfur and Oxygen Anionic Species

Rate constants and product ions were determined for a series of anions reacting with singlet molecular oxygen O₂ (a ¹Δ_g) at thermal energy using an electrospray ionization-selected ion flow tube. The 20 naturally occurring amino acids were used to produce corresponding deprotonated anions; only [Cys-H]^- and [Pro-H]^- were found to be reactive with O₂ (a ¹Δ_g), generating OSCH₂CH(NH₂)CO₂^- + HO and C₅H₆NO₂^- + H₂O₂, respectively. The reaction of O₂ (a ¹Δ_g) with [Cys-H]^- has a rate constant more than ten times larger than the reaction of O₂ (a ¹Δ_g) with [Pro-H]^-. Furthermore, reactions of O₂ (a ¹Δ_g) with carboxylic acid and thiol anions were carried out to elucidate the reactivity of the sulfur-containing functional groups. Potential energy surfaces and overall reaction exothermicities were calculated for representative reactions using density functional theory. Reactions in which attack occurs at the sulfur produce HCSO^- as an ionic product. Reactions of several carboxylic acid anions likely proceed through a hydroperoxide intermediate that is analogous to that calculated for reactions with amino acid anions at a higher collision energy. Overall, rate constants for reactions of carboxylic acid anions RC(O)O^- were found to be smaller for larger R groups.

Ignition Delay Reduction with Sodium Addition to Imidazolium-Based Dicyanamide Ionic Liquid

A range of ionic liquids (ILs) have been synthesized and modeled to better understand the role of the cation in the ignition of hypergolic ionic liquids. Vogelhuber et al. have shown by density functional theory methods that the addition of sodium cations to an ionic liquid promotes ignition with white fuming nitric acid (WFNA) by lowering energy barriers. To validate this prediction, solid sodium dicyanamide (Na⁺DCA^-) was added at various weight percents to 1-butyl-3-methylimidazolium dicyanamide (BMIM⁺DCA^-). The ignition delay was measured for each mixture with WFNA. Overall, it was found that the Na⁺DCA^- lowered the ignition delay by 11 ms at 7 wt %. The calculations done by Vogelhuber et al. appear to be consistent with this observation. The sodium cation may play a role by orienting the anion with the WFNA resulting in the favorable reaction energetics observed.

IL-oxidizer/IL-fuel combinations as greener hypergols

Two ionic liquids were ignited spontaneously upon contact, suggesting an effective route to develop greener hypergolic IL–IL bipropellants.

Theoretical Investigation of the Reactivity of Sodium Dicyanamide with Nitric Acid

There is a need to replace current hydrazine fuels with safer propellants, and dicyanamide (DCA^-)-based systems have emerged as promising alternatives because they autoignite when mixed with some oxidizers. Previous studies of the hypergolic reaction mechanism have focused on the reaction between DCA^- and the oxidizer HNO₃; here, we compare the calculated pathway of DCA^- + HNO₃ with the reaction coordinate of the ion pair sodium dicyanamide with nitric acid, Na[DCA] + HNO₃. Enthalpies and free energies are calculated in the gas phase and in solution using a quantum mechanical continuum solvation model, SMD-GIL. The barriers to the Na[DCA] + HNO₃ reaction are dramatically lowered relative to those of the reaction with the bare anion, and an exothermic exit channel to produce NaNO₃ and the reactive intermediate HDCA appears. These results suggest that Na[DCA] may accelerate the ignition reaction.

Kinetics of CO+ and CO2+ with N and O atoms

We have measured reaction rate constants for CO⁺ and CO₂⁺ reacting with N and O atoms using a selected ion flow tube apparatus equipped with a microwave discharge atom source. Experimental work was supplemented by molecular structure calculations. Calculated pathways show the sensitivity of kinetic barriers to theoretical methods and imply that high-level ab initio methods are required for accurate energetics. We report room-temperature rate constants of 1.0 ± 0.4 × 10^-11 cm³ s^-1 and 4.0 ± 1.6 × 10^-11 cm³ s^-1 for the reactions of CO⁺ with N and O atoms, respectively, and 8.0 ± 3.0 × 10^-12 cm³ s^-1 and 2.0 ± 0.8 × 10^-11 cm³ s^-1 for the reactions of CO₂⁺ with N and O atoms, respectively. The reaction of CO₂⁺ + O is observed to yield O₂⁺ exclusively. These values help resolve discrepancies in the literature and are important for modeling of the Martian atmosphere.

Reactivity of amino acid anions with nitrogen and oxygen atoms

Gas-phase reaction of deprotonated tyrosine with a ground state O atom generates five ionic products.

Reactions of Sulfur- and Oxygen-Containing Anions with Hydrogen Atoms: A Comparative Study

Reactions of hydrogen atoms with small sulfur-containing anions, SCN^-, CH₃COS^-, C₆H₅COS^-, ^-SCH₂COOH, C₆H₅S^-, 2-HOOCC₆H₄S^-, and related oxygen-containing anions, OCN^-, CH₃COO^-, C₆H₅COO^-, HOCH₂COO^-, C₆H₅O^-, 2-HOOCC₆H₄O^-, have been studied both experimentally and computationally. The experimental results show that associative electron detachment (AED) is the only channel for the reactions. The rate constants for reactions between sulfur-containing anions and H atoms are generally higher than for the related oxygen-containing anions with the exception of the reaction of SCN^-. The generally higher reactivity of the sulfur anions contrasts with previous results where AED reactivity was found to correlate with reaction exothermicity. Density functional theory calculations indicate that the reaction enthalpies, the characteristics of the reaction potential energy surfaces, and other structural and electronic factors can influence the reaction rate constants. This study indicates that organic sulfur anions can be more reactive than related oxygen anions in the interstellar medium where hydrogen atoms are abundant.

Effect of Boron Clusters on the Ignition Reaction of HNO3 and Dicynanamide-Based Ionic Liquids

Many ionic liquids containing the dicynamide anion (DCA^-, formula N(CN)₂^-) exhibit hypergolic ignition when exposed to the common oxidizer nitric acid. However, the ignition delay is often about 10 times longer than the desired 5 ms for rocket applications, so that improvements are desired. Experiments in the past decade have suggested both a mechanism for the early reaction steps and also that additives such as decaborane can reduce the ignition delay. The mechanisms for reactions of nitric acid with both DCA^- and protonated DCAH are considered here, using accurate wave function methods. Complexation of DCA^- or DCAH with borane clusters B₁₀H₁₄ or B₉H₁₄^- is found to modify these mechanisms slightly by changing the nature of some of the intermediate saddle points and by small reductions in the reaction barriers.

Triggering the Chemical Instability of an Ionic Liquid under High Pressure

Ionic liquids are an interesting class of materials due to their distinguished properties, allowing their use in an impressive range of applications, from catalysis to hypergolic fuels. However, the reactivity triggered by the application of high pressure can give rise to a new class of materials, which is not achieved under normal conditions. Here, we report on the high-pressure chemical instability of the ionic liquid 1-allyl-3-methylimidazolium dicyanamide, [allylC1im][N(CN)2], probed by both Raman and IR techniques and supported by quantum chemical calculations. Our results show a reaction occurring above 8 GPa, involving the terminal double bond of the allyl group, giving rise to an oligomeric product. The results presented herein contribute to our understanding of the stability of ionic liquids, which is of paramount interest for engineering applications. Moreover, gaining insight into this peculiar kind of reactivity could lead to the development of new or alternative synthetic routes to achieve, for example, poly(ionic liquids).