Thermophysical properties of energetic ionic liquids/nitric acid mixtures: insights from molecular dynamics simulations

Structural Properties of HEHN- and HAN-Based Ionic Liquid Mixtures: A Polarizable Molecular Dynamics Study

Molecular dynamics simulations of binary mixtures comprising 2-hydroxyethylhydrazinium nitrate (HEHN) and hydroxylammonium nitrate (HAN) were conducted using the polarizable APPLE&P force field to investigate fundamental properties of multimode propulsion (MMP) propellants. Calculated densities as a function of temperature were in good agreement with experiments and similar simulations. The structural properties of neat HEHN and HAN-HEHN provided insights into their inherent, protic nature. Radial distribution functions (RDFs) identified key hydrogen bonding sites located at N-H···O and O-H···O within a first solvation shell of approximately 2 Å. Angular distribution functions further affirmed the relatively strong nature of the hydrogen bonds with nearly linear directionality. The increased hydroxylammonium cation (HA+) mole fraction shows the influence of competitively strong hydrogen bonds on the overall hydrogen bond network. Dominant spatial motifs via three-dimensional distribution functions along with nearly nanosecond-long hydrogen bond lifetimes highlight the local bonding environment that may precede proton transfer reactions.

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

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.