Effect of alkyl chain length in protic ionic liquids: an AIMD perspective

Crystallization Kinetics of Phosphonium Ionic Liquids: Effect of Cation Alkyl Chain Length and Thermal History

The effects of alkyl chain length on the crystallization kinetics and ion mobility of tetraalkylphosphonium, [P666,n][TFSI], (n = 2, 6, 8, and 12) ionic liquids were studied by differential scanning calorimetry (DSC) and broadband dielectric spectroscopy (BDS) over a wide temperature range. The liquid-glass transition temperature (Tg) and ion dynamics examined over a broad T range were almost insensitive to structural modifications of the phosphonium cation. In contrast, the crystallization kinetics were strongly affected by the length of the fourth alkyl chain. Furthermore, the thermal history of the sample (cold vs melt crystallization) significantly impacted the crystallization rate. It has been found that the nature of crystallization phenomena is the same across the homologous series, while the kinetic aspect differs. Finally, electric conductivity in supercooled liquid and crystalline solid phases was measured for all samples, revealing significant ionic conductivity, largely independent of the cation structure.

Bulk Nanostructure of Mixtures of Choline Arginate, Choline Lysinate, and Water

Neutron diffraction with empirical potential structure refinement was used to investigate the bulk liquid nanostructure of mixtures of choline arginate (Ch[Arg]), choline lysinate (Ch[Lys]), and water at mole ratios of 1Ch[Arg]:1Ch[Lys]:6H2O (balanced), 1Ch[Arg]:1Ch[Lys]:20H2O (balanced dilute), 3Ch[Arg]:1Ch[Lys]:12H2O (Arg- rich), and 1Ch[Arg]:3Ch[Lys]:12H2O (Lys- rich). The Arg- and Lys- anions tend not to associate due to electrostatic repulsion between charge groups and weak anion-anion attractions. This means that the local ion structures around the anions in these mixtures resemble the parent single-component systems. The bulk liquid nanostructure varies with the Arg-:Lys- ratio. In the Lys--rich mixture (1Ch[Arg]:3Ch[Lys]:12H2O), Lys- side chains cluster into a continuous apolar domain separated from a charged domain of polar groups. In the balanced mixture (1Ch[Arg]:1Ch[Lys]:6H2O), Lys- side chains form discrete apolar aggregates within a continuous polar domain of Arg-, Ch+, and water, and in the Arg--rich mixture (3Ch[Arg]:1Ch[Lys]:12H2O), the distribution of Lys- and Arg- is nearly homogeneous. Finally, in the balance dilute system (1Ch[Arg]:1Ch[Lys]:20H2O), a percolating water domain forms.

Liquid Structure and Hydrogen Bonding in Aqueous Hydroxylammonium Nitrate

Hydroxylammonium nitrate (HAN) has emerged as a promising component in ionic liquid-based spacecraft propellants. However, the physicochemical and structural properties of aqueous HAN have been largely overlooked. The purpose of this study is to investigate the hydrogen bonding in aqueous HAN and understand its implications on these properties and the proton transfer mechanism as a function of concentration. Classical polarizable molecular dynamics simulations have been employed with the APPLE&P force field to analyze the geometry of individual hydrogen bonds and the overall hydrogen-bonding network in various concentrations of aqueous HAN. Radial distribution functions (RDFs) and spatial distribution functions (SDFs) indicate the structural arrangement of the species and their hydrogen bonds. Projections of water density and the orientation of its electric dipole moment near the ions provide insight into the hydrogen-bonding network. The incorporation of water into the hydrogen-bonding network at high ion concentrations occurs via interstitial accommodation around the ions immediately outside the first solvation shell. While ion pairs are observed at all concentrations considered, the frequency of Ha···On hydrogen bonds increases substantially with the ion concentration. The findings contribute to a better fundamental understanding of HAN and the precursors of reactivity, crucial to the development of "green" spacecraft propellants.

Molecular simulations of understanding the Na+ ion structure, dynamic and thermodynamic behavior in ionic liquids: Butyl ammonium hydrogen bisulfate and tri-butyl ammonium hydrogen bisulfate

