Monte Carlo QM/MM simulation studies of the Cannizzaro reaction in ionic liquids for improved biofuel production

The conversion of biomass to 5-hydroxymethylfurfural (HMF) holds substantial promise as a renewable energy source. Notably, HMF can be transformed into 2,5-bis(hydroxymethyl)furan (BHMF), a crucial reactant in biofuel production, but requires harsh operating conditions, H2, and precious metal catalysts. A recently reported Cannizzaro reaction of HMF to BHMF, characterized by its efficiency, mild conditions, and eco-friendliness, instead employed ionic liquids (ILs) to achieve high yields. In this study, combined quantum mechanical and molecular mechanical (QM/MM) simulations in conjunction with Metropolis Monte Carlo statistical mechanics and free-energy perturbation theory utilized M06-2X/6-31+G(d), PDDG/PM3, and the OPLS-VSIL force field to uncover important solute-solvent interactions present in the HMF to BHMF reaction pathway. The Cannizzaro reaction was examined in water and in five ILs composed of the 1-butyl-3-methylimidazolium [BMIM] cation coupled to hexafluorophosphate, tetrafluoroborate, thiocyanate, chloride, and bromide. Energetic and structural analysis of the rate-determining hydride transfer between HMF and the hydride-donor anion HMFOH- attributed the enhanced reactivity to highly organized solvent interactions featuring (1) hydrogen bonding between the ring protons of [BMIM] and the negatively charged carbonyl oxygen atoms on the transition structure, (2) favorable electrostatic interactions between the IL anions and solute hydroxyl groups, and (3) beneficial π-π stacking interactions between [BMIM] and the two aromatic rings present in HMF and HMFOH-.

Wave mechanics in an ionic liquid mixture

Experimental measurements of interactions in ionic liquids and concentrated electrolytes over the past decade or so have revealed simultaneous monotonic and oscillatory decay modes. These observations have been hard to interpret using classical theories, which typically allow for just one electrostatic decay mode in electrolytes. Meanwhile, substantial progress in the theoretical description of dielectric response and ion correlations in electrolytes has illuminated the deep connection between density and charge correlations and the multiplicity of decay modes characterising a liquid electrolyte. The challenge in front of us is to build connections between the theoretical expressions for a pair of correlation functions and the directly measured free energy of interaction between macroscopic surfaces in experiments. Towards this aim, we here present measurements and analysis of the interactions between macroscopic bodies across a fluid mixture of two ionic liquids of widely diverging ionic size. The measured oscillatory interaction forces in the liquid mixtures are significantly more complex than for either of the pure ionic liquids, but can be fitted to a superposition of two oscillatory and one monotonic mode with parameters matching those of the pure liquids. We discuss this empirical finding, which hints at a kind of wave mechanics for interactions in liquid matter.

Unveiling the Rational Development of Stimuli-Responsive Silk Fibroin-Based Ionogel Formulations

We present an approach for the rational development of stimuli-responsive ionogels which can be formulated for precise control of multiple unique ionogel features and fill niche pharmaceutical applications. Ionogels are captivating materials, exhibiting self-healing characteristics, tunable mechanical and structural properties, high thermal stability, and electroconductivity. However, the majority of ionogels developed require complex chemistry, exhibit high viscosity, poor biocompatibility, and low biodegradability. In our work, we overcome these limitations. We employ a facile production process and strategically integrate silk fibroin, the biocompatible ionic liquids (ILs) choline acetate ([Cho][OAc]), choline dihydrogen phosphate ([Cho][DHP]), and choline chloride ([Cho][Cl]), traditional pharmaceutical excipients, and the model antiepileptic drug phenobarbital. In the absence of ILs, we failed to observe gel formation; yet in the presence of ILs, thermoresponsive ionogels formed. Systems were assessed via visual tests, transmission electron microscopy, confocal reflection microscopy, dynamic light scattering, zeta potential and rheology measurements. We formed diverse ionogels of strengths ranging between 18 and 642 Pa. Under 25 °C storage, formulations containing polyvinylpyrrolidone (PVP) showed an ionogel formation period ranging over 14 days, increasing in the order of [Cho][DHP], [Cho][OAc], and [Cho][Cl]. Formulations lacking PVP showed an ionogel formation period ranging over 32 days, increasing in the order of [Cho][OAc], [Cho][DHP] and [Cho][Cl]. By heating from 25 to 60 °C, immediately following preparation, thermoresponsive ionogels formed below 41 °C in the absence of PVP. Based on our experimental results and density functional theory calculations, we attribute ionogel formation to macromolecular crowding and confinement effects, further enhanced upon PVP inclusion. Holistically, applying our rational development strategy enables the production of ionogels of tunable physicochemical and rheological properties, enhanced drug solubility, and structural and energetic stability. We believe our rational development approach will advance the design of biomaterials and smart platforms for diverse drug delivery applications.

