Flow-Tube Investigations of Hypergolic Reactions of a Dicyanamide Ionic Liquid Via Tunable Vacuum Ultraviolet Aerosol Mass Spectrometry

Unraveling the initial steps of the ignition chemistry of the hypergolic ionic liquid 1-ethyl-3-methylimidazolium cyanoborohydride ([EMIM][CBH]) with nitric acid (HNO3) exploiting chirped pulse triggered droplet merging

The composition of the products and the mechanistic routes for the reaction of the hypergolic ionic liquid (HIL) 1-ethyl-3-methylimidazolium cyanoborohydride ([EMIM][CBH]) and nitric acid (HNO₃) at various concentrations from 10% to 70% were explored using a contactless single droplet merging within an ultrasonic levitation setup in an inert atmosphere of argon to reveal the initial steps that cause hypergolicity. The reactions were initiated through controlled droplet-merging manipulation triggered by a frequency chirp pulse amplitude modulation. Utilizing the high-speed optical and infrared cameras surrounding the levitation process chamber, intriguing visual images were unveiled: (i) extensive gas release and (ii) temperature rises of up to 435 K in the merged droplets. The gas development was validated qualitatively and quantitatively with Fourier Transform Infrared Spectroscopy (FTIR) indicating the major gas-phase products to be hydrogen cyanide (HCN) and nitrous oxide (N₂O). The merged droplet was also probed by pulsed Raman spectroscopy which deciphered features for key functional groups of the reaction products and intermediates (-BH, -BH₂, -BH₃, -NCO); reaction kinetics revealed that the reaction was initiated by the interaction of the [CBH]^- anion of the HIL with the oxidizer (HNO₃) through proton transfer. Computations indicate the formation of a van-der-Waals complex between the [CBH]^- anion and HNO₃ initially, followed by proton transfer from the acid to the anion and subsequent extensive isomerization; these rearrangements were found to be essential for the formation of HCN and N₂O. The exoergicity observed during the merging process provides a molar enthalpy change up to 10 kJ mol^-1 to the system, which could be sufficient for a significant fraction of the reactants of about 11% to overcome the reaction barriers in the individual steps of the computationally determined minimum energy pathways.

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.

Coordination polymerization of nitrogen-rich linkers and dicyanamide anions toward energetic coordination polymers with low sensitivities

The design and synthesis of energetic materials (EMs) with high energy and reliable stabilities has attracted much attention in the field of EMs. In this work, we employed a strategy of the coordination polymerization of mild dicyanamide ions (DCA^-), two isomeric ligands 1-methyl-5-aminotetrazole (1-MAT) and 2-methyl-5-aminotetrazole (2-MAT) to construct energetic coordination polymers (ECPs). Four new ECPs {[Co(DCA)₂(1-MAT)₂]·H₂O}_n1, [Cu(DCA)₂(1-MAT)]_n2, [Cd(DCA)₂(1-MAT)₂]_n3 and [Cd(DCA)₂(2-MAT)₂]_n4 were successfully synthesized through solvent evaporation routes. Compounds 1 and 4 display 1D chains, while 2 and 3 exhibit 2D-layered structures. Compounds 1-3 with the 1-MAT ligand all exhibit reliable thermal stabilities (> 200 °C). The calculated heats of detonation (ΔH_det) of 1-3 are all higher than 1.4 kJ g^-1, which are higher than traditional explosive TNT (1.22 kJ g^-1) and the reported ECP AgMtta (HMtta = 5-methyl-1H-tetrazole, ΔH_det = 1.32 kJ g^-1). Furthermore, sensitivity testing demonstrates that 1-4 features low mechanical sensitivity to external mechanical action in contrast with the extremely sensitive azide-based ECPs [Cu₃(2-MAT)₂(N₃)₆]_n. In addition, compound 2 shows hypergolic properties via an 'oxidizer-fuel' drop experiment, demonstrating its application prospects in the field of propellants. This work details an approach of synthesizing multipurpose ECPs with reliable stabilities by introducing mild dicyanamide anions into nitrogen-rich skeletons.

Catalytic decomposition of hydroxylamine nitrate and hydrazine nitrate using Ru/ZSM-5 catalyst under mild reaction conditions

