Understanding the mechanism of cellulose dissolution in 1-butyl-3-methylimidazolium chloride ionic liquid via quantum chemistry calculations and molecular dynamics simulations

While N,N′-dialkylimidazolium ionic liquids (ILs) have been well-established as effective solvents for dissolution and processing of cellulose, the detailed mechanism at the molecular level still remains unclear. In this work, we present a combined quantum chemistry and molecular dynamics simulation study on how the ILs dissolve cellulose. On the basis of calculations on 1-butyl-3-methylimidazolium chloride, one of the most effective ILs dissolving cellulose, we further studied the molecular behavior of cellulose models (i.e. cellulose oligomers with degrees of polymerization n = 2, 4, and 6) in the IL, including the structural features and hydrogen bonding patterns. The collected data indicate that both chloride anions and imidazolium cations of the IL interact with the oligomer via hydrogen bonds. However, the anions occupy the first coordination shell of the oligomer, and the strength and number of hydrogen bonds and the interaction energy between anions and the oligomer are much larger than those between cations and the oligomer. It is observed that the intramolecular hydrogen bond in the oligomer is broken under the combined effect of anions and cations. The present results emphasize that the chloride anions play a critically important role and the imidazolium cations also present a remarkable contribution in the cellulose dissolution. This point of view is different from previous one that only underlines the importance of the chloride anions in the cellulose dissolution. The present results improve our understanding for the cellulose dissolution in imidazolium chloride ILs.     

Ionic liquid processing of cellulose

Utilization of natural polymers has attracted increasing attention because of the consumption and over-exploitation of non-renewable resources, such as coal and oil. The development of green processing of cellulose, the most abundant biorenewable material on Earth, is urgent from the viewpoints of both sustainability and environmental protection. The discovery of the dissolution of cellulose in ionic liquids (ILs, salts which melt below 100 °C) provides new opportunities for the processing of this biopolymer, however, many fundamental and practical questions need to be answered in order to determine if this will ultimately be a green or sustainable strategy. In this critical review, the open fundamental questions regarding the interactions of cellulose with both the IL cations and anions in the dissolution process are discussed. Investigations have shown that the interactions between the anion and cellulose play an important role in the solvation of cellulose, however, opinions on the role of the cation are conflicting. Some researchers have concluded that the cations are hydrogen bonding to this biopolymer, while others suggest they are not. Our review of the available data has led us to urge the use of more chemical units of solubility, such as 'g cellulose per mole of IL' or 'mol IL per mol hydroxyl in cellulose' to provide more consistency in data reporting and more insight into the dissolution mechanism. This review will also assess the greenness and sustainability of IL processing of biomass, where it would seem that the choices of cation and anion are critical not only to the science of the dissolution, but to the ultimate 'greenness' of any process (142 references).     

Solvatochromic studies of ionic liquid/organic mixtures

Room-temperature ionic liquids (ILs) have potential for many different applications, including catalysis and synthesis. Organics are often present during IL applications; therefore, a more fundamental understanding of the interactions between IL and organics is necessary. A systematic study of the effects of organic cosolvents, cations, and anions on the solvent strength of IL/organic mixtures will allow for a greater understanding and potential for tuning of ILs for specific purposes. Solvent strength is commonly quantified using spectroscopic probes. We report the solvent strength of IL/organic mixtures using Reichardt's dyes 30 and 33, Kamlet-Taft parameters, and phenol blue. The results show that the polarity of ILs is largely unaffected by the organic cosolvent; that is, the probes are preferentially solvated by the ILs. However, more specific solvation forces, such as hydrogen bonding, can be influenced indirectly by the strength of the anion/cation interaction, giving counterintuitive results.     

Dual effects of dimethylsulfoxide on cellulose solvating ability of 1-allyl-3-methylimidazolium chloride