This manuscript presents the all-atom molecular dynamics simulations to investigate intermolecular structure and solvation thermodynamics of Na+ ion in two different ammonium-based protic ionic liquids (1) Butyl Ammonium hydrogen bisulfate [BA+][HSO4-], (2) Tri-butyl ammonium hydrogen bisulfate [TBA+][HSO4-]. The ionic liquid [BA+][HSO4-] show a more coordinated behavior when compared to [TBA+][HSO4-], which is observed over the temperature range from 278 K to 348 K. Hydrogens of the cations show a hydrogen bonding interaction with oxygens of anions. The cationic [TBA+] molecules show more solvation behavior with anions when compared to the [BA+]. The Na+ ion show a strong coordination structure with [HSO4-] in [TBA+][HSO4-] when compared to the [BA+][HSO4-]. We further calculate the detailed solvation free energy (ΔG) calculations using thermodynamic integration. We found that the ΔG of Na+ is more favorable in [TBA+][HSO4-] when compared to [BA+][HSO4-] in the temperature range varying from 278 K to 348 K. With the temperature rise, we observe the more favorable solvation of Na+ in both ionic liquids. On the other hand, the solvation of Cl- becomes less favorable. Overall, this manuscript provides detailed molecular level structural and thermodynamic origins of Na+ in protic ionic liquids useful for designing and developing sustainable electrolytes for Na+ battery applications.

Understanding X-ray Photoelectron Spectra of Ionic Liquids: Experiments and Simulations of 1-Butyl-3-methylimidazolium Thiocyanate

We demonstrate a combined experimental and computational approach to probe the electronic structure and atomic environment of an ionic liquid, based on core level binding energies. The 1-butyl-3-methylimidazolium thiocyanate [C₄C₁Im][SCN] ionic liquid was studied using ab initio molecular dynamics, and results were compared against previously published and new experimental X-ray photoelectron spectroscopy (XPS) data. The long-held assumption that initial-state effects in XPS dominate the measured binding energies is proven correct, which validates the established premise that the ground state electronic structure of the ionic liquid can be inferred directly from XPS measurements. A regression model based upon site electrostatic potentials and intramolecular bond lengths is shown to account accurately for variations in core-level binding energies within the ionic liquid, demonstrating the important effect of long-range interactions on the core levels and throwing into question the validity of traditional single ion pair ionic liquid calculations for interpreting XPS data.

Comparison between Ab Initio Molecular Dynamics and OPLS-Based Force Fields for Ionic Liquid Solvent Organization

OPLS-based force fields (FFs) have been shown to provide accurate bulk-phase properties for a wide variety of imidazolium-based ionic liquids (ILs). However, the ability of OPLS to reproduce an IL solvent structure is not as well-validated given the relative lack of high-level theoretical or experimental data available for comparison. In this study, ab initio molecular dynamics (AIMD) simulations were performed for three widely used ILs: the 1-butyl-3-methylimidazolium cation with chloride, tetrafluoroborate, or hexafluorophosphate anions, that is, [BMIM][Cl], [BMIM][BF₄], and [BMIM][PF₆], respectively, as a basis for further assessment of two unique IL FFs: the ±0.8 charge-scaled OPLS-2009IL FF and the OPLS-VSIL FF. The OPLS-2009IL FF employs a traditional all-atom functional form, whereas the OPLS-VSIL FF was developed using a virtual site that offloads negative charge to inside the plane of the ring with careful attention given to reproducing hydrogen bonding. Detailed comparisons between AIMD and the OPLS FFs were made based on radial distribution functions (RDFs), combined distribution functions (CDFs), and spatial distribution functions (SDFs) to examine cation-anion interactions and π⁺-π⁺ stacking between the imidazolium rings. While both FFs were able to correctly capture the general solvent structure of these popular ILs, the OPLS-VSIL FF quantitatively reproduced interaction distances more accurately. In addition, this work provides further insights into the different short- and long-range structure patterns of these popular ILs.

New insights into chloromethyl-oxirane and chloromethyl-thiirane in liquid and solid phase from low-temperature infrared spectroscopy and ab initio modeling