Uncovering the Critical Factors that Enable Extractive Desulfurization of Fuels in Ionic Liquids and Deep Eutectic Solvents from Simulations

Environmental regulatory agencies have implemented stringent restrictions on the permissible levels of sulfur compounds in fuel to reduce harmful emissions and improve air quality. Problematically, traditional desulfurization methods have shown low effectiveness in the removal of refractory sulfur compounds, e.g., thiophene (TS), dibenzothiophene (DBT), and 4-methyldibenzothiophene (MDBT). In this work, molecular dynamics (MD) simulations and free energy perturbation (FEP) have been applied to investigate the use of ionic liquids (ILs) and deep eutectic solvents (DESs) as efficient TS/DBT/MDBT extractants. For the IL simulations, the selected cation was 1-butyl-3-methylimidazolium [BMIM] and the anions included chloride [Cl], thiocyanate [SCN], tetrafluoroborate [BF₄], hexafluorophosphate [PF₆], and bis(trifluoromethylsulfonyl)amide [NTf₂]. The DESs were composed of choline chloride with ethylene glycol (CCEtg) or with glycerol (CCGly). Calculation of excess chemical potentials predicted the ILs to be more promising extractants with energies lower by 1-3 kcal/mol compared to DESs. Increasing IL anion size was positively correlated to enhanced solvation of S-compounds, which was influenced by energetically dominant solute-anion interactions and favorable solute-[BMIM] π-π stacking. For the DESs, the solvent components offered a range of synergistic, yet comparatively weaker, electrostatic interactions that included hydrogen bonding and cation-π interactions. An in-depth analysis of the structure of IL and DES systems is presented, along with a discussion of the critical factors behind experimental trends of S-compound extraction efficiency.

Effects of Methylating Imidazolium-Based Ionic Liquids on Viscosity: New Insights from the Compensated Arrhenius Formalism

Methylation of the C(2) carbon on imidazolium-based room temperature ionic liquids (RTILs) causes an unexpected increase in viscosity when paired with the anion bis(trifluoromethylsulfonamide) [Tf₂N]^-, but the viscosity decreases when the methylated imidazolium is paired with a tetracyanoborate [B(CN)₄]^- anion. This paper investigates these different observations in viscosity using the compensated Arrhenius formalism (CAF) for fluidity (inverse viscosity), which assumes fluidity to be a thermally activated process. CAF activation energies are determined for imidazolium [Tf₂N]^- and methylated imidazolium [Tf₂N]^- and compared to imidazolium [B(CN)₄]^- and methylated imidazolium [B(CN)₄]^-. The results show that the activation energy increases with methylation for [Tf₂N]^-, but it decreases with methylation for [B(CN)₄]^-. The CAF results also yield information concerning the entropy of activation, which are compared for the two systems.

Insights into structure-property relationships in ionic liquids using cyclic perfluoroalkylsulfonylimides

Room temperature ionic liquids of cyclic sulfonimide anions ncPFSI (ring size: n = 4-6) with the cations [EMIm]+ (1-ethyl-3-methylimidazolium), [BMIm]+ (1-butyl-3-methylimidazolium) and [BMPL]+ (BMPL = 1-butyl-1-methylpyrrolidinium) have been synthesized. Their solid-state structures have been elucidated by single-crystal X-ray diffraction and their physicochemical properties (thermal behaviour and stability, dynamic viscosity and specific conductivity) have been assessed. In addition, the ion diffusion was studied by pulsed field gradient stimulated echo (PFGSTE) NMR spectroscopy. The decisive influence of the ring size of the cyclic sulfonimide anions on the physicochemical properties of the ILs has been revealed. All ILs show different properties compared to those of the non-cyclic TFSI anion. While these differences are especially distinct for ILs with the very rigid 6cPFSI anion, the 5-membered ring anion 5cPFSI was found to result in ILs with relatively similar properties. The difference between the properties of the TFSI anion and the cyclic sulfonimide anions has been rationalized by the rigidity (conformational lock) of the cyclic sulfonimide anions. The comparison of selected IL properties was augmented by MD simulations. These highlight the importance of π+-π+ interactions between pairs of [EMIm]+ cations in the liquid phase. The π+-π+ interactions are evident for the solid state from the molecular structures of the [EMIm]+-ILs with the three cyclic imide anions determined by single-crystal X-ray diffraction.

Cumulative Millisecond-Long Sampling for a Comprehensive Energetic Evaluation of Aqueous Ionic Liquid Effects on Amino Acid Interactions