Polar aprotic solvents are considered to act as cosolvents with ionic liquids (ILs) for cellulose, strengthening the solvating ability of ILs by improving their cellulose solvating kinetics without influencing the solubility of cellulose in ILs. In this work, it was found that dimethylsulfoxide (DMSO) at low concentration improves the cellulose solvating ability of [AMIM][Cl], but weakens it at high concentration. To clarify the mechanism of these dual effects of DMSO on the cellulose solvating ability of [AMIM][Cl], the [AMIM][Cl]/DMSO system was investigated using excess infrared spectroscopy, nuclear magnetic resonance (NMR) T ₂ relaxometry, ¹H NMR, ³⁵Cl NMR, and dynamic light scattering. The results indicate that the tight association between the cation and anion in the [AMIM][Cl] network is loosened at low DMSO concentration. As a result, mass transport is accelerated due to the enhanced dynamics of [AMIM][Cl], promoting the cellulose solvating kinetics of [AMIM][Cl]. However, ion clusters of [AMIM][Cl] start to form when the molar fraction of DMSO (x _DMSO) exceeds 0.5. The hydrogen bonds between cations and anions in the ion clusters become much stronger than in pure [AMIM][Cl], leading to decreased ability of [AMIM][Cl] to form hydrogen bonds with cellulose and thus decreased cellulose solubility in the [AMIM][Cl]/DMSO mixture.     

Hydrogen bonding mediated ion pairs of some aprotic ionic liquids and their structural transition in aqueous solution

Ion pair speciation of ionic liquids(ILs) has an important effect on the physical and chemical properties of ILs and recognition of the structure of ion pairs in solution is essential. It has been reported that ion pairs of some ILs can be formed by hydrogen bonding interactions between cations and anions of them. Considering the fact that far-IR(FIR) spectroscopy is a powerful tool in indicating the intermolecular and intramolecular hydrogen bonding, in this work, this spectroscopic technique has been combined with molecular dynamic(MD) simulation and nuclear magnetic resonance hydrogen spectroscopy(~1H NMR) to investigate ion pairs of aprotic ILs [Bmim][NO_3], [BuPy][NO_3], [Pyr_(14)][NO_3], [PP_(14)][NO_3] and [Bu-choline][NO_3] in aqueous IL mixtures. The FIR spectra have been assigned with the aid of density functional theory(DFT) calculations, and the results are used to understand the effect of cationic nature on the structure of ion pairs. It is found that contact ion pairs formed in the neat aprotic ILs by hydrogen bonding interactions between cation and anion, were still maintained in aqueous solutions up to high water mole fraction(say 0.80 for [BuPy][NO3]). When water content was increased to a critical mole fraction of water(say 0.83 for [BuPy][NO3]), the contact ion pairs could be transformed into solvent-separated ion pairs due to the formation of the hydrogen bonding between ions and water. With the further dilution of the aqueous ILs solution, the solvent-separated ion pairs was finally turned into free cations and free anions(fully hydrated cations or anions). The concentrations of the ILs at which the contact ion pairs were transformed into solvent-separated ion pairs and solvent-separated ion pairs were transformed into free ions(fully hydrated ion) were dependent on the cationic structures. These information provides direct spectral evidence for ion pair structures of the aprotic ILs in aqueous solution. MD simulation and ~1H NMR results support the conclusion drawn from FIR spectra investigations.     

Glucose solvation by the ionic liquid 1,3-dimethylimidazolium chloride: a simulation study

Simulations of beta-glucose in the ionic liquid 1,3-dimethylimidazolium chloride have been performed in order to examine the solvation environment of the carbohydrate. Both single molecule and 1:5 glucose:ionic liquid (16.7 wt %) solutions are studied, and the hydrogen bonding between sugar and solvent is examined. The primary solvation shell around the perimeter of the glucose ring consists predominantly of chloride anions which hydrogen bond to the hydroxyl groups. A small presence of the cation is also found, with the association occurring through the weakly acidic hydrogen at the 2-position of the imidazolium ring interacting with the oxygen atoms of the sugar secondary hydroxyls. An average chloride coordination number of 4 is found around the glucose for both the single molecule and high concentration simulations, despite the reduced chloride:glucose ratio in the latter case. In relation to the cation, the glucose molecules occupy positions above and below the plane of the imidazolium ring. Importantly, even at high glucose concentrations, no significant change in the anion-cation interactions and overall liquid structure of the ionic liquid is found, indicating that the glucose is readily accommodated by the solvent at this concentration. Dominant contributions to the sugar-ionic liquid interaction energy come from favorable hydrogen bonding (electrostatic) interactions between hydroxyls and chlorides, although a small favorable van der Waals energy contribution is also seen between the sugar and cations suggesting that the cation could be tailored in order to further improve the dissolution of glucose/cellulose in ionic liquid systems.     

