The Nature of Hydrogen Bonding in Protic Ionic Liquids

Small-Angle X-ray Scattering Study of the Amphiphilic Bulk Nanostructure of Tetraalkylammonium Deep Eutectic Solvents

Deep eutectic solvents (DESs) are low-melting mixtures, often prepared from a salt and a molecular hydrogen bond donor. Like ionic liquids, DESs that contain at least one sufficiently amphiphilic component can form bicontinuous nanostructures consisting of polar and nonpolar domains, although this has not been widely explored for many DES combinations. Here, the bulk nanostructures of DESs comprising tetraalkylammonium bromide salts (tetrabutylammonium bromide, tetraoctylammonium bromide, and methyltrioctylammonium bromide) with alkanols and alkanoic acids of systematically varied chain lengths (C2, C6, C8, and C10) as hydrogen bond donors have been studied. Small-angle X-ray scattering techniques were used to identify the relationship between the alkyl chain length and functionality of the hydrogen bond donor on the nature of the amphiphilic nanostructures formed. These findings demonstrated that the amphiphilic nanostructures of the DESs were not affected by the functional group on the hydrogen bond donor, with these nanostructures influenced primarily by both the absolute and relative alkyl chain lengths of the salt and hydrogen bond donor.

Exploiting Cation Structure and Water Content in Modulating the Acidity of Ammonium Hydrogen Sulfate Protic Ionic Liquids

In this paper, we investigated the effect of cation structure and water content on proton dissociation in alkylammonium [HSO4]- protic ionic liquids (ILs) doped with 20 wt % water and correlated this with experimental Hammett acidities. For pure systems, increased cation substitution resulted in a reduction in the number of direct anion-anion neighbors leading to larger numbers of small aggregates, which is further enhanced with addition of water. We also observed spontaneous proton dissociation from [HSO4]- to water only for primary amine-based protic ILs, preceded by the formation of an anion trimer motif. Investigation using DFT calculations revealed spontaneous proton dissociation from [HSO4]- to water can occur for each of the protic ILs investigated; however, this is dependent on the size of the anion aggregates. These findings are important in the fields of catalysis and lignocellulosic biomass, where solvent acidity is a crucial parameter in biomass fractionation and lignin chemistry.

H-bond network, interfacial tension and chain melting temperature govern phospholipid self-assembly in ionic liquids

HYPOTHESIS

The self-assembly structures and phase behaviour of phospholipids in protic ionic liquids (ILs) depend on intermolecular forces that can be controlled through changes in the size, polarity, and H-bond capacity of the solvent.

EXPERIMENTS

The structure and temperature stability of the self-assembled phases formed by four phospholipids in three ILs was determined by a combination of small- and wide-angle X-ray scattering (SAXS and WAXS) and small-angle neutron scattering (SANS). The phospholipids have identical phosphocholine head groups but different alkyl tail lengths and saturations (DOPC, POPC, DPPC and DSPC), while the ILs' amphiphilicity, H-bond network density and polarity are varied between propylammonium nitrate (PAN) to ethylammonium nitrate (EAN) to ethanolammonium nitrate (EtAN).

FINDINGS

The observed structures and phase behaviour of the lipids becomes more surfactant-like with decreasing average solvent polarity, H-bond network density and surface tension. In PAN, all the investigated phospholipids behave like surfactants in water. In EAN they exhibit anomalous phase sequences and unexpected transitions as a function of temperature, while EtAN supports structures that share characteristics with water and EAN. Structures formed are also sensitive to proximity to the lipid chain melting temperature.

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.

Intermolecular interactions in clusters of ethylammonium nitrate and 1-amino-1,2,3-triazole