The interactions of amino acid side-chains confer diverse energetic contributions and physical properties to a protein's stability and function. Various computational tools estimate the effect of changing a given amino acid on the protein's stability based on parametrized (free) energy functions. When parametrized for the prediction of protein stability in water, such energy functions can lead to suboptimal results for other solvents, such as ionic liquids (IL), aqueous ionic liquids (aIL), or salt solutions. However, to our knowledge, no comprehensive data are available describing the energetic effects of aIL on intramolecular protein interactions. Here, we present the most comprehensive set of potential of mean force (PMF) profiles of pairwise protein-residue interactions to date, covering 50 relevant interactions in water, the two biotechnologically relevant aIL [BMIM/Cl] and [BMIM/TfO], and [Na/Cl]. These results are based on a cumulated simulation time of >1 ms. aIL and salt ions can weaken, but also strengthen, specific residue interactions by more than 3 kcal mol^-1, depending on the residue pair, residue-residue configuration, participating ions, and concentration, necessitating considering such interactions specifically. These changes originate from a complex interplay of competitive or cooperative noncovalent ion-residue interactions, changes in solvent structural dynamics, or unspecific charge screening effects and occur at the contact distance but also at larger, solvent-separated distances. This data provide explanations at the atomistic and energetic levels for complex IL effects on protein stability and should help improve the prediction accuracies of computational tools that estimate protein stability based on (free) energy functions.

Placing gold on a π+-surface: ligand design and impact on reactivity

We describe a novel gold chloride complex supported by an ambiphilic phosphine/xanthylium ligand in which the AuCl moiety interacts with the π⁺ surface of the xanthylium unit as indicated by structural studies. Energy decomposition analyses carried out on a model system indicates the prevalence of non-covalent interactions in which the electrostatic and dispersion terms cumulatively dominate. The presence of these AuCl-π⁺ interactions correlates with the high catalytic activity of this complex in the cyclisation of 2-(phenylethynyl)phenylboronic acid, N-propargyl-t-butylamide, and 2-allyl-2-(2-propynyl)malonate. Comparison with the significantly less active acridinium and the 9-oxa-10-boraanthracene analogues reinforces this conclusion.

Infrared and Terahertz Spectroscopic Investigation of Imidazolium, Pyridinium, and Tetraalkylammonium Tetrafluoroborate Ionic Liquids

We performed terahertz time-domain spectroscopy and infrared spectroscopy of imidazolium-based, pyridinium-based, and tetraalkylammonium-based tetrafluoroborate ionic liquids to study their characteristic intermolecular and intramolecular vibrational modes to clarify interactions between various cations and the tetrafluoroborate anion. It was found that the central frequency of the intermolecular vibrational band for these ionic liquids has a relatively high frequency, ranging from 90 to 100 cm^-1. In the 900-1150 cm^-1 range, the intramolecular vibrational absorption band of the 3-fold degenerate mode of tetrafluoroborate anions in the ionic liquids was observed. Although the tetrafluoroborate anion is attributable to one of the weakly coordinated anions, the spectroscopic splitting behavior of the 3-fold degenerate mode differs depending on the cation species. It was revealed that the degenerate mode is very sensitive to local interactions between the tetrafluoroborate anion and each cation.

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.

Molecular-Level Understanding of Surface Roughness Boosting Segregation Behavior at the ZIF-8/Ionic Liquid Interfaces

Here, we perform a series of classical molecular dynamics simulations for two different [HEMIM][DCA] and [BMIM][BF₄] ionic liquids (ILs) on the ZIF-8 surface to explore the interfacial properties of metal-organic framework (MOFs)/IL composite materials at the molecular level. Our simulation results reveal that the interfacial structures of anions and cations on the ZIF-8 surface are dominated by the surface roughness due to the steric hindrance, which is extremely different from the driving mechanism based on solid-ion interactions of ILs on flat solid surfaces. At the ZIF-8/IL interfaces, the open sodalite (SOD) cages of the ZIF-8 surface can block most of the large-size cations outside and significantly boost the segregation behavior of anions and cations. In comparison with the [BMIM][BF₄] IL, the [HEMIM][DCA] IL has much more anions entering into the open SOD cages owing to the combination of stronger ZIF-8-[DCA]^- interactions and more ordered arrangement of [DCA]^- anions on the ZIF-8 surface. Furthermore, more and stronger ZIF-8-[BF₄]^- hydrogen bonds (HBs) are found to exist on the cage edges than the ZIF-8-[DCA]^- HBs, further preventing [BF₄]^- anions from entering into SOD cages. By more detailed analyses, we find that the hydrophobic interaction has an important influence on the interfacial structures of the side chains of [HEMIM]⁺ and [BMIM]⁺ cations, while the π-π stacking interaction plays a key role in determining the interfacial structures of the imidazolium rings of both cations. Our simulation results in this work provide a molecular-level understanding of the underlying driving mechanism on segregation behavior at the ZIF-8/IL interfaces.