Hydroxylamine nitrate and hydrazine nitrate are dangerous explosives and toxic chemicals. Catalytic decomposition is an efficient way for disposal of these chemicals. In the current work, a Ru/ZSM-5 catalyst has been fabricated and evaluated for the decomposition of hydroxylamine nitrate and hydrazine nitrate in 1.0 mol L⁻¹ HNO₃. The hydroxylamine nitrate and hydrazine nitrate can be thoroughly decomposed under 80 °C. And the Ru/ZSM-5 catalyst can be separated from the reaction mixture and reused at least 130 times with stable catalytic performance. Easy operation, less solid waste generation, and a simple catalytic device make the strategy reported here practical, environmentally friendly, and economically attractive.

The Ru/ZSM-5 catalyst was shown to be an efficient heterogeneous catalyst for decomposition of hydroxylamine nitrate and hydrazine nitrate. The catalyst could be recovered by filtration and reused, showing good potential for industrial application.

Thermal and Catalytic Decomposition of 2-Hydroxyethylhydrazine and 2-Hydroxyethylhydrazinium Nitrate Ionic Liquid

To develop chemical kinetics models for the combustion of ionic liquid-based monopropellants, identification of the elementary steps in the thermal and catalytic decomposition of components such as 2-hydroxyethylhydrazinium nitrate (HEHN) is needed but is currently not well understood. The first decomposition step in protic ionic liquids such as HEHN is typically the proton transfer from the cation to the anion, resulting in the formation of 2-hydroxyethylhydrazine (HEH) and HNO₃. In the first part of this investigation, the high-temperature thermal decomposition of HEH is probed with flash pyrolysis (<1400 K) and vacuum ultraviolet (10.45 eV) photoionization time-of-flight mass spectrometry (VUV-PI-TOFMS). Next, the investigation into the thermal and catalytic decomposition of HEHN includes two mass spectrometric techniques: (1) tunable VUV-PI-TOFMS (7.4-15 eV) and (2) ambient ionization mass spectrometry utilizing both plasma and laser ionization techniques whereby HEHN is introduced onto a heated inert or iridium catalytic surface and the products are probed. The products can be identified by their masses, their ionization energies, and their collision-induced fragmentation patterns. Formation of product species indicates that catalytic surface recombination is an important reaction process in the decomposition mechanism of HEHN. The products and their possible elementary reaction mechanisms are discussed.

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₃.

Controlled Chemistry via Contactless Manipulation and Merging of Droplets in an Acoustic Levitator

A unique, versatile, and material-independent approach to manipulate contactlessly and merge two chemically distinct droplets suspended in an acoustic levitator is reported. Large-amplitude axial oscillations are induced in the top droplet by low-frequency amplitude modulation of the ultrasonic carrier wave, which causes the top sample to merge with the sample in the pressure minimum below. The levitator is enclosed within a pressure-compatible process chamber to enable control of the environmental conditions. The merging technique permits precise control of the substances affecting the chemical reactions, the sample temperature, the volumes of the liquid reactants down to the picoliter range, and the mixing locations in space and time. The performance of this approach is demonstrated by merging droplets of water (H₂O) and ethanol (C₂H₅OH), conducting an acid-base reaction between aqueous droplets of sodium hydroxycarbonate (NaHCO₃) and acetic acid (CH₃COOH), the hypergolic explosion produced via merging a droplet of an ionic liquid with nitric acid (HNO₃), and the coalescence of a solid particle (CuSO₄·5H₂O) and a water droplet followed by dehydration using a carbon dioxide laser. The physical and chemical changes produced by the merging are traced in real time via complementary Raman, Fourier-transform infrared, and ultraviolet-visible spectroscopies. The concept of the contactless manipulation of liquid droplets and solid particles may fundamentally change how scientists control and study chemical reactions relevant to, for example, combustion systems, material sciences, medicinal chemistry, planetary sciences, and biochemistry.

Molecular Dynamics Simulations and Product Vibrational Spectral Analysis for the Reactions of NO2 with 1-Ethyl-3-methylimidazolium Dicyanamide (EMIM+DCA-), 1-Butyl-3-methylimidazolium Dicyanamide (BMIM+DCA-), and 1-Allyl-3-methylimidazolium Dicyanamide (AMIM+DCA-)