The intermolecular interaction energies, including hydrogen bonds (H-bonds), of clusters of the ionic liquid ethylammonium nitrate (EAN) and 1-amino-1,2,3-triazole (1-AT) based deep eutectic propellants (DeEP) are examined. 1-AT is introduced as a neutral hydrogen bond donor (HBD) to EAN in order to form a eutectic mixture. The effective fragment potential (EFP) is used to examine the bonding interactions in the DeEP clusters. The resolution of the Identity (RI) approximated second order Møller-Plesset perturbation theory (RI-MP2) and coupled cluster theory (RI-CCSD(T)) are used to validate the EFP results. The EFP method predicts that there are significant polarization and charge transfer effects in the EAN:1-AT complexes, along with Coulombic, dispersion and exchange repulsion interactions. The EFP interaction energies are in good agreement with the RI-MP2 and RI-CCSD(T) results. The quasi-atomic orbital (QUAO) bonding and kinetic bond order (KBO) analyses are additionally used to develop a conceptual and semi-quantitative understanding of the H-bonding interactions as a function of the size of the system. The QUAO and KBO analyses suggest that the H-bonds in the examined clusters follow the characteristic hydrogen bonding three-center four electron interactions. The strongest H-bonding interactions between the (EAN)1:(1-AT)n and (EAN)2:(1-AT)n (n = 1-5) complexes are observed internally within EAN; that is, between the ethylammonium cation [EA]+ and the nitrate anion ([NO3]-). The weakest H-bonding interactions occur between [NO3]- and 1-AT. Consequently, the average strengths of the H-bonds within a given (EAN)x:(1-AT)n complex decrease as more 1-AT molecules are introduced into the EAN monomer and EAN dimer. The QUAO bonding analysis suggests that 1-AT in (EAN)x:(1-AT)n can act as both a HBD and a hydrogen bond acceptor simultaneously. It is observed that two 1-AT molecules can form H-bonds to each other. Although the KBOs that correspond to H-bonding interactions in [EA]+:1-AT, [NO3]-:1-AT and between two 1-AT molecules are weaker than the H-bonds in EAN, those weak H-bond networks with 1-AT could be important to form a stable DeEP.

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.

Electron Density Learning of Z-Bonds in Ionic Liquids and Its Application

Ionic liquids (ILs) exhibit fascinating properties due to special Z-bonds and have been widely used in electrochemical systems. The local Z-bond networks potentially cause a discrepancy in electrochemical properties. Understanding the correlations between the Z-bond energy (EZ-bond) and the electrochemical properties is helpful to identify appropriate ILs. It is difficult to estimate the correlations from single density functional theory calculations or molecular dynamic simulations. In this work, a machine learning model targeting the electronic density (ρBCP) of Z-bonds has been trained successfully, as expected for use in systems above the nanoscale size. The connection between the EZ-bond and the electrochemical potential window in ILs@TiO2, as well as that between the EZ-bond and the charge carrier mobility in ILs-PEDOT:Tos@SiO2, was separately investigated. This study highlights an efficient model for predicting ρBCP in nanoscale systems and anticipates exploring the connection between Z-bonds and the electrochemical properties of IL-based systems.

Unusual phosphatidylcholine lipid phase behavior in the ionic liquid ethylammonium nitrate

HYPOTHESIS

The forces that govern lipid self-assembly ionic liquids are similar to water, but their different balance can result in unexpected behaviour.

EXPERIMENTS

The self-assembly behaviour and phase equilibria of two phospholipids, 1,2-distearoyl-sn-glycero-3-phosphocholine (DSPC) and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), in the most common protic ionic liquid, ethylammonium nitrate (EAN) have been investigated as function of composition and temperature by small- and wide-angle X-ray scattering (SAXS/WAXS) and small-angle neutron scattering (SANS).

FINDINGS

Both lipids form unusual self-assembly structures and show complex and unexpected phase behaviour unlike that seen in water; DSPC undergoes a gel L_β to crystalline L_c phase transition on warming, while POPC forms worm-like micelles L₁ upon dilution. This surprising phase behaviour is attributed to the large size of the EAN ions that solvate the lipid headgroup compared to water changing amphiphile packing. Weaker H-bonding between EAN and lipid headgroups also contributes. These results provide new insight for the design of lipid based nanostructured materials in ionic liquids with atypical properties.

Proton transfer between sulfonic acids and various propylamines by density functional theory calculations

CONTEXT