Enhancement of charge-assisted hydrogen bond capabilities due to O-alkylation proximity in alkoxy cationic polythiophenes: solution- and solid-state evidence via EPR, AFM and surface free energy

Despite the array of applications for cationic polythiophenes (CPTs), there is still a need for structure-function guidelines and mechanistic understanding of their solution- and solid-state properties. This work presents a solution- and solid-state investigation of the effect of O-alkylation proximity on the hydrogen bonding (H-bonding) capabilities of alkoxy-CPTs, based on comparing an imidazolium alkoxy CPT with strong cation-pi, pi+ and positive charge-assisted hydrogen bonding (+CAHB) capabilities (PIMa), with two isothiouronium alkoxy CPTs with two-point +CAHB capabilities (PT1 & PT2), which have short and long alkoxy side chains, respectively. Our results show that a closer proximity of O-alkylation strengthens the +CAHB capabilities of PT1: in aqueous solutions, PT2 aggregates have a stronger interaction with cationic EPR spin probes than aggregates of PIMa and PT1, which in turn show a similar extent of repulsion towards the cationic spin probes. In solid-state, atomic force microscopy (AFM) shows that PIMa generates dendritic structures onto mica, with features of diffusion-limited aggregation (DLA), indicating strong interactions with the anionic substrate due to a high configurational entropy during spreading, regardless of being drop-casted from water or 1,4-dioxane-water (W-DI), despite the latter disturbing H-bonding due to selective solvation. PT1 is also capable of generating dendritic structures resembling ballistic aggregation (BA). However, this occurs only when casting from water, since W-DI generates island-like aggregates resembling attachment limited aggregation (ALA), which is the morphology generated by PT2 regardless of the solvent. Finally, spin-coated films of PIMa and PT1 show similar dispersivity of the surface free energy (SFE), which in turn is larger than that in PT2 films, which are also more affected when casted from W-DI, presenting much larger decreases of dispersivity. These results constitute a novel empirical structure-function guideline that could be useful for optimal design and/or processing of alkoxy CPTs. For example, dendritic patterns have recently gained attention since the colloidal droplet drying is related to engineering applications including inkjet printing, biosensing, and functional material design, while the SFE is relevant for opto- and bio-electronic applications of conjugated polyelectrolytes (CPEs). This information could also be useful when analyzing previous results obtained from alkoxy CPTs with different side chain lengths.

Isolation and reactivity of a gold(I) hydroxytrifluoroborate complex stabilized by anion-π+ interactions

A 9,9-dimethylxanthene-based ligand substituted at the 4- and 5-positions by a phosphine and a xanthylium unit, respectively, has been prepared and converted into an AuCl complex, the structure of which reveals an intramolecular Au-Cl⋯π⁺ interaction. This new ligand platform was also found to support the formation of an unprecedented hydroxytrifluoroborate derivative featuring a "hard/soft" mismatched Au-μ(OH)-BF₃ motif. Despite its surprising stability, this gold hydroxytrifluoroborate complex is a remarkably potent carbophilic catalyst which readily activates alkynes, without activator.

Pressure-Dependent Clustering in Ionic-Liquid-Poly (Vinylidene Fluoride) Mixtures: An Infrared Spectroscopic Study

The nanostructures of ionic liquids (ILs) have been the focus of considerable research attention in recent years. Nevertheless, the nanoscale structures of ILs in the presence of polymers have not been described in detail at present. In this study, nanostructures of ILs disturbed by poly(vinylidene fluoride) (PVdF) were investigated via high-pressure infrared spectra. For 1-(2-hydroxyethyl)-3-methylimidazolium bis(trifluoromethanesulfonyl)imide ([HEMIm][TFSI])-PVdF mixtures, non-monotonic frequency shifts of the C^4,5-H vibrations upon dilution were observed under ambient pressure. The experimental results suggest the presence of microheterogeneity in the [HEMIm][TFSI] systems. Upon compression, PVdF further influenced the local structure of C^4,5–H via pressure-enhanced IL–PVdF interactions; however, the local structures of C²–H and hydrogen-bonded O–H were not affected by PVdF under high pressures. For choline [TFSI]–PVdF mixtures, PVdF may disturb the local structures of hydrogen-bonded O–H. In the absence of the C^4,5–H⋯anion and C²–H⋯anion in choline [TFSI]–PVdF mixtures, the O–H group becomes a favorable moiety for pressure-enhanced IL–PVdF interactions. Our results indicate the potential of high-pressure application for designing pressure-dependent electronic switches based on the possible changes in the microheterogeneity and electrical conductivity in IL-PVdF systems under various pressures.