Direct dynamics trajectory simulations were carried out for the NO₂ oxidation of 1-ethyl-3-methylimidazolium dicyanamide (EMIM⁺DCA^-), which were aimed at probing the nature of the primary and secondary reactions in the system. Guided by trajectory results, reaction coordinates and potential energy diagrams were mapped out for NO₂ with EMIM⁺DCA^-, as well as with its analogues 1-butyl-3-methylimidazolium dicyanamide (BMIM⁺DCA^-) and 1-allyl-3-methylimidazolium dicyanamide (AMIM⁺DCA^-). Reactions of the dialkylimidazolium-dicyanamide (DCA) ionic liquids (ILs) are all initiated by proton transfer and/or alkyl abstraction between 1,3-dialkylimidazolium cations and DCA^- anion, of which two exoergic pathways are particularly relevant to their oxidation activities. One pathway is the transfer of a H^β-proton from the ethyl, butyl, or allyl group of the dialkylimidazolium cation to DCA^- that results in the concomitant elimination of the corresponding alkyl as a neutral alkene, and the other pathway is the alkyl abstraction by DCA^- via a second order nucleophilic substitution (SN₂) mechanism. The intra-ion-pair reaction products, including [dialkylimidazolium⁺ - H_C2⁺], alkylimidazole, alkene, alkyl-DCA, HDCA, and DCA^-, react with NO₂ and favor the formation of nitrite (-ONO) complexes over nitro (-NO₂) complexes, albeit the two complex structures have similar formation energies. The exoergic intra-ion-pair reactions in the dialkylimidazolium-DCA ILs account for their significantly higher oxidation activities over the previously reported 1-methyl-4-amino-1,2,4-triazolium dicyanamide [Liu, J.; J. Phys. Chem. B 2019, 123, 2956-2970] and for the relatively higher reactivity of BMIM⁺DCA^- vs AMIM⁺DCA^- as BMIM⁺ has a higher reaction path degeneracy for intra-ion-pair H^β-proton transfer and its H^β-transfer is more energetically favorable. To validate and directly compare our computational results with spectral measurements in the ILs, infrared and Raman spectra of BMIM⁺DCA^- and AMIM⁺DCA^- and their products with NO₂ were calculated using an ionic liquid solvation model. The simulated spectra reproduced all of the vibrational frequencies detected in the reactions of BMIM⁺DCA^- and AMIM⁺DCA^- IL droplets with NO₂ (as reported by Brotton et al. [ J. Phys. Chem. A 2018, 122, 7351-7377] and Lucas et al. [ J. Phys. Chem. A 2019, 123, 400-416]).

Thermal Decomposition and Hypergolic Reaction of a Dicyanoborohydride Ionic Liquid

In this study, in situ infrared spectroscopy techniques and thermogravimetric analysis coupled with mass spectrometry (TGA-MS) are employed to characterize the reactivity of the ionic liquid, 1-butyl-3-methylimidazolium dicyanoborohydride (BMIM⁺DCBH^-), in comparison to the well-characterized 1-butyl-3-methylimidazolium dicyanamide (BMIM⁺DCA^-) ionic liquid. TGA measurements determined the enthalpy of vaporization (ΔH_vap) to be 112.7 ± 12.3 kJ/mol at 298 K. A rapid scan Fourier transform infrared spectrometer was used to obtain vibrational information useful in tracking the appearance and disappearance of species in the hypergolic reactions of BMIM⁺DCBH^- and BMIM⁺DCA^- with white fuming nitric acid (WFNA) and in the thermal decomposition of these energetic ionic liquids. Attenuated total reflectance measurements recorded the infrared spectra of the reactant sample (BMIM⁺DCBH^-) and the liquid reaction products after reacting with WFNA. Computational chemistry efforts, aided by the experimental results, were used to propose key reaction pathways leading to the hypergolic ignition of BMIM⁺DCBH^- + WFNA. Experimental results indicate that the hypergolic reaction of BMIM⁺DCBH^- with WFNA generates both common and unique intermediates as compared to previous BMIM⁺DCA^- + WFNA investigations: nitrous oxide was generated during both hypergolic reactions indicating that it may play a crucial role in the hypergolic ignition process, NO₂ was generated in significantly higher concentrations for BMIM⁺DCBH^- than for BMIM⁺DCA^-, CO₂ was only generated for BMIM⁺DCA^-, and HCN was only generated during thermal decomposition and hypergolic ignition of BMIM⁺DCBH^-.