Proton transfer in acid-base systems is not well understood. Some acid-base reactions do not proceed to the extent that is expected from the difference in the pK_a values between the base and acid in aqueous solutions, yet some do. In that regard, we have computationally studied the process of proton transfer from the acids of varying strength (benzenesulfonic acid (BSu), methansulfonic acid (MsO), and sulfuric acid (SA)) to the amines with different numbers of propyl substituents on the nitrogen atom (propylamine (PrA), dipropylamine (DPrA), and tripropylamine (TPrA)) upon complexation. Density functional theory calculations were used to thoroughly examine the energetic and structural aspects of the molecular complexes and/or ionic pairs resulting from the acid-base interaction. The potential energy curves along the proton transfer coordinate in these acid-amine systems were analyzed. The change in free energies accompanying the molecular complexes and ionic pair formations was calculated, and the relationship between the energy values and the ΔРА parameter (difference in proton affinity of the acid anion and amine) was established. The larger ΔРА values were found to be unfavorable for the formation of ionic pairs. Using structural, energy, QTAIM, and NBO analyses, we determined that the hydrogen bonds in the molecular complexes PrA-MsO and PrA-BSu are stronger than those in their corresponding ionic pairs. The ionic pairs with the TPrA cation possess the strongest hydrogen bonds of all the ionic pairs being studied, regardless of the anion. The results showed that hydrogen bonding interactions in the molecular complexes contribute significantly to the energies of the acid-base interaction, while in the ionic pairs, the most important energy contribution comes from Coulomb interactions, followed by hydrogen bonding and dispersion forces. The ionic pairs with propylammonium, dipropylammonium, and tripropylammonium cations have stronger ion-ion interactions than tetrapropylammonium (TetPrA)-containing ionic pairs with the same anions. This effect rises with the order of the cation: TetPrA → TPrA → DPrA → PrA, and the sequence of anions is SA → BSu → MsO. The results obtained here expand the concept of acid-base interaction and provide an alternative to experimental searches for suitable acids and bases to obtain new types of protic ionic liquids.

METHODS

All quantum-chemical calculations were carried out by using the DFT/B3LYP-GD3/6-31++G(d,p) level as implemented in the Gaussian 09 software package. For the resulting structures, the electron density distribution was analyzed by the "atoms in molecules" (QTAIM) and the natural bond orbital (NBO) methods on the wave functions obtained at the same level of theory by AIMAll Version 10.05.04 and Gaussian NBO Version 3.1 programs, respectively.

Quantification and Distribution of Three Types of Hydrogen Bonds in Mixtures of an Ionic Liquid with the Hydrogen-Bond-Accepting Molecular Solvent DMSO Explored by Neutron Diffraction and Molecular Dynamics Simulations

The concept of hydrogen bonding is celebrating its 100th birthday. Hydrogen bonds (H-bonds) play a key role in the structure and function of biological molecules, the strength of materials, and molecular binding. Herein, we study H-bonding in mixtures of a hydroxyl-functionalized ionic liquid with the neutral, H-bond-accepting molecular liquid dimethylsulfoxide (DMSO) using neutron diffraction experiments and molecular dynamics simulations. We report the geometry, strength, and distribution of three different types of H-bond OH···O, formed between the hydroxyl group of the cation and either the oxygen atom of another cation, the counteranion, or the neutral molecule. Such a variety of different strengths and distributions of H-bonds in one single mixture could hold the promise of providing solvents with potential applications in H-bond-related chemistry, for example, to alter the natural selectivity patterns of catalytic reactions or the conformation of catalysts.

New Protic Ionic Liquids as Potential Additives to Lubricate Si-Based MEMS/NEMS

The motivation for this work was to develop new protic ionic liquids (PILs) as additives for the lubrication of micro and nanoelectromechanical systems (MEMS and NEMS). Ten PILs based on the combination of methylimidazolium ([MIMH]), 4-picolinium ([4-picH]), pyridinium ([PyrH]), 1,8-diazabicyclo[5.4.0]-undec-7-ene-8-ium ([DBUH]) and tetramethylguanidinium ([TMGH]) cations with hydrogen sulfate([HSO4]) and mesylate ([MeSO3]) anions were tested as additives in polyethylene glycol (PEG200) to lubricate steel/silicon and silicon/silicon contacts. The best additive was [4-picH][HSO4], which adsorbed strongly on the Si surface, leading to a protective film that reduced wear by up to 15 times compared to PEG200.

The structure of protic ionic liquids based on sulfuric acid, doped with excess of sulfuric acid or with water