A force field for bio-polymers in ionic liquids (BILFF) - part 1: [EMIm][OAc]/water mixtures

We present BILFF, a novel force field for bio-polymers in ionic liquids. In the first part of our study, we introduce optimized force field parameters for mixtures of the ionic liquid (IL) 1-ethyl-3-methylimidazolium acetate ([EMIm][OAc]) with water. This imidazolium-based IL is of particular practical importance as it can dissolve significant amounts of cellulose even at room temperature. An understanding of this dissolution process via molecular dynamics simulations requires a quantitative description of the microscopic structure and the strong hydrogen bonds with a method able of simulating at least several dozen nanoseconds, which is the main aim of our novel force field. To reach this goal, we optimize the force field parameters to reproduce radial, spatial, and combined distribution functions, hydrogen bond lifetimes, diffusion coefficients, and several other quantities from reference ab initio molecular dynamics (AIMD) simulations. Non-trivial effects such as dispersion interactions between the side chains and π-π stacking of the cations are reproduced very well. We further validate the force field by comparison to experimental data such as thermal expansion coefficients, bulk modulus, and density at different temperatures, which yields good agreement and correct trends. No other force field with optimized parameters for mixtures of [EMIm][OAc] and water has been presented in the literature yet. Optimized force field parameters for cellulose and other ILs will be published in upcoming articles.

Masking specific effects of ionic liquid constituents at the solid-liquid interface by surface functionalization

Ion specific effects of ionic liquid (IL) constituents on the surface charge and aggregation properties of two types of particles (positively charged amidine (AL) and polyimidazolium-functionalized sulfate (SL-IP-2) latexes) were investigated in IL solutions containing different anions and the 1-butyl-3-methylimidazolium cation. For the AL systems, the affinity of IL anions to the particle surface followed the sequence chloride < bromide < nitrate < acetate. The critical coagulation concentration values decreased in the same order indicating that ion specific adsorption determines the surface charge density and the extent of the repulsive interparticle forces. In contrast, no tendencies were observed for the SL-IP-2 particles, i.e., both charge and aggregation features were insensitive to the type of anions. This surprising behavior sheds light on that surface functionalization with the polyimidazolium compound effectively masks interfacial ion specific effects. These results indicate new possible routes to the design of processable particle dispersions in ILs irrespective of their composition.

A Systematic Study of DFT Performance for Geometry Optimizations of Ionic Liquid Clusters

Clusters of two ion pairs of imidazolium-based ionic liquids were optimized with 43 different levels of theory, including DFT functionals and MP2-based methods combined with varying Dunning's basis sets, and added dispersion corrections. Better preforming DFT functionals were then applied to clusters consisting of four ion pairs. Excellent performance of some DFT functionals for the two ion pair clusters did not always match that of the four ion-paired clusters despite interionic distances remaining constant between the optimized two and four ion-paired clusters of the same ionic liquid. Combinations of DFT functional and basis set such as ωB97X-D/cc-pVDZ, M06-2X/aug-cc-pVDZ, B3LYP-D3/cc-pVTZ, and TPSS-D3/cc-pVTZ gave excellent results for geometry optimization of two ion-paired clusters of imidazolium ionic liquids but gave larger deviations when applied to the four ion-paired clusters of varying ionic liquids. Empirical dispersion corrections were seen to be crucial in correctly capturing correlation effects in the studied ionic liquid clusters, becoming more important in larger clusters. Dunning's double-ζ basis set, cc-pVDZ, is associated with the smallest root mean squared deviations for geometries; however, it also produces the largest deviations in total electronic energies. ωB97X-D and M06-2X produced the best performance with the augmented version of this basis set. The triple-ζ basis set, cc-pVTZ, leads to the best performance of most of the DFT functionals (especially the dispersion-corrected ones) used, whereas its augmented version, aug-cc-pVTZ, was not seen to improve results. The combinations of functional and basis set that gave the best geometry and energetics in both two and four ion-paired clusters were PBE-D3/cc-pVTZ, ωB97X-D/aug-cc-pVDZ, and BLYP-D3/cc-pVTZ. All three combinations are recommended for geometry optimizations of larger clusters of ionic liquids. PBE-D3/cc-pVTZ performed the best with an average deviation of 2.3 kJ mol^-1 and a standard deviation of 3.4 kJ mol^-1 for total electronic energy when applied to four ion-paired clusters. Geometries optimized with FMO2-SRS-MP2/cc-pVTZ produced total energy within 2.0 kJ mol^-1 off the benchmark in two ion-paired clusters, with the cc-pVDZ basis set performing unsurprisingly poorly with the same method. The error increased to 4.8 kJ mol^-1 on average in four ion-paired clusters, with the smallest RMSD deviations in geometries when compared to the benchmark ones. This study is the first report that investigated the performance of DFT functionals for two and four ion-paired clusters of a wide range of ionic liquids consisting of commonly used cations such as pyrrolidinium, imidazolium, pyridinium, and ammonium. It also identified the importance of assessing the performance of quantum chemical methods for ionic liquids on a variety of cation-anion combinations.