First-Principles Molecular Dynamics of Monomethylhydrazine and Nitrogen Dioxide

The exploration of chemical reactions preceding ignition is essential for the development of ideal hypergolic propellants. Unexpected reaction pathways of a hypergolic mixture composed of monomethylhydrazine and nitrogen dioxide are predicted through a cooperative combination of (i) spin-unrestricted ab initio molecular dynamics (AIMD) and (ii) wave packet dynamics of protons. Ensembles of AIMD trajectories reveal a sequence of reaction steps for proton transfer and rupture of the C-N bond. The details of proton transfer are explored by wave packet dynamics on the basis of ab initio potential energy surfaces from AIMD trajectories. The possibility of spontaneous ignition of this hypergolic mixture at room temperature is predicted as a quantized feature of proton-transfer dynamics.

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 ].

Oxidation of a Levitated Droplet of 1-Allyl-3-methylimidazolium Dicyanamide by Nitrogen Dioxide

Understanding the reaction mechanisms of ionic liquids and their oxidizers is necessary to develop the next generation of hypergolic, ionic-liquid-based fuels. We studied reactions between a levitated droplet of 1-allyl-3-methylimidazolium dicyanamide ([AMIM][DCA]), with and without hydrogen-capped boron nanoparticles, and nitrogen dioxide (NO₂). The reactions were monitored with Fourier-transform infrared (FTIR) and Raman spectroscopy. The emergence of new structures in the FTIR and Raman spectra is consistent with the formation of functional groups including organic nitrites (RONO), nitroamines (R¹R²NNO₂), and carbonitrates (R¹R²C=NO₂^-). Possible reaction mechanisms based on these new functional groups are discussed. The reaction rates were deduced at various temperatures by heating the levitated droplets with a carbon dioxide laser. We thereby determined an overall activation energy of 38.5 ± 2.3 kJ mol^-1 for the oxidation of [AMIM][DCA] for the first time.

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.

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.

Spectroscopic Investigation of the Primary Reaction Intermediates in the Oxidation of Levitated Droplets of Energetic Ionic Liquids

The production of the next generation of hypergolic, ionic-liquid-based fuels requires an understanding of the reaction mechanisms between the ionic liquid and oxidizer. We probed reactions between a levitated droplet of 1-methyl-4-amino-1,2,4-triazolium dicyanamide ([MAT][DCA]), with and without hydrogen-capped boron nanoparticles, and the nitrogen dioxide (NO₂) oxidizer. The apparatus exploits an ultrasonic levitator enclosed within a pressure-compatible process chamber equipped with complementary Raman, ultraviolet-visible, and Fourier-transform infrared (FTIR) spectroscopic probes. Vibrational modes were first assigned to the FTIR and Raman spectra of droplets levitated in argon. Spectra were subsequently collected for pure and boron-doped [MAT][DCA] exposed to nitrogen dioxide. By comparison with electronic structure calculations, some of the newly formed modes suggest that the N atom of the NO₂ molecule bonds to a terminal N on the dicyanamide anion yielding [O₂N-NCNCN]^-. This represents the first spectroscopic evidence of a key reaction intermediate in the oxidation of levitated ionic liquid droplets.

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.

Catalytic Decomposition of Hydroxylammonium Nitrate Ionic Liquid: Enhancement of NO Formation

Hydroxylammonium nitrate (HAN) is a promising candidate to replace highly toxic hydrazine in monopropellant thruster space applications. The reactivity of HAN aerosols on heated copper and iridium targets was investigated using tunable vacuum ultraviolet photoionization time-of-flight aerosol mass spectrometry. The reaction products were identified by their mass-to-charge ratios and their ionization energies. Products include NH₃, H₂O, NO, hydroxylamine (HA), HNO₃, and a small amount of NO₂ at high temperature. No N₂O was detected under these experimental conditions, despite the fact that N₂O is one of the expected products according to the generally accepted thermal decomposition mechanism of HAN. Upon introduction of iridium catalyst, a significant enhancement of the NO/HA ratio was observed. This observation indicates that the formation of NO via decomposition of HA is an important pathway in the catalytic decomposition of HAN.