Neutron scattering with isotopic substitution was used to study the structure of concentrated sulfuric acid, and two protic ionic liquids (PILs): a Brønsted-acidic PIL, synthesised using pyridine and excess of sulfuric acid, [Hpy][HSO₄]·H₂SO₄, and a hydrated PIL, in which an equimolar mixture of sulfuric acid and pyridine has been doped with water, [Hpy][HSO₄]·2H₂O. Brønsted acidic PILs are excellent solvents/catalysts for esterifications, driving reaction to completion by phase-separating water and ester products. Water-doped PILs are efficient solvents/antisolvents in biomass fractionation. This study was carried out to provide an insight into the relationship between the performance of PILs in the two respective processes and their liquid structure. It was found that a persistent sulfate/sulfuric acid/water network structure was retained through the transition from sulfuric acid to PILs, even in the presence of 2 moles (∼17 wt%) of water. Hydrogen sulfate PILs have the propensity to incorporate water into hydrogen-bonded anionic chains, with strong and directional hydrogen bonds, which essentially form a new water-in-salt solvent system, with its own distinct structure and physico-chemical properties. It is the properties of this hydrated PIL that can be credited both for the good performance in esterification and beneficial solvent/antisolvent behaviour in biomass fractionation.

Multiple Ensembles of the Hydrogen-bonded Network in Ethylammonium Nitrate versus Water from Vibrational Spectral Dynamics of SCN- Probe

We performed classical molecular dynamics simulations to monitor the structural interactions and ultrafast dynamical and spectral response in the protic ionic liquid, ethylammonium nitrate (EAN) and water using the nitrile stretching mode of thiocyanate ion (SCN-) as the vibrational probe. The normalized stretch frequency distribution of nitrile stretch of SCN- attains an asymmetric shape in EAN, indicating the existence of more than one hydrogen-bonding environment in EAN. We computed the 2D IR spectrum from classical trajectories, applying the response function formalism. Spectral diffusion dynamics in EAN undergo an initial rattling of the SCN - inside the local ion-cage occurring at a timescale of 0.10 ps, followed by the breakup of the ion-cage activating molecular diffusion at 7.86 ps timescale. In contrast, the dynamics of structural reorganization occur at a timescale of 0.58 ps in H 2 O. Hence, the time dependence of the frequency-frequency correlation function decay hints at the local molecular structure and ultrafast ion dynamics of the SCN - probe. The loss of frequency correlation read from the peak shape changes in the 2D correlation spectrum as a function of waiting time is faster in H 2 O than in EAN due to the enhanced structural ordering and higher viscosity of the latter. We provide an atomic-level interpretation of the solvation environment around SCN - in EAN and water, which indicates the multiple ensembles of the hydrogen bond network in EAN.

A DLPNO-CCSD(T) benchmarking study of intermolecular interactions of ionic liquids

The accuracy of correlation energy recovered by coupled cluster single-, double-, and perturbative triple-excitations, CCSD(T), has led to the method being considered the gold standard of computational chemistry. The application of CCSD(T) has been limited to medium-sized molecular systems due to its steep scaling with molecular size. The recent development of alternative domain-based local pair natural orbital coupled-cluster method, DLPNO-CCSD(T), has significantly broadened the range of chemical systems to which CCSD(T) level calculations can be applied. Condensed systems such as ionic liquids (ILs) have a large contribution from London dispersion forces of up to 150 kJ mol^-1 in large-scale clusters. Ionic liquids show appreciable charge transfer effects that result in the increased valence orbital delocalization over the entire ionic network, raising the question whether the application of methods based on localized orbitals is reliable for these semi-Coulombic materials. Here the performance of DLPNO-CCSD(T) is validated for the prediction of correlation interaction energies of two data sets incorporating single-ion pairs of protic and aprotic ILs. DLPNO-CCSD(T) produced results within chemical accuracy with tight parameter settings and a non-iterative treatment of triple excitations. To achieve spectroscopic accuracy of 1 kJ mol^-1 , especially for hydrogen-bonded ILs and those containing halides, the DLPNO settings had to be increased by two orders of magnitude and include the iterative treatment of triple excitations, resulting in a 2.5-fold increase in computational cost. Two new sets of parameters are put forward to produce the performance of DLPNO-CCSD(T) within chemical and spectroscopic accuracy.

Effect of Water and Salt on the Colloidal Stability of Latex Particles in Ionic Liquid Solutions