Dissolving Cellulose in 1,2,3-Triazolium- and Imidazolium-Based Ionic Liquids with Aromatic Anions

We present 1,2,3-triazolium- and imidazolium-based ionic liquids (ILs) with aromatic anions as a new class of cellulose solvents. The two anions in our study, benzoate and salicylate, possess a lower basicity when compared to acetate and therefore should lead to a lower amount of N-heterocyclic carbenes (NHCs) in the ILs. We characterize their physicochemical properties and find that all of them are liquids at room temperature. By applying force field molecular dynamics (MD) simulations, we investigate the structure and dynamics of the liquids and find strong and long-lived hydrogen bonds, as well as significant π-π stacking between the aromatic anion and cation. Our ILs dissolve up to 8.5 wt.-% cellulose. Via NMR spectroscopy of the solution, we rule out chain degradation or derivatization, even after several weeks at elevated temperature. Based on our MD simulations, we estimate the enthalpy of solvation and derive a simple model for semi-quantitative prediction of cellulose solubility in ILs. With the help of Sankey diagrams, we illustrate the hydrogen bond network topology of the solutions, which is characterized by competing hydrogen bond donors and acceptors. The hydrogen bonds between cellulose and the anions possess average lifetimes in the nanosecond range, which is longer than found in common pure ILs.

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.

The Effect of Phenyl Substitutions on Microstructures and Dynamics of Tetraalkylphosphonium Bis(trifluoro- methylsulfonyl)imide Ionic Liquids

Extensive atomistic simulations demonstrated that a gradual substitution of hexyl chains with phenyl groups in tetraalkylphosphonium cations results in remarkable changes in hydrogen bonding interactions, liquid structures and scattering structural functions, and rotational dynamics of hexyl chains and phenyl groups in tetraalkylphosphonium bis(trifluoromethylsulfonyl)imide ionic liquids. Hydrogen donor sites in hexyl chains present competitive characteristics with those in phenyl groups in coordinating anions, as well as their continuous and intermittent hydrogen bonding dynamics. Cation-cation and anion-anion spatial correlations show concomitant shift to short distances with decreased peak intensities with variations of cation structures, whereas cation-anion correlations have a distinct shift to large radial distances due to decreased associations of anions with neighboring cations. These microstructural changes are qualitatively manifested in shifts of prominent peaks for prevalent charge alternations and adjacency correlations between ion species in scattering structural functions. Meanwhile, rotational dynamics of hexyl chains speed up, which, in turn, slow down rotations of phenyl groups, whereas anions exhibit imperceptible changes in their rotational dynamics. These computational results are intrinsically correlated with conformational flexibilities, molecular sizes, and steric hindrance effects of phenyl groups in comparison with hexyl chains, and constrained distributions of anions around cations in heterogeneous ionic environments.

One-pot synthesis of symmetric imidazolium ionic liquids N,N-disubstituted with long alkyl chains

The Debus–Radziszewski imidazole synthesis was adapted to directly yield long-chain imidazolium ionic liquids. Imidazolium acetate ionic liquids with side-chains up to sixteen carbon atoms were synthesised in excellent yields via an on-water, one-pot reaction. The imidazolium acetate ILs acted as surfactants when dissolved in various solvents. The imidazolium acetate ionic liquids were also derivatised via an acid metathesis to the chloride, nitrate, and hydrogen oxalate derivatives. The thermal behaviour of all the ionic liquids was determined via thermogravimetric and calorimetric analysis.

The modified Debus–Radziszewski reaction was used as a one-pot on-water reaction to allow a greener synthesis of long-chain 1,3-dialkylimidazolium acetate ionic liquids in high yield from long-chain linear amines.

Synthesis and Performance of Aromatic Polyamide Ionenes as Gas Separation Membranes