A detailed study of the conformational landscape of chloromethyl-oxirane and chloromethyl-thiirane is here reported. The equilibrium of the three different conformers of the two molecules was assessed, using a joint approach of experimental and theoretical methods. High quality infrared spectroscopy measurements of the liquid and of the crystalline phases were interpreted with the aid of ab initio Molecular Dynamics (AIMD) simulations, anharmonic frequencies and free energy calculations, obtaining a very good reproduction of the experimental data. The modulation of the conformational equilibrium upon the addition of polar and non-polar solvents was computationally evaluated and results found a confirmation in experimental measures.

Choline Hydrogen Dicarboxylate Ionic Liquids by X-ray Scattering, Vibrational Spectroscopy and Molecular Dynamics: H-Fumarate and H-Maleate and Their Conformations

We explore the structure of two ionic liquids based on the choline cation and the monoanion of the maleic acid. We consider two isomers of the anion (H-maleate, the cis-isomer and H-fumarate, the trans-isomer) having different physical chemical properties. H-maleate assumes a closed structure and forms a strong intramolecular hydrogen bond whereas H-fumarate has an open structure. X-ray diffraction, infrared and Raman spectroscopy and molecular dynamics have been used to provide a reliable picture of the interactions which characterize the structure of the fluids. All calculations indicate that the choline cation prefers to connect mainly to the carboxylate group through OH⋯O interactions in both the compounds and orient the charged head N(CH3)3+ toward the negative portion of the anion. However, the different structure of the two anions affects the distribution of the ionic components in the fluid. The trans conformation of H-fumarate allows further interactions between anions through COOH and CO2- groups whereas intramolecular hydrogen bonding in H-maleate prevents this association. Our theoretical findings have been validated by comparing them with experimental X-ray data and infrared and Raman spectra.

Microstructural and Dynamical Heterogeneities in Ionic Liquids

Ionic liquids (ILs) are a special category of molten salts solely composed of ions with varied molecular symmetry and charge delocalization. The versatility in combining varied cation-anion moieties and in functionalizing ions with different atoms and molecular groups contributes to their peculiar interactions ranging from weak isotropic associations to strong, specific, and anisotropic forces. A delicate interplay among intra- and intermolecular interactions facilitates the formation of heterogeneous microstructures and liquid morphologies, which further contributes to their striking dynamical properties. Microstructural and dynamical heterogeneities of ILs lead to their multifaceted properties described by an inherent designer feature, which makes ILs important candidates for novel solvents, electrolytes, and functional materials in academia and industrial applications. Due to a massive number of combinations of ion pairs with ion species having distinct molecular structures and IL mixtures containing varied molecular solvents, a comprehensive understanding of their hierarchical structural and dynamical quantities is of great significance for a rational selection of ILs with appropriate properties and thereafter advancing their macroscopic functionalities in applications. In this review, we comprehensively trace recent advances in understanding delicate interplay of strong and weak interactions that underpin their complex phase behaviors with a particular emphasis on understanding heterogeneous microstructures and dynamics of ILs in bulk liquids, in mixtures with cosolvents, and in interfacial regions.

Coupled hydroxyl and ether functionalisation in EAN derivatives: the effect of hydrogen bond donor/acceptor groups on the structural heterogeneity studied with X-ray diffractions and fixed charge/polarizable simulations

We present a study by energy-dispersive X-ray diffraction of liquid 2-(2-hydroxyethoxy)ethan-1-ammonium nitrate, NH₃CH₂CH₂(OCH₂CH₂OH)⁺NO₃^- (22HHEAN). This ionic liquid is derived from the parent ethylammonium nitrate (EAN) with an ether link in the chain and a hydroxyl group in the terminal position. The absence of peaks at low-q values in the experimental diffraction curve indicates that the added polar groups and the high conformational isomerism of the cations alter strongly the nanosegregation of the parent EAN liquid. Aggregation between ionic species may involve hydrogen bonding between cations and anions and a variety of intermolecular hydrogen bonds between cations. Diffraction patterns are compared with the results of molecular dynamics simulations with two different force fields: the fixed point charge force field (GAFF) with different charge scaling protocols and the polarizable AMOEBA force field. Most point charge models lead to the appearance of a quite evident low q-peak which decreases gradually, when the percentage and type of the scaling (uniform vs. non-uniform) are increased. In the polarisable model and in the model where only anion charges are scaled to 20%, instead, the pre-peak is absent in agreement with our experiments.