Molecular dynamics simulation of dual amino-functionalized imidazolium-based ionic liquids

In this work, the all-atom (AA) force fields were set up for three kinds of dual amino-functionalized imidazolium-based ionic liquids (ILs), composed by cations with different alkyl chain length and amino acid anion [Gly]⁻. The force field was based on our previous work and the default parameters were developed in this study. Molecular dynamics simulations were performed. Validation was carried out by comparing simulation densities with experimental data, and good agreement was obtained. Molar volume and heat capacity at constant pressure were predicted. Mean square displacements for these ILs were computed and these ILs were proved to move very slowly. It may be caused by hydrogen-bonded network between ions and the terminal azyl. To depict the microscopic structures of the ILs, many types of radial distribution functions were investigated. It is interesting to find that not only the cation and anion, but also the anions themselves will form hydrogen bonds.     

Cluster Formation through Hydrogen Bond Bridges across Chloride Anions in a Hydroxyl-Functionalized Ionic Liquid

Several recent studies of hydroxyl-functionalized ionic liquids (ILs) have shown that cation-cation interactions can be dominating these materials at the molecular level when the anion involved is weakly interacting. The hydrogen bonds between the like ions led to the formation of interesting chain-like, ring-like, or distinct dimeric (i. e. two ion pairs) supermolecular clusters. In the present work, vibrational spectroscopy (ATR-IR and Raman) and density functional theory (DFT) calculations of the hydroxyl-functionalized imidazolium ionic liquid C₂OHmimCl indicate that anion-cation hydrogen bonding interactions are dominating, leading to the formation of distinct dimeric ion pair clusters. In this arrangement, the Cl⁻ anions function as a bridge between the cations by establishing bifurcated hydrogen bonds with the OH group of one cation and the C(2)-H of another cation. Cation–cation interactions, on the other hand, do not play a significant role in the observed clusters.     

Dispersion and Hydrogen Bonding Rule: Why the Vaporization Enthalpies of Aprotic Ionic Liquids Are Significantly Larger than those of Protic Ionic liquids

It is well known that gas‐phase experiments and computational methods point to the dominance of dispersion forces in the molecular association of hydrocarbons. Estimates or even quantification of these weak forces are complicated due to solvent effects in solution. The dissection of interaction energies and quantification of dispersion interactions is particularly challenging for polar systems such as ionic liquids (ILs) which are characterized by a subtle balance between Coulomb interactions, hydrogen bonding, and dispersion forces. Here, we have used vaporization enthalpies, far‐infrared spectroscopy, and dispersion‐corrected calculations to dissect the interaction energies between cations and anions in aprotic (AILs), and protic (PILs) ionic liquids. It was found that the higher total interaction energy in PILs results from the strong and directional hydrogen bonds between cation and anion, whereas the larger vaporization enthalpies of AILs clearly arise from increasing dispersion forces between ion pairs.     

Ionic liquid based on α-amino acid anion and N7,N9-dimethylguaninium cation ([dMG][AA]): theoretical study on the structure and electronic properties

The interactions between five amino acid based anions ([AA](-) (AA = Gly, Phe, His, Try, and Tyr)) and N7,N9-dimethylguaninium cation ([dMG](+)) have been investigated by the hybrid density functional theory method B3LYP together with the basis set 6-311++G(d,p). The calculated interaction energy was found to decrease in magnitude with increasing side-chain length in the amino acid anion. The interaction between the [dMG](+) cation and [AA](-) anion in the most stable configurations of ion pairs is a hydrogen bonding interaction. These hydrogen bonds (H bonds) were analyzed by the quantum theory of atoms in molecules (QTAIM) and natural bond orbital (NBO) analysis. Finally, several correlations between electron densities in bond critical points of hydrogen bonds and interaction energy as well as vibrational frequencies in the most stable configurations of ion pairs have been checked.     