Here, we report the synthesis and thermophysical properties of seven primarily aromatic, imidazolium-based polyamide ionenes. The effects of varied para-, meta-, and ortho-connectivity, and spacing of ionic and amide functional groups, on structural and thermophysical properties were analyzed. Suitable, robust derivatives were cast into thin films, neat, or with stoichiometric equivalents of the ionic liquid (IL) 1-benzy-3-methylimidazolium bistriflimide ([Bnmim][Tf2N]), and the gas transport properties of these membranes were measured. Pure gas permeabilities and permselectivities for N2, CH4, and CO2 are reported. Consistent para-connectivity in the backbone was shown to yield the highest CO2 permeability and suitability for casting as a very thin, flexible film. Derivatives containing terephthalamide segments exhibited the highest CO2/CH4 and CO2/N2 selectivities, yet CO2 permeability decreased with further deviation from consistent para-linkages.

Cationic Imidazolium Polythiophenes: Effects of Imidazolium-Methylation on Solution Concentration-Driven Aggregation and Surface Free Energy of Films Processed from Solvents with Different Polarity

Cationic imidazolium-functionalized polythiophenes with single- or double-methylation of the imidazolium ring were used to study the impact of imidazolium-methylation on (i) the solution concentration-driven aggregation in the presence of paramagnetic probes with different ionic and hydrophobic constituents and (ii) their surface free energy (SFE) as spin-coated films deposited on plasma-activated glass. Electron paramagnetic resonance spectroscopy shows that the differences in film structuration between the polymers with different methylations originate from the early stages of aggregation. In the solid state, higher degree of imidazolium-methylation generates smaller values of total SFE, γS, (by around 2 mN/m), which could be relevant in optoelectronic applications. Methylation also causes a decrease in the polar contribution of γS (γSp), suggesting that methylation decreases the polar nature of the imidazolium ring, probably due to the blocking of its H-bonding capabilities. The values of γS obtained in the present work are similar to the values obtained for doped films of neutral conjugated polymers, such as polyaniline, poly(3-hexylthiophene), and polypyrrole. However, imidazolium-polythiophenes generate films with a larger predominance of the dispersive component of γS (γSd), probably due to the motion restriction in the ionic functionalities in a conjugated polyelectrolyte, in comparison to regular dopants. The presence of 1,4-dioxane increases γSp, especially, in the polymer with larger imidazolium-methylation (and therefore unable to interact through H-bonding), probably by a decrease of the imidazolium-glass interactions. Singly-methylated imidazolium polythiophenes have been applied as electrode selective ("buffer") interlayers in conventional and inverted organic solar cells, improving their performance. However, clear structure-function guidelines are still needed for designing high-performance polythiophene-based interlayer materials. Therefore, the information reported in this work could be useful for such applications.

How Molecular Chiralities of Bis(mandelato)borate Anions Affect Their Binding Structures With Alkali Metal Ions and Microstructural Properties in Tetraalkylphosphonium Ionic Liquids

Spiroborate anion-based inorganic electrolytes and ionic liquids (ILs) have fascinating electrochemical and tribological properties and have received widespread attention in industrial applications. The molecular chiralities of spiroborate anions have a significant effect on the microstructures and macroscopic functionalities of these ionic materials in application and thus deserve fundamental consideration. In the current work, we performed quantum chemistry calculations to address the binding strength and coordination structures of chiral bis(mandelato)borate ([BMB]) anions with representative alkali metal ions, as well as the electronic properties of alkali metal ion-[BMB] ion pair complexes. The optimized [BMB] conformers are categorized into V-shaped, bent, and twisted structures with varied electrostatic potential contours and conformational energies and distinct alkali metal ion-[BMB] binding structures. Alkali metal ions have additional associations with phenyl groups in V-shaped [BMB] conformers owing to preferential cation-π interactions. Furthermore, the effects of the molecular chiralities of [BMB] anions on the thermodynamics and microstructural properties of tetraalkylphosphonium [BMB] ILs were studied by performing extensive atomistic interactions. Oxygen atoms in [BMB] anions have competitive hydrogen bonding interactions with hydrogen atoms in cations depending on the molecular chiralities and steric hindrance effects of [BMB] anions. However, the molecular chiralities of [BMB] anions have a negligible effect on the liquid densities of tetraalkylphosphonium [BMB] ILs and the spatial distributions of boron atoms in anions around phosphorous atoms in cations. Enlarging tetraalkylphosphonium cation sizes leads to enhanced cation-anion intermolecular hydrogen bonding and Coulombic interactions due to enhanced segregation of polar groups in apolar networks in heterogeneous IL matrices, as verified by scattering structural functions.

Studies on the Solubility of Terephthalic Acid in Ionic Liquids