The colloidal stability of sulfate (SL) and polyimidazolium-modified sulfate (SL-IP-2) latex particles was studied in an ionic liquid (IL) of ethylammonium nitrate (EAN) and its water mixtures. Aggregation rates were found to vary systematically as a function of the IL-to-water ratio. Repulsive electrostatic interactions between particles dominated at low IL concentrations, while they were significantly screened at intermediate IL concentrations, leading to destabilization of the dispersions. When the IL concentration was further increased, the aggregation of latex particles slowed down due to the increased viscosity and finally, a striking stabilization was observed in the IL-rich regime close to the pure IL solvent. The latter stabilization is due to the formation of IL layers at the interface between particles and IL, which induce repulsive oscillatory forces. The presence of the added salt in the system affected differently the structure of the interfaces around SL and SL-IP-2 particles. The sign of the charge and the composition of the particle surfaces were found to be the most important parameters affecting the colloidal stability. The nature of the counterions also plays an important role in the interfacial properties due to their influence on the structure of the IL surface layers. No evidence was observed for the presence of long-range electrostatic interactions between the particles in pure ILs. The results indicate that the presence of even low concentrations of water and salt in the system (as undesirable impurities) can strongly alter the interfacial structure and thus, the aggregation mechanism in particle IL dispersions.

Alkylimidazolium Protic Ionic Liquids: Structural Features and Physicochemical Properties

We focus on a series of protic ionic liquids (PILs) with imidazolium and alkylimidazolium (1R3HIm, R = methyl, ethyl, propyl, and butyl) cations. Using the literature data and our experimental results on the thermal and transport properties, we analyze the effects of the anion nature and the alkyl radical length in the cation structure on the above properties. DFT calculations in gas and solvent phase have resulted in microscopic insights into the structure and cation-anion binding in these PILs. We show that the higher thermodynamic stability of an ion pair raises the PIL decomposition temperature. The melting points of the salts with the same cation decrease as the hydrocarbon radical in the cation becomes longer, which correlates with the weaker ion-ion interaction in the ion pairs. A comparative analysis of the protic ILs and corresponding ILs (1R3MeIm) with the same radical (R) in the cation structure and the same anion has been performed. The lower melting points of the ILs with 1R3MeIm cations are assumed to result from the weakening both of the ion-ion interaction and hydrogen bond.

Spectroscopic Signatures of Hydrogen-Bonding Motifs in Protonic Ionic Liquid Systems: Insights from Diethylammonium Nitrate in the Solid State

Diethylammonium nitrate, [N0 0 2 2][NO3], and its perdeuterated analogue, [N D D 2 2] [NO3], were structurally characterized and studied by infrared, Raman, and inelastic neutron scattering (INS) spectroscopy. Using these experimental data along with state-of-the-art computational materials modeling, we report unambiguous spectroscopic signatures of hydrogen-bonding interactions between the two counterions. An exhaustive assignment of the spectral features observed with each technique has been provided, and a number of distinct modes related to NH···O dynamics have been identified. We put a particular emphasis on a detailed interpretation of the high-resolution, broadband INS experiments. In particular, the INS data highlight the importance of conformational degrees of freedom within the alkyl chains, a ubiquitous feature of ionic liquid (IL) systems. These findings also enable an in-depth physicochemical understanding of protonic IL systems, a first and necessary step to the tailoring of hydrogen-bonding networks in this important class of materials.

An ultrafast vibrational study of dynamical heterogeneity in the protic ionic liquid ethyl-ammonium nitrate. I. Room temperature dynamics

Using ultrafast two-dimensional infrared spectroscopy (2D-IR), a vibrational probe (thiocyanate, SCN^-) was used to investigate the hydrogen bonding network of the protic ionic liquid ethyl-ammonium nitrate (EAN) in comparison to H₂O. The 2D-IR experiments were performed in both parallel (⟨ZZZZ⟩) and perpendicular (⟨ZZXX⟩) polarizations at room temperature. In EAN, the non-Gaussian lineshape in the FTIR spectrum of SCN^- suggests two sub-ensembles. Vibrational relaxation rates extracted from the 2D-IR spectra provide evidence of the dynamical differences between the two sub-ensembles. We support the interpretation of two sub-ensembles with response function simulations of two overlapping bands with different vibrational relaxation rates and, otherwise, similar dynamics. The measured rates for spectral diffusion depend on polarization, indicating reorientation-induced spectral diffusion (RISD). A model of restricted molecular rotation (wobbling in a cone) fully describes the observed spectral diffusion in EAN. In H₂O, both RISD and structural spectral diffusion contribute with similar timescales. This complete characterization of the dynamics at room temperature provides the basis for the temperature-dependent measurements in Paper II of this series.