Vibrational spectroscopy and dynamics of small anions in ionic liquid solutions

Fourier-transform infrared (FTIR) and time-resolved IR spectroscopies have been used to study vibrational band positions, vibrational energy relaxation (VER) rates, and reorientation times of anions in several ionic liquid (IL) solutions. The ILs primarily investigated are based on the 1-butyl-2,3-dimethylimidazolium ([BM(2)IM]) cation with thiocyanate (NCS-), dicyanamide (N(CN)2-), and tetrafluoroborate (BF4-) anions. Spectroscopic studies are carried out near 2000 cm-1 for the C[Triple Bond]N stretching bands of NCS- and N(CN)2- as the IL anion as well as for NCS-, N(CN)2-, and azide (N3-) anions dissolved in [BM2IM][BF4]. The VER studies of N(CN)2- are reported for the first time. VER of N3-, NCS-, and N(CN)2- is measured in normal solvents, such as N-methylformamide, to compare with the IL solutions. The spectral shifts and VER rates of the anions in IL solution are quite similar to those in polar aprotic, conventional organic solvents, i.e., dimethylsulfoxide, and significantly different than those in methanol, in which there is hydrogen bonding. Similar studies were also carried out for the anions in another IL, 1-butyl-3-methylimidazolium tetrafluoroborate ([BMIM][BF4]), in which the C2 hydrogen is present. The results for the anions are similar to those in the [BM2IM] containing ILs, in which the C2 hydrogen is methyl substituted. This suggests that substituting this hydrogen has, at most, a minor effect on the degree of hydrogen bonding in the anion-IL solvation interaction based on the infrared spectra and dynamics.     

The electrochemical stability of ionic liquids and deep eutectic solvents

Room temperature ionic liquids (ILs) composed of cations and anions, as well as deep eutectic solvents (DESs) composed of hydrogen bond donors (HBDs) and hydrogen bond acceptors (HBAs), are regarded as green solvents due to their low volatility. They have been used widely for electrochemically driven reactions because they exhibit high conductivity and excellent electrochemical stability. However, no systematic investigations on the electrochemical potential windows (EPWs), which could be used to characterize the electrochemical stability, have been reported. In this regard, the EPWs of 33 ILs and 23 DESs have been studied utilizing cyclic voltammetry (CV) method and the effects of structural factors (cations and anions of ILs, and HBDs and HBAs of DESs) and external factors (electrode, water content) on the EPWs have been comprehensively investigated. The electrochemical stability of selected ILs comprising five traditional cations, namely imidazolium, pyridinium, pyrrolidinium, piperidinium and ammonium and 13 kinds of versatile anions was studied. The results show that for ILs, both cation and anion play an important role on the reductive and oxidative potential limit. For a same IL at different working electrode, for example, glassy carbon (GC), gold (Au) and platinum (Pt) electrode, the largest potential window is almost observed on the GC working electrode. The investigations on the EPWs of choline chloride (ChCl), choline bromide (ChBr), choline iodide (ChI), and methyl urea based DESs show that the DES composed of ChCl and methyl urea has the largest potential window. This work may aid the selection of ILs or DESs for use as a direct electrolyte or a solvent in electrochemical applications.     

Three in One: The Versatility of Hydrogen Bonding Interaction in Halide Salts with Hydroxy-Functionalized Pyridinium Cations