Protic Ionic Liquids Based on the Alkyl-Imidazolium Cation: Effect of the Alkyl Chain Length on Structure and Dynamics

Protic ionic liquids are known to form extended hydrogen-bonded networks that can lead to properties different from those encountered in the aprotic analogous liquids, in particular with respect to the structure and transport behavior. In this context, the present paper focuses on a wide series of 1-alkyl-imidazolium bis(trifluoromethylsulfonyl)imide ionic liquids, [HC _nIm][TFSI], with the alkyl chain length ( n) on the imidazolium cation varying from ethyl ( n = 2) to dodecyl ( n = 12). A combination of several methods, such as vibrational spectroscopy, wide-angle X-ray scattering (WAXS), broadband dielectric spectroscopy, and ¹H NMR spectroscopy, is used to understand the correlation between local cation-anion coordination, nature of nanosegregation, and transport properties. The results indicate the propensity of the -NH site on the cation to form stronger H-bonds with the anion as the alkyl chain length increases. In addition, the position and width of the scattering peak q₁ (or the pre-peak), resolved by WAXS and due to the nanosegregation of the polar from the nonpolar domains, are clearly dependent on the alkyl chain length. However, we find no evidence from pulsed-field gradient NMR of a proton motion decoupled from molecular diffusion, hypothesized to be facilitated by the longer N-H bonds localized in the segregated ionic domains. Finally, for all protic ionic liquids investigated, the ionic conductivity displays a Vogel-Fulcher-Tammann dependence on inverse temperature, with an activation energy E_a that also depends on the alkyl chain length, although not strictly linearly.

Structure and dynamics of propylammonium nitrate-acetonitrile mixtures: An intricate multi-scale system probed with experimental and theoretical techniques

In this article, we report the study of structural and dynamical properties for a series of acetonitrile/propylammonium nitrate mixtures as a function of their composition. These systems display an unusual increase in intensity in their X-ray diffraction patterns in the low-q regime, and their ¹H-NMR diffusion-ordered NMR spectroscopy (DOSY) spectra display unusual diffusivities. However, the magnitude of both phenomena for mixtures of propylammonium nitrate is smaller than those observed for ethylammonium nitrate mixtures with the same cosolvent, suggesting that the cation alkyl tail plays an important role in these observations. The experimental X-ray scattering data are compared with the results of molecular dynamics simulations, including both ab initio studies used to interpret short-range interactions and classical simulations to describe longer range interactions. The higher level calculations highlight the presence of a strong hydrogen bond network within the ionic liquid, only slightly perturbed even at high acetonitrile concentration. These strong interactions lead to the symmetry breaking of the NO₃^- vibrations, with a splitting of about 88 cm^-1 in the ν₃ antisymmetric stretch. The classical force field simulations use a greater number of ion pairs, but are not capable of fully describing the longest range interactions, although they do successfully account for the observed concentration trend, and the analysis of the models confirms the nano-inhomogeneity of these kinds of samples.

Mapping the Free Energy of Lithium Solvation in the Protic Ionic Liquid Ethylammonuim Nitrate: A Metadynamics Study

Understanding lithium solvation and transport in ionic liquids is important due to their possible application in electrochemical devices. Using first-principles simulations aided by a metadynamics approach we study the free-energy landscape for lithium ions at infinite dilution in ethylammonium nitrate, a protic ionic liquid. We analyze the local structure of the liquid around the lithium cation and obtain a quantitative picture in agreement with experimental findings. Our simulations show that the lowest two free energy minima correspond to conformations with the lithium ion being solvated either by three or four nitrate ions with a transition barrier between them of 0.2 eV. Other less probable conformations having different solvation pattern are also investigated.