Protic Ionic Liquid as Reagent, Catalyst, and Solvent: 1-Methylimidazolium Thiocyanate

We propose a new concept of the triple role of protic ionic liquids with nucleophilic anions: a) a regenerable solvent, b) a Brønsted acid inducing diverse transformations via general acid catalysis, and c) a source of a nucleophile. The efficiency of this strategy was demonstrated using thiocyanate-based protic ionic liquids for the ring-opening of donor-acceptor cyclopropanes. A wide variety of activated cyclopropanes were found to react with 1-methylimidazolium thiocyanate under mild metal-free conditions via unusual nitrogen attack of the ambident thiocyanate ion on the electrophilic center of the three-membered ring affording pyrrolidine-2-thiones bearing donor and acceptor substituents at the C(5) and C(3) atoms, respectively, in a single time-efficient step. The ability of 1-methylimidazolium thiocyanate to serve as a triplex reagent was exemplarily illustrated by (4+2)-annulation with 1-acyl-2-(2-hydroxyphenyl)cyclopropane, epoxide ring-opening and other organic transformations.

Extension of the CL&Pol Polarizable Force Field to Electrolytes, Protic Ionic Liquids, and Deep Eutectic Solvents

The polarizable CL&Pol force field presented in our previous study, Transferable, Polarizable Force Field for Ionic Liquids (J. Chem. Theory Comput. 2019, 15, 5858, DOI: http://doi.org/10.1021/acs.jctc.9b0068910.1021/acs.jctc.9b00689), is extended to electrolytes, protic ionic liquids (PIL), deep eutectic solvents (DES), and glycols. These systems are problematic in polarizable simulations because they contain either small, highly charged ions or strong hydrogen bonds, which cause trajectory instabilities due to the pull exerted on the induced dipoles. We use a Tang-Toennies (TT) function to dampen, or smear, the interactions between charges and induced dipole at a short range involving small, highly charged atoms (such as hydrogen or lithium), thus preventing the "polarization catastrophe". The new force field gives stable trajectories and is validated through comparison with experimental data on density, viscosity, and ion diffusion coefficients of liquid systems of the above-mentioned classes. The results also shed light on the hydrogen-bonding pattern in ethylammonium nitrate, a PIL, for which the literature contains conflicting views. We describe the implementation of the TT damping function, of the temperature-grouped Nosé-Hoover thermostat for polarizable molecular dynamics (MD) and of the periodic perturbation method for viscosity evaluation from non-equilibrium trajectories in the LAMMPS MD code. The main result of this work is the wider applicability of the CL&Pol polarizable force field to new, important classes of fluids, achieving robust trajectories and a good description of equilibrium and transport properties in challenging systems. The fragment-based approach of CL&Pol will allow ready extension to a wide variety of PILs, DES, and electrolytes.

Effect of Brønsted Acidity on Ion Conduction in Fluorinated Acetic Acid and N-Methylimidazole Equimolar Mixtures as Pseudo-protic Ionic Liquids

To clarify proton conduction mechanism in protic ionic liquids (PILs) and pseudo-PILs (pPILs), equimolar mixtures of N-methylimidazole (C₁Im) with fluorinated acetic acids were investigated by Raman spectroscopy, X-ray scattering, and dielectric relaxation spectroscopy (DRS). Only the ionic species exist in the equimolar mixture of C₁Im and HTFA (HTFA: trifluoroacetic acid). On the other hand, the equimolar mixture of C₁Im and HDFA (HDFA: difluoroacetic acid) consists of both ionic and electrically neutral species. In particular, not only the electrostatic but also van der Waals interactions with the F atoms were observed in the liquid structures of both [C₁hIm⁺][TFA^-] and [C₁hIm⁺][DFA^-]. The concept for proton conduction mechanism that we have proposed in previous study was revisited; the proton conduction mechanism could be classified with two linear free energy relationship lines for proton exchange reaction and translation/rotation of proton carriers. Our results exhibit that the proton conduction mechanism changes from proton hopping to vehicle mechanism with increasing acidity of an acid HA in PILs.

Liquid Structure of Single and Mixed Cation Alkylammonium Bromide Urea Deep Eutectic Solvents

The liquid structures of three alkyl ammonium bromide and urea DESs, ethylammonium bromide:urea (1:1), butylammonium bromide:urea (1:1), and ethylammonium bromide/butylammonium bromide:urea (0.5:0.5:1), have been studied using small-angle neutron diffraction with H/D substituted sample contrasts. The diffraction data was fit using empirical potential structure refinement (EPSR). An amphiphilic nanostructure was found in all DESs due to cation alkyl chains being solvophobically excluded from charged domains, and due to clustering together. The polar domain was continuous in all three DESs, whereas the apolar domain was continuous for the butylammonium DES and in the mixed DES, but not the ethylammonium DES. This is attributed to solvophobic interactions being weaker for the short ethyl chain. Surprisingly, the urea also forms large clusters in all three DESs. In ethylammonium bromide:urea (1:1), urea-urea orientations are mainly perpendicular, but in butylammonium bromide:urea (1:1) and the mixed system in-plane and perpendicular arrangements are found. The liquid nanostructures found in this work, especially for the ethylammonium DES, are different from those found previously for the corresponding DESs formed using glycerol, revealing that the DES amphiphilic nanostructure is sensitive to the nature of the HBD (hydrogen bond donor).

Catanionic Surfactant Self-Assembly in Protic Ionic Liquids

Mixing of cationic and anionic surfactants in water can result in pseudo-double-tailed catanionic surfactant ion pairs that form lamellar phases or vesicles that are unstable toward electrolyte addition. Here we show that despite the very high ionic strengths, catanionic surfactants counterintuitively form a wider variety of self-assembled aggregates in pure ionic liquids (ILs, pure salts in a liquid phase) than in water, including micelles, vesicles, and lyotropic phases. Self-assembled structures only form when the IL is sufficiently polar to drive self-assembly through electrostatic interactions and/or H-bond networks, but the catanionic effect is manifested only when the IL does not itself exhibit pronounced amphiphilic nanostructure. This enables the type of catanionic aggregate formed to be designed by changing the hydrogen bonds between the ions through variation of the structures of the cation and anion. These results reveal an entirely new way of controlling catanionic surfactant self-assembly under nonaqueous and high-salt conditions.

Hydrogen bonding and charge transport in a protic polymerized ionic liquid

Hydrogen bonding and charge transport in the protic polymerized ionic liquid poly[tris(2-(2-methoxyethoxy)ethyl)ammoniumacryloxypropyl sulfonate] (PAAPS) are studied by combining Fourier transform infrared (FTIR) and broadband dielectric spectroscopy (BDS) in a wide temperature range from 170 to 300 K. While the former enables to determine precisely the formation of hydrogen bonds and other moiety-specific quantized vibrational states, the latter allows for recording the complex conductivity in a spectral range from 10-2 to 10+9  Hz. A pronounced thermal hysteresis is observed for the H-bond network formation in distinct contrast to the reversibility of the effective conductivity measured by BDS. On the basis of this finding and the fact that the conductivity changes with temperature by orders of magnitude, whereas the integrated absorbance of the N-H stretching vibration (being proportional to the number density of protons in the hydrogen bond network) changes only by a factor of 4, it is concluded that charge transport takes place predominantly due to hopping conduction assisted by glassy dynamics (dynamic glass transition assisted hopping) and is not significantly affected by the establishment of H-bonds.

Ion Pair Structures and Hydrogen Bonding in RnNH4-n Alkylammonium Ionic Liquids with Hydrogen Sulfate and Mesylate Anions by DFT Computations

Density function theory calculations are employed to study the interaction of amines bearing different numbers of alkyl substituents of different sizes on the nitrogen atom with sulfuric and methanesulfonic acids. The proton affinities of the studied amines are calculated, and it is shown that the higher the value is, the more probable is its protonation. The most stable structures of the ion pairs resulting from the acid-base interaction are obtained and characterized. The geometric parameters of the ion pairs and the characteristics derived from the NBO and QTAIM analysis show that there are hydrogen bonding interactions between the cation and the anion. The hydrogen bonding character of the ion pairs and the strength of the interaction between the ions strongly depend on the nature of the cation itself. The interaction between the ions in the ion pairs weakens with the increase in the cation size. The trend of change in the structural parameters of the H-bonds and energetic characteristics in the cation series for the studied ion pairs is not dependent on the nature of the anion.

Ion hydration and association in aqueous potassium tetrahydroxyborate solutions

Potassium tetrahydroxyborate solution is a significant material in the borate solution family, but there is limited knowledge about hydration structures and interactions of K⁺, [B(OH)₄^-], and water. In this study, the X-ray diffraction measurements of potassium tetrahydroxyborate solutions have been made. The experimental structure factors are subjected to Empirical Potential Structure Refinement (EPSR) modeling to reveal the details of ion hydration and association in the aqueous solutions. This study shows that the O(W)-O(W) distance of water in the studied solutions ranges from 2.82 to 2.76 Å with a coordination number that ranges from 4.7 ± 1.4 to 3.1 ± 1.3 when the value of the water-salt molar ratio (WSR) is decreased from 30 to 6. The addition of ions slightly affects the tetrahedral structure of water even when the concentration of ions is high. The first hydration distance of K⁺ remained at ∼2.67 Å, whereas the value of the coordination number (CN) decreased from 5.4 ± 1.3 to 3.9 ± 1.5 when the concentration of the borate solution was increased. The hydration ability of K⁺ was weak and almost did not have a fixed local hydration structure. The pair distribution function (PDF) of g_B-O(W)(r) shows that [B(OH)₄^-] has a broad hydration distance from 2.9 to 5.4 Å because of the complex interactive relationship between K⁺, [B(OH)₄^-] and water. There is a competitive hydration between K⁺ and [B(OH)₄^-]. Both the X-ray diffraction and DFT-based calculations confirm that the main species is monodentate contact ion pairs when WSR = 30, bidentate contact ion pairs when WSR = 14, and triple contact ion pairs when WSR = 6. These results will provide a new understanding about potassium tetrahydroxyborate solution.

Unexpected Strong Acidity Enhancing the Effect in Protic Ionic Liquids Quantified by Equilibrium Acidity Studies: A Crucial Role of Cation Structures on Dictating the Solvation Properties

The equilibrium acidities for several series of structural and electronic different organic acids were measured in 3-pyrrolidinium-based aprotic and protic ionic liquid (IL) analogues, that is, [Bmpy][NTf₂], [BpyH][NTf₂], and [PyH₂][NTf₂], by the UV-Vis method. The acidities of neutral acids are found to be much stronger in the protic ILs (PILs) than aprotic ILs (AILs), and the acidifying effect in the two PILs roughly increases proportionally to the number of protons in the cation of the PIL. On the other hand, interestingly, the cationic N⁺-H acids exhibit similar acidities in [Bmpy][NTf₂] and [BpyH][NTf₂] but much weaker than those in [PyH₂][NTf₂]. The Hammett ρ values for the acidic dissociation of para-substituted benzoic acids in two PILs are about the same as that in water (1.00) but significantly smaller than that in the AIL [Bmpy][NTf₂] (2.66). The correlations between the acidities in the PILs and water show double-linear relationships with different slopes and intercepts for the neutral and cationic acids. These, together with previous observations in the PILs [DBUH][OTf] (DBU = 1,5-diazabicyclo(5.4.0)-5-undecene) and EAN (ethylammonium nitrate), clearly indicate that the structure of the cation plays a subtle but a crucial role on sensing the electronic nature of solutes and strongly affects the solvation behaviors of PILs.

Amphiphilic nanostructure in choline carboxylate and amino acid ionic liquids and solutions

The liquid structures of six choline carboxylate/amino acid ionic liquids (bio-ILs) and their mixtures with water and various n-alkanols have been investigated by small-angle X-ray scattering (SAXS). The ILs exhibit long-range amphiphilic nanostructure comprised of polar and apolar domains that can be controlled by choice of anion, and which is tolerant to water dilution. Mixtures with n-alkanols can lead to marked changes in domain size and ordering. Utilising the Teubner-Strey model, we find amphiphilicity factors in many of these mixtures are comparable to those observed in conventional microemulsions, and that cooperative assembly in bio-IL/alkanol mixtures can enhance amphiphilicity, with potential to improve performance in a range of applications.