The paradigm of supramolecular chemistry relies on the delicate balance of noncovalent forces. Here we present a systematic approach for controlling the structural versatility of halide salts by the nature of hydrogen bonding interactions. We synthesized halide salts with hydroxy-functionalized pyridinium cations [HOC_nPy]⁺ (n=2, 3, 4) and chloride, bromide and iodide anions, which are typically used as precursor material for synthesizing ionic liquids by anion metathesis reaction. The X-ray structures of these omnium halides show two types of hydrogen bonding: ‘intra-ionic’ H-bonds, wherein the anion interacts with the hydroxy group and the positively charged ring at the same cation, and ‘inter-ionic’ H-bonds, wherein the anion also interacts with the hydroxy group and the ring system but of different cations. We show that hydrogen bonding is controllable by the length of the hydroxyalkyl chain and the interaction strength of the anion. Some molten halide salts exhibit a third type of hydrogen bonding. IR spectra reveal elusive H-bonds between the OH groups of cations, showing interaction between ions of like charge. They are formed despite the repulsive interaction between the like-charged ions and compete with the favored cation-anion H-bonds. All types of H-bonding are analyzed by quantum chemical methods and the natural bond orbital approach, emphasizing the importance of charge transfer in these interactions. For simple omnium salts, we evidenced three distinct types of hydrogen bonds: Three in one!     

Dirubidium hexaaquacobalt(II) tetrakis(hydrogen phthalate) tetrahydrate and coordination modes of the hydrogen phthalate anion

The title compound, Rb₂[Co(H₂O)₆](C₈H₅O₄)₄·4H₂O, consists of nearly regular octahedral [Co(H₂O)₆]²⁺ cations with the Co^II cations on the inversion centre (special position 2a), Rb⁺ cations, hydrogen phthalate (Hpht⁻) anions and disordered water molecules. The Rb⁺ cation is surrounded by nine O atoms from Hpht⁻ anions and water molecules, with a strongly deformed pentagonal–bipyramidal geometry and one apex split into three positions. The crystal packing is governed by numerous hydrogen bonds involving all water molecules and Hpht⁻ anions. In this way, layers parallel to the ab plane are formed, with the aromatic rings of the Hpht⁻ anions esentially directed along the c axis. While Hpht⁻ anions form the outer part of the layers, disordered water molecules and Rb⁺ cations alternate with [Co(H₂O)₆]²⁺ cations in the inner parts. The only interactions between the layers are van der Waals forces between the atoms of the aromatic rings. A search of the Cambridge Structural Database for coordination modes and types of hydrogen‐bonding interaction of the Hpht⁻ anion showed that, when uncoordinated Hpht⁻ anions are present, compounds with intermolecular hydrogen bonds are more numerous than compounds with intramolecular hydrogen bonds. For coordinated Hpht⁻ anions, chelating and bridging anions are almost equally common, while monodentate anions are relatively scarce. The same coordination modes appear for Hpht⁻ anions with or without intramolecular hydrogen bonds, although intramolecular hydrogen bonds are less common.     

Hydrogen Bonding Between Ions of Like Charge in Ionic Liquids Characterized by NMR Deuteron Quadrupole Coupling Constants—Comparison with Salt Bridges and Molecular Systems

We present deuteron quadrupole coupling constants (DQCC) for hydroxyl‐functionalized ionic liquids (ILs) in the crystalline or glassy states characterizing two types of hydrogen bonding: The regular Coulomb‐enhanced hydrogen bonds between cation and anion (c–a), and the unusual hydrogen bonds between cation and cation (c–c), which are present despite repulsive Coulomb forces. We measure these sensitive probes of hydrogen bonding by means of solid‐state NMR spectroscopy. The DQCCs of (c–a) ion pairs and (c–c) H‐bonds are compared to those of salt bridges in supramolecular complexes and those present in molecular liquids. At low temperatures, the (c–c) species successfully compete with the (c–a) ion pairs and dominate the cluster populations. Equilibrium constants obtained from molecular‐dynamics (MD) simulations show van't Hoff behavior with small transition enthalpies between the differently H‐bonded species. We show that cationic‐cluster formation prevents these ILs from crystallizing. With cooling, the (c–c) hydrogen bonds persist, resulting in supercooling and glass formation.     

The relationship between intramolecular hydrogen bonds and certain physical properties of regioselectively substituted cellulose derivatives