Low solubility of terephthalic acid in common solvents makes its industrial production very difficult and not environmentally benign. Ionic liquids are known for their extraordinary solvent properties, with capability to dissolve a wide variety of materials, from common solvents to cellulose, opening new possibilities to find more suitable solvents for terephthalic acid. This work presents studies on the solubility of terephthalic acid in ionic liquids, and demonstrates that terephthalic acid is soluble in ionic liquids, such as 1-ethyl-3-methylimidazolium diethylphosphate, 1-butyl-3-methylimidazolium acetate, and dialkylimidazolium chlorides up to four times higher than in DMSO. Additionally, the temperature effect and correlation of ionic liquid structure with solubility efficiency are discussed.

The interplay and strength of the π⋯HF, C⋯HF, F⋯HF and F⋯HC hydrogen bonds upon the formation of multimolecular complexes based on C2H2⋯HF and C2H4⋯HF small dimers

The conception of this theoretical research was idealized aiming to unveil the intermolecular structures of complexes formed by acetylene or ethylene and hydrofluoric acid. At light of computational calculations by using the B3LYP/6-311++G(d,p) method, the geometries of the C₂H₂⋯(HF), C₂H₂⋯2(HF), C₂H₂⋯4(HF), C₂H₄⋯(HF), C₂H₄⋯2(HF) and C₂H₄⋯4(HF) hydrogen-bonded complexes were fully optimized. Moreover, the Post-Hartree-Fock calculations MP2/6-311++G(d,p), MP2/aug-cc-pVTZ, MP4(SDQ)/6-311++G(d,p) and CCSD/6-311++G(d,p) also were also used. The infrared spectra were analyzed in order to identify the new vibrational modes and frequencies of the proton donors shifted to red region. Through the modeling of charge-fluxes on the basis of the Quantum Theory of Atoms In Molecules (QTAIM) and, by contradicting the expectation of the hydrofluorination mechanisms of acetylene or ethylene, C⋯HF was recognized as a new type of hydrogen bond instead of the already well known π⋯H. The calculations of the Natural Bonding Orbital (NBO) and Charges derived from the Electrostatic Potential Grid-based (ChElPG) were also applied to interpret the shifting frequencies as well as measuring of the punctual charge-transfer after the formation of the complexes. Finally, the determination of the stabilization energy was carried out through the arguments of the Fock matrix in NBO basis and through the supermolecule approach. Also it is worthwhile to notice that some algebraic formulations were used for determining the electronic cooperative effect (CE).

Ionic liquid-based nanofluids (ionanofluids) for thermal applications: an experimental thermophysical characterization

Heat transfer fluids materials are manufactured for the purpose of transfer, distribution and storage of heat. Several of their important properties can be listed (for example flash point, thermal expansivity or technical safety). However, to assess the thermal exchange performance of these fluids, a prior knowledge of their heat capacity, density, viscosity and thermal conductivity is obligatory. The most popular heat transfer fluids are based on organic liquids, such as ethylene glycol. However, new technologies and development require more efficient materials. Ionanofluids, mixtures of ionic liquids and nanoparticles, were proposed as a viable replacement for those commonly used fluids due to the properties of ionic liquids (wide liquid range or low vapour pressure and flammability) combined with enhanced thermophysical properties of nanofluids caused by the dispersion of nanoparticles (mainly thermal conductivity and heat capacity). Very few authors reported the extensive analysis of those systems thermophysical properties and impact on the heat exchange efficiency. Moreover, the availability of published data is very limited. The aim of this work is to investigate ionanofluids based on the trihexyl(tetradecyl)phosphonium cation paired with the acetate, butanoate, hexanoate, octanoate or decanoate anion, mixed with carbon nanotubes, boron nitride, graphite or mesoporous carbon as nanoparticles with concentration up to 3 wt %. The density, heat capacity, thermal stability, thermal conductivity and viscosity of selected ionanofluids were determined experimentally as functions of the temperature (up to 363.15 K) and compared with theoretical tools to evaluate the predictive capability. Based on the experimental results, lubrication, heat storage potential and economic analysis were also discussed and compared to commercial heat transfer fluids.

A concerted addition mechanism in [Hmim]Br-triggered thiol–ene reactions: a typical “ionic liquid effect” revealed by DFT and experimental studies

The π⁺–π and H-bond interactions between [Hmim]Br and substrates promote a special one-step addition mechanism in thiol–ene reactions.

Systematic study on anion–cation interactions via doubly ionic H-bonds in 1,3-dimethylimidazolium salts comprising chalcogenolate anions MMIm [ER] (E = S, Se; R = H, tBu, SiMe3)

We present convenient syntheses of so far inaccessible, crystalline and highly pure 1,3-dialkylimidazolium salts with extremely nucleophilic thiolate and selenolate anions [ER]⁻ (R = H, tBu, SiMe₃).