This article tries to provide some direct evidence about the relationship between the intramolecular hydrogen bonds in cellulose and their corresponding effect on physical properties. The formation of intramolecular hydrogen bonds has been proved to contribute directly to certain physical properties of cellulose, such as its solubility in solvents having different polarities, the relative reactivities of the hydroxyls in a repeating unit and its crystallinity, using a 6-O-methylcellulose (6MC) film that was known¹ to have intramolecular hydrogen bonds. The excellent solubility of 6MC when compared with other cellulose derivatives indicated a lack of interchain hydrogen bonds. A comparison of the relative reactivities between the C-2 and C-3 position hydroxyls in 6MC also indicates that intramolecular hydrogen bonds once formed in 6MC films are possibly maintained even after dissolution in solvents. In addition, the poor crystallinity exhibited by 6MC supports the idea that crystallization in cellulosics may be dependent more upon preferencial interchain hydrogen bonding at the C-6 position hydroxyls than upon a uniform structure such as that found in 6MC, where every structural unit is completely and regioselectively substituted, distinguishing it from other synthetic polymers such as polyolefins and polyesters. © 1997 John Wiley & Sons, Inc. J Polym Sci B: Polym Phys, 35: 717–723, 1997     

Effects of Cationic Structure on Cellulose Dissolution in Ionic Liquids: A Molecular Dynamics Study

In recent years, great progress has been made in the dissolution of cellulose with ionic liquids (ILs). However, the mechanism of cellulose dissolution, especially the role the IL cation played in the dissolution process, has not been clearly understood. Herein, the mixtures of cellulose with a series of imidazolium‐based chloride ionic liquids and 1‐butyl‐3‐methyl pyridinium chloride ([C₄mpy]Cl) were simulated to study the effect that varying the heterocyclic structure and alkyl chain length of the IL cation has on the dissolution of cellulose. It was shown that the dissolution of cellulose in [C₄mpy]Cl is better than that in [C₄mim]Cl. For imidazolium‐based ILs, the shorter the alkyl chain is, the higher the solubility will be. In addition, an all‐atom force field for 1‐allyl‐3‐methyl imidazolium cation ([Amim]⁺) was developed, for the first time, to investigate the effect the electron‐withdrawing group within the alkyl chain of the IL cation has on the dissolution of cellulose. It was found that the interaction energy between [Amim]⁺ and cellulose was greater than that between [C₃mim]⁺ and cellulose, indicating that the presence of electron‐withdrawing group in alkyl chain of the cation enhanced the interaction between the cation and cellulose due to the increase of electronegativity of the cations. These findings are used to assess the cationic effect on the dissolution of cellulose in ILs. They are also expected to be important for rational design of novel ILs for efficient dissolution of cellulose.     

A Bond Path and an Attractive Ehrenfest Force Do Not Necessarily Indicate Bonding Interactions: Case Study on M2X2 (M=Li,Na, K; X=H,OH, F,Cl)

Interactions in dimers of model alkali metal derivatives M₂X₂ (M=Li or Na or K; X=H or F, Cl, OH) are studied in the frame of the quantum theory of atoms in molecules (QTAIM) using the interacting quantum atoms approach (IQA). Contrary to opinion prevalent in QTAIM studies, the interaction between two anions linked by a bond path is demonstrated to be strongly repulsive. One may therefore say that a bond path does not necessarily indicate bonding interactions. The interactions between two anions or two cations that are not linked by a bond path are also strongly repulsive. The repulsive anion–anion and cation–cation interactions are outweighed by much stronger attractive anion–cation interactions, and the model molecules are therefore in a stable state. The attractive Ehrenfest forces (calculated in the frame of the QTAIM) acting across interatomic surfaces shared by anions in the dimers do not reflect the repulsive interactions between anions. Probable reasons of this disagreement are discussed. The force exerted on the nucleus and the electrons of a particular atom by the nucleus and the electrons of any another atom in question is proposed. It is assumed that this force unambiguously exposes whether basins of two atoms are attracted or repelled by each other in a polyatomic